CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
HcalIsoTrkSimAnalyzer Class Reference
Inheritance diagram for HcalIsoTrkSimAnalyzer:
edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources > edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

 HcalIsoTrkSimAnalyzer (edm::ParameterSet const &)
 
 ~HcalIsoTrkSimAnalyzer () override=default
 
- Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
 EDAnalyzer ()=default
 
 EDAnalyzer (const EDAnalyzer &)=delete
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
const EDAnalyzeroperator= (const EDAnalyzer &)=delete
 
bool wantsGlobalLuminosityBlocks () const noexcept final
 
bool wantsGlobalRuns () const noexcept final
 
bool wantsInputProcessBlocks () const noexcept final
 
bool wantsProcessBlocks () const noexcept final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const noexcept
 
bool wantsStreamRuns () const noexcept
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESResolverIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESResolverIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Member Functions

void analyze (edm::Event const &, edm::EventSetup const &) override
 
void beginJob () override
 
void beginRun (edm::Run const &, edm::EventSetup const &) override
 
double dR (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
 
void endRun (edm::Run const &, edm::EventSetup const &) override
 
double eThreshold (const DetId &id, const CaloGeometry *geo) const
 
std::array< int, 3 > fillTree (std::vector< math::XYZTLorentzVector > &vecL1, std::vector< math::XYZTLorentzVector > &vecL3, math::XYZPoint &leadPV, std::vector< spr::propagatedGenParticleID > &trackIDs, const CaloGeometry *geo, const CaloTopology *topo, const HcalTopology *theHBHETopology, const EcalChannelStatus *theEcalChStatus, const EcalSeverityLevelAlgo *theEcalSevlv, edm::Handle< EcalRecHitCollection > &barrelRecHitsHandle, edm::Handle< EcalRecHitCollection > &endcapRecHitsHandle, edm::Handle< HBHERecHitCollection > &hbhe, const edm::Handle< CaloTowerCollection > &towerHandle, edm::Handle< reco::GenParticleCollection > &genParticles, const HcalRespCorrs *respCorrs)
 
DetId newId (const DetId &)
 
bool notaMuon (const reco::GenParticle *pTrack)
 
double rhoh (const edm::Handle< CaloTowerCollection > &)
 
void storeEnergy (int indx, const HcalRespCorrs *respCorrs, const std::vector< DetId > &ids, std::vector< double > &edet, double &eHcal, std::vector< unsigned int > *detIds, std::vector< double > *hitEnergies)
 

Private Attributes

double a_charIsoR_
 
double a_coneR1_
 
double a_coneR2_
 
const double a_coneR_
 
const double a_mipR2_
 
const double a_mipR3_
 
const double a_mipR4_
 
const double a_mipR5_
 
const double a_mipR_
 
const edm::InputTag algTag_
 
const bool collapseDepth_
 
const int dataType_
 
const double eEcalMax_
 
const double eIsolate1_
 
const double eIsolate2_
 
std::vector< double > etabins_
 
double etadist_
 
double etahalfdist_
 
const double etaMax_
 
const EcalPFRecHitThresholdseThresholds_
 
const edm::InputTag extTag_
 
const bool fillTreeRange_
 
const double hcalScale_
 
const HcalDDDRecConstantshdc_
 
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::L1TGlobalUtill1GtUtils_
 
const std::string l1TrigName_
 
const std::string l2Filter_
 
const std::string l3Filter_
 
const std::string labelBS_
 
const std::string labelEB_
 
const std::string labelEE_
 
const std::string labelGenTrack_
 
const std::string labelHBHE_
 
const std::string labelRecVtx_
 
const std::string labelTower_
 
const double maxRestrictionP_
 
const int mode_
 
const std::string modnam_
 
const std::vector< int > newDepth_
 
unsigned int nHigh_
 
unsigned int nLow_
 
unsigned int nRun_
 
std::vector< int > oldDepth_
 
std::vector< int > oldDet_
 
std::vector< int > oldEta_
 
const std::vector< int > oldID_
 
std::vector< double > phibins_
 
double phidist_
 
double phihalfdist_
 
const std::string prdnam_
 
const int prescaleHigh_
 
const int prescaleLow_
 
const std::string processName_
 
const double ptMin_
 
const double pTrackHigh_
 
const double pTrackLow_
 
const double pTrackMin_
 
const double slopeRestrictionP_
 
int t_allvertex
 
int t_DataType
 
std::vector< unsigned int > * t_DetIds
 
std::vector< unsigned int > * t_DetIds1
 
std::vector< unsigned int > * t_DetIds3
 
double t_eAnnular
 
double t_eHcal
 
double t_eHcal10
 
double t_eHcal30
 
double t_emaxNearP
 
double t_eMipDR
 
double t_eMipDR2
 
double t_eMipDR3
 
double t_eMipDR4
 
double t_eMipDR5
 
int t_Event
 
unsigned int t_EventNo
 
double t_EventWeight
 
double t_gentrackP
 
int t_goodPV
 
double t_hAnnular
 
std::vector< double > * t_HitEnergies
 
std::vector< double > * t_HitEnergies1
 
std::vector< double > * t_HitEnergies3
 
std::vector< bool > * t_hltbits
 
double t_hmaxNearP
 
int t_ieta
 
std::vector< int > * t_ietaAll
 
std::vector< int > * t_ietaGood
 
int t_iphi
 
bool t_L1Bit
 
double t_l1eta
 
double t_l1phi
 
double t_l1pt
 
double t_l3eta
 
double t_l3phi
 
double t_l3pt
 
double t_mindR1
 
double t_mindR2
 
int t_nTrk
 
int t_nVtx
 
double t_p
 
double t_phi
 
double t_pt
 
bool t_qltyFlag
 
bool t_qltyMissFlag
 
bool t_qltyPVFlag
 
double t_rhoh
 
int t_Run
 
unsigned int t_RunNo
 
bool t_selectTk
 
int t_Tracks
 
int t_TracksLoose
 
int t_TracksProp
 
int t_TracksSaved
 
int t_TracksTight
 
std::vector< int > * t_trackType
 
std::vector< bool > * t_trgbits
 
bool t_TrigPass
 
bool t_TrigPassSel
 
const edm::InputTag theTriggerResultsLabel_
 
edm::EDGetTokenT< BXVector< GlobalAlgBlk > > tok_alg_
 
edm::ESGetToken< MagneticField, IdealMagneticFieldRecordtok_bFieldH_
 
edm::EDGetTokenT< reco::BeamSpottok_bs_
 
edm::EDGetTokenT< CaloTowerCollectiontok_cala_
 
edm::ESGetToken< CaloTopology, CaloTopologyRecordtok_caloTopology_
 
edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecordtok_ddrec_
 
edm::EDGetTokenT< EcalRecHitCollectiontok_EB_
 
edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcdtok_ecalChStatus_
 
edm::ESGetToken< EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcdtok_ecalPFRecHitThresholds_
 
edm::EDGetTokenT< EcalRecHitCollectiontok_EE_
 
edm::EDGetTokenT< GenEventInfoProducttok_ew_
 
edm::ESGetToken< CaloGeometry, CaloGeometryRecordtok_geom_
 
edm::EDGetTokenT< HBHERecHitCollectiontok_hbhe_
 
edm::ESGetToken< HcalTopology, HcalRecNumberingRecordtok_htopo_
 
edm::EDGetTokenT< reco::GenParticleCollectiontok_parts_
 
edm::ESGetToken< HepPDT::ParticleDataTable, PDTRecordtok_pdt_
 
edm::EDGetTokenT< reco::VertexCollectiontok_recVtx_
 
edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcdtok_resp_
 
edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcdtok_sevlv_
 
edm::EDGetTokenT< trigger::TriggerEventtok_trigEvt_
 
edm::EDGetTokenT< edm::TriggerResultstok_trigRes_
 
TTree * tree
 
TTree * tree2
 
const edm::InputTag triggerEvent_
 
const std::vector< std::string > trigNames_
 
const bool unCorrect_
 
const bool useL1Trigger_
 
const bool usePFThresh_
 
const int useRaw_
 

Additional Inherited Members

- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< Bconsumes (edm::InputTag tag) noexcept
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Definition at line 77 of file HcalIsoTrkSimAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalIsoTrkSimAnalyzer()

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

Definition at line 200 of file HcalIsoTrkSimAnalyzer.cc.

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

201  : trigNames_(iConfig.getParameter<std::vector<std::string> >("triggers")),
202  processName_(iConfig.getParameter<std::string>("processName")),
203  l1Filter_(iConfig.getParameter<std::string>("l1Filter")),
204  l2Filter_(iConfig.getParameter<std::string>("l2Filter")),
205  l3Filter_(iConfig.getParameter<std::string>("l3Filter")),
206  ptMin_(iConfig.getParameter<double>("pTMin")),
207  etaMax_(iConfig.getParameter<double>("maxChargedHadronEta")),
208  a_coneR_(iConfig.getParameter<double>("coneRadius")),
209  a_mipR_(iConfig.getParameter<double>("coneRadiusMIP")),
210  a_mipR2_(iConfig.getParameter<double>("coneRadiusMIP2")),
211  a_mipR3_(iConfig.getParameter<double>("coneRadiusMIP3")),
212  a_mipR4_(iConfig.getParameter<double>("coneRadiusMIP4")),
213  a_mipR5_(iConfig.getParameter<double>("coneRadiusMIP5")),
214  pTrackMin_(iConfig.getParameter<double>("minimumTrackP")),
215  eEcalMax_(iConfig.getParameter<double>("maximumEcalEnergy")),
216  maxRestrictionP_(iConfig.getParameter<double>("maxTrackP")),
217  slopeRestrictionP_(iConfig.getParameter<double>("slopeTrackP")),
218  hcalScale_(iConfig.getUntrackedParameter<double>("hHcalScale", 1.0)),
219  eIsolate1_(iConfig.getParameter<double>("isolationEnergyTight")),
220  eIsolate2_(iConfig.getParameter<double>("isolationEnergyLoose")),
221  pTrackLow_(iConfig.getParameter<double>("momentumLow")),
222  pTrackHigh_(iConfig.getParameter<double>("momentumHigh")),
223  prescaleLow_(iConfig.getParameter<int>("prescaleLow")),
224  prescaleHigh_(iConfig.getParameter<int>("prescaleHigh")),
225  useRaw_(iConfig.getUntrackedParameter<int>("useRaw", 0)),
226  dataType_(iConfig.getUntrackedParameter<int>("dataType", 0)),
227  mode_(iConfig.getUntrackedParameter<int>("outMode", 11)),
228  ignoreTrigger_(iConfig.getUntrackedParameter<bool>("ignoreTriggers", false)),
229  useL1Trigger_(iConfig.getUntrackedParameter<bool>("useL1Trigger", false)),
230  unCorrect_(iConfig.getUntrackedParameter<bool>("unCorrect", false)),
231  collapseDepth_(iConfig.getUntrackedParameter<bool>("collapseDepth", false)),
232  fillTreeRange_(iConfig.getUntrackedParameter<bool>("fillTreeRange", false)),
233  hitEthrEB_(iConfig.getParameter<double>("EBHitEnergyThreshold")),
234  hitEthrEE0_(iConfig.getParameter<double>("EEHitEnergyThreshold0")),
235  hitEthrEE1_(iConfig.getParameter<double>("EEHitEnergyThreshold1")),
236  hitEthrEE2_(iConfig.getParameter<double>("EEHitEnergyThreshold2")),
237  hitEthrEE3_(iConfig.getParameter<double>("EEHitEnergyThreshold3")),
238  hitEthrEELo_(iConfig.getParameter<double>("EEHitEnergyThresholdLow")),
239  hitEthrEEHi_(iConfig.getParameter<double>("EEHitEnergyThresholdHigh")),
240  triggerEvent_(iConfig.getParameter<edm::InputTag>("labelTriggerEvent")),
241  theTriggerResultsLabel_(iConfig.getParameter<edm::InputTag>("labelTriggerResult")),
242  labelGenTrack_(iConfig.getParameter<std::string>("labelTrack")),
243  labelRecVtx_(iConfig.getParameter<std::string>("labelVertex")),
244  labelEB_(iConfig.getParameter<std::string>("labelEBRecHit")),
245  labelEE_(iConfig.getParameter<std::string>("labelEERecHit")),
246  labelHBHE_(iConfig.getParameter<std::string>("labelHBHERecHit")),
247  labelTower_(iConfig.getParameter<std::string>("labelCaloTower")),
248  l1TrigName_(iConfig.getUntrackedParameter<std::string>("l1TrigName", "L1_SingleJet60")),
249  oldID_(iConfig.getUntrackedParameter<std::vector<int> >("oldID")),
250  newDepth_(iConfig.getUntrackedParameter<std::vector<int> >("newDepth")),
251  hep17_(iConfig.getUntrackedParameter<bool>("hep17")),
252  usePFThresh_(iConfig.getParameter<bool>("usePFThreshold")),
253  labelBS_(iConfig.getParameter<std::string>("labelBeamSpot")),
254  modnam_(iConfig.getUntrackedParameter<std::string>("moduleName", "")),
255  prdnam_(iConfig.getUntrackedParameter<std::string>("producerName", "")),
256  algTag_(iConfig.getParameter<edm::InputTag>("algInputTag")),
257  extTag_(iConfig.getParameter<edm::InputTag>("extInputTag")),
258  tok_trigEvt_(consumes<trigger::TriggerEvent>(triggerEvent_)),
259  tok_trigRes_(consumes<edm::TriggerResults>(theTriggerResultsLabel_)),
260  tok_bs_(consumes<reco::BeamSpot>(labelBS_)),
261  tok_recVtx_((modnam_.empty()) ? consumes<reco::VertexCollection>(labelRecVtx_)
283  nRun_(0),
284  nLow_(0),
285  nHigh_(0),
286  hdc_(nullptr) {
287  usesResource(TFileService::kSharedResource);
288 
289  //now do whatever initialization is needed
290  const double isolationRadius(28.9), innerR(10.0), outerR(30.0);
291  a_charIsoR_ = a_coneR_ + isolationRadius;
292  a_coneR1_ = a_coneR_ + innerR;
293  a_coneR2_ = a_coneR_ + outerR;
294  // Different isolation cuts are described in DN-2016/029
295  // Tight cut uses 2 GeV; Loose cut uses 10 GeV
296  // Eta dependent cut uses (maxRestrictionP_ * exp(|ieta|*log(2.5)/18))
297  // with the factor for exponential slopeRestrictionP_ = log(2.5)/18
298  // maxRestrictionP_ = 8 GeV as came from a study
299 
300  for (unsigned int k = 0; k < oldID_.size(); ++k) {
301  oldDet_.emplace_back((oldID_[k] / 10000) % 10);
302  oldEta_.emplace_back((oldID_[k] / 100) % 100);
303  oldDepth_.emplace_back(oldID_[k] % 100);
304  }
305 
306  l1GtUtils_ =
308 
309  if (modnam_.empty()) {
310  edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
311  << labelBS_ << " " << labelRecVtx_ << " " << labelGenTrack_ << " "
312  << edm::InputTag("ecalRecHit", labelEB_) << " "
313  << edm::InputTag("ecalRecHit", labelEE_) << " " << labelHBHE_ << " "
314  << labelTower_;
315  } else {
316  edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
317  << labelBS_ << " " << edm::InputTag(modnam_, labelRecVtx_, prdnam_) << " "
318  << labelGenTrack_ << " " << edm::InputTag(modnam_, labelEB_, prdnam_) << " "
321  }
322 
323  edm::LogVerbatim("HcalIsoTrack") << "Parameters read from config file \n"
324  << "\t minPt " << ptMin_ << "\t etaMax " << etaMax_ << "\t a_coneR " << a_coneR_
325  << ":" << a_coneR1_ << ":" << a_coneR2_ << "\t a_charIsoR " << a_charIsoR_
326  << "\t a_mipR " << a_mipR_ << "\t a_mipR2 " << a_mipR2_ << "\t a_mipR3 " << a_mipR3_
327  << "\t a_mipR4 " << a_mipR4_ << "\t a_mipR5 " << a_mipR5_ << "\n pTrackMin_ "
328  << pTrackMin_ << "\t eEcalMax_ " << eEcalMax_ << "\t maxRestrictionP_ "
329  << maxRestrictionP_ << "\t slopeRestrictionP_ " << slopeRestrictionP_
330  << "\t eIsolateStrong_ " << eIsolate1_ << "\t eIsolateSoft_ " << eIsolate2_
331  << "\t hcalScale_ " << hcalScale_ << "\n\t momentumLow_ " << pTrackLow_
332  << "\t prescaleLow_ " << prescaleLow_ << "\t momentumHigh_ " << pTrackHigh_
333  << "\t prescaleHigh_ " << prescaleHigh_ << "\n\t useRaw_ " << useRaw_
334  << "\t ignoreTrigger_ " << ignoreTrigger_ << "\n\t useL1Trigegr_ " << useL1Trigger_
335  << "\t dataType_ " << dataType_ << "\t mode_ " << mode_
336  << "\t unCorrect_ " << unCorrect_ << "\t collapseDepth_ " << collapseDepth_
337  << "\t fillTreeRange " << fillTreeRange_ << "\t L1TrigName_ " << l1TrigName_
338  << "\nThreshold flag used " << usePFThresh_ << " value for EB " << hitEthrEB_
339  << " EE " << hitEthrEE0_ << ":" << hitEthrEE1_ << ":" << hitEthrEE2_ << ":"
340  << hitEthrEE3_ << ":" << hitEthrEELo_ << ":" << hitEthrEEHi_;
341  edm::LogVerbatim("HcalIsoTrack") << "Process " << processName_ << " L1Filter:" << l1Filter_
342  << " L2Filter:" << l2Filter_ << " L3Filter:" << l3Filter_;
343  for (unsigned int k = 0; k < trigNames_.size(); ++k) {
344  edm::LogVerbatim("HcalIsoTrack") << "Trigger[" << k << "] " << trigNames_[k];
345  }
346  edm::LogVerbatim("HcalIsoTrack") << oldID_.size() << " DetIDs to be corrected with HEP17 flag:" << hep17_;
347  for (unsigned int k = 0; k < oldID_.size(); ++k)
348  edm::LogVerbatim("HcalIsoTrack") << "[" << k << "] Det " << oldDet_[k] << " EtaAbs " << oldEta_[k] << " Depth "
349  << oldDepth_[k] << ":" << newDepth_[k];
350 
351  for (int i = 0; i < 10; i++)
352  phibins_.push_back(-M_PI + 0.1 * (2 * i + 1) * M_PI);
353  for (int i = 0; i < 8; ++i)
354  etabins_.push_back(-2.1 + 0.6 * i);
355  etadist_ = etabins_[1] - etabins_[0];
356  phidist_ = phibins_[1] - phibins_[0];
357  etahalfdist_ = 0.5 * etadist_;
358  phihalfdist_ = 0.5 * phidist_;
359  edm::LogVerbatim("HcalIsoTrack") << "EtaDist " << etadist_ << " " << etahalfdist_ << " PhiDist " << phidist_ << " "
360  << phihalfdist_;
361  unsigned int k1(0), k2(0);
362  for (auto phi : phibins_) {
363  edm::LogVerbatim("HcalIsoTrack") << "phibin_[" << k1 << "] " << phi;
364  ++k1;
365  }
366  for (auto eta : etabins_) {
367  edm::LogVerbatim("HcalIsoTrack") << "etabin_[" << k2 << "] " << eta;
368  ++k2;
369  }
370 }
static const std::string kSharedResource
Definition: TFileService.h:76
Log< level::Info, true > LogVerbatim
std::vector< GenParticle > GenParticleCollection
collection of GenParticles
const std::vector< int > newDepth_
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
HepPDT::ParticleDataTable ParticleDataTable
const std::string labelRecVtx_
const std::vector< std::string > trigNames_
l1t::L1TGlobalUtil * l1GtUtils_
const edm::InputTag algTag_
std::vector< double > etabins_
std::vector< double > phibins_
const edm::InputTag triggerEvent_
std::vector< Vertex > VertexCollection
Definition: Vertex.h:31
edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd > tok_sevlv_
edm::ESGetToken< CaloTopology, CaloTopologyRecord > tok_caloTopology_
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
const std::vector< int > oldID_
std::vector< int > oldDepth_
edm::EDGetTokenT< BXVector< GlobalAlgBlk > > tok_alg_
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
const HcalDDDRecConstants * hdc_
Transition
Definition: Transition.h:12
edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tok_htopo_
const edm::InputTag extTag_
edm::EDGetTokenT< reco::GenParticleCollection > tok_parts_
edm::EDGetTokenT< GenEventInfoProduct > tok_ew_
#define M_PI
edm::EDGetTokenT< CaloTowerCollection > tok_cala_
edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > tok_ecalChStatus_
edm::EDGetTokenT< edm::TriggerResults > tok_trigRes_
edm::ESGetToken< EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd > tok_ecalPFRecHitThresholds_
edm::EDGetTokenT< trigger::TriggerEvent > tok_trigEvt_
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
const std::string processName_
edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord > tok_ddrec_
fixed size matrix
HLT enums.
edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_resp_
edm::ESGetToken< HepPDT::ParticleDataTable, PDTRecord > tok_pdt_
const std::string labelGenTrack_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord > tok_geom_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > tok_bFieldH_
const edm::InputTag theTriggerResultsLabel_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_

◆ ~HcalIsoTrkSimAnalyzer()

HcalIsoTrkSimAnalyzer::~HcalIsoTrkSimAnalyzer ( )
overridedefault

Member Function Documentation

◆ analyze()

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

Implements edm::one::EDAnalyzerBase.

Definition at line 372 of file HcalIsoTrkSimAnalyzer.cc.

References Calorimetry_cff::bField, dataType_, l1t::L1TGlobalUtil::decisionsFinal(), fileCollector::done, l1ctLayer1_cff::dr, dR(), trigger::TriggerObject::energy(), trigger::TriggerObject::eta(), etaMax_, eThresholds_, fillTree(), spr::find(), genParticles2HepMC_cfi::genEventInfo, AJJGenJetFilter_cfi::genParticles, edm::EventSetup::getData(), photonIsolationHIProducer_cfi::hbhe, ValidationMatrix::hlt, hltConfig_, mps_fire::i, trigger::TriggerObject::id(), iEvent, ignoreTrigger_, edm::HandleBase::isValid(), dqmdumpme::k, l1Filter_, l1GtUtils_, l1TrigName_, l2Filter_, l3Filter_, label, edm::InputTag::label(), labelEB_, labelEE_, labelHBHE_, trigger::TriggerObject::mass(), HLTConfigProvider::moduleLabels(), LogMessageMonitor_cff::modules, names, convertSQLiteXML::ok, trigger::TriggerObject::phi(), reco::BeamSpot::position(), edm::Handle< T >::product(), spr::propagateCALO(), trigger::TriggerObject::pt(), trigger::TriggerObject::px(), trigger::TriggerObject::py(), trigger::TriggerObject::pz(), l1t::L1TGlobalUtil::retrieveL1(), HcalCondObjectContainerBase::setTopo(), findQualityFiles::size, AlCaHLTBitMon_QueryRunRegistry::string, t_allvertex, t_DataType, t_Event, t_EventNo, t_EventWeight, t_goodPV, t_hltbits, t_ietaAll, t_ietaGood, t_L1Bit, t_l1eta, t_l1phi, t_l1pt, t_l3eta, t_l3phi, t_l3pt, t_nVtx, t_Run, t_RunNo, t_Tracks, t_TracksLoose, t_TracksProp, t_TracksSaved, t_TracksTight, t_trackType, t_trgbits, t_TrigPass, t_TrigPassSel, tok_alg_, tok_bFieldH_, tok_bs_, tok_cala_, tok_caloTopology_, tok_EB_, tok_ecalChStatus_, tok_ecalPFRecHitThresholds_, tok_EE_, tok_ew_, tok_geom_, tok_hbhe_, tok_htopo_, tok_parts_, tok_pdt_, tok_recVtx_, tok_resp_, tok_sevlv_, tok_trigEvt_, tok_trigRes_, tree2, PDWG_DiPhoton_SD_cff::triggerEvent, triggerEvent_, HLTConfigProvider::triggerIndex(), L1TEGammaOffline_cfi::triggerNames, triggerResults, trigNames_, useL1Trigger_, x, y, and z.

372  {
373  t_Run = iEvent.id().run();
374  t_Event = iEvent.id().event();
376 #ifdef EDM_ML_DEBUG
377  edm::LogVerbatim("HcalIsoTrack") << "Run " << t_Run << " Event " << t_Event << " type " << t_DataType
378  << " Luminosity " << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing();
379 #endif
380  //Get magnetic field and ECAL channel status
381  const MagneticField* bField = &iSetup.getData(tok_bFieldH_);
382  const EcalChannelStatus* theEcalChStatus = &iSetup.getData(tok_ecalChStatus_);
383  const EcalSeverityLevelAlgo* theEcalSevlv = &iSetup.getData(tok_sevlv_);
384 
385  // get calogeometry and calotopology
386  const CaloGeometry* geo = &iSetup.getData(tok_geom_);
387  const CaloTopology* caloTopology = &iSetup.getData(tok_caloTopology_);
388  const HcalTopology* theHBHETopology = &iSetup.getData(tok_htopo_);
389 
390  // get response correction
391  const HcalRespCorrs* resp = &iSetup.getData(tok_resp_);
392  HcalRespCorrs* respCorrs = new HcalRespCorrs(*resp);
393  respCorrs->setTopo(theHBHETopology);
394 
395  // get ECAL thresholds
396  eThresholds_ = &iSetup.getData(tok_ecalPFRecHitThresholds_);
397 
398  // get particle data table
399  const HepPDT::ParticleDataTable* pdt = &iSetup.getData(tok_pdt_);
400 
401  //=== genParticle information
402  bool okC(true);
404  if (!genParticles.isValid()) {
405  edm::LogWarning("HcalIsoTrack") << "Cannot access the genParticles collection";
406  okC = false;
407  }
408 
409  //event weight for FLAT sample
410  t_EventWeight = 1.0;
412  if (genEventInfo.isValid())
413  t_EventWeight = genEventInfo->weight();
414 
415  //Define the best vertex and the beamspot
416  const edm::Handle<reco::VertexCollection> recVtxs = iEvent.getHandle(tok_recVtx_);
417  const edm::Handle<reco::BeamSpot> beamSpotH = iEvent.getHandle(tok_bs_);
418  math::XYZPoint leadPV(0, 0, 0);
419  t_goodPV = t_nVtx = 0;
420  if (recVtxs.isValid() && !(recVtxs->empty())) {
421  t_nVtx = recVtxs->size();
422  for (unsigned int k = 0; k < recVtxs->size(); ++k) {
423  if (!((*recVtxs)[k].isFake()) && ((*recVtxs)[k].ndof() > 4)) {
424  if (t_goodPV == 0)
425  leadPV = math::XYZPoint((*recVtxs)[k].x(), (*recVtxs)[k].y(), (*recVtxs)[k].z());
426  t_goodPV++;
427  }
428  }
429  }
430  if (t_goodPV == 0 && beamSpotH.isValid()) {
431  leadPV = beamSpotH->position();
432  }
434 #ifdef EDM_ML_DEBUG
435  edm::LogVerbatim("HcalIsoTrack") << "Primary Vertex " << leadPV << " out of " << t_goodPV << " vertex";
436  if (beamSpotH.isValid()) {
437  edm::LogVerbatim("HcalIsoTrack") << " Beam Spot " << beamSpotH->position();
438  }
439 #endif
440  // RecHits
441  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle = iEvent.getHandle(tok_EB_);
442  if (!barrelRecHitsHandle.isValid()) {
443  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEB_;
444  okC = false;
445  }
446  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle = iEvent.getHandle(tok_EE_);
447  if (!endcapRecHitsHandle.isValid()) {
448  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEE_;
449  okC = false;
450  }
452  if (!hbhe.isValid()) {
453  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelHBHE_;
454  okC = false;
455  }
456  const edm::Handle<CaloTowerCollection> caloTower = iEvent.getHandle(tok_cala_);
457 
458  //Propagate tracks to calorimeter surface)
459  std::vector<spr::propagatedGenParticleID> trackIDs =
460  spr::propagateCALO(genParticles, pdt, geo, bField, etaMax_, false);
461  std::vector<math::XYZTLorentzVector> vecL1, vecL3;
462  t_RunNo = iEvent.id().run();
463  t_EventNo = iEvent.id().event();
464  t_Tracks = (genParticles.product())->size();
465  t_TracksProp = trackIDs.size();
466  t_ietaAll->clear();
467  t_ietaGood->clear();
468  t_trackType->clear();
469  t_trgbits->clear();
470  t_hltbits->clear();
471 #ifdef EDM_ML_DEBUG
472  edm::LogVerbatim("HcalIsoTrack") << "# of propagated tracks " << t_TracksProp << " out of " << t_Tracks
473  << " with Trigger " << ignoreTrigger_;
474 #endif
475 
476  //Trigger
477  t_trgbits->assign(trigNames_.size(), false);
478  t_hltbits->assign(trigNames_.size(), false);
480  t_L1Bit = true;
481  t_TrigPass = false;
482 
483  if (!ignoreTrigger_) {
484  //L1
486  const std::vector<std::pair<std::string, bool> >& finalDecisions = l1GtUtils_->decisionsFinal();
487  for (const auto& decision : finalDecisions) {
488  if (decision.first.find(l1TrigName_) != std::string::npos) {
489  t_L1Bit = decision.second;
490  break;
491  }
492  }
493 #ifdef EDM_ML_DEBUG
494  edm::LogVerbatim("HcalIsoTrack") << "Trigger Information for " << l1TrigName_ << " is " << t_L1Bit
495  << " from a list of " << finalDecisions.size() << " decisions";
496 #endif
497 
498  //HLT
500  if (triggerResults.isValid()) {
501  const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
502  const std::vector<std::string>& names = triggerNames.triggerNames();
503  if (!trigNames_.empty()) {
504  for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
505  int hlt = triggerResults->accept(iHLT);
506  for (unsigned int i = 0; i < trigNames_.size(); ++i) {
507  if (names[iHLT].find(trigNames_[i]) != std::string::npos) {
508  t_trgbits->at(i) = (hlt > 0);
509  t_hltbits->at(i) = (hlt > 0);
510  if (hlt > 0)
511  t_TrigPass = true;
512 #ifdef EDM_ML_DEBUG
513  edm::LogVerbatim("HcalIsoTrack")
514  << "This trigger " << names[iHLT] << " Flag " << hlt << ":" << t_trgbits->at(i);
515 #endif
516  }
517  }
518  }
519  }
520  }
521 #ifdef EDM_ML_DEBUG
522  edm::LogVerbatim("HcalIsoTrack") << "HLT Information shows " << t_TrigPass << ":" << trigNames_.empty() << ":"
523  << okC;
524 #endif
525  }
526 
527  std::array<int, 3> ntksave{{0, 0, 0}};
528  if (ignoreTrigger_ || useL1Trigger_) {
529  t_l1pt = t_l1eta = t_l1phi = 0;
530  t_l3pt = t_l3eta = t_l3phi = 0;
531  if (ignoreTrigger_ || t_L1Bit)
532  ntksave = fillTree(vecL1,
533  vecL3,
534  leadPV,
535  trackIDs,
536  geo,
537  caloTopology,
538  theHBHETopology,
539  theEcalChStatus,
540  theEcalSevlv,
541  barrelRecHitsHandle,
542  endcapRecHitsHandle,
543  hbhe,
544  caloTower,
545  genParticles,
546  respCorrs);
547  t_TracksSaved = ntksave[0];
548  t_TracksLoose = ntksave[1];
549  t_TracksTight = ntksave[2];
550  } else {
552  edm::Handle<trigger::TriggerEvent> triggerEventHandle = iEvent.getHandle(tok_trigEvt_);
553  if (!triggerEventHandle.isValid()) {
554  edm::LogWarning("HcalIsoTrack") << "Error! Can't get the product " << triggerEvent_.label();
555  } else if (okC) {
556  triggerEvent = *(triggerEventHandle.product());
558  const trigger::TriggerObjectCollection& TOC(triggerEvent.getObjects());
559  bool done(false);
560  if (triggerResults.isValid()) {
561  std::vector<std::string> modules;
562  const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
563  const std::vector<std::string>& names = triggerNames.triggerNames();
564  for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
565  bool ok = (t_TrigPass) || (trigNames_.empty());
566  if (ok) {
567  unsigned int triggerindx = hltConfig_.triggerIndex(names[iHLT]);
568  const std::vector<std::string>& moduleLabels(hltConfig_.moduleLabels(triggerindx));
569  std::vector<math::XYZTLorentzVector> vecL2;
570  vecL1.clear();
571  vecL3.clear();
572  //loop over all trigger filters in event (i.e. filters passed)
573  for (unsigned int ifilter = 0; ifilter < triggerEvent.sizeFilters(); ++ifilter) {
574  std::vector<int> Keys;
575  std::string label = triggerEvent.filterTag(ifilter).label();
576  //loop over keys to objects passing this filter
577  for (unsigned int imodule = 0; imodule < moduleLabels.size(); imodule++) {
578  if (label.find(moduleLabels[imodule]) != std::string::npos) {
579 #ifdef EDM_ML_DEBUG
580  edm::LogVerbatim("HcalIsoTrack") << "FilterName " << label;
581 #endif
582  for (unsigned int ifiltrKey = 0; ifiltrKey < triggerEvent.filterKeys(ifilter).size(); ++ifiltrKey) {
583  Keys.push_back(triggerEvent.filterKeys(ifilter)[ifiltrKey]);
584  const trigger::TriggerObject& TO(TOC[Keys[ifiltrKey]]);
585  math::XYZTLorentzVector v4(TO.px(), TO.py(), TO.pz(), TO.energy());
586  if (label.find(l2Filter_) != std::string::npos) {
587  vecL2.push_back(v4);
588  } else if (label.find(l3Filter_) != std::string::npos) {
589  vecL3.push_back(v4);
590  } else if ((label.find(l1Filter_) != std::string::npos) || (l1Filter_.empty())) {
591  vecL1.push_back(v4);
592  }
593 #ifdef EDM_ML_DEBUG
594  edm::LogVerbatim("HcalIsoTrack")
595  << "key " << ifiltrKey << " : pt " << TO.pt() << " eta " << TO.eta() << " phi " << TO.phi()
596  << " mass " << TO.mass() << " Id " << TO.id();
597 #endif
598  }
599 #ifdef EDM_ML_DEBUG
600  edm::LogVerbatim("HcalIsoTrack")
601  << "sizes " << vecL1.size() << ":" << vecL2.size() << ":" << vecL3.size();
602 #endif
603  }
604  }
605  }
607  math::XYZTLorentzVector mindRvec1;
608  double mindR1(999);
609  for (unsigned int i = 0; i < vecL2.size(); i++) {
610  double dr = dR(vecL1[0], vecL2[i]);
611 #ifdef EDM_ML_DEBUG
612  edm::LogVerbatim("HcalIsoTrack") << "lvl2[" << i << "] dR " << dr;
613 #endif
614  if (dr < mindR1) {
615  mindR1 = dr;
616  mindRvec1 = vecL2[i];
617  }
618  }
619 #ifdef EDM_ML_DEBUG
620  edm::LogVerbatim("HcalIsoTrack") << "L2 object closest to L1 " << mindRvec1 << " at Dr " << mindR1;
621 #endif
622 
623  if (!vecL1.empty()) {
624  t_l1pt = vecL1[0].pt();
625  t_l1eta = vecL1[0].eta();
626  t_l1phi = vecL1[0].phi();
627  } else {
628  t_l1pt = t_l1eta = t_l1phi = 0;
629  }
630  if (!vecL3.empty()) {
631  t_l3pt = vecL3[0].pt();
632  t_l3eta = vecL3[0].eta();
633  t_l3phi = vecL3[0].phi();
634  } else {
635  t_l3pt = t_l3eta = t_l3phi = 0;
636  }
637  // Now fill in the tree for each selected track
638  if (!done) {
639  ntksave = fillTree(vecL1,
640  vecL3,
641  leadPV,
642  trackIDs,
643  geo,
644  caloTopology,
645  theHBHETopology,
646  theEcalChStatus,
647  theEcalSevlv,
648  barrelRecHitsHandle,
649  endcapRecHitsHandle,
650  hbhe,
651  caloTower,
652  genParticles,
653  respCorrs);
654  t_TracksSaved += ntksave[0];
655  t_TracksLoose += ntksave[1];
656  t_TracksTight += ntksave[2];
657  done = true;
658  }
659  }
660  }
661  }
662  }
663  }
664 #ifdef EDM_ML_DEBUG
665  edm::LogVerbatim("HcalIsoTrack") << "Final results on selected tracks " << t_TracksSaved << ":" << t_TracksLoose
666  << ":" << t_TracksTight;
667 #endif
669  tree2->Fill();
670 }
size
Write out results.
Log< level::Info, true > LogVerbatim
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
std::vector< spr::propagatedTrackID > propagateCALO(edm::Handle< reco::TrackCollection > &trkCollection, const CaloGeometry *geo, const MagneticField *bField, const std::string &theTrackQuality, bool debug=false)
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
The single EDProduct to be saved for each event (AOD case)
Definition: TriggerEvent.h:26
const Point & position() const
position
Definition: BeamSpot.h:59
HepPDT::ParticleDataTable ParticleDataTable
double dR(math::XYZTLorentzVector &, math::XYZTLorentzVector &)
T const * product() const
Definition: Handle.h:70
const std::vector< std::string > trigNames_
l1t::L1TGlobalUtil * l1GtUtils_
std::vector< int > * t_ietaGood
std::string const & label() const
Definition: InputTag.h:36
const edm::InputTag triggerEvent_
std::vector< bool > * t_hltbits
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
const std::string names[nVars_]
edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd > tok_sevlv_
edm::ESGetToken< CaloTopology, CaloTopologyRecord > tok_caloTopology_
Single trigger physics object (e.g., an isolated muon)
Definition: TriggerObject.h:21
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
char const * label
int iEvent
Definition: GenABIO.cc:224
const EcalPFRecHitThresholds * eThresholds_
edm::EDGetTokenT< BXVector< GlobalAlgBlk > > tok_alg_
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
std::vector< int > * t_ietaAll
std::vector< int > * t_trackType
edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tok_htopo_
unsigned int triggerIndex(const std::string &triggerName) const
slot position of trigger path in trigger table (0 to size-1)
std::array< int, 3 > fillTree(std::vector< math::XYZTLorentzVector > &vecL1, std::vector< math::XYZTLorentzVector > &vecL3, math::XYZPoint &leadPV, std::vector< spr::propagatedGenParticleID > &trackIDs, const CaloGeometry *geo, const CaloTopology *topo, const HcalTopology *theHBHETopology, const EcalChannelStatus *theEcalChStatus, const EcalSeverityLevelAlgo *theEcalSevlv, edm::Handle< EcalRecHitCollection > &barrelRecHitsHandle, edm::Handle< EcalRecHitCollection > &endcapRecHitsHandle, edm::Handle< HBHERecHitCollection > &hbhe, const edm::Handle< CaloTowerCollection > &towerHandle, edm::Handle< reco::GenParticleCollection > &genParticles, const HcalRespCorrs *respCorrs)
static std::string const triggerResults
Definition: EdmProvDump.cc:47
edm::EDGetTokenT< reco::GenParticleCollection > tok_parts_
edm::EDGetTokenT< GenEventInfoProduct > tok_ew_
edm::EDGetTokenT< CaloTowerCollection > tok_cala_
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:75
edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > tok_ecalChStatus_
edm::EDGetTokenT< edm::TriggerResults > tok_trigRes_
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
std::vector< size_type > Keys
edm::ESGetToken< EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd > tok_ecalPFRecHitThresholds_
edm::EDGetTokenT< trigger::TriggerEvent > tok_trigEvt_
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
bool isValid() const
Definition: HandleBase.h:70
std::vector< bool > * t_trgbits
edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_resp_
const std::vector< std::pair< std::string, bool > > & decisionsFinal()
edm::ESGetToken< HepPDT::ParticleDataTable, PDTRecord > tok_pdt_
Log< level::Warning, false > LogWarning
edm::ESGetToken< CaloGeometry, CaloGeometryRecord > tok_geom_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > tok_bFieldH_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
void setTopo(const HcalTopology *topo)
void retrieveL1(const edm::Event &iEvent, const edm::EventSetup &evSetup)
initialize the class (mainly reserve)

◆ beginJob()

void HcalIsoTrkSimAnalyzer::beginJob ( )
overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 672 of file HcalIsoTrkSimAnalyzer.cc.

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

672  {
674  tree = fs->make<TTree>("CalibTree", "CalibTree");
675 
676  tree->Branch("t_Run", &t_Run, "t_Run/I");
677  tree->Branch("t_Event", &t_Event, "t_Event/I");
678  tree->Branch("t_DataType", &t_DataType, "t_DataType/I");
679  tree->Branch("t_ieta", &t_ieta, "t_ieta/I");
680  tree->Branch("t_iphi", &t_iphi, "t_iphi/I");
681  tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
682  tree->Branch("t_nVtx", &t_nVtx, "t_nVtx/I");
683  tree->Branch("t_nTrk", &t_nTrk, "t_nTrk/I");
684  tree->Branch("t_goodPV", &t_goodPV, "t_goodPV/I");
685  tree->Branch("t_l1pt", &t_l1pt, "t_l1pt/D");
686  tree->Branch("t_l1eta", &t_l1eta, "t_l1eta/D");
687  tree->Branch("t_l1phi", &t_l1phi, "t_l1phi/D");
688  tree->Branch("t_l3pt", &t_l3pt, "t_l3pt/D");
689  tree->Branch("t_l3eta", &t_l3eta, "t_l3eta/D");
690  tree->Branch("t_l3phi", &t_l3phi, "t_l3phi/D");
691  tree->Branch("t_p", &t_p, "t_p/D");
692  tree->Branch("t_pt", &t_pt, "t_pt/D");
693  tree->Branch("t_phi", &t_phi, "t_phi/D");
694  tree->Branch("t_mindR1", &t_mindR1, "t_mindR1/D");
695  tree->Branch("t_mindR2", &t_mindR2, "t_mindR2/D");
696  tree->Branch("t_eMipDR", &t_eMipDR, "t_eMipDR/D");
697  tree->Branch("t_eMipDR2", &t_eMipDR2, "t_eMipDR2/D");
698  tree->Branch("t_eMipDR3", &t_eMipDR3, "t_eMipDR3/D");
699  tree->Branch("t_eMipDR4", &t_eMipDR4, "t_eMipDR4/D");
700  tree->Branch("t_eMipDR5", &t_eMipDR5, "t_eMipDR5/D");
701  tree->Branch("t_eHcal", &t_eHcal, "t_eHcal/D");
702  tree->Branch("t_eHcal10", &t_eHcal10, "t_eHcal10/D");
703  tree->Branch("t_eHcal30", &t_eHcal30, "t_eHcal30/D");
704  tree->Branch("t_hmaxNearP", &t_hmaxNearP, "t_hmaxNearP/D");
705  tree->Branch("t_emaxNearP", &t_emaxNearP, "t_emaxNearP/D");
706  tree->Branch("t_eAnnular", &t_eAnnular, "t_eAnnular/D");
707  tree->Branch("t_hAnnular", &t_hAnnular, "t_hAnnular/D");
708  tree->Branch("t_rhoh", &t_rhoh, "t_rhoh/D");
709  tree->Branch("t_selectTk", &t_selectTk, "t_selectTk/O");
710  tree->Branch("t_qltyFlag", &t_qltyFlag, "t_qltyFlag/O");
711  tree->Branch("t_qltyMissFlag", &t_qltyMissFlag, "t_qltyMissFlag/O");
712  tree->Branch("t_qltyPVFlag", &t_qltyPVFlag, "t_qltyPVFlag/O");
713  tree->Branch("t_gentrackP", &t_gentrackP, "t_gentrackP/D");
714 
715  t_DetIds = new std::vector<unsigned int>();
716  t_DetIds1 = new std::vector<unsigned int>();
717  t_DetIds3 = new std::vector<unsigned int>();
718  t_HitEnergies = new std::vector<double>();
719  t_HitEnergies1 = new std::vector<double>();
720  t_HitEnergies3 = new std::vector<double>();
721  t_trgbits = new std::vector<bool>();
722  tree->Branch("t_DetIds", "std::vector<unsigned int>", &t_DetIds);
723  tree->Branch("t_HitEnergies", "std::vector<double>", &t_HitEnergies);
724  tree->Branch("t_trgbits", "std::vector<bool>", &t_trgbits);
725  tree->Branch("t_DetIds1", "std::vector<unsigned int>", &t_DetIds1);
726  tree->Branch("t_DetIds3", "std::vector<unsigned int>", &t_DetIds3);
727  tree->Branch("t_HitEnergies1", "std::vector<double>", &t_HitEnergies1);
728  tree->Branch("t_HitEnergies3", "std::vector<double>", &t_HitEnergies3);
729 
730  tree2 = fs->make<TTree>("EventInfo", "Event Information");
731 
732  tree2->Branch("t_RunNo", &t_RunNo, "t_RunNo/i");
733  tree2->Branch("t_EventNo", &t_EventNo, "t_EventNo/i");
734  tree2->Branch("t_Tracks", &t_Tracks, "t_Tracks/I");
735  tree2->Branch("t_TracksProp", &t_TracksProp, "t_TracksProp/I");
736  tree2->Branch("t_TracksSaved", &t_TracksSaved, "t_TracksSaved/I");
737  tree2->Branch("t_TracksLoose", &t_TracksLoose, "t_TracksLoose/I");
738  tree2->Branch("t_TracksTight", &t_TracksTight, "t_TracksTight/I");
739  tree2->Branch("t_TrigPass", &t_TrigPass, "t_TrigPass/O");
740  tree2->Branch("t_TrigPassSel", &t_TrigPassSel, "t_TrigPassSel/O");
741  tree2->Branch("t_L1Bit", &t_L1Bit, "t_L1Bit/O");
742  tree2->Branch("t_allvertex", &t_allvertex, "t_allvertex/I");
743  t_hltbits = new std::vector<bool>();
744  t_ietaAll = new std::vector<int>();
745  t_ietaGood = new std::vector<int>();
746  t_trackType = new std::vector<int>();
747  tree2->Branch("t_ietaAll", "std::vector<int>", &t_ietaAll);
748  tree2->Branch("t_ietaGood", "std::vector<int>", &t_ietaGood);
749  tree2->Branch("t_trackType", "std::vector<int>", &t_trackType);
750  tree2->Branch("t_hltbits", "std::vector<bool>", &t_hltbits);
751 }
std::vector< double > * t_HitEnergies3
std::vector< double > * t_HitEnergies
std::vector< int > * t_ietaGood
std::vector< bool > * t_hltbits
std::vector< unsigned int > * t_DetIds1
std::vector< int > * t_ietaAll
std::vector< int > * t_trackType
std::vector< unsigned int > * t_DetIds3
std::vector< unsigned int > * t_DetIds
std::vector< double > * t_HitEnergies1
std::vector< bool > * t_trgbits
Definition: tree.py:1

◆ beginRun()

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

Definition at line 754 of file HcalIsoTrkSimAnalyzer.cc.

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

754  {
755  hdc_ = &iSetup.getData(tok_ddrec_);
756 
757  bool changed_(true);
758  bool flag = hltConfig_.init(iRun, iSetup, processName_, changed_);
759  edm::LogVerbatim("HcalIsoTrack") << "Run[" << nRun_ << "] " << iRun.run() << " process " << processName_
760  << " init flag " << flag << " change flag " << changed_;
761  // check if trigger names in (new) config
762  if (changed_) {
763 #ifdef EDM_ML_DEBUG
764  edm::LogVerbatim("HcalIsoTrack") << "New trigger menu found !!!";
765 #endif
766  const unsigned int n(hltConfig_.size());
767  for (unsigned itrig = 0; itrig < trigNames_.size(); itrig++) {
768  unsigned int triggerindx = hltConfig_.triggerIndex(trigNames_[itrig]);
769  if (triggerindx >= n) {
770  edm::LogWarning("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " does not exist in "
771  << "the current menu";
772 #ifdef EDM_ML_DEBUG
773  } else {
774  edm::LogVerbatim("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " exists";
775 #endif
776  }
777  }
778  }
779 }
Log< level::Info, true > LogVerbatim
const std::vector< std::string > trigNames_
unsigned int size() const
number of trigger paths in trigger table
const HcalDDDRecConstants * hdc_
unsigned int triggerIndex(const std::string &triggerName) const
slot position of trigger path in trigger table (0 to size-1)
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d&#39;tor
const std::string processName_
edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord > tok_ddrec_
Log< level::Warning, false > LogWarning

◆ dR()

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

Definition at line 1223 of file HcalIsoTrkSimAnalyzer.cc.

References reco::deltaR().

Referenced by analyze(), and fillTree().

1223  {
1224  return reco::deltaR(vec1.eta(), vec1.phi(), vec2.eta(), vec2.phi());
1225 }
std::vector< double > vec1
Definition: HCALResponse.h:15
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
std::vector< vec1 > vec2
Definition: HCALResponse.h:16

◆ endRun()

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

Definition at line 782 of file HcalIsoTrkSimAnalyzer.cc.

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

782  {
783  nRun_++;
784  edm::LogVerbatim("HcalIsoTrack") << "endRun[" << nRun_ << "] " << iRun.run();
785 }
Log< level::Info, true > LogVerbatim

◆ eThreshold()

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

Definition at line 1258 of file HcalIsoTrkSimAnalyzer.cc.

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

Referenced by fillTree().

1258  {
1259  double eThr(hitEthrEB_);
1260  if (usePFThresh_) {
1261  eThr = static_cast<double>((*eThresholds_)[id]);
1262  } else {
1263  const GlobalPoint& pos = geo->getPosition(id);
1264  double eta = std::abs(pos.eta());
1265  if (id.subdetId() != EcalBarrel) {
1266  eThr = (((eta * hitEthrEE3_ + hitEthrEE2_) * eta + hitEthrEE1_) * eta + hitEthrEE0_);
1267  if (eThr < hitEthrEELo_)
1268  eThr = hitEthrEELo_;
1269  else if (eThr > hitEthrEEHi_)
1270  eThr = hitEthrEEHi_;
1271  }
1272  }
1273  return eThr;
1274 }
GlobalPoint getPosition(const DetId &id) const
Get the position of a given detector id.
Definition: CaloGeometry.cc:50
Abs< T >::type abs(const T &t)
Definition: Abs.h:22

◆ fillDescriptions()

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

Definition at line 787 of file HcalIsoTrkSimAnalyzer.cc.

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

787  {
789  std::vector<std::string> trig = {"HLT_PFJet40",
790  "HLT_PFJet60",
791  "HLT_PFJet80",
792  "HLT_PFJet140",
793  "HLT_PFJet200",
794  "HLT_PFJet260",
795  "HLT_PFJet320",
796  "HLT_PFJet400",
797  "HLT_PFJet450",
798  "HLT_PFJet500"};
799  desc.add<std::vector<std::string> >("triggers", trig);
800  desc.add<std::string>("processName", "HLT");
801  desc.add<std::string>("l1Filter", "");
802  desc.add<std::string>("l2Filter", "L2Filter");
803  desc.add<std::string>("l3Filter", "Filter");
804  desc.add<double>("pTMin", 1.0), desc.add<double>("maxChargedHadronEta", 3.0);
805  // Minimum momentum of selected isolated track and signal zone
806  desc.add<double>("minimumTrackP", 10.0);
807  desc.add<double>("coneRadius", 34.98);
808  // signal zone in ECAL and MIP energy cutoff
809  desc.add<double>("coneRadiusMIP", 14.0);
810  desc.add<double>("coneRadiusMIP2", 18.0);
811  desc.add<double>("coneRadiusMIP3", 20.0);
812  desc.add<double>("coneRadiusMIP4", 22.0);
813  desc.add<double>("coneRadiusMIP5", 24.0);
814  desc.add<double>("maximumEcalEnergy", 2.0);
815  // following 4 parameters are for isolation cuts and described in the code
816  desc.add<double>("maxTrackP", 8.0);
817  desc.add<double>("slopeTrackP", 0.05090504066);
818  desc.add<double>("isolationEnergyTight", 2.0);
819  desc.add<double>("isolationEnergyLoose", 10.0);
820  // energy thershold for ECAL (from Egamma group)
821  desc.add<double>("EBHitEnergyThreshold", 0.08);
822  desc.add<double>("EEHitEnergyThreshold0", 0.30);
823  desc.add<double>("EEHitEnergyThreshold1", 0.00);
824  desc.add<double>("EEHitEnergyThreshold2", 0.00);
825  desc.add<double>("EEHitEnergyThreshold3", 0.00);
826  desc.add<double>("EEHitEnergyThresholdLow", 0.30);
827  desc.add<double>("EEHitEnergyThresholdHigh", 0.30);
828  // prescale factors
829  desc.add<double>("momentumLow", 40.0);
830  desc.add<double>("momentumHigh", 60.0);
831  desc.add<int>("prescaleLow", 1);
832  desc.add<int>("prescaleHigh", 1);
833  // various labels for collections used in the code
834  desc.add<edm::InputTag>("labelTriggerEvent", edm::InputTag("hltTriggerSummaryAOD", "", "HLT"));
835  desc.add<edm::InputTag>("labelTriggerResult", edm::InputTag("TriggerResults", "", "HLT"));
836  desc.add<std::string>("labelTrack", "generalTracks");
837  desc.add<std::string>("labelVertex", "offlinePrimaryVertices");
838  desc.add<std::string>("labelEBRecHit", "EcalRecHitsEB");
839  desc.add<std::string>("labelEERecHit", "EcalRecHitsEE");
840  desc.add<std::string>("labelHBHERecHit", "hbhereco");
841  desc.add<std::string>("labelBeamSpot", "offlineBeamSpot");
842  desc.add<std::string>("labelCaloTower", "towerMaker");
843  desc.add<edm::InputTag>("algInputTag", edm::InputTag("gtStage2Digis"));
844  desc.add<edm::InputTag>("extInputTag", edm::InputTag("gtStage2Digis"));
845  desc.addUntracked<std::string>("moduleName", "");
846  desc.addUntracked<std::string>("producerName", "");
847  // Various flags used for selecting tracks, choice of energy Method2/0
848  // Data type 0/1 for single jet trigger or others
849  desc.addUntracked<int>("useRaw", 0);
850  desc.addUntracked<bool>("ignoreTriggers", false);
851  desc.addUntracked<bool>("useL1Trigger", false);
852  desc.addUntracked<double>("hcalScale", 1.0);
853  desc.addUntracked<int>("dataType", 0);
854  desc.addUntracked<bool>("unCorrect", false);
855  desc.addUntracked<bool>("collapseDepth", false);
856  desc.addUntracked<bool>("fillTreeRange", false);
857  desc.addUntracked<std::string>("l1TrigName", "L1_SingleJet60");
858  desc.addUntracked<int>("outMode", 11);
859  std::vector<int> dummy;
860  desc.addUntracked<std::vector<int> >("oldID", dummy);
861  desc.addUntracked<std::vector<int> >("newDepth", dummy);
862  desc.addUntracked<bool>("hep17", false);
863  desc.add<bool>("usePFThreshold", true);
864  descriptions.add("hcalIsoTrkSimAnalyzer", desc);
865 }
void add(std::string const &label, ParameterSetDescription const &psetDescription)

◆ fillTree()

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

Definition at line 867 of file HcalIsoTrkSimAnalyzer.cc.

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

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

881  {
882  int nSave(0), nLoose(0), nTight(0);
883  //Loop over tracks
884  unsigned int nTracks(0);
885  t_nTrk = trackIDs.size();
886  t_rhoh = (tower.isValid()) ? rhoh(tower) : 0;
887  for (auto const& trkDetItr : trackIDs) {
888  const reco::GenParticle* pTrack = &(*(trkDetItr.trkItr));
889  math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
890 #ifdef EDM_ML_DEBUG
891  edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) :" << pTrack->pt() << "|"
892  << pTrack->eta() << "|" << pTrack->phi() << "|" << pTrack->p();
893 #endif
894  t_mindR2 = 999;
895  for (unsigned int k = 0; k < vecL3.size(); ++k) {
896  double dr = dR(vecL3[k], v4);
897  if (dr < t_mindR2) {
898  t_mindR2 = dr;
899  }
900  }
901  t_mindR1 = (!vecL1.empty()) ? dR(vecL1[0], v4) : 999;
902 #ifdef EDM_ML_DEBUG
903  edm::LogVerbatim("HcalIsoTrack") << "Closest L3 object at dr :" << t_mindR2 << " and from L1 " << t_mindR1;
904 #endif
905  t_ieta = t_iphi = 0;
906  if (trkDetItr.okHCAL) {
907  HcalDetId detId = (HcalDetId)(trkDetItr.detIdHCAL);
908  t_ieta = detId.ieta();
909  t_iphi = detId.iphi();
910  if (t_p > 40.0 && t_p <= 60.0)
911  t_ietaAll->emplace_back(t_ieta);
912  }
913  //Selection of good track
915  double eIsolation = maxRestrictionP_ * exp(slopeRestrictionP_ * std::abs((double)t_ieta));
916  if (eIsolation < eIsolate1_)
917  eIsolation = eIsolate1_;
918  if (eIsolation < eIsolate2_)
919  eIsolation = eIsolate2_;
920 #ifdef EDM_ML_DEBUG
921  edm::LogVerbatim("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " << t_qltyFlag << "|" << trkDetItr.okECAL << "|"
922  << trkDetItr.okHCAL << " eIsolation " << eIsolation;
923 #endif
924  t_qltyFlag = (t_selectTk && trkDetItr.okECAL && trkDetItr.okHCAL);
925  bool notMuon = notaMuon(pTrack);
926  if (t_qltyFlag && notMuon) {
927  int nNearTRKs(0);
929  std::vector<DetId> eIds;
930  std::vector<double> eHit;
932  barrelRecHitsHandle,
933  endcapRecHitsHandle,
934  trkDetItr.pointHCAL,
935  trkDetItr.pointECAL,
936  a_mipR_,
937  trkDetItr.directionECAL,
938  eIds,
939  eHit);
940  double eEcal(0);
941  for (unsigned int k = 0; k < eIds.size(); ++k) {
942  if (eHit[k] > eThreshold(eIds[k], geo))
943  eEcal += eHit[k];
944  }
945 #ifdef EDM_ML_DEBUG
946  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR << ":" << eEcal;
947 #endif
948  t_eMipDR = eEcal;
951  std::vector<DetId> eIds2;
952  std::vector<double> eHit2;
954  barrelRecHitsHandle,
955  endcapRecHitsHandle,
956  trkDetItr.pointHCAL,
957  trkDetItr.pointECAL,
958  a_mipR2_,
959  trkDetItr.directionECAL,
960  eIds2,
961  eHit2);
962  double eEcal2(0);
963  for (unsigned int k = 0; k < eIds2.size(); ++k) {
964  if (eHit2[k] > eThreshold(eIds2[k], geo))
965  eEcal2 += eHit2[k];
966  }
967 #ifdef EDM_ML_DEBUG
968  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR2 << ":" << eEcal2;
969 #endif
970  t_eMipDR2 = eEcal2;
973  std::vector<DetId> eIds3;
974  std::vector<double> eHit3;
976  barrelRecHitsHandle,
977  endcapRecHitsHandle,
978  trkDetItr.pointHCAL,
979  trkDetItr.pointECAL,
980  a_mipR3_,
981  trkDetItr.directionECAL,
982  eIds3,
983  eHit3);
984  double eEcal3(0);
985  for (unsigned int k = 0; k < eIds3.size(); ++k) {
986  if (eHit3[k] > eThreshold(eIds3[k], geo))
987  eEcal3 += eHit3[k];
988  }
989 #ifdef EDM_ML_DEBUG
990  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR3 << ":" << eEcal3;
991 #endif
992  t_eMipDR3 = eEcal3;
995  std::vector<DetId> eIds4;
996  std::vector<double> eHit4;
998  barrelRecHitsHandle,
999  endcapRecHitsHandle,
1000  trkDetItr.pointHCAL,
1001  trkDetItr.pointECAL,
1002  a_mipR4_,
1003  trkDetItr.directionECAL,
1004  eIds4,
1005  eHit4);
1006  double eEcal4(0);
1007  for (unsigned int k = 0; k < eIds4.size(); ++k) {
1008  if (eHit4[k] > eThreshold(eIds4[k], geo))
1009  eEcal4 += eHit4[k];
1010  }
1011 #ifdef EDM_ML_DEBUG
1012  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR4 << ":" << eEcal4;
1013 #endif
1014  t_eMipDR4 = eEcal4;
1017  std::vector<DetId> eIds5;
1018  std::vector<double> eHit5;
1019  t_eMipDR5 = spr::eCone_ecal(geo,
1020  barrelRecHitsHandle,
1021  endcapRecHitsHandle,
1022  trkDetItr.pointHCAL,
1023  trkDetItr.pointECAL,
1024  a_mipR5_,
1025  trkDetItr.directionECAL,
1026  eIds5,
1027  eHit5);
1028  double eEcal5(0);
1029  for (unsigned int k = 0; k < eIds5.size(); ++k) {
1030  if (eHit5[k] > eThreshold(eIds5[k], geo))
1031  eEcal5 += eHit5[k];
1032  }
1033 #ifdef EDM_ML_DEBUG
1034  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR5 << ":" << eEcal5;
1035 #endif
1036  t_eMipDR5 = eEcal5;
1038 
1039  t_emaxNearP = spr::chargeIsolationGenEcal(nTracks, trackIDs, geo, caloTopology, 15, 15);
1040  const DetId cellE(trkDetItr.detIdECAL);
1041  std::pair<double, bool> e11x11P = spr::eECALmatrix(cellE,
1042  barrelRecHitsHandle,
1043  endcapRecHitsHandle,
1044  *theEcalChStatus,
1045  geo,
1046  caloTopology,
1047  theEcalSevlv,
1048  5,
1049  5,
1050  -100.0,
1051  -100.0,
1052  -100.0,
1053  100.0);
1054  std::pair<double, bool> e15x15P = spr::eECALmatrix(cellE,
1055  barrelRecHitsHandle,
1056  endcapRecHitsHandle,
1057  *theEcalChStatus,
1058  geo,
1059  caloTopology,
1060  theEcalSevlv,
1061  7,
1062  7,
1063  -100.0,
1064  -100.0,
1065  -100.0,
1066  100.0);
1067  if (e11x11P.second && e15x15P.second) {
1068  t_eAnnular = (e15x15P.first - e11x11P.first);
1069  } else {
1070  t_eAnnular = -(e15x15P.first - e11x11P.first);
1071  }
1072  t_hmaxNearP = spr::chargeIsolationGenCone(nTracks, trackIDs, a_charIsoR_, nNearTRKs, false);
1073  const DetId cellH(trkDetItr.detIdHCAL);
1074  double h5x5 = spr::eHCALmatrix(
1075  theHBHETopology, cellH, hbhe, 2, 2, false, true, -100.0, -100.0, -100.0, -100.0, -100.0, 100.0);
1076  double h7x7 = spr::eHCALmatrix(
1077  theHBHETopology, cellH, hbhe, 3, 3, false, true, -100.0, -100.0, -100.0, -100.0, -100.0, 100.0);
1078  t_hAnnular = h7x7 - h5x5;
1079 #ifdef EDM_ML_DEBUG
1080  edm::LogVerbatim("HcalIsoTrack") << "max p Near (Ecal) " << t_emaxNearP << " (Hcal) " << t_hmaxNearP
1081  << " Annular E (Ecal) " << e11x11P.first << ":" << e15x15P.first << ":"
1082  << t_eAnnular << " (Hcal) " << h5x5 << ":" << h7x7 << ":" << t_hAnnular;
1083 #endif
1084  t_gentrackP = pTrack->p();
1085  if (t_eMipDR < eEcalMax_ && t_hmaxNearP < eIsolation) {
1086  t_DetIds->clear();
1087  t_HitEnergies->clear();
1088  t_DetIds1->clear();
1089  t_HitEnergies1->clear();
1090  t_DetIds3->clear();
1091  t_HitEnergies3->clear();
1092  int nRecHits(-999), nRecHits1(-999), nRecHits3(-999);
1093  std::vector<DetId> ids, ids1, ids3;
1094  std::vector<double> edet0, edet1, edet3;
1095  t_eHcal = spr::eCone_hcal(geo,
1096  hbhe,
1097  trkDetItr.pointHCAL,
1098  trkDetItr.pointECAL,
1099  a_coneR_,
1100  trkDetItr.directionHCAL,
1101  nRecHits,
1102  ids,
1103  edet0,
1104  useRaw_);
1105  if (!oldID_.empty()) {
1106  for (unsigned k = 0; k < ids.size(); ++k)
1107  ids[k] = newId(ids[k]);
1108  }
1109  storeEnergy(0, respCorrs, ids, edet0, t_eHcal, t_DetIds, t_HitEnergies);
1110 
1111  //----- hcal energy in the extended cone 1 (a_coneR+10) --------------
1112  t_eHcal10 = spr::eCone_hcal(geo,
1113  hbhe,
1114  trkDetItr.pointHCAL,
1115  trkDetItr.pointECAL,
1116  a_coneR1_,
1117  trkDetItr.directionHCAL,
1118  nRecHits1,
1119  ids1,
1120  edet1,
1121  useRaw_);
1122  if (!oldID_.empty()) {
1123  for (unsigned k = 0; k < ids1.size(); ++k)
1124  ids1[k] = newId(ids1[k]);
1125  }
1126  storeEnergy(1, respCorrs, ids1, edet1, t_eHcal10, t_DetIds1, t_HitEnergies1);
1127 
1128  //----- hcal energy in the extended cone 3 (a_coneR+30) --------------
1129  t_eHcal30 = spr::eCone_hcal(geo,
1130  hbhe,
1131  trkDetItr.pointHCAL,
1132  trkDetItr.pointECAL,
1133  a_coneR2_,
1134  trkDetItr.directionHCAL,
1135  nRecHits3,
1136  ids3,
1137  edet3,
1138  useRaw_);
1139  if (!oldID_.empty()) {
1140  for (unsigned k = 0; k < ids3.size(); ++k)
1141  ids3[k] = newId(ids3[k]);
1142  }
1143  storeEnergy(3, respCorrs, ids3, edet3, t_eHcal30, t_DetIds3, t_HitEnergies3);
1144 
1145  t_p = pTrack->p();
1146  t_pt = pTrack->pt();
1147  t_phi = pTrack->phi();
1148 
1149 #ifdef EDM_ML_DEBUG
1150  edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) :" << t_pt << "|"
1151  << pTrack->eta() << "|" << t_phi << "|" << t_p << " Generator Level p "
1152  << t_gentrackP;
1153  edm::LogVerbatim("HcalIsoTrack") << "e_MIP " << t_eMipDR << " Chg Isolation " << t_hmaxNearP << " eHcal"
1154  << t_eHcal << " ieta " << t_ieta << " Quality " << t_qltyMissFlag << ":"
1155  << t_qltyPVFlag << ":" << t_selectTk;
1156  for (unsigned int ll = 0; ll < t_DetIds->size(); ll++) {
1157  edm::LogVerbatim("HcalIsoTrack")
1158  << "det id is = " << t_DetIds->at(ll) << " hit enery is = " << t_HitEnergies->at(ll);
1159  }
1160  for (unsigned int ll = 0; ll < t_DetIds1->size(); ll++) {
1161  edm::LogVerbatim("HcalIsoTrack")
1162  << "det id is = " << t_DetIds1->at(ll) << " hit enery is = " << t_HitEnergies1->at(ll);
1163  }
1164  for (unsigned int ll = 0; ll < t_DetIds3->size(); ll++) {
1165  edm::LogVerbatim("HcalIsoTrack")
1166  << "det id is = " << t_DetIds3->at(ll) << " hit enery is = " << t_HitEnergies3->at(ll);
1167  }
1168 #endif
1169  bool accept(false);
1170  if (t_p > pTrackMin_) {
1171  if (t_p < pTrackLow_) {
1172  ++nLow_;
1173  if (prescaleLow_ <= 1)
1174  accept = true;
1175  else if (nLow_ % prescaleLow_ == 1)
1176  accept = true;
1177  } else if (t_p > pTrackHigh_) {
1178  ++nHigh_;
1179  if (prescaleHigh_ <= 1)
1180  accept = true;
1181  else if (nHigh_ % prescaleHigh_ == 1)
1182  accept = true;
1183  } else {
1184  accept = true;
1185  }
1186  }
1187  if (fillTreeRange_) {
1188  if ((t_p < pTrackLow_) || (t_p > pTrackHigh_))
1189  accept = false;
1190  }
1191  if (accept) {
1192  tree->Fill();
1193  nSave++;
1194  int type(0);
1195  if (t_eMipDR < 1.0) {
1196  if (t_hmaxNearP < eIsolate2_) {
1197  ++nLoose;
1198  type = 1;
1199  }
1200  if (t_hmaxNearP < eIsolate1_) {
1201  ++nTight;
1202  type = 2;
1203  }
1204  }
1205  if (t_p > 40.0 && t_p <= 60.0 && t_selectTk) {
1206  t_ietaGood->emplace_back(t_ieta);
1207  t_trackType->emplace_back(type);
1208  }
1209 #ifdef EDM_ML_DEBUG
1210  for (unsigned int k = 0; k < t_trgbits->size(); k++) {
1211  edm::LogVerbatim("HcalIsoTrack") << "trigger bit is = " << t_trgbits->at(k);
1212  }
1213 #endif
1214  }
1215  }
1216  }
1217  ++nTracks;
1218  }
1219  std::array<int, 3> i3{{nSave, nLoose, nTight}};
1220  return i3;
1221 }
double eCone_hcal(const CaloGeometry *geo, edm::Handle< T > &hits, const GlobalPoint &hpoint1, const GlobalPoint &point1, double dR, const GlobalVector &trackMom, int &nRecHits, double hbThr=-100, double heThr=-100, double hfThr=-100, double hoThr=-100, double tMin=-500, double tMax=500, int detOnly=-1, int useRaw=0, bool debug=false)
Log< level::Info, true > LogVerbatim
std::vector< double > * t_HitEnergies3
double pz() const final
z coordinate of momentum vector
double pt() const final
transverse momentum
std::vector< double > * t_HitEnergies
double dR(math::XYZTLorentzVector &, math::XYZTLorentzVector &)
double chargeIsolationGenCone(unsigned int trkIndex, std::vector< spr::propagatedGenParticleID > &trackIDs, double dR, int &nNearTRKs, bool debug=false)
std::vector< int > * t_ietaGood
double chargeIsolationGenEcal(unsigned int trkIndex, std::vector< spr::propagatedGenParticleID > &trackIDs, const CaloGeometry *geo, const CaloTopology *caloTopology, int ieta, int iphi, bool debug=false)
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
std::vector< unsigned int > * t_DetIds1
DetId newId(const DetId &)
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
double rhoh(const edm::Handle< CaloTowerCollection > &)
const std::vector< int > oldID_
bool notaMuon(const reco::GenParticle *pTrack)
double px() const final
x coordinate of momentum vector
double p() const final
magnitude of momentum vector
std::vector< int > * t_ietaAll
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
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)
std::vector< int > * t_trackType
double eCone_ecal(const CaloGeometry *geo, edm::Handle< T > &barrelhits, edm::Handle< T > &endcaphits, const GlobalPoint &hpoint1, const GlobalPoint &point1, double dR, const GlobalVector &trackMom, int &nRecHits, double ebThr=-100, double eeThr=-100, double tMin=-500, double tMax=500, bool debug=false)
double py() const final
y coordinate of momentum vector
std::vector< unsigned int > * t_DetIds3
Definition: DetId.h:17
std::vector< unsigned int > * t_DetIds
std::vector< double > * t_HitEnergies1
std::vector< bool > * t_trgbits
Definition: tree.py:1
double eThreshold(const DetId &id, const CaloGeometry *geo) const
double phi() const final
momentum azimuthal angle
double eECALmatrix(const DetId &detId, edm::Handle< T > &hitsEB, edm::Handle< T > &hitsEE, const CaloGeometry *geo, const CaloTopology *caloTopology, int ieta, int iphi, double ebThr=-100, double eeThr=-100, double tMin=-500, double tMax=500, bool debug=false)
double eHCALmatrix(const HcalTopology *topology, const DetId &det, edm::Handle< T > &hits, int ieta, int iphi, bool includeHO=false, bool algoNew=true, double hbThr=-100, double heThr=-100, double hfThr=-100, double hoThr=-100, double tMin=-500, double tMax=500, int useRaw=0, bool debug=false)
double eta() const final
momentum pseudorapidity

◆ newId()

DetId HcalIsoTrkSimAnalyzer::newId ( const DetId id)
private

Definition at line 1276 of file HcalIsoTrkSimAnalyzer.cc.

References HcalDetId::depth(), hep17_, l1ctLayer2EG_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().

1276  {
1277  HcalDetId hid(id);
1278  if (hep17_ && ((hid.iphi() < 63) || (hid.iphi() > 66) || (hid.zside() < 0)))
1279  return id;
1280  for (unsigned int k = 0; k < oldID_.size(); ++k) {
1281  if ((hid.subdetId() == oldDet_[k]) && (hid.ietaAbs() == oldEta_[k]) && (hid.depth() == oldDepth_[k])) {
1282  return static_cast<DetId>(HcalDetId(hid.subdet(), hid.ieta(), hid.iphi(), newDepth_[k]));
1283  }
1284  }
1285  return id;
1286 }
const std::vector< int > newDepth_
const std::vector< int > oldID_
std::vector< int > oldDepth_
Definition: DetId.h:17

◆ notaMuon()

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

Definition at line 1347 of file HcalIsoTrkSimAnalyzer.cc.

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

Referenced by fillTree().

1347  {
1348  int id = pTrack->pdgId();
1349  bool flag = ((id != 13) && (id != -13) && (pTrack->charge() != 0));
1350 #ifdef EDM_ML_DEBUG
1351  edm::LogVerbatim("HcalIsoTrack") << "notaMuon: ID " << id << " charge " << pTrack->charge() << " Flag " << flag;
1352 #endif
1353  return flag;
1354 }
Log< level::Info, true > LogVerbatim
int pdgId() const final
PDG identifier.
int charge() const final
electric charge

◆ rhoh()

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

Definition at line 1227 of file HcalIsoTrkSimAnalyzer.cc.

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

Referenced by fillTree().

1227  {
1228  std::vector<double> sumPFNallSMDQH2;
1229  sumPFNallSMDQH2.reserve(phibins_.size() * etabins_.size());
1230 
1231  for (auto eta : etabins_) {
1232  for (auto phi : phibins_) {
1233  double hadder = 0;
1234  for (const auto& pf_it : (*tower)) {
1235  if (fabs(eta - pf_it.eta()) > etahalfdist_)
1236  continue;
1237  if (fabs(reco::deltaPhi(phi, pf_it.phi())) > phihalfdist_)
1238  continue;
1239  hadder += pf_it.hadEt();
1240  }
1241  sumPFNallSMDQH2.emplace_back(hadder);
1242  }
1243  }
1244 
1245  double evt_smdq(0);
1246  std::sort(sumPFNallSMDQH2.begin(), sumPFNallSMDQH2.end());
1247  if (sumPFNallSMDQH2.size() % 2)
1248  evt_smdq = sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 1) / 2];
1249  else
1250  evt_smdq = (sumPFNallSMDQH2[sumPFNallSMDQH2.size() / 2] + sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 2) / 2]) / 2.;
1251  double rhoh = evt_smdq / (etadist_ * phidist_);
1252 #ifdef EDM_ML_DEBUG
1253  edm::LogVerbatim("HcalIsoTrack") << "Rho " << evt_smdq << ":" << rhoh;
1254 #endif
1255  return rhoh;
1256 }
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:26
Log< level::Info, true > LogVerbatim
std::vector< double > etabins_
std::vector< double > phibins_
double rhoh(const edm::Handle< CaloTowerCollection > &)

◆ storeEnergy()

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

Definition at line 1288 of file HcalIsoTrkSimAnalyzer.cc.

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

Referenced by fillTree().

1294  {
1295  double ehcal(0);
1296  if (unCorrect_) {
1297  for (unsigned int k = 0; k < ids.size(); ++k) {
1298  double corr = (respCorrs->getValues(ids[k]))->getValue();
1299  if (corr != 0)
1300  edet[k] /= corr;
1301  ehcal += edet[k];
1302  }
1303  } else {
1304  for (const auto& en : edet)
1305  ehcal += en;
1306  }
1307  if ((std::abs(ehcal - eHcal) > 0.001) && (!unCorrect_))
1308  edm::LogWarning("HcalIsoTrack") << "Check inconsistent energies: " << indx << " " << eHcal << ":" << ehcal
1309  << " from " << ids.size() << " cells";
1310  eHcal = hcalScale_ * ehcal;
1311 
1312  if (collapseDepth_) {
1313  std::map<HcalDetId, double> hitMap;
1314  for (unsigned int k = 0; k < ids.size(); ++k) {
1315  HcalDetId id = hdc_->mergedDepthDetId(HcalDetId(ids[k]));
1316  auto itr = hitMap.find(id);
1317  if (itr == hitMap.end()) {
1318  hitMap[id] = edet[k];
1319  } else {
1320  (itr->second) += edet[k];
1321  }
1322  }
1323  detIds->reserve(hitMap.size());
1324  hitEnergies->reserve(hitMap.size());
1325  for (const auto& hit : hitMap) {
1326  detIds->emplace_back(hit.first.rawId());
1327  hitEnergies->emplace_back(hit.second);
1328  }
1329  } else {
1330  detIds->reserve(ids.size());
1331  hitEnergies->reserve(ids.size());
1332  for (unsigned int k = 0; k < ids.size(); ++k) {
1333  detIds->emplace_back(ids[k].rawId());
1334  hitEnergies->emplace_back(edet[k]);
1335  }
1336  }
1337 #ifdef EDM_ML_DEBUG
1338  edm::LogVerbatim("HcalIsoTrack") << "Input to storeEnergy with " << ids.size() << " cells";
1339  for (unsigned int k = 0; k < ids.size(); ++k)
1340  edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId(ids[k]) << " E " << edet[k];
1341  edm::LogVerbatim("HcalIsoTrack") << "Output of storeEnergy with " << detIds->size() << " cells and Etot " << eHcal;
1342  for (unsigned int k = 0; k < detIds->size(); ++k)
1343  edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId((*detIds)[k]) << " E " << (*hitEnergies)[k];
1344 #endif
1345 }
Log< level::Info, true > LogVerbatim
HcalDetId mergedDepthDetId(const HcalDetId &id) const
const Item * getValues(DetId fId, bool throwOnFail=true) const
dictionary corr
const HcalDDDRecConstants * hdc_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
const JetExtendedData & getValue(const Container &, const reco::JetBaseRef &)
get value for the association. Throw exception if no association found
Log< level::Warning, false > LogWarning

Member Data Documentation

◆ a_charIsoR_

double HcalIsoTrkSimAnalyzer::a_charIsoR_
private

Definition at line 169 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_coneR1_

double HcalIsoTrkSimAnalyzer::a_coneR1_
private

Definition at line 169 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_coneR2_

double HcalIsoTrkSimAnalyzer::a_coneR2_
private

Definition at line 169 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_coneR_

const double HcalIsoTrkSimAnalyzer::a_coneR_
private

Definition at line 124 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_mipR2_

const double HcalIsoTrkSimAnalyzer::a_mipR2_
private

Definition at line 124 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_mipR3_

const double HcalIsoTrkSimAnalyzer::a_mipR3_
private

Definition at line 124 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_mipR4_

const double HcalIsoTrkSimAnalyzer::a_mipR4_
private

Definition at line 125 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_mipR5_

const double HcalIsoTrkSimAnalyzer::a_mipR5_
private

Definition at line 125 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ a_mipR_

const double HcalIsoTrkSimAnalyzer::a_mipR_
private

Definition at line 124 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ algTag_

const edm::InputTag HcalIsoTrkSimAnalyzer::algTag_
private

Definition at line 143 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ collapseDepth_

const bool HcalIsoTrkSimAnalyzer::collapseDepth_
private

Definition at line 132 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and storeEnergy().

◆ dataType_

const int HcalIsoTrkSimAnalyzer::dataType_
private

Definition at line 130 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ eEcalMax_

const double HcalIsoTrkSimAnalyzer::eEcalMax_
private

Definition at line 125 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ eIsolate1_

const double HcalIsoTrkSimAnalyzer::eIsolate1_
private

Definition at line 127 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ eIsolate2_

const double HcalIsoTrkSimAnalyzer::eIsolate2_
private

Definition at line 127 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ etabins_

std::vector<double> HcalIsoTrkSimAnalyzer::etabins_
private

Definition at line 173 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ etadist_

double HcalIsoTrkSimAnalyzer::etadist_
private

Definition at line 175 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ etahalfdist_

double HcalIsoTrkSimAnalyzer::etahalfdist_
private

Definition at line 175 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ etaMax_

const double HcalIsoTrkSimAnalyzer::etaMax_
private

Definition at line 123 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ eThresholds_

const EcalPFRecHitThresholds* HcalIsoTrkSimAnalyzer::eThresholds_
private

Definition at line 171 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ extTag_

const edm::InputTag HcalIsoTrkSimAnalyzer::extTag_
private

Definition at line 143 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ fillTreeRange_

const bool HcalIsoTrkSimAnalyzer::fillTreeRange_
private

Definition at line 132 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ hcalScale_

const double HcalIsoTrkSimAnalyzer::hcalScale_
private

Definition at line 127 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and storeEnergy().

◆ hdc_

const HcalDDDRecConstants* HcalIsoTrkSimAnalyzer::hdc_
private

Definition at line 170 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginRun(), and storeEnergy().

◆ hep17_

const bool HcalIsoTrkSimAnalyzer::hep17_
private

Definition at line 140 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and newId().

◆ hitEthrEB_

const double HcalIsoTrkSimAnalyzer::hitEthrEB_
private

Definition at line 133 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEE0_

const double HcalIsoTrkSimAnalyzer::hitEthrEE0_
private

Definition at line 133 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEE1_

const double HcalIsoTrkSimAnalyzer::hitEthrEE1_
private

Definition at line 133 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEE2_

const double HcalIsoTrkSimAnalyzer::hitEthrEE2_
private

Definition at line 134 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEE3_

const double HcalIsoTrkSimAnalyzer::hitEthrEE3_
private

Definition at line 134 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEEHi_

const double HcalIsoTrkSimAnalyzer::hitEthrEEHi_
private

Definition at line 135 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hitEthrEELo_

const double HcalIsoTrkSimAnalyzer::hitEthrEELo_
private

Definition at line 135 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ hltConfig_

HLTConfigProvider HcalIsoTrkSimAnalyzer::hltConfig_
private

Definition at line 119 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginRun().

◆ ignoreTrigger_

const bool HcalIsoTrkSimAnalyzer::ignoreTrigger_
private

Definition at line 131 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ l1Filter_

const std::string HcalIsoTrkSimAnalyzer::l1Filter_
private

Definition at line 121 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ l1GtUtils_

l1t::L1TGlobalUtil* HcalIsoTrkSimAnalyzer::l1GtUtils_
private

Definition at line 118 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ l1TrigName_

const std::string HcalIsoTrkSimAnalyzer::l1TrigName_
private

Definition at line 138 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ l2Filter_

const std::string HcalIsoTrkSimAnalyzer::l2Filter_
private

Definition at line 122 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ l3Filter_

const std::string HcalIsoTrkSimAnalyzer::l3Filter_
private

Definition at line 122 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ labelBS_

const std::string HcalIsoTrkSimAnalyzer::labelBS_
private

Definition at line 142 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ labelEB_

const std::string HcalIsoTrkSimAnalyzer::labelEB_
private

Definition at line 137 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ labelEE_

const std::string HcalIsoTrkSimAnalyzer::labelEE_
private

Definition at line 138 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ labelGenTrack_

const std::string HcalIsoTrkSimAnalyzer::labelGenTrack_
private

Definition at line 137 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ labelHBHE_

const std::string HcalIsoTrkSimAnalyzer::labelHBHE_
private

Definition at line 138 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ labelRecVtx_

const std::string HcalIsoTrkSimAnalyzer::labelRecVtx_
private

Definition at line 137 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ labelTower_

const std::string HcalIsoTrkSimAnalyzer::labelTower_
private

Definition at line 138 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ maxRestrictionP_

const double HcalIsoTrkSimAnalyzer::maxRestrictionP_
private

Definition at line 126 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ mode_

const int HcalIsoTrkSimAnalyzer::mode_
private

Definition at line 130 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ modnam_

const std::string HcalIsoTrkSimAnalyzer::modnam_
private

Definition at line 142 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ newDepth_

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

Definition at line 139 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and newId().

◆ nHigh_

unsigned int HcalIsoTrkSimAnalyzer::nHigh_
private

Definition at line 168 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree().

◆ nLow_

unsigned int HcalIsoTrkSimAnalyzer::nLow_
private

Definition at line 168 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree().

◆ nRun_

unsigned int HcalIsoTrkSimAnalyzer::nRun_
private

Definition at line 168 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginRun(), and endRun().

◆ oldDepth_

std::vector<int> HcalIsoTrkSimAnalyzer::oldDepth_
private

Definition at line 174 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and newId().

◆ oldDet_

std::vector<int> HcalIsoTrkSimAnalyzer::oldDet_
private

Definition at line 174 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and newId().

◆ oldEta_

std::vector<int> HcalIsoTrkSimAnalyzer::oldEta_
private

Definition at line 174 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and newId().

◆ oldID_

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

Definition at line 139 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ phibins_

std::vector<double> HcalIsoTrkSimAnalyzer::phibins_
private

Definition at line 173 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ phidist_

double HcalIsoTrkSimAnalyzer::phidist_
private

Definition at line 175 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ phihalfdist_

double HcalIsoTrkSimAnalyzer::phihalfdist_
private

Definition at line 175 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and rhoh().

◆ prdnam_

const std::string HcalIsoTrkSimAnalyzer::prdnam_
private

Definition at line 142 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ prescaleHigh_

const int HcalIsoTrkSimAnalyzer::prescaleHigh_
private

Definition at line 129 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ prescaleLow_

const int HcalIsoTrkSimAnalyzer::prescaleLow_
private

Definition at line 129 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ processName_

const std::string HcalIsoTrkSimAnalyzer::processName_
private

Definition at line 121 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkSimAnalyzer().

◆ ptMin_

const double HcalIsoTrkSimAnalyzer::ptMin_
private

Definition at line 123 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ pTrackHigh_

const double HcalIsoTrkSimAnalyzer::pTrackHigh_
private

Definition at line 128 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ pTrackLow_

const double HcalIsoTrkSimAnalyzer::pTrackLow_
private

Definition at line 128 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ pTrackMin_

const double HcalIsoTrkSimAnalyzer::pTrackMin_
private

Definition at line 125 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ slopeRestrictionP_

const double HcalIsoTrkSimAnalyzer::slopeRestrictionP_
private

Definition at line 126 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().

◆ t_allvertex

int HcalIsoTrkSimAnalyzer::t_allvertex
private

Definition at line 196 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_DataType

int HcalIsoTrkSimAnalyzer::t_DataType
private

Definition at line 179 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_DetIds

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

Definition at line 191 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_DetIds1

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

Definition at line 191 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_DetIds3

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

Definition at line 191 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eAnnular

double HcalIsoTrkSimAnalyzer::t_eAnnular
private

Definition at line 187 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal

double HcalIsoTrkSimAnalyzer::t_eHcal
private

Definition at line 188 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal10

double HcalIsoTrkSimAnalyzer::t_eHcal10
private

Definition at line 188 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal30

double HcalIsoTrkSimAnalyzer::t_eHcal30
private

Definition at line 188 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_emaxNearP

double HcalIsoTrkSimAnalyzer::t_emaxNearP
private

Definition at line 187 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR

double HcalIsoTrkSimAnalyzer::t_eMipDR
private

Definition at line 185 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR2

double HcalIsoTrkSimAnalyzer::t_eMipDR2
private

Definition at line 185 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR3

double HcalIsoTrkSimAnalyzer::t_eMipDR3
private

Definition at line 185 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR4

double HcalIsoTrkSimAnalyzer::t_eMipDR4
private

Definition at line 185 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR5

double HcalIsoTrkSimAnalyzer::t_eMipDR5
private

Definition at line 186 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Event

int HcalIsoTrkSimAnalyzer::t_Event
private

Definition at line 179 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_EventNo

unsigned int HcalIsoTrkSimAnalyzer::t_EventNo
private

Definition at line 178 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_EventWeight

double HcalIsoTrkSimAnalyzer::t_EventWeight
private

Definition at line 181 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_gentrackP

double HcalIsoTrkSimAnalyzer::t_gentrackP
private

Definition at line 186 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_goodPV

int HcalIsoTrkSimAnalyzer::t_goodPV
private

Definition at line 180 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_hAnnular

double HcalIsoTrkSimAnalyzer::t_hAnnular
private

Definition at line 187 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies

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

Definition at line 192 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies1

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

Definition at line 192 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies3

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

Definition at line 192 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_hltbits

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

Definition at line 193 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_hmaxNearP

double HcalIsoTrkSimAnalyzer::t_hmaxNearP
private

Definition at line 186 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_ieta

int HcalIsoTrkSimAnalyzer::t_ieta
private

Definition at line 179 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_ietaAll

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

Definition at line 197 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ t_ietaGood

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

Definition at line 197 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ t_iphi

int HcalIsoTrkSimAnalyzer::t_iphi
private

Definition at line 179 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_L1Bit

bool HcalIsoTrkSimAnalyzer::t_L1Bit
private

Definition at line 194 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1eta

double HcalIsoTrkSimAnalyzer::t_l1eta
private

Definition at line 182 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1phi

double HcalIsoTrkSimAnalyzer::t_l1phi
private

Definition at line 182 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1pt

double HcalIsoTrkSimAnalyzer::t_l1pt
private

Definition at line 182 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3eta

double HcalIsoTrkSimAnalyzer::t_l3eta
private

Definition at line 183 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3phi

double HcalIsoTrkSimAnalyzer::t_l3phi
private

Definition at line 183 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3pt

double HcalIsoTrkSimAnalyzer::t_l3pt
private

Definition at line 183 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_mindR1

double HcalIsoTrkSimAnalyzer::t_mindR1
private

Definition at line 184 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_mindR2

double HcalIsoTrkSimAnalyzer::t_mindR2
private

Definition at line 184 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_nTrk

int HcalIsoTrkSimAnalyzer::t_nTrk
private

Definition at line 180 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_nVtx

int HcalIsoTrkSimAnalyzer::t_nVtx
private

Definition at line 180 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_p

double HcalIsoTrkSimAnalyzer::t_p
private

Definition at line 181 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_phi

double HcalIsoTrkSimAnalyzer::t_phi
private

Definition at line 181 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_pt

double HcalIsoTrkSimAnalyzer::t_pt
private

Definition at line 181 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyFlag

bool HcalIsoTrkSimAnalyzer::t_qltyFlag
private

Definition at line 189 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyMissFlag

bool HcalIsoTrkSimAnalyzer::t_qltyMissFlag
private

Definition at line 189 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyPVFlag

bool HcalIsoTrkSimAnalyzer::t_qltyPVFlag
private

Definition at line 190 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_rhoh

double HcalIsoTrkSimAnalyzer::t_rhoh
private

Definition at line 188 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Run

int HcalIsoTrkSimAnalyzer::t_Run
private

Definition at line 179 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_RunNo

unsigned int HcalIsoTrkSimAnalyzer::t_RunNo
private

Definition at line 178 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_selectTk

bool HcalIsoTrkSimAnalyzer::t_selectTk
private

Definition at line 189 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Tracks

int HcalIsoTrkSimAnalyzer::t_Tracks
private

Definition at line 195 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksLoose

int HcalIsoTrkSimAnalyzer::t_TracksLoose
private

Definition at line 196 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksProp

int HcalIsoTrkSimAnalyzer::t_TracksProp
private

Definition at line 195 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksSaved

int HcalIsoTrkSimAnalyzer::t_TracksSaved
private

Definition at line 195 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksTight

int HcalIsoTrkSimAnalyzer::t_TracksTight
private

Definition at line 196 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_trackType

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

Definition at line 197 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ t_trgbits

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

Definition at line 193 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ t_TrigPass

bool HcalIsoTrkSimAnalyzer::t_TrigPass
private

Definition at line 190 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TrigPassSel

bool HcalIsoTrkSimAnalyzer::t_TrigPassSel
private

Definition at line 190 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ theTriggerResultsLabel_

const edm::InputTag HcalIsoTrkSimAnalyzer::theTriggerResultsLabel_
private

Definition at line 136 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer().

◆ tok_alg_

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

Definition at line 155 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_bFieldH_

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

Definition at line 158 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_bs_

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

Definition at line 147 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_cala_

edm::EDGetTokenT<CaloTowerCollection> HcalIsoTrkSimAnalyzer::tok_cala_
private

Definition at line 154 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_caloTopology_

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

Definition at line 162 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_ddrec_

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

Definition at line 157 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by beginRun().

◆ tok_EB_

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkSimAnalyzer::tok_EB_
private

Definition at line 149 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_ecalChStatus_

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

Definition at line 159 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_ecalPFRecHitThresholds_

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

Definition at line 166 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_EE_

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkSimAnalyzer::tok_EE_
private

Definition at line 150 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_ew_

edm::EDGetTokenT<GenEventInfoProduct> HcalIsoTrkSimAnalyzer::tok_ew_
private

Definition at line 152 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_geom_

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

Definition at line 161 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> HcalIsoTrkSimAnalyzer::tok_hbhe_
private

Definition at line 151 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_htopo_

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

Definition at line 163 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_parts_

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

Definition at line 153 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_pdt_

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

Definition at line 165 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_recVtx_

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

Definition at line 148 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_resp_

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

Definition at line 164 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_sevlv_

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

Definition at line 160 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_trigEvt_

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

Definition at line 145 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tok_trigRes_

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

Definition at line 146 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze().

◆ tree

TTree* HcalIsoTrkSimAnalyzer::tree
private

◆ tree2

TTree * HcalIsoTrkSimAnalyzer::tree2
private

Definition at line 177 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ triggerEvent_

const edm::InputTag HcalIsoTrkSimAnalyzer::triggerEvent_
private

Definition at line 136 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ trigNames_

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

Definition at line 120 of file HcalIsoTrkSimAnalyzer.cc.

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

◆ unCorrect_

const bool HcalIsoTrkSimAnalyzer::unCorrect_
private

Definition at line 132 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by HcalIsoTrkSimAnalyzer(), and storeEnergy().

◆ useL1Trigger_

const bool HcalIsoTrkSimAnalyzer::useL1Trigger_
private

Definition at line 131 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkSimAnalyzer().

◆ usePFThresh_

const bool HcalIsoTrkSimAnalyzer::usePFThresh_
private

Definition at line 141 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkSimAnalyzer().

◆ useRaw_

const int HcalIsoTrkSimAnalyzer::useRaw_
private

Definition at line 130 of file HcalIsoTrkSimAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkSimAnalyzer().