Inheritance diagram for HcalIsoTrkAnalyzer:

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

	~HcalIsoTrkAnalyzer ()

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &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 ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

virtual void	beginJob ()

virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	beginRun (edm::Run const &, edm::EventSetup const &) override

double	dR (math::XYZTLorentzVector &, math::XYZTLorentzVector &)

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	endRun (edm::Run const &, edm::EventSetup const &) override

int	fillTree (std::vector< math::XYZTLorentzVector > &vecL1, std::vector< math::XYZTLorentzVector > &vecL3, math::XYZPoint &leadPV, std::vector< spr::propagatedTrackDirection > &trkCaloDirections, const CaloGeometry *geo, edm::Handle< EcalRecHitCollection > &barrelRecHitsHandle, edm::Handle< EcalRecHitCollection > &endcapRecHitsHandle, edm::Handle< HBHERecHitCollection > &hbhe)

Private Attributes
double	a_charIsoR_

double	a_coneR_

double	a_mipR_

int	dataType_

double	eEcalMax_

double	eIsolation_

edm::Service< TFileService >	fs

double	hcalScale_

HLTConfigProvider	hltConfig_

bool	ignoreTrigger_

std::string	l1Filter_

std::string	l2Filter_

std::string	l3Filter_

std::string	labelEB_

std::string	labelEE_

std::string	labelGenTrack_

std::string	labelHBHE_

std::string	labelRecVtx_

double	maxRestrictionPt_

unsigned int	nRun_

std::vector< double >	pbin

std::string	processName_

double	pTrackMin_

spr::trackSelectionParameters	selectionParameter_

double	slopeRestrictionPt_

int	t_DataType

std::vector< unsigned int > *	t_DetIds

double	t_eHcal

double	t_eMipDR

int	t_Event

double	t_EventWeight

int	t_goodPV

std::vector< double > *	t_HitEnergies

double	t_hmaxNearP

int	t_ieta

std::vector< int > *	t_ietaAll

std::vector< int > *	t_ietaGood

double	t_l1eta

double	t_l1phi

double	t_l1pt

double	t_l3eta

double	t_l3phi

double	t_l3pt

double	t_mindR1

double	t_mindR2

double	t_p

bool	t_qltyFlag

bool	t_qltyMissFlag

bool	t_qltyPVFlag

int	t_Run

bool	t_selectTk

int	t_Tracks

int	t_TracksProp

int	t_TracksSaved

std::vector< bool > *	t_trgbits

std::string	theTrackQuality_

edm::InputTag	theTriggerResultsLabel_

edm::EDGetTokenT< reco::BeamSpot >	tok_bs_

edm::EDGetTokenT < EcalRecHitCollection >	tok_EB_

edm::EDGetTokenT < EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT < reco::TrackCollection >	tok_genTrack_

edm::EDGetTokenT < HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT < reco::VertexCollection >	tok_recVtx_

edm::EDGetTokenT < trigger::TriggerEvent >	tok_trigEvt_

edm::EDGetTokenT < edm::TriggerResults >	tok_trigRes_

TTree *	tree

TTree *	tree2

edm::InputTag	triggerEvent_

std::vector< std::string >	trigNames_

bool	useRaw_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Definition at line 68 of file HcalIsoTrkAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 130 of file HcalIsoTrkAnalyzer.cc.

                                                                      : 
   nRun_(0) {
 
   usesResource("TFileService");
 
   //now do whatever initialization is needed
   const double isolationRadius(28.9);
   trigNames_                          = iConfig.getParameter<std::vector<std::string> >("Triggers");
   theTrackQuality_                    = iConfig.getParameter<std::string>("TrackQuality");
   processName_                        = iConfig.getParameter<std::string>("ProcessName");
   l1Filter_                           = iConfig.getParameter<std::string>("L1Filter");
   l2Filter_                           = iConfig.getParameter<std::string>("L2Filter");
   l3Filter_                           = iConfig.getParameter<std::string>("L3Filter");
   reco::TrackBase::TrackQuality trackQuality_=reco::TrackBase::qualityByName(theTrackQuality_);
   maxRestrictionPt_                   = iConfig.getParameter<double>("MinTrackPt");
   slopeRestrictionPt_                 = iConfig.getParameter<double>("SlopeTrackPt");
   selectionParameter_.minPt           = (slopeRestrictionPt_ > 0) ? (maxRestrictionPt_ - 2.5/slopeRestrictionPt_) : maxRestrictionPt_;
   selectionParameter_.minQuality      = trackQuality_;
   selectionParameter_.maxDxyPV        = iConfig.getParameter<double>("MaxDxyPV");
   selectionParameter_.maxDzPV         = iConfig.getParameter<double>("MaxDzPV");
   selectionParameter_.maxChi2         = iConfig.getParameter<double>("MaxChi2");
   selectionParameter_.maxDpOverP      = iConfig.getParameter<double>("MaxDpOverP");
   selectionParameter_.minOuterHit     = iConfig.getParameter<int>("MinOuterHit");
   selectionParameter_.minLayerCrossed = iConfig.getParameter<int>("MinLayerCrossed");
   selectionParameter_.maxInMiss       = iConfig.getParameter<int>("MaxInMiss");
   selectionParameter_.maxOutMiss      = iConfig.getParameter<int>("MaxOutMiss");
   a_coneR_                            = iConfig.getParameter<double>("ConeRadius");
   a_charIsoR_                         = a_coneR_ + isolationRadius;
   a_mipR_                             = iConfig.getParameter<double>("ConeRadiusMIP");
   pTrackMin_                          = iConfig.getParameter<double>("MinimumTrackP");
   eEcalMax_                           = iConfig.getParameter<double>("MaximumEcalEnergy");
   eIsolation_                         = iConfig.getParameter<double>("IsolationEnergy");
   triggerEvent_                       = iConfig.getParameter<edm::InputTag>("TriggerEventLabel");
   theTriggerResultsLabel_             = iConfig.getParameter<edm::InputTag>("TriggerResultLabel");
   labelGenTrack_                      = iConfig.getParameter<std::string>("TrackLabel");
   labelRecVtx_                        = iConfig.getParameter<std::string>("VertexLabel");
   labelEB_                            = iConfig.getParameter<std::string>("EBRecHitLabel");
   labelEE_                            = iConfig.getParameter<std::string>("EERecHitLabel");
   labelHBHE_                          = iConfig.getParameter<std::string>("HBHERecHitLabel");
   std::string labelBS                 = iConfig.getParameter<std::string>("BeamSpotLabel");
   std::string modnam                  = iConfig.getUntrackedParameter<std::string>("ModuleName","");
   std::string prdnam                  = iConfig.getUntrackedParameter<std::string>("ProducerName","");
   ignoreTrigger_                      = iConfig.getUntrackedParameter<bool>("IgnoreTriggers", false);
   useRaw_                             = iConfig.getUntrackedParameter<bool>("UseRaw", false);
   hcalScale_                          = iConfig.getUntrackedParameter<double>("HcalScale", 1.0);
   dataType_                           = iConfig.getUntrackedParameter<int>("DataType", 0);
 
   // define tokens for access
   tok_trigEvt_  = consumes<trigger::TriggerEvent>(triggerEvent_);
   tok_trigRes_  = consumes<edm::TriggerResults>(theTriggerResultsLabel_);
   tok_bs_       = consumes<reco::BeamSpot>(labelBS);
   tok_genTrack_ = consumes<reco::TrackCollection>(labelGenTrack_);
   if (modnam == "") {
     tok_recVtx_   = consumes<reco::VertexCollection>(labelRecVtx_);
     tok_EB_       = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit",labelEB_));
     tok_EE_       = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit",labelEE_));
     tok_hbhe_     = consumes<HBHERecHitCollection>(labelHBHE_);
     edm::LogInfo("HcalIsoTrack") << "Labels used " << triggerEvent_ << " "
                  << theTriggerResultsLabel_ << " "
                  << labelBS << " " << labelRecVtx_ << " " 
                  << labelGenTrack_ << " " 
                  << edm::InputTag("ecalRecHit",labelEB_) << " " 
                  << edm::InputTag("ecalRecHit",labelEE_) << " "
                  << labelHBHE_;
   } else {
     tok_recVtx_   = consumes<reco::VertexCollection>(edm::InputTag(modnam,labelRecVtx_,prdnam));
     tok_EB_       = consumes<EcalRecHitCollection>(edm::InputTag(modnam,labelEB_,prdnam));
     tok_EE_       = consumes<EcalRecHitCollection>(edm::InputTag(modnam,labelEE_,prdnam));
     tok_hbhe_     = consumes<HBHERecHitCollection>(edm::InputTag(modnam,labelHBHE_,prdnam));
     edm::LogInfo("HcalIsoTrack") << "Labels used "   << triggerEvent_ 
                  << "\n            " << theTriggerResultsLabel_ 
                  << "\n            " << labelBS 
                  << "\n            " << edm::InputTag(modnam,labelRecVtx_,prdnam)
                  << "\n            " << labelGenTrack_
                  << "\n            " << edm::InputTag(modnam,labelEB_,prdnam)
                  << "\n            " << edm::InputTag(modnam,labelEE_,prdnam)
                  << "\n            " << edm::InputTag(modnam,labelHBHE_,prdnam);
   }
 
   edm::LogInfo("HcalIsoTrack") <<"Parameters read from config file \n" 
                    <<"\t minPt "           << maxRestrictionPt_ <<":" << slopeRestrictionPt_
                    <<"\t theTrackQuality " << theTrackQuality_
                    <<"\t minQuality "      << selectionParameter_.minQuality
                    <<"\t maxDxyPV "        << selectionParameter_.maxDxyPV          
                    <<"\t maxDzPV "         << selectionParameter_.maxDzPV          
                    <<"\t maxChi2 "         << selectionParameter_.maxChi2          
                    <<"\t maxDpOverP "      << selectionParameter_.maxDpOverP
                    <<"\t minOuterHit "     << selectionParameter_.minOuterHit
                    <<"\t minLayerCrossed " << selectionParameter_.minLayerCrossed
                    <<"\t maxInMiss "       << selectionParameter_.maxInMiss
                    <<"\t maxOutMiss "      << selectionParameter_.maxOutMiss
                    <<"\t a_coneR "         << a_coneR_
                    <<"\t a_charIsoR "      << a_charIsoR_
                    <<"\t a_mipR "          << a_mipR_
                    <<"\t hcalScale_ "      << hcalScale_
                    <<"\t useRaw_ "         << useRaw_
                    <<"\t ignoreTrigger_ "  << ignoreTrigger_
                    <<"\t dataType_      "  << dataType_;
   edm::LogInfo("HcalIsoTrack") << "Process " << processName_ << " L1Filter:" 
                    << l1Filter_ << " L2Filter:" << l2Filter_
                    << " L3Filter:" << l3Filter_;
   for (unsigned int k=0; k<trigNames_.size(); ++k) {
     edm::LogInfo("HcalIsoTrack") << "Trigger[" << k << "] " << trigNames_[k];
   }
 }

HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer ( )

Definition at line 236 of file HcalIsoTrkAnalyzer.cc.

236 { }

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 238 of file HcalIsoTrkAnalyzer.cc.

                                                                                     {
 
   t_Run      = iEvent.id().run();
   t_Event    = iEvent.id().event();
   t_DataType = dataType_;
 #ifdef DebugLog
   edm::LogInfo("HcalIsoTrack") << "Run " << t_Run << " Event " << t_Event 
                    << " type " << t_DataType << " Luminosity " 
                    << iEvent.luminosityBlock() << " Bunch " 
                    << iEvent.bunchCrossing();
 #endif
   //Get magnetic field and ECAL channel status
   edm::ESHandle<MagneticField> bFieldH;
   iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
   const MagneticField *bField = bFieldH.product();
 
   // get handles to calogeometry and calotopology
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   const CaloGeometry* geo = pG.product();
   
   bool okC(true);
   //Get track collection
   edm::Handle<reco::TrackCollection> trkCollection;
   iEvent.getByToken(tok_genTrack_, trkCollection);
   reco::TrackCollection::const_iterator trkItr;
   if (!trkCollection.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelGenTrack_;
     okC = false;
   }
  
   //event weight for FLAT sample
   t_EventWeight = 1.0;
 
   //Define the best vertex and the beamspot
   edm::Handle<reco::VertexCollection> recVtxs;
   iEvent.getByToken(tok_recVtx_, recVtxs);  
   edm::Handle<reco::BeamSpot> beamSpotH;
   iEvent.getByToken(tok_bs_, beamSpotH);
   math::XYZPoint leadPV(0,0,0);
   t_goodPV = 0;
   if (recVtxs.isValid() && recVtxs->size()>0) {
     for (unsigned int k=0; k<recVtxs->size(); ++k) {
       if (!((*recVtxs)[k].isFake()) && ((*recVtxs)[k].ndof() > 4)) {
     if (t_goodPV == 0) leadPV = math::XYZPoint((*recVtxs)[k].x(),(*recVtxs)[k].y(),(*recVtxs)[k].z());
     t_goodPV++;
       }
     }
   } 
   if (t_goodPV == 0 && beamSpotH.isValid()) {
     leadPV = beamSpotH->position();
   }
 #ifdef DebugLog
   edm::LogInfo("HcalIsoTrack") << "Primary Vertex " << leadPV << " out of "
                    << t_goodPV << " vertex";
   if (beamSpotH.isValid()) {
     edm::LogInfo("HcalIsoTrack") << " Beam Spot " << beamSpotH->position();
   }
 #endif  
   // RecHits
   edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
   if (!barrelRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEB_;
     okC = false;
   }
   edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
   if (!endcapRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEE_;
     okC = false;
   }
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(tok_hbhe_, hbhe);
   if (!hbhe.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelHBHE_;
     okC = false;
   }
       
   //Propagate tracks to calorimeter surface)
   std::vector<spr::propagatedTrackDirection> trkCaloDirections;
   spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_,
              trkCaloDirections, false);
   std::vector<math::XYZTLorentzVector> vecL1, vecL3;
   t_Tracks     = trkCollection->size();
   t_TracksProp = trkCaloDirections.size();
   t_ietaAll->clear(); t_ietaGood->clear();
 #ifdef DebugLog
   edm::LogInfo("HcalIsoTrack") << "# of propagated tracks " << t_TracksProp
                    << " out of " << t_Tracks << " with Trigger "
                    << ignoreTrigger_;
 #endif
   //Trigger
   t_trgbits->clear();
   t_TracksSaved = 0;
   for (unsigned int i=0; i<trigNames_.size(); ++i) t_trgbits->push_back(false);
   if (ignoreTrigger_) {
     t_l1pt  = t_l1eta = t_l1phi = 0;
     t_l3pt  = t_l3eta = t_l3phi = 0;
     t_TracksSaved = fillTree(vecL1, vecL3, leadPV, trkCaloDirections, geo, 
                  barrelRecHitsHandle, endcapRecHitsHandle, hbhe);
   } else {
     trigger::TriggerEvent triggerEvent;
     edm::Handle<trigger::TriggerEvent> triggerEventHandle;
     iEvent.getByToken(tok_trigEvt_, triggerEventHandle);
     if (!triggerEventHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Error! Can't get the product "
                     << triggerEvent_.label() ;
     } else if (okC) {
       triggerEvent = *(triggerEventHandle.product());
       const trigger::TriggerObjectCollection& TOC(triggerEvent.getObjects());
       edm::Handle<edm::TriggerResults> triggerResults;
       iEvent.getByToken(tok_trigRes_, triggerResults);
       bool done(false);
       if (triggerResults.isValid()) {
     std::vector<std::string> modules;
     const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults);
     const std::vector<std::string> & triggerNames_ = triggerNames.triggerNames();
     int ntksave(0);
     for (unsigned int iHLT=0; iHLT<triggerResults->size(); iHLT++) {
       bool ok(false);
       int hlt    = triggerResults->accept(iHLT);
       if (trigNames_.size() > 0) {
         for (unsigned int i=0; i<trigNames_.size(); ++i) {
           if (triggerNames_[iHLT].find(trigNames_[i].c_str())!=std::string::npos) {
         t_trgbits->at(i) = (hlt>0);
         if (hlt > 0) ok = true;
 #ifdef DebugLog
         edm::LogInfo("HcalIsoTrack") << "This trigger "
                          << triggerNames_[iHLT] << " Flag " 
                          << hlt << ":" << ok;
 #endif
           }
         }
       } else {
         ok = (hlt>0);
       }
 
       if (ok) {
         unsigned int triggerindx = hltConfig_.triggerIndex(triggerNames_[iHLT]);
         const std::vector<std::string>& moduleLabels(hltConfig_.moduleLabels(triggerindx));
         std::vector<math::XYZTLorentzVector> vecL2;
         vecL1.clear(); vecL3.clear();
         //loop over all trigger filters in event (i.e. filters passed)
         for (unsigned int ifilter=0; ifilter<triggerEvent.sizeFilters();
          ++ifilter) {  
           std::vector<int> Keys;
           std::string label = triggerEvent.filterTag(ifilter).label();
           //loop over keys to objects passing this filter
           for (unsigned int imodule=0; imodule<moduleLabels.size(); imodule++) {
         if (label.find(moduleLabels[imodule]) != std::string::npos) {
 #ifdef DebugLog
           edm::LogInfo("HcalIsoTrack") << "FilterName " << label;
 #endif
           for (unsigned int ifiltrKey=0; ifiltrKey<triggerEvent.filterKeys(ifilter).size(); ++ifiltrKey) {
             Keys.push_back(triggerEvent.filterKeys(ifilter)[ifiltrKey]);
             const trigger::TriggerObject& TO(TOC[Keys[ifiltrKey]]);
             math::XYZTLorentzVector v4(TO.px(), TO.py(), TO.pz(), TO.energy());
             if (label.find(l2Filter_) != std::string::npos) {
               vecL2.push_back(v4);
             } else if (label.find(l3Filter_) != std::string::npos) {
               vecL3.push_back(v4);
             } else if (label.find(l1Filter_) != std::string::npos || l1Filter_ == "") {
               vecL1.push_back(v4);
             }
 #ifdef DebugLog
             edm::LogInfo("HcalIsoTrack") << "key " <<ifiltrKey<<" : pt "
                          << TO.pt() << " eta "<<TO.eta()
                          << " phi " <<TO.phi()<<" mass "
                          << TO.mass() <<" Id "<<TO.id();
 #endif
           }
 #ifdef DebugLog
           edm::LogInfo("HcalIsoTrack") << "sizes " << vecL1.size() <<":"
                            << vecL2.size() <<":" <<vecL3.size();
 #endif
         }
           }
         }
         math::XYZTLorentzVector mindRvec1;
         double mindR1(999);
         for (unsigned int i=0; i<vecL2.size(); i++) {
           double dr   = dR(vecL1[0],vecL2[i]);
 #ifdef DebugLog
           edm::LogInfo("HcalIsoTrack") << "lvl2[" << i << "] dR " << dr;
 #endif
           if (dr<mindR1) {
         mindR1    = dr;
         mindRvec1 = vecL2[i];
           }
         }
 #ifdef DebugLog
         edm::LogInfo("HcalIsoTrack") << "L2 object closest to L1 "
                      << mindRvec1 << " at Dr " << mindR1;
 #endif
       
         if (vecL1.size()>0) {
           t_l1pt  = vecL1[0].pt();
           t_l1eta = vecL1[0].eta();
           t_l1phi = vecL1[0].phi();
         } else {
           t_l1pt  = t_l1eta = t_l1phi = 0;
         }
         if (vecL3.size()>0) {
           t_l3pt  = vecL3[0].pt();
           t_l3eta = vecL3[0].eta();
           t_l3phi = vecL3[0].phi();
         } else {
           t_l3pt  = t_l3eta = t_l3phi = 0;
         }
         // Now fill in the tree for each selected track
         if (!done) {
           ntksave = fillTree(vecL1, vecL3, leadPV, trkCaloDirections, geo, 
                  barrelRecHitsHandle,endcapRecHitsHandle, hbhe);
           done = true;
         } 
         t_TracksSaved += ntksave;
       }
     }
       }
     }
   }
   tree2->Fill();
 }

void HcalIsoTrkAnalyzer::beginJob ( void )

privatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 465 of file HcalIsoTrkAnalyzer.cc.

References fs, TFileService::make(), t_DataType, t_DetIds, t_eHcal, t_eMipDR, t_Event, t_EventWeight, t_goodPV, t_HitEnergies, t_hmaxNearP, t_ieta, t_ietaAll, t_ietaGood, t_l1eta, t_l1phi, t_l1pt, t_l3eta, t_l3phi, t_l3pt, t_mindR1, t_mindR2, t_p, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_Run, t_selectTk, t_Tracks, t_TracksProp, t_TracksSaved, t_trgbits, tree, and tree2.

                                   {
 
   tree = fs->make<TTree>("CalibTree", "CalibTree");
      
   tree->Branch("t_Run",         &t_Run,         "t_Run/I");
   tree->Branch("t_Event",       &t_Event,       "t_Event/I");
   tree->Branch("t_DataType",    &t_DataType,    "t_DataType/I");
   tree->Branch("t_ieta",        &t_ieta,        "t_ieta/I");
   tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
   tree->Branch("t_goodPV",      &t_goodPV,      "t_goodPV/I");
   tree->Branch("t_l1pt",        &t_l1pt,        "t_l1pt/D");
   tree->Branch("t_l1eta",       &t_l1eta,       "t_l1eta/D");
   tree->Branch("t_l1phi",       &t_l1phi,       "t_l1phi/D"); 
   tree->Branch("t_l3pt",        &t_l3pt,        "t_l3pt/D");
   tree->Branch("t_l3eta",       &t_l3eta,       "t_l3eta/D"); 
   tree->Branch("t_l3phi",       &t_l3phi,       "t_l3phi/D");
   tree->Branch("t_p",           &t_p,           "t_p/D");
   tree->Branch("t_mindR1",      &t_mindR1,      "t_mindR1/D");
   tree->Branch("t_mindR2",      &t_mindR2,      "t_mindR2/D");
   tree->Branch("t_eMipDR",      &t_eMipDR,      "t_eMipDR/D");
   tree->Branch("t_eHcal",       &t_eHcal,       "t_eHcal/D");
   tree->Branch("t_hmaxNearP",   &t_hmaxNearP,   "t_hmaxNearP/D");
   tree->Branch("t_selectTk",    &t_selectTk,    "t_selectTk/O");
   tree->Branch("t_qltyFlag",    &t_qltyFlag,    "t_qltyFlag/O");
   tree->Branch("t_qltyMissFlag",&t_qltyMissFlag,"t_qltyMissFlag/O");
   tree->Branch("t_qltyPVFlag",  &t_qltyPVFlag,  "t_qltyPVFlag/O)");
 
   t_DetIds      = new std::vector<unsigned int>();
   t_HitEnergies = new std::vector<double>();
   t_trgbits     = new std::vector<bool>();
   tree->Branch("t_DetIds",      "std::vector<unsigned int>", &t_DetIds);
   tree->Branch("t_HitEnergies", "std::vector<double>",       &t_HitEnergies);
   tree->Branch("t_trgbits",     "std::vector<bool>",         &t_trgbits); 
 
   tree2 = fs->make<TTree>("EventInfo", "Event Information");
      
   tree2->Branch("t_Tracks",      &t_Tracks,      "t_Tracks/I");
   tree2->Branch("t_TracksProp",  &t_TracksProp,  "t_TracksProp/I");
   tree2->Branch("t_TracksSaved", &t_TracksSaved, "t_TracksSaved/I");
   t_ietaAll     = new std::vector<int>();
   t_ietaGood    = new std::vector<int>();
   tree2->Branch("t_ietaAll",    "std::vector<int>",          &t_ietaAll);
   tree2->Branch("t_ietaGood",   "std::vector<int>",          &t_ietaGood);
 }

void HcalIsoTrkAnalyzer::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

privatevirtual

Definition at line 549 of file HcalIsoTrkAnalyzer.cc.

549 {}

void HcalIsoTrkAnalyzer::beginRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

overrideprivatevirtual

Definition at line 512 of file HcalIsoTrkAnalyzer.cc.

References hltConfig_, HLTConfigProvider::init(), gen::n, nRun_, processName_, edm::RunBase::run(), HLTConfigProvider::size(), HLTConfigProvider::triggerIndex(), and trigNames_.

                                                                                  {
   bool changed_(true);
   bool flag =  hltConfig_.init(iRun,iSetup,processName_,changed_);
   edm::LogInfo("HcalIsoTrack") << "Run[" << nRun_ << "] " << iRun.run() 
                    << " process " << processName_ << " init flag "
                    << flag << " change flag " << changed_;
 
   // check if trigger names in (new) config                       
   if (changed_) {
     changed_ = false;
 #ifdef DebugLog
     edm::LogInfo("HcalIsoTrack") << "New trigger menu found !!!";
 #endif
     const unsigned int n(hltConfig_.size());
     for (unsigned itrig=0; itrig<trigNames_.size(); itrig++) {
       unsigned int triggerindx = hltConfig_.triggerIndex(trigNames_[itrig]);
       if (triggerindx >= n) {
     edm::LogWarning("HcalIsoTrack") << trigNames_[itrig] << " " 
                     << triggerindx << " does not exist in "
                     << "the current menu";
 #ifdef DebugLog
       } else {
     edm::LogInfo("HcalIsoTrack") << trigNames_[itrig] << " " 
                      << triggerindx << " exists";
 #endif
       }
     }
   }
 }

double HcalIsoTrkAnalyzer::dR	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 680 of file HcalIsoTrkAnalyzer.cc.

References reco::deltaR().

Referenced by analyze(), and fillTree().

                                                                                       {
   return reco::deltaR(vec1.eta(),vec1.phi(),vec2.eta(),vec2.phi());
 }

void HcalIsoTrkAnalyzer::endLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

privatevirtual

Definition at line 551 of file HcalIsoTrkAnalyzer.cc.

551 {}

void HcalIsoTrkAnalyzer::endRun	(	edm::Run const &	iRun,
		edm::EventSetup const &
	)

overrideprivatevirtual

Definition at line 543 of file HcalIsoTrkAnalyzer.cc.

References nRun_, and edm::RunBase::run().

                                                                         {
   nRun_++;
   edm::LogInfo("HcalIsoTrack") << "endRun[" << nRun_ << "] " << iRun.run();
 }

void HcalIsoTrkAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 553 of file HcalIsoTrkAnalyzer.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                     {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

int HcalIsoTrkAnalyzer::fillTree	(	std::vector< math::XYZTLorentzVector > &	vecL1,
		std::vector< math::XYZTLorentzVector > &	vecL3,
		math::XYZPoint &	leadPV,
		std::vector< spr::propagatedTrackDirection > &	trkCaloDirections,
		const CaloGeometry *	geo,
		edm::Handle< EcalRecHitCollection > &	barrelRecHitsHandle,
		edm::Handle< EcalRecHitCollection > &	endcapRecHitsHandle,
		edm::Handle< HBHERecHitCollection > &	hbhe
	)

private

Definition at line 561 of file HcalIsoTrkAnalyzer.cc.

References a_charIsoR_, a_coneR_, a_mipR_, funct::abs(), spr::chargeIsolationCone(), dR(), spr::eCone_ecal(), spr::eCone_hcal(), eEcalMax_, eIsolation_, reco::TrackBase::eta(), spr::goodTrack(), hcalScale_, HcalDetId::ieta(), relval_2017::k, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, maxRestrictionPt_, spr::trackSelectionParameters::minPt, nTracks(), reco::TrackBase::p(), reco::TrackBase::phi(), reco::TrackBase::pt(), pTrackMin_, reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), selectionParameter_, slopeRestrictionPt_, t_DetIds, t_eHcal, t_eMipDR, t_HitEnergies, t_hmaxNearP, t_ieta, t_ietaAll, t_ietaGood, t_mindR1, t_mindR2, t_p, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_selectTk, t_trgbits, tree, and useRaw_.

Referenced by analyze(), and core.AutoFillTreeProducer.AutoFillTreeProducer::process().

                                                         {
 
   int nsave = 0;
   //Loop over tracks
   std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
   unsigned int nTracks(0), nselTracks(0);
   for (trkDetItr = trkCaloDirections.begin(),nTracks=0; 
        trkDetItr != trkCaloDirections.end(); trkDetItr++,nTracks++) {
     const reco::Track* pTrack = &(*(trkDetItr->trkItr));
     math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), 
                    pTrack->pz(), pTrack->p());
 #ifdef DebugLog
     edm::LogInfo("HcalIsoTrack") << "This track : " << nTracks 
                  << " (pt|eta|phi|p) :" << pTrack->pt()
                  << "|" << pTrack->eta() << "|" 
                  << pTrack->phi() << "|" <<pTrack->p();
 #endif
     t_mindR2 = 999;
     for (unsigned int k=0; k<vecL3.size(); ++k) {
       double dr   = dR(vecL3[k],v4); 
       if (dr<t_mindR2) {
     t_mindR2  = dr;
       }
     }
     t_mindR1 = (vecL1.size() > 0) ? dR(vecL1[0],v4) : 999;
 #ifdef DebugLog
     edm::LogInfo("HcalIsoTrack") << "Closest L3 object at dr :" 
                  << t_mindR2 << " and from L1 " << t_mindR1;
 #endif      
     //Selection of good track
     selectionParameter_.minPt = maxRestrictionPt_ - std::abs(pTrack->eta())/slopeRestrictionPt_;
     t_selectTk = spr::goodTrack(pTrack,leadPV,selectionParameter_,false);
     spr::trackSelectionParameters oneCutParameters = selectionParameter_;
     oneCutParameters.maxDxyPV  = 10;
     oneCutParameters.maxDzPV   = 100;
     oneCutParameters.maxInMiss = 2;
     oneCutParameters.maxOutMiss= 2;
     bool qltyFlag  = spr::goodTrack(pTrack,leadPV,oneCutParameters,false);
     oneCutParameters           = selectionParameter_;
     oneCutParameters.maxDxyPV  = 10;
     oneCutParameters.maxDzPV   = 100;
     t_qltyMissFlag = spr::goodTrack(pTrack,leadPV,oneCutParameters,false);
     oneCutParameters           = selectionParameter_;
     oneCutParameters.maxInMiss = 2;
     oneCutParameters.maxOutMiss= 2;
     t_qltyPVFlag   = spr::goodTrack(pTrack,leadPV,oneCutParameters,false);
     t_ieta         = 0;
     t_p            = pTrack->p();
     if (trkDetItr->okHCAL) {
       HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
       t_ieta = detId.ieta();
       if (t_p > 40.0 && t_p <= 60.0) t_ietaAll->push_back(t_ieta);
     }
 #ifdef DebugLog
     edm::LogInfo("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " 
                  << qltyFlag << "|" << trkDetItr->okECAL
                  << "|" << trkDetItr->okHCAL;
 #endif
     t_qltyFlag = (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL);
     if (t_qltyFlag) {
       nselTracks++;
       int nRH_eMipDR(0), nNearTRKs(0), nRecHits(-999);
       t_eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle, 
                  endcapRecHitsHandle, trkDetItr->pointHCAL,
                  trkDetItr->pointECAL, a_mipR_, 
                  trkDetItr->directionECAL, nRH_eMipDR);
       t_DetIds->clear(); t_HitEnergies->clear();
       std::vector<DetId> ids;
       t_eHcal = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL, 
                 trkDetItr->pointECAL, a_coneR_, 
                 trkDetItr->directionHCAL,nRecHits, 
                 ids, *t_HitEnergies, useRaw_);
       t_eHcal *= hcalScale_;
       for (unsigned int k=0; k<ids.size(); ++k) {
     t_DetIds->push_back(ids[k].rawId());
       }
       t_hmaxNearP = spr::chargeIsolationCone(nTracks, trkCaloDirections,
                          a_charIsoR_, nNearTRKs, false);
 #ifdef DebugLog
       edm::LogInfo("HcalIsoTrack") << "This track : " << nTracks 
                    << " (pt|eta|phi|p) :"  << pTrack->pt() 
                    << "|" << pTrack->eta() << "|" 
                    << pTrack->phi() << "|" << t_p;
       edm::LogInfo("HcalIsoTrack") << "e_MIP " << t_eMipDR 
                    << " Chg Isolation " << t_hmaxNearP 
                    << " eHcal" << t_eHcal << " ieta " 
                    << t_ieta << " Quality " 
                    << t_qltyMissFlag << ":" 
                    << t_qltyPVFlag << ":" << t_selectTk;
       for (unsigned int lll=0;lll<t_DetIds->size();lll++) {
     edm::LogInfo("HcalIsoTrack") << "det id is = " <<t_DetIds->at(lll)
                      << "   hit enery is  = "  
                      << t_HitEnergies->at(lll) ;
       }
 #endif
       if (t_p>pTrackMin_ && t_eMipDR<eEcalMax_ && 
       t_hmaxNearP<eIsolation_) {
     tree->Fill();
     nsave++;
     if (t_p > 40.0 && t_p <= 60.0 && t_selectTk) t_ietaGood->push_back(t_ieta);
 #ifdef DebugLog
     for (unsigned int k=0; k<t_trgbits->size(); k++) {
       edm::LogInfo("HcalIsoTrack") << "trigger bit is  = " 
                        << t_trgbits->at(k);
     }
 #endif
       }
     }
   }
   return nsave;
 }

Member Data Documentation

double HcalIsoTrkAnalyzer::a_charIsoR_

private

Definition at line 101 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::a_coneR_

private

Definition at line 101 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::a_mipR_

private

Definition at line 101 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

int HcalIsoTrkAnalyzer::dataType_

private

Definition at line 99 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::eEcalMax_

private

Definition at line 102 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::eIsolation_

private

Definition at line 102 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

edm::Service<TFileService> HcalIsoTrkAnalyzer::fs

private

Definition at line 95 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob().

double HcalIsoTrkAnalyzer::hcalScale_

private

Definition at line 102 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

HLTConfigProvider HcalIsoTrkAnalyzer::hltConfig_

private

Definition at line 96 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginRun().

bool HcalIsoTrkAnalyzer::ignoreTrigger_

private

Definition at line 107 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::string HcalIsoTrkAnalyzer::l1Filter_

private

Definition at line 106 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::string HcalIsoTrkAnalyzer::l2Filter_

private

Definition at line 106 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::string HcalIsoTrkAnalyzer::l3Filter_

private

Definition at line 106 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::string HcalIsoTrkAnalyzer::labelEB_

private

Definition at line 104 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::string HcalIsoTrkAnalyzer::labelEE_

private

Definition at line 104 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::string HcalIsoTrkAnalyzer::labelGenTrack_

private

Definition at line 104 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::string HcalIsoTrkAnalyzer::labelHBHE_

private

Definition at line 105 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::string HcalIsoTrkAnalyzer::labelRecVtx_

private

Definition at line 104 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::maxRestrictionPt_

private

Definition at line 100 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

unsigned int HcalIsoTrkAnalyzer::nRun_

private

Definition at line 94 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and endRun().

std::vector<double> HcalIsoTrkAnalyzer::pbin

private

Definition at line 124 of file HcalIsoTrkAnalyzer.cc.

std::string HcalIsoTrkAnalyzer::processName_

private

Definition at line 105 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::pTrackMin_

private

Definition at line 102 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

spr::trackSelectionParameters HcalIsoTrkAnalyzer::selectionParameter_

private

Definition at line 98 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::slopeRestrictionPt_

private

Definition at line 100 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

int HcalIsoTrkAnalyzer::t_DataType

private

Definition at line 118 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::vector<unsigned int>* HcalIsoTrkAnalyzer::t_DetIds

private

Definition at line 123 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eHcal

private

Definition at line 121 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR

private

Definition at line 121 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_Event

private

Definition at line 118 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_EventWeight

private

Definition at line 119 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_goodPV

private

Definition at line 118 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::vector<double>* HcalIsoTrkAnalyzer::t_HitEnergies

private

Definition at line 124 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_hmaxNearP

private

Definition at line 121 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_ieta

private

Definition at line 118 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

std::vector<int>* HcalIsoTrkAnalyzer::t_ietaAll

private

Definition at line 127 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), beginJob(), and fillTree().

std::vector<int> * HcalIsoTrkAnalyzer::t_ietaGood

private

Definition at line 127 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_l1eta

private

Definition at line 119 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1phi

private

Definition at line 119 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1pt

private

Definition at line 119 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3eta

private

Definition at line 120 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3phi

private

Definition at line 120 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3pt

private

Definition at line 120 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_mindR1

private

Definition at line 120 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_mindR2

private

Definition at line 121 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_p

private

Definition at line 120 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_qltyFlag

private

Definition at line 122 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_qltyMissFlag

private

Definition at line 122 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_qltyPVFlag

private

Definition at line 122 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_Run

private

Definition at line 118 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

bool HcalIsoTrkAnalyzer::t_selectTk

private

Definition at line 122 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_Tracks

private

Definition at line 126 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_TracksProp

private

Definition at line 126 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_TracksSaved

private

Definition at line 126 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::vector<bool>* HcalIsoTrkAnalyzer::t_trgbits

private

Definition at line 125 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), beginJob(), and fillTree().

std::string HcalIsoTrkAnalyzer::theTrackQuality_

private

Definition at line 105 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::InputTag HcalIsoTrkAnalyzer::theTriggerResultsLabel_

private

Definition at line 103 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

edm::EDGetTokenT<reco::BeamSpot> HcalIsoTrkAnalyzer::tok_bs_

private

Definition at line 112 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_EB_

private

Definition at line 113 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_EE_

private

Definition at line 114 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

private

Definition at line 110 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<HBHERecHitCollection> HcalIsoTrkAnalyzer::tok_hbhe_

private

Definition at line 115 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<reco::VertexCollection> HcalIsoTrkAnalyzer::tok_recVtx_

private

Definition at line 111 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<trigger::TriggerEvent> HcalIsoTrkAnalyzer::tok_trigEvt_

private

Definition at line 108 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<edm::TriggerResults> HcalIsoTrkAnalyzer::tok_trigRes_

private

Definition at line 109 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

TTree* HcalIsoTrkAnalyzer::tree

private

Definition at line 117 of file HcalIsoTrkAnalyzer.cc.

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

TTree * HcalIsoTrkAnalyzer::tree2

private

Definition at line 117 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

edm::InputTag HcalIsoTrkAnalyzer::triggerEvent_

private

Definition at line 103 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

std::vector<std::string> HcalIsoTrkAnalyzer::trigNames_

private

Definition at line 97 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), beginRun(), and HcalIsoTrkAnalyzer().

bool HcalIsoTrkAnalyzer::useRaw_

private

Definition at line 107 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation