Inheritance diagram for HcalIsoTrkSimAnalyzer:

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

	~HcalIsoTrkSimAnalyzer () override=default

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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
void	analyze (edm::Event const &, edm::EventSetup const &) override

void	beginJob () override

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

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

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

double	eThreshold (const DetId &id, const CaloGeometry *geo) const

std::array< int, 3 >	fillTree (std::vector< math::XYZTLorentzVector > &vecL1, std::vector< math::XYZTLorentzVector > &vecL3, math::XYZPoint &leadPV, std::vector< spr::propagatedGenParticleID > &trackIDs, const CaloGeometry geo, const CaloTopology topo, const HcalTopology theHBHETopology, const EcalChannelStatus theEcalChStatus, const EcalSeverityLevelAlgo theEcalSevlv, edm::Handle< EcalRecHitCollection > &barrelRecHitsHandle, edm::Handle< EcalRecHitCollection > &endcapRecHitsHandle, edm::Handle< HBHERecHitCollection > &hbhe, const edm::Handle< CaloTowerCollection > &towerHandle, edm::Handle< reco::GenParticleCollection > &genParticles, const HcalRespCorrs respCorrs)

DetId	newId (const DetId &)

bool	notaMuon (const reco::GenParticle *pTrack)

double	rhoh (const edm::Handle< CaloTowerCollection > &)

void	storeEnergy (int indx, const HcalRespCorrs respCorrs, const std::vector< DetId > &ids, std::vector< double > &edet, double &eHcal, std::vector< unsigned int > detIds, std::vector< double > *hitEnergies)

Private Attributes
double	a_charIsoR_

double	a_coneR1_

double	a_coneR2_

const double	a_coneR_

const double	a_mipR2_

const double	a_mipR3_

const double	a_mipR4_

const double	a_mipR5_

const double	a_mipR_

const edm::InputTag	algTag_

const bool	collapseDepth_

const int	dataType_

const double	eEcalMax_

const double	eIsolate1_

const double	eIsolate2_

std::vector< double >	etabins_

double	etadist_

double	etahalfdist_

const double	etaMax_

const EcalPFRecHitThresholds *	eThresholds_

const edm::InputTag	extTag_

const bool	fillTreeRange_

const double	hcalScale_

const HcalDDDRecConstants *	hdc_

const bool	hep17_

const double	hitEthrEB_

const double	hitEthrEE0_

const double	hitEthrEE1_

const double	hitEthrEE2_

const double	hitEthrEE3_

const double	hitEthrEEHi_

const double	hitEthrEELo_

HLTConfigProvider	hltConfig_

const bool	ignoreTrigger_

const std::string	l1Filter_

l1t::L1TGlobalUtil *	l1GtUtils_

const std::string	l1TrigName_

const std::string	l2Filter_

const std::string	l3Filter_

const std::string	labelBS_

const std::string	labelEB_

const std::string	labelEE_

const std::string	labelGenTrack_

const std::string	labelHBHE_

const std::string	labelRecVtx_

const std::string	labelTower_

const double	maxRestrictionP_

const int	mode_

const std::string	modnam_

const std::vector< int >	newDepth_

unsigned int	nHigh_

unsigned int	nLow_

unsigned int	nRun_

std::vector< int >	oldDepth_

std::vector< int >	oldDet_

std::vector< int >	oldEta_

const std::vector< int >	oldID_

std::vector< double >	phibins_

double	phidist_

double	phihalfdist_

const std::string	prdnam_

const int	prescaleHigh_

const int	prescaleLow_

const std::string	processName_

const double	ptMin_

const double	pTrackHigh_

const double	pTrackLow_

const double	pTrackMin_

const double	slopeRestrictionP_

int	t_allvertex

int	t_DataType

std::vector< unsigned int > *	t_DetIds

std::vector< unsigned int > *	t_DetIds1

std::vector< unsigned int > *	t_DetIds3

double	t_eAnnular

double	t_eHcal

double	t_eHcal10

double	t_eHcal30

double	t_emaxNearP

double	t_eMipDR

double	t_eMipDR2

double	t_eMipDR3

double	t_eMipDR4

double	t_eMipDR5

int	t_Event

unsigned int	t_EventNo

double	t_EventWeight

double	t_gentrackP

int	t_goodPV

double	t_hAnnular

std::vector< double > *	t_HitEnergies

std::vector< double > *	t_HitEnergies1

std::vector< double > *	t_HitEnergies3

std::vector< bool > *	t_hltbits

double	t_hmaxNearP

int	t_ieta

std::vector< int > *	t_ietaAll

std::vector< int > *	t_ietaGood

int	t_iphi

bool	t_L1Bit

double	t_l1eta

double	t_l1phi

double	t_l1pt

double	t_l3eta

double	t_l3phi

double	t_l3pt

double	t_mindR1

double	t_mindR2

int	t_nTrk

int	t_nVtx

double	t_p

double	t_phi

double	t_pt

bool	t_qltyFlag

bool	t_qltyMissFlag

bool	t_qltyPVFlag

double	t_rhoh

int	t_Run

unsigned int	t_RunNo

bool	t_selectTk

int	t_Tracks

int	t_TracksLoose

int	t_TracksProp

int	t_TracksSaved

int	t_TracksTight

std::vector< int > *	t_trackType

std::vector< bool > *	t_trgbits

bool	t_TrigPass

bool	t_TrigPassSel

const edm::InputTag	theTriggerResultsLabel_

edm::EDGetTokenT< BXVector< GlobalAlgBlk > >	tok_alg_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	tok_bFieldH_

edm::EDGetTokenT< reco::BeamSpot >	tok_bs_

edm::EDGetTokenT< CaloTowerCollection >	tok_cala_

edm::ESGetToken< CaloTopology, CaloTopologyRecord >	tok_caloTopology_

edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord >	tok_ddrec_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd >	tok_ecalChStatus_

edm::ESGetToken< EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd >	tok_ecalPFRecHitThresholds_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT< GenEventInfoProduct >	tok_ew_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	tok_htopo_

edm::EDGetTokenT< reco::GenParticleCollection >	tok_parts_

edm::ESGetToken< HepPDT::ParticleDataTable, PDTRecord >	tok_pdt_

edm::EDGetTokenT< reco::VertexCollection >	tok_recVtx_

edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd >	tok_resp_

edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	tok_sevlv_

edm::EDGetTokenT< trigger::TriggerEvent >	tok_trigEvt_

edm::EDGetTokenT< edm::TriggerResults >	tok_trigRes_

TTree *	tree

TTree *	tree2

const edm::InputTag	triggerEvent_

const std::vector< std::string >	trigNames_

const bool	unCorrect_

const bool	useL1Trigger_

const bool	usePFThresh_

const int	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)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 77 of file HcalIsoTrkSimAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalIsoTrkSimAnalyzer()

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

explicit

Definition at line 200 of file HcalIsoTrkSimAnalyzer.cc.

References a_charIsoR_, a_coneR1_, a_coneR2_, a_coneR_, a_mipR2_, a_mipR3_, a_mipR4_, a_mipR5_, a_mipR_, algTag_, collapseDepth_, edm::EDConsumerBase::consumesCollector(), dataType_, eEcalMax_, eIsolate1_, eIsolate2_, PVValHelper::eta, etabins_, etadist_, etahalfdist_, etaMax_, l1t::Event, extTag_, fillTreeRange_, hcalScale_, hep17_, hitEthrEB_, hitEthrEE0_, hitEthrEE1_, hitEthrEE2_, hitEthrEE3_, hitEthrEEHi_, hitEthrEELo_, mps_fire::i, ignoreTrigger_, ProducerED_cfi::InputTag, isotrackApplyRegressor::k, relval_steps::k2, TFileService::kSharedResource, l1Filter_, l1GtUtils_, l1TrigName_, l2Filter_, l3Filter_, labelBS_, labelEB_, labelEE_, labelGenTrack_, labelHBHE_, labelRecVtx_, labelTower_, M_PI, maxRestrictionP_, mode_, modnam_, newDepth_, oldDepth_, oldDet_, oldEta_, oldID_, phi, phibins_, phidist_, phihalfdist_, prdnam_, prescaleHigh_, prescaleLow_, processName_, ptMin_, pTrackHigh_, pTrackLow_, pTrackMin_, slopeRestrictionP_, theTriggerResultsLabel_, triggerEvent_, trigNames_, unCorrect_, useL1Trigger_, usePFThresh_, and useRaw_.

     : trigNames_(iConfig.getParameter<std::vector<std::string> >("triggers")),
       processName_(iConfig.getParameter<std::string>("processName")),
       l1Filter_(iConfig.getParameter<std::string>("l1Filter")),
       l2Filter_(iConfig.getParameter<std::string>("l2Filter")),
       l3Filter_(iConfig.getParameter<std::string>("l3Filter")),
       ptMin_(iConfig.getParameter<double>("pTMin")),
       etaMax_(iConfig.getParameter<double>("maxChargedHadronEta")),
       a_coneR_(iConfig.getParameter<double>("coneRadius")),
       a_mipR_(iConfig.getParameter<double>("coneRadiusMIP")),
       a_mipR2_(iConfig.getParameter<double>("coneRadiusMIP2")),
       a_mipR3_(iConfig.getParameter<double>("coneRadiusMIP3")),
       a_mipR4_(iConfig.getParameter<double>("coneRadiusMIP4")),
       a_mipR5_(iConfig.getParameter<double>("coneRadiusMIP5")),
       pTrackMin_(iConfig.getParameter<double>("minimumTrackP")),
       eEcalMax_(iConfig.getParameter<double>("maximumEcalEnergy")),
       maxRestrictionP_(iConfig.getParameter<double>("maxTrackP")),
       slopeRestrictionP_(iConfig.getParameter<double>("slopeTrackP")),
       hcalScale_(iConfig.getUntrackedParameter<double>("hHcalScale", 1.0)),
       eIsolate1_(iConfig.getParameter<double>("isolationEnergyTight")),
       eIsolate2_(iConfig.getParameter<double>("isolationEnergyLoose")),
       pTrackLow_(iConfig.getParameter<double>("momentumLow")),
       pTrackHigh_(iConfig.getParameter<double>("momentumHigh")),
       prescaleLow_(iConfig.getParameter<int>("prescaleLow")),
       prescaleHigh_(iConfig.getParameter<int>("prescaleHigh")),
       useRaw_(iConfig.getUntrackedParameter<int>("useRaw", 0)),
       dataType_(iConfig.getUntrackedParameter<int>("dataType", 0)),
       mode_(iConfig.getUntrackedParameter<int>("outMode", 11)),
       ignoreTrigger_(iConfig.getUntrackedParameter<bool>("ignoreTriggers", false)),
       useL1Trigger_(iConfig.getUntrackedParameter<bool>("useL1Trigger", false)),
       unCorrect_(iConfig.getUntrackedParameter<bool>("unCorrect", false)),
       collapseDepth_(iConfig.getUntrackedParameter<bool>("collapseDepth", false)),
       fillTreeRange_(iConfig.getUntrackedParameter<bool>("fillTreeRange", false)),
       hitEthrEB_(iConfig.getParameter<double>("EBHitEnergyThreshold")),
       hitEthrEE0_(iConfig.getParameter<double>("EEHitEnergyThreshold0")),
       hitEthrEE1_(iConfig.getParameter<double>("EEHitEnergyThreshold1")),
       hitEthrEE2_(iConfig.getParameter<double>("EEHitEnergyThreshold2")),
       hitEthrEE3_(iConfig.getParameter<double>("EEHitEnergyThreshold3")),
       hitEthrEELo_(iConfig.getParameter<double>("EEHitEnergyThresholdLow")),
       hitEthrEEHi_(iConfig.getParameter<double>("EEHitEnergyThresholdHigh")),
       triggerEvent_(iConfig.getParameter<edm::InputTag>("labelTriggerEvent")),
       theTriggerResultsLabel_(iConfig.getParameter<edm::InputTag>("labelTriggerResult")),
       labelGenTrack_(iConfig.getParameter<std::string>("labelTrack")),
       labelRecVtx_(iConfig.getParameter<std::string>("labelVertex")),
       labelEB_(iConfig.getParameter<std::string>("labelEBRecHit")),
       labelEE_(iConfig.getParameter<std::string>("labelEERecHit")),
       labelHBHE_(iConfig.getParameter<std::string>("labelHBHERecHit")),
       labelTower_(iConfig.getParameter<std::string>("labelCaloTower")),
       l1TrigName_(iConfig.getUntrackedParameter<std::string>("l1TrigName", "L1_SingleJet60")),
       oldID_(iConfig.getUntrackedParameter<std::vector<int> >("oldID")),
       newDepth_(iConfig.getUntrackedParameter<std::vector<int> >("newDepth")),
       hep17_(iConfig.getUntrackedParameter<bool>("hep17")),
       usePFThresh_(iConfig.getParameter<bool>("usePFThreshold")),
       labelBS_(iConfig.getParameter<std::string>("labelBeamSpot")),
       modnam_(iConfig.getUntrackedParameter<std::string>("moduleName", "")),
       prdnam_(iConfig.getUntrackedParameter<std::string>("producerName", "")),
       algTag_(iConfig.getParameter<edm::InputTag>("algInputTag")),
       extTag_(iConfig.getParameter<edm::InputTag>("extInputTag")),
       tok_trigEvt_(consumes<trigger::TriggerEvent>(triggerEvent_)),
       tok_trigRes_(consumes<edm::TriggerResults>(theTriggerResultsLabel_)),
       tok_bs_(consumes<reco::BeamSpot>(labelBS_)),
       tok_recVtx_((modnam_.empty()) ? consumes<reco::VertexCollection>(labelRecVtx_)
                                     : consumes<reco::VertexCollection>(edm::InputTag(modnam_, labelRecVtx_, prdnam_))),
       tok_EB_((modnam_.empty()) ? consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEB_))
                                 : consumes<EcalRecHitCollection>(edm::InputTag(modnam_, labelEB_, prdnam_))),
       tok_EE_((modnam_.empty()) ? consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEE_))
                                 : consumes<EcalRecHitCollection>(edm::InputTag(modnam_, labelEE_, prdnam_))),
       tok_hbhe_((modnam_.empty()) ? consumes<HBHERecHitCollection>(labelHBHE_)
                                   : consumes<HBHERecHitCollection>(edm::InputTag(modnam_, labelHBHE_, prdnam_))),
       tok_ew_(consumes<GenEventInfoProduct>(edm::InputTag("generator"))),
       tok_parts_(consumes<reco::GenParticleCollection>(edm::InputTag("genParticles"))),
       tok_cala_(consumes<CaloTowerCollection>(labelTower_)),
       tok_alg_(consumes<BXVector<GlobalAlgBlk> >(algTag_)),
       tok_ddrec_(esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>()),
       tok_bFieldH_(esConsumes<MagneticField, IdealMagneticFieldRecord>()),
       tok_ecalChStatus_(esConsumes<EcalChannelStatus, EcalChannelStatusRcd>()),
       tok_sevlv_(esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>()),
       tok_geom_(esConsumes<CaloGeometry, CaloGeometryRecord>()),
       tok_caloTopology_(esConsumes<CaloTopology, CaloTopologyRecord>()),
       tok_htopo_(esConsumes<HcalTopology, HcalRecNumberingRecord>()),
       tok_resp_(esConsumes<HcalRespCorrs, HcalRespCorrsRcd>()),
       tok_pdt_(esConsumes<HepPDT::ParticleDataTable, PDTRecord>()),
       tok_ecalPFRecHitThresholds_(esConsumes<EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd>()),
       nRun_(0),
       nLow_(0),
       nHigh_(0),
       hdc_(nullptr) {
   usesResource(TFileService::kSharedResource);
 
   //now do whatever initialization is needed
   const double isolationRadius(28.9), innerR(10.0), outerR(30.0);
   a_charIsoR_ = a_coneR_ + isolationRadius;
   a_coneR1_ = a_coneR_ + innerR;
   a_coneR2_ = a_coneR_ + outerR;
   // Different isolation cuts are described in DN-2016/029
   // Tight cut uses 2 GeV; Loose cut uses 10 GeV
   // Eta dependent cut uses (maxRestrictionP_ * exp(|ieta|*log(2.5)/18))
   // with the factor for exponential slopeRestrictionP_ = log(2.5)/18
   // maxRestrictionP_ = 8 GeV as came from a study
 
   for (unsigned int k = 0; k < oldID_.size(); ++k) {
     oldDet_.emplace_back((oldID_[k] / 10000) % 10);
     oldEta_.emplace_back((oldID_[k] / 100) % 100);
     oldDepth_.emplace_back(oldID_[k] % 100);
   }
 
   l1GtUtils_ =
       new l1t::L1TGlobalUtil(iConfig, consumesCollector(), *this, algTag_, extTag_, l1t::UseEventSetupIn::Event);
 
   if (modnam_.empty()) {
     edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
                                      << labelBS_ << " " << labelRecVtx_ << " " << labelGenTrack_ << " "
                                      << edm::InputTag("ecalRecHit", labelEB_) << " "
                                      << edm::InputTag("ecalRecHit", labelEE_) << " " << labelHBHE_ << " "
                                      << labelTower_;
   } else {
     edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
                                      << labelBS_ << " " << edm::InputTag(modnam_, labelRecVtx_, prdnam_) << " "
                                      << labelGenTrack_ << " " << edm::InputTag(modnam_, labelEB_, prdnam_) << " "
                                      << edm::InputTag(modnam_, labelEE_, prdnam_) << " "
                                      << edm::InputTag(modnam_, labelHBHE_, prdnam_) << " " << labelTower_;
   }
 
   edm::LogVerbatim("HcalIsoTrack") << "Parameters read from config file \n"
                                    << "\t minPt " << ptMin_ << "\t etaMax " << etaMax_ << "\t a_coneR " << a_coneR_
                                    << ":" << a_coneR1_ << ":" << a_coneR2_ << "\t a_charIsoR " << a_charIsoR_
                                    << "\t a_mipR " << a_mipR_ << "\t a_mipR2 " << a_mipR2_ << "\t a_mipR3 " << a_mipR3_
                                    << "\t a_mipR4 " << a_mipR4_ << "\t a_mipR5 " << a_mipR5_ << "\n pTrackMin_ "
                                    << pTrackMin_ << "\t eEcalMax_ " << eEcalMax_ << "\t maxRestrictionP_ "
                                    << maxRestrictionP_ << "\t slopeRestrictionP_ " << slopeRestrictionP_
                                    << "\t eIsolateStrong_ " << eIsolate1_ << "\t eIsolateSoft_ " << eIsolate2_
                                    << "\t hcalScale_ " << hcalScale_ << "\n\t momentumLow_ " << pTrackLow_
                                    << "\t prescaleLow_ " << prescaleLow_ << "\t momentumHigh_ " << pTrackHigh_
                                    << "\t prescaleHigh_ " << prescaleHigh_ << "\n\t useRaw_ " << useRaw_
                                    << "\t ignoreTrigger_ " << ignoreTrigger_ << "\n\t useL1Trigegr_ " << useL1Trigger_
                                    << "\t dataType_      " << dataType_ << "\t mode_          " << mode_
                                    << "\t unCorrect_     " << unCorrect_ << "\t collapseDepth_ " << collapseDepth_
                                    << "\t fillTreeRange " << fillTreeRange_ << "\t L1TrigName_    " << l1TrigName_
                                    << "\nThreshold flag used " << usePFThresh_ << " value for EB " << hitEthrEB_
                                    << " EE " << hitEthrEE0_ << ":" << hitEthrEE1_ << ":" << hitEthrEE2_ << ":"
                                    << hitEthrEE3_ << ":" << hitEthrEELo_ << ":" << hitEthrEEHi_;
   edm::LogVerbatim("HcalIsoTrack") << "Process " << processName_ << " L1Filter:" << l1Filter_
                                    << " L2Filter:" << l2Filter_ << " L3Filter:" << l3Filter_;
   for (unsigned int k = 0; k < trigNames_.size(); ++k) {
     edm::LogVerbatim("HcalIsoTrack") << "Trigger[" << k << "] " << trigNames_[k];
   }
   edm::LogVerbatim("HcalIsoTrack") << oldID_.size() << " DetIDs to be corrected with HEP17 flag:" << hep17_;
   for (unsigned int k = 0; k < oldID_.size(); ++k)
     edm::LogVerbatim("HcalIsoTrack") << "[" << k << "] Det " << oldDet_[k] << " EtaAbs " << oldEta_[k] << " Depth "
                                      << oldDepth_[k] << ":" << newDepth_[k];
 
   for (int i = 0; i < 10; i++)
     phibins_.push_back(-M_PI + 0.1 * (2 * i + 1) * M_PI);
   for (int i = 0; i < 8; ++i)
     etabins_.push_back(-2.1 + 0.6 * i);
   etadist_ = etabins_[1] - etabins_[0];
   phidist_ = phibins_[1] - phibins_[0];
   etahalfdist_ = 0.5 * etadist_;
   phihalfdist_ = 0.5 * phidist_;
   edm::LogVerbatim("HcalIsoTrack") << "EtaDist " << etadist_ << " " << etahalfdist_ << " PhiDist " << phidist_ << " "
                                    << phihalfdist_;
   unsigned int k1(0), k2(0);
   for (auto phi : phibins_) {
     edm::LogVerbatim("HcalIsoTrack") << "phibin_[" << k1 << "] " << phi;
     ++k1;
   }
   for (auto eta : etabins_) {
     edm::LogVerbatim("HcalIsoTrack") << "etabin_[" << k2 << "] " << eta;
     ++k2;
   }
 }

◆ ~HcalIsoTrkSimAnalyzer()

HcalIsoTrkSimAnalyzer::~HcalIsoTrkSimAnalyzer ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 372 of file HcalIsoTrkSimAnalyzer.cc.

                                                                                        {
   t_Run = iEvent.id().run();
   t_Event = iEvent.id().event();
   t_DataType = dataType_;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "Run " << t_Run << " Event " << t_Event << " type " << t_DataType
                                    << " Luminosity " << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing();
 #endif
   //Get magnetic field and ECAL channel status
   const MagneticField* bField = &iSetup.getData(tok_bFieldH_);
   const EcalChannelStatus* theEcalChStatus = &iSetup.getData(tok_ecalChStatus_);
   const EcalSeverityLevelAlgo* theEcalSevlv = &iSetup.getData(tok_sevlv_);
 
   // get calogeometry and calotopology
   const CaloGeometry* geo = &iSetup.getData(tok_geom_);
   const CaloTopology* caloTopology = &iSetup.getData(tok_caloTopology_);
   const HcalTopology* theHBHETopology = &iSetup.getData(tok_htopo_);
 
   // get response correction
   const HcalRespCorrs* resp = &iSetup.getData(tok_resp_);
   HcalRespCorrs* respCorrs = new HcalRespCorrs(*resp);
   respCorrs->setTopo(theHBHETopology);
 
   // get ECAL thresholds
   eThresholds_ = &iSetup.getData(tok_ecalPFRecHitThresholds_);
 
   // get particle data table
   const HepPDT::ParticleDataTable* pdt = &iSetup.getData(tok_pdt_);
 
   //=== genParticle information
   bool okC(true);
   edm::Handle<reco::GenParticleCollection> genParticles = iEvent.getHandle(tok_parts_);
   if (!genParticles.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the genParticles collection";
     okC = false;
   }
 
   //event weight for FLAT sample
   t_EventWeight = 1.0;
   const edm::Handle<GenEventInfoProduct> genEventInfo = iEvent.getHandle(tok_ew_);
   if (genEventInfo.isValid())
     t_EventWeight = genEventInfo->weight();
 
   //Define the best vertex and the beamspot
   const edm::Handle<reco::VertexCollection> recVtxs = iEvent.getHandle(tok_recVtx_);
   const edm::Handle<reco::BeamSpot> beamSpotH = iEvent.getHandle(tok_bs_);
   math::XYZPoint leadPV(0, 0, 0);
   t_goodPV = t_nVtx = 0;
   if (recVtxs.isValid() && !(recVtxs->empty())) {
     t_nVtx = recVtxs->size();
     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();
   }
   t_allvertex = t_goodPV;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "Primary Vertex " << leadPV << " out of " << t_goodPV << " vertex";
   if (beamSpotH.isValid()) {
     edm::LogVerbatim("HcalIsoTrack") << " Beam Spot " << beamSpotH->position();
   }
 #endif
   // RecHits
   edm::Handle<EcalRecHitCollection> barrelRecHitsHandle = iEvent.getHandle(tok_EB_);
   if (!barrelRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEB_;
     okC = false;
   }
   edm::Handle<EcalRecHitCollection> endcapRecHitsHandle = iEvent.getHandle(tok_EE_);
   if (!endcapRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEE_;
     okC = false;
   }
   edm::Handle<HBHERecHitCollection> hbhe = iEvent.getHandle(tok_hbhe_);
   if (!hbhe.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelHBHE_;
     okC = false;
   }
   const edm::Handle<CaloTowerCollection> caloTower = iEvent.getHandle(tok_cala_);
 
   //Propagate tracks to calorimeter surface)
   std::vector<spr::propagatedGenParticleID> trackIDs =
       spr::propagateCALO(genParticles, pdt, geo, bField, etaMax_, false);
   std::vector<math::XYZTLorentzVector> vecL1, vecL3;
   t_RunNo = iEvent.id().run();
   t_EventNo = iEvent.id().event();
   t_Tracks = (genParticles.product())->size();
   t_TracksProp = trackIDs.size();
   t_ietaAll->clear();
   t_ietaGood->clear();
   t_trackType->clear();
   t_trgbits->clear();
   t_hltbits->clear();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "# of propagated tracks " << t_TracksProp << " out of " << t_Tracks
                                    << " with Trigger " << ignoreTrigger_;
 #endif
 
   //Trigger
   t_trgbits->assign(trigNames_.size(), false);
   t_hltbits->assign(trigNames_.size(), false);
   t_TracksSaved = t_TracksLoose = t_TracksTight = 0;
   t_L1Bit = true;
   t_TrigPass = false;
 
   if (!ignoreTrigger_) {
     //L1
     l1GtUtils_->retrieveL1(iEvent, iSetup, tok_alg_);
     const std::vector<std::pair<std::string, bool> >& finalDecisions = l1GtUtils_->decisionsFinal();
     for (const auto& decision : finalDecisions) {
       if (decision.first.find(l1TrigName_) != std::string::npos) {
         t_L1Bit = decision.second;
         break;
       }
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalIsoTrack") << "Trigger Information for " << l1TrigName_ << " is " << t_L1Bit
                                      << " from a list of " << finalDecisions.size() << " decisions";
 #endif
 
     //HLT
     const edm::Handle<edm::TriggerResults> triggerResults = iEvent.getHandle(tok_trigRes_);
     if (triggerResults.isValid()) {
       const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
       const std::vector<std::string>& names = triggerNames.triggerNames();
       if (!trigNames_.empty()) {
         for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
           int hlt = triggerResults->accept(iHLT);
           for (unsigned int i = 0; i < trigNames_.size(); ++i) {
             if (names[iHLT].find(trigNames_[i]) != std::string::npos) {
               t_trgbits->at(i) = (hlt > 0);
               t_hltbits->at(i) = (hlt > 0);
               if (hlt > 0)
                 t_TrigPass = true;
 #ifdef EDM_ML_DEBUG
               edm::LogVerbatim("HcalIsoTrack")
                   << "This trigger " << names[iHLT] << " Flag " << hlt << ":" << t_trgbits->at(i);
 #endif
             }
           }
         }
       }
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalIsoTrack") << "HLT Information shows " << t_TrigPass << ":" << trigNames_.empty() << ":"
                                      << okC;
 #endif
   }
 
   std::array<int, 3> ntksave{{0, 0, 0}};
   if (ignoreTrigger_ || useL1Trigger_) {
     t_l1pt = t_l1eta = t_l1phi = 0;
     t_l3pt = t_l3eta = t_l3phi = 0;
     if (ignoreTrigger_ || t_L1Bit)
       ntksave = fillTree(vecL1,
                          vecL3,
                          leadPV,
                          trackIDs,
                          geo,
                          caloTopology,
                          theHBHETopology,
                          theEcalChStatus,
                          theEcalSevlv,
                          barrelRecHitsHandle,
                          endcapRecHitsHandle,
                          hbhe,
                          caloTower,
                          genParticles,
                          respCorrs);
     t_TracksSaved = ntksave[0];
     t_TracksLoose = ntksave[1];
     t_TracksTight = ntksave[2];
   } else {
     trigger::TriggerEvent triggerEvent;
     edm::Handle<trigger::TriggerEvent> triggerEventHandle = iEvent.getHandle(tok_trigEvt_);
     if (!triggerEventHandle.isValid()) {
       edm::LogWarning("HcalIsoTrack") << "Error! Can't get the product " << triggerEvent_.label();
     } else if (okC) {
       triggerEvent = *(triggerEventHandle.product());
       const edm::Handle<edm::TriggerResults> triggerResults = iEvent.getHandle(tok_trigRes_);
       const trigger::TriggerObjectCollection& TOC(triggerEvent.getObjects());
       bool done(false);
       if (triggerResults.isValid()) {
         std::vector<std::string> modules;
         const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
         const std::vector<std::string>& names = triggerNames.triggerNames();
         for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
           bool ok = (t_TrigPass) || (trigNames_.empty());
           if (ok) {
             unsigned int triggerindx = hltConfig_.triggerIndex(names[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 EDM_ML_DEBUG
                   edm::LogVerbatim("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_.empty())) {
                       vecL1.push_back(v4);
                     }
 #ifdef EDM_ML_DEBUG
                     edm::LogVerbatim("HcalIsoTrack")
                         << "key " << ifiltrKey << " : pt " << TO.pt() << " eta " << TO.eta() << " phi " << TO.phi()
                         << " mass " << TO.mass() << " Id " << TO.id();
 #endif
                   }
 #ifdef EDM_ML_DEBUG
                   edm::LogVerbatim("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 EDM_ML_DEBUG
               edm::LogVerbatim("HcalIsoTrack") << "lvl2[" << i << "] dR " << dr;
 #endif
               if (dr < mindR1) {
                 mindR1 = dr;
                 mindRvec1 = vecL2[i];
               }
             }
 #ifdef EDM_ML_DEBUG
             edm::LogVerbatim("HcalIsoTrack") << "L2 object closest to L1 " << mindRvec1 << " at Dr " << mindR1;
 #endif
 
             if (!vecL1.empty()) {
               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.empty()) {
               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,
                                  trackIDs,
                                  geo,
                                  caloTopology,
                                  theHBHETopology,
                                  theEcalChStatus,
                                  theEcalSevlv,
                                  barrelRecHitsHandle,
                                  endcapRecHitsHandle,
                                  hbhe,
                                  caloTower,
                                  genParticles,
                                  respCorrs);
               t_TracksSaved += ntksave[0];
               t_TracksLoose += ntksave[1];
               t_TracksTight += ntksave[2];
               done = true;
             }
           }
         }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "Final results on selected tracks " << t_TracksSaved << ":" << t_TracksLoose
                                    << ":" << t_TracksTight;
 #endif
   t_TrigPassSel = (t_TracksSaved > 0);
   tree2->Fill();
 }

◆ beginJob()

void HcalIsoTrkSimAnalyzer::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 672 of file HcalIsoTrkSimAnalyzer.cc.

References compareTotals::fs, t_allvertex, t_DataType, t_DetIds, t_DetIds1, t_DetIds3, t_eAnnular, t_eHcal, t_eHcal10, t_eHcal30, t_emaxNearP, t_eMipDR, t_eMipDR2, t_eMipDR3, t_eMipDR4, t_eMipDR5, t_Event, t_EventNo, t_EventWeight, t_gentrackP, t_goodPV, t_hAnnular, t_HitEnergies, t_HitEnergies1, t_HitEnergies3, t_hltbits, t_hmaxNearP, t_ieta, t_ietaAll, t_ietaGood, t_iphi, t_L1Bit, t_l1eta, t_l1phi, t_l1pt, t_l3eta, t_l3phi, t_l3pt, t_mindR1, t_mindR2, t_nTrk, t_nVtx, t_p, t_phi, t_pt, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_rhoh, t_Run, t_RunNo, t_selectTk, t_Tracks, t_TracksLoose, t_TracksProp, t_TracksSaved, t_TracksTight, t_trackType, t_trgbits, t_TrigPass, t_TrigPassSel, and tree2.

                                      {
   edm::Service<TFileService> fs;
   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_iphi", &t_iphi, "t_iphi/I");
   tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
   tree->Branch("t_nVtx", &t_nVtx, "t_nVtx/I");
   tree->Branch("t_nTrk", &t_nTrk, "t_nTrk/I");
   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_pt", &t_pt, "t_pt/D");
   tree->Branch("t_phi", &t_phi, "t_phi/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_eMipDR2", &t_eMipDR2, "t_eMipDR2/D");
   tree->Branch("t_eMipDR3", &t_eMipDR3, "t_eMipDR3/D");
   tree->Branch("t_eMipDR4", &t_eMipDR4, "t_eMipDR4/D");
   tree->Branch("t_eMipDR5", &t_eMipDR5, "t_eMipDR5/D");
   tree->Branch("t_eHcal", &t_eHcal, "t_eHcal/D");
   tree->Branch("t_eHcal10", &t_eHcal10, "t_eHcal10/D");
   tree->Branch("t_eHcal30", &t_eHcal30, "t_eHcal30/D");
   tree->Branch("t_hmaxNearP", &t_hmaxNearP, "t_hmaxNearP/D");
   tree->Branch("t_emaxNearP", &t_emaxNearP, "t_emaxNearP/D");
   tree->Branch("t_eAnnular", &t_eAnnular, "t_eAnnular/D");
   tree->Branch("t_hAnnular", &t_hAnnular, "t_hAnnular/D");
   tree->Branch("t_rhoh", &t_rhoh, "t_rhoh/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");
   tree->Branch("t_gentrackP", &t_gentrackP, "t_gentrackP/D");
 
   t_DetIds = new std::vector<unsigned int>();
   t_DetIds1 = new std::vector<unsigned int>();
   t_DetIds3 = new std::vector<unsigned int>();
   t_HitEnergies = new std::vector<double>();
   t_HitEnergies1 = new std::vector<double>();
   t_HitEnergies3 = 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);
   tree->Branch("t_DetIds1", "std::vector<unsigned int>", &t_DetIds1);
   tree->Branch("t_DetIds3", "std::vector<unsigned int>", &t_DetIds3);
   tree->Branch("t_HitEnergies1", "std::vector<double>", &t_HitEnergies1);
   tree->Branch("t_HitEnergies3", "std::vector<double>", &t_HitEnergies3);
 
   tree2 = fs->make<TTree>("EventInfo", "Event Information");
 
   tree2->Branch("t_RunNo", &t_RunNo, "t_RunNo/i");
   tree2->Branch("t_EventNo", &t_EventNo, "t_EventNo/i");
   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");
   tree2->Branch("t_TracksLoose", &t_TracksLoose, "t_TracksLoose/I");
   tree2->Branch("t_TracksTight", &t_TracksTight, "t_TracksTight/I");
   tree2->Branch("t_TrigPass", &t_TrigPass, "t_TrigPass/O");
   tree2->Branch("t_TrigPassSel", &t_TrigPassSel, "t_TrigPassSel/O");
   tree2->Branch("t_L1Bit", &t_L1Bit, "t_L1Bit/O");
   tree2->Branch("t_allvertex", &t_allvertex, "t_allvertex/I");
   t_hltbits = new std::vector<bool>();
   t_ietaAll = new std::vector<int>();
   t_ietaGood = new std::vector<int>();
   t_trackType = new std::vector<int>();
   tree2->Branch("t_ietaAll", "std::vector<int>", &t_ietaAll);
   tree2->Branch("t_ietaGood", "std::vector<int>", &t_ietaGood);
   tree2->Branch("t_trackType", "std::vector<int>", &t_trackType);
   tree2->Branch("t_hltbits", "std::vector<bool>", &t_hltbits);
 }

◆ beginRun()

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

overrideprivate

Definition at line 754 of file HcalIsoTrkSimAnalyzer.cc.

References RemoveAddSevLevel::flag, edm::EventSetup::getData(), hdc_, hltConfig_, HLTConfigProvider::init(), create_idmaps::n, nRun_, processName_, edm::RunBase::run(), HLTConfigProvider::size(), tok_ddrec_, HLTConfigProvider::triggerIndex(), and trigNames_.

                                                                                     {
   hdc_ = &iSetup.getData(tok_ddrec_);
 
   bool changed_(true);
   bool flag = hltConfig_.init(iRun, iSetup, processName_, changed_);
   edm::LogVerbatim("HcalIsoTrack") << "Run[" << nRun_ << "] " << iRun.run() << " process " << processName_
                                    << " init flag " << flag << " change flag " << changed_;
   // check if trigger names in (new) config
   if (changed_) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("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 EDM_ML_DEBUG
       } else {
         edm::LogVerbatim("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " exists";
 #endif
       }
     }
   }
 }

◆ dR()

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

private

Definition at line 1223 of file HcalIsoTrkSimAnalyzer.cc.

References reco::deltaR().

Referenced by analyze(), and fillTree().

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

◆ endRun()

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

overrideprivate

Definition at line 782 of file HcalIsoTrkSimAnalyzer.cc.

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

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

◆ eThreshold()

double HcalIsoTrkSimAnalyzer::eThreshold	(	const DetId &	id,
		const CaloGeometry *	geo
	)		const

private

Definition at line 1258 of file HcalIsoTrkSimAnalyzer.cc.

References funct::abs(), EcalBarrel, PVValHelper::eta, CaloGeometry::getPosition(), hitEthrEB_, hitEthrEE0_, hitEthrEE1_, hitEthrEE2_, hitEthrEE3_, hitEthrEEHi_, hitEthrEELo_, EcalPhiSymFlatTableProducers_cfi::id, and usePFThresh_.

Referenced by fillTree().

                                                                                        {
   double eThr(hitEthrEB_);
   if (usePFThresh_) {
     eThr = static_cast<double>((*eThresholds_)[id]);
   } else {
     const GlobalPoint& pos = geo->getPosition(id);
     double eta = std::abs(pos.eta());
     if (id.subdetId() != EcalBarrel) {
       eThr = (((eta * hitEthrEE3_ + hitEthrEE2_) * eta + hitEthrEE1_) * eta + hitEthrEE0_);
       if (eThr < hitEthrEELo_)
         eThr = hitEthrEELo_;
       else if (eThr > hitEthrEEHi_)
         eThr = hitEthrEEHi_;
     }
   }
   return eThr;
 }

◆ fillDescriptions()

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

static

Definition at line 787 of file HcalIsoTrkSimAnalyzer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, AlCaHLTBitMon_QueryRunRegistry::string, and ClusterTask_cfi::trig.

                                                                                        {
   edm::ParameterSetDescription desc;
   std::vector<std::string> trig = {"HLT_PFJet40",
                                    "HLT_PFJet60",
                                    "HLT_PFJet80",
                                    "HLT_PFJet140",
                                    "HLT_PFJet200",
                                    "HLT_PFJet260",
                                    "HLT_PFJet320",
                                    "HLT_PFJet400",
                                    "HLT_PFJet450",
                                    "HLT_PFJet500"};
   desc.add<std::vector<std::string> >("triggers", trig);
   desc.add<std::string>("processName", "HLT");
   desc.add<std::string>("l1Filter", "");
   desc.add<std::string>("l2Filter", "L2Filter");
   desc.add<std::string>("l3Filter", "Filter");
   desc.add<double>("pTMin", 1.0), desc.add<double>("maxChargedHadronEta", 3.0);
   // Minimum momentum of selected isolated track and signal zone
   desc.add<double>("minimumTrackP", 10.0);
   desc.add<double>("coneRadius", 34.98);
   // signal zone in ECAL and MIP energy cutoff
   desc.add<double>("coneRadiusMIP", 14.0);
   desc.add<double>("coneRadiusMIP2", 18.0);
   desc.add<double>("coneRadiusMIP3", 20.0);
   desc.add<double>("coneRadiusMIP4", 22.0);
   desc.add<double>("coneRadiusMIP5", 24.0);
   desc.add<double>("maximumEcalEnergy", 2.0);
   // following 4 parameters are for isolation cuts and described in the code
   desc.add<double>("maxTrackP", 8.0);
   desc.add<double>("slopeTrackP", 0.05090504066);
   desc.add<double>("isolationEnergyTight", 2.0);
   desc.add<double>("isolationEnergyLoose", 10.0);
   // energy thershold for ECAL (from Egamma group)
   desc.add<double>("EBHitEnergyThreshold", 0.08);
   desc.add<double>("EEHitEnergyThreshold0", 0.30);
   desc.add<double>("EEHitEnergyThreshold1", 0.00);
   desc.add<double>("EEHitEnergyThreshold2", 0.00);
   desc.add<double>("EEHitEnergyThreshold3", 0.00);
   desc.add<double>("EEHitEnergyThresholdLow", 0.30);
   desc.add<double>("EEHitEnergyThresholdHigh", 0.30);
   // prescale factors
   desc.add<double>("momentumLow", 40.0);
   desc.add<double>("momentumHigh", 60.0);
   desc.add<int>("prescaleLow", 1);
   desc.add<int>("prescaleHigh", 1);
   // various labels for collections used in the code
   desc.add<edm::InputTag>("labelTriggerEvent", edm::InputTag("hltTriggerSummaryAOD", "", "HLT"));
   desc.add<edm::InputTag>("labelTriggerResult", edm::InputTag("TriggerResults", "", "HLT"));
   desc.add<std::string>("labelTrack", "generalTracks");
   desc.add<std::string>("labelVertex", "offlinePrimaryVertices");
   desc.add<std::string>("labelEBRecHit", "EcalRecHitsEB");
   desc.add<std::string>("labelEERecHit", "EcalRecHitsEE");
   desc.add<std::string>("labelHBHERecHit", "hbhereco");
   desc.add<std::string>("labelBeamSpot", "offlineBeamSpot");
   desc.add<std::string>("labelCaloTower", "towerMaker");
   desc.add<edm::InputTag>("algInputTag", edm::InputTag("gtStage2Digis"));
   desc.add<edm::InputTag>("extInputTag", edm::InputTag("gtStage2Digis"));
   desc.addUntracked<std::string>("moduleName", "");
   desc.addUntracked<std::string>("producerName", "");
   //  Various flags used for selecting tracks, choice of energy Method2/0
   //  Data type 0/1 for single jet trigger or others
   desc.addUntracked<int>("useRaw", 0);
   desc.addUntracked<bool>("ignoreTriggers", false);
   desc.addUntracked<bool>("useL1Trigger", false);
   desc.addUntracked<double>("hcalScale", 1.0);
   desc.addUntracked<int>("dataType", 0);
   desc.addUntracked<bool>("unCorrect", false);
   desc.addUntracked<bool>("collapseDepth", false);
   desc.addUntracked<bool>("fillTreeRange", false);
   desc.addUntracked<std::string>("l1TrigName", "L1_SingleJet60");
   desc.addUntracked<int>("outMode", 11);
   std::vector<int> dummy;
   desc.addUntracked<std::vector<int> >("oldID", dummy);
   desc.addUntracked<std::vector<int> >("newDepth", dummy);
   desc.addUntracked<bool>("hep17", false);
   desc.add<bool>("usePFThreshold", true);
   descriptions.add("hcalIsoTrkSimAnalyzer", desc);
 }

◆ fillTree()

std::array< int, 3 > HcalIsoTrkSimAnalyzer::fillTree	(	std::vector< math::XYZTLorentzVector > &	vecL1,
		std::vector< math::XYZTLorentzVector > &	vecL3,
		math::XYZPoint &	leadPV,
		std::vector< spr::propagatedGenParticleID > &	trackIDs,
		const CaloGeometry *	geo,
		const CaloTopology *	topo,
		const HcalTopology *	theHBHETopology,
		const EcalChannelStatus *	theEcalChStatus,
		const EcalSeverityLevelAlgo *	theEcalSevlv,
		edm::Handle< EcalRecHitCollection > &	barrelRecHitsHandle,
		edm::Handle< EcalRecHitCollection > &	endcapRecHitsHandle,
		edm::Handle< HBHERecHitCollection > &	hbhe,
		const edm::Handle< CaloTowerCollection > &	towerHandle,
		edm::Handle< reco::GenParticleCollection > &	genParticles,
		const HcalRespCorrs *	respCorrs
	)

private

Definition at line 867 of file HcalIsoTrkSimAnalyzer.cc.

References a_charIsoR_, a_coneR1_, a_coneR2_, a_coneR_, a_mipR2_, a_mipR3_, a_mipR4_, a_mipR5_, a_mipR_, funct::abs(), accept(), spr::chargeIsolationGenCone(), spr::chargeIsolationGenEcal(), hcalRecHitTable_cff::detId, l1ctLayer1_cff::dr, dR(), spr::eCone_ecal(), spr::eCone_hcal(), spr::eECALmatrix(), eEcalMax_, spr::eHCALmatrix(), eIsolate1_, eIsolate2_, reco::LeafCandidate::eta(), eThreshold(), JetChargeProducer_cfi::exp, fillTreeRange_, photonIsolationHIProducer_cfi::hbhe, testProducerWithPsetDescEmpty_cfi::i3, isotrackApplyRegressor::k, maxRestrictionP_, newId(), nHigh_, jetSelector_cfi::nLoose, nLow_, notaMuon(), beamSpotPI::nTracks, oldID_, reco::LeafCandidate::p(), reco::LeafCandidate::phi(), prescaleHigh_, prescaleLow_, reco::LeafCandidate::pt(), pTrackHigh_, pTrackLow_, pTrackMin_, reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), rhoh(), slopeRestrictionP_, storeEnergy(), t_DetIds, t_DetIds1, t_DetIds3, t_eAnnular, t_eHcal, t_eHcal10, t_eHcal30, t_emaxNearP, t_eMipDR, t_eMipDR2, t_eMipDR3, t_eMipDR4, t_eMipDR5, t_gentrackP, t_hAnnular, t_HitEnergies, t_HitEnergies1, t_HitEnergies3, t_hmaxNearP, t_ieta, t_ietaAll, t_ietaGood, t_iphi, t_mindR1, t_mindR2, t_nTrk, t_p, t_phi, t_pt, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_rhoh, t_selectTk, t_trackType, t_trgbits, l1tHGCalTowerProducer_cfi::tower, and useRaw_.

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

                                                                                    {
   int nSave(0), nLoose(0), nTight(0);
   //Loop over tracks
   unsigned int nTracks(0);
   t_nTrk = trackIDs.size();
   t_rhoh = (tower.isValid()) ? rhoh(tower) : 0;
   for (auto const& trkDetItr : trackIDs) {
     const reco::GenParticle* pTrack = &(*(trkDetItr.trkItr));
     math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("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.empty()) ? dR(vecL1[0], v4) : 999;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalIsoTrack") << "Closest L3 object at dr :" << t_mindR2 << " and from L1 " << t_mindR1;
 #endif
     t_ieta = t_iphi = 0;
     if (trkDetItr.okHCAL) {
       HcalDetId detId = (HcalDetId)(trkDetItr.detIdHCAL);
       t_ieta = detId.ieta();
       t_iphi = detId.iphi();
       if (t_p > 40.0 && t_p <= 60.0)
         t_ietaAll->emplace_back(t_ieta);
     }
     //Selection of good track
     t_selectTk = t_qltyMissFlag = t_qltyPVFlag = true;
     double eIsolation = maxRestrictionP_ * exp(slopeRestrictionP_ * std::abs((double)t_ieta));
     if (eIsolation < eIsolate1_)
       eIsolation = eIsolate1_;
     if (eIsolation < eIsolate2_)
       eIsolation = eIsolate2_;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " << t_qltyFlag << "|" << trkDetItr.okECAL << "|"
                                      << trkDetItr.okHCAL << " eIsolation " << eIsolation;
 #endif
     t_qltyFlag = (t_selectTk && trkDetItr.okECAL && trkDetItr.okHCAL);
     bool notMuon = notaMuon(pTrack);
     if (t_qltyFlag && notMuon) {
       int nNearTRKs(0);
       std::vector<DetId> eIds;
       std::vector<double> eHit;
       t_eMipDR = spr::eCone_ecal(geo,
                                  barrelRecHitsHandle,
                                  endcapRecHitsHandle,
                                  trkDetItr.pointHCAL,
                                  trkDetItr.pointECAL,
                                  a_mipR_,
                                  trkDetItr.directionECAL,
                                  eIds,
                                  eHit);
       double eEcal(0);
       for (unsigned int k = 0; k < eIds.size(); ++k) {
         if (eHit[k] > eThreshold(eIds[k], geo))
           eEcal += eHit[k];
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR << ":" << eEcal;
 #endif
       t_eMipDR = eEcal;
       std::vector<DetId> eIds2;
       std::vector<double> eHit2;
       t_eMipDR2 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr.pointHCAL,
                                   trkDetItr.pointECAL,
                                   a_mipR2_,
                                   trkDetItr.directionECAL,
                                   eIds2,
                                   eHit2);
       double eEcal2(0);
       for (unsigned int k = 0; k < eIds2.size(); ++k) {
         if (eHit2[k] > eThreshold(eIds2[k], geo))
           eEcal2 += eHit2[k];
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR2 << ":" << eEcal2;
 #endif
       t_eMipDR2 = eEcal2;
       std::vector<DetId> eIds3;
       std::vector<double> eHit3;
       t_eMipDR3 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr.pointHCAL,
                                   trkDetItr.pointECAL,
                                   a_mipR3_,
                                   trkDetItr.directionECAL,
                                   eIds3,
                                   eHit3);
       double eEcal3(0);
       for (unsigned int k = 0; k < eIds3.size(); ++k) {
         if (eHit3[k] > eThreshold(eIds3[k], geo))
           eEcal3 += eHit3[k];
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR3 << ":" << eEcal3;
 #endif
       t_eMipDR3 = eEcal3;
       std::vector<DetId> eIds4;
       std::vector<double> eHit4;
       t_eMipDR4 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr.pointHCAL,
                                   trkDetItr.pointECAL,
                                   a_mipR4_,
                                   trkDetItr.directionECAL,
                                   eIds4,
                                   eHit4);
       double eEcal4(0);
       for (unsigned int k = 0; k < eIds4.size(); ++k) {
         if (eHit4[k] > eThreshold(eIds4[k], geo))
           eEcal4 += eHit4[k];
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR4 << ":" << eEcal4;
 #endif
       t_eMipDR4 = eEcal4;
       std::vector<DetId> eIds5;
       std::vector<double> eHit5;
       t_eMipDR5 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr.pointHCAL,
                                   trkDetItr.pointECAL,
                                   a_mipR5_,
                                   trkDetItr.directionECAL,
                                   eIds5,
                                   eHit5);
       double eEcal5(0);
       for (unsigned int k = 0; k < eIds5.size(); ++k) {
         if (eHit5[k] > eThreshold(eIds5[k], geo))
           eEcal5 += eHit5[k];
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR5 << ":" << eEcal5;
 #endif
       t_eMipDR5 = eEcal5;
 
       t_emaxNearP = spr::chargeIsolationGenEcal(nTracks, trackIDs, geo, caloTopology, 15, 15);
       const DetId cellE(trkDetItr.detIdECAL);
       std::pair<double, bool> e11x11P = spr::eECALmatrix(cellE,
                                                          barrelRecHitsHandle,
                                                          endcapRecHitsHandle,
                                                          *theEcalChStatus,
                                                          geo,
                                                          caloTopology,
                                                          theEcalSevlv,
                                                          5,
                                                          5,
                                                          -100.0,
                                                          -100.0,
                                                          -100.0,
                                                          100.0);
       std::pair<double, bool> e15x15P = spr::eECALmatrix(cellE,
                                                          barrelRecHitsHandle,
                                                          endcapRecHitsHandle,
                                                          *theEcalChStatus,
                                                          geo,
                                                          caloTopology,
                                                          theEcalSevlv,
                                                          7,
                                                          7,
                                                          -100.0,
                                                          -100.0,
                                                          -100.0,
                                                          100.0);
       if (e11x11P.second && e15x15P.second) {
         t_eAnnular = (e15x15P.first - e11x11P.first);
       } else {
         t_eAnnular = -(e15x15P.first - e11x11P.first);
       }
       t_hmaxNearP = spr::chargeIsolationGenCone(nTracks, trackIDs, a_charIsoR_, nNearTRKs, false);
       const DetId cellH(trkDetItr.detIdHCAL);
       double h5x5 = spr::eHCALmatrix(
           theHBHETopology, cellH, hbhe, 2, 2, false, true, -100.0, -100.0, -100.0, -100.0, -100.0, 100.0);
       double h7x7 = spr::eHCALmatrix(
           theHBHETopology, cellH, hbhe, 3, 3, false, true, -100.0, -100.0, -100.0, -100.0, -100.0, 100.0);
       t_hAnnular = h7x7 - h5x5;
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "max p Near (Ecal) " << t_emaxNearP << " (Hcal) " << t_hmaxNearP
                                        << " Annular E (Ecal) " << e11x11P.first << ":" << e15x15P.first << ":"
                                        << t_eAnnular << " (Hcal) " << h5x5 << ":" << h7x7 << ":" << t_hAnnular;
 #endif
       t_gentrackP = pTrack->p();
       if (t_eMipDR < eEcalMax_ && t_hmaxNearP < eIsolation) {
         t_DetIds->clear();
         t_HitEnergies->clear();
         t_DetIds1->clear();
         t_HitEnergies1->clear();
         t_DetIds3->clear();
         t_HitEnergies3->clear();
         int nRecHits(-999), nRecHits1(-999), nRecHits3(-999);
         std::vector<DetId> ids, ids1, ids3;
         std::vector<double> edet0, edet1, edet3;
         t_eHcal = spr::eCone_hcal(geo,
                                   hbhe,
                                   trkDetItr.pointHCAL,
                                   trkDetItr.pointECAL,
                                   a_coneR_,
                                   trkDetItr.directionHCAL,
                                   nRecHits,
                                   ids,
                                   edet0,
                                   useRaw_);
         if (!oldID_.empty()) {
           for (unsigned k = 0; k < ids.size(); ++k)
             ids[k] = newId(ids[k]);
         }
         storeEnergy(0, respCorrs, ids, edet0, t_eHcal, t_DetIds, t_HitEnergies);
 
         //----- hcal energy in the extended cone 1 (a_coneR+10) --------------
         t_eHcal10 = spr::eCone_hcal(geo,
                                     hbhe,
                                     trkDetItr.pointHCAL,
                                     trkDetItr.pointECAL,
                                     a_coneR1_,
                                     trkDetItr.directionHCAL,
                                     nRecHits1,
                                     ids1,
                                     edet1,
                                     useRaw_);
         if (!oldID_.empty()) {
           for (unsigned k = 0; k < ids1.size(); ++k)
             ids1[k] = newId(ids1[k]);
         }
         storeEnergy(1, respCorrs, ids1, edet1, t_eHcal10, t_DetIds1, t_HitEnergies1);
 
         //----- hcal energy in the extended cone 3 (a_coneR+30) --------------
         t_eHcal30 = spr::eCone_hcal(geo,
                                     hbhe,
                                     trkDetItr.pointHCAL,
                                     trkDetItr.pointECAL,
                                     a_coneR2_,
                                     trkDetItr.directionHCAL,
                                     nRecHits3,
                                     ids3,
                                     edet3,
                                     useRaw_);
         if (!oldID_.empty()) {
           for (unsigned k = 0; k < ids3.size(); ++k)
             ids3[k] = newId(ids3[k]);
         }
         storeEnergy(3, respCorrs, ids3, edet3, t_eHcal30, t_DetIds3, t_HitEnergies3);
 
         t_p = pTrack->p();
         t_pt = pTrack->pt();
         t_phi = pTrack->phi();
 
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) :" << t_pt << "|"
                                          << pTrack->eta() << "|" << t_phi << "|" << t_p << " Generator Level p "
                                          << t_gentrackP;
         edm::LogVerbatim("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 ll = 0; ll < t_DetIds->size(); ll++) {
           edm::LogVerbatim("HcalIsoTrack")
               << "det id is = " << t_DetIds->at(ll) << "   hit enery is  = " << t_HitEnergies->at(ll);
         }
         for (unsigned int ll = 0; ll < t_DetIds1->size(); ll++) {
           edm::LogVerbatim("HcalIsoTrack")
               << "det id is = " << t_DetIds1->at(ll) << "   hit enery is  = " << t_HitEnergies1->at(ll);
         }
         for (unsigned int ll = 0; ll < t_DetIds3->size(); ll++) {
           edm::LogVerbatim("HcalIsoTrack")
               << "det id is = " << t_DetIds3->at(ll) << "   hit enery is  = " << t_HitEnergies3->at(ll);
         }
 #endif
         bool accept(false);
         if (t_p > pTrackMin_) {
           if (t_p < pTrackLow_) {
             ++nLow_;
             if (prescaleLow_ <= 1)
               accept = true;
             else if (nLow_ % prescaleLow_ == 1)
               accept = true;
           } else if (t_p > pTrackHigh_) {
             ++nHigh_;
             if (prescaleHigh_ <= 1)
               accept = true;
             else if (nHigh_ % prescaleHigh_ == 1)
               accept = true;
           } else {
             accept = true;
           }
         }
         if (fillTreeRange_) {
           if ((t_p < pTrackLow_) || (t_p > pTrackHigh_))
             accept = false;
         }
         if (accept) {
           tree->Fill();
           nSave++;
           int type(0);
           if (t_eMipDR < 1.0) {
             if (t_hmaxNearP < eIsolate2_) {
               ++nLoose;
               type = 1;
             }
             if (t_hmaxNearP < eIsolate1_) {
               ++nTight;
               type = 2;
             }
           }
           if (t_p > 40.0 && t_p <= 60.0 && t_selectTk) {
             t_ietaGood->emplace_back(t_ieta);
             t_trackType->emplace_back(type);
           }
 #ifdef EDM_ML_DEBUG
           for (unsigned int k = 0; k < t_trgbits->size(); k++) {
             edm::LogVerbatim("HcalIsoTrack") << "trigger bit is  = " << t_trgbits->at(k);
           }
 #endif
         }
       }
     }
     ++nTracks;
   }
   std::array<int, 3> i3{{nSave, nLoose, nTight}};
   return i3;
 }

◆ newId()

DetId HcalIsoTrkSimAnalyzer::newId ( const DetId & id )

private

Definition at line 1276 of file HcalIsoTrkSimAnalyzer.cc.

References HcalDetId::depth(), hep17_, EcalPhiSymFlatTableProducers_cfi::id, HcalDetId::ieta(), HcalDetId::ietaAbs(), HcalDetId::iphi(), isotrackApplyRegressor::k, newDepth_, oldDepth_, oldDet_, oldEta_, oldID_, HcalDetId::subdet(), DetId::subdetId(), and HcalDetId::zside().

Referenced by fillTree().

                                                   {
   HcalDetId hid(id);
   if (hep17_ && ((hid.iphi() < 63) || (hid.iphi() > 66) || (hid.zside() < 0)))
     return id;
   for (unsigned int k = 0; k < oldID_.size(); ++k) {
     if ((hid.subdetId() == oldDet_[k]) && (hid.ietaAbs() == oldEta_[k]) && (hid.depth() == oldDepth_[k])) {
       return static_cast<DetId>(HcalDetId(hid.subdet(), hid.ieta(), hid.iphi(), newDepth_[k]));
     }
   }
   return id;
 }

◆ notaMuon()

bool HcalIsoTrkSimAnalyzer::notaMuon ( const reco::GenParticle * pTrack )

private

Definition at line 1347 of file HcalIsoTrkSimAnalyzer.cc.

References reco::LeafCandidate::charge(), RemoveAddSevLevel::flag, and reco::LeafCandidate::pdgId().

Referenced by fillTree().

                                                                   {
   int id = pTrack->pdgId();
   bool flag = ((id != 13) && (id != -13) && (pTrack->charge() != 0));
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "notaMuon: ID " << id << " charge " << pTrack->charge() << " Flag " << flag;
 #endif
   return flag;
 }

◆ rhoh()

double HcalIsoTrkSimAnalyzer::rhoh ( const edm::Handle< CaloTowerCollection > & tower )

private

Definition at line 1227 of file HcalIsoTrkSimAnalyzer.cc.

References reco::deltaPhi(), PVValHelper::eta, etabins_, etadist_, etahalfdist_, phi, phibins_, phidist_, phihalfdist_, jetUpdater_cfi::sort, and l1tHGCalTowerProducer_cfi::tower.

Referenced by fillTree().

                                                                               {
   std::vector<double> sumPFNallSMDQH2;
   sumPFNallSMDQH2.reserve(phibins_.size() * etabins_.size());
 
   for (auto eta : etabins_) {
     for (auto phi : phibins_) {
       double hadder = 0;
       for (const auto& pf_it : (*tower)) {
         if (fabs(eta - pf_it.eta()) > etahalfdist_)
           continue;
         if (fabs(reco::deltaPhi(phi, pf_it.phi())) > phihalfdist_)
           continue;
         hadder += pf_it.hadEt();
       }
       sumPFNallSMDQH2.emplace_back(hadder);
     }
   }
 
   double evt_smdq(0);
   std::sort(sumPFNallSMDQH2.begin(), sumPFNallSMDQH2.end());
   if (sumPFNallSMDQH2.size() % 2)
     evt_smdq = sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 1) / 2];
   else
     evt_smdq = (sumPFNallSMDQH2[sumPFNallSMDQH2.size() / 2] + sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 2) / 2]) / 2.;
   double rhoh = evt_smdq / (etadist_ * phidist_);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "Rho " << evt_smdq << ":" << rhoh;
 #endif
   return rhoh;
 }

◆ storeEnergy()

void HcalIsoTrkSimAnalyzer::storeEnergy	(	int	indx,
		const HcalRespCorrs *	respCorrs,
		const std::vector< DetId > &	ids,
		std::vector< double > &	edet,
		double &	eHcal,
		std::vector< unsigned int > *	detIds,
		std::vector< double > *	hitEnergies
	)

private

Definition at line 1288 of file HcalIsoTrkSimAnalyzer.cc.

References funct::abs(), collapseDepth_, alignCSCRings::corr, reco::JetExtendedAssociation::getValue(), HcalCondObjectContainer< Item >::getValues(), hcalScale_, hdc_, EcalPhiSymFlatTableProducers_cfi::id, isotrackApplyRegressor::k, HcalDDDRecConstants::mergedDepthDetId(), nano_mu_digi_cff::rawId, and unCorrect_.

Referenced by fillTree().

                                                                         {
   double ehcal(0);
   if (unCorrect_) {
     for (unsigned int k = 0; k < ids.size(); ++k) {
       double corr = (respCorrs->getValues(ids[k]))->getValue();
       if (corr != 0)
         edet[k] /= corr;
       ehcal += edet[k];
     }
   } else {
     for (const auto& en : edet)
       ehcal += en;
   }
   if ((std::abs(ehcal - eHcal) > 0.001) && (!unCorrect_))
     edm::LogWarning("HcalIsoTrack") << "Check inconsistent energies: " << indx << " " << eHcal << ":" << ehcal
                                     << " from " << ids.size() << " cells";
   eHcal = hcalScale_ * ehcal;
 
   if (collapseDepth_) {
     std::map<HcalDetId, double> hitMap;
     for (unsigned int k = 0; k < ids.size(); ++k) {
       HcalDetId id = hdc_->mergedDepthDetId(HcalDetId(ids[k]));
       auto itr = hitMap.find(id);
       if (itr == hitMap.end()) {
         hitMap[id] = edet[k];
       } else {
         (itr->second) += edet[k];
       }
     }
     detIds->reserve(hitMap.size());
     hitEnergies->reserve(hitMap.size());
     for (const auto& hit : hitMap) {
       detIds->emplace_back(hit.first.rawId());
       hitEnergies->emplace_back(hit.second);
     }
   } else {
     detIds->reserve(ids.size());
     hitEnergies->reserve(ids.size());
     for (unsigned int k = 0; k < ids.size(); ++k) {
       detIds->emplace_back(ids[k].rawId());
       hitEnergies->emplace_back(edet[k]);
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "Input to storeEnergy with " << ids.size() << " cells";
   for (unsigned int k = 0; k < ids.size(); ++k)
     edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId(ids[k]) << " E " << edet[k];
   edm::LogVerbatim("HcalIsoTrack") << "Output of storeEnergy with " << detIds->size() << " cells and Etot " << eHcal;
   for (unsigned int k = 0; k < detIds->size(); ++k)
     edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId((*detIds)[k]) << " E " << (*hitEnergies)[k];
 #endif
 }

Member Data Documentation

◆ a_charIsoR_

double HcalIsoTrkSimAnalyzer::a_charIsoR_

private

Definition at line 169 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_coneR1_

double HcalIsoTrkSimAnalyzer::a_coneR1_

private

Definition at line 169 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_coneR2_

double HcalIsoTrkSimAnalyzer::a_coneR2_

private

Definition at line 169 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_coneR_

const double HcalIsoTrkSimAnalyzer::a_coneR_

private

Definition at line 124 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_mipR2_

const double HcalIsoTrkSimAnalyzer::a_mipR2_

private

Definition at line 124 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_mipR3_

const double HcalIsoTrkSimAnalyzer::a_mipR3_

private

Definition at line 124 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_mipR4_

const double HcalIsoTrkSimAnalyzer::a_mipR4_

private

Definition at line 125 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_mipR5_

const double HcalIsoTrkSimAnalyzer::a_mipR5_

private

Definition at line 125 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_mipR_

const double HcalIsoTrkSimAnalyzer::a_mipR_

private

Definition at line 124 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ algTag_

const edm::InputTag HcalIsoTrkSimAnalyzer::algTag_

private

Definition at line 143 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ collapseDepth_

const bool HcalIsoTrkSimAnalyzer::collapseDepth_

private

Definition at line 132 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and storeEnergy().

◆ dataType_

const int HcalIsoTrkSimAnalyzer::dataType_

private

Definition at line 130 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ eEcalMax_

const double HcalIsoTrkSimAnalyzer::eEcalMax_

private

Definition at line 125 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ eIsolate1_

const double HcalIsoTrkSimAnalyzer::eIsolate1_

private

Definition at line 127 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ eIsolate2_

const double HcalIsoTrkSimAnalyzer::eIsolate2_

private

Definition at line 127 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ etabins_

std::vector<double> HcalIsoTrkSimAnalyzer::etabins_

private

Definition at line 173 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ etadist_

double HcalIsoTrkSimAnalyzer::etadist_

private

Definition at line 175 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ etahalfdist_

double HcalIsoTrkSimAnalyzer::etahalfdist_

private

Definition at line 175 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ etaMax_

const double HcalIsoTrkSimAnalyzer::etaMax_

private

Definition at line 123 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ eThresholds_

const EcalPFRecHitThresholds* HcalIsoTrkSimAnalyzer::eThresholds_

private

Definition at line 171 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ extTag_

const edm::InputTag HcalIsoTrkSimAnalyzer::extTag_

private

Definition at line 143 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ fillTreeRange_

const bool HcalIsoTrkSimAnalyzer::fillTreeRange_

private

Definition at line 132 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ hcalScale_

const double HcalIsoTrkSimAnalyzer::hcalScale_

private

Definition at line 127 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and storeEnergy().

◆ hdc_

const HcalDDDRecConstants* HcalIsoTrkSimAnalyzer::hdc_

private

Definition at line 170 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginRun(), and storeEnergy().

◆ hep17_

const bool HcalIsoTrkSimAnalyzer::hep17_

private

Definition at line 140 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and newId().

◆ hitEthrEB_

const double HcalIsoTrkSimAnalyzer::hitEthrEB_

private

Definition at line 133 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEE0_

const double HcalIsoTrkSimAnalyzer::hitEthrEE0_

private

Definition at line 133 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEE1_

const double HcalIsoTrkSimAnalyzer::hitEthrEE1_

private

Definition at line 133 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEE2_

const double HcalIsoTrkSimAnalyzer::hitEthrEE2_

private

Definition at line 134 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEE3_

const double HcalIsoTrkSimAnalyzer::hitEthrEE3_

private

Definition at line 134 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEEHi_

const double HcalIsoTrkSimAnalyzer::hitEthrEEHi_

private

Definition at line 135 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEELo_

const double HcalIsoTrkSimAnalyzer::hitEthrEELo_

private

Definition at line 135 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hltConfig_

HLTConfigProvider HcalIsoTrkSimAnalyzer::hltConfig_

private

Definition at line 119 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginRun().

◆ ignoreTrigger_

const bool HcalIsoTrkSimAnalyzer::ignoreTrigger_

private

Definition at line 131 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ l1Filter_

const std::string HcalIsoTrkSimAnalyzer::l1Filter_

private

Definition at line 121 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ l1GtUtils_

l1t::L1TGlobalUtil* HcalIsoTrkSimAnalyzer::l1GtUtils_

private

Definition at line 118 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ l1TrigName_

const std::string HcalIsoTrkSimAnalyzer::l1TrigName_

private

Definition at line 138 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ l2Filter_

const std::string HcalIsoTrkSimAnalyzer::l2Filter_

private

Definition at line 122 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ l3Filter_

const std::string HcalIsoTrkSimAnalyzer::l3Filter_

private

Definition at line 122 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ labelBS_

const std::string HcalIsoTrkSimAnalyzer::labelBS_

private

Definition at line 142 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ labelEB_

const std::string HcalIsoTrkSimAnalyzer::labelEB_

private

Definition at line 137 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ labelEE_

const std::string HcalIsoTrkSimAnalyzer::labelEE_

private

Definition at line 138 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ labelGenTrack_

const std::string HcalIsoTrkSimAnalyzer::labelGenTrack_

private

Definition at line 137 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ labelHBHE_

const std::string HcalIsoTrkSimAnalyzer::labelHBHE_

private

Definition at line 138 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ labelRecVtx_

const std::string HcalIsoTrkSimAnalyzer::labelRecVtx_

private

Definition at line 137 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ labelTower_

const std::string HcalIsoTrkSimAnalyzer::labelTower_

private

Definition at line 138 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ maxRestrictionP_

const double HcalIsoTrkSimAnalyzer::maxRestrictionP_

private

Definition at line 126 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ mode_

const int HcalIsoTrkSimAnalyzer::mode_

private

Definition at line 130 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ modnam_

const std::string HcalIsoTrkSimAnalyzer::modnam_

private

Definition at line 142 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ newDepth_

const std::vector<int> HcalIsoTrkSimAnalyzer::newDepth_

private

Definition at line 139 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and newId().

◆ nHigh_

unsigned int HcalIsoTrkSimAnalyzer::nHigh_

private

Definition at line 168 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree().

◆ nLow_

unsigned int HcalIsoTrkSimAnalyzer::nLow_

private

Definition at line 168 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree().

◆ nRun_

unsigned int HcalIsoTrkSimAnalyzer::nRun_

private

Definition at line 168 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginRun(), and endRun().

◆ oldDepth_

std::vector<int> HcalIsoTrkSimAnalyzer::oldDepth_

private

Definition at line 174 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and newId().

◆ oldDet_

std::vector<int> HcalIsoTrkSimAnalyzer::oldDet_

private

Definition at line 174 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and newId().

◆ oldEta_

std::vector<int> HcalIsoTrkSimAnalyzer::oldEta_

private

Definition at line 174 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and newId().

◆ oldID_

const std::vector<int> HcalIsoTrkSimAnalyzer::oldID_

private

Definition at line 139 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), HcalIsoTrkSimAnalyzer(), and newId().

◆ phibins_

std::vector<double> HcalIsoTrkSimAnalyzer::phibins_

private

Definition at line 173 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ phidist_

double HcalIsoTrkSimAnalyzer::phidist_

private

Definition at line 175 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ phihalfdist_

double HcalIsoTrkSimAnalyzer::phihalfdist_

private

Definition at line 175 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ prdnam_

const std::string HcalIsoTrkSimAnalyzer::prdnam_

private

Definition at line 142 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ prescaleHigh_

const int HcalIsoTrkSimAnalyzer::prescaleHigh_

private

Definition at line 129 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ prescaleLow_

const int HcalIsoTrkSimAnalyzer::prescaleLow_

private

Definition at line 129 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ processName_

const std::string HcalIsoTrkSimAnalyzer::processName_

private

Definition at line 121 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkSimAnalyzer().

◆ ptMin_

const double HcalIsoTrkSimAnalyzer::ptMin_

private

Definition at line 123 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ pTrackHigh_

const double HcalIsoTrkSimAnalyzer::pTrackHigh_

private

Definition at line 128 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ pTrackLow_

const double HcalIsoTrkSimAnalyzer::pTrackLow_

private

Definition at line 128 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ pTrackMin_

const double HcalIsoTrkSimAnalyzer::pTrackMin_

private

Definition at line 125 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ slopeRestrictionP_

const double HcalIsoTrkSimAnalyzer::slopeRestrictionP_

private

Definition at line 126 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ t_allvertex

int HcalIsoTrkSimAnalyzer::t_allvertex

private

Definition at line 196 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_DataType

int HcalIsoTrkSimAnalyzer::t_DataType

private

Definition at line 179 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_DetIds

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

private

Definition at line 191 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_DetIds1

std::vector<unsigned int> * HcalIsoTrkSimAnalyzer::t_DetIds1

private

Definition at line 191 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_DetIds3

std::vector<unsigned int> * HcalIsoTrkSimAnalyzer::t_DetIds3

private

Definition at line 191 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eAnnular

double HcalIsoTrkSimAnalyzer::t_eAnnular

private

Definition at line 187 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal

double HcalIsoTrkSimAnalyzer::t_eHcal

private

Definition at line 188 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal10

double HcalIsoTrkSimAnalyzer::t_eHcal10

private

Definition at line 188 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal30

double HcalIsoTrkSimAnalyzer::t_eHcal30

private

Definition at line 188 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_emaxNearP

double HcalIsoTrkSimAnalyzer::t_emaxNearP

private

Definition at line 187 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR

double HcalIsoTrkSimAnalyzer::t_eMipDR

private

Definition at line 185 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR2

double HcalIsoTrkSimAnalyzer::t_eMipDR2

private

Definition at line 185 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR3

double HcalIsoTrkSimAnalyzer::t_eMipDR3

private

Definition at line 185 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR4

double HcalIsoTrkSimAnalyzer::t_eMipDR4

private

Definition at line 185 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR5

double HcalIsoTrkSimAnalyzer::t_eMipDR5

private

Definition at line 186 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Event

int HcalIsoTrkSimAnalyzer::t_Event

private

Definition at line 179 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_EventNo

unsigned int HcalIsoTrkSimAnalyzer::t_EventNo

private

Definition at line 178 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_EventWeight

double HcalIsoTrkSimAnalyzer::t_EventWeight

private

Definition at line 181 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_gentrackP

double HcalIsoTrkSimAnalyzer::t_gentrackP

private

Definition at line 186 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_goodPV

int HcalIsoTrkSimAnalyzer::t_goodPV

private

Definition at line 180 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_hAnnular

double HcalIsoTrkSimAnalyzer::t_hAnnular

private

Definition at line 187 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies

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

private

Definition at line 192 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies1

std::vector<double> * HcalIsoTrkSimAnalyzer::t_HitEnergies1

private

Definition at line 192 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies3

std::vector<double> * HcalIsoTrkSimAnalyzer::t_HitEnergies3

private

Definition at line 192 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_hltbits

std::vector<bool> * HcalIsoTrkSimAnalyzer::t_hltbits

private

Definition at line 193 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_hmaxNearP

double HcalIsoTrkSimAnalyzer::t_hmaxNearP

private

Definition at line 186 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_ieta

int HcalIsoTrkSimAnalyzer::t_ieta

private

Definition at line 179 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_ietaAll

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

private

Definition at line 197 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ t_ietaGood

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

private

Definition at line 197 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ t_iphi

int HcalIsoTrkSimAnalyzer::t_iphi

private

Definition at line 179 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_L1Bit

bool HcalIsoTrkSimAnalyzer::t_L1Bit

private

Definition at line 194 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1eta

double HcalIsoTrkSimAnalyzer::t_l1eta

private

Definition at line 182 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1phi

double HcalIsoTrkSimAnalyzer::t_l1phi

private

Definition at line 182 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1pt

double HcalIsoTrkSimAnalyzer::t_l1pt

private

Definition at line 182 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3eta

double HcalIsoTrkSimAnalyzer::t_l3eta

private

Definition at line 183 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3phi

double HcalIsoTrkSimAnalyzer::t_l3phi

private

Definition at line 183 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3pt

double HcalIsoTrkSimAnalyzer::t_l3pt

private

Definition at line 183 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_mindR1

double HcalIsoTrkSimAnalyzer::t_mindR1

private

Definition at line 184 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_mindR2

double HcalIsoTrkSimAnalyzer::t_mindR2

private

Definition at line 184 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_nTrk

int HcalIsoTrkSimAnalyzer::t_nTrk

private

Definition at line 180 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_nVtx

int HcalIsoTrkSimAnalyzer::t_nVtx

private

Definition at line 180 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_p

double HcalIsoTrkSimAnalyzer::t_p

private

Definition at line 181 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_phi

double HcalIsoTrkSimAnalyzer::t_phi

private

Definition at line 181 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_pt

double HcalIsoTrkSimAnalyzer::t_pt

private

Definition at line 181 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyFlag

bool HcalIsoTrkSimAnalyzer::t_qltyFlag

private

Definition at line 189 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyMissFlag

bool HcalIsoTrkSimAnalyzer::t_qltyMissFlag

private

Definition at line 189 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyPVFlag

bool HcalIsoTrkSimAnalyzer::t_qltyPVFlag

private

Definition at line 190 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_rhoh

double HcalIsoTrkSimAnalyzer::t_rhoh

private

Definition at line 188 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Run

int HcalIsoTrkSimAnalyzer::t_Run

private

Definition at line 179 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_RunNo

unsigned int HcalIsoTrkSimAnalyzer::t_RunNo

private

Definition at line 178 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_selectTk

bool HcalIsoTrkSimAnalyzer::t_selectTk

private

Definition at line 189 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Tracks

int HcalIsoTrkSimAnalyzer::t_Tracks

private

Definition at line 195 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksLoose

int HcalIsoTrkSimAnalyzer::t_TracksLoose

private

Definition at line 196 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksProp

int HcalIsoTrkSimAnalyzer::t_TracksProp

private

Definition at line 195 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksSaved

int HcalIsoTrkSimAnalyzer::t_TracksSaved

private

Definition at line 195 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksTight

int HcalIsoTrkSimAnalyzer::t_TracksTight

private

Definition at line 196 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_trackType

std::vector<int> * HcalIsoTrkSimAnalyzer::t_trackType

private

Definition at line 197 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ t_trgbits

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

private

Definition at line 193 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ t_TrigPass

bool HcalIsoTrkSimAnalyzer::t_TrigPass

private

Definition at line 190 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TrigPassSel

bool HcalIsoTrkSimAnalyzer::t_TrigPassSel

private

Definition at line 190 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ theTriggerResultsLabel_

const edm::InputTag HcalIsoTrkSimAnalyzer::theTriggerResultsLabel_

private

Definition at line 136 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ tok_alg_

edm::EDGetTokenT<BXVector<GlobalAlgBlk> > HcalIsoTrkSimAnalyzer::tok_alg_

private

Definition at line 155 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_bFieldH_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> HcalIsoTrkSimAnalyzer::tok_bFieldH_

private

Definition at line 158 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_bs_

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

private

Definition at line 147 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_cala_

edm::EDGetTokenT<CaloTowerCollection> HcalIsoTrkSimAnalyzer::tok_cala_

private

Definition at line 154 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_caloTopology_

edm::ESGetToken<CaloTopology, CaloTopologyRecord> HcalIsoTrkSimAnalyzer::tok_caloTopology_

private

Definition at line 162 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_ddrec_

edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> HcalIsoTrkSimAnalyzer::tok_ddrec_

private

Definition at line 157 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginRun().

◆ tok_EB_

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkSimAnalyzer::tok_EB_

private

Definition at line 149 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_ecalChStatus_

edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> HcalIsoTrkSimAnalyzer::tok_ecalChStatus_

private

Definition at line 159 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_ecalPFRecHitThresholds_

edm::ESGetToken<EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd> HcalIsoTrkSimAnalyzer::tok_ecalPFRecHitThresholds_

private

Definition at line 166 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_EE_

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkSimAnalyzer::tok_EE_

private

Definition at line 150 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_ew_

edm::EDGetTokenT<GenEventInfoProduct> HcalIsoTrkSimAnalyzer::tok_ew_

private

Definition at line 152 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_geom_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HcalIsoTrkSimAnalyzer::tok_geom_

private

Definition at line 161 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> HcalIsoTrkSimAnalyzer::tok_hbhe_

private

Definition at line 151 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_htopo_

edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> HcalIsoTrkSimAnalyzer::tok_htopo_

private

Definition at line 163 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_parts_

edm::EDGetTokenT<reco::GenParticleCollection> HcalIsoTrkSimAnalyzer::tok_parts_

private

Definition at line 153 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_pdt_

edm::ESGetToken<HepPDT::ParticleDataTable, PDTRecord> HcalIsoTrkSimAnalyzer::tok_pdt_

private

Definition at line 165 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_recVtx_

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

private

Definition at line 148 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_resp_

edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> HcalIsoTrkSimAnalyzer::tok_resp_

private

Definition at line 164 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_sevlv_

edm::ESGetToken<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd> HcalIsoTrkSimAnalyzer::tok_sevlv_

private

Definition at line 160 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_trigEvt_

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

private

Definition at line 145 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_trigRes_

edm::EDGetTokenT<edm::TriggerResults> HcalIsoTrkSimAnalyzer::tok_trigRes_

private

Definition at line 146 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tree

TTree* HcalIsoTrkSimAnalyzer::tree

private

Definition at line 177 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eventstfile.Events::__iter__(), svgfig.SVG::__str__(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::beginLoop(), tree.Tree::branch_(), MetTreeProducer.MetTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::fillTree(), SimpleTreeProducer.SimpleTreeProducer::process(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::process(), MetTreeProducer.MetTreeProducer::process(), python.cmstools.EventTree::SetAlias(), eventstfile.Events::size(), eventstfile.Events::to(), tree.Tree::var(), and tree.Tree::vector().

◆ tree2

TTree * HcalIsoTrkSimAnalyzer::tree2

private

Definition at line 177 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ triggerEvent_

const edm::InputTag HcalIsoTrkSimAnalyzer::triggerEvent_

private

Definition at line 136 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ trigNames_

const std::vector<std::string> HcalIsoTrkSimAnalyzer::trigNames_

private

Definition at line 120 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ unCorrect_

const bool HcalIsoTrkSimAnalyzer::unCorrect_

private

Definition at line 132 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and storeEnergy().

◆ useL1Trigger_

const bool HcalIsoTrkSimAnalyzer::useL1Trigger_

private

Definition at line 131 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ usePFThresh_

const bool HcalIsoTrkSimAnalyzer::usePFThresh_

private

Definition at line 141 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ useRaw_

const int HcalIsoTrkSimAnalyzer::useRaw_

private

Definition at line 130 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HcalIsoTrkSimAnalyzer()

◆ ~HcalIsoTrkSimAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginRun()

◆ dR()

◆ endRun()

◆ eThreshold()

◆ fillDescriptions()

◆ fillTree()

◆ newId()

◆ notaMuon()

◆ rhoh()

◆ storeEnergy()

Member Data Documentation

◆ a_charIsoR_

◆ a_coneR1_

◆ a_coneR2_

◆ a_coneR_

◆ a_mipR2_

◆ a_mipR3_

◆ a_mipR4_

◆ a_mipR5_

◆ a_mipR_

◆ algTag_

◆ collapseDepth_

◆ dataType_

◆ eEcalMax_

◆ eIsolate1_

◆ eIsolate2_

◆ etabins_

◆ etadist_

◆ etahalfdist_

◆ etaMax_

◆ eThresholds_

◆ extTag_

◆ fillTreeRange_

◆ hcalScale_

◆ hdc_

◆ hep17_

◆ hitEthrEB_

◆ hitEthrEE0_

◆ hitEthrEE1_

◆ hitEthrEE2_

◆ hitEthrEE3_

◆ hitEthrEEHi_

◆ hitEthrEELo_

◆ hltConfig_

◆ ignoreTrigger_

◆ l1Filter_

◆ l1GtUtils_

◆ l1TrigName_

◆ l2Filter_

◆ l3Filter_

◆ labelBS_

◆ labelEB_

◆ labelEE_

◆ labelGenTrack_

◆ labelHBHE_

◆ labelRecVtx_

◆ labelTower_

◆ maxRestrictionP_

◆ mode_

◆ modnam_

◆ newDepth_

◆ nHigh_

◆ nLow_

◆ nRun_

◆ oldDepth_

◆ oldDet_

◆ oldEta_

◆ oldID_