Inheritance diagram for HcalIsoTrkAnalyzer:

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

	~HcalIsoTrkAnalyzer () override

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::propagatedTrackDirection > &trkCaloDirections, std::vector< spr::propagatedTrackID > &trkCaloDets, const CaloGeometry geo, const CaloTopology topo, const EcalChannelStatus theEcalChStatus, const EcalSeverityLevelAlgo theEcalSevlv, edm::Handle< EcalRecHitCollection > &barrelRecHitsHandle, edm::Handle< EcalRecHitCollection > &endcapRecHitsHandle, edm::Handle< HBHERecHitCollection > &hbhe, edm::Handle< CaloTowerCollection > &towerHandle, edm::Handle< reco::GenParticleCollection > &genParticles, const HcalDbService *conditions, const edm::Handle< reco::MuonCollection > &muonh)

double	gainFactor (const HcalDbService *dbserv, const HcalDetId &id)

DetId	newId (const DetId &)

bool	notaMuon (const reco::Track *pTrack0, const edm::Handle< reco::MuonCollection > &muonh)

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

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

double	trackP (const reco::Track *, const edm::Handle< reco::GenParticleCollection > &)

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 std::vector< int >	debEvents_

bool	debug_

const double	eEcalMax_

const double	eIsolate1_

const double	eIsolate2_

std::vector< double >	etabins_

double	etadist_

double	etahalfdist_

const EcalPFRecHitThresholds *	eThresholds_

const edm::InputTag	extTag_

const bool	fillInRange_

const bool	getCharge_

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	labelMuon_

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	pTrackHigh_

const double	pTrackLow_

const double	pTrackMin_

HcalRespCorrs *	respCorrs_

spr::trackSelectionParameters	selectionParameter_

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 HcalTopology *	theHBHETopology_

const std::string	theTrackQuality_

const edm::InputTag	theTriggerResultsLabel_

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

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	tok_bFieldH_

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

const edm::EDGetTokenT< CaloTowerCollection >	tok_cala_

const edm::ESGetToken< CaloTopology, CaloTopologyRecord >	tok_caloTopology_

const edm::ESGetToken< HcalDbService, HcalDbRecord >	tok_dbservice_

const edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord >	tok_ddrec_

const edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

const edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd >	tok_ecalChStatus_

const edm::ESGetToken< EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd >	tok_ecalPFRecHitThresholds_

const edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

const edm::EDGetTokenT< GenEventInfoProduct >	tok_ew_

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

const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_

const edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	tok_htopo_

const edm::EDGetTokenT< reco::MuonCollection >	tok_Muon_

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

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

const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd >	tok_resp_

const edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	tok_sevlv_

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

const 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 89 of file HcalIsoTrkAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalIsoTrkAnalyzer()

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

explicit

Definition at line 223 of file HcalIsoTrkAnalyzer.cc.

     : trigNames_(iConfig.getParameter<std::vector<std::string>>("triggers")),
       theTrackQuality_(iConfig.getParameter<std::string>("trackQuality")),
       processName_(iConfig.getParameter<std::string>("processName")),
       l1Filter_(iConfig.getParameter<std::string>("l1Filter")),
       l2Filter_(iConfig.getParameter<std::string>("l2Filter")),
       l3Filter_(iConfig.getParameter<std::string>("l3Filter")),
       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)),
       getCharge_(iConfig.getUntrackedParameter<bool>("getCharge", false)),
       collapseDepth_(iConfig.getUntrackedParameter<bool>("collapseDepth", false)),
       fillInRange_(iConfig.getUntrackedParameter<bool>("fillInRange", 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")),
       debEvents_(iConfig.getParameter<std::vector<int>>("debugEvents")),
       usePFThresh_(iConfig.getParameter<bool>("usePFThreshold")),
       labelBS_(iConfig.getParameter<std::string>("labelBeamSpot")),
       modnam_(iConfig.getUntrackedParameter<std::string>("moduleName", "")),
       prdnam_(iConfig.getUntrackedParameter<std::string>("producerName", "")),
       labelMuon_(iConfig.getParameter<std::string>("labelMuon")),
       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_parts_(consumes<reco::GenParticleCollection>(edm::InputTag("genParticles"))),
       tok_genTrack_(consumes<reco::TrackCollection>(labelGenTrack_)),
       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_cala_(consumes<CaloTowerCollection>(labelTower_)),
       tok_ew_(consumes<GenEventInfoProduct>(edm::InputTag("generator"))),
       tok_alg_(consumes<BXVector<GlobalAlgBlk>>(algTag_)),
       tok_Muon_(consumes<reco::MuonCollection>(labelMuon_)),
       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, edm::Transition::BeginRun>()),
       tok_resp_(esConsumes<HcalRespCorrs, HcalRespCorrsRcd, edm::Transition::BeginRun>()),
       tok_dbservice_(esConsumes<HcalDbService, HcalDbRecord>()),
       tok_ecalPFRecHitThresholds_(esConsumes<EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd>()),
       nRun_(0),
       nLow_(0),
       nHigh_(0),
       theHBHETopology_(nullptr),
       respCorrs_(nullptr),
       hdc_(nullptr) {
   usesResource(TFileService::kSharedResource);
 
   //now do whatever initialization is needed
   const double isolationRadius(28.9), innerR(10.0), outerR(30.0);
   reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality_);
   selectionParameter_.minPt = iConfig.getParameter<double>("minTrackPt");
   ;
   selectionParameter_.minQuality = trackQuality_;
   selectionParameter_.maxDxyPV = iConfig.getParameter<double>("maxDxyPV");
   selectionParameter_.maxDzPV = iConfig.getParameter<double>("maxDzPV");
   selectionParameter_.maxChi2 = iConfig.getParameter<double>("maxChi2");
   selectionParameter_.maxDpOverP = iConfig.getParameter<double>("maxDpOverP");
   selectionParameter_.minOuterHit = iConfig.getParameter<int>("minOuterHit");
   selectionParameter_.minLayerCrossed = iConfig.getParameter<int>("minLayerCrossed");
   selectionParameter_.maxInMiss = iConfig.getParameter<int>("maxInMiss");
   selectionParameter_.maxOutMiss = iConfig.getParameter<int>("maxOutMiss");
   a_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);
   // define tokens for access
 
   if (modnam_.empty()) {
     edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
                                      << labelBS_ << " " << labelRecVtx_ << " " << labelGenTrack_ << " "
                                      << edm::InputTag("ecalRecHit", labelEB_) << " "
                                      << edm::InputTag("ecalRecHit", labelEE_) << " " << labelHBHE_ << " " << labelTower_
                                      << " " << labelMuon_;
   } 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_ << " "
                                      << labelMuon_;
   }
 
   edm::LogVerbatim("HcalIsoTrack")
       << "Parameters read from config file \n"
       << "\t minPt " << selectionParameter_.minPt << "\t theTrackQuality " << theTrackQuality_ << "\t minQuality "
       << selectionParameter_.minQuality << "\t maxDxyPV " << selectionParameter_.maxDxyPV << "\t maxDzPV "
       << selectionParameter_.maxDzPV << "\t maxChi2 " << selectionParameter_.maxChi2 << "\t maxDpOverP "
       << selectionParameter_.maxDpOverP << "\t minOuterHit " << selectionParameter_.minOuterHit << "\t minLayerCrossed "
       << selectionParameter_.minLayerCrossed << "\t maxInMiss " << selectionParameter_.maxInMiss << "\t maxOutMiss "
       << selectionParameter_.maxOutMiss << "\t a_coneR " << a_coneR_ << ":" << 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 GetCharge " << getCharge_ << "\t FillInRange " << fillInRange_ << "\t L1TrigName_    "
       << l1TrigName_ << "\nThreshold flag used " << usePFThresh_ << " value for EB " << hitEthrEB_ << " EE "
       << hitEthrEE0_ << ":" << hitEthrEE1_ << ":" << hitEthrEE2_ << ":" << hitEthrEE3_ << ":" << hitEthrEELo_ << ":"
       << hitEthrEEHi_ << " and " << debEvents_.size() << " events to be debugged";
   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;
   }
 }

◆ ~HcalIsoTrkAnalyzer()

HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer ( )

override

Definition at line 417 of file HcalIsoTrkAnalyzer.cc.

References respCorrs_.

                                         {
   if (respCorrs_)
     delete respCorrs_;
 }

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 422 of file HcalIsoTrkAnalyzer.cc.

                                                                                     {
   t_Run = iEvent.id().run();
   t_Event = iEvent.id().event();
   t_DataType = dataType_;
   debug_ = (debEvents_.empty())
                ? true
                : (std::find(debEvents_.begin(), debEvents_.end(), iEvent.id().event()) != debEvents_.end());
 #ifdef EDM_ML_DEBUG
   if (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_);
   eThresholds_ = &iSetup.getData(tok_ecalPFRecHitThresholds_);
 
   // get calogeometry and calotopology
   const CaloGeometry* geo = &iSetup.getData(tok_geom_);
   const CaloTopology* caloTopology = &iSetup.getData(tok_caloTopology_);
 
   // get Hcal response corrections
   const HcalDbService* conditions = &iSetup.getData(tok_dbservice_);
 
   //=== genParticle information
   edm::Handle<reco::GenParticleCollection> genParticles = iEvent.getHandle(tok_parts_);
 
   bool okC(true);
   //Get track collection
   edm::Handle<reco::TrackCollection> trkCollection = iEvent.getHandle(tok_genTrack_);
   if (!trkCollection.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelGenTrack_;
     okC = false;
   }
 
   //Get muon collection
   const edm::Handle<reco::MuonCollection> muonh = iEvent.getHandle(tok_Muon_);
 
   //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
   if (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;
   }
   edm::Handle<CaloTowerCollection> caloTower = iEvent.getHandle(tok_cala_);
 
   //Propagate tracks to calorimeter surface)
   std::vector<spr::propagatedTrackDirection> trkCaloDirections;
   spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDirections, false);
   std::vector<spr::propagatedTrackID> trkCaloDets;
   spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDets, false);
   std::vector<math::XYZTLorentzVector> vecL1, vecL3;
   t_RunNo = iEvent.id().run();
   t_EventNo = iEvent.id().event();
   t_Tracks = trkCollection->size();
   t_TracksProp = trkCaloDirections.size();
   t_ietaAll->clear();
   t_ietaGood->clear();
   t_trackType->clear();
   t_trgbits->clear();
   t_hltbits->clear();
 #ifdef EDM_ML_DEBUG
   if (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 auto& 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
     if (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
               if (debug_)
                 edm::LogVerbatim("HcalIsoTrack")
                     << "This trigger " << names[iHLT] << " Flag " << hlt << ":" << t_trgbits->at(i);
 #endif
             }
           }
         }
       }
     }
 #ifdef EDM_ML_DEBUG
     if (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,
                          trkCaloDirections,
                          trkCaloDets,
                          geo,
                          caloTopology,
                          theEcalChStatus,
                          theEcalSevlv,
                          barrelRecHitsHandle,
                          endcapRecHitsHandle,
                          hbhe,
                          caloTower,
                          genParticles,
                          conditions,
                          muonh);
     t_TracksSaved = ntksave[0];
     t_TracksLoose = ntksave[1];
     t_TracksTight = ntksave[2];
   } else {
     trigger::TriggerEvent triggerEvent;
     const 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) {
       const edm::Handle<edm::TriggerResults> triggerResults = iEvent.getHandle(tok_trigRes_);
       triggerEvent = *(triggerEventHandle.product());
       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
                   if (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
                     if (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
                   if (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
               if (debug_)
                 edm::LogVerbatim("HcalIsoTrack") << "lvl2[" << i << "] dR " << dr;
 #endif
               if (dr < mindR1) {
                 mindR1 = dr;
                 mindRvec1 = vecL2[i];
               }
             }
 #ifdef EDM_ML_DEBUG
             if (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,
                                  trkCaloDirections,
                                  trkCaloDets,
                                  geo,
                                  caloTopology,
                                  theEcalChStatus,
                                  theEcalSevlv,
                                  barrelRecHitsHandle,
                                  endcapRecHitsHandle,
                                  hbhe,
                                  caloTower,
                                  genParticles,
                                  conditions,
                                  muonh);
               t_TracksSaved += ntksave[0];
               t_TracksLoose += ntksave[1];
               t_TracksTight += ntksave[2];
               done = true;
             }
           }
         }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   if (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 HcalIsoTrkAnalyzer::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 740 of file HcalIsoTrkAnalyzer.cc.

References make_classfiles::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 HcalIsoTrkAnalyzer::beginRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

overrideprivate

Definition at line 822 of file HcalIsoTrkAnalyzer.cc.

References RemoveAddSevLevel::flag, edm::EventSetup::getData(), hdc_, hltConfig_, ignoreTrigger_, HLTConfigProvider::init(), create_idmaps::n, nRun_, processName_, respCorrs_, edm::RunBase::run(), HcalCondObjectContainerBase::setTopo(), HLTConfigProvider::size(), theHBHETopology_, tok_ddrec_, tok_htopo_, tok_resp_, HLTConfigProvider::triggerIndex(), and trigNames_.

                                                                                  {
   theHBHETopology_ = &iSetup.getData(tok_htopo_);
   const HcalRespCorrs* resp = &iSetup.getData(tok_resp_);
   respCorrs_ = new HcalRespCorrs(*resp);
   respCorrs_->setTopo(theHBHETopology_);
   edm::LogVerbatim("HcalIsoTrack") << "beginRun " << iRun.run() << " get responseCoorection " << respCorrs_;
 
   hdc_ = &iSetup.getData(tok_ddrec_);
 
   if (!ignoreTrigger_) {
     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 HcalIsoTrkAnalyzer::dR	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 1349 of file HcalIsoTrkAnalyzer.cc.

References reco::deltaR().

Referenced by analyze(), fillTree(), notaMuon(), and trackP().

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

◆ endRun()

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

overrideprivate

Definition at line 858 of file HcalIsoTrkAnalyzer.cc.

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

                                                                         {
   nRun_++;
   edm::LogVerbatim("HcalIsoTrack") << "endRun[" << nRun_ << "] " << iRun.run();
 
   if (respCorrs_) {
     delete respCorrs_;
     respCorrs_ = nullptr;
   }
 }

◆ eThreshold()

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

private

Definition at line 1401 of file HcalIsoTrkAnalyzer.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 HcalIsoTrkAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 868 of file HcalIsoTrkAnalyzer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, events, 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");
   // following 10 parameters are parameters to select good tracks
   desc.add<std::string>("trackQuality", "highPurity");
   desc.add<double>("minTrackPt", 1.0);
   desc.add<double>("maxDxyPV", 0.02);
   desc.add<double>("maxDzPV", 0.02);
   desc.add<double>("maxChi2", 5.0);
   desc.add<double>("maxDpOverP", 0.1);
   desc.add<int>("minOuterHit", 4);
   desc.add<int>("minLayerCrossed", 8);
   desc.add<int>("maxInMiss", 0);
   desc.add<int>("maxOutMiss", 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<std::string>("labelMuon", "muons");
   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>("getCharge", false);
   desc.addUntracked<bool>("collapseDepth", false);
   desc.addUntracked<bool>("fillInRange", 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);
   std::vector<int> events;
   desc.add<std::vector<int>>("debugEvents", events);
   desc.add<bool>("usePFThreshold", true);
   descriptions.add("hcalIsoTrkAnalyzer", desc);
 }

◆ fillTree()

std::array< int, 3 > HcalIsoTrkAnalyzer::fillTree	(	std::vector< math::XYZTLorentzVector > &	vecL1,
		std::vector< math::XYZTLorentzVector > &	vecL3,
		math::XYZPoint &	leadPV,
		std::vector< spr::propagatedTrackDirection > &	trkCaloDirections,
		std::vector< spr::propagatedTrackID > &	trkCaloDets,
		const CaloGeometry *	geo,
		const CaloTopology *	topo,
		const EcalChannelStatus *	theEcalChStatus,
		const EcalSeverityLevelAlgo *	theEcalSevlv,
		edm::Handle< EcalRecHitCollection > &	barrelRecHitsHandle,
		edm::Handle< EcalRecHitCollection > &	endcapRecHitsHandle,
		edm::Handle< HBHERecHitCollection > &	hbhe,
		edm::Handle< CaloTowerCollection > &	towerHandle,
		edm::Handle< reco::GenParticleCollection > &	genParticles,
		const HcalDbService *	conditions,
		const edm::Handle< reco::MuonCollection > &	muonh
	)

private

Definition at line 962 of file HcalIsoTrkAnalyzer.cc.

References a_charIsoR_, a_coneR1_, a_coneR2_, a_coneR_, a_mipR2_, a_mipR3_, a_mipR4_, a_mipR5_, a_mipR_, funct::abs(), accept(), spr::chargeIsolationCone(), spr::chargeIsolationEcal(), PixelBaryCentreAnalyzer_cfg::conditions, debug_, hcalRecHitTable_cff::detId, l1ctLayer1_cff::dr, dR(), spr::eCone_ecal(), spr::eCone_hcal(), spr::eECALmatrix(), eEcalMax_, spr::eHCALmatrix(), eIsolate1_, eIsolate2_, reco::TrackBase::eta(), eThreshold(), JetChargeProducer_cfi::exp, fillInRange_, AJJGenJetFilter_cfi::genParticles, spr::goodTrack(), photonIsolationHIProducer_cfi::hbhe, testProducerWithPsetDescEmpty_cfi::i3, edm::HandleBase::isValid(), isotrackApplyRegressor::k, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, maxRestrictionP_, newId(), nHigh_, jetSelector_cfi::nLoose, nLow_, notaMuon(), beamSpotPI::nTracks, oldID_, reco::TrackBase::p(), reco::TrackBase::phi(), prescaleHigh_, prescaleLow_, reco::TrackBase::pt(), pTrackHigh_, pTrackLow_, pTrackMin_, reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), rhoh(), selectionParameter_, 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_EventNo, 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_RunNo, t_selectTk, t_trackType, t_trgbits, theHBHETopology_, l1tHGCalTowerProducer_cfi::tower, trackP(), 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 = trkCaloDirections.size();
   t_rhoh = (tower.isValid()) ? rhoh(tower) : 0;
   for (const auto& trkDetItr : trkCaloDirections) {
     const reco::Track* pTrack = &(*(trkDetItr.trkItr));
     math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
     t_p = pTrack->p();
     t_pt = pTrack->pt();
     t_phi = pTrack->phi();
 
 #ifdef EDM_ML_DEBUG
     if (debug_)
       edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) : " << t_pt << "|"
                                        << pTrack->eta() << "|" << t_phi << "|" << t_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
     if (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 = spr::goodTrack(pTrack, leadPV, selectionParameter_, false);
     spr::trackSelectionParameters oneCutParameters = selectionParameter_;
     oneCutParameters.maxDxyPV = 10;
     oneCutParameters.maxDzPV = 100;
     oneCutParameters.maxInMiss = 2;
     oneCutParameters.maxOutMiss = 2;
     bool qltyFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
     oneCutParameters = selectionParameter_;
     oneCutParameters.maxDxyPV = 10;
     oneCutParameters.maxDzPV = 100;
     t_qltyMissFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
     oneCutParameters = selectionParameter_;
     oneCutParameters.maxInMiss = 2;
     oneCutParameters.maxOutMiss = 2;
     t_qltyPVFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
     double eIsolation = maxRestrictionP_ * exp(slopeRestrictionP_ * std::abs((double)t_ieta));
     if (eIsolation < eIsolate1_)
       eIsolation = eIsolate1_;
     if (eIsolation < eIsolate2_)
       eIsolation = eIsolate2_;
 #ifdef EDM_ML_DEBUG
     if (debug_)
       edm::LogVerbatim("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " << qltyFlag << "|" << trkDetItr.okECAL << "|"
                                        << trkDetItr.okHCAL << " eIsolation " << eIsolation;
 #endif
     t_qltyFlag = (qltyFlag && trkDetItr.okECAL && trkDetItr.okHCAL);
     bool notMuon = (muonh.isValid()) ? notaMuon(pTrack, muonh) : true;
     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
       if (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
       if (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
       if (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
       if (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
       if (debug_)
         edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR5 << ":" << eEcal5;
 #endif
       t_eMipDR5 = eEcal5;
 
       t_emaxNearP = spr::chargeIsolationEcal(nTracks, trkCaloDets, 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::chargeIsolationCone(nTracks, trkCaloDirections, 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
       if (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 = trackP(pTrack, genParticles);
       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, conditions, 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, conditions, 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, conditions, ids3, edet3, t_eHcal30, t_DetIds3, t_HitEnergies3);
 
 #ifdef EDM_ML_DEBUG
         if (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 = " << HcalDetId(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 = " << HcalDetId(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 = " << HcalDetId(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 (fillInRange_) {
           if ((t_p < pTrackLow_) || (t_p > pTrackHigh_))
             accept = false;
         }
         if (accept) {
           tree->Fill();
           edm::LogVerbatim("HcalIsoTrackX")
               << "Run " << t_RunNo << " Event " << t_EventNo << " Track " << nTracks << " p " << t_p;
           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
           if (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;
 }

◆ gainFactor()

double HcalIsoTrkAnalyzer::gainFactor	(	const HcalDbService *	dbserv,
		const HcalDetId &	id
	)

private

Definition at line 1524 of file HcalIsoTrkAnalyzer.cc.

References PixelBaryCentreAnalyzer_cfg::conditions, Hcal_Conditions_forGlobalTag_cff::gain, and HcalCalibrations::respcorrgain().

Referenced by storeEnergy().

                                                                                           {
   double gain(0.0);
   const HcalCalibrations& calibs = conditions->getHcalCalibrations(id);
   for (int capid = 0; capid < 4; ++capid)
     gain += (0.25 * calibs.respcorrgain(capid));
   return gain;
 }

◆ newId()

DetId HcalIsoTrkAnalyzer::newId ( const DetId & id )

private

Definition at line 1419 of file HcalIsoTrkAnalyzer.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 HcalIsoTrkAnalyzer::notaMuon	(	const reco::Track *	pTrack0,
		const edm::Handle< reco::MuonCollection > &	muonh
	)

private

Definition at line 1499 of file HcalIsoTrkAnalyzer.cc.

References reco::deltaR(), dR(), reco::TrackBase::eta(), RemoveAddSevLevel::flag, reco::TrackBase::phi(), and muon::segmentCompatibility().

Referenced by fillTree().

                                                                                                       {
   bool flag(true);
   for (reco::MuonCollection::const_iterator recMuon = muonh->begin(); recMuon != muonh->end(); ++recMuon) {
     if (recMuon->innerTrack().isNonnull()) {
       const reco::Track* pTrack = (recMuon->innerTrack()).get();
       bool mediumMuon = (((recMuon->isPFMuon()) && (recMuon->isGlobalMuon() || recMuon->isTrackerMuon())) &&
                          (recMuon->innerTrack()->validFraction() > 0.49));
       if (mediumMuon) {
         double chiGlobal = ((recMuon->globalTrack().isNonnull()) ? recMuon->globalTrack()->normalizedChi2() : 999);
         bool goodGlob = (recMuon->isGlobalMuon() && chiGlobal < 3 &&
                          recMuon->combinedQuality().chi2LocalPosition < 12 && recMuon->combinedQuality().trkKink < 20);
         mediumMuon = muon::segmentCompatibility(*recMuon) > (goodGlob ? 0.303 : 0.451);
       }
       if (mediumMuon) {
         double dR = reco::deltaR(pTrack->eta(), pTrack->phi(), pTrack0->eta(), pTrack0->phi());
         if (dR < 0.1) {
           flag = false;
           break;
         }
       }
     }
   }
   return flag;
 }

◆ rhoh()

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

private

Definition at line 1369 of file HcalIsoTrkAnalyzer.cc.

References debug_, 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
   if (debug_)
     edm::LogVerbatim("HcalIsoTrack") << "Rho " << evt_smdq << ":" << rhoh;
 #endif
   return rhoh;
 }

◆ storeEnergy()

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

private

Definition at line 1431 of file HcalIsoTrkAnalyzer.cc.

References funct::abs(), collapseDepth_, PixelBaryCentreAnalyzer_cfg::conditions, alignCSCRings::corr, debug_, Hcal_Conditions_forGlobalTag_cff::gain, gainFactor(), getCharge_, reco::JetExtendedAssociation::getValue(), HcalCondObjectContainer< Item >::getValues(), hcalScale_, hdc_, EcalPhiSymFlatTableProducers_cfi::id, isotrackApplyRegressor::k, HcalDDDRecConstants::mergedDepthDetId(), nano_mu_digi_cff::rawId, respCorrs_, 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 if (getCharge_) {
     for (unsigned int k = 0; k < ids.size(); ++k) {
       double gain = gainFactor(conditions, HcalDetId(ids[k]));
       if (gain != 0)
         edet[k] /= gain;
       ehcal += edet[k];
     }
   } else {
     for (const auto& en : edet)
       ehcal += en;
   }
   if ((std::abs(ehcal - eHcal) > 0.001) && (!unCorrect_) && (!getCharge_))
     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
   if (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
 }

◆ trackP()

double HcalIsoTrkAnalyzer::trackP	(	const reco::Track *	pTrack,
		const edm::Handle< reco::GenParticleCollection > &	genParticles
	)

private

Definition at line 1353 of file HcalIsoTrkAnalyzer.cc.

References reco::deltaR(), dR(), reco::TrackBase::eta(), AJJGenJetFilter_cfi::genParticles, AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::phi(), and isotrackTrainRegressor::pmom.

Referenced by fillTree().

                                                                                             {
   double pmom = -1.0;
   if (genParticles.isValid()) {
     double mindR(999.9);
     for (const auto& p : (*genParticles)) {
       double dR = reco::deltaR(pTrack->eta(), pTrack->phi(), p.momentum().Eta(), p.momentum().Phi());
       if (dR < mindR) {
         mindR = dR;
         pmom = p.momentum().R();
       }
     }
   }
   return pmom;
 }

Member Data Documentation

◆ a_charIsoR_

double HcalIsoTrkAnalyzer::a_charIsoR_

private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_coneR1_

double HcalIsoTrkAnalyzer::a_coneR1_

private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_coneR2_

double HcalIsoTrkAnalyzer::a_coneR2_

private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_coneR_

const double HcalIsoTrkAnalyzer::a_coneR_

private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR2_

const double HcalIsoTrkAnalyzer::a_mipR2_

private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR3_

const double HcalIsoTrkAnalyzer::a_mipR3_

private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR4_

const double HcalIsoTrkAnalyzer::a_mipR4_

private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR5_

const double HcalIsoTrkAnalyzer::a_mipR5_

private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR_

const double HcalIsoTrkAnalyzer::a_mipR_

private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ algTag_

const edm::InputTag HcalIsoTrkAnalyzer::algTag_

private

Definition at line 160 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ collapseDepth_

const bool HcalIsoTrkAnalyzer::collapseDepth_

private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

◆ dataType_

const int HcalIsoTrkAnalyzer::dataType_

private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ debEvents_

const std::vector<int> HcalIsoTrkAnalyzer::debEvents_

private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ debug_

bool HcalIsoTrkAnalyzer::debug_

private

Definition at line 220 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), fillTree(), rhoh(), and storeEnergy().

◆ eEcalMax_

const double HcalIsoTrkAnalyzer::eEcalMax_

private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ eIsolate1_

const double HcalIsoTrkAnalyzer::eIsolate1_

private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ eIsolate2_

const double HcalIsoTrkAnalyzer::eIsolate2_

private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ etabins_

std::vector<double> HcalIsoTrkAnalyzer::etabins_

private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ etadist_

double HcalIsoTrkAnalyzer::etadist_

private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ etahalfdist_

double HcalIsoTrkAnalyzer::etahalfdist_

private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ eThresholds_

const EcalPFRecHitThresholds* HcalIsoTrkAnalyzer::eThresholds_

private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ extTag_

const edm::InputTag HcalIsoTrkAnalyzer::extTag_

private

Definition at line 160 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ fillInRange_

const bool HcalIsoTrkAnalyzer::fillInRange_

private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ getCharge_

const bool HcalIsoTrkAnalyzer::getCharge_

private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

◆ hcalScale_

const double HcalIsoTrkAnalyzer::hcalScale_

private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

◆ hdc_

const HcalDDDRecConstants* HcalIsoTrkAnalyzer::hdc_

private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and storeEnergy().

◆ hep17_

const bool HcalIsoTrkAnalyzer::hep17_

private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ hitEthrEB_

const double HcalIsoTrkAnalyzer::hitEthrEB_

private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE0_

const double HcalIsoTrkAnalyzer::hitEthrEE0_

private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE1_

const double HcalIsoTrkAnalyzer::hitEthrEE1_

private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE2_

const double HcalIsoTrkAnalyzer::hitEthrEE2_

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE3_

const double HcalIsoTrkAnalyzer::hitEthrEE3_

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEEHi_

const double HcalIsoTrkAnalyzer::hitEthrEEHi_

private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEELo_

const double HcalIsoTrkAnalyzer::hitEthrEELo_

private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hltConfig_

HLTConfigProvider HcalIsoTrkAnalyzer::hltConfig_

private

Definition at line 134 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginRun().

◆ ignoreTrigger_

const bool HcalIsoTrkAnalyzer::ignoreTrigger_

private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

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

◆ l1Filter_

const std::string HcalIsoTrkAnalyzer::l1Filter_

private

Definition at line 138 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l1GtUtils_

l1t::L1TGlobalUtil* HcalIsoTrkAnalyzer::l1GtUtils_

private

Definition at line 133 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l1TrigName_

const std::string HcalIsoTrkAnalyzer::l1TrigName_

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l2Filter_

const std::string HcalIsoTrkAnalyzer::l2Filter_

private

Definition at line 139 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l3Filter_

const std::string HcalIsoTrkAnalyzer::l3Filter_

private

Definition at line 139 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelBS_

const std::string HcalIsoTrkAnalyzer::labelBS_

private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ labelEB_

const std::string HcalIsoTrkAnalyzer::labelEB_

private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelEE_

const std::string HcalIsoTrkAnalyzer::labelEE_

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelGenTrack_

const std::string HcalIsoTrkAnalyzer::labelGenTrack_

private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelHBHE_

const std::string HcalIsoTrkAnalyzer::labelHBHE_

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelMuon_

const std::string HcalIsoTrkAnalyzer::labelMuon_

private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ labelRecVtx_

const std::string HcalIsoTrkAnalyzer::labelRecVtx_

private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ labelTower_

const std::string HcalIsoTrkAnalyzer::labelTower_

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ maxRestrictionP_

const double HcalIsoTrkAnalyzer::maxRestrictionP_

private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ mode_

const int HcalIsoTrkAnalyzer::mode_

private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ modnam_

const std::string HcalIsoTrkAnalyzer::modnam_

private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ newDepth_

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

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ nHigh_

unsigned int HcalIsoTrkAnalyzer::nHigh_

private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree().

◆ nLow_

unsigned int HcalIsoTrkAnalyzer::nLow_

private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree().

◆ nRun_

unsigned int HcalIsoTrkAnalyzer::nRun_

private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and endRun().

◆ oldDepth_

std::vector<int> HcalIsoTrkAnalyzer::oldDepth_

private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ oldDet_

std::vector<int> HcalIsoTrkAnalyzer::oldDet_

private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ oldEta_

std::vector<int> HcalIsoTrkAnalyzer::oldEta_

private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ oldID_

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

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

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

◆ phibins_

std::vector<double> HcalIsoTrkAnalyzer::phibins_

private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ phidist_

double HcalIsoTrkAnalyzer::phidist_

private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ phihalfdist_

double HcalIsoTrkAnalyzer::phihalfdist_

private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ prdnam_

const std::string HcalIsoTrkAnalyzer::prdnam_

private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ prescaleHigh_

const int HcalIsoTrkAnalyzer::prescaleHigh_

private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ prescaleLow_

const int HcalIsoTrkAnalyzer::prescaleLow_

private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ processName_

const std::string HcalIsoTrkAnalyzer::processName_

private

Definition at line 138 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkAnalyzer().

◆ pTrackHigh_

const double HcalIsoTrkAnalyzer::pTrackHigh_

private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ pTrackLow_

const double HcalIsoTrkAnalyzer::pTrackLow_

private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ pTrackMin_

const double HcalIsoTrkAnalyzer::pTrackMin_

private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ respCorrs_

HcalRespCorrs* HcalIsoTrkAnalyzer::respCorrs_

private

Definition at line 190 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), endRun(), storeEnergy(), and ~HcalIsoTrkAnalyzer().

◆ selectionParameter_

spr::trackSelectionParameters HcalIsoTrkAnalyzer::selectionParameter_

private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ slopeRestrictionP_

const double HcalIsoTrkAnalyzer::slopeRestrictionP_

private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ t_allvertex

int HcalIsoTrkAnalyzer::t_allvertex

private

Definition at line 217 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_DataType

int HcalIsoTrkAnalyzer::t_DataType

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_DetIds

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

private

Definition at line 212 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_DetIds1

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

private

Definition at line 212 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_DetIds3

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

private

Definition at line 212 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eAnnular

double HcalIsoTrkAnalyzer::t_eAnnular

private

Definition at line 208 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal

double HcalIsoTrkAnalyzer::t_eHcal

private

Definition at line 209 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal10

double HcalIsoTrkAnalyzer::t_eHcal10

private

Definition at line 209 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal30

double HcalIsoTrkAnalyzer::t_eHcal30

private

Definition at line 209 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_emaxNearP

double HcalIsoTrkAnalyzer::t_emaxNearP

private

Definition at line 208 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR

double HcalIsoTrkAnalyzer::t_eMipDR

private

Definition at line 206 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR2

double HcalIsoTrkAnalyzer::t_eMipDR2

private

Definition at line 206 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR3

double HcalIsoTrkAnalyzer::t_eMipDR3

private

Definition at line 206 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR4

double HcalIsoTrkAnalyzer::t_eMipDR4

private

Definition at line 206 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR5

double HcalIsoTrkAnalyzer::t_eMipDR5

private

Definition at line 207 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Event

int HcalIsoTrkAnalyzer::t_Event

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_EventNo

unsigned int HcalIsoTrkAnalyzer::t_EventNo

private

Definition at line 199 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_EventWeight

double HcalIsoTrkAnalyzer::t_EventWeight

private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_gentrackP

double HcalIsoTrkAnalyzer::t_gentrackP

private

Definition at line 207 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_goodPV

int HcalIsoTrkAnalyzer::t_goodPV

private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_hAnnular

double HcalIsoTrkAnalyzer::t_hAnnular

private

Definition at line 208 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies

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

private

Definition at line 213 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies1

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

private

Definition at line 213 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies3

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

private

Definition at line 213 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_hltbits

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

private

Definition at line 214 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_hmaxNearP

double HcalIsoTrkAnalyzer::t_hmaxNearP

private

Definition at line 207 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_ieta

int HcalIsoTrkAnalyzer::t_ieta

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_ietaAll

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

private

Definition at line 218 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_ietaGood

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

private

Definition at line 218 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_iphi

int HcalIsoTrkAnalyzer::t_iphi

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_L1Bit

bool HcalIsoTrkAnalyzer::t_L1Bit

private

Definition at line 215 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1eta

double HcalIsoTrkAnalyzer::t_l1eta

private

Definition at line 203 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1phi

double HcalIsoTrkAnalyzer::t_l1phi

private

Definition at line 203 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1pt

double HcalIsoTrkAnalyzer::t_l1pt

private

Definition at line 203 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3eta

double HcalIsoTrkAnalyzer::t_l3eta

private

Definition at line 204 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3phi

double HcalIsoTrkAnalyzer::t_l3phi

private

Definition at line 204 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3pt

double HcalIsoTrkAnalyzer::t_l3pt

private

Definition at line 204 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_mindR1

double HcalIsoTrkAnalyzer::t_mindR1

private

Definition at line 205 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_mindR2

double HcalIsoTrkAnalyzer::t_mindR2

private

Definition at line 205 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_nTrk

int HcalIsoTrkAnalyzer::t_nTrk

private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_nVtx

int HcalIsoTrkAnalyzer::t_nVtx

private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_p

double HcalIsoTrkAnalyzer::t_p

private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_phi

double HcalIsoTrkAnalyzer::t_phi

private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_pt

double HcalIsoTrkAnalyzer::t_pt

private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyFlag

bool HcalIsoTrkAnalyzer::t_qltyFlag

private

Definition at line 210 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyMissFlag

bool HcalIsoTrkAnalyzer::t_qltyMissFlag

private

Definition at line 210 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyPVFlag

bool HcalIsoTrkAnalyzer::t_qltyPVFlag

private

Definition at line 211 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_rhoh

double HcalIsoTrkAnalyzer::t_rhoh

private

Definition at line 209 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Run

int HcalIsoTrkAnalyzer::t_Run

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_RunNo

unsigned int HcalIsoTrkAnalyzer::t_RunNo

private

Definition at line 199 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_selectTk

bool HcalIsoTrkAnalyzer::t_selectTk

private

Definition at line 210 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Tracks

int HcalIsoTrkAnalyzer::t_Tracks

private

Definition at line 216 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksLoose

int HcalIsoTrkAnalyzer::t_TracksLoose

private

Definition at line 217 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksProp

int HcalIsoTrkAnalyzer::t_TracksProp

private

Definition at line 216 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksSaved

int HcalIsoTrkAnalyzer::t_TracksSaved

private

Definition at line 216 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksTight

int HcalIsoTrkAnalyzer::t_TracksTight

private

Definition at line 217 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_trackType

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

private

Definition at line 218 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_trgbits

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

private

Definition at line 214 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_TrigPass

bool HcalIsoTrkAnalyzer::t_TrigPass

private

Definition at line 211 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TrigPassSel

bool HcalIsoTrkAnalyzer::t_TrigPassSel

private

Definition at line 211 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ theHBHETopology_

const HcalTopology* HcalIsoTrkAnalyzer::theHBHETopology_

private

Definition at line 189 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and fillTree().

◆ theTrackQuality_

const std::string HcalIsoTrkAnalyzer::theTrackQuality_

private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ theTriggerResultsLabel_

const edm::InputTag HcalIsoTrkAnalyzer::theTriggerResultsLabel_

private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ tok_alg_

const edm::EDGetTokenT<BXVector<GlobalAlgBlk> > HcalIsoTrkAnalyzer::tok_alg_

private

Definition at line 173 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_bFieldH_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> HcalIsoTrkAnalyzer::tok_bFieldH_

private

Definition at line 177 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_bs_

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

private

Definition at line 166 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_cala_

const edm::EDGetTokenT<CaloTowerCollection> HcalIsoTrkAnalyzer::tok_cala_

private

Definition at line 171 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_caloTopology_

const edm::ESGetToken<CaloTopology, CaloTopologyRecord> HcalIsoTrkAnalyzer::tok_caloTopology_

private

Definition at line 181 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_dbservice_

const edm::ESGetToken<HcalDbService, HcalDbRecord> HcalIsoTrkAnalyzer::tok_dbservice_

private

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_ddrec_

const edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> HcalIsoTrkAnalyzer::tok_ddrec_

private

Definition at line 176 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun().

◆ tok_EB_

const edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_EB_

private

Definition at line 168 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_ecalChStatus_

const edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> HcalIsoTrkAnalyzer::tok_ecalChStatus_

private

Definition at line 178 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_ecalPFRecHitThresholds_

const edm::ESGetToken<EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd> HcalIsoTrkAnalyzer::tok_ecalPFRecHitThresholds_

private

Definition at line 185 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_EE_

const edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_EE_

private

Definition at line 169 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_ew_

const edm::EDGetTokenT<GenEventInfoProduct> HcalIsoTrkAnalyzer::tok_ew_

private

Definition at line 172 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_genTrack_

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

private

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_geom_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HcalIsoTrkAnalyzer::tok_geom_

private

Definition at line 180 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_hbhe_

const edm::EDGetTokenT<HBHERecHitCollection> HcalIsoTrkAnalyzer::tok_hbhe_

private

Definition at line 170 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_htopo_

const edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> HcalIsoTrkAnalyzer::tok_htopo_

private

Definition at line 182 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun().

◆ tok_Muon_

const edm::EDGetTokenT<reco::MuonCollection> HcalIsoTrkAnalyzer::tok_Muon_

private

Definition at line 174 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_parts_

const edm::EDGetTokenT<reco::GenParticleCollection> HcalIsoTrkAnalyzer::tok_parts_

private

Definition at line 164 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_recVtx_

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

private

Definition at line 167 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_resp_

const edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> HcalIsoTrkAnalyzer::tok_resp_

private

Definition at line 183 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun().

◆ tok_sevlv_

const edm::ESGetToken<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd> HcalIsoTrkAnalyzer::tok_sevlv_

private

Definition at line 179 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_trigEvt_

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

private

Definition at line 162 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tok_trigRes_

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

private

Definition at line 163 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze().

◆ tree

TTree* HcalIsoTrkAnalyzer::tree

private

Definition at line 198 of file HcalIsoTrkAnalyzer.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 * HcalIsoTrkAnalyzer::tree2

private

Definition at line 198 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ triggerEvent_

const edm::InputTag HcalIsoTrkAnalyzer::triggerEvent_

private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ trigNames_

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

private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

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

◆ unCorrect_

const bool HcalIsoTrkAnalyzer::unCorrect_

private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

◆ useL1Trigger_

const bool HcalIsoTrkAnalyzer::useL1Trigger_

private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ usePFThresh_

const bool HcalIsoTrkAnalyzer::usePFThresh_

private

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ useRaw_

const int HcalIsoTrkAnalyzer::useRaw_

private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HcalIsoTrkAnalyzer()

◆ ~HcalIsoTrkAnalyzer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginRun()

◆ dR()

◆ endRun()

◆ eThreshold()

◆ fillDescriptions()

◆ fillTree()

◆ gainFactor()

◆ newId()

◆ notaMuon()

◆ rhoh()

◆ storeEnergy()

◆ trackP()

Member Data Documentation

◆ a_charIsoR_

◆ a_coneR1_

◆ a_coneR2_

◆ a_coneR_

◆ a_mipR2_

◆ a_mipR3_

◆ a_mipR4_

◆ a_mipR5_

◆ a_mipR_

◆ algTag_

◆ collapseDepth_

◆ dataType_

◆ debEvents_

◆ debug_

◆ eEcalMax_

◆ eIsolate1_

◆ eIsolate2_

◆ etabins_

◆ etadist_

◆ etahalfdist_

◆ eThresholds_

◆ extTag_

◆ fillInRange_

◆ getCharge_

◆ hcalScale_

◆ hdc_

◆ hep17_

◆ hitEthrEB_

◆ hitEthrEE0_

◆ hitEthrEE1_

◆ hitEthrEE2_

◆ hitEthrEE3_

◆ hitEthrEEHi_

◆ hitEthrEELo_

◆ hltConfig_

◆ ignoreTrigger_

◆ l1Filter_

◆ l1GtUtils_

◆ l1TrigName_

◆ l2Filter_

◆ l3Filter_

◆ labelBS_

◆ labelEB_

◆ labelEE_

◆ labelGenTrack_

◆ labelHBHE_

◆ labelMuon_

◆ labelRecVtx_

◆ labelTower_

◆ maxRestrictionP_

◆ mode_

◆ modnam_

◆ newDepth_

◆ nHigh_

◆ nLow_