|
|
|
Definition at line 67 of file ValidIsoTrkCalib.cc.
| ValidIsoTrkCalib::ValidIsoTrkCalib | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 198 of file ValidIsoTrkCalib.cc.
References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and Parameters::parameters.
{
//takeAllRecHits_=iConfig.getUntrackedParameter<bool>("takeAllRecHits");
takeGenTracks_=iConfig.getUntrackedParameter<bool>("takeGenTracks");
genTracksLabel_ = iConfig.getParameter<edm::InputTag>("genTracksLabel");
genhbheLabel_= iConfig.getParameter<edm::InputTag>("genHBHE");
//genhoLabel_=iConfig.getParameter<edm::InputTag>("genHO");
//genecalLabel_=iConfig.getParameter<std::vector<edm::InputTag> >("genECAL");
// m_hcalLabel = iConfig.getUntrackedParameter<std::string> ("hcalRecHitsLabel","hbhereco");
hbheLabel_= iConfig.getParameter<edm::InputTag>("hbheInput");
hoLabel_=iConfig.getParameter<edm::InputTag>("hoInput");
//eLabel_=iConfig.getParameter<edm::InputTag>("eInput");
trackLabel_ = iConfig.getParameter<edm::InputTag>("HcalIsolTrackInput");
trackLabel1_ = iConfig.getParameter<edm::InputTag>("trackInput");
associationConeSize_=iConfig.getParameter<double>("associationConeSize");
allowMissingInputs_=iConfig.getUntrackedParameter<bool>("allowMissingInputs", true);
// outputFileName_=iConfig.getParameter<std::string>("outputFileName");
// calibFactorsFileName_=iConfig.getParameter<std::string>("calibFactorsFileName");
AxB_=iConfig.getParameter<std::string>("AxB");
calibrationConeSize_=iConfig.getParameter<double>("calibrationConeSize");
MinNTrackHitsBarrel = iConfig.getParameter<int>("MinNTrackHitsBarrel");
MinNTECHitsEndcap = iConfig.getParameter<int>("MinNTECHitsEndcap");
energyECALmip = iConfig.getParameter<double>("energyECALmip");
energyMinIso = iConfig.getParameter<double>("energyMinIso");
energyMaxIso = iConfig.getParameter<double>("energyMaxIso");
maxPNear = iConfig.getParameter<double>("maxPNear");
edm::ParameterSet parameters = iConfig.getParameter<edm::ParameterSet>("TrackAssociatorParameters");
parameters_.loadParameters( parameters );
trackAssociator_.useDefaultPropagator();
// taECALCone_=iConfig.getUntrackedParameter<double>("TrackAssociatorECALCone",0.5);
//taHCALCone_=iConfig.getUntrackedParameter<double>("TrackAssociatorHCALCone",0.6);
}
| ValidIsoTrkCalib::~ValidIsoTrkCalib | ( | ) |
Definition at line 244 of file ValidIsoTrkCalib.cc.
{
// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)
}
| void ValidIsoTrkCalib::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 255 of file ValidIsoTrkCalib.cc.
References abs, gather_cfg::cout, MaxHit_struct::depthhit, MaxHit_struct::dr, alignCSCRings::e, edm::EventID::event(), HcalObjRepresent::Fill(), cmsDownloadME::gen, GeneralTracks_cfi::generalTracks, edm::EventSetup::get(), edm::Event::getByLabel(), getDistInPlaneSimple(), reco::JetExtendedAssociation::getValue(), DetId::Hcal, HcalBarrel, MaxHit_struct::hitenergy, i, edm::EventBase::id(), HcalDetId::ieta(), MaxHit_struct::ietahitm, info, HcalDetId::iphi(), MaxHit_struct::iphihitm, n, pos, edm::ESHandle< T >::product(), edm::EventID::run(), convertSQLiteXML::runNumber, mathSSE::sqrt(), reco::btau::trackEta, reco::btau::trackPhi, TrackDetMatchInfo::trkGlobPosAtEcal, and TrackDetMatchInfo::trkGlobPosAtHcal.
{
using namespace edm;
try{
edm::ESHandle <HcalRespCorrs> recalibCorrs;
iSetup.get<HcalRespCorrsRcd>().get("recalibrate",recalibCorrs);
respRecalib = recalibCorrs.product();
LogInfo("CalibConstants")<<" Loaded: OK ";
}catch(const cms::Exception & e) {
LogWarning("CalibConstants")<<" Not Found!! ";
}
edm::Handle<reco::TrackCollection> generalTracks;
iEvent.getByLabel(genTracksLabel_, generalTracks);
edm::Handle<reco::TrackCollection> isoProdTracks;
iEvent.getByLabel(trackLabel1_,isoProdTracks);
edm::Handle<reco::IsolatedPixelTrackCandidateCollection> isoPixelTracks;
//edm::Handle<reco::TrackCollection> isoPixelTracks;
iEvent.getByLabel(trackLabel_,isoPixelTracks);
/*
edm::Handle<EcalRecHitCollection> ecal;
iEvent.getByLabel(eLabel_,ecal);
const EcalRecHitCollection Hitecal = *(ecal.product());
*/
edm::Handle<HBHERecHitCollection> hbhe;
iEvent.getByLabel(hbheLabel_,hbhe);
const HBHERecHitCollection Hithbhe = *(hbhe.product());
edm::ESHandle<CaloGeometry> pG;
iSetup.get<CaloGeometryRecord>().get(pG);
geo = pG.product();
const CaloSubdetectorGeometry* gHcal = geo->getSubdetectorGeometry(DetId::Hcal,HcalBarrel);
//Note: even though it says HcalBarrel, we actually get the whole Hcal detector geometry!
// Lumi_n=iEvent.luminosityBlock();
parameters_.useEcal = true;
parameters_.useHcal = true;
parameters_.useCalo = false;
parameters_.useMuon = false;
//parameters_.dREcal = taECALCone_;
//parameters_.dRHcal = taHCALCone_;
//cout<<"Hello World. TrackCollectionSize: "<< isoPixelTracks->size()<<endl;
if (isoPixelTracks->size()==0) return;
for (reco::TrackCollection::const_iterator trit=isoProdTracks->begin(); trit!=isoProdTracks->end(); trit++)
{
reco::IsolatedPixelTrackCandidateCollection::const_iterator isoMatched=isoPixelTracks->begin();
//reco::TrackCollection::const_iterator isoMatched=isoPixelTracks->begin();
bool matched=false;
//for (reco::IsolatedPixelTrackCandidateCollection::const_iterator trit = isoPixelTracks->begin(); trit!=isoPixelTracks->end(); trit++)
for (reco::IsolatedPixelTrackCandidateCollection::const_iterator it = isoPixelTracks->begin(); it!=isoPixelTracks->end(); it++)
//for (reco::TrackCollection::const_iterator it = isoPixelTracks->begin(); it!=isoPixelTracks->end(); it++)
{
if (abs((trit->pt() - it->pt())/it->pt()) < 0.005 && abs(trit->eta() - it->eta()) < 0.01)
{
isoMatched=it;
matched=true;
break;
}
}
// CUT
if (!matched) continue;
if(isoMatched->maxPtPxl() > maxPNear) continue;
ptNear = isoMatched->maxPtPxl();
//cout<<"Point 0.1 isoMatch. ptnear: "<<ptNear<<endl;
// CUT
if (trit->hitPattern().numberOfValidHits()<MinNTrackHitsBarrel) continue;
if (fabs(trit->eta())>1.47&&trit->hitPattern().numberOfValidStripTECHits()<MinNTECHitsEndcap) continue;
//cout<<"Point 0.2.1 after numofvalidhits HB: "<<trit->hitPattern().numberOfValidHits()<<endl;
//cout<<"Point 0.2.2 after numofvalidstrips HE: "<<trit->hitPattern().numberOfValidStripTECHits()<<endl;
numVH = trit->hitPattern().numberOfValidHits();
numVS = trit->hitPattern().numberOfValidStripTECHits();
trackE = sqrt(trit->px()*trit->px()+trit->py()*trit->py()+trit->pz()*trit->pz()+0.14*0.14);
trackPt = trit->pt();
trackEta = trit->eta();
trackPhi = trit->phi();
emEnergy = isoMatched->energyIn();
//cout<<"Point 0.3. Matched :: pt: "<<trit->pt()<<" wholeEnergy: "<<trackE<<" emEnergy: "<<emEnergy<<" eta: "<<etahcal<<" phi: "<<phihcal<<endl;
//cout<<"Point 0.4. EM energy in cone: "<<emEnergy<<" EtaHcal: "<<etahcal<<" PhiHcal: "<<phihcal<<endl;
TrackDetMatchInfo info = trackAssociator_.associate(iEvent, iSetup,trackAssociator_.getFreeTrajectoryState(iSetup, *trit), parameters_);
//float etaecal=info.trkGlobPosAtEcal.eta();
//float phiecal=info.trkGlobPosAtEcal.phi();
// float etahcal=info.trkGlobPosAtHcal.eta();
// float phihcal=info.trkGlobPosAtHcal.phi();
xTrkEcal=info.trkGlobPosAtEcal.x();
yTrkEcal=info.trkGlobPosAtEcal.y();
zTrkEcal=info.trkGlobPosAtEcal.z();
xTrkHcal=info.trkGlobPosAtHcal.x();
yTrkHcal=info.trkGlobPosAtHcal.y();
zTrkHcal=info.trkGlobPosAtHcal.z();
if (xTrkEcal==0 && yTrkEcal==0&& zTrkEcal==0) {cout<<"zero point at Ecal"<<endl; continue;}
if (xTrkHcal==0 && yTrkHcal==0&& zTrkHcal==0) {cout<<"zero point at Hcal"<<endl; continue;}
/*GlobalVector trackMomAtEcal = info.trkMomAtEcal;
GlobalVector trackMomAtHcal = info.trkMomAtHcal;
PxTrkHcal = trackMomAtHcal.x();
PyTrkHcal = trackMomAtHcal.y();
PzTrkHcal = trackMomAtHcal.z();
*/
GlobalPoint gPointEcal(xTrkEcal,yTrkEcal,zTrkEcal);
GlobalPoint gPointHcal(xTrkHcal,yTrkHcal,zTrkHcal);
int iphitrue = -10;
int ietatrue = 100;
const HcalDetId tempId = gHcal->getClosestCell(gPointHcal);
ietatrue = tempId.ieta();
iphitrue = tempId.iphi();
MaxHit_struct MaxHit;
MaxHit.hitenergy=-100.;
//container for used recHits
std::vector<int> usedHits;
//
usedHits.clear();
//cout <<"Point 1. Entrance to HBHECollection"<<endl;
//float dddeta = 1000.;
//float dddphi = 1000.;
//int iphitrue = 1234;
//int ietatrue = 1234;
GlobalPoint gPhot;
for (HBHERecHitCollection::const_iterator hhit=Hithbhe.begin(); hhit!=Hithbhe.end(); hhit++)
{
//check that this hit was not considered before and push it into usedHits
bool hitIsUsed=false;
int hitHashedIndex=hhit->id().hashed_index();
for (uint32_t i=0; i<usedHits.size(); i++)
{
if (usedHits[i]==hitHashedIndex) hitIsUsed=true;
}
if (hitIsUsed) continue;
usedHits.push_back(hitHashedIndex);
//
// rof 16.05.2008 start: include the possibility for recalibration
float recal = 1;
// rof end
GlobalPoint pos = geo->getPosition(hhit->detid());
//float phihit = pos.phi();
//float etahit = pos.eta();
int iphihitm = (hhit->id()).iphi();
int ietahitm = (hhit->id()).ieta();
int depthhit = (hhit->id()).depth();
float enehit = hhit->energy() * recal;
if (depthhit!=1) continue;
/*
float dphi = fabs(phihcal - phihit);
if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
float deta = fabs(etahcal - etahit);
float dr = sqrt(dphi*dphi + deta*deta);
*/
//double distAtHcal = getDistInPlaneTrackDir(gPointHcal, trackMomAtHcal, pos);
double distAtHcal = getDistInPlaneSimple(gPointHcal, pos);
if(distAtHcal < associationConeSize_)
{
for (HBHERecHitCollection::const_iterator hhit2=Hithbhe.begin(); hhit2!=Hithbhe.end(); hhit2++)
{
int iphihitm2 = (hhit2->id()).iphi();
int ietahitm2 = (hhit2->id()).ieta();
int depthhit2 = (hhit2->id()).depth();
float enehit2 = hhit2->energy() * recal;
if (iphihitm==iphihitm2 && ietahitm==ietahitm2 && depthhit!=depthhit2) enehit = enehit+enehit2;
}
//cout<<"IN CONE ieta: "<<ietahitm<<" iphi: "<<iphihitm<<" depthhit: "<<depthhit<<" dr: "<<dr<<" energy: "<<enehit<<endl;
//Find a Hit with Maximum Energy
if(enehit > MaxHit.hitenergy)
{
MaxHit.hitenergy = enehit;
MaxHit.ietahitm = (hhit->id()).ieta();
MaxHit.iphihitm = (hhit->id()).iphi();
MaxHit.dr = distAtHcal;
//MaxHit.depthhit = (hhit->id()).depth();
MaxHit.depthhit = 1;
//gPhot = geo->getPosition(hhit->detid());
}
}
} //end of all HBHE hits cycle
usedHits.clear();
//cout<<"Hottest ieta: "<<MaxHit.ietahitm<<" iphi: "<<MaxHit.iphihitm<<" dr: "<<MaxHit.dr<<endl;
//cout<<"Track ieta: "<<ietatrue<<" iphi: "<<iphitrue<<endl;
//cout<<"Point 3. MaxHit :::En "<<MaxHit.hitenergy<<" ieta: "<<MaxHit.ietahitm<<" iphi: "<<MaxHit.iphihitm<<endl;
Bool_t passCuts = kFALSE;
if(trackE > energyMinIso && trackE < energyMaxIso && emEnergy < energyECALmip && MaxHit.hitenergy > 0. && abs(MaxHit.ietahitm)<29) passCuts = kTRUE;
//cout<<"Pont 0.1.1. trackE:"<<trackE<<" emEn: "<<emEnergy<<endl;
numHits=0;
eClustBefore = 0.0;
eClustAfter = 0.0;
eBeforeDepth1 = 0.0;
eAfterDepth1 = 0.0;
eBeforeDepth2 = 0.0;
eAfterDepth2 = 0.0;
CentHitFactor = 0.0;
e3x3After = 0.0;
e3x3Before = 0.0;
e5x5After = 0.0;
e5x5Before = 0.0;
for (HBHERecHitCollection::const_iterator hhit=Hithbhe.begin(); hhit!=Hithbhe.end(); hhit++)
{
//check that this hit was not considered before and push it into usedHits
bool hitIsUsed=false;
int hitHashedIndex=hhit->id().hashed_index();
for (uint32_t i=0; i<usedHits.size(); i++)
{
if (usedHits[i]==hitHashedIndex) hitIsUsed=true;
}
if (hitIsUsed) continue;
usedHits.push_back(hitHashedIndex);
int DIETA = 100;
if(MaxHit.ietahitm*(hhit->id()).ieta()>0)
{
DIETA = MaxHit.ietahitm - (hhit->id()).ieta();
}
if(MaxHit.ietahitm*(hhit->id()).ieta()<0)
{
DIETA = MaxHit.ietahitm - (hhit->id()).ieta();
DIETA = DIETA>0 ? DIETA-1 : DIETA+1;
}
int DIPHI = abs(MaxHit.iphihitm - (hhit->id()).iphi());
DIPHI = DIPHI>36 ? 72-DIPHI : DIPHI;
int numbercell=100; //always collect Wide clastor!
//if(AxB_=="5x5" || AxB_=="3x3" || AxB_=="7x7"|| AxB_=="Cone")
//if(AxB_=="3x3") numbercell = 1;
//if(AxB_=="5x5") numbercell = 2;
//if(AxB_=="Cone") numbercell = 1000;
if( abs(DIETA)<=numbercell && (abs(DIPHI)<=numbercell || ( abs(MaxHit.ietahitm)>=20 && abs(DIPHI)<=numbercell+1)) )
{
const GlobalPoint pos2 = geo->getPosition(hhit->detid());
if(passCuts && hhit->energy()>0)
{
float factor = 0.0;
// factor = CalibFactors[hhit->id()];
factor = respRecalib -> getValues(hhit->id())->getValue();
//if(i<5){cout<<" calib factors: "<<factor<<" ij "<<100*i+j<<endl;}
if (hhit->id().ieta() == MaxHit.ietahitm && hhit->id().iphi() == MaxHit.iphihitm) CentHitFactor=factor;
if (hhit->id().ieta() == MaxHit.ietahitm && hhit->id().iphi() == MaxHit.iphihitm) iTime = hhit->time();
if(AxB_!="3x3" && AxB_!="5x5" && AxB_!="Cone") LogWarning(" AxB ")<<" Not supported: "<< AxB_;
if (abs(DIETA)<=2 && (abs(DIPHI)<=2 || ((abs(MaxHit.ietahitm)>20 && abs(DIPHI)<=4) && !( (abs(MaxHit.ietahitm)==21 || abs(MaxHit.ietahitm)==22) && abs((hhit->id()).ieta())<=20 && abs(DIPHI)>2) )) )
{
e5x5Before += hhit->energy();
e5x5After += hhit->energy()*factor;
}
if (abs(DIETA)<=1 && (abs(DIPHI)<=1 || ((abs(MaxHit.ietahitm)>20 && abs(DIPHI)<=2) && !(abs(MaxHit.ietahitm)==21 && abs((hhit->id()).ieta())<=20 && abs(DIPHI)>1) )) )
{
e3x3Before += hhit->energy();
e3x3After += hhit->energy()*factor;
}
if(AxB_=="5x5")
{
if (abs(DIETA)<=2 && (abs(DIPHI)<=2 || ( abs(MaxHit.ietahitm)>20 && abs(DIPHI)<=4) ) )
{
if (abs(MaxHit.ietahitm)==21 && abs((hhit->id()).ieta())<=20 && abs(DIPHI)>3) continue;
HTime[numHits]= hhit->time();
numHits++;
eClustBefore += hhit->energy();
eClustAfter += hhit->energy()*factor;
if ((hhit->id().depth() == 1) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29))
{
eBeforeDepth1 += hhit->energy();
eAfterDepth1 += hhit->energy()*factor;
}
else if((hhit->id().depth() == 2) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29))
{
eBeforeDepth2 += hhit->energy();
eAfterDepth2 += hhit->energy()*factor;
}
}
}//end of 5x5
if(AxB_=="3x3")
{
if (abs(DIETA)<=1 && (abs(DIPHI)<=1 || ( abs(MaxHit.ietahitm)>20 && abs(DIPHI)<=2) ) )
{
if (abs(MaxHit.ietahitm)==21 && abs((hhit->id()).ieta())<=20 && abs(DIPHI)>2) continue;
HTime[numHits]= hhit->time();
numHits++;
eClustBefore += hhit->energy();
eClustAfter += hhit->energy() * factor;
if ((hhit->id().depth() == 1) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29))
{
eBeforeDepth1 += hhit->energy();
eAfterDepth1 += hhit->energy() * factor;
}
else if((hhit->id().depth() == 2) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29))
{
eBeforeDepth2 += hhit->energy();
eAfterDepth2 += hhit->energy() * factor;
}
}
}//end of 3x3
if (AxB_=="Cone" && getDistInPlaneSimple(gPointHcal, pos2) < calibrationConeSize_) {
HTime[numHits]= hhit->time();
numHits++;
eClustBefore += hhit->energy();
eClustAfter += hhit->energy() * factor;
if ((hhit->id().depth() == 1) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29))
{
eBeforeDepth1 += hhit->energy();
eAfterDepth1 += hhit->energy() * factor;
}
else if((hhit->id().depth() == 2) && (abs(hhit->id().ieta()) > 17) && (abs(hhit->id().ieta()) < 29))
{
eBeforeDepth2 += hhit->energy();
eAfterDepth2 += hhit->energy() * factor;
}
}//end of Cone
}//end of passCuts
}//end of DIETA DIPHI
} //end of associatedcone HBHE hits cycle
int dieta_M_P = 100;
int diphi_M_P = 100;
if(MaxHit.ietahitm*ietatrue>0) {dieta_M_P = abs (MaxHit.ietahitm-ietatrue);}
if(MaxHit.ietahitm*ietatrue<0) {dieta_M_P = abs(MaxHit.ietahitm-ietatrue)-1;}
diphi_M_P = abs(MaxHit.iphihitm-iphitrue);
diphi_M_P = diphi_M_P>36 ? 72-diphi_M_P : diphi_M_P;
if(passCuts)
{
eventNumber = iEvent.id().event();
runNumber = iEvent.id().run();
eCentHitBefore = MaxHit.hitenergy;
eCentHitAfter = MaxHit.hitenergy*CentHitFactor;
eECAL = emEnergy;
numValidTrkHits = numVH;
numValidTrkStrips = numVS;
PtNearBy = ptNear;
eTrack = trackE;
phiTrack = trackPhi;
etaTrack = trackEta;
iEta = MaxHit.ietahitm;
iPhi = MaxHit.iphihitm;
iEtaTr = ietatrue;
iPhiTr = iphitrue;
iDr = sqrt(diphi_M_P*diphi_M_P+dieta_M_P*dieta_M_P);
delR = MaxHit.dr;
dietatr = dieta_M_P;
diphitr = diphi_M_P;
fTree -> Fill();
}
} //end of isoProdTracks cycle
/* ------------------ Some stuff for general tracks ---------- ----*/
//cout<<" generalTracks Size: "<< generalTracks->size()<<endl;
int n = generalTracks->size();
nTracks->Fill(n);
if(takeGenTracks_ && iEvent.id().event()%10==1)
{
gen = generalTracks->size();
iso = isoProdTracks->size();
pix = isoPixelTracks->size();
genPt[0] = -33;
genPhi[0] = -33;
genEta[0] = -33;
isoPt[0] = -33;
isoPhi[0] = -33;
isoEta[0] = -33;
pixPt[0] = -33;
pixPhi[0] = -33;
pixEta[0] = -33;
Int_t gencount=0, isocount=0, pixcount=0;
for (reco::TrackCollection::const_iterator gentr=generalTracks->begin(); gentr!=generalTracks->end(); gentr++)
{
genPt[gencount] = gentr->pt();
genPhi[gencount] = gentr->phi();
genEta[gencount] = gentr->eta();
gencount++;
}
for (reco::TrackCollection::const_iterator isotr=isoProdTracks->begin(); isotr!=isoProdTracks->end(); isotr++)
{
isoPt[isocount] = isotr->pt();
isoPhi[isocount] = isotr->phi();
isoEta[isocount] = isotr->eta();
isocount++;
}
for (reco::IsolatedPixelTrackCandidateCollection::const_iterator pixtr=isoPixelTracks->begin(); pixtr!=isoPixelTracks->end(); pixtr++)
{
pixPt[pixcount] = pixtr->pt();
pixPhi[pixcount] = pixtr->phi();
pixEta[pixcount] = pixtr->eta();
pixcount++;
}
}
tTree -> Fill();
}
| void ValidIsoTrkCalib::beginJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 776 of file ValidIsoTrkCalib.cc.
References cmsDownloadME::gen, and convertSQLiteXML::runNumber.
{
// if(!ReadCalibFactors(calibFactorsFileName_.c_str() )) {cout<<"Cant read file with cailib coefficients!! ---"<<endl;}
// try{
// edm::ESHandle <HcalRespCorrs> recalibCorrs;
// iSetup.get<HcalRespCorrsRcd>().get("recalibrate",recalibCorrs);
// respRecalib = recalibCorrs.product();
//
// LogInfo("CalibConstants")<<" Loaded: OK ";
//
// }catch(const cms::Exception & e) {
// LogWarning("CalibConstants")<<" Not Found!! ";
// }
// rootFile = new TFile(outputFileName_.c_str(),"RECREATE");
//@@@@@@@@@@@@@
//TFileDirectory ValDir = fs->mkdir("Validation");
nTracks = fs->make<TH1F>("nTracks","general;number of general tracks",11,-0.5,10.5);
tTree = fs->make<TTree>("tTree", "Tree for gen info");
fTree = fs->make<TTree>("fTree", "Tree for IsoTrack Calibration");
fTree->Branch("eventNumber", &eventNumber, "eventNumber/I");
fTree->Branch("runNumber", &runNumber, "runNumber/I");
fTree->Branch("eClustBefore", &eClustBefore,"eClustBefore/F");
fTree->Branch("eClustAfter", &eClustAfter,"eClustAfter/F");
fTree->Branch("eTrack", &eTrack, "eTrack/F");
fTree->Branch("etaTrack", &etaTrack, "etaTrack/F");
fTree->Branch("phiTrack", &phiTrack, "phiTrack/F");
fTree->Branch("numHits", &numHits, "numHits/I");
fTree->Branch("eECAL", &eECAL, "eECAL/F");
fTree->Branch("PtNearBy", &PtNearBy, "PtNearBy/F");
fTree->Branch("numValidTrkHits", &numValidTrkHits, "numValidTrkHits/F");
fTree->Branch("numValidTrkStrips", &numValidTrkStrips, "numValidTrkStrips/F");
fTree->Branch("eBeforeDepth1", &eBeforeDepth1, "eBeforeDepth1/F");
fTree->Branch("eBeforeDepth2", &eBeforeDepth2, "eBeforeDepth2/F");
fTree->Branch("eAfterDepth1", &eAfterDepth1, "eAfterDepth1/F");
fTree->Branch("eAfterDepth2", &eAfterDepth2, "eAfterDepth2/F");
fTree->Branch("e3x3Before", &e3x3Before, "e3x3Before/F");
fTree->Branch("e3x3After", &e3x3After, "e3x3After/F");
fTree->Branch("e5x5Before", &e5x5Before, "e5x5Before/F");
fTree->Branch("e5x5After", &e5x5After, "e5x5After/F");
fTree->Branch("eCentHitAfter", &eCentHitAfter, "eCentHitAfter/F");
fTree->Branch("eCentHitBefore", &eCentHitBefore, "eCentHitBefore/F");
fTree->Branch("iEta", &iEta, "iEta/I");
fTree->Branch("iPhi", &iPhi, "iPhi/I");
fTree->Branch("iEtaTr", &iEtaTr, "iEtaTr/I");
fTree->Branch("iPhiTr", &iPhiTr, "iPhiTr/I");
fTree->Branch("dietatr", &dietatr, "dietatr/I");
fTree->Branch("diphitr", &diphitr, "diphitr/I");
fTree->Branch("iDr", &iDr, "iDr/F");
fTree->Branch("delR", &delR, "delR/F");
fTree->Branch("iTime", &iTime, "iTime/F");
fTree->Branch("HTime", HTime, "HTime[numHits]/F");
fTree->Branch("xTrkEcal", &xTrkEcal, "xTrkEcal/F");
fTree->Branch("yTrkEcal", &yTrkEcal, "yTrkEcal/F");
fTree->Branch("zTrkEcal", &zTrkEcal, "zTrkEcal/F");
fTree->Branch("xTrkHcal", &xTrkHcal, "xTrkHcal/F");
fTree->Branch("yTrkHcal", &yTrkHcal, "yTrkHcal/F");
fTree->Branch("zTrkHcal", &zTrkHcal, "zTrkHcal/F");
if(takeGenTracks_) {
tTree->Branch("gen", &gen, "gen/I");
tTree->Branch("iso", &iso, "iso/I");
tTree->Branch("pix", &pix, "pix/I");
tTree->Branch("genPt", genPt, "genPt[gen]/F");
tTree->Branch("genPhi", genPhi, "genPhi[gen]/F");
tTree->Branch("genEta", genEta, "genEta[gen]/F");
tTree->Branch("isoPt", isoPt, "isoPt[iso]/F");
tTree->Branch("isoPhi", isoPhi, "isoPhi[iso]/F");
tTree->Branch("isoEta", isoEta, "isoEta[iso]/F");
tTree->Branch("pixPt", pixPt, "pixPt[pix]/F");
tTree->Branch("pixPhi", pixPhi, "pixPhi[pix]/F");
tTree->Branch("pixEta", pixEta, "pixEta[pix]/F");
}
}
| void ValidIsoTrkCalib::endJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 877 of file ValidIsoTrkCalib.cc.
{
// rootFile->Write();
//rootFile->Close();
}
bool ValidIsoTrkCalib::allowMissingInputs_ [private] |
Definition at line 108 of file ValidIsoTrkCalib.cc.
double ValidIsoTrkCalib::associationConeSize_ [private] |
Definition at line 104 of file ValidIsoTrkCalib.cc.
string ValidIsoTrkCalib::AxB_ [private] |
Definition at line 105 of file ValidIsoTrkCalib.cc.
double ValidIsoTrkCalib::calibrationConeSize_ [private] |
Definition at line 106 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::CentHitFactor [private] |
Definition at line 165 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::delR [private] |
Definition at line 170 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::dietatr [private] |
Definition at line 171 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::diphitr [private] |
Definition at line 172 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::e3x3After [private] |
Definition at line 177 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::e3x3Before [private] |
Definition at line 176 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::e5x5After [private] |
Definition at line 179 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::e5x5Before [private] |
Definition at line 178 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::eAfterDepth1 [private] |
Definition at line 160 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::eAfterDepth2 [private] |
Definition at line 162 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::eBeforeDepth1 [private] |
Definition at line 159 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::eBeforeDepth2 [private] |
Definition at line 161 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::eCentHitAfter [private] |
Definition at line 164 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::eCentHitBefore [private] |
Definition at line 163 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::eClustAfter [private] |
Definition at line 151 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::eClustBefore [private] |
Definition at line 150 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::eECAL [private] |
Definition at line 155 of file ValidIsoTrkCalib.cc.
Float_t ValidIsoTrkCalib::emEnergy [private] |
Definition at line 134 of file ValidIsoTrkCalib.cc.
double ValidIsoTrkCalib::energyECALmip [private] |
Definition at line 130 of file ValidIsoTrkCalib.cc.
double ValidIsoTrkCalib::energyMaxIso [private] |
Definition at line 131 of file ValidIsoTrkCalib.cc.
double ValidIsoTrkCalib::energyMinIso [private] |
Definition at line 131 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::etaTrack [private] |
Definition at line 153 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::eTrack [private] |
Definition at line 152 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::eventNumber [private] |
Definition at line 180 of file ValidIsoTrkCalib.cc.
edm::Service<TFileService> ValidIsoTrkCalib::fs [private] |
Definition at line 191 of file ValidIsoTrkCalib.cc.
TTree * ValidIsoTrkCalib::fTree [private] |
Definition at line 139 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::gen [private] |
Definition at line 115 of file ValidIsoTrkCalib.cc.
std::vector<edm::InputTag> ValidIsoTrkCalib::genecalLabel_ [private] |
Definition at line 95 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::genEta[500] [private] |
Definition at line 116 of file ValidIsoTrkCalib.cc.
InputTag ValidIsoTrkCalib::genhbheLabel_ [private] |
Definition at line 93 of file ValidIsoTrkCalib.cc.
InputTag ValidIsoTrkCalib::genhoLabel_ [private] |
Definition at line 94 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::genPhi[500] [private] |
Definition at line 116 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::genPt[500] [private] |
Definition at line 116 of file ValidIsoTrkCalib.cc.
InputTag ValidIsoTrkCalib::genTracksLabel_ [private] |
Definition at line 92 of file ValidIsoTrkCalib.cc.
const CaloGeometry* ValidIsoTrkCalib::geo [private] |
Definition at line 91 of file ValidIsoTrkCalib.cc.
InputTag ValidIsoTrkCalib::hbheLabel_ [private] |
Definition at line 96 of file ValidIsoTrkCalib.cc.
InputTag ValidIsoTrkCalib::hoLabel_ [private] |
Definition at line 97 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::HTime[100] [private] |
Definition at line 175 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::iDr [private] |
Definition at line 170 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::iEta [private] |
Definition at line 166 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::iEtaTr [private] |
Definition at line 168 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::iPhi [private] |
Definition at line 167 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::iPhiTr [private] |
Definition at line 169 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::iso [private] |
Definition at line 115 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::isoEta[500] [private] |
Definition at line 117 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::isoPhi[500] [private] |
Definition at line 117 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::isoPt[500] [private] |
Definition at line 117 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::iTime [private] |
Definition at line 174 of file ValidIsoTrkCalib.cc.
double ValidIsoTrkCalib::maxPNear [private] |
Definition at line 130 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::MinNTECHitsEndcap [private] |
Definition at line 129 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::MinNTrackHitsBarrel [private] |
Definition at line 129 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::Nhits [private] |
Definition at line 149 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::NisoTrk [private] |
Definition at line 123 of file ValidIsoTrkCalib.cc.
TH1F* ValidIsoTrkCalib::nTracks [private] |
Definition at line 189 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::numHits [private] |
Definition at line 156 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::numValidTrkHits [private] |
Definition at line 183 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::numValidTrkStrips [private] |
Definition at line 183 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::numVH [private] |
Definition at line 183 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::numVS [private] |
Definition at line 183 of file ValidIsoTrkCalib.cc.
string ValidIsoTrkCalib::outputFileName_ [private] |
Definition at line 109 of file ValidIsoTrkCalib.cc.
Definition at line 87 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::phiTrack [private] |
Definition at line 154 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::pix [private] |
Definition at line 115 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::pixEta[500] [private] |
Definition at line 118 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::pixPhi[500] [private] |
Definition at line 118 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::pixPt[500] [private] |
Definition at line 118 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::ptNear [private] |
Definition at line 125 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::PtNearBy [private] |
Definition at line 182 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::ptrack [private] |
Definition at line 126 of file ValidIsoTrkCalib.cc.
const HcalRespCorrs* ValidIsoTrkCalib::respRecalib [private] |
Definition at line 185 of file ValidIsoTrkCalib.cc.
int ValidIsoTrkCalib::runNumber [private] |
Definition at line 181 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::rvert [private] |
Definition at line 126 of file ValidIsoTrkCalib.cc.
double ValidIsoTrkCalib::taECALCone_ [private] |
Definition at line 88 of file ValidIsoTrkCalib.cc.
double ValidIsoTrkCalib::taHCALCone_ [private] |
Definition at line 89 of file ValidIsoTrkCalib.cc.
bool ValidIsoTrkCalib::takeGenTracks_ [private] |
Definition at line 113 of file ValidIsoTrkCalib.cc.
Float_t ValidIsoTrkCalib::targetE [private] |
Definition at line 135 of file ValidIsoTrkCalib.cc.
Definition at line 86 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::trackE [private] |
Definition at line 124 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::trackEta [private] |
Definition at line 124 of file ValidIsoTrkCalib.cc.
InputTag ValidIsoTrkCalib::trackLabel1_ [private] |
Definition at line 99 of file ValidIsoTrkCalib.cc.
InputTag ValidIsoTrkCalib::trackLabel_ [private] |
Definition at line 98 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::trackPhi [private] |
Definition at line 124 of file ValidIsoTrkCalib.cc.
float ValidIsoTrkCalib::trackPt [private] |
Definition at line 124 of file ValidIsoTrkCalib.cc.
TTree* ValidIsoTrkCalib::tTree [private] |
Definition at line 139 of file ValidIsoTrkCalib.cc.
Float_t ValidIsoTrkCalib::xTrkEcal [private] |
Definition at line 141 of file ValidIsoTrkCalib.cc.
Float_t ValidIsoTrkCalib::xTrkHcal [private] |
Definition at line 145 of file ValidIsoTrkCalib.cc.
Float_t ValidIsoTrkCalib::yTrkEcal [private] |
Definition at line 142 of file ValidIsoTrkCalib.cc.
Float_t ValidIsoTrkCalib::yTrkHcal [private] |
Definition at line 146 of file ValidIsoTrkCalib.cc.
Float_t ValidIsoTrkCalib::zTrkEcal [private] |
Definition at line 143 of file ValidIsoTrkCalib.cc.
Float_t ValidIsoTrkCalib::zTrkHcal [private] |
Definition at line 147 of file ValidIsoTrkCalib.cc.
1.7.3