#include <Calibration/HcalCalibAlgos/src/HcalIsoTrkAnalyzer.cc>

Inheritance diagram for HcalIsoTrkAnalyzer:

Public Member Functions
	HcalIsoTrkAnalyzer (const edm::ParameterSet &)

	~HcalIsoTrkAnalyzer ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &) override

virtual void	beginJob () override

virtual void	endJob () override

Private Attributes
bool	allowMissingInputs_

double	associationConeSize_

string	AxB_

double	calibrationConeSize_

UInt_t	cell

Float_t	cellEnergy

TClonesArray *	cells

TRefArray *	cells3x3

TRefArray *	cellsPF

double	EcalCone

double	EcalConeOuter

Float_t	emEnergy

Float_t	emRingEnergy

double	energyECALmip

double	energyMaxIso

double	energyMinIso

Float_t	etVetoJet

UInt_t	eventNumber

float	eventWeight

const CaloGeometry *	geo

double	hottestHitDistance

Int_t	iEtaHit

std::ofstream	input_to_L3

UInt_t	iPhiHit

std::string	m_ebInstance

std::string	m_ecalLabel

std::string	m_eeInstance

std::string	m_hcalLabel

int	m_histoFlag

std::string	m_inputTrackLabel

double	maxPNear

int	MinNTECHitsEndcap

int	MinNTrackHitsBarrel

int	nIterations

UInt_t	numberOfCells

string	outputFileName_

TrackAssociatorParameters	parameters_

Float_t	probeJetEmFrac

TLorentzVector *	probeJetP4

Float_t	PxTrkHcal

Float_t	PyTrkHcal

Float_t	PzTrkHcal

TFile *	rootFile

UInt_t	runNumber

double	rvert

Float_t	tagJetEmFrac

TLorentzVector *	tagJetP4

Float_t	targetE

edm::EDGetTokenT < EcalRecHitCollection >	tok_ecal_

edm::EDGetTokenT < HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT < HORecHitCollection >	tok_ho_

edm::EDGetTokenT < reco::TrackCollection >	tok_track1_

edm::EDGetTokenT < reco::IsolatedPixelTrackCandidateCollection >	tok_track_

TrackDetectorAssociator	trackAssociator_

double	trackEta

double	trackPhi

TTree *	tree

Float_t	xTrkEcal

Float_t	xTrkHcal

Float_t	yTrkEcal

Float_t	yTrkHcal

Float_t	zTrkEcal

Float_t	zTrkHcal

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: see twiki for details: https://twiki.cern.ch/twiki/bin/view/CMS/IsoTrackAnalyzer

Implementation: <Notes on="" implementation>="">

Definition at line 92 of file HcalIsoTrkAnalyzer.cc.

Constructor & Destructor Documentation

HcalIsoTrkAnalyzer::HcalIsoTrkAnalyzer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 172 of file HcalIsoTrkAnalyzer.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), Parameters::parameters, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   
   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");
  
   tok_hbhe_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInput"));
   tok_ho_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("hoInput"));
   tok_ecal_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("eInput"));
   tok_track_ = consumes<reco::IsolatedPixelTrackCandidateCollection>(iConfig.getParameter<edm::InputTag>("HcalIsolTrackInput"));
   tok_track1_ = consumes<reco::TrackCollection>(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");
   hottestHitDistance = iConfig.getParameter<double>("hottestHitDistance");
   EcalCone = iConfig.getParameter<double>("EcalCone");
   EcalConeOuter = iConfig.getParameter<double>("EcalConeOuter");
 
   edm::ParameterSet parameters = iConfig.getParameter<edm::ParameterSet>("TrackAssociatorParameters");
   edm::ConsumesCollector iC = consumesCollector();
   parameters_.loadParameters( parameters, iC );
   trackAssociator_.useDefaultPropagator();
 
 }

HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer ( )

Definition at line 211 of file HcalIsoTrkAnalyzer.cc.

212 {

213 }

Member Function Documentation

void HcalIsoTrkAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 217 of file HcalIsoTrkAnalyzer.cc.

References funct::abs(), edm::SortedCollection< T, SORT >::begin(), gather_cfg::cout, HLT_25ns14e33_v1_cff::depth, MaxHit_struct::depthhit, MaxHit_struct::dr, alignCSCRings::e, patCandidatesForDimuonsSequences_cff::ecal, ecalEnergyInCone(), edm::SortedCollection< T, SORT >::end(), edm::EventID::event(), edm::EventSetup::get(), edm::Event::getByToken(), getDistInPlaneTrackDir(), DetId::Hcal, HcalBarrel, MaxHit_struct::hitenergy, i, edm::EventBase::id(), HcalDetId::ieta(), MaxHit_struct::ietahitm, info(), HcalDetId::iphi(), MaxHit_struct::iphihitm, MaxHit_struct::posMax, edm::Handle< T >::product(), edm::EventID::run(), convertSQLiteXML::runNumber, mathSSE::sqrt(), reco::btau::trackEta, reco::btau::trackPhi, MainPageGenerator::tree, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

 {
   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> isoProdTracks;
   iEvent.getByToken(tok_track1_,isoProdTracks);  
 
   edm::Handle<reco::IsolatedPixelTrackCandidateCollection> isoPixelTracks;
   iEvent.getByToken(tok_track_,isoPixelTracks);
     
   edm::Handle<EcalRecHitCollection> ecal;
   iEvent.getByToken(tok_ecal_,ecal);
   const EcalRecHitCollection Hitecal = *(ecal.product());
     
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(tok_hbhe_,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!
 
 
 // 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 (isoPixelTracks->size()==0) return;
 
   for (reco::TrackCollection::const_iterator trit=isoProdTracks->begin(); trit!=isoProdTracks->end(); trit++)
     {
       reco::IsolatedPixelTrackCandidateCollection::const_iterator isoMatched=isoPixelTracks->begin();
       bool matched=false;
 
       //Note: this part needs to be changed: obtain isoPixelTracks from isoProdTracks (should be a way)
       for (reco::IsolatedPixelTrackCandidateCollection::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;
         }
     }
       if (!matched) continue;
       
       if (trit->hitPattern().numberOfValidHits()<MinNTrackHitsBarrel) continue;
       if (fabs(trit->eta())>1.47&&trit->hitPattern().numberOfValidStripTECHits()<MinNTECHitsEndcap) continue;
       
 
       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_);
 
         
       xTrkHcal=info.trkGlobPosAtHcal.x();
       yTrkHcal=info.trkGlobPosAtHcal.y();
       zTrkHcal=info.trkGlobPosAtHcal.z();
       
       xTrkEcal=info.trkGlobPosAtEcal.x();
       yTrkEcal=info.trkGlobPosAtEcal.y();
       zTrkEcal=info.trkGlobPosAtEcal.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 trackMomAtHcal = info.trkMomAtHcal;
     
       PxTrkHcal = trackMomAtHcal.x();
       PyTrkHcal = trackMomAtHcal.y();
       PzTrkHcal = trackMomAtHcal.z();
       
       //PxTrkHcal=0;
       //PyTrkHcal=0;
       //PzTrkHcal=0;
 
       GlobalPoint gPointEcal(xTrkEcal,yTrkEcal,zTrkEcal);
       GlobalPoint gPointHcal(xTrkHcal,yTrkHcal,zTrkHcal);
       
       //        emEnergy = isoMatched->energyIn();
       emEnergy =  ecalEnergyInCone(gPointEcal, EcalCone, Hitecal, geo);
       
       emRingEnergy = ecalEnergyInCone(gPointEcal, EcalConeOuter, Hitecal, geo) - ecalEnergyInCone(gPointEcal, EcalCone, Hitecal, geo);
 
       
     int iphitrue = -10;
     int ietatrue = 100;
     
     const HcalDetId tempId = gHcal->getClosestCell(gPointHcal);
     ietatrue = tempId.ieta();
     iphitrue = tempId.iphi();
 
 
       //container for used recHits
     std::vector<DetId> usedHits; 
       //      
       //clear usedHits
       usedHits.clear();
 
 
       // Find Hcal RecHit with maximum energy and collect other information
       MaxHit_struct MaxHit;
       MaxHit.hitenergy=-100;
       
 
       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;
       DetId hitId=hhit->id();
       for (uint32_t i=0; i<usedHits.size(); i++)
         {
           if (usedHits[i]==hitId) hitIsUsed=true;
         }
       if (hitIsUsed) continue;
       usedHits.push_back(hitId);
       //
 
       // rof 16.05.2008 start: include the possibility for recalibration
       float recal = 1;
       // rof end
 
       GlobalPoint pos = geo->getPosition(hhit->detid());
 
       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 distAtHcal =  getDistInPlaneTrackDir(gPointHcal, trackMomAtHcal, pos);
       //double distAtHcal =  getDistInPlaneSimple(gPointHcal, pos);
 
       //if(dr<associationConeSize_)      
       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();
           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.dr   = distAtHcal;
 
           MaxHit.posMax = geo->getPosition(hhit->detid());
           
         }
         }
     }
       
 
       Bool_t passCuts = kFALSE;
       if(calEnergy > energyMinIso && calEnergy < energyMaxIso && emEnergy < energyECALmip && 
      isoMatched->maxPtPxl() < maxPNear && abs(MaxHit.ietahitm)<30 && MaxHit.hitenergy > 0. && 
           MaxHit.dr < hottestHitDistance && emRingEnergy<8.){ 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;
           DetId hitId=hhit->id();
           for (uint32_t i=0; i<usedHits.size(); i++)
         {
           if (usedHits[i]==hitId) hitIsUsed=true;
         }
           if (hitIsUsed) 
         {
           //cout<<"Hit is Used!   ieta: "<< (hhit->id()).ieta()<<"  iphi: "<<(hhit->id()).iphi()<<"  depth: "<<(hhit->id()).depth()<<"  hashed_index: "<<hitHashedIndex<<endl;
           continue;}
           
           usedHits.push_back(hitId);
           //
 
           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" && getDistInPlaneTrackDir(gPointHcal, trackMomAtHcal, pos2) < calibrationConeSize_) {
             //if (AxB_=="Cone" && getDistInPlaneSimple(gPointHcal,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;
 
     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.getByToken(tok_ho_,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;
         DetId hitId=hoItr->id();
         for (uint32_t i=0; i<usedHits.size(); i++)
           {
         if (usedHits[i]==hitId) hitIsUsed=true;
           }
         if (hitIsUsed) continue;
 
         usedHits.push_back(hitId);
 
         /*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 )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 598 of file HcalIsoTrkAnalyzer.cc.

References convertSQLiteXML::runNumber, and MainPageGenerator::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("PxTrkHcal", &PxTrkHcal, "PxTrkHcal/F");
     tree->Branch("PyTrkHcal", &PyTrkHcal, "PyTrkHcal/F");
     tree->Branch("PzTrkHcal", &PzTrkHcal, "PzTrkHcal/F");
 
 
     tree->Branch("xTrkHcal", &xTrkHcal, "xTrkHcal/F");
     tree->Branch("yTrkHcal", &yTrkHcal, "yTrkHcal/F");
     tree->Branch("zTrkHcal", &zTrkHcal, "zTrkHcal/F");
 
     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 )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 642 of file HcalIsoTrkAnalyzer.cc.

                            {
 
     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
 
 }

Member Data Documentation

bool HcalIsoTrkAnalyzer::allowMissingInputs_

private

Definition at line 126 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::associationConeSize_

private

Definition at line 124 of file HcalIsoTrkAnalyzer.cc.

string HcalIsoTrkAnalyzer::AxB_

private

Definition at line 125 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::calibrationConeSize_

private

Definition at line 124 of file HcalIsoTrkAnalyzer.cc.

UInt_t HcalIsoTrkAnalyzer::cell

private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::cellEnergy

private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

TClonesArray* HcalIsoTrkAnalyzer::cells

private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

TRefArray* HcalIsoTrkAnalyzer::cells3x3

private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

TRefArray* HcalIsoTrkAnalyzer::cellsPF

private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::EcalCone

private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::EcalConeOuter

private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::emEnergy

private

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::emRingEnergy

private

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::energyECALmip

private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::energyMaxIso

private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::energyMinIso

private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::etVetoJet

private

Definition at line 163 of file HcalIsoTrkAnalyzer.cc.

UInt_t HcalIsoTrkAnalyzer::eventNumber

private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by Vispa.Plugins.EdmBrowser.EdmDataAccessor.EdmDataAccessor::setFilterBranches().

float HcalIsoTrkAnalyzer::eventWeight

private

Definition at line 134 of file HcalIsoTrkAnalyzer.cc.

const CaloGeometry* HcalIsoTrkAnalyzer::geo

private

Definition at line 110 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::hottestHitDistance

private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Int_t HcalIsoTrkAnalyzer::iEtaHit

private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

std::ofstream HcalIsoTrkAnalyzer::input_to_L3

private

Definition at line 167 of file HcalIsoTrkAnalyzer.cc.

UInt_t HcalIsoTrkAnalyzer::iPhiHit

private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::m_ebInstance

private

Definition at line 119 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::m_ecalLabel

private

Definition at line 118 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::m_eeInstance

private

Definition at line 120 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::m_hcalLabel

private

Definition at line 121 of file HcalIsoTrkAnalyzer.cc.

int HcalIsoTrkAnalyzer::m_histoFlag

private

Definition at line 122 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::m_inputTrackLabel

private

Definition at line 117 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::maxPNear

private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

int HcalIsoTrkAnalyzer::MinNTECHitsEndcap

private

Definition at line 133 of file HcalIsoTrkAnalyzer.cc.

int HcalIsoTrkAnalyzer::MinNTrackHitsBarrel

private

Definition at line 133 of file HcalIsoTrkAnalyzer.cc.

int HcalIsoTrkAnalyzer::nIterations

private

Definition at line 133 of file HcalIsoTrkAnalyzer.cc.

UInt_t HcalIsoTrkAnalyzer::numberOfCells

private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

string HcalIsoTrkAnalyzer::outputFileName_

private

Definition at line 127 of file HcalIsoTrkAnalyzer.cc.

TrackAssociatorParameters HcalIsoTrkAnalyzer::parameters_

private

Definition at line 108 of file HcalIsoTrkAnalyzer.cc.

Referenced by Mixins._TypedParameterizable::clone(), Types.PSet::clone(), and Mixins._TypedParameterizable::copy().

Float_t HcalIsoTrkAnalyzer::probeJetEmFrac

private

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

TLorentzVector* HcalIsoTrkAnalyzer::probeJetP4

private

Definition at line 162 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::PxTrkHcal

private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::PyTrkHcal

private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::PzTrkHcal

private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

TFile* HcalIsoTrkAnalyzer::rootFile

private

Definition at line 139 of file HcalIsoTrkAnalyzer.cc.

UInt_t HcalIsoTrkAnalyzer::runNumber

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::rvert

private

Definition at line 130 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::tagJetEmFrac

private

Definition at line 164 of file HcalIsoTrkAnalyzer.cc.

TLorentzVector* HcalIsoTrkAnalyzer::tagJetP4

private

Definition at line 161 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::targetE

private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_ecal_

private

Definition at line 113 of file HcalIsoTrkAnalyzer.cc.

edm::EDGetTokenT<HBHERecHitCollection> HcalIsoTrkAnalyzer::tok_hbhe_

private

Definition at line 111 of file HcalIsoTrkAnalyzer.cc.

edm::EDGetTokenT<HORecHitCollection> HcalIsoTrkAnalyzer::tok_ho_

private

Definition at line 112 of file HcalIsoTrkAnalyzer.cc.

edm::EDGetTokenT<reco::TrackCollection> HcalIsoTrkAnalyzer::tok_track1_

private

Definition at line 115 of file HcalIsoTrkAnalyzer.cc.

edm::EDGetTokenT<reco::IsolatedPixelTrackCandidateCollection> HcalIsoTrkAnalyzer::tok_track_

private

Definition at line 114 of file HcalIsoTrkAnalyzer.cc.

TrackDetectorAssociator HcalIsoTrkAnalyzer::trackAssociator_

private

Definition at line 107 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::trackEta

private

Definition at line 129 of file HcalIsoTrkAnalyzer.cc.

double HcalIsoTrkAnalyzer::trackPhi

private

Definition at line 129 of file HcalIsoTrkAnalyzer.cc.

TTree* HcalIsoTrkAnalyzer::tree

private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by eventstfile.Events::__iter__(), svgfig.SVG::__str__(), MetTreeProducer.MetTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::fillTree(), MetTreeProducer.MetTreeProducer::process(), python.cmstools.EventTree::SetAlias(), and eventstfile.Events::to().

Float_t HcalIsoTrkAnalyzer::xTrkEcal

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::xTrkHcal

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::yTrkEcal

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::yTrkHcal

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::zTrkEcal

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Float_t HcalIsoTrkAnalyzer::zTrkHcal

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation