|
|
|
#include <Calibration/HcalCalibAlgos/src/HcalIsoTrkAnalyzer.cc>
Description: see twiki for details: https://twiki.cern.ch/twiki/bin/view/CMS/IsoTrackAnalyzer
Implementation: <Notes on="" implementation>="">
Definition at line 98 of file HcalIsoTrkAnalyzer.cc.
| HcalIsoTrkAnalyzer::HcalIsoTrkAnalyzer | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 239 of file HcalIsoTrkAnalyzer.cc.
References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and Parameters::parameters.
{
m_ecalLabel = iConfig.getUntrackedParameter<std::string> ("ecalRecHitsLabel","ecalRecHit");
m_ebInstance = iConfig.getUntrackedParameter<std::string> ("ebRecHitsInstance","EcalRecHitsEB");
m_eeInstance = iConfig.getUntrackedParameter<std::string> ("eeRecHitsInstance","EcalRecHitsEE");
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");
AxB_=iConfig.getParameter<std::string>("AxB");
calibrationConeSize_=iConfig.getParameter<double>("calibrationConeSize");
nIterations = iConfig.getParameter<int>("noOfIterations");
eventWeight = iConfig.getParameter<double>("eventWeight");
energyMinIso = iConfig.getParameter<double>("energyMinIso");
energyMaxIso = iConfig.getParameter<double>("energyMaxIso");
energyECALmip = iConfig.getParameter<double>("energyECALmip");
maxPNear = iConfig.getParameter<double>("maxPNear");
MinNTrackHitsBarrel = iConfig.getParameter<int>("MinNTrackHitsBarrel");
MinNTECHitsEndcap = iConfig.getParameter<int>("MinNTECHitsEndcap");
edm::ParameterSet parameters = iConfig.getParameter<edm::ParameterSet>("TrackAssociatorParameters");
parameters_.loadParameters( parameters );
trackAssociator_.useDefaultPropagator();
}
| HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer | ( | ) |
Definition at line 274 of file HcalIsoTrkAnalyzer.cc.
{
}
| void HcalIsoTrkAnalyzer::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 280 of file HcalIsoTrkAnalyzer.cc.
References abs, edm::SortedCollection< T, SORT >::begin(), MaxHit_struct::depthhit, patCandidatesForDimuonsSequences_cff::ecal, EcalBarrel, EcalEndcap, edm::SortedCollection< T, SORT >::end(), PV3DBase< T, PVType, FrameType >::eta(), edm::EventID::event(), GeneralTracks_cfi::generalTracks, edm::EventSetup::get(), edm::Event::getByLabel(), CaloSubdetectorGeometry::getClosestCell(), getDistInPlaneSimple(), EEDetId::hashedIndex(), EBDetId::hashedIndex(), DetId::Hcal, HcalBarrel, HcalEndcap, MaxHit_struct::hitenergy, i, edm::EventBase::id(), HcalDetId::ieta(), MaxHit_struct::ietahitm, info, HcalDetId::iphi(), MaxHit_struct::iphihitm, PV3DBase< T, PVType, FrameType >::phi(), pos, MaxHit_struct::posMax, edm::Handle< T >::product(), edm::SortedCollection< T, SORT >::push_back(), edm::EventID::run(), inputsource_file_cfi::runNumber, mathSSE::sqrt(), reco::btau::trackEta, reco::btau::trackPhi, and diffTreeTool::tree.
{
using namespace edm;
using namespace std;
vector<float> rawEnergyVec;
vector<int> detiphi;
vector<int> detieta;
vector<int> i3i5;
vector<HcalDetId> detidvec;
float calEnergy;
edm::Handle<reco::TrackCollection> generalTracks;
iEvent.getByLabel(trackLabel1_,generalTracks);
edm::Handle<reco::IsolatedPixelTrackCandidateCollection> trackCollection;
iEvent.getByLabel(trackLabel_,trackCollection);
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();
// rof 16.05.2008 start: include the possibility for recalibration (use "recalibrate" label for safety)
/*
edm::ESHandle <HcalRespCorrs> recalibCorrs;
iSetup.get<HcalRespCorrsRcd>().get("recalibrate",recalibCorrs);
const HcalRespCorrs* myRecalib = recalibCorrs.product();
*/
// rof end
parameters_.useEcal = true;
parameters_.useHcal = true;
parameters_.useCalo = false;
parameters_.useMuon = false;
parameters_.dREcal = 0.5;
parameters_.dRHcal = 0.6;
if (trackCollection->size()==0) return;
for (reco::TrackCollection::const_iterator trit=generalTracks->begin(); trit!=generalTracks->end(); trit++)
{
reco::IsolatedPixelTrackCandidateCollection::const_iterator isoMatched=trackCollection->begin();
bool matched=false;
for (reco::IsolatedPixelTrackCandidateCollection::const_iterator it = trackCollection->begin(); it!=trackCollection->end(); it++)
{
if (abs((trit->pt() - it->pt())/it->pt()) < 0.005 && abs(trit->eta() - it->eta()) < 0.01)
{
isoMatched=it;
matched=true;
break;
}
}
if (!matched) continue;
if (trit->hitPattern().numberOfValidHits()<MinNTrackHitsBarrel) continue;
if (fabs(trit->eta())>1.47&&trit->hitPattern().numberOfValidStripTECHits()<MinNTECHitsEndcap) continue;
//container for used recHits
std::vector<int> usedHits;
//
calEnergy = sqrt(trit->px()*trit->px()+trit->py()*trit->py()+trit->pz()*trit->pz()+0.14*0.14);
trackEta = trit->eta();
trackPhi = trit->phi();
double corrHCAL = 1.; //another possibility for correction. - why?
// cout << endl << " ISO TRACK E = "<< calEnergy << " ETA = " << trackEta<< " PHI = " << trackPhi << " Correction " << corrHCAL<< endl;
rvert = sqrt(trit->vx()*trit->vx()+trit->vy()*trit->vy()+trit->vz()*trit->vz());
//Associate track with a calorimeter
TrackDetMatchInfo info = trackAssociator_.associate(iEvent, iSetup,trackAssociator_.getFreeTrajectoryState(iSetup, *trit),parameters_);
//*(it->track().get())
double etaecal=info.trkGlobPosAtEcal.eta();
double phiecal=info.trkGlobPosAtEcal.phi();
double etahcal=info.trkGlobPosAtHcal.eta();
//double 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();
GlobalPoint gP(xTrkHcal,yTrkHcal,zTrkHcal);
int iphitrue = -10;
int ietatrue = 100;
if (etahcal<1.392)
{
const CaloSubdetectorGeometry* gHB = geo->getSubdetectorGeometry(DetId::Hcal,HcalBarrel);
// const GlobalPoint tempPoint(newx, newy, newz);
//const DetId tempId = gHB->getClosestCell(tempPoint);
const HcalDetId tempId = gHB->getClosestCell(gP);
ietatrue = tempId.ieta();
iphitrue = tempId.iphi();
}
if (etahcal>1.392 && etahcal<3.0)
{
const CaloSubdetectorGeometry* gHE = geo->getSubdetectorGeometry(DetId::Hcal,HcalEndcap);
const HcalDetId tempId = gHE->getClosestCell(gP);
ietatrue = tempId.ieta();
iphitrue = tempId.iphi();
}
// eecal=info.coneEnergy(parameters_.dREcal, TrackDetMatchInfo::EcalRecHits);
// ehcal=info.coneEnergy(parameters_.dRHcal, TrackDetMatchInfo::HcalRecHits);
double rmin = 0.07;
if( fabs(etaecal) > 1.47 ) rmin = 0.07*(fabs(etaecal)-0.47)*1.2;
if( fabs(etaecal) > 2.2 ) rmin = 0.07*(fabs(etaecal)-0.47)*1.4;
MaxHit_struct MaxHit;
// Find Ecal RecHit with maximum energy and collect other information
MaxHit.hitenergy=-100;
double econus = 0.;
float ecal_cluster = 0.;
//clear usedHits
usedHits.clear();
//
for (std::vector<EcalRecHit>::const_iterator ehit=Hitecal.begin(); ehit!=Hitecal.end(); ehit++)
{
//check that this hit was not considered before and push it into usedHits
bool hitIsUsed=false;
int hitHashedIndex=-10000;
if (ehit->id().subdetId()==EcalBarrel)
{
EBDetId did(ehit->id());
hitHashedIndex=did.hashedIndex();
}
if (ehit->id().subdetId()==EcalEndcap)
{
EEDetId did(ehit->id());
hitHashedIndex=did.hashedIndex();
}
for (uint32_t i=0; i<usedHits.size(); i++)
{
if (usedHits[i]==hitHashedIndex) hitIsUsed=true;
}
if (hitIsUsed) continue;
usedHits.push_back(hitHashedIndex);
//
if((*ehit).energy() > MaxHit.hitenergy)
{
MaxHit.hitenergy = (*ehit).energy();
}
GlobalPoint pos = geo->getPosition((*ehit).detid());
double phihit = pos.phi();
double etahit = pos.eta();
double dphi = fabs(phiecal - phihit);
if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
double deta = fabs(etaecal - etahit);
double dr = sqrt(dphi*dphi + deta*deta);
if (dr < rmin) {
econus = econus + (*ehit).energy();
}
if (dr < 0.13) ecal_cluster += (*ehit).energy();
}
MaxHit.hitenergy=-100;
//clear usedHits
usedHits.clear();
//
//float dddeta = 1000.;
//float dddphi = 1000.;
//int iphitrue = 1234;
//int ietatrue = 1234;
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());
double phihit = pos.phi();
double etahit = pos.eta();
int iphihitm = (hhit->id()).iphi();
int ietahitm = (hhit->id()).ieta();
int depthhit = (hhit->id()).depth();
double enehit = hhit->energy() * recal;
if (depthhit!=1) continue;
double dphi = fabs(info.trkGlobPosAtHcal.phi() - phihit);
if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
double deta = fabs(info.trkGlobPosAtHcal.eta() - etahit);
double dr = sqrt(dphi*dphi + deta*deta);
/* Commented. Another way of finding true projection is implemented
if (deta<dddeta) {
ietatrue = ietahitm;
dddeta=deta;
}
if (dphi<dddphi) {
iphitrue = iphihitm;
dddphi=dphi;
}*/
if(dr<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();
double enehit2 = hhit2->energy() * recal;
if ( iphihitm==iphihitm2 && ietahitm==ietahitm2 && depthhit!=depthhit2){
enehit = enehit+enehit2;
}
}
if(enehit > MaxHit.hitenergy)
{
MaxHit.hitenergy = enehit;
MaxHit.ietahitm = (hhit->id()).ieta();
MaxHit.iphihitm = (hhit->id()).iphi();
MaxHit.depthhit = (hhit->id()).depth();
MaxHit.posMax = geo->getPosition(hhit->detid());
}
}
}
Bool_t passCuts = kFALSE;
if(calEnergy > energyMinIso && calEnergy < energyMaxIso && isoMatched->energyIn() < energyECALmip &&
isoMatched->maxPtPxl() < maxPNear && abs(MaxHit.ietahitm)<30 && MaxHit.hitenergy > 0.){ passCuts = kTRUE; }
if(AxB_=="5x5" || AxB_=="3x3" || AxB_=="7x7"|| AxB_=="Cone")
{
//clear usedHits
usedHits.clear();
//
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;
/*AP DIPHI = DIPHI<-36 ? 72+DIPHI : DIPHI; */
int numbercell=0;
if(AxB_=="3x3") numbercell = 1;
if(AxB_=="5x5") numbercell = 2;
if(AxB_=="7x7") numbercell = 3;
if(AxB_=="Cone") numbercell = 1000;
if( abs(DIETA)<=numbercell && (abs(DIPHI)<=numbercell || ( abs(MaxHit.ietahitm)>=20 && abs(DIPHI)<=numbercell+1)) ) {
// rof 16.05.2008 start: include the possibility for recalibration
float recal = 1;
int iii3i5 = 0;
const GlobalPoint pos2 = geo->getPosition(hhit->detid());
if(passCuts){
if(AxB_=="5x5" || AxB_=="3x3" || AxB_=="7x7") {
rawEnergyVec.push_back(hhit->energy() * recal * corrHCAL);
detidvec.push_back(hhit->id());
detiphi.push_back((hhit->id()).iphi());
detieta.push_back((hhit->id()).ieta());
i3i5.push_back(iii3i5);
}
if (AxB_=="Cone" && getDistInPlaneSimple(gP,pos2) < calibrationConeSize_) {
rawEnergyVec.push_back(hhit->energy() * recal * corrHCAL);
detidvec.push_back(hhit->id());
detiphi.push_back((hhit->id()).iphi());
detieta.push_back((hhit->id()).ieta());
i3i5.push_back(iii3i5);
}
}
}
}
}
if(AxB_!="3x3" && AxB_!="5x5" && AxB_!="7x7" && AxB_!="Cone") LogWarning(" AxB ")<<" Not supported: "<< AxB_;
if(passCuts){
input_to_L3 << rawEnergyVec.size() << " " << calEnergy;
for (unsigned int i=0; i<rawEnergyVec.size(); i++)
{
input_to_L3 << " " << rawEnergyVec.at(i) << " " << detidvec.at(i).rawId() ;
}
input_to_L3 <<endl;
eventNumber = iEvent.id().event();
runNumber = iEvent.id().run();
iEtaHit = ietatrue;
iPhiHit = iphitrue;
emEnergy = isoMatched->energyIn();
// exampleP4->SetPxPyPzE(2, 1, 1, 10);
numberOfCells=rawEnergyVec.size();
targetE = calEnergy;
for (unsigned int ia=0; ia<numberOfCells; ++ia) {
cellEnergy = rawEnergyVec.at(ia);
cell = detidvec.at(ia).rawId();
new((*cells)[ia]) TCell(cell, cellEnergy);
}
tree->Fill();
cells->Clear();
}
rawEnergyVec.clear();
detidvec.clear();
detiphi.clear();
detieta.clear();
i3i5.clear();
try {
Handle<HORecHitCollection> ho;
iEvent.getByLabel(hoLabel_,ho);
const HORecHitCollection Hitho = *(ho.product());
//clear usedHits
usedHits.clear();
//
for(HORecHitCollection::const_iterator hoItr=Hitho.begin(); hoItr!=Hitho.end(); hoItr++)
{
//check that this hit was not considered before and push it into usedHits
bool hitIsUsed=false;
int hitHashedIndex=hoItr->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);
/*AP
GlobalPoint pos = geo->getPosition(hoItr->detid());
double phihit = pos.phi();
double etahit = pos.eta();
int iphihitm = (hoItr->id()).iphi();
int ietahitm = (hoItr->id()).ieta();
int depthhit = (hoItr->id()).depth();
double dphi = fabs(trackPhi - phihit);
if(dphi > 4.*atan(1.)) dphi = 8.*atan(1.) - dphi;
double deta = fabs(trackEta - etahit);
double dr = sqrt(dphi*dphi + deta*deta);
*/
}
} catch (cms::Exception& e) { // can't find it!
if (!allowMissingInputs_) throw e;
}
}
}
| void HcalIsoTrkAnalyzer::beginJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 736 of file HcalIsoTrkAnalyzer.cc.
References BuildWebpage::rootFile, inputsource_file_cfi::runNumber, and diffTreeTool::tree.
{
// MyL3Algo = new MinL3AlgoUniv<HcalDetId>(eventWeight);
input_to_L3.open("input_to_l3.txt");
rootFile = new TFile("rootFile.root", "RECREATE");
tree = new TTree("hcalCalibTree", "Tree for IsoTrack Calibration");
cells = new TClonesArray("TCell", 10000);
// cells3x3 = new TRefArray;
// cellsPF = new TRefArray;
// exampleP4 = new TLorentzVector();
tagJetP4 = new TLorentzVector(); // dijet
probeJetP4 = new TLorentzVector(); // dijet
tree->Branch("cells", &cells, 64000, 0);
tree->Branch("targetE", &targetE, "targetE/F");
tree->Branch("emEnergy", &emEnergy, "emEnergy/F");
tree->Branch("xTrkEcal", &xTrkEcal, "xTrkEcal/F");
tree->Branch("yTrkEcal", &yTrkEcal, "yTrkEcal/F");
tree->Branch("zTrkEcal", &zTrkEcal, "zTrkEcal/F");
tree->Branch("xTrkHcal", &xTrkHcal, "xTrkHcal/F");
tree->Branch("yTrkHcal", &yTrkHcal, "yTrkHcal/F");
tree->Branch("zTrkHcal", &zTrkHcal, "zTrkHcal/F");
// tree->Branch("exampleP4", "TLorentzVector", &exampleP4);
// tree->Branch("cells3x3", &cells3x3, 64000, 0);
// tree->Branch("cellsPF", &cellsPF, 64000, 0);
tree->Branch("iEtaHit", &iEtaHit, "iEtaHit/I");
tree->Branch("iPhiHit", &iPhiHit, "iPhiHit/i");
tree->Branch("eventNumber", &eventNumber, "eventNumber/i");
tree->Branch("runNumber", &runNumber, "runNumber/i");
tree->Branch("tagJetP4", "TLorentzVector", &tagJetP4); // dijet
tree->Branch("probeJetP4", "TLorentzVector", &probeJetP4); // dijet
tree->Branch("etVetoJet", &etVetoJet, "etVetoJet/F"); // dijet
tree->Branch("tagJetEmFrac", &tagJetEmFrac,"tagJetEmFrac/F"); // dijet
tree->Branch("probeJetEmFrac", &probeJetEmFrac,"probeJetEmFrac/F"); // dijet
}
| void HcalIsoTrkAnalyzer::endJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 781 of file HcalIsoTrkAnalyzer.cc.
References BuildWebpage::rootFile.
{
input_to_L3.close();
rootFile->Write();
rootFile->Close();
if (cells) delete cells;
// if (cells3x3) delete cells3x3;
// if (cellsPF) delete cellsPF;
// if (exampleP4) delete exampleP4;
if (tagJetP4) delete tagJetP4; // dijet
if (probeJetP4) delete probeJetP4; // dijet
}
| double HcalIsoTrkAnalyzer::getDistInPlaneSimple | ( | const GlobalPoint | caloPoint, |
| const GlobalPoint | rechitPoint | ||
| ) |
Definition at line 194 of file HcalIsoTrkAnalyzer.cc.
References Vector3DBase< T, FrameTag >::dot(), PV3DBase< T, PVType, FrameType >::mag(), Vector3DBase< T, FrameTag >::unit(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
// Simplified version of getDistInPlane
// Assume track direction is origin -> point of hcal intersection
const GlobalVector caloIntersectVector(caloPoint.x(),
caloPoint.y(),
caloPoint.z());
const GlobalVector caloIntersectUnitVector = caloIntersectVector.unit();
const GlobalVector rechitVector(rechitPoint.x(),
rechitPoint.y(),
rechitPoint.z());
const GlobalVector rechitUnitVector = rechitVector.unit();
double dotprod = caloIntersectUnitVector.dot(rechitUnitVector);
double rechitdist = caloIntersectVector.mag()/dotprod;
const GlobalVector effectiveRechitVector = rechitdist*rechitUnitVector;
const GlobalPoint effectiveRechitPoint(effectiveRechitVector.x(),
effectiveRechitVector.y(),
effectiveRechitVector.z());
GlobalVector distance_vector = effectiveRechitPoint-caloPoint;
if (dotprod > 0.)
{
return distance_vector.mag();
}
else
{
return 999999.;
}
}
bool HcalIsoTrkAnalyzer::allowMissingInputs_ [private] |
Definition at line 132 of file HcalIsoTrkAnalyzer.cc.
double HcalIsoTrkAnalyzer::associationConeSize_ [private] |
Definition at line 130 of file HcalIsoTrkAnalyzer.cc.
string HcalIsoTrkAnalyzer::AxB_ [private] |
Definition at line 131 of file HcalIsoTrkAnalyzer.cc.
double HcalIsoTrkAnalyzer::calibrationConeSize_ [private] |
Definition at line 130 of file HcalIsoTrkAnalyzer.cc.
UInt_t HcalIsoTrkAnalyzer::cell [private] |
Definition at line 163 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::cellEnergy [private] |
Definition at line 164 of file HcalIsoTrkAnalyzer.cc.
TClonesArray* HcalIsoTrkAnalyzer::cells [private] |
Definition at line 166 of file HcalIsoTrkAnalyzer.cc.
TRefArray* HcalIsoTrkAnalyzer::cells3x3 [private] |
Definition at line 167 of file HcalIsoTrkAnalyzer.cc.
TRefArray* HcalIsoTrkAnalyzer::cellsPF [private] |
Definition at line 168 of file HcalIsoTrkAnalyzer.cc.
char HcalIsoTrkAnalyzer::dirname[50] [private] |
Definition at line 155 of file HcalIsoTrkAnalyzer.cc.
InputTag HcalIsoTrkAnalyzer::eLabel_ [private] |
Definition at line 119 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::emEnergy [private] |
Definition at line 182 of file HcalIsoTrkAnalyzer.cc.
double HcalIsoTrkAnalyzer::energyECALmip [private] |
Definition at line 153 of file HcalIsoTrkAnalyzer.cc.
double HcalIsoTrkAnalyzer::energyMaxIso [private] |
Definition at line 152 of file HcalIsoTrkAnalyzer.cc.
double HcalIsoTrkAnalyzer::energyMinIso [private] |
Definition at line 152 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::etVetoJet [private] |
Definition at line 186 of file HcalIsoTrkAnalyzer.cc.
UInt_t HcalIsoTrkAnalyzer::eventNumber [private] |
Definition at line 169 of file HcalIsoTrkAnalyzer.cc.
float HcalIsoTrkAnalyzer::eventWeight [private] |
Definition at line 151 of file HcalIsoTrkAnalyzer.cc.
const CaloGeometry* HcalIsoTrkAnalyzer::geo [private] |
Definition at line 116 of file HcalIsoTrkAnalyzer.cc.
InputTag HcalIsoTrkAnalyzer::hbheLabel_ [private] |
Definition at line 117 of file HcalIsoTrkAnalyzer.cc.
char HcalIsoTrkAnalyzer::hname[20] [private] |
Definition at line 156 of file HcalIsoTrkAnalyzer.cc.
InputTag HcalIsoTrkAnalyzer::hoLabel_ [private] |
Definition at line 118 of file HcalIsoTrkAnalyzer.cc.
char HcalIsoTrkAnalyzer::htitle[80] [private] |
Definition at line 157 of file HcalIsoTrkAnalyzer.cc.
Int_t HcalIsoTrkAnalyzer::iEtaHit [private] |
Definition at line 171 of file HcalIsoTrkAnalyzer.cc.
ofstream HcalIsoTrkAnalyzer::input_to_L3 [private] |
Definition at line 190 of file HcalIsoTrkAnalyzer.cc.
UInt_t HcalIsoTrkAnalyzer::iPhiHit [private] |
Definition at line 172 of file HcalIsoTrkAnalyzer.cc.
std::string HcalIsoTrkAnalyzer::m_ebInstance [private] |
Definition at line 125 of file HcalIsoTrkAnalyzer.cc.
std::string HcalIsoTrkAnalyzer::m_ecalLabel [private] |
Definition at line 124 of file HcalIsoTrkAnalyzer.cc.
std::string HcalIsoTrkAnalyzer::m_eeInstance [private] |
Definition at line 126 of file HcalIsoTrkAnalyzer.cc.
std::string HcalIsoTrkAnalyzer::m_hcalLabel [private] |
Definition at line 127 of file HcalIsoTrkAnalyzer.cc.
int HcalIsoTrkAnalyzer::m_histoFlag [private] |
Definition at line 128 of file HcalIsoTrkAnalyzer.cc.
std::string HcalIsoTrkAnalyzer::m_inputTrackLabel [private] |
Definition at line 123 of file HcalIsoTrkAnalyzer.cc.
double HcalIsoTrkAnalyzer::maxPNear [private] |
Definition at line 153 of file HcalIsoTrkAnalyzer.cc.
int HcalIsoTrkAnalyzer::MinNTECHitsEndcap [private] |
Definition at line 150 of file HcalIsoTrkAnalyzer.cc.
int HcalIsoTrkAnalyzer::MinNTrackHitsBarrel [private] |
Definition at line 150 of file HcalIsoTrkAnalyzer.cc.
int HcalIsoTrkAnalyzer::nIterations [private] |
Definition at line 150 of file HcalIsoTrkAnalyzer.cc.
UInt_t HcalIsoTrkAnalyzer::numberOfCells [private] |
Definition at line 162 of file HcalIsoTrkAnalyzer.cc.
string HcalIsoTrkAnalyzer::outputFileName_ [private] |
Definition at line 133 of file HcalIsoTrkAnalyzer.cc.
Definition at line 114 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::probeJetEmFrac [private] |
Definition at line 188 of file HcalIsoTrkAnalyzer.cc.
TLorentzVector* HcalIsoTrkAnalyzer::probeJetP4 [private] |
Definition at line 185 of file HcalIsoTrkAnalyzer.cc.
TFile* HcalIsoTrkAnalyzer::rootFile [private] |
Definition at line 159 of file HcalIsoTrkAnalyzer.cc.
UInt_t HcalIsoTrkAnalyzer::runNumber [private] |
Definition at line 170 of file HcalIsoTrkAnalyzer.cc.
double HcalIsoTrkAnalyzer::rvert [private] |
Definition at line 136 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::tagJetEmFrac [private] |
Definition at line 187 of file HcalIsoTrkAnalyzer.cc.
TLorentzVector* HcalIsoTrkAnalyzer::tagJetP4 [private] |
Definition at line 184 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::targetE [private] |
Definition at line 161 of file HcalIsoTrkAnalyzer.cc.
Definition at line 113 of file HcalIsoTrkAnalyzer.cc.
double HcalIsoTrkAnalyzer::trackEta [private] |
Definition at line 135 of file HcalIsoTrkAnalyzer.cc.
InputTag HcalIsoTrkAnalyzer::trackLabel1_ [private] |
Definition at line 121 of file HcalIsoTrkAnalyzer.cc.
InputTag HcalIsoTrkAnalyzer::trackLabel_ [private] |
Definition at line 120 of file HcalIsoTrkAnalyzer.cc.
double HcalIsoTrkAnalyzer::trackPhi [private] |
Definition at line 135 of file HcalIsoTrkAnalyzer.cc.
TTree* HcalIsoTrkAnalyzer::tree [private] |
Definition at line 160 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::xTrkEcal [private] |
Definition at line 174 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::xTrkHcal [private] |
Definition at line 178 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::yTrkEcal [private] |
Definition at line 175 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::yTrkHcal [private] |
Definition at line 179 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::zTrkEcal [private] |
Definition at line 176 of file HcalIsoTrkAnalyzer.cc.
Float_t HcalIsoTrkAnalyzer::zTrkHcal [private] |
Definition at line 180 of file HcalIsoTrkAnalyzer.cc.
1.7.3