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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const 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

bool	wantsStreamRuns () const

	~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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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

std::vector< ESProxyIndex > 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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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 HcalTopology theHBHETopology, 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 HcalRespCorrs respCorrs)

DetId	newId (const DetId &)

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

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

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 bool	collapseDepth_

const int	dataType_

const double	eEcalMax_

const double	eIsolate1_

const double	eIsolate2_

std::vector< double >	etabins_

double	etadist_

double	etahalfdist_

edm::Service< TFileService >	fs

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	labelEB_

const std::string	labelEE_

const std::string	labelGenTrack_

const std::string	labelHBHE_

const std::string	labelRecVtx_

const std::string	labelTower_

const double	maxRestrictionP_

const int	mode_

const std::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 int	prescaleHigh_

const int	prescaleLow_

const std::string	processName_

const double	pTrackHigh_

const double	pTrackLow_

const double	pTrackMin_

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 std::string	theTrackQuality_

const edm::InputTag	theTriggerResultsLabel_

edm::EDGetTokenT< BXVector< GlobalAlgBlk > >	tok_alg_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	tok_bFieldH_

edm::EDGetTokenT< reco::BeamSpot >	tok_bs_

edm::EDGetTokenT< CaloTowerCollection >	tok_cala_

edm::ESGetToken< CaloTopology, CaloTopologyRecord >	tok_caloTopology_

edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord >	tok_ddrec_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd >	tok_ecalChStatus_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT< GenEventInfoProduct >	tok_ew_

edm::EDGetTokenT< reco::TrackCollection >	tok_genTrack_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	tok_htopo_

edm::EDGetTokenT< reco::GenParticleCollection >	tok_parts_

edm::EDGetTokenT< reco::VertexCollection >	tok_recVtx_

edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd >	tok_resp_

edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	tok_sevlv_

edm::EDGetTokenT< trigger::TriggerEvent >	tok_trigEvt_

edm::EDGetTokenT< edm::TriggerResults >	tok_trigRes_

TTree *	tree

TTree *	tree2

const edm::InputTag	triggerEvent_

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

const bool	unCorrect_

const bool	useL1Trigger_

const 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
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

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

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

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

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

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

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

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

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

Detailed Description

Definition at line 87 of file HcalIsoTrkAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalIsoTrkAnalyzer()

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

explicit

Definition at line 205 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)),
       collapseDepth_(iConfig.getUntrackedParameter<bool>("collapseDepth", 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")),
       nRun_(0),
       nLow_(0),
       nHigh_(0),
       hdc_(nullptr) {
   usesResource(TFileService::kSharedResource);
  
   //now do whatever initialization is needed
   const double isolationRadius(28.9), innerR(10.0), outerR(30.0);
   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
   std::string labelBS = iConfig.getParameter<std::string>("labelBeamSpot");
   std::string modnam = iConfig.getUntrackedParameter<std::string>("moduleName", "");
   std::string prdnam = iConfig.getUntrackedParameter<std::string>("producerName", "");
   edm::InputTag algTag = iConfig.getParameter<edm::InputTag>("algInputTag");
   edm::InputTag extTag = iConfig.getParameter<edm::InputTag>("extInputTag");
   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
   tok_trigEvt_ = consumes<trigger::TriggerEvent>(triggerEvent_);
   tok_trigRes_ = consumes<edm::TriggerResults>(theTriggerResultsLabel_);
   tok_bs_ = consumes<reco::BeamSpot>(labelBS);
   tok_genTrack_ = consumes<reco::TrackCollection>(labelGenTrack_);
   tok_ew_ = consumes<GenEventInfoProduct>(edm::InputTag("generator"));
   tok_parts_ = consumes<reco::GenParticleCollection>(edm::InputTag("genParticles"));
   tok_cala_ = consumes<CaloTowerCollection>(labelTower_);
   tok_alg_ = consumes<BXVector<GlobalAlgBlk>>(algTag);
  
   if (modnam.empty()) {
     tok_recVtx_ = consumes<reco::VertexCollection>(labelRecVtx_);
     tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEB_));
     tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEE_));
     tok_hbhe_ = consumes<HBHERecHitCollection>(labelHBHE_);
     edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
                                      << labelBS << " " << labelRecVtx_ << " " << labelGenTrack_ << " "
                                      << edm::InputTag("ecalRecHit", labelEB_) << " "
                                      << edm::InputTag("ecalRecHit", labelEE_) << " " << labelHBHE_ << " "
                                      << labelTower_;
   } else {
     tok_recVtx_ = consumes<reco::VertexCollection>(edm::InputTag(modnam, labelRecVtx_, prdnam));
     tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag(modnam, labelEB_, prdnam));
     tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag(modnam, labelEE_, prdnam));
     tok_hbhe_ = consumes<HBHERecHitCollection>(edm::InputTag(modnam, labelHBHE_, prdnam));
     edm::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_;
   }
  
   tok_ddrec_ = esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>();
   tok_bFieldH_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
   tok_ecalChStatus_ = esConsumes<EcalChannelStatus, EcalChannelStatusRcd>();
   tok_sevlv_ = esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>();
   tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
   tok_caloTopology_ = esConsumes<CaloTopology, CaloTopologyRecord>();
   tok_htopo_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
   tok_resp_ = esConsumes<HcalRespCorrs, HcalRespCorrsRcd>();
  
   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 L1TrigName_    " << l1TrigName_ << "\nThreshold for EB " << hitEthrEB_ << " EE "
       << hitEthrEE0_ << ":" << hitEthrEE1_ << ":" << hitEthrEE2_ << ":" << hitEthrEE3_ << ":" << hitEthrEELo_ << ":"
       << hitEthrEEHi_;
   edm::LogVerbatim("HcalIsoTrack") << "Process " << processName_ << " L1Filter:" << l1Filter_
                                    << " L2Filter:" << l2Filter_ << " L3Filter:" << l3Filter_;
   for (unsigned int k = 0; k < trigNames_.size(); ++k) {
     edm::LogVerbatim("HcalIsoTrack") << "Trigger[" << k << "] " << trigNames_[k];
   }
   edm::LogVerbatim("HcalIsoTrack") << oldID_.size() << " DetIDs to be corrected with HEP17 flag:" << hep17_;
   for (unsigned int k = 0; k < oldID_.size(); ++k)
     edm::LogVerbatim("HcalIsoTrack") << "[" << k << "] Det " << oldDet_[k] << " EtaAbs " << oldEta_[k] << " Depth "
                                      << oldDepth_[k] << ":" << newDepth_[k];
  
   for (int i = 0; i < 10; i++)
     phibins_.push_back(-M_PI + 0.1 * (2 * i + 1) * M_PI);
   for (int i = 0; i < 8; ++i)
     etabins_.push_back(-2.1 + 0.6 * i);
   etadist_ = etabins_[1] - etabins_[0];
   phidist_ = phibins_[1] - phibins_[0];
   etahalfdist_ = 0.5 * etadist_;
   phihalfdist_ = 0.5 * phidist_;
   edm::LogVerbatim("HcalIsoTrack") << "EtaDist " << etadist_ << " " << etahalfdist_ << " PhiDist " << phidist_ << " "
                                    << phihalfdist_;
   unsigned int k1(0), k2(0);
   for (auto phi : phibins_) {
     edm::LogVerbatim("HcalIsoTrack") << "phibin_[" << k1 << "] " << phi;
     ++k1;
   }
   for (auto eta : etabins_) {
     edm::LogVerbatim("HcalIsoTrack") << "etabin_[" << k2 << "] " << eta;
     ++k2;
   }
 }

◆ ~HcalIsoTrkAnalyzer()

HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer ( )

inlineoverride

Definition at line 90 of file HcalIsoTrkAnalyzer.cc.

90 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 386 of file HcalIsoTrkAnalyzer.cc.

                                                                                     {
   t_Run = iEvent.id().run();
   t_Event = iEvent.id().event();
   t_DataType = dataType_;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "Run " << t_Run << " Event " << t_Event << " type " << t_DataType
                                    << " Luminosity " << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing();
 #endif
   //Get magnetic field and ECAL channel status
   const MagneticField* bField = &iSetup.getData(tok_bFieldH_);
   const EcalChannelStatus* theEcalChStatus = &iSetup.getData(tok_ecalChStatus_);
   const EcalSeverityLevelAlgo* theEcalSevlv = &iSetup.getData(tok_sevlv_);
  
   // get handles to calogeometry and calotopology
   const CaloGeometry* geo = &iSetup.getData(tok_geom_);
   const CaloTopology* caloTopology = &iSetup.getData(tok_caloTopology_);
   const HcalTopology* theHBHETopology = &iSetup.getData(tok_htopo_);
   const HcalRespCorrs* resp = &iSetup.getData(tok_resp_);
   HcalRespCorrs* respCorrs = new HcalRespCorrs(*resp);
   respCorrs->setTopo(theHBHETopology);
  
   //=== genParticle information
   edm::Handle<reco::GenParticleCollection> genParticles;
   iEvent.getByToken(tok_parts_, genParticles);
  
   bool okC(true);
   //Get track collection
   edm::Handle<reco::TrackCollection> trkCollection;
   iEvent.getByToken(tok_genTrack_, trkCollection);
   if (!trkCollection.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelGenTrack_;
     okC = false;
   }
  
   //event weight for FLAT sample
   t_EventWeight = 1.0;
   edm::Handle<GenEventInfoProduct> genEventInfo;
   iEvent.getByToken(tok_ew_, genEventInfo);
   if (genEventInfo.isValid())
     t_EventWeight = genEventInfo->weight();
  
   //Define the best vertex and the beamspot
   edm::Handle<reco::VertexCollection> recVtxs;
   iEvent.getByToken(tok_recVtx_, recVtxs);
   edm::Handle<reco::BeamSpot> beamSpotH;
   iEvent.getByToken(tok_bs_, beamSpotH);
   math::XYZPoint leadPV(0, 0, 0);
   t_goodPV = t_nVtx = 0;
   if (recVtxs.isValid() && !(recVtxs->empty())) {
     t_nVtx = recVtxs->size();
     for (unsigned int k = 0; k < recVtxs->size(); ++k) {
       if (!((*recVtxs)[k].isFake()) && ((*recVtxs)[k].ndof() > 4)) {
         if (t_goodPV == 0)
           leadPV = math::XYZPoint((*recVtxs)[k].x(), (*recVtxs)[k].y(), (*recVtxs)[k].z());
         t_goodPV++;
       }
     }
   }
   if (t_goodPV == 0 && beamSpotH.isValid()) {
     leadPV = beamSpotH->position();
   }
   t_allvertex = t_goodPV;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "Primary Vertex " << leadPV << " out of " << t_goodPV << " vertex";
   if (beamSpotH.isValid()) {
     edm::LogVerbatim("HcalIsoTrack") << " Beam Spot " << beamSpotH->position();
   }
 #endif
   // RecHits
   edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
   if (!barrelRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEB_;
     okC = false;
   }
   edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
   if (!endcapRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEE_;
     okC = false;
   }
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(tok_hbhe_, hbhe);
   if (!hbhe.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelHBHE_;
     okC = false;
   }
   edm::Handle<CaloTowerCollection> caloTower;
   iEvent.getByToken(tok_cala_, caloTower);
  
   //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
   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;
  
   //L1
   l1GtUtils_->retrieveL1(iEvent, iSetup, tok_alg_);
   const std::vector<std::pair<std::string, bool>>& finalDecisions = l1GtUtils_->decisionsFinal();
   for (const auto& decision : finalDecisions) {
     if (decision.first.find(l1TrigName_) != std::string::npos) {
       t_L1Bit = decision.second;
       break;
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "Trigger Information for " << l1TrigName_ << " is " << t_L1Bit
                                    << " from a list of " << finalDecisions.size() << " decisions";
 #endif
  
   //HLT
   edm::Handle<edm::TriggerResults> triggerResults;
   iEvent.getByToken(tok_trigRes_, triggerResults);
   if (triggerResults.isValid()) {
     const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
     const std::vector<std::string>& names = triggerNames.triggerNames();
     if (!trigNames_.empty()) {
       for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
         int hlt = triggerResults->accept(iHLT);
         for (unsigned int i = 0; i < trigNames_.size(); ++i) {
           if (names[iHLT].find(trigNames_[i]) != std::string::npos) {
             t_trgbits->at(i) = (hlt > 0);
             t_hltbits->at(i) = (hlt > 0);
             if (hlt > 0)
               t_TrigPass = true;
 #ifdef EDM_ML_DEBUG
             edm::LogVerbatim("HcalIsoTrack")
                 << "This trigger " << names[iHLT] << " Flag " << hlt << ":" << t_trgbits->at(i);
 #endif
           }
         }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "HLT Information shows " << t_TrigPass << ":" << trigNames_.empty() << ":" << okC;
 #endif
  
   std::array<int, 3> ntksave{{0, 0, 0}};
   if (ignoreTrigger_ || useL1Trigger_) {
     t_l1pt = t_l1eta = t_l1phi = 0;
     t_l3pt = t_l3eta = t_l3phi = 0;
     if (ignoreTrigger_ || t_L1Bit)
       ntksave = fillTree(vecL1,
                          vecL3,
                          leadPV,
                          trkCaloDirections,
                          trkCaloDets,
                          geo,
                          caloTopology,
                          theHBHETopology,
                          theEcalChStatus,
                          theEcalSevlv,
                          barrelRecHitsHandle,
                          endcapRecHitsHandle,
                          hbhe,
                          caloTower,
                          genParticles,
                          respCorrs);
     t_TracksSaved = ntksave[0];
     t_TracksLoose = ntksave[1];
     t_TracksTight = ntksave[2];
   } else {
     trigger::TriggerEvent triggerEvent;
     edm::Handle<trigger::TriggerEvent> triggerEventHandle;
     iEvent.getByToken(tok_trigEvt_, triggerEventHandle);
     if (!triggerEventHandle.isValid()) {
       edm::LogWarning("HcalIsoTrack") << "Error! Can't get the product " << triggerEvent_.label();
     } else if (okC) {
       triggerEvent = *(triggerEventHandle.product());
       const trigger::TriggerObjectCollection& TOC(triggerEvent.getObjects());
       bool done(false);
       if (triggerResults.isValid()) {
         std::vector<std::string> modules;
         const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
         const std::vector<std::string>& names = triggerNames.triggerNames();
         for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
           bool ok = (t_TrigPass) || (trigNames_.empty());
           if (ok) {
             unsigned int triggerindx = hltConfig_.triggerIndex(names[iHLT]);
             const std::vector<std::string>& moduleLabels(hltConfig_.moduleLabels(triggerindx));
             std::vector<math::XYZTLorentzVector> vecL2;
             vecL1.clear();
             vecL3.clear();
             //loop over all trigger filters in event (i.e. filters passed)
             for (unsigned int ifilter = 0; ifilter < triggerEvent.sizeFilters(); ++ifilter) {
               std::vector<int> Keys;
               std::string label = triggerEvent.filterTag(ifilter).label();
               //loop over keys to objects passing this filter
               for (unsigned int imodule = 0; imodule < moduleLabels.size(); imodule++) {
                 if (label.find(moduleLabels[imodule]) != std::string::npos) {
 #ifdef EDM_ML_DEBUG
                   edm::LogVerbatim("HcalIsoTrack") << "FilterName " << label;
 #endif
                   for (unsigned int ifiltrKey = 0; ifiltrKey < triggerEvent.filterKeys(ifilter).size(); ++ifiltrKey) {
                     Keys.push_back(triggerEvent.filterKeys(ifilter)[ifiltrKey]);
                     const trigger::TriggerObject& TO(TOC[Keys[ifiltrKey]]);
                     math::XYZTLorentzVector v4(TO.px(), TO.py(), TO.pz(), TO.energy());
                     if (label.find(l2Filter_) != std::string::npos) {
                       vecL2.push_back(v4);
                     } else if (label.find(l3Filter_) != std::string::npos) {
                       vecL3.push_back(v4);
                     } else if ((label.find(l1Filter_) != std::string::npos) || (l1Filter_.empty())) {
                       vecL1.push_back(v4);
                     }
 #ifdef EDM_ML_DEBUG
                     edm::LogVerbatim("HcalIsoTrack")
                         << "key " << ifiltrKey << " : pt " << TO.pt() << " eta " << TO.eta() << " phi " << TO.phi()
                         << " mass " << TO.mass() << " Id " << TO.id();
 #endif
                   }
 #ifdef EDM_ML_DEBUG
                   edm::LogVerbatim("HcalIsoTrack")
                       << "sizes " << vecL1.size() << ":" << vecL2.size() << ":" << vecL3.size();
 #endif
                 }
               }
             }
             math::XYZTLorentzVector mindRvec1;
             double mindR1(999);
             for (unsigned int i = 0; i < vecL2.size(); i++) {
               double dr = dR(vecL1[0], vecL2[i]);
 #ifdef EDM_ML_DEBUG
               edm::LogVerbatim("HcalIsoTrack") << "lvl2[" << i << "] dR " << dr;
 #endif
               if (dr < mindR1) {
                 mindR1 = dr;
                 mindRvec1 = vecL2[i];
               }
             }
 #ifdef EDM_ML_DEBUG
             edm::LogVerbatim("HcalIsoTrack") << "L2 object closest to L1 " << mindRvec1 << " at Dr " << mindR1;
 #endif
  
             if (!vecL1.empty()) {
               t_l1pt = vecL1[0].pt();
               t_l1eta = vecL1[0].eta();
               t_l1phi = vecL1[0].phi();
             } else {
               t_l1pt = t_l1eta = t_l1phi = 0;
             }
             if (!vecL3.empty()) {
               t_l3pt = vecL3[0].pt();
               t_l3eta = vecL3[0].eta();
               t_l3phi = vecL3[0].phi();
             } else {
               t_l3pt = t_l3eta = t_l3phi = 0;
             }
             // Now fill in the tree for each selected track
             if (!done) {
               ntksave = fillTree(vecL1,
                                  vecL3,
                                  leadPV,
                                  trkCaloDirections,
                                  trkCaloDets,
                                  geo,
                                  caloTopology,
                                  theHBHETopology,
                                  theEcalChStatus,
                                  theEcalSevlv,
                                  barrelRecHitsHandle,
                                  endcapRecHitsHandle,
                                  hbhe,
                                  caloTower,
                                  genParticles,
                                  respCorrs);
               t_TracksSaved += ntksave[0];
               t_TracksLoose += ntksave[1];
               t_TracksTight += ntksave[2];
               done = true;
             }
           }
         }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalIsoTrack") << "Final results on selected tracks " << t_TracksSaved << ":" << t_TracksLoose
                                    << ":" << t_TracksTight;
 #endif
   t_TrigPassSel = (t_TracksSaved > 0);
   tree2->Fill();
 }

◆ beginJob()

void HcalIsoTrkAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 692 of file HcalIsoTrkAnalyzer.cc.

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

References fs, TFileService::make(), 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.

◆ beginRun()

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

overrideprivate

Definition at line 773 of file HcalIsoTrkAnalyzer.cc.

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

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

◆ dR()

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

private

Definition at line 1262 of file HcalIsoTrkAnalyzer.cc.

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

References reco::deltaR().

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

◆ endRun()

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

overrideprivate

Definition at line 801 of file HcalIsoTrkAnalyzer.cc.

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

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

◆ eThreshold()

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

private

Definition at line 1313 of file HcalIsoTrkAnalyzer.cc.

                                                                                     {
   const GlobalPoint& pos = geo->getPosition(id);
   double eta = std::abs(pos.eta());
   double eThr(hitEthrEB_);
   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;
 }

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

Referenced by fillTree().

◆ fillDescriptions()

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

static

Definition at line 806 of file HcalIsoTrkAnalyzer.cc.

                                                                                     {
   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", 20.0);
   desc.add<double>("coneRadius", 34.98);
   // signal zone in ECAL and MIP energy cutoff
   desc.add<double>("coneRadiusMIP", 14.0);
   desc.add<double>("coneRadiusMIP2", 18.0);
   desc.add<double>("coneRadiusMIP3", 20.0);
   desc.add<double>("coneRadiusMIP4", 22.0);
   desc.add<double>("coneRadiusMIP5", 24.0);
   desc.add<double>("maximumEcalEnergy", 2.0);
   // following 4 parameters are for isolation cuts and described in the code
   desc.add<double>("maxTrackP", 8.0);
   desc.add<double>("slopeTrackP", 0.05090504066);
   desc.add<double>("isolationEnergyTight", 2.0);
   desc.add<double>("isolationEnergyLoose", 10.0);
   // energy thershold for ECAL (from Egamma group)
   desc.add<double>("EBHitEnergyThreshold", 0.08);
   desc.add<double>("EEHitEnergyThreshold0", 0.30);
   desc.add<double>("EEHitEnergyThreshold1", 0.00);
   desc.add<double>("EEHitEnergyThreshold2", 0.00);
   desc.add<double>("EEHitEnergyThreshold3", 0.00);
   desc.add<double>("EEHitEnergyThresholdLow", 0.30);
   desc.add<double>("EEHitEnergyThresholdHigh", 0.30);
   // prescale factors
   desc.add<double>("momentumLow", 40.0);
   desc.add<double>("momentumHigh", 60.0);
   desc.add<int>("prescaleLow", 1);
   desc.add<int>("prescaleHigh", 1);
   // various labels for collections used in the code
   desc.add<edm::InputTag>("labelTriggerEvent", edm::InputTag("hltTriggerSummaryAOD", "", "HLT"));
   desc.add<edm::InputTag>("labelTriggerResult", edm::InputTag("TriggerResults", "", "HLT"));
   desc.add<std::string>("labelTrack", "generalTracks");
   desc.add<std::string>("labelVertex", "offlinePrimaryVertices");
   desc.add<std::string>("labelEBRecHit", "EcalRecHitsEB");
   desc.add<std::string>("labelEERecHit", "EcalRecHitsEE");
   desc.add<std::string>("labelHBHERecHit", "hbhereco");
   desc.add<std::string>("labelBeamSpot", "offlineBeamSpot");
   desc.add<std::string>("labelCaloTower", "towerMaker");
   desc.add<edm::InputTag>("algInputTag", edm::InputTag("gtStage2Digis"));
   desc.add<edm::InputTag>("extInputTag", edm::InputTag("gtStage2Digis"));
   desc.addUntracked<std::string>("moduleName", "");
   desc.addUntracked<std::string>("producerName", "");
   //  Various flags used for selecting tracks, choice of energy Method2/0
   //  Data type 0/1 for single jet trigger or others
   desc.addUntracked<int>("useRaw", 0);
   desc.addUntracked<bool>("ignoreTriggers", false);
   desc.addUntracked<bool>("useL1Trigger", false);
   desc.addUntracked<double>("hcalScale", 1.0);
   desc.addUntracked<int>("dataType", 0);
   desc.addUntracked<bool>("unCorrect", false);
   desc.addUntracked<bool>("collapseDepth", false);
   desc.addUntracked<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);
   descriptions.add("HcalIsoTrkAnalyzer", desc);
 }

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

◆ 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 HcalTopology *	theHBHETopology,
		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 HcalRespCorrs *	respCorrs
	)

private

Definition at line 894 of file HcalIsoTrkAnalyzer.cc.

                                                                                 {
   int nSave(0), nLoose(0), nTight(0);
   //Loop over tracks
   std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
   unsigned int nTracks(0), nselTracks(0);
   t_nTrk = trkCaloDirections.size();
   t_rhoh = (tower.isValid()) ? rhoh(tower) : 0;
   for (trkDetItr = trkCaloDirections.begin(), nTracks = 0; trkDetItr != trkCaloDirections.end();
        trkDetItr++, nTracks++) {
     const reco::Track* pTrack = &(*(trkDetItr->trkItr));
     math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) :" << pTrack->pt() << "|"
                                      << pTrack->eta() << "|" << pTrack->phi() << "|" << pTrack->p();
 #endif
     t_mindR2 = 999;
     for (unsigned int k = 0; k < vecL3.size(); ++k) {
       double dr = dR(vecL3[k], v4);
       if (dr < t_mindR2) {
         t_mindR2 = dr;
       }
     }
     t_mindR1 = (!vecL1.empty()) ? dR(vecL1[0], v4) : 999;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalIsoTrack") << "Closest L3 object at dr :" << t_mindR2 << " and from L1 " << t_mindR1;
 #endif
     t_ieta = t_iphi = 0;
     if (trkDetItr->okHCAL) {
       HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
       t_ieta = detId.ieta();
       t_iphi = detId.iphi();
       if (t_p > 40.0 && t_p <= 60.0)
         t_ietaAll->emplace_back(t_ieta);
     }
     //Selection of good track
     t_selectTk = 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
     edm::LogVerbatim("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " << qltyFlag << "|" << trkDetItr->okECAL << "|"
                                      << trkDetItr->okHCAL << " eIsolation " << eIsolation;
 #endif
     t_qltyFlag = (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL);
     if (t_qltyFlag) {
       nselTracks++;
       int nNearTRKs(0);
       std::vector<DetId> eIds;
       std::vector<double> eHit;
       t_eMipDR = spr::eCone_ecal(geo,
                                  barrelRecHitsHandle,
                                  endcapRecHitsHandle,
                                  trkDetItr->pointHCAL,
                                  trkDetItr->pointECAL,
                                  a_mipR_,
                                  trkDetItr->directionECAL,
                                  eIds,
                                  eHit);
       double eEcal(0);
       for (unsigned int k = 0; k < eIds.size(); ++k) {
         if (eHit[k] > eThreshold(eIds[k], geo))
           eEcal += eHit[k];
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR << ":" << eEcal;
 #endif
       t_eMipDR = eEcal;
       std::vector<DetId> eIds2;
       std::vector<double> eHit2;
       t_eMipDR2 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_mipR2_,
                                   trkDetItr->directionECAL,
                                   eIds2,
                                   eHit2);
       double eEcal2(0);
       for (unsigned int k = 0; k < eIds2.size(); ++k) {
         if (eHit2[k] > eThreshold(eIds2[k], geo))
           eEcal2 += eHit2[k];
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR2 << ":" << eEcal2;
 #endif
       t_eMipDR2 = eEcal2;
       std::vector<DetId> eIds3;
       std::vector<double> eHit3;
       t_eMipDR3 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_mipR3_,
                                   trkDetItr->directionECAL,
                                   eIds3,
                                   eHit3);
       double eEcal3(0);
       for (unsigned int k = 0; k < eIds3.size(); ++k) {
         if (eHit3[k] > eThreshold(eIds3[k], geo))
           eEcal3 += eHit3[k];
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR3 << ":" << eEcal3;
 #endif
       t_eMipDR3 = eEcal3;
       std::vector<DetId> eIds4;
       std::vector<double> eHit4;
       t_eMipDR4 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_mipR4_,
                                   trkDetItr->directionECAL,
                                   eIds4,
                                   eHit4);
       double eEcal4(0);
       for (unsigned int k = 0; k < eIds4.size(); ++k) {
         if (eHit4[k] > eThreshold(eIds4[k], geo))
           eEcal4 += eHit4[k];
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR4 << ":" << eEcal4;
 #endif
       t_eMipDR4 = eEcal4;
       std::vector<DetId> eIds5;
       std::vector<double> eHit5;
       t_eMipDR5 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_mipR5_,
                                   trkDetItr->directionECAL,
                                   eIds5,
                                   eHit5);
       double eEcal5(0);
       for (unsigned int k = 0; k < eIds5.size(); ++k) {
         if (eHit5[k] > eThreshold(eIds5[k], geo))
           eEcal5 += eHit5[k];
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR5 << ":" << eEcal5;
 #endif
       t_eMipDR5 = eEcal5;
  
       t_emaxNearP = spr::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
       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, respCorrs, ids, edet0, t_eHcal, t_DetIds, t_HitEnergies);
  
         //----- hcal energy in the extended cone 1 (a_coneR+10) --------------
         t_eHcal10 = spr::eCone_hcal(geo,
                                     hbhe,
                                     trkDetItr->pointHCAL,
                                     trkDetItr->pointECAL,
                                     a_coneR1_,
                                     trkDetItr->directionHCAL,
                                     nRecHits1,
                                     ids1,
                                     edet1,
                                     useRaw_);
         if (!oldID_.empty()) {
           for (unsigned k = 0; k < ids1.size(); ++k)
             ids1[k] = newId(ids1[k]);
         }
         storeEnergy(1, respCorrs, ids1, edet1, t_eHcal10, t_DetIds1, t_HitEnergies1);
  
         //----- hcal energy in the extended cone 3 (a_coneR+30) --------------
         t_eHcal30 = spr::eCone_hcal(geo,
                                     hbhe,
                                     trkDetItr->pointHCAL,
                                     trkDetItr->pointECAL,
                                     a_coneR2_,
                                     trkDetItr->directionHCAL,
                                     nRecHits3,
                                     ids3,
                                     edet3,
                                     useRaw_);
         if (!oldID_.empty()) {
           for (unsigned k = 0; k < ids3.size(); ++k)
             ids3[k] = newId(ids3[k]);
         }
         storeEnergy(3, respCorrs, ids3, edet3, t_eHcal30, t_DetIds3, t_HitEnergies3);
  
         t_p = pTrack->p();
         t_pt = pTrack->pt();
         t_phi = pTrack->phi();
  
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) :" << t_pt << "|"
                                          << pTrack->eta() << "|" << t_phi << "|" << t_p << " Generator Level p "
                                          << t_gentrackP;
         edm::LogVerbatim("HcalIsoTrack") << "e_MIP " << t_eMipDR << " Chg Isolation " << t_hmaxNearP << " eHcal"
                                          << t_eHcal << " ieta " << t_ieta << " Quality " << t_qltyMissFlag << ":"
                                          << t_qltyPVFlag << ":" << t_selectTk;
         for (unsigned int ll = 0; ll < t_DetIds->size(); ll++) {
           edm::LogVerbatim("HcalIsoTrack")
               << "det id is = " << t_DetIds->at(ll) << "   hit enery is  = " << t_HitEnergies->at(ll);
         }
         for (unsigned int ll = 0; ll < t_DetIds1->size(); ll++) {
           edm::LogVerbatim("HcalIsoTrack")
               << "det id is = " << t_DetIds1->at(ll) << "   hit enery is  = " << t_HitEnergies1->at(ll);
         }
         for (unsigned int ll = 0; ll < t_DetIds3->size(); ll++) {
           edm::LogVerbatim("HcalIsoTrack")
               << "det id is = " << t_DetIds3->at(ll) << "   hit enery is  = " << t_HitEnergies3->at(ll);
         }
 #endif
         bool accept(false);
         if (t_p > pTrackMin_) {
           if (t_p < pTrackLow_) {
             ++nLow_;
             if (prescaleLow_ <= 1)
               accept = true;
             else if (nLow_ % prescaleLow_ == 1)
               accept = true;
           } else if (t_p > pTrackHigh_) {
             ++nHigh_;
             if (prescaleHigh_ <= 1)
               accept = true;
             else if (nHigh_ % prescaleHigh_ == 1)
               accept = true;
           } else {
             accept = true;
           }
         }
         if (accept) {
           tree->Fill();
           nSave++;
           int type(0);
           if (t_eMipDR < 1.0) {
             if (t_hmaxNearP < eIsolate2_) {
               ++nLoose;
               type = 1;
             }
             if (t_hmaxNearP < eIsolate1_) {
               ++nTight;
               type = 2;
             }
           }
           if (t_p > 40.0 && t_p <= 60.0 && t_selectTk) {
             t_ietaGood->emplace_back(t_ieta);
             t_trackType->emplace_back(type);
           }
 #ifdef EDM_ML_DEBUG
           for (unsigned int k = 0; k < t_trgbits->size(); k++) {
             edm::LogVerbatim("HcalIsoTrack") << "trigger bit is  = " << t_trgbits->at(k);
           }
 #endif
         }
       }
     }
   }
   std::array<int, 3> i3{{nSave, nLoose, nTight}};
   return i3;
 }

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(), flavorHistoryFilter_cfi::dr, dR(), spr::eCone_ecal(), spr::eCone_hcal(), spr::eECALmatrix(), eEcalMax_, spr::eHCALmatrix(), eIsolate1_, eIsolate2_, reco::TrackBase::eta(), eThreshold(), JetChargeProducer_cfi::exp, genParticles2HepMC_cfi::genParticles, spr::goodTrack(), photonIsolationHIProducer_cfi::hbhe, testProducerWithPsetDescEmpty_cfi::i3, HcalDetId::ieta(), HcalDetId::iphi(), dqmdumpme::k, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, maxRestrictionP_, newId(), nHigh_, jetSelector_cfi::nLoose, nLow_, 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_gentrackP, t_hAnnular, t_HitEnergies, t_HitEnergies1, t_HitEnergies3, t_hmaxNearP, t_ieta, t_ietaAll, t_ietaGood, t_iphi, t_mindR1, t_mindR2, t_nTrk, t_p, t_phi, t_pt, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_rhoh, t_selectTk, t_trackType, t_trgbits, hgcalTowerProducer_cfi::tower, trackP(), and useRaw_.

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

◆ newId()

DetId HcalIsoTrkAnalyzer::newId ( const DetId & id )

private

Definition at line 1327 of file HcalIsoTrkAnalyzer.cc.

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

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

Referenced by fillTree().

◆ rhoh()

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

private

Definition at line 1282 of file HcalIsoTrkAnalyzer.cc.

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

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

Referenced by fillTree().

◆ storeEnergy()

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

private

Definition at line 1339 of file HcalIsoTrkAnalyzer.cc.

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

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

Referenced by fillTree().

◆ trackP()

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

private

Definition at line 1266 of file HcalIsoTrkAnalyzer.cc.

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

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

Referenced by fillTree().

Member Data Documentation

◆ a_charIsoR_

double HcalIsoTrkAnalyzer::a_charIsoR_

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_coneR1_

double HcalIsoTrkAnalyzer::a_coneR1_

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_coneR2_

double HcalIsoTrkAnalyzer::a_coneR2_

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_coneR_

const double HcalIsoTrkAnalyzer::a_coneR_

private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR2_

const double HcalIsoTrkAnalyzer::a_mipR2_

private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR3_

const double HcalIsoTrkAnalyzer::a_mipR3_

private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR4_

const double HcalIsoTrkAnalyzer::a_mipR4_

private

Definition at line 138 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR5_

const double HcalIsoTrkAnalyzer::a_mipR5_

private

Definition at line 138 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR_

const double HcalIsoTrkAnalyzer::a_mipR_

private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ collapseDepth_

const bool HcalIsoTrkAnalyzer::collapseDepth_

private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

◆ dataType_

const int HcalIsoTrkAnalyzer::dataType_

private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ eEcalMax_

const double HcalIsoTrkAnalyzer::eEcalMax_

private

Definition at line 138 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ eIsolate1_

const double HcalIsoTrkAnalyzer::eIsolate1_

private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ eIsolate2_

const double HcalIsoTrkAnalyzer::eIsolate2_

private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ etabins_

std::vector<double> HcalIsoTrkAnalyzer::etabins_

private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ etadist_

double HcalIsoTrkAnalyzer::etadist_

private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ etahalfdist_

double HcalIsoTrkAnalyzer::etahalfdist_

private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ fs

edm::Service<TFileService> HcalIsoTrkAnalyzer::fs

private

Definition at line 130 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob().

◆ hcalScale_

const double HcalIsoTrkAnalyzer::hcalScale_

private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

◆ hdc_

const HcalDDDRecConstants* HcalIsoTrkAnalyzer::hdc_

private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and storeEnergy().

◆ hep17_

const bool HcalIsoTrkAnalyzer::hep17_

private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ hitEthrEB_

const double HcalIsoTrkAnalyzer::hitEthrEB_

private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE0_

const double HcalIsoTrkAnalyzer::hitEthrEE0_

private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE1_

const double HcalIsoTrkAnalyzer::hitEthrEE1_

private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE2_

const double HcalIsoTrkAnalyzer::hitEthrEE2_

private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE3_

const double HcalIsoTrkAnalyzer::hitEthrEE3_

private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEEHi_

const double HcalIsoTrkAnalyzer::hitEthrEEHi_

private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEELo_

const double HcalIsoTrkAnalyzer::hitEthrEELo_

private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hltConfig_

HLTConfigProvider HcalIsoTrkAnalyzer::hltConfig_

private

Definition at line 131 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginRun().

◆ ignoreTrigger_

const bool HcalIsoTrkAnalyzer::ignoreTrigger_

private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l1Filter_

const std::string HcalIsoTrkAnalyzer::l1Filter_

private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l1GtUtils_

l1t::L1TGlobalUtil* HcalIsoTrkAnalyzer::l1GtUtils_

private

Definition at line 129 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l1TrigName_

const std::string HcalIsoTrkAnalyzer::l1TrigName_

private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l2Filter_

const std::string HcalIsoTrkAnalyzer::l2Filter_

private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l3Filter_

const std::string HcalIsoTrkAnalyzer::l3Filter_

private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelEB_

const std::string HcalIsoTrkAnalyzer::labelEB_

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelEE_

const std::string HcalIsoTrkAnalyzer::labelEE_

private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelGenTrack_

const std::string HcalIsoTrkAnalyzer::labelGenTrack_

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelHBHE_

const std::string HcalIsoTrkAnalyzer::labelHBHE_

private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelRecVtx_

const std::string HcalIsoTrkAnalyzer::labelRecVtx_

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ labelTower_

const std::string HcalIsoTrkAnalyzer::labelTower_

private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ maxRestrictionP_

const double HcalIsoTrkAnalyzer::maxRestrictionP_

private

Definition at line 139 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ mode_

const int HcalIsoTrkAnalyzer::mode_

private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ newDepth_

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

private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ nHigh_

unsigned int HcalIsoTrkAnalyzer::nHigh_

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree().

◆ nLow_

unsigned int HcalIsoTrkAnalyzer::nLow_

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree().

◆ nRun_

unsigned int HcalIsoTrkAnalyzer::nRun_

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and endRun().

◆ oldDepth_

std::vector<int> HcalIsoTrkAnalyzer::oldDepth_

private

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ oldDet_

std::vector<int> HcalIsoTrkAnalyzer::oldDet_

private

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ oldEta_

std::vector<int> HcalIsoTrkAnalyzer::oldEta_

private

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ oldID_

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

private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

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

◆ phibins_

std::vector<double> HcalIsoTrkAnalyzer::phibins_

private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ phidist_

double HcalIsoTrkAnalyzer::phidist_

private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ phihalfdist_

double HcalIsoTrkAnalyzer::phihalfdist_

private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ prescaleHigh_

const int HcalIsoTrkAnalyzer::prescaleHigh_

private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ prescaleLow_

const int HcalIsoTrkAnalyzer::prescaleLow_

private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ processName_

const std::string HcalIsoTrkAnalyzer::processName_

private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkAnalyzer().

◆ pTrackHigh_

const double HcalIsoTrkAnalyzer::pTrackHigh_

private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ pTrackLow_

const double HcalIsoTrkAnalyzer::pTrackLow_

private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ pTrackMin_

const double HcalIsoTrkAnalyzer::pTrackMin_

private

Definition at line 138 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ selectionParameter_

spr::trackSelectionParameters HcalIsoTrkAnalyzer::selectionParameter_

private

Definition at line 133 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ slopeRestrictionP_

const double HcalIsoTrkAnalyzer::slopeRestrictionP_

private

Definition at line 139 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ t_allvertex

int HcalIsoTrkAnalyzer::t_allvertex

private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_DataType

int HcalIsoTrkAnalyzer::t_DataType

private

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_DetIds

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

private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_DetIds1

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

private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_DetIds3

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

private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eAnnular

double HcalIsoTrkAnalyzer::t_eAnnular

private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal

double HcalIsoTrkAnalyzer::t_eHcal

private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal10

double HcalIsoTrkAnalyzer::t_eHcal10

private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal30

double HcalIsoTrkAnalyzer::t_eHcal30

private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_emaxNearP

double HcalIsoTrkAnalyzer::t_emaxNearP

private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR

double HcalIsoTrkAnalyzer::t_eMipDR

private

Definition at line 190 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR2

double HcalIsoTrkAnalyzer::t_eMipDR2

private

Definition at line 190 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR3

double HcalIsoTrkAnalyzer::t_eMipDR3

private

Definition at line 190 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR4

double HcalIsoTrkAnalyzer::t_eMipDR4

private

Definition at line 190 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR5

double HcalIsoTrkAnalyzer::t_eMipDR5

private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Event

int HcalIsoTrkAnalyzer::t_Event

private

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_EventNo

unsigned int HcalIsoTrkAnalyzer::t_EventNo

private

Definition at line 183 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_EventWeight

double HcalIsoTrkAnalyzer::t_EventWeight

private

Definition at line 186 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_gentrackP

double HcalIsoTrkAnalyzer::t_gentrackP

private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_goodPV

int HcalIsoTrkAnalyzer::t_goodPV

private

Definition at line 185 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_hAnnular

double HcalIsoTrkAnalyzer::t_hAnnular

private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies

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

private

Definition at line 197 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies1

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

private

Definition at line 197 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies3

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

private

Definition at line 197 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_hltbits

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

private

Definition at line 198 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_hmaxNearP

double HcalIsoTrkAnalyzer::t_hmaxNearP

private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_ieta

int HcalIsoTrkAnalyzer::t_ieta

private

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_ietaAll

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

private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_ietaGood

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

private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_iphi

int HcalIsoTrkAnalyzer::t_iphi

private

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_L1Bit

bool HcalIsoTrkAnalyzer::t_L1Bit

private

Definition at line 199 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1eta

double HcalIsoTrkAnalyzer::t_l1eta

private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1phi

double HcalIsoTrkAnalyzer::t_l1phi

private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1pt

double HcalIsoTrkAnalyzer::t_l1pt

private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3eta

double HcalIsoTrkAnalyzer::t_l3eta

private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3phi

double HcalIsoTrkAnalyzer::t_l3phi

private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3pt

double HcalIsoTrkAnalyzer::t_l3pt

private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_mindR1

double HcalIsoTrkAnalyzer::t_mindR1

private

Definition at line 189 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_mindR2

double HcalIsoTrkAnalyzer::t_mindR2

private

Definition at line 189 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_nTrk

int HcalIsoTrkAnalyzer::t_nTrk

private

Definition at line 185 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_nVtx

int HcalIsoTrkAnalyzer::t_nVtx

private

Definition at line 185 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_p

double HcalIsoTrkAnalyzer::t_p

private

Definition at line 186 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_phi

double HcalIsoTrkAnalyzer::t_phi

private

Definition at line 186 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_pt

double HcalIsoTrkAnalyzer::t_pt

private

Definition at line 186 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyFlag

bool HcalIsoTrkAnalyzer::t_qltyFlag

private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyMissFlag

bool HcalIsoTrkAnalyzer::t_qltyMissFlag

private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyPVFlag

bool HcalIsoTrkAnalyzer::t_qltyPVFlag

private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_rhoh

double HcalIsoTrkAnalyzer::t_rhoh

private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Run

int HcalIsoTrkAnalyzer::t_Run

private

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_RunNo

unsigned int HcalIsoTrkAnalyzer::t_RunNo

private

Definition at line 183 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_selectTk

bool HcalIsoTrkAnalyzer::t_selectTk

private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Tracks

int HcalIsoTrkAnalyzer::t_Tracks

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksLoose

int HcalIsoTrkAnalyzer::t_TracksLoose

private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksProp

int HcalIsoTrkAnalyzer::t_TracksProp

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksSaved

int HcalIsoTrkAnalyzer::t_TracksSaved

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksTight

int HcalIsoTrkAnalyzer::t_TracksTight

private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_trackType

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

private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_trgbits

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

private

Definition at line 198 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_TrigPass

bool HcalIsoTrkAnalyzer::t_TrigPass

private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TrigPassSel

bool HcalIsoTrkAnalyzer::t_TrigPassSel

private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ theTrackQuality_

const std::string HcalIsoTrkAnalyzer::theTrackQuality_

private

Definition at line 134 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ theTriggerResultsLabel_

const edm::InputTag HcalIsoTrkAnalyzer::theTriggerResultsLabel_

private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ tok_alg_

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

private

Definition at line 171 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_bFieldH_

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

private

Definition at line 174 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_bs_

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

private

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_cala_

edm::EDGetTokenT<CaloTowerCollection> HcalIsoTrkAnalyzer::tok_cala_

private

Definition at line 169 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_caloTopology_

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

private

Definition at line 178 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_ddrec_

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

private

Definition at line 173 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkAnalyzer().

◆ tok_EB_

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_EB_

private

Definition at line 166 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_ecalChStatus_

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

private

Definition at line 175 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_EE_

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_EE_

private

Definition at line 167 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_ew_

edm::EDGetTokenT<GenEventInfoProduct> HcalIsoTrkAnalyzer::tok_ew_

private

Definition at line 170 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_genTrack_

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

private

Definition at line 163 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_geom_

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

private

Definition at line 177 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> HcalIsoTrkAnalyzer::tok_hbhe_

private

Definition at line 168 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_htopo_

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

private

Definition at line 179 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_parts_

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

private

Definition at line 162 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_recVtx_

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

private

Definition at line 164 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_resp_

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

private

Definition at line 180 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_sevlv_

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

private

Definition at line 176 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_trigEvt_

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

private

Definition at line 160 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_trigRes_

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

private

Definition at line 161 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tree

TTree* HcalIsoTrkAnalyzer::tree

private

Definition at line 182 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 182 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ triggerEvent_

const edm::InputTag HcalIsoTrkAnalyzer::triggerEvent_

private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ trigNames_

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

private

Definition at line 132 of file HcalIsoTrkAnalyzer.cc.

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

◆ unCorrect_

const bool HcalIsoTrkAnalyzer::unCorrect_

private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

◆ useL1Trigger_

const bool HcalIsoTrkAnalyzer::useL1Trigger_

private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ useRaw_

const int HcalIsoTrkAnalyzer::useRaw_

private

Definition at line 143 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()

◆ newId()

◆ rhoh()

◆ storeEnergy()

◆ trackP()

Member Data Documentation

◆ a_charIsoR_

◆ a_coneR1_

◆ a_coneR2_

◆ a_coneR_

◆ a_mipR2_

◆ a_mipR3_

◆ a_mipR4_

◆ a_mipR5_

◆ a_mipR_

◆ collapseDepth_

◆ dataType_

◆ eEcalMax_

◆ eIsolate1_

◆ eIsolate2_

◆ etabins_

◆ etadist_

◆ etahalfdist_

◆ fs

◆ hcalScale_

◆ hdc_

◆ hep17_

◆ hitEthrEB_

◆ hitEthrEE0_

◆ hitEthrEE1_

◆ hitEthrEE2_

◆ hitEthrEE3_

◆ hitEthrEEHi_

◆ hitEthrEELo_

◆ hltConfig_

◆ ignoreTrigger_

◆ l1Filter_

◆ l1GtUtils_

◆ l1TrigName_

◆ l2Filter_

◆ l3Filter_

◆ labelEB_

◆ labelEE_

◆ labelGenTrack_

◆ labelHBHE_

◆ labelRecVtx_

◆ labelTower_

◆ maxRestrictionP_

◆ mode_

◆ newDepth_

◆ nHigh_

◆ nLow_

◆ nRun_

◆ oldDepth_

◆ oldDet_

◆ oldEta_

◆ oldID_

◆ phibins_

◆ phidist_

◆ phihalfdist_

◆ prescaleHigh_

◆ prescaleLow_

◆ processName_