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

Public Member Functions

 HcalIsoTrkAnalyzer (edm::ParameterSet const &)
 
 ~HcalIsoTrkAnalyzer () override
 
- Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
 EDAnalyzer ()=default
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

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

Private Member Functions

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

Private Attributes

double a_charIsoR_
 
double a_coneR1_
 
double a_coneR2_
 
const double a_coneR_
 
const double a_mipR2_
 
const double a_mipR3_
 
const double a_mipR4_
 
const double a_mipR5_
 
const double a_mipR_
 
const bool collapseDepth_
 
const int dataType_
 
const double eEcalMax_
 
const double eIsolate1_
 
const double eIsolate2_
 
std::vector< double > etabins_
 
double etadist_
 
double etahalfdist_
 
edm::Service< TFileService > fs
 
const double hcalScale_
 
const HcalDDDRecConstants * hdc_
 
const double hitEthrEB_
 
const double hitEthrEE0_
 
const double hitEthrEE1_
 
const double hitEthrEE2_
 
const double hitEthrEE3_
 
const double hitEthrEEHi_
 
const double hitEthrEELo_
 
HLTConfigProvider hltConfig_
 
const bool ignoreTrigger_
 
const std::string l1Filter_
 
l1t::L1TGlobalUtil * l1GtUtils_
 
const std::string l1TrigName_
 
const std::string l2Filter_
 
const std::string l3Filter_
 
const std::string labelEB_
 
const std::string labelEE_
 
const std::string labelGenTrack_
 
const std::string labelHBHE_
 
const std::string labelRecVtx_
 
const std::string labelTower_
 
const double maxRestrictionP_
 
unsigned int nHigh_
 
unsigned int nLow_
 
unsigned int nRun_
 
std::vector< double > phibins_
 
double phidist_
 
double phihalfdist_
 
const int prescaleHigh_
 
const int prescaleLow_
 
const std::string processName_
 
const double pTrackHigh_
 
const double pTrackLow_
 
const double pTrackMin_
 
spr::trackSelectionParameters selectionParameter_
 
const double slopeRestrictionP_
 
int t_allvertex
 
int t_DataType
 
std::vector< unsigned int > * t_DetIds
 
std::vector< unsigned int > * t_DetIds1
 
std::vector< unsigned int > * t_DetIds3
 
double t_eAnnular
 
double t_eHcal
 
double t_eHcal10
 
double t_eHcal30
 
double t_emaxNearP
 
double t_eMipDR
 
double t_eMipDR2
 
double t_eMipDR3
 
double t_eMipDR4
 
double t_eMipDR5
 
int t_Event
 
unsigned int t_EventNo
 
double t_EventWeight
 
double t_gentrackP
 
int t_goodPV
 
double t_hAnnular
 
std::vector< double > * t_HitEnergies
 
std::vector< double > * t_HitEnergies1
 
std::vector< double > * t_HitEnergies3
 
std::vector< bool > * t_hltbits
 
double t_hmaxNearP
 
int t_ieta
 
std::vector< int > * t_ietaAll
 
std::vector< int > * t_ietaGood
 
int t_iphi
 
bool t_L1Bit
 
double t_l1eta
 
double t_l1phi
 
double t_l1pt
 
double t_l3eta
 
double t_l3phi
 
double t_l3pt
 
double t_mindR1
 
double t_mindR2
 
int t_nTrk
 
int t_nVtx
 
double t_p
 
double t_phi
 
double t_pt
 
bool t_qltyFlag
 
bool t_qltyMissFlag
 
bool t_qltyPVFlag
 
double t_rhoh
 
int t_Run
 
unsigned int t_RunNo
 
bool t_selectTk
 
int t_Tracks
 
int t_TracksLoose
 
int t_TracksProp
 
int t_TracksSaved
 
int t_TracksTight
 
std::vector< int > * t_trackType
 
std::vector< bool > * t_trgbits
 
bool t_TrigPass
 
bool t_TrigPassSel
 
const std::string theTrackQuality_
 
const edm::InputTag theTriggerResultsLabel_
 
edm::EDGetTokenT< BXVector< GlobalAlgBlk > > tok_alg_
 
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
 
edm::EDGetTokenT< CaloTowerCollection > tok_cala_
 
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
 
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
 
edm::EDGetTokenT< GenEventInfoProduct > tok_ew_
 
edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
 
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
 
edm::EDGetTokenT< reco::GenParticleCollection > tok_parts_
 
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
 
edm::EDGetTokenT< trigger::TriggerEvent > tok_trigEvt_
 
edm::EDGetTokenT< edm::TriggerResults > tok_trigRes_
 
TTree * tree
 
TTree * tree2
 
const edm::InputTag triggerEvent_
 
const std::vector< std::string > trigNames_
 
const bool unCorrect_
 
const bool useL1Trigger_
 
const int useRaw_
 

Additional Inherited Members

- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 87 of file HcalIsoTrkAnalyzer.cc.

Constructor & Destructor Documentation

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

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

References a_charIsoR_, a_coneR1_, a_coneR2_, a_coneR_, a_mipR2_, a_mipR3_, a_mipR4_, a_mipR5_, a_mipR_, collapseDepth_, edm::EDConsumerBase::consumesCollector(), dataType_, eEcalMax_, eIsolate1_, eIsolate2_, PVValHelper::eta, etabins_, etadist_, etahalfdist_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcalScale_, hitEthrEB_, hitEthrEE0_, hitEthrEE1_, hitEthrEE2_, hitEthrEE3_, hitEthrEEHi_, hitEthrEELo_, mps_fire::i, ignoreTrigger_, HLT_2018_cff::InputTag, dqmdumpme::k, TFileService::kSharedResource, l1Filter_, l1GtUtils_, l1TrigName_, l2Filter_, l3Filter_, labelEB_, labelEE_, labelGenTrack_, labelHBHE_, labelRecVtx_, labelTower_, M_PI, spr::trackSelectionParameters::maxChi2, spr::trackSelectionParameters::maxDpOverP, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, maxRestrictionP_, spr::trackSelectionParameters::minLayerCrossed, spr::trackSelectionParameters::minOuterHit, spr::trackSelectionParameters::minPt, spr::trackSelectionParameters::minQuality, phi, phibins_, phidist_, phihalfdist_, prescaleHigh_, prescaleLow_, processName_, pTrackHigh_, pTrackLow_, pTrackMin_, reco::TrackBase::qualityByName(), selectionParameter_, slopeRestrictionP_, AlCaHLTBitMon_QueryRunRegistry::string, theTrackQuality_, theTriggerResultsLabel_, tok_alg_, tok_bs_, tok_cala_, tok_EB_, tok_EE_, tok_ew_, tok_genTrack_, tok_hbhe_, tok_parts_, tok_recVtx_, tok_trigEvt_, tok_trigRes_, triggerEvent_, trigNames_, unCorrect_, useL1Trigger_, and useRaw_.

193  : trigNames_(iConfig.getParameter<std::vector<std::string>>("triggers")),
194  theTrackQuality_(iConfig.getParameter<std::string>("trackQuality")),
195  processName_(iConfig.getParameter<std::string>("processName")),
196  l1Filter_(iConfig.getParameter<std::string>("l1Filter")),
197  l2Filter_(iConfig.getParameter<std::string>("l2Filter")),
198  l3Filter_(iConfig.getParameter<std::string>("l3Filter")),
199  a_coneR_(iConfig.getParameter<double>("coneRadius")),
200  a_mipR_(iConfig.getParameter<double>("coneRadiusMIP")),
201  a_mipR2_(iConfig.getParameter<double>("coneRadiusMIP2")),
202  a_mipR3_(iConfig.getParameter<double>("coneRadiusMIP3")),
203  a_mipR4_(iConfig.getParameter<double>("coneRadiusMIP4")),
204  a_mipR5_(iConfig.getParameter<double>("coneRadiusMIP5")),
205  pTrackMin_(iConfig.getParameter<double>("minimumTrackP")),
206  eEcalMax_(iConfig.getParameter<double>("maximumEcalEnergy")),
207  maxRestrictionP_(iConfig.getParameter<double>("maxTrackP")),
208  slopeRestrictionP_(iConfig.getParameter<double>("slopeTrackP")),
209  hcalScale_(iConfig.getUntrackedParameter<double>("hHcalScale", 1.0)),
210  eIsolate1_(iConfig.getParameter<double>("isolationEnergyTight")),
211  eIsolate2_(iConfig.getParameter<double>("isolationEnergyLoose")),
212  pTrackLow_(iConfig.getParameter<double>("momentumLow")),
213  pTrackHigh_(iConfig.getParameter<double>("momentumHigh")),
214  prescaleLow_(iConfig.getParameter<int>("prescaleLow")),
215  prescaleHigh_(iConfig.getParameter<int>("prescaleHigh")),
216  useRaw_(iConfig.getUntrackedParameter<int>("useRaw", 0)),
217  dataType_(iConfig.getUntrackedParameter<int>("dataType", 0)),
218  ignoreTrigger_(iConfig.getUntrackedParameter<bool>("ignoreTriggers", false)),
219  useL1Trigger_(iConfig.getUntrackedParameter<bool>("useL1Trigger", false)),
220  unCorrect_(iConfig.getUntrackedParameter<bool>("unCorrect", false)),
221  collapseDepth_(iConfig.getUntrackedParameter<bool>("collapseDepth", false)),
222  hitEthrEB_(iConfig.getParameter<double>("EBHitEnergyThreshold")),
223  hitEthrEE0_(iConfig.getParameter<double>("EEHitEnergyThreshold0")),
224  hitEthrEE1_(iConfig.getParameter<double>("EEHitEnergyThreshold1")),
225  hitEthrEE2_(iConfig.getParameter<double>("EEHitEnergyThreshold2")),
226  hitEthrEE3_(iConfig.getParameter<double>("EEHitEnergyThreshold3")),
227  hitEthrEELo_(iConfig.getParameter<double>("EEHitEnergyThresholdLow")),
228  hitEthrEEHi_(iConfig.getParameter<double>("EEHitEnergyThresholdHigh")),
229  triggerEvent_(iConfig.getParameter<edm::InputTag>("labelTriggerEvent")),
230  theTriggerResultsLabel_(iConfig.getParameter<edm::InputTag>("labelTriggerResult")),
231  labelGenTrack_(iConfig.getParameter<std::string>("labelTrack")),
232  labelRecVtx_(iConfig.getParameter<std::string>("labelVertex")),
233  labelEB_(iConfig.getParameter<std::string>("labelEBRecHit")),
234  labelEE_(iConfig.getParameter<std::string>("labelEERecHit")),
235  labelHBHE_(iConfig.getParameter<std::string>("labelHBHERecHit")),
236  labelTower_(iConfig.getParameter<std::string>("labelCaloTower")),
237  l1TrigName_(iConfig.getUntrackedParameter<std::string>("l1TrigName", "L1_SingleJet60")),
238  nRun_(0),
239  nLow_(0),
240  nHigh_(0),
241  hdc_(nullptr) {
242  usesResource(TFileService::kSharedResource);
243 
244  //now do whatever initialization is needed
245  const double isolationRadius(28.9), innerR(10.0), outerR(30.0);
246  reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality_);
247  selectionParameter_.minPt = iConfig.getParameter<double>("minTrackPt");
248  ;
249  selectionParameter_.minQuality = trackQuality_;
250  selectionParameter_.maxDxyPV = iConfig.getParameter<double>("maxDxyPV");
251  selectionParameter_.maxDzPV = iConfig.getParameter<double>("maxDzPV");
252  selectionParameter_.maxChi2 = iConfig.getParameter<double>("maxChi2");
253  selectionParameter_.maxDpOverP = iConfig.getParameter<double>("maxDpOverP");
254  selectionParameter_.minOuterHit = iConfig.getParameter<int>("minOuterHit");
255  selectionParameter_.minLayerCrossed = iConfig.getParameter<int>("minLayerCrossed");
256  selectionParameter_.maxInMiss = iConfig.getParameter<int>("maxInMiss");
257  selectionParameter_.maxOutMiss = iConfig.getParameter<int>("maxOutMiss");
258  a_charIsoR_ = a_coneR_ + isolationRadius;
259  a_coneR1_ = a_coneR_ + innerR;
260  a_coneR2_ = a_coneR_ + outerR;
261  // Different isolation cuts are described in DN-2016/029
262  // Tight cut uses 2 GeV; Loose cut uses 10 GeV
263  // Eta dependent cut uses (maxRestrictionP_ * exp(|ieta|*log(2.5)/18))
264  // with the factor for exponential slopeRestrictionP_ = log(2.5)/18
265  // maxRestrictionP_ = 8 GeV as came from a study
266  std::string labelBS = iConfig.getParameter<std::string>("labelBeamSpot");
267  std::string modnam = iConfig.getUntrackedParameter<std::string>("moduleName", "");
268  std::string prdnam = iConfig.getUntrackedParameter<std::string>("producerName", "");
269  edm::InputTag algTag = iConfig.getParameter<edm::InputTag>("algInputTag");
270  edm::InputTag extTag = iConfig.getParameter<edm::InputTag>("extInputTag");
271  l1GtUtils_ = new l1t::L1TGlobalUtil(iConfig, consumesCollector(), *this, algTag, extTag);
272  // define tokens for access
273  tok_trigEvt_ = consumes<trigger::TriggerEvent>(triggerEvent_);
274  tok_trigRes_ = consumes<edm::TriggerResults>(theTriggerResultsLabel_);
275  tok_bs_ = consumes<reco::BeamSpot>(labelBS);
276  tok_genTrack_ = consumes<reco::TrackCollection>(labelGenTrack_);
277  tok_ew_ = consumes<GenEventInfoProduct>(edm::InputTag("generator"));
278  tok_parts_ = consumes<reco::GenParticleCollection>(edm::InputTag("genParticles"));
279  tok_cala_ = consumes<CaloTowerCollection>(labelTower_);
280  tok_alg_ = consumes<BXVector<GlobalAlgBlk>>(algTag);
281 
282  if (modnam.empty()) {
283  tok_recVtx_ = consumes<reco::VertexCollection>(labelRecVtx_);
284  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEB_));
285  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEE_));
286  tok_hbhe_ = consumes<HBHERecHitCollection>(labelHBHE_);
287  edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
288  << labelBS << " " << labelRecVtx_ << " " << labelGenTrack_ << " "
289  << edm::InputTag("ecalRecHit", labelEB_) << " "
290  << edm::InputTag("ecalRecHit", labelEE_) << " " << labelHBHE_ << " "
291  << labelTower_;
292  } else {
293  tok_recVtx_ = consumes<reco::VertexCollection>(edm::InputTag(modnam, labelRecVtx_, prdnam));
294  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag(modnam, labelEB_, prdnam));
295  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag(modnam, labelEE_, prdnam));
296  tok_hbhe_ = consumes<HBHERecHitCollection>(edm::InputTag(modnam, labelHBHE_, prdnam));
297  edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
298  << labelBS << " " << edm::InputTag(modnam, labelRecVtx_, prdnam) << " "
299  << labelGenTrack_ << " " << edm::InputTag(modnam, labelEB_, prdnam) << " "
300  << edm::InputTag(modnam, labelEE_, prdnam) << " "
301  << edm::InputTag(modnam, labelHBHE_, prdnam) << " " << labelTower_;
302  }
303 
304  edm::LogVerbatim("HcalIsoTrack")
305  << "Parameters read from config file \n"
306  << "\t minPt " << selectionParameter_.minPt << "\t theTrackQuality " << theTrackQuality_ << "\t minQuality "
307  << selectionParameter_.minQuality << "\t maxDxyPV " << selectionParameter_.maxDxyPV << "\t maxDzPV "
308  << selectionParameter_.maxDzPV << "\t maxChi2 " << selectionParameter_.maxChi2 << "\t maxDpOverP "
309  << selectionParameter_.maxDpOverP << "\t minOuterHit " << selectionParameter_.minOuterHit << "\t minLayerCrossed "
310  << selectionParameter_.minLayerCrossed << "\t maxInMiss " << selectionParameter_.maxInMiss << "\t maxOutMiss "
311  << selectionParameter_.maxOutMiss << "\t a_coneR " << a_coneR_ << ":" << a_coneR1_ << ":" << a_coneR2_
312  << "\t a_charIsoR " << a_charIsoR_ << "\t a_mipR " << a_mipR_ << "\t a_mipR2 " << a_mipR2_ << "\t a_mipR3 "
313  << a_mipR3_ << "\t a_mipR4 " << a_mipR4_ << "\t a_mipR5 " << a_mipR5_ << "\n pTrackMin_ " << pTrackMin_
314  << "\t eEcalMax_ " << eEcalMax_ << "\t maxRestrictionP_ " << maxRestrictionP_ << "\t slopeRestrictionP_ "
315  << slopeRestrictionP_ << "\t eIsolateStrong_ " << eIsolate1_ << "\t eIsolateSoft_ " << eIsolate2_
316  << "\t hcalScale_ " << hcalScale_ << "\n\t momentumLow_ " << pTrackLow_ << "\t prescaleLow_ " << prescaleLow_
317  << "\t momentumHigh_ " << pTrackHigh_ << "\t prescaleHigh_ " << prescaleHigh_ << "\n\t useRaw_ " << useRaw_
318  << "\t ignoreTrigger_ " << ignoreTrigger_ << "\n\t useL1Trigegr_ " << useL1Trigger_ << "\t dataType_ "
319  << dataType_ << "\t unCorrect_ " << unCorrect_ << "\t collapseDepth_ " << collapseDepth_
320  << "\t L1TrigName_ " << l1TrigName_ << "\nThreshold for EB " << hitEthrEB_ << " EE " << hitEthrEE0_ << ":"
321  << hitEthrEE1_ << ":" << hitEthrEE2_ << ":" << hitEthrEE3_ << ":" << hitEthrEELo_ << ":" << hitEthrEEHi_;
322  edm::LogVerbatim("HcalIsoTrack") << "Process " << processName_ << " L1Filter:" << l1Filter_
323  << " L2Filter:" << l2Filter_ << " L3Filter:" << l3Filter_;
324  for (unsigned int k = 0; k < trigNames_.size(); ++k) {
325  edm::LogVerbatim("HcalIsoTrack") << "Trigger[" << k << "] " << trigNames_[k];
326  }
327 
328  for (int i = 0; i < 10; i++)
329  phibins_.push_back(-M_PI + 0.1 * (2 * i + 1) * M_PI);
330  for (int i = 0; i < 8; ++i)
331  etabins_.push_back(-2.1 + 0.6 * i);
332  etadist_ = etabins_[1] - etabins_[0];
333  phidist_ = phibins_[1] - phibins_[0];
334  etahalfdist_ = 0.5 * etadist_;
335  phihalfdist_ = 0.5 * phidist_;
336  edm::LogVerbatim("HcalIsoTrack") << "EtaDist " << etadist_ << " " << etahalfdist_ << " PhiDist " << phidist_ << " "
337  << phihalfdist_;
338  unsigned int k1(0), k2(0);
339  for (auto phi : phibins_) {
340  edm::LogVerbatim("HcalIsoTrack") << "phibin_[" << k1 << "] " << phi;
341  ++k1;
342  }
343  for (auto eta : etabins_) {
344  edm::LogVerbatim("HcalIsoTrack") << "etabin_[" << k2 << "] " << eta;
345  ++k2;
346  }
347 }
TFileService::kSharedResource
static const std::string kSharedResource
Definition: TFileService.h:76
HcalIsoTrkAnalyzer::labelHBHE_
const std::string labelHBHE_
Definition: HcalIsoTrkAnalyzer.cc:150
HcalIsoTrkAnalyzer::trigNames_
const std::vector< std::string > trigNames_
Definition: HcalIsoTrkAnalyzer.cc:131
HcalIsoTrkAnalyzer::selectionParameter_
spr::trackSelectionParameters selectionParameter_
Definition: HcalIsoTrkAnalyzer.cc:132
HLT_2018_cff.InputTag
InputTag
Definition: HLT_2018_cff.py:78913
HcalIsoTrkAnalyzer::l1TrigName_
const std::string l1TrigName_
Definition: HcalIsoTrkAnalyzer.cc:150
HcalIsoTrkAnalyzer::maxRestrictionP_
const double maxRestrictionP_
Definition: HcalIsoTrkAnalyzer.cc:138
mps_fire.i
i
Definition: mps_fire.py:341
HcalIsoTrkAnalyzer::tok_EE_
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: HcalIsoTrkAnalyzer.cc:163
HcalIsoTrkAnalyzer::labelRecVtx_
const std::string labelRecVtx_
Definition: HcalIsoTrkAnalyzer.cc:149
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:150
HcalIsoTrkAnalyzer::labelTower_
const std::string labelTower_
Definition: HcalIsoTrkAnalyzer.cc:150
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HcalIsoTrkAnalyzer::tok_alg_
edm::EDGetTokenT< BXVector< GlobalAlgBlk > > tok_alg_
Definition: HcalIsoTrkAnalyzer.cc:167
HcalIsoTrkAnalyzer::tok_EB_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: HcalIsoTrkAnalyzer.cc:162
HcalIsoTrkAnalyzer::tok_genTrack_
edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
Definition: HcalIsoTrkAnalyzer.cc:159
HcalIsoTrkAnalyzer::etabins_
std::vector< double > etabins_
Definition: HcalIsoTrkAnalyzer.cc:154
HcalIsoTrkAnalyzer::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: HcalIsoTrkAnalyzer.cc:164
HcalIsoTrkAnalyzer::tok_parts_
edm::EDGetTokenT< reco::GenParticleCollection > tok_parts_
Definition: HcalIsoTrkAnalyzer.cc:158
HcalIsoTrkAnalyzer::tok_trigEvt_
edm::EDGetTokenT< trigger::TriggerEvent > tok_trigEvt_
Definition: HcalIsoTrkAnalyzer.cc:156
HcalIsoTrkAnalyzer::labelGenTrack_
const std::string labelGenTrack_
Definition: HcalIsoTrkAnalyzer.cc:149
HcalIsoTrkAnalyzer::tok_bs_
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
Definition: HcalIsoTrkAnalyzer.cc:161
HcalIsoTrkAnalyzer::tok_ew_
edm::EDGetTokenT< GenEventInfoProduct > tok_ew_
Definition: HcalIsoTrkAnalyzer.cc:166
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:71
HcalIsoTrkAnalyzer::l1GtUtils_
l1t::L1TGlobalUtil * l1GtUtils_
Definition: HcalIsoTrkAnalyzer.cc:128
HcalIsoTrkAnalyzer::l3Filter_
const std::string l3Filter_
Definition: HcalIsoTrkAnalyzer.cc:135
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:50
HcalIsoTrkAnalyzer::slopeRestrictionP_
const double slopeRestrictionP_
Definition: HcalIsoTrkAnalyzer.cc:138
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
HcalIsoTrkAnalyzer::l2Filter_
const std::string l2Filter_
Definition: HcalIsoTrkAnalyzer.cc:135
HcalIsoTrkAnalyzer::processName_
const std::string processName_
Definition: HcalIsoTrkAnalyzer.cc:134
HcalIsoTrkAnalyzer::labelEB_
const std::string labelEB_
Definition: HcalIsoTrkAnalyzer.cc:149
HcalIsoTrkAnalyzer::phibins_
std::vector< double > phibins_
Definition: HcalIsoTrkAnalyzer.cc:154
HcalIsoTrkAnalyzer::theTriggerResultsLabel_
const edm::InputTag theTriggerResultsLabel_
Definition: HcalIsoTrkAnalyzer.cc:148
HcalIsoTrkAnalyzer::hdc_
const HcalDDDRecConstants * hdc_
Definition: HcalIsoTrkAnalyzer.cc:153
spr::trackSelectionParameters::minQuality
reco::TrackBase::TrackQuality minQuality
Definition: TrackSelection.h:21
HcalIsoTrkAnalyzer::tok_trigRes_
edm::EDGetTokenT< edm::TriggerResults > tok_trigRes_
Definition: HcalIsoTrkAnalyzer.cc:157
HcalIsoTrkAnalyzer::triggerEvent_
const edm::InputTag triggerEvent_
Definition: HcalIsoTrkAnalyzer.cc:148
HcalIsoTrkAnalyzer::tok_recVtx_
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
Definition: HcalIsoTrkAnalyzer.cc:160
HcalIsoTrkAnalyzer::tok_cala_
edm::EDGetTokenT< CaloTowerCollection > tok_cala_
Definition: HcalIsoTrkAnalyzer.cc:165
HcalIsoTrkAnalyzer::labelEE_
const std::string labelEE_
Definition: HcalIsoTrkAnalyzer.cc:150
HcalIsoTrkAnalyzer::l1Filter_
const std::string l1Filter_
Definition: HcalIsoTrkAnalyzer.cc:134
HcalIsoTrkAnalyzer::theTrackQuality_
const std::string theTrackQuality_
Definition: HcalIsoTrkAnalyzer.cc:133
HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer ( )
inlineoverride

Definition at line 90 of file HcalIsoTrkAnalyzer.cc.

References analyze(), beginJob(), beginRun(), dR(), endRun(), eThreshold(), fillDescriptions(), fillTree(), genParticles2HepMC_cfi::genParticles, photonIsolationHIProducer_cfi::hbhe, photons_cff::ids, rhoh(), storeEnergy(), and trackP().

90 {}

Member Function Documentation

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

Implements edm::one::EDAnalyzerBase.

Definition at line 349 of file HcalIsoTrkAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), Calorimetry_cff::bField, edm::EventBase::bunchCrossing(), dataType_, l1t::L1TGlobalUtil::decisionsFinal(), fileCollector::done, flavorHistoryFilter_cfi::dr, dR(), trigger::TriggerObject::energy(), trigger::TriggerObject::eta(), edm::EventID::event(), fillTree(), trigger::TriggerEvent::filterKeys(), trigger::TriggerEvent::filterTag(), spr::find(), genParticles2HepMC_cfi::genEventInfo, genParticles2HepMC_cfi::genParticles, edm::EventSetup::get(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), photonIsolationHIProducer_cfi::hbhe, ValidationMatrix::hlt, hltConfig_, mps_fire::i, trigger::TriggerObject::id(), edm::EventBase::id(), ignoreTrigger_, edm::HandleBase::isValid(), dqmdumpme::k, l1Filter_, l1GtUtils_, l1TrigName_, l2Filter_, l3Filter_, label, edm::InputTag::label(), labelEB_, labelEE_, labelGenTrack_, labelHBHE_, edm::EventBase::luminosityBlock(), trigger::TriggerObject::mass(), HLTConfigProvider::moduleLabels(), LogMessageMonitor_cff::modules, names, convertSQLiteXML::ok, trigger::TriggerObject::phi(), reco::BeamSpot::position(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), spr::propagateCALO(), trigger::TriggerObject::pt(), trigger::TriggerObject::px(), trigger::TriggerObject::py(), trigger::TriggerObject::pz(), l1t::L1TGlobalUtil::retrieveL1(), edm::EventID::run(), HcalCondObjectContainerBase::setTopo(), edm::HLTGlobalStatus::size(), trigger::TriggerEvent::sizeFilters(), 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, theTrackQuality_, tok_alg_, tok_bs_, tok_cala_, tok_EB_, tok_EE_, tok_ew_, tok_genTrack_, tok_hbhe_, tok_parts_, tok_recVtx_, tok_trigEvt_, tok_trigRes_, tree2, PDWG_DiPhoton_SD_cff::triggerEvent, triggerEvent_, HLTConfigProvider::triggerIndex(), L1TEGammaOffline_cfi::triggerNames, edm::TriggerNames::triggerNames(), edm::Event::triggerNames(), triggerResults, trigNames_, useL1Trigger_, GenEventInfoProduct::weight(), x, y, and z.

Referenced by ~HcalIsoTrkAnalyzer().

349  {
350  t_Run = iEvent.id().run();
351  t_Event = iEvent.id().event();
352  t_DataType = dataType_;
353 #ifdef EDM_ML_DEBUG
354  edm::LogVerbatim("HcalIsoTrack") << "Run " << t_Run << " Event " << t_Event << " type " << t_DataType
355  << " Luminosity " << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing();
356 #endif
357  //Get magnetic field and ECAL channel status
358  edm::ESHandle<MagneticField> bFieldH;
359  iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
360  const MagneticField* bField = bFieldH.product();
361 
362  edm::ESHandle<EcalChannelStatus> ecalChStatus;
363  iSetup.get<EcalChannelStatusRcd>().get(ecalChStatus);
364  const EcalChannelStatus* theEcalChStatus = ecalChStatus.product();
365 
366  edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
367  iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
368  const EcalSeverityLevelAlgo* theEcalSevlv = sevlv.product();
369 
370  // get handles to calogeometry and calotopology
371  edm::ESHandle<CaloGeometry> pG;
372  iSetup.get<CaloGeometryRecord>().get(pG);
373  const CaloGeometry* geo = pG.product();
374 
375  edm::ESHandle<CaloTopology> theCaloTopology;
376  iSetup.get<CaloTopologyRecord>().get(theCaloTopology);
377  const CaloTopology* caloTopology = theCaloTopology.product();
378 
379  edm::ESHandle<HcalTopology> htopo;
380  iSetup.get<HcalRecNumberingRecord>().get(htopo);
381  const HcalTopology* theHBHETopology = htopo.product();
382 
383  edm::ESHandle<HcalRespCorrs> resp;
384  iSetup.get<HcalRespCorrsRcd>().get(resp);
385  HcalRespCorrs* respCorrs = new HcalRespCorrs(*resp.product());
386  respCorrs->setTopo(theHBHETopology);
387 
388  //=== genParticle information
389  edm::Handle<reco::GenParticleCollection> genParticles;
390  iEvent.getByToken(tok_parts_, genParticles);
391 
392  bool okC(true);
393  //Get track collection
394  edm::Handle<reco::TrackCollection> trkCollection;
395  iEvent.getByToken(tok_genTrack_, trkCollection);
396  if (!trkCollection.isValid()) {
397  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelGenTrack_;
398  okC = false;
399  }
400 
401  //event weight for FLAT sample
402  t_EventWeight = 1.0;
403  edm::Handle<GenEventInfoProduct> genEventInfo;
404  iEvent.getByToken(tok_ew_, genEventInfo);
405  if (genEventInfo.isValid())
406  t_EventWeight = genEventInfo->weight();
407 
408  //Define the best vertex and the beamspot
409  edm::Handle<reco::VertexCollection> recVtxs;
410  iEvent.getByToken(tok_recVtx_, recVtxs);
411  edm::Handle<reco::BeamSpot> beamSpotH;
412  iEvent.getByToken(tok_bs_, beamSpotH);
413  math::XYZPoint leadPV(0, 0, 0);
414  t_goodPV = t_nVtx = 0;
415  if (recVtxs.isValid() && !(recVtxs->empty())) {
416  t_nVtx = recVtxs->size();
417  for (unsigned int k = 0; k < recVtxs->size(); ++k) {
418  if (!((*recVtxs)[k].isFake()) && ((*recVtxs)[k].ndof() > 4)) {
419  if (t_goodPV == 0)
420  leadPV = math::XYZPoint((*recVtxs)[k].x(), (*recVtxs)[k].y(), (*recVtxs)[k].z());
421  t_goodPV++;
422  }
423  }
424  }
425  if (t_goodPV == 0 && beamSpotH.isValid()) {
426  leadPV = beamSpotH->position();
427  }
428  t_allvertex = t_goodPV;
429 #ifdef EDM_ML_DEBUG
430  edm::LogVerbatim("HcalIsoTrack") << "Primary Vertex " << leadPV << " out of " << t_goodPV << " vertex";
431  if (beamSpotH.isValid()) {
432  edm::LogVerbatim("HcalIsoTrack") << " Beam Spot " << beamSpotH->position();
433  }
434 #endif
435  // RecHits
436  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
437  iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
438  if (!barrelRecHitsHandle.isValid()) {
439  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEB_;
440  okC = false;
441  }
442  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
443  iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
444  if (!endcapRecHitsHandle.isValid()) {
445  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEE_;
446  okC = false;
447  }
448  edm::Handle<HBHERecHitCollection> hbhe;
449  iEvent.getByToken(tok_hbhe_, hbhe);
450  if (!hbhe.isValid()) {
451  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelHBHE_;
452  okC = false;
453  }
454  edm::Handle<CaloTowerCollection> caloTower;
455  iEvent.getByToken(tok_cala_, caloTower);
456 
457  //Propagate tracks to calorimeter surface)
458  std::vector<spr::propagatedTrackDirection> trkCaloDirections;
459  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDirections, false);
460  std::vector<spr::propagatedTrackID> trkCaloDets;
461  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDets, false);
462  std::vector<math::XYZTLorentzVector> vecL1, vecL3;
463  t_RunNo = iEvent.id().run();
464  t_EventNo = iEvent.id().event();
465  t_Tracks = trkCollection->size();
466  t_TracksProp = trkCaloDirections.size();
467  t_ietaAll->clear();
468  t_ietaGood->clear();
469  t_trackType->clear();
470  t_trgbits->clear();
471  t_hltbits->clear();
472 #ifdef EDM_ML_DEBUG
473  edm::LogVerbatim("HcalIsoTrack") << "# of propagated tracks " << t_TracksProp << " out of " << t_Tracks
474  << " with Trigger " << ignoreTrigger_;
475 #endif
476 
477  //Trigger
478  t_trgbits->assign(trigNames_.size(), false);
479  t_hltbits->assign(trigNames_.size(), false);
480  t_TracksSaved = t_TracksLoose = t_TracksTight = 0;
481  t_L1Bit = true;
482  t_TrigPass = false;
483 
484  //L1
485  l1GtUtils_->retrieveL1(iEvent, iSetup, tok_alg_);
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
499  edm::Handle<edm::TriggerResults> triggerResults;
500  iEvent.getByToken(tok_trigRes_, triggerResults);
501  if (triggerResults.isValid()) {
502  const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
503  const std::vector<std::string>& names = triggerNames.triggerNames();
504  if (!trigNames_.empty()) {
505  for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
506  int hlt = triggerResults->accept(iHLT);
507  for (unsigned int i = 0; i < trigNames_.size(); ++i) {
508  if (names[iHLT].find(trigNames_[i]) != std::string::npos) {
509  t_trgbits->at(i) = (hlt > 0);
510  t_hltbits->at(i) = (hlt > 0);
511  if (hlt > 0)
512  t_TrigPass = true;
513 #ifdef EDM_ML_DEBUG
514  edm::LogVerbatim("HcalIsoTrack")
515  << "This trigger " << names[iHLT] << " Flag " << hlt << ":" << t_trgbits->at(i);
516 #endif
517  }
518  }
519  }
520  }
521  }
522 #ifdef EDM_ML_DEBUG
523  edm::LogVerbatim("HcalIsoTrack") << "HLT Information shows " << t_TrigPass << ":" << trigNames_.empty() << ":" << okC;
524 #endif
525 
526  std::array<int, 3> ntksave{{0, 0, 0}};
527  if (ignoreTrigger_ || useL1Trigger_) {
528  t_l1pt = t_l1eta = t_l1phi = 0;
529  t_l3pt = t_l3eta = t_l3phi = 0;
530  if (ignoreTrigger_ || t_L1Bit)
531  ntksave = fillTree(vecL1,
532  vecL3,
533  leadPV,
534  trkCaloDirections,
535  trkCaloDets,
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 {
551  trigger::TriggerEvent triggerEvent;
552  edm::Handle<trigger::TriggerEvent> triggerEventHandle;
553  iEvent.getByToken(tok_trigEvt_, triggerEventHandle);
554  if (!triggerEventHandle.isValid()) {
555  edm::LogWarning("HcalIsoTrack") << "Error! Can't get the product " << triggerEvent_.label();
556  } else if (okC) {
557  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  trkCaloDirections,
643  trkCaloDets,
644  geo,
645  caloTopology,
646  theHBHETopology,
647  theEcalChStatus,
648  theEcalSevlv,
649  barrelRecHitsHandle,
650  endcapRecHitsHandle,
651  hbhe,
652  caloTower,
653  genParticles,
654  respCorrs);
655  t_TracksSaved += ntksave[0];
656  t_TracksLoose += ntksave[1];
657  t_TracksTight += ntksave[2];
658  done = true;
659  }
660  }
661  }
662  }
663  }
664  }
665 #ifdef EDM_ML_DEBUG
666  edm::LogVerbatim("HcalIsoTrack") << "Final results on selected tracks " << t_TracksSaved << ":" << t_TracksLoose
667  << ":" << t_TracksTight;
668 #endif
669  t_TrigPassSel = (t_TracksSaved > 0);
670  tree2->Fill();
671 }
HcalIsoTrkAnalyzer::labelHBHE_
const std::string labelHBHE_
Definition: HcalIsoTrkAnalyzer.cc:150
HcalIsoTrkAnalyzer::trigNames_
const std::vector< std::string > trigNames_
Definition: HcalIsoTrkAnalyzer.cc:131
spr::propagateCALO
std::vector< spr::propagatedTrackID > propagateCALO(edm::Handle< reco::TrackCollection > &trkCollection, const CaloGeometry *geo, const MagneticField *bField, const std::string &theTrackQuality, bool debug=false)
Definition: CaloPropagateTrack.cc:141
HcalIsoTrkAnalyzer::l1TrigName_
const std::string l1TrigName_
Definition: HcalIsoTrkAnalyzer.cc:150
mps_fire.i
i
Definition: mps_fire.py:341
HcalIsoTrkAnalyzer::tok_EE_
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: HcalIsoTrkAnalyzer.cc:163
trigger::TriggerEvent
The single EDProduct to be saved for each event (AOD case)
Definition: TriggerEvent.h:25
trigger::TriggerEvent::sizeFilters
trigger::size_type sizeFilters() const
Definition: TriggerEvent.h:146
L1TEGammaOffline_cfi.triggerNames
triggerNames
Definition: L1TEGammaOffline_cfi.py:40
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HcalIsoTrkAnalyzer::t_ietaAll
std::vector< int > * t_ietaAll
Definition: HcalIsoTrkAnalyzer.cc:189
edm::HLTGlobalStatus::accept
bool accept() const
Has at least one path accepted the event?
Definition: HLTGlobalStatus.h:49
trigger::TriggerEvent::filterKeys
const Keys & filterKeys(trigger::size_type index) const
Definition: TriggerEvent.h:118
HcalIsoTrkAnalyzer::tok_alg_
edm::EDGetTokenT< BXVector< GlobalAlgBlk > > tok_alg_
Definition: HcalIsoTrkAnalyzer.cc:167
HcalIsoTrkAnalyzer::tok_EB_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: HcalIsoTrkAnalyzer.cc:162
HcalIsoTrkAnalyzer::tok_genTrack_
edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
Definition: HcalIsoTrkAnalyzer.cc:159
GenEventInfoProduct::weight
double weight() const
Definition: GenEventInfoProduct.h:33
HcalIsoTrkAnalyzer::fillTree
std::array< int, 3 > fillTree(std::vector< math::XYZTLorentzVector > &vecL1, std::vector< math::XYZTLorentzVector > &vecL3, math::XYZPoint &leadPV, std::vector< spr::propagatedTrackDirection > &trkCaloDirections, std::vector< spr::propagatedTrackID > &trkCaloDets, const CaloGeometry *geo, const CaloTopology *topo, const HcalTopology *theHBHETopology, const EcalChannelStatus *theEcalChStatus, const EcalSeverityLevelAlgo *theEcalSevlv, edm::Handle< EcalRecHitCollection > &barrelRecHitsHandle, edm::Handle< EcalRecHitCollection > &endcapRecHitsHandle, edm::Handle< HBHERecHitCollection > &hbhe, edm::Handle< CaloTowerCollection > &towerHandle, edm::Handle< reco::GenParticleCollection > &genParticles, const HcalRespCorrs *respCorrs)
Definition: HcalIsoTrkAnalyzer.cc:873
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
edm::TriggerNames::triggerNames
Strings const & triggerNames() const
Definition: TriggerNames.cc:20
names
const std::string names[nVars_]
Definition: PhotonIDValueMapProducer.cc:122
HcalIsoTrkAnalyzer::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: HcalIsoTrkAnalyzer.cc:164
Calorimetry_cff.bField
bField
Definition: Calorimetry_cff.py:292
HcalIsoTrkAnalyzer::tok_parts_
edm::EDGetTokenT< reco::GenParticleCollection > tok_parts_
Definition: HcalIsoTrkAnalyzer.cc:158
HLTConfigProvider::triggerIndex
unsigned int triggerIndex(const std::string &triggerName) const
slot position of trigger path in trigger table (0 to size-1)
Definition: HLTConfigProvider.h:72
trigger::TriggerObject
Single trigger physics object (e.g., an isolated muon)
Definition: TriggerObject.h:21
photonIsolationHIProducer_cfi.hbhe
hbhe
Definition: photonIsolationHIProducer_cfi.py:8
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
label
char const * label
Definition: PFTauDecayModeTools.cc:11
iEvent
int iEvent
Definition: GenABIO.cc:224
HcalIsoTrkAnalyzer::tok_trigEvt_
edm::EDGetTokenT< trigger::TriggerEvent > tok_trigEvt_
Definition: HcalIsoTrkAnalyzer.cc:156
HcalIsoTrkAnalyzer::labelGenTrack_
const std::string labelGenTrack_
Definition: HcalIsoTrkAnalyzer.cc:149
HcalIsoTrkAnalyzer::tok_bs_
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
Definition: HcalIsoTrkAnalyzer.cc:161
HcalIsoTrkAnalyzer::tok_ew_
edm::EDGetTokenT< GenEventInfoProduct > tok_ew_
Definition: HcalIsoTrkAnalyzer.cc:166
LogMessageMonitor_cff.modules
modules
Definition: LogMessageMonitor_cff.py:7
trigger::TriggerEvent::getObjects
const TriggerObjectCollection & getObjects() const
Definition: TriggerEvent.h:101
HcalIsoTrkAnalyzer::l1GtUtils_
l1t::L1TGlobalUtil * l1GtUtils_
Definition: HcalIsoTrkAnalyzer.cc:128
edm::HLTGlobalStatus::size
unsigned int size() const
Get number of paths stored.
Definition: HLTGlobalStatus.h:35
HcalIsoTrkAnalyzer::t_hltbits
std::vector< bool > * t_hltbits
Definition: HcalIsoTrkAnalyzer.cc:185
triggerResults
static std::string const triggerResults
Definition: EdmProvDump.cc:45
HcalIsoTrkAnalyzer::hltConfig_
HLTConfigProvider hltConfig_
Definition: HcalIsoTrkAnalyzer.cc:130
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
HcalIsoTrkAnalyzer::l3Filter_
const std::string l3Filter_
Definition: HcalIsoTrkAnalyzer.cc:135
HLTConfigProvider::moduleLabels
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
Definition: HLTConfigProvider.h:75
PDWG_DiPhoton_SD_cff.triggerEvent
triggerEvent
Definition: PDWG_DiPhoton_SD_cff.py:39
HcalIsoTrkAnalyzer::t_trgbits
std::vector< bool > * t_trgbits
Definition: HcalIsoTrkAnalyzer.cc:185
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:75
trigger::TriggerEvent::filterTag
const edm::InputTag filterTag(trigger::size_type index) const
Definition: TriggerEvent.h:108
HcalIsoTrkAnalyzer::l2Filter_
const std::string l2Filter_
Definition: HcalIsoTrkAnalyzer.cc:135
edm::Handle::product
T const * product() const
Definition: Handle.h:69
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:19
HcalIsoTrkAnalyzer::labelEB_
const std::string labelEB_
Definition: HcalIsoTrkAnalyzer.cc:149
HcalIsoTrkAnalyzer::dR
double dR(math::XYZTLorentzVector &, math::XYZTLorentzVector &)
Definition: HcalIsoTrkAnalyzer.cc:1229
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
genParticles2HepMC_cfi.genEventInfo
genEventInfo
Definition: genParticles2HepMC_cfi.py:6
l1t::L1TGlobalUtil::decisionsFinal
const std::vector< std::pair< std::string, bool > > & decisionsFinal()
Definition: L1TGlobalUtil.h:136
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:59
genParticles2HepMC_cfi.genParticles
genParticles
Definition: genParticles2HepMC_cfi.py:4
HcalIsoTrkAnalyzer::tok_trigRes_
edm::EDGetTokenT< edm::TriggerResults > tok_trigRes_
Definition: HcalIsoTrkAnalyzer.cc:157
HcalIsoTrkAnalyzer::triggerEvent_
const edm::InputTag triggerEvent_
Definition: HcalIsoTrkAnalyzer.cc:148
HcalIsoTrkAnalyzer::t_ietaGood
std::vector< int > * t_ietaGood
Definition: HcalIsoTrkAnalyzer.cc:189
HcalIsoTrkAnalyzer::tok_recVtx_
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
Definition: HcalIsoTrkAnalyzer.cc:160
HcalIsoTrkAnalyzer::tok_cala_
edm::EDGetTokenT< CaloTowerCollection > tok_cala_
Definition: HcalIsoTrkAnalyzer.cc:165
flavorHistoryFilter_cfi.dr
dr
Definition: flavorHistoryFilter_cfi.py:37
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
HcalIsoTrkAnalyzer::labelEE_
const std::string labelEE_
Definition: HcalIsoTrkAnalyzer.cc:150
HcalCondObjectContainerBase::setTopo
void setTopo(const HcalTopology *topo)
Definition: HcalCondObjectContainerBase.cc:17
l1t::L1TGlobalUtil::retrieveL1
void retrieveL1(const edm::Event &iEvent, const edm::EventSetup &evSetup)
initialize the class (mainly reserve)
Definition: L1TGlobalUtil.cc:68
HcalIsoTrkAnalyzer::l1Filter_
const std::string l1Filter_
Definition: HcalIsoTrkAnalyzer.cc:134
HcalIsoTrkAnalyzer::theTrackQuality_
const std::string theTrackQuality_
Definition: HcalIsoTrkAnalyzer.cc:133
HcalIsoTrkAnalyzer::t_trackType
std::vector< int > * t_trackType
Definition: HcalIsoTrkAnalyzer.cc:189
void HcalIsoTrkAnalyzer::beginJob ( void  )
overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 673 of file HcalIsoTrkAnalyzer.cc.

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

Referenced by ~HcalIsoTrkAnalyzer().

673  {
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 }
HcalIsoTrkAnalyzer::t_ietaAll
std::vector< int > * t_ietaAll
Definition: HcalIsoTrkAnalyzer.cc:189
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
HcalIsoTrkAnalyzer::t_HitEnergies
std::vector< double > * t_HitEnergies
Definition: HcalIsoTrkAnalyzer.cc:184
HcalIsoTrkAnalyzer::t_DetIds
std::vector< unsigned int > * t_DetIds
Definition: HcalIsoTrkAnalyzer.cc:183
HcalIsoTrkAnalyzer::t_HitEnergies1
std::vector< double > * t_HitEnergies1
Definition: HcalIsoTrkAnalyzer.cc:184
HcalIsoTrkAnalyzer::t_HitEnergies3
std::vector< double > * t_HitEnergies3
Definition: HcalIsoTrkAnalyzer.cc:184
HcalIsoTrkAnalyzer::fs
edm::Service< TFileService > fs
Definition: HcalIsoTrkAnalyzer.cc:129
HcalIsoTrkAnalyzer::t_hltbits
std::vector< bool > * t_hltbits
Definition: HcalIsoTrkAnalyzer.cc:185
HcalIsoTrkAnalyzer::t_trgbits
std::vector< bool > * t_trgbits
Definition: HcalIsoTrkAnalyzer.cc:185
HcalIsoTrkAnalyzer::t_DetIds1
std::vector< unsigned int > * t_DetIds1
Definition: HcalIsoTrkAnalyzer.cc:183
tree
Definition: tree.py:1
HcalIsoTrkAnalyzer::t_ietaGood
std::vector< int > * t_ietaGood
Definition: HcalIsoTrkAnalyzer.cc:189
HcalIsoTrkAnalyzer::t_DetIds3
std::vector< unsigned int > * t_DetIds3
Definition: HcalIsoTrkAnalyzer.cc:183
HcalIsoTrkAnalyzer::t_trackType
std::vector< int > * t_trackType
Definition: HcalIsoTrkAnalyzer.cc:189
void HcalIsoTrkAnalyzer::beginRun ( edm::Run const &  iRun,
edm::EventSetup const &  iSetup 
)
overrideprivate

Definition at line 754 of file HcalIsoTrkAnalyzer.cc.

References RemoveAddSevLevel::flag, edm::EventSetup::get(), hdc_, hltConfig_, HLTConfigProvider::init(), dqmiodumpmetadata::n, nRun_, processName_, edm::ESHandle< T >::product(), edm::RunBase::run(), HLTConfigProvider::size(), HLTConfigProvider::triggerIndex(), and trigNames_.

Referenced by ~HcalIsoTrkAnalyzer().

754  {
755  edm::ESHandle<HcalDDDRecConstants> pHRNDC;
756  iSetup.get<HcalRecNumberingRecord>().get(pHRNDC);
757  hdc_ = pHRNDC.product();
758 
759  bool changed_(true);
760  bool flag = hltConfig_.init(iRun, iSetup, processName_, changed_);
761  edm::LogVerbatim("HcalIsoTrack") << "Run[" << nRun_ << "] " << iRun.run() << " process " << processName_
762  << " init flag " << flag << " change flag " << changed_;
763  // check if trigger names in (new) config
764  if (changed_) {
765  changed_ = false;
766 #ifdef EDM_ML_DEBUG
767  edm::LogVerbatim("HcalIsoTrack") << "New trigger menu found !!!";
768 #endif
769  const unsigned int n(hltConfig_.size());
770  for (unsigned itrig = 0; itrig < trigNames_.size(); itrig++) {
771  unsigned int triggerindx = hltConfig_.triggerIndex(trigNames_[itrig]);
772  if (triggerindx >= n) {
773  edm::LogWarning("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " does not exist in "
774  << "the current menu";
775 #ifdef EDM_ML_DEBUG
776  } else {
777  edm::LogVerbatim("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " exists";
778 #endif
779  }
780  }
781  }
782 }
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:62
HcalIsoTrkAnalyzer::trigNames_
const std::vector< std::string > trigNames_
Definition: HcalIsoTrkAnalyzer.cc:131
HLTConfigProvider::triggerIndex
unsigned int triggerIndex(const std::string &triggerName) const
slot position of trigger path in trigger table (0 to size-1)
Definition: HLTConfigProvider.h:72
HcalIsoTrkAnalyzer::hltConfig_
HLTConfigProvider hltConfig_
Definition: HcalIsoTrkAnalyzer.cc:130
HcalIsoTrkAnalyzer::processName_
const std::string processName_
Definition: HcalIsoTrkAnalyzer.cc:134
HLTConfigProvider::init
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d&#39;tor
Definition: HLTConfigProvider.cc:36
HcalIsoTrkAnalyzer::hdc_
const HcalDDDRecConstants * hdc_
Definition: HcalIsoTrkAnalyzer.cc:153
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
double HcalIsoTrkAnalyzer::dR ( math::XYZTLorentzVector &  vec1,
math::XYZTLorentzVector &  vec2 
)
private

Definition at line 1229 of file HcalIsoTrkAnalyzer.cc.

References reco::deltaR().

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

1229  {
1230  return reco::deltaR(vec1.eta(), vec1.phi(), vec2.eta(), vec2.phi());
1231 }
vec1
std::vector< double > vec1
Definition: HCALResponse.h:15
reco::deltaR
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
vec2
std::vector< vec1 > vec2
Definition: HCALResponse.h:16
void HcalIsoTrkAnalyzer::endRun ( edm::Run const &  iRun,
edm::EventSetup const &   
)
overrideprivate

Definition at line 785 of file HcalIsoTrkAnalyzer.cc.

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

Referenced by ~HcalIsoTrkAnalyzer().

785  {
786  nRun_++;
787  edm::LogVerbatim("HcalIsoTrack") << "endRun[" << nRun_ << "] " << iRun.run();
788 }
double HcalIsoTrkAnalyzer::eThreshold ( const DetId &  id,
const CaloGeometry *  geo 
) const
private

Definition at line 1280 of file HcalIsoTrkAnalyzer.cc.

References funct::abs(), EcalBarrel, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), CaloGeometry::getPosition(), hitEthrEB_, hitEthrEE0_, hitEthrEE1_, hitEthrEE2_, hitEthrEE3_, hitEthrEEHi_, and hitEthrEELo_.

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

1280  {
1281  const GlobalPoint& pos = geo->getPosition(id);
1282  double eta = std::abs(pos.eta());
1283  double eThr(hitEthrEB_);
1284  if (id.subdetId() != EcalBarrel) {
1285  eThr = (((eta * hitEthrEE3_ + hitEthrEE2_) * eta + hitEthrEE1_) * eta + hitEthrEE0_);
1286  if (eThr < hitEthrEELo_)
1287  eThr = hitEthrEELo_;
1288  else if (eThr > hitEthrEEHi_)
1289  eThr = hitEthrEEHi_;
1290  }
1291  return eThr;
1292 }
CaloGeometry::getPosition
GlobalPoint getPosition(const DetId &id) const
Get the position of a given detector id.
Definition: CaloGeometry.cc:50
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:73
void HcalIsoTrkAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 790 of file HcalIsoTrkAnalyzer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), HLT_2018_cff::InputTag, AlCaHLTBitMon_QueryRunRegistry::string, and ClusterTask_cfi::trig.

Referenced by ~HcalIsoTrkAnalyzer().

790  {
791  edm::ParameterSetDescription desc;
792  std::vector<std::string> trig = {"HLT_PFJet40",
793  "HLT_PFJet60",
794  "HLT_PFJet80",
795  "HLT_PFJet140",
796  "HLT_PFJet200",
797  "HLT_PFJet260",
798  "HLT_PFJet320",
799  "HLT_PFJet400",
800  "HLT_PFJet450",
801  "HLT_PFJet500"};
802  desc.add<std::vector<std::string>>("triggers", trig);
803  desc.add<std::string>("processName", "HLT");
804  desc.add<std::string>("l1Filter", "");
805  desc.add<std::string>("l2Filter", "L2Filter");
806  desc.add<std::string>("l3Filter", "Filter");
807  // following 10 parameters are parameters to select good tracks
808  desc.add<std::string>("trackQuality", "highPurity");
809  desc.add<double>("minTrackPt", 1.0);
810  desc.add<double>("maxDxyPV", 0.02);
811  desc.add<double>("maxDzPV", 0.02);
812  desc.add<double>("maxChi2", 5.0);
813  desc.add<double>("maxDpOverP", 0.1);
814  desc.add<int>("minOuterHit", 4);
815  desc.add<int>("minLayerCrossed", 8);
816  desc.add<int>("maxInMiss", 0);
817  desc.add<int>("maxOutMiss", 0);
818  // Minimum momentum of selected isolated track and signal zone
819  desc.add<double>("minimumTrackP", 20.0);
820  desc.add<double>("coneRadius", 34.98);
821  // signal zone in ECAL and MIP energy cutoff
822  desc.add<double>("coneRadiusMIP", 14.0);
823  desc.add<double>("coneRadiusMIP2", 18.0);
824  desc.add<double>("coneRadiusMIP3", 20.0);
825  desc.add<double>("coneRadiusMIP4", 22.0);
826  desc.add<double>("coneRadiusMIP5", 24.0);
827  desc.add<double>("maximumEcalEnergy", 2.0);
828  // following 4 parameters are for isolation cuts and described in the code
829  desc.add<double>("maxTrackP", 8.0);
830  desc.add<double>("slopeTrackP", 0.05090504066);
831  desc.add<double>("isolationEnergyTight", 2.0);
832  desc.add<double>("isolationEnergyLoose", 10.0);
833  // energy thershold for ECAL (from Egamma group)
834  desc.add<double>("EBHitEnergyThreshold", 0.08);
835  desc.add<double>("EEHitEnergyThreshold0", 0.30);
836  desc.add<double>("EEHitEnergyThreshold1", 0.00);
837  desc.add<double>("EEHitEnergyThreshold2", 0.00);
838  desc.add<double>("EEHitEnergyThreshold3", 0.00);
839  desc.add<double>("EEHitEnergyThresholdLow", 0.30);
840  desc.add<double>("EEHitEnergyThresholdHigh", 0.30);
841  // prescale factors
842  desc.add<double>("momentumLow", 40.0);
843  desc.add<double>("momentumHigh", 60.0);
844  desc.add<int>("prescaleLow", 1);
845  desc.add<int>("prescaleHigh", 1);
846  // various labels for collections used in the code
847  desc.add<edm::InputTag>("labelTriggerEvent", edm::InputTag("hltTriggerSummaryAOD", "", "HLT"));
848  desc.add<edm::InputTag>("labelTriggerResult", edm::InputTag("TriggerResults", "", "HLT"));
849  desc.add<std::string>("labelTrack", "generalTracks");
850  desc.add<std::string>("labelVertex", "offlinePrimaryVertices");
851  desc.add<std::string>("labelEBRecHit", "EcalRecHitsEB");
852  desc.add<std::string>("labelEERecHit", "EcalRecHitsEE");
853  desc.add<std::string>("labelHBHERecHit", "hbhereco");
854  desc.add<std::string>("labelBeamSpot", "offlineBeamSpot");
855  desc.add<std::string>("labelCaloTower", "towerMaker");
856  desc.add<edm::InputTag>("algInputTag", edm::InputTag("gtStage2Digis"));
857  desc.add<edm::InputTag>("extInputTag", edm::InputTag("gtStage2Digis"));
858  desc.addUntracked<std::string>("moduleName", "");
859  desc.addUntracked<std::string>("producerName", "");
860  // Various flags used for selecting tracks, choice of energy Method2/0
861  // Data type 0/1 for single jet trigger or others
862  desc.addUntracked<int>("useRaw", 0);
863  desc.addUntracked<bool>("ignoreTriggers", false);
864  desc.addUntracked<bool>("useL1Trigger", false);
865  desc.addUntracked<double>("hcalScale", 1.0);
866  desc.addUntracked<int>("dataType", 0);
867  desc.addUntracked<bool>("unCorrect", false);
868  desc.addUntracked<bool>("collapseDepth", false);
869  desc.addUntracked<std::string>("l1TrigName", "L1_SingleJet60");
870  descriptions.add("HcalIsoTrkAnalyzer", desc);
871 }
HLT_2018_cff.InputTag
InputTag
Definition: HLT_2018_cff.py:78913
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
ClusterTask_cfi.trig
trig
Definition: ClusterTask_cfi.py:336
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
std::array< int, 3 > HcalIsoTrkAnalyzer::fillTree ( std::vector< math::XYZTLorentzVector > &  vecL1,
std::vector< math::XYZTLorentzVector > &  vecL3,
math::XYZPoint &  leadPV,
std::vector< spr::propagatedTrackDirection > &  trkCaloDirections,
std::vector< spr::propagatedTrackID > &  trkCaloDets,
const CaloGeometry *  geo,
const CaloTopology *  topo,
const HcalTopology *  theHBHETopology,
const EcalChannelStatus *  theEcalChStatus,
const EcalSeverityLevelAlgo *  theEcalSevlv,
edm::Handle< EcalRecHitCollection > &  barrelRecHitsHandle,
edm::Handle< EcalRecHitCollection > &  endcapRecHitsHandle,
edm::Handle< HBHERecHitCollection > &  hbhe,
edm::Handle< CaloTowerCollection > &  towerHandle,
edm::Handle< reco::GenParticleCollection > &  genParticles,
const HcalRespCorrs *  respCorrs 
)
private

Definition at line 873 of file HcalIsoTrkAnalyzer.cc.

References a_charIsoR_, a_coneR1_, a_coneR2_, a_coneR_, a_mipR2_, a_mipR3_, a_mipR4_, a_mipR5_, a_mipR_, funct::abs(), accept(), spr::chargeIsolationCone(), spr::chargeIsolationEcal(), flavorHistoryFilter_cfi::dr, dR(), spr::eCone_ecal(), spr::eCone_hcal(), spr::eECALmatrix(), eEcalMax_, spr::eHCALmatrix(), eIsolate1_, eIsolate2_, reco::TrackBase::eta(), eThreshold(), JetChargeProducer_cfi::exp, spr::goodTrack(), testProducerWithPsetDescEmpty_cfi::i3, photons_cff::ids, HcalDetId::ieta(), HcalDetId::iphi(), edm::HandleBase::isValid(), dqmdumpme::k, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, maxRestrictionP_, nHigh_, jetSelector_cfi::nLoose, nLow_, nTracks(), reco::TrackBase::p(), reco::TrackBase::phi(), prescaleHigh_, prescaleLow_, reco::TrackBase::pt(), pTrackHigh_, pTrackLow_, pTrackMin_, reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), rhoh(), selectionParameter_, slopeRestrictionP_, storeEnergy(), t_DetIds, t_DetIds1, t_DetIds3, t_eAnnular, t_eHcal, t_eHcal10, t_eHcal30, t_emaxNearP, t_eMipDR, t_eMipDR2, t_eMipDR3, t_eMipDR4, t_eMipDR5, t_gentrackP, t_hAnnular, t_HitEnergies, t_HitEnergies1, t_HitEnergies3, t_hmaxNearP, t_ieta, t_ietaAll, t_ietaGood, t_iphi, t_mindR1, t_mindR2, t_nTrk, t_p, t_phi, t_pt, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_rhoh, t_selectTk, t_trackType, t_trgbits, trackP(), and useRaw_.

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

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

Definition at line 1249 of file HcalIsoTrkAnalyzer.cc.

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

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

1249  {
1250  std::vector<double> sumPFNallSMDQH2;
1251  sumPFNallSMDQH2.reserve(phibins_.size() * etabins_.size());
1252 
1253  for (auto eta : etabins_) {
1254  for (auto phi : phibins_) {
1255  double hadder = 0;
1256  for (const auto& pf_it : (*tower)) {
1257  if (fabs(eta - pf_it.eta()) > etahalfdist_)
1258  continue;
1259  if (fabs(reco::deltaPhi(phi, pf_it.phi())) > phihalfdist_)
1260  continue;
1261  hadder += pf_it.hadEt();
1262  }
1263  sumPFNallSMDQH2.emplace_back(hadder);
1264  }
1265  }
1266 
1267  double evt_smdq(0);
1268  std::sort(sumPFNallSMDQH2.begin(), sumPFNallSMDQH2.end());
1269  if (sumPFNallSMDQH2.size() % 2)
1270  evt_smdq = sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 1) / 2];
1271  else
1272  evt_smdq = (sumPFNallSMDQH2[sumPFNallSMDQH2.size() / 2] + sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 2) / 2]) / 2.;
1273  double rhoh = evt_smdq / (etadist_ * phidist_);
1274 #ifdef EDM_ML_DEBUG
1275  edm::LogVerbatim("HcalIsoTrack") << "Rho " << evt_smdq << ":" << rhoh;
1276 #endif
1277  return rhoh;
1278 }
reco::deltaPhi
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:26
HcalIsoTrkAnalyzer::etabins_
std::vector< double > etabins_
Definition: HcalIsoTrkAnalyzer.cc:154
HcalIsoTrkAnalyzer::phibins_
std::vector< double > phibins_
Definition: HcalIsoTrkAnalyzer.cc:154
HcalIsoTrkAnalyzer::rhoh
double rhoh(const edm::Handle< CaloTowerCollection > &)
Definition: HcalIsoTrkAnalyzer.cc:1249
void HcalIsoTrkAnalyzer::storeEnergy ( int  indx,
const HcalRespCorrs *  respCorrs,
const std::vector< DetId > &  ids,
std::vector< double > &  edet,
double &  eHcal,
std::vector< unsigned int > *  detIds,
std::vector< double > *  hitEnergies 
)
private

Definition at line 1294 of file HcalIsoTrkAnalyzer.cc.

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

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

1300  {
1301  double ehcal(0);
1302  if (unCorrect_) {
1303  for (unsigned int k = 0; k < ids.size(); ++k) {
1304  double corr = (respCorrs->getValues(ids[k]))->getValue();
1305  if (corr != 0)
1306  edet[k] /= corr;
1307  ehcal += edet[k];
1308  }
1309  } else {
1310  for (const auto& en : edet)
1311  ehcal += en;
1312  }
1313  if (std::abs(ehcal - eHcal) > 0.001)
1314  edm::LogWarning("HcalIsoTrack") << "Check inconsistent energies: " << indx << " " << eHcal << ":" << ehcal
1315  << " from " << ids.size() << " cells";
1316  eHcal = hcalScale_ * ehcal;
1317 
1318  if (collapseDepth_) {
1319  std::map<HcalDetId, double> hitMap;
1320  for (unsigned int k = 0; k < ids.size(); ++k) {
1321  HcalDetId id = hdc_->mergedDepthDetId(HcalDetId(ids[k]));
1322  auto itr = hitMap.find(id);
1323  if (itr == hitMap.end()) {
1324  hitMap[id] = edet[k];
1325  } else {
1326  (itr->second) += edet[k];
1327  }
1328  }
1329  detIds->reserve(hitMap.size());
1330  hitEnergies->reserve(hitMap.size());
1331  for (const auto& hit : hitMap) {
1332  detIds->emplace_back(hit.first.rawId());
1333  hitEnergies->emplace_back(hit.second);
1334  }
1335  } else {
1336  detIds->reserve(ids.size());
1337  hitEnergies->reserve(ids.size());
1338  for (unsigned int k = 0; k < ids.size(); ++k) {
1339  detIds->emplace_back(ids[k].rawId());
1340  hitEnergies->emplace_back(edet[k]);
1341  }
1342  }
1343 #ifdef EDM_ML_DEBUG
1344  edm::LogVerbatim("HcalIsoTrack") << "Input to storeEnergy with " << ids.size() << " cells";
1345  for (unsigned int k = 0; k < ids.size(); ++k)
1346  edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId(ids[k]) << " E " << edet[k];
1347  edm::LogVerbatim("HcalIsoTrack") << "Output of storeEnergy with " << detIds->size() << " cells and Etot " << eHcal;
1348  for (unsigned int k = 0; k < detIds->size(); ++k)
1349  edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId((*detIds)[k]) << " E " << (*hitEnergies)[k];
1350 #endif
1351 }
HcalDDDRecConstants::mergedDepthDetId
HcalDetId mergedDepthDetId(const HcalDetId &id) const
Definition: HcalDDDRecConstants.cc:647
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:157
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:31
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
photons_cff.ids
ids
Definition: photons_cff.py:50
corr
JetCorrectorParameters corr
Definition: classes.h:5
HcalIsoTrkAnalyzer::hdc_
const HcalDDDRecConstants * hdc_
Definition: HcalIsoTrkAnalyzer.cc:153
reco::JetExtendedAssociation::getValue
const JetExtendedData & getValue(const Container &, const reco::JetBaseRef &)
get value for the association. Throw exception if no association found
Definition: JetExtendedAssociation.cc:50
double HcalIsoTrkAnalyzer::trackP ( const reco::Track *  pTrack,
const edm::Handle< reco::GenParticleCollection > &  genParticles 
)
private

Definition at line 1233 of file HcalIsoTrkAnalyzer.cc.

References reco::deltaR(), dR(), reco::TrackBase::eta(), edm::HandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, and reco::TrackBase::phi().

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

1234  {
1235  double pmom = -1.0;
1236  if (genParticles.isValid()) {
1237  double mindR(999.9);
1238  for (const auto& p : (*genParticles)) {
1239  double dR = reco::deltaR(pTrack->eta(), pTrack->phi(), p.momentum().Eta(), p.momentum().Phi());
1240  if (dR < mindR) {
1241  mindR = dR;
1242  pmom = p.momentum().R();
1243  }
1244  }
1245  }
1246  return pmom;
1247 }
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:614
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:617
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
reco::deltaR
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
HcalIsoTrkAnalyzer::dR
double dR(math::XYZTLorentzVector &, math::XYZTLorentzVector &)
Definition: HcalIsoTrkAnalyzer.cc:1229

Member Data Documentation

double HcalIsoTrkAnalyzer::a_charIsoR_
private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::a_coneR1_
private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::a_coneR2_
private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_coneR_
private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_mipR2_
private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_mipR3_
private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_mipR4_
private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_mipR5_
private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_mipR_
private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const bool HcalIsoTrkAnalyzer::collapseDepth_
private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

const int HcalIsoTrkAnalyzer::dataType_
private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::eEcalMax_
private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::eIsolate1_
private

Definition at line 139 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::eIsolate2_
private

Definition at line 139 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

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

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

double HcalIsoTrkAnalyzer::etadist_
private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

double HcalIsoTrkAnalyzer::etahalfdist_
private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

edm::Service<TFileService> HcalIsoTrkAnalyzer::fs
private

Definition at line 129 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob().

const double HcalIsoTrkAnalyzer::hcalScale_
private

Definition at line 139 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

const HcalDDDRecConstants* HcalIsoTrkAnalyzer::hdc_
private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and storeEnergy().

const double HcalIsoTrkAnalyzer::hitEthrEB_
private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEE0_
private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEE1_
private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEE2_
private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEE3_
private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEEHi_
private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEELo_
private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

HLTConfigProvider HcalIsoTrkAnalyzer::hltConfig_
private

Definition at line 130 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginRun().

const bool HcalIsoTrkAnalyzer::ignoreTrigger_
private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::l1Filter_
private

Definition at line 134 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

l1t::L1TGlobalUtil* HcalIsoTrkAnalyzer::l1GtUtils_
private

Definition at line 128 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::l1TrigName_
private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::l2Filter_
private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::l3Filter_
private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelEB_
private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelEE_
private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelGenTrack_
private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelHBHE_
private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelRecVtx_
private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelTower_
private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::maxRestrictionP_
private

Definition at line 138 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

unsigned int HcalIsoTrkAnalyzer::nHigh_
private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree().

unsigned int HcalIsoTrkAnalyzer::nLow_
private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree().

unsigned int HcalIsoTrkAnalyzer::nRun_
private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and endRun().

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

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

double HcalIsoTrkAnalyzer::phidist_
private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

double HcalIsoTrkAnalyzer::phihalfdist_
private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

const int HcalIsoTrkAnalyzer::prescaleHigh_
private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const int HcalIsoTrkAnalyzer::prescaleLow_
private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::processName_
private

Definition at line 134 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::pTrackHigh_
private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::pTrackLow_
private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::pTrackMin_
private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

spr::trackSelectionParameters HcalIsoTrkAnalyzer::selectionParameter_
private

Definition at line 132 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::slopeRestrictionP_
private

Definition at line 138 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

int HcalIsoTrkAnalyzer::t_allvertex
private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_DataType
private

Definition at line 171 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

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

Definition at line 183 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

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

Definition at line 183 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

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

Definition at line 183 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eAnnular
private

Definition at line 179 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eHcal
private

Definition at line 180 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eHcal10
private

Definition at line 180 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eHcal30
private

Definition at line 180 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_emaxNearP
private

Definition at line 179 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR
private

Definition at line 177 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR2
private

Definition at line 177 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR3
private

Definition at line 177 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR4
private

Definition at line 177 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR5
private

Definition at line 178 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_Event
private

Definition at line 171 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

unsigned int HcalIsoTrkAnalyzer::t_EventNo
private

Definition at line 170 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_EventWeight
private

Definition at line 173 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_gentrackP
private

Definition at line 178 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_goodPV
private

Definition at line 172 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_hAnnular
private

Definition at line 179 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

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

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

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

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

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

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

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

Definition at line 185 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_hmaxNearP
private

Definition at line 178 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_ieta
private

Definition at line 171 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

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

Definition at line 189 of file HcalIsoTrkAnalyzer.cc.

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

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

Definition at line 189 of file HcalIsoTrkAnalyzer.cc.

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

int HcalIsoTrkAnalyzer::t_iphi
private

Definition at line 171 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_L1Bit
private

Definition at line 186 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1eta
private

Definition at line 174 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1phi
private

Definition at line 174 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1pt
private

Definition at line 174 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3eta
private

Definition at line 175 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3phi
private

Definition at line 175 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3pt
private

Definition at line 175 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_mindR1
private

Definition at line 176 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_mindR2
private

Definition at line 176 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_nTrk
private

Definition at line 172 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_nVtx
private

Definition at line 172 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_p
private

Definition at line 173 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_phi
private

Definition at line 173 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_pt
private

Definition at line 173 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_qltyFlag
private

Definition at line 181 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_qltyMissFlag
private

Definition at line 181 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_qltyPVFlag
private

Definition at line 182 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_rhoh
private

Definition at line 180 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_Run
private

Definition at line 171 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

unsigned int HcalIsoTrkAnalyzer::t_RunNo
private

Definition at line 170 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

bool HcalIsoTrkAnalyzer::t_selectTk
private

Definition at line 181 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_Tracks
private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_TracksLoose
private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_TracksProp
private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_TracksSaved
private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_TracksTight
private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

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

Definition at line 189 of file HcalIsoTrkAnalyzer.cc.

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

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

Definition at line 185 of file HcalIsoTrkAnalyzer.cc.

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

bool HcalIsoTrkAnalyzer::t_TrigPass
private

Definition at line 182 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

bool HcalIsoTrkAnalyzer::t_TrigPassSel
private

Definition at line 182 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

const std::string HcalIsoTrkAnalyzer::theTrackQuality_
private

Definition at line 133 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const edm::InputTag HcalIsoTrkAnalyzer::theTriggerResultsLabel_
private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

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

Definition at line 167 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 161 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 162 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 163 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 166 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 164 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 160 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

TTree* HcalIsoTrkAnalyzer::tree
private

Definition at line 169 of file HcalIsoTrkAnalyzer.cc.

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

TTree * HcalIsoTrkAnalyzer::tree2
private

Definition at line 169 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

const edm::InputTag HcalIsoTrkAnalyzer::triggerEvent_
private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

Definition at line 131 of file HcalIsoTrkAnalyzer.cc.

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

const bool HcalIsoTrkAnalyzer::unCorrect_
private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

const bool HcalIsoTrkAnalyzer::useL1Trigger_
private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const int HcalIsoTrkAnalyzer::useRaw_
private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

Generated for CMSSW Reference Manual by doxygen 1.8.11