Page History

    CalRecon *calRec = rec->getCalRecon();
    
    if (!calRec) return;
 
    TObjArray*  xtalCol = calRec->getCalXtalRecCol();
    Long64_t numXtal = xtalCol->GetEntries();
    ((TH1F*)GetObjectPtr("CALXTALCOUNT"))->Fill((Float_t)numXtal);
    float totXE = 0.;
    for (int xc=0;xc<numXtal; xc++) {
      CalXtalRecData* xtal = (CalXtalRecData*)xtalCol->At(xc);
      float xtalEnergy = xtal->getEnergy();
      ((TH1F*)GetObjectPtr("CALXTALE"))->Fill(xtalEnergy);
      if (xtalEnergy > 2000) {
	CalXtalId id = xtal->getPackedId();
	int lyr = id.getLayer();
	int twr = id.getTower();
	int col = id.getColumn();
	CalRangeRecData* rData = xtal->getRangeRecData(0);
	int range = rData->getRange(0);
	double ph0 = xtal->getEnergySelectedRange(range,0);
	std::cout << "twr, lyr, col = " << twr << " " << lyr << " " << col <<
	  " range " << range << " ph0 " << ph0 << " E " << xtalEnergy << std::endl;

	continue;
      }
      totXE += xtalEnergy;
    }
    ((TH1F*)GetObjectPtr("CALXTALTOTE"))->Fill(totXE);