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
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase &&)=default
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESProxyIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

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

Private Member Functions

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

Private Attributes

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

Additional Inherited Members

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

Detailed Description

Definition at line 91 of file HcalIsoTrkAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalIsoTrkAnalyzer()

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

Definition at line 212 of file HcalIsoTrkAnalyzer.cc.

213  : trigNames_(iConfig.getParameter<std::vector<std::string> >("triggers")),
214  theTrackQuality_(iConfig.getParameter<std::string>("trackQuality")),
215  processName_(iConfig.getParameter<std::string>("processName")),
216  l1Filter_(iConfig.getParameter<std::string>("l1Filter")),
217  l2Filter_(iConfig.getParameter<std::string>("l2Filter")),
218  l3Filter_(iConfig.getParameter<std::string>("l3Filter")),
219  a_coneR_(iConfig.getParameter<double>("coneRadius")),
220  a_mipR_(iConfig.getParameter<double>("coneRadiusMIP")),
221  a_mipR2_(iConfig.getParameter<double>("coneRadiusMIP2")),
222  a_mipR3_(iConfig.getParameter<double>("coneRadiusMIP3")),
223  a_mipR4_(iConfig.getParameter<double>("coneRadiusMIP4")),
224  a_mipR5_(iConfig.getParameter<double>("coneRadiusMIP5")),
225  pTrackMin_(iConfig.getParameter<double>("minimumTrackP")),
226  eEcalMax_(iConfig.getParameter<double>("maximumEcalEnergy")),
227  maxRestrictionP_(iConfig.getParameter<double>("maxTrackP")),
228  slopeRestrictionP_(iConfig.getParameter<double>("slopeTrackP")),
229  hcalScale_(iConfig.getUntrackedParameter<double>("hHcalScale", 1.0)),
230  eIsolate1_(iConfig.getParameter<double>("isolationEnergyTight")),
231  eIsolate2_(iConfig.getParameter<double>("isolationEnergyLoose")),
232  pTrackLow_(iConfig.getParameter<double>("momentumLow")),
233  pTrackHigh_(iConfig.getParameter<double>("momentumHigh")),
234  prescaleLow_(iConfig.getParameter<int>("prescaleLow")),
235  prescaleHigh_(iConfig.getParameter<int>("prescaleHigh")),
236  useRaw_(iConfig.getUntrackedParameter<int>("useRaw", 0)),
237  dataType_(iConfig.getUntrackedParameter<int>("dataType", 0)),
238  mode_(iConfig.getUntrackedParameter<int>("outMode", 11)),
239  ignoreTrigger_(iConfig.getUntrackedParameter<bool>("ignoreTriggers", false)),
240  useL1Trigger_(iConfig.getUntrackedParameter<bool>("useL1Trigger", false)),
241  unCorrect_(iConfig.getUntrackedParameter<bool>("unCorrect", false)),
242  collapseDepth_(iConfig.getUntrackedParameter<bool>("collapseDepth", false)),
243  hitEthrEB_(iConfig.getParameter<double>("EBHitEnergyThreshold")),
244  hitEthrEE0_(iConfig.getParameter<double>("EEHitEnergyThreshold0")),
245  hitEthrEE1_(iConfig.getParameter<double>("EEHitEnergyThreshold1")),
246  hitEthrEE2_(iConfig.getParameter<double>("EEHitEnergyThreshold2")),
247  hitEthrEE3_(iConfig.getParameter<double>("EEHitEnergyThreshold3")),
248  hitEthrEELo_(iConfig.getParameter<double>("EEHitEnergyThresholdLow")),
249  hitEthrEEHi_(iConfig.getParameter<double>("EEHitEnergyThresholdHigh")),
250  triggerEvent_(iConfig.getParameter<edm::InputTag>("labelTriggerEvent")),
251  theTriggerResultsLabel_(iConfig.getParameter<edm::InputTag>("labelTriggerResult")),
252  labelGenTrack_(iConfig.getParameter<std::string>("labelTrack")),
253  labelRecVtx_(iConfig.getParameter<std::string>("labelVertex")),
254  labelEB_(iConfig.getParameter<std::string>("labelEBRecHit")),
255  labelEE_(iConfig.getParameter<std::string>("labelEERecHit")),
256  labelHBHE_(iConfig.getParameter<std::string>("labelHBHERecHit")),
257  labelTower_(iConfig.getParameter<std::string>("labelCaloTower")),
258  l1TrigName_(iConfig.getUntrackedParameter<std::string>("l1TrigName", "L1_SingleJet60")),
259  oldID_(iConfig.getUntrackedParameter<std::vector<int> >("oldID")),
260  newDepth_(iConfig.getUntrackedParameter<std::vector<int> >("newDepth")),
261  hep17_(iConfig.getUntrackedParameter<bool>("hep17")),
262  nRun_(0),
263  nLow_(0),
264  nHigh_(0),
265  hdc_(nullptr) {
266  usesResource(TFileService::kSharedResource);
267 
268  //now do whatever initialization is needed
269  const double isolationRadius(28.9), innerR(10.0), outerR(30.0);
270  reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality_);
271  selectionParameter_.minPt = iConfig.getParameter<double>("minTrackPt");
272  ;
273  selectionParameter_.minQuality = trackQuality_;
274  selectionParameter_.maxDxyPV = iConfig.getParameter<double>("maxDxyPV");
275  selectionParameter_.maxDzPV = iConfig.getParameter<double>("maxDzPV");
276  selectionParameter_.maxChi2 = iConfig.getParameter<double>("maxChi2");
277  selectionParameter_.maxDpOverP = iConfig.getParameter<double>("maxDpOverP");
278  selectionParameter_.minOuterHit = iConfig.getParameter<int>("minOuterHit");
279  selectionParameter_.minLayerCrossed = iConfig.getParameter<int>("minLayerCrossed");
280  selectionParameter_.maxInMiss = iConfig.getParameter<int>("maxInMiss");
281  selectionParameter_.maxOutMiss = iConfig.getParameter<int>("maxOutMiss");
282  a_charIsoR_ = a_coneR_ + isolationRadius;
283  a_coneR1_ = a_coneR_ + innerR;
284  a_coneR2_ = a_coneR_ + outerR;
285  // Different isolation cuts are described in DN-2016/029
286  // Tight cut uses 2 GeV; Loose cut uses 10 GeV
287  // Eta dependent cut uses (maxRestrictionP_ * exp(|ieta|*log(2.5)/18))
288  // with the factor for exponential slopeRestrictionP_ = log(2.5)/18
289  // maxRestrictionP_ = 8 GeV as came from a study
290  std::string labelBS = iConfig.getParameter<std::string>("labelBeamSpot");
291  std::string modnam = iConfig.getUntrackedParameter<std::string>("moduleName", "");
292  std::string prdnam = iConfig.getUntrackedParameter<std::string>("producerName", "");
293  edm::InputTag algTag = iConfig.getParameter<edm::InputTag>("algInputTag");
294  edm::InputTag extTag = iConfig.getParameter<edm::InputTag>("extInputTag");
295  std::string labelMuon = iConfig.getParameter<std::string>("labelMuon");
296 
297  for (unsigned int k = 0; k < oldID_.size(); ++k) {
298  oldDet_.emplace_back((oldID_[k] / 10000) % 10);
299  oldEta_.emplace_back((oldID_[k] / 100) % 100);
300  oldDepth_.emplace_back(oldID_[k] % 100);
301  }
302 
303  l1GtUtils_ = new l1t::L1TGlobalUtil(iConfig, consumesCollector(), *this, algTag, extTag, l1t::UseEventSetupIn::Event);
304  // define tokens for access
305  tok_trigEvt_ = consumes<trigger::TriggerEvent>(triggerEvent_);
306  tok_trigRes_ = consumes<edm::TriggerResults>(theTriggerResultsLabel_);
307  tok_bs_ = consumes<reco::BeamSpot>(labelBS);
308  tok_genTrack_ = consumes<reco::TrackCollection>(labelGenTrack_);
309  tok_ew_ = consumes<GenEventInfoProduct>(edm::InputTag("generator"));
310  tok_parts_ = consumes<reco::GenParticleCollection>(edm::InputTag("genParticles"));
311  tok_cala_ = consumes<CaloTowerCollection>(labelTower_);
312  tok_alg_ = consumes<BXVector<GlobalAlgBlk> >(algTag);
313  tok_Muon_ = consumes<reco::MuonCollection>(labelMuon);
314 
315  if (modnam.empty()) {
316  tok_recVtx_ = consumes<reco::VertexCollection>(labelRecVtx_);
317  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEB_));
318  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEE_));
319  tok_hbhe_ = consumes<HBHERecHitCollection>(labelHBHE_);
320  edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
321  << labelBS << " " << labelRecVtx_ << " " << labelGenTrack_ << " "
322  << edm::InputTag("ecalRecHit", labelEB_) << " "
323  << edm::InputTag("ecalRecHit", labelEE_) << " " << labelHBHE_ << " " << labelTower_
324  << " " << labelMuon;
325  } else {
326  tok_recVtx_ = consumes<reco::VertexCollection>(edm::InputTag(modnam, labelRecVtx_, prdnam));
327  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag(modnam, labelEB_, prdnam));
328  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag(modnam, labelEE_, prdnam));
329  tok_hbhe_ = consumes<HBHERecHitCollection>(edm::InputTag(modnam, labelHBHE_, prdnam));
330  edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
331  << labelBS << " " << edm::InputTag(modnam, labelRecVtx_, prdnam) << " "
332  << labelGenTrack_ << " " << edm::InputTag(modnam, labelEB_, prdnam) << " "
333  << edm::InputTag(modnam, labelEE_, prdnam) << " "
334  << edm::InputTag(modnam, labelHBHE_, prdnam) << " " << labelTower_ << " "
335  << labelMuon;
336  }
337 
338  tok_ddrec_ = esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>();
339  tok_bFieldH_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
340  tok_ecalChStatus_ = esConsumes<EcalChannelStatus, EcalChannelStatusRcd>();
341  tok_sevlv_ = esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>();
342  tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
343  tok_caloTopology_ = esConsumes<CaloTopology, CaloTopologyRecord>();
344  tok_htopo_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
345  tok_resp_ = esConsumes<HcalRespCorrs, HcalRespCorrsRcd>();
346 
347  edm::LogVerbatim("HcalIsoTrack")
348  << "Parameters read from config file \n"
349  << "\t minPt " << selectionParameter_.minPt << "\t theTrackQuality " << theTrackQuality_ << "\t minQuality "
350  << selectionParameter_.minQuality << "\t maxDxyPV " << selectionParameter_.maxDxyPV << "\t maxDzPV "
351  << selectionParameter_.maxDzPV << "\t maxChi2 " << selectionParameter_.maxChi2 << "\t maxDpOverP "
352  << selectionParameter_.maxDpOverP << "\t minOuterHit " << selectionParameter_.minOuterHit << "\t minLayerCrossed "
353  << selectionParameter_.minLayerCrossed << "\t maxInMiss " << selectionParameter_.maxInMiss << "\t maxOutMiss "
354  << selectionParameter_.maxOutMiss << "\t a_coneR " << a_coneR_ << ":" << a_coneR1_ << ":" << a_coneR2_
355  << "\t a_charIsoR " << a_charIsoR_ << "\t a_mipR " << a_mipR_ << "\t a_mipR2 " << a_mipR2_ << "\t a_mipR3 "
356  << a_mipR3_ << "\t a_mipR4 " << a_mipR4_ << "\t a_mipR5 " << a_mipR5_ << "\n pTrackMin_ " << pTrackMin_
357  << "\t eEcalMax_ " << eEcalMax_ << "\t maxRestrictionP_ " << maxRestrictionP_ << "\t slopeRestrictionP_ "
358  << slopeRestrictionP_ << "\t eIsolateStrong_ " << eIsolate1_ << "\t eIsolateSoft_ " << eIsolate2_
359  << "\t hcalScale_ " << hcalScale_ << "\n\t momentumLow_ " << pTrackLow_ << "\t prescaleLow_ " << prescaleLow_
360  << "\t momentumHigh_ " << pTrackHigh_ << "\t prescaleHigh_ " << prescaleHigh_ << "\n\t useRaw_ " << useRaw_
361  << "\t ignoreTrigger_ " << ignoreTrigger_ << "\n\t useL1Trigegr_ " << useL1Trigger_ << "\t dataType_ "
362  << dataType_ << "\t mode_ " << mode_ << "\t unCorrect_ " << unCorrect_ << "\t collapseDepth_ "
363  << collapseDepth_ << "\t L1TrigName_ " << l1TrigName_ << "\nThreshold for EB " << hitEthrEB_ << " EE "
364  << hitEthrEE0_ << ":" << hitEthrEE1_ << ":" << hitEthrEE2_ << ":" << hitEthrEE3_ << ":" << hitEthrEELo_ << ":"
365  << hitEthrEEHi_;
366  edm::LogVerbatim("HcalIsoTrack") << "Process " << processName_ << " L1Filter:" << l1Filter_
367  << " L2Filter:" << l2Filter_ << " L3Filter:" << l3Filter_;
368  for (unsigned int k = 0; k < trigNames_.size(); ++k) {
369  edm::LogVerbatim("HcalIsoTrack") << "Trigger[" << k << "] " << trigNames_[k];
370  }
371  edm::LogVerbatim("HcalIsoTrack") << oldID_.size() << " DetIDs to be corrected with HEP17 flag:" << hep17_;
372  for (unsigned int k = 0; k < oldID_.size(); ++k)
373  edm::LogVerbatim("HcalIsoTrack") << "[" << k << "] Det " << oldDet_[k] << " EtaAbs " << oldEta_[k] << " Depth "
374  << oldDepth_[k] << ":" << newDepth_[k];
375 
376  for (int i = 0; i < 10; i++)
377  phibins_.push_back(-M_PI + 0.1 * (2 * i + 1) * M_PI);
378  for (int i = 0; i < 8; ++i)
379  etabins_.push_back(-2.1 + 0.6 * i);
380  etadist_ = etabins_[1] - etabins_[0];
381  phidist_ = phibins_[1] - phibins_[0];
382  etahalfdist_ = 0.5 * etadist_;
383  phihalfdist_ = 0.5 * phidist_;
384  edm::LogVerbatim("HcalIsoTrack") << "EtaDist " << etadist_ << " " << etahalfdist_ << " PhiDist " << phidist_ << " "
385  << phihalfdist_;
386  unsigned int k1(0), k2(0);
387  for (auto phi : phibins_) {
388  edm::LogVerbatim("HcalIsoTrack") << "phibin_[" << k1 << "] " << phi;
389  ++k1;
390  }
391  for (auto eta : etabins_) {
392  edm::LogVerbatim("HcalIsoTrack") << "etabin_[" << k2 << "] " << eta;
393  ++k2;
394  }
395 }

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_, l1t::Event, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcalScale_, hep17_, hitEthrEB_, hitEthrEE0_, hitEthrEE1_, hitEthrEE2_, hitEthrEE3_, hitEthrEEHi_, hitEthrEELo_, mps_fire::i, ignoreTrigger_, HLT_FULL_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, mode_, newDepth_, oldDepth_, oldDet_, oldEta_, oldID_, phi, phibins_, phidist_, phihalfdist_, prescaleHigh_, prescaleLow_, processName_, pTrackHigh_, pTrackLow_, pTrackMin_, reco::TrackBase::qualityByName(), selectionParameter_, slopeRestrictionP_, AlCaHLTBitMon_QueryRunRegistry::string, theTrackQuality_, theTriggerResultsLabel_, tok_alg_, tok_bFieldH_, tok_bs_, tok_cala_, tok_caloTopology_, tok_ddrec_, tok_EB_, tok_ecalChStatus_, tok_EE_, tok_ew_, tok_genTrack_, tok_geom_, tok_hbhe_, tok_htopo_, tok_Muon_, tok_parts_, tok_recVtx_, tok_resp_, tok_sevlv_, tok_trigEvt_, tok_trigRes_, triggerEvent_, trigNames_, unCorrect_, useL1Trigger_, and useRaw_.

◆ ~HcalIsoTrkAnalyzer()

HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer ( )
inlineoverride

Definition at line 94 of file HcalIsoTrkAnalyzer.cc.

94 {}

Member Function Documentation

◆ analyze()

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

Implements edm::one::EDAnalyzerBase.

Definition at line 397 of file HcalIsoTrkAnalyzer.cc.

397  {
398  t_Run = iEvent.id().run();
399  t_Event = iEvent.id().event();
400  t_DataType = dataType_;
401 #ifdef EDM_ML_DEBUG
402  edm::LogVerbatim("HcalIsoTrack") << "Run " << t_Run << " Event " << t_Event << " type " << t_DataType
403  << " Luminosity " << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing();
404 #endif
405  //Get magnetic field and ECAL channel status
406  const MagneticField* bField = &iSetup.getData(tok_bFieldH_);
407  const EcalChannelStatus* theEcalChStatus = &iSetup.getData(tok_ecalChStatus_);
408  const EcalSeverityLevelAlgo* theEcalSevlv = &iSetup.getData(tok_sevlv_);
409 
410  // get handles to calogeometry and calotopology
411  const CaloGeometry* geo = &iSetup.getData(tok_geom_);
412  const CaloTopology* caloTopology = &iSetup.getData(tok_caloTopology_);
413  const HcalTopology* theHBHETopology = &iSetup.getData(tok_htopo_);
414  const HcalRespCorrs* resp = &iSetup.getData(tok_resp_);
415  HcalRespCorrs* respCorrs = new HcalRespCorrs(*resp);
416  respCorrs->setTopo(theHBHETopology);
417 
418  //=== genParticle information
419  edm::Handle<reco::GenParticleCollection> genParticles;
420  iEvent.getByToken(tok_parts_, genParticles);
421 
422  bool okC(true);
423  //Get track collection
424  edm::Handle<reco::TrackCollection> trkCollection;
425  iEvent.getByToken(tok_genTrack_, trkCollection);
426  if (!trkCollection.isValid()) {
427  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelGenTrack_;
428  okC = false;
429  }
430 
431  //Get muon collection
432  edm::Handle<reco::MuonCollection> muonh;
433  iEvent.getByToken(tok_Muon_, muonh);
434 
435  //event weight for FLAT sample
436  t_EventWeight = 1.0;
437  edm::Handle<GenEventInfoProduct> genEventInfo;
438  iEvent.getByToken(tok_ew_, genEventInfo);
439  if (genEventInfo.isValid())
440  t_EventWeight = genEventInfo->weight();
441 
442  //Define the best vertex and the beamspot
443  edm::Handle<reco::VertexCollection> recVtxs;
444  iEvent.getByToken(tok_recVtx_, recVtxs);
445  edm::Handle<reco::BeamSpot> beamSpotH;
446  iEvent.getByToken(tok_bs_, beamSpotH);
447  math::XYZPoint leadPV(0, 0, 0);
448  t_goodPV = t_nVtx = 0;
449  if (recVtxs.isValid() && !(recVtxs->empty())) {
450  t_nVtx = recVtxs->size();
451  for (unsigned int k = 0; k < recVtxs->size(); ++k) {
452  if (!((*recVtxs)[k].isFake()) && ((*recVtxs)[k].ndof() > 4)) {
453  if (t_goodPV == 0)
454  leadPV = math::XYZPoint((*recVtxs)[k].x(), (*recVtxs)[k].y(), (*recVtxs)[k].z());
455  t_goodPV++;
456  }
457  }
458  }
459  if (t_goodPV == 0 && beamSpotH.isValid()) {
460  leadPV = beamSpotH->position();
461  }
462  t_allvertex = t_goodPV;
463 #ifdef EDM_ML_DEBUG
464  edm::LogVerbatim("HcalIsoTrack") << "Primary Vertex " << leadPV << " out of " << t_goodPV << " vertex";
465  if (beamSpotH.isValid()) {
466  edm::LogVerbatim("HcalIsoTrack") << " Beam Spot " << beamSpotH->position();
467  }
468 #endif
469  // RecHits
470  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
471  iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
472  if (!barrelRecHitsHandle.isValid()) {
473  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEB_;
474  okC = false;
475  }
476  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
477  iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
478  if (!endcapRecHitsHandle.isValid()) {
479  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEE_;
480  okC = false;
481  }
482  edm::Handle<HBHERecHitCollection> hbhe;
483  iEvent.getByToken(tok_hbhe_, hbhe);
484  if (!hbhe.isValid()) {
485  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelHBHE_;
486  okC = false;
487  }
488  edm::Handle<CaloTowerCollection> caloTower;
489  iEvent.getByToken(tok_cala_, caloTower);
490 
491  //Propagate tracks to calorimeter surface)
492  std::vector<spr::propagatedTrackDirection> trkCaloDirections;
493  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDirections, false);
494  std::vector<spr::propagatedTrackID> trkCaloDets;
495  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDets, false);
496  std::vector<math::XYZTLorentzVector> vecL1, vecL3;
497  t_RunNo = iEvent.id().run();
498  t_EventNo = iEvent.id().event();
499  t_Tracks = trkCollection->size();
500  t_TracksProp = trkCaloDirections.size();
501  t_ietaAll->clear();
502  t_ietaGood->clear();
503  t_trackType->clear();
504  t_trgbits->clear();
505  t_hltbits->clear();
506 #ifdef EDM_ML_DEBUG
507  edm::LogVerbatim("HcalIsoTrack") << "# of propagated tracks " << t_TracksProp << " out of " << t_Tracks
508  << " with Trigger " << ignoreTrigger_;
509 #endif
510 
511  //Trigger
512  t_trgbits->assign(trigNames_.size(), false);
513  t_hltbits->assign(trigNames_.size(), false);
514  t_TracksSaved = t_TracksLoose = t_TracksTight = 0;
515  t_L1Bit = true;
516  t_TrigPass = false;
517 
518  //L1
519  l1GtUtils_->retrieveL1(iEvent, iSetup, tok_alg_);
520  const std::vector<std::pair<std::string, bool> >& finalDecisions = l1GtUtils_->decisionsFinal();
521  for (const auto& decision : finalDecisions) {
522  if (decision.first.find(l1TrigName_) != std::string::npos) {
523  t_L1Bit = decision.second;
524  break;
525  }
526  }
527 #ifdef EDM_ML_DEBUG
528  edm::LogVerbatim("HcalIsoTrack") << "Trigger Information for " << l1TrigName_ << " is " << t_L1Bit
529  << " from a list of " << finalDecisions.size() << " decisions";
530 #endif
531 
532  //HLT
533  edm::Handle<edm::TriggerResults> triggerResults;
534  iEvent.getByToken(tok_trigRes_, triggerResults);
535  if (triggerResults.isValid()) {
536  const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
537  const std::vector<std::string>& names = triggerNames.triggerNames();
538  if (!trigNames_.empty()) {
539  for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
540  int hlt = triggerResults->accept(iHLT);
541  for (unsigned int i = 0; i < trigNames_.size(); ++i) {
542  if (names[iHLT].find(trigNames_[i]) != std::string::npos) {
543  t_trgbits->at(i) = (hlt > 0);
544  t_hltbits->at(i) = (hlt > 0);
545  if (hlt > 0)
546  t_TrigPass = true;
547 #ifdef EDM_ML_DEBUG
548  edm::LogVerbatim("HcalIsoTrack")
549  << "This trigger " << names[iHLT] << " Flag " << hlt << ":" << t_trgbits->at(i);
550 #endif
551  }
552  }
553  }
554  }
555  }
556 #ifdef EDM_ML_DEBUG
557  edm::LogVerbatim("HcalIsoTrack") << "HLT Information shows " << t_TrigPass << ":" << trigNames_.empty() << ":" << okC;
558 #endif
559 
560  std::array<int, 3> ntksave{{0, 0, 0}};
561  if (ignoreTrigger_ || useL1Trigger_) {
562  t_l1pt = t_l1eta = t_l1phi = 0;
563  t_l3pt = t_l3eta = t_l3phi = 0;
564  if (ignoreTrigger_ || t_L1Bit)
565  ntksave = fillTree(vecL1,
566  vecL3,
567  leadPV,
568  trkCaloDirections,
569  trkCaloDets,
570  geo,
571  caloTopology,
572  theHBHETopology,
573  theEcalChStatus,
574  theEcalSevlv,
575  barrelRecHitsHandle,
576  endcapRecHitsHandle,
577  hbhe,
578  caloTower,
579  genParticles,
580  respCorrs,
581  muonh);
582  t_TracksSaved = ntksave[0];
583  t_TracksLoose = ntksave[1];
584  t_TracksTight = ntksave[2];
585  } else {
586  trigger::TriggerEvent triggerEvent;
587  edm::Handle<trigger::TriggerEvent> triggerEventHandle;
588  iEvent.getByToken(tok_trigEvt_, triggerEventHandle);
589  if (!triggerEventHandle.isValid()) {
590  edm::LogWarning("HcalIsoTrack") << "Error! Can't get the product " << triggerEvent_.label();
591  } else if (okC) {
592  triggerEvent = *(triggerEventHandle.product());
593  const trigger::TriggerObjectCollection& TOC(triggerEvent.getObjects());
594  bool done(false);
595  if (triggerResults.isValid()) {
596  std::vector<std::string> modules;
597  const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
598  const std::vector<std::string>& names = triggerNames.triggerNames();
599  for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
600  bool ok = (t_TrigPass) || (trigNames_.empty());
601  if (ok) {
602  unsigned int triggerindx = hltConfig_.triggerIndex(names[iHLT]);
603  const std::vector<std::string>& moduleLabels(hltConfig_.moduleLabels(triggerindx));
604  std::vector<math::XYZTLorentzVector> vecL2;
605  vecL1.clear();
606  vecL3.clear();
607  //loop over all trigger filters in event (i.e. filters passed)
608  for (unsigned int ifilter = 0; ifilter < triggerEvent.sizeFilters(); ++ifilter) {
609  std::vector<int> Keys;
610  std::string label = triggerEvent.filterTag(ifilter).label();
611  //loop over keys to objects passing this filter
612  for (unsigned int imodule = 0; imodule < moduleLabels.size(); imodule++) {
613  if (label.find(moduleLabels[imodule]) != std::string::npos) {
614 #ifdef EDM_ML_DEBUG
615  edm::LogVerbatim("HcalIsoTrack") << "FilterName " << label;
616 #endif
617  for (unsigned int ifiltrKey = 0; ifiltrKey < triggerEvent.filterKeys(ifilter).size(); ++ifiltrKey) {
618  Keys.push_back(triggerEvent.filterKeys(ifilter)[ifiltrKey]);
619  const trigger::TriggerObject& TO(TOC[Keys[ifiltrKey]]);
620  math::XYZTLorentzVector v4(TO.px(), TO.py(), TO.pz(), TO.energy());
621  if (label.find(l2Filter_) != std::string::npos) {
622  vecL2.push_back(v4);
623  } else if (label.find(l3Filter_) != std::string::npos) {
624  vecL3.push_back(v4);
625  } else if ((label.find(l1Filter_) != std::string::npos) || (l1Filter_.empty())) {
626  vecL1.push_back(v4);
627  }
628 #ifdef EDM_ML_DEBUG
629  edm::LogVerbatim("HcalIsoTrack")
630  << "key " << ifiltrKey << " : pt " << TO.pt() << " eta " << TO.eta() << " phi " << TO.phi()
631  << " mass " << TO.mass() << " Id " << TO.id();
632 #endif
633  }
634 #ifdef EDM_ML_DEBUG
635  edm::LogVerbatim("HcalIsoTrack")
636  << "sizes " << vecL1.size() << ":" << vecL2.size() << ":" << vecL3.size();
637 #endif
638  }
639  }
640  }
642  math::XYZTLorentzVector mindRvec1;
643  double mindR1(999);
644  for (unsigned int i = 0; i < vecL2.size(); i++) {
645  double dr = dR(vecL1[0], vecL2[i]);
646 #ifdef EDM_ML_DEBUG
647  edm::LogVerbatim("HcalIsoTrack") << "lvl2[" << i << "] dR " << dr;
648 #endif
649  if (dr < mindR1) {
650  mindR1 = dr;
651  mindRvec1 = vecL2[i];
652  }
653  }
654 #ifdef EDM_ML_DEBUG
655  edm::LogVerbatim("HcalIsoTrack") << "L2 object closest to L1 " << mindRvec1 << " at Dr " << mindR1;
656 #endif
657 
658  if (!vecL1.empty()) {
659  t_l1pt = vecL1[0].pt();
660  t_l1eta = vecL1[0].eta();
661  t_l1phi = vecL1[0].phi();
662  } else {
663  t_l1pt = t_l1eta = t_l1phi = 0;
664  }
665  if (!vecL3.empty()) {
666  t_l3pt = vecL3[0].pt();
667  t_l3eta = vecL3[0].eta();
668  t_l3phi = vecL3[0].phi();
669  } else {
670  t_l3pt = t_l3eta = t_l3phi = 0;
671  }
672  // Now fill in the tree for each selected track
673  if (!done) {
674  ntksave = fillTree(vecL1,
675  vecL3,
676  leadPV,
677  trkCaloDirections,
678  trkCaloDets,
679  geo,
680  caloTopology,
681  theHBHETopology,
682  theEcalChStatus,
683  theEcalSevlv,
684  barrelRecHitsHandle,
685  endcapRecHitsHandle,
686  hbhe,
687  caloTower,
688  genParticles,
689  respCorrs,
690  muonh);
691  t_TracksSaved += ntksave[0];
692  t_TracksLoose += ntksave[1];
693  t_TracksTight += ntksave[2];
694  done = true;
695  }
696  }
697  }
698  }
699  }
700  }
701 #ifdef EDM_ML_DEBUG
702  edm::LogVerbatim("HcalIsoTrack") << "Final results on selected tracks " << t_TracksSaved << ":" << t_TracksLoose
703  << ":" << t_TracksTight;
704 #endif
705  t_TrigPassSel = (t_TracksSaved > 0);
706  tree2->Fill();
707 }

References Calorimetry_cff::bField, dataType_, l1t::L1TGlobalUtil::decisionsFinal(), fileCollector::done, flavorHistoryFilter_cfi::dr, dR(), trigger::TriggerObject::energy(), trigger::TriggerObject::eta(), fillTree(), spr::find(), genParticles2HepMC_cfi::genEventInfo, genParticles2HepMC_cfi::genParticles, edm::EventSetup::getData(), photonIsolationHIProducer_cfi::hbhe, ValidationMatrix::hlt, hltConfig_, mps_fire::i, trigger::TriggerObject::id(), iEvent, ignoreTrigger_, edm::HandleBase::isValid(), dqmdumpme::k, l1Filter_, l1GtUtils_, l1TrigName_, l2Filter_, l3Filter_, label, edm::InputTag::label(), labelEB_, labelEE_, labelGenTrack_, labelHBHE_, trigger::TriggerObject::mass(), HLTConfigProvider::moduleLabels(), LogMessageMonitor_cff::modules, names, convertSQLiteXML::ok, trigger::TriggerObject::phi(), reco::BeamSpot::position(), edm::Handle< T >::product(), spr::propagateCALO(), trigger::TriggerObject::pt(), trigger::TriggerObject::px(), trigger::TriggerObject::py(), trigger::TriggerObject::pz(), l1t::L1TGlobalUtil::retrieveL1(), HcalCondObjectContainerBase::setTopo(), 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_bFieldH_, tok_bs_, tok_cala_, tok_caloTopology_, tok_EB_, tok_ecalChStatus_, tok_EE_, tok_ew_, tok_genTrack_, tok_geom_, tok_hbhe_, tok_htopo_, tok_Muon_, tok_parts_, tok_recVtx_, tok_resp_, tok_sevlv_, tok_trigEvt_, tok_trigRes_, tree2, PDWG_DiPhoton_SD_cff::triggerEvent, triggerEvent_, HLTConfigProvider::triggerIndex(), L1TEGammaOffline_cfi::triggerNames, triggerResults, trigNames_, useL1Trigger_, x, y, and z.

◆ beginJob()

void HcalIsoTrkAnalyzer::beginJob ( void  )
overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 709 of file HcalIsoTrkAnalyzer.cc.

709  {
710  tree = fs->make<TTree>("CalibTree", "CalibTree");
711 
712  tree->Branch("t_Run", &t_Run, "t_Run/I");
713  tree->Branch("t_Event", &t_Event, "t_Event/I");
714  tree->Branch("t_DataType", &t_DataType, "t_DataType/I");
715  tree->Branch("t_ieta", &t_ieta, "t_ieta/I");
716  tree->Branch("t_iphi", &t_iphi, "t_iphi/I");
717  tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
718  tree->Branch("t_nVtx", &t_nVtx, "t_nVtx/I");
719  tree->Branch("t_nTrk", &t_nTrk, "t_nTrk/I");
720  tree->Branch("t_goodPV", &t_goodPV, "t_goodPV/I");
721  tree->Branch("t_l1pt", &t_l1pt, "t_l1pt/D");
722  tree->Branch("t_l1eta", &t_l1eta, "t_l1eta/D");
723  tree->Branch("t_l1phi", &t_l1phi, "t_l1phi/D");
724  tree->Branch("t_l3pt", &t_l3pt, "t_l3pt/D");
725  tree->Branch("t_l3eta", &t_l3eta, "t_l3eta/D");
726  tree->Branch("t_l3phi", &t_l3phi, "t_l3phi/D");
727  tree->Branch("t_p", &t_p, "t_p/D");
728  tree->Branch("t_pt", &t_pt, "t_pt/D");
729  tree->Branch("t_phi", &t_phi, "t_phi/D");
730  tree->Branch("t_mindR1", &t_mindR1, "t_mindR1/D");
731  tree->Branch("t_mindR2", &t_mindR2, "t_mindR2/D");
732  tree->Branch("t_eMipDR", &t_eMipDR, "t_eMipDR/D");
733  tree->Branch("t_eMipDR2", &t_eMipDR2, "t_eMipDR2/D");
734  tree->Branch("t_eMipDR3", &t_eMipDR3, "t_eMipDR3/D");
735  tree->Branch("t_eMipDR4", &t_eMipDR4, "t_eMipDR4/D");
736  tree->Branch("t_eMipDR5", &t_eMipDR5, "t_eMipDR5/D");
737  tree->Branch("t_eHcal", &t_eHcal, "t_eHcal/D");
738  tree->Branch("t_eHcal10", &t_eHcal10, "t_eHcal10/D");
739  tree->Branch("t_eHcal30", &t_eHcal30, "t_eHcal30/D");
740  tree->Branch("t_hmaxNearP", &t_hmaxNearP, "t_hmaxNearP/D");
741  tree->Branch("t_emaxNearP", &t_emaxNearP, "t_emaxNearP/D");
742  tree->Branch("t_eAnnular", &t_eAnnular, "t_eAnnular/D");
743  tree->Branch("t_hAnnular", &t_hAnnular, "t_hAnnular/D");
744  tree->Branch("t_rhoh", &t_rhoh, "t_rhoh/D");
745  tree->Branch("t_selectTk", &t_selectTk, "t_selectTk/O");
746  tree->Branch("t_qltyFlag", &t_qltyFlag, "t_qltyFlag/O");
747  tree->Branch("t_qltyMissFlag", &t_qltyMissFlag, "t_qltyMissFlag/O");
748  tree->Branch("t_qltyPVFlag", &t_qltyPVFlag, "t_qltyPVFlag/O");
749  tree->Branch("t_gentrackP", &t_gentrackP, "t_gentrackP/D");
750 
751  t_DetIds = new std::vector<unsigned int>();
752  t_DetIds1 = new std::vector<unsigned int>();
753  t_DetIds3 = new std::vector<unsigned int>();
754  t_HitEnergies = new std::vector<double>();
755  t_HitEnergies1 = new std::vector<double>();
756  t_HitEnergies3 = new std::vector<double>();
757  t_trgbits = new std::vector<bool>();
758  tree->Branch("t_DetIds", "std::vector<unsigned int>", &t_DetIds);
759  tree->Branch("t_HitEnergies", "std::vector<double>", &t_HitEnergies);
760  tree->Branch("t_trgbits", "std::vector<bool>", &t_trgbits);
761  tree->Branch("t_DetIds1", "std::vector<unsigned int>", &t_DetIds1);
762  tree->Branch("t_DetIds3", "std::vector<unsigned int>", &t_DetIds3);
763  tree->Branch("t_HitEnergies1", "std::vector<double>", &t_HitEnergies1);
764  tree->Branch("t_HitEnergies3", "std::vector<double>", &t_HitEnergies3);
765 
766  tree2 = fs->make<TTree>("EventInfo", "Event Information");
767 
768  tree2->Branch("t_RunNo", &t_RunNo, "t_RunNo/i");
769  tree2->Branch("t_EventNo", &t_EventNo, "t_EventNo/i");
770  tree2->Branch("t_Tracks", &t_Tracks, "t_Tracks/I");
771  tree2->Branch("t_TracksProp", &t_TracksProp, "t_TracksProp/I");
772  tree2->Branch("t_TracksSaved", &t_TracksSaved, "t_TracksSaved/I");
773  tree2->Branch("t_TracksLoose", &t_TracksLoose, "t_TracksLoose/I");
774  tree2->Branch("t_TracksTight", &t_TracksTight, "t_TracksTight/I");
775  tree2->Branch("t_TrigPass", &t_TrigPass, "t_TrigPass/O");
776  tree2->Branch("t_TrigPassSel", &t_TrigPassSel, "t_TrigPassSel/O");
777  tree2->Branch("t_L1Bit", &t_L1Bit, "t_L1Bit/O");
778  tree2->Branch("t_allvertex", &t_allvertex, "t_allvertex/I");
779  t_hltbits = new std::vector<bool>();
780  t_ietaAll = new std::vector<int>();
781  t_ietaGood = new std::vector<int>();
782  t_trackType = new std::vector<int>();
783  tree2->Branch("t_ietaAll", "std::vector<int>", &t_ietaAll);
784  tree2->Branch("t_ietaGood", "std::vector<int>", &t_ietaGood);
785  tree2->Branch("t_trackType", "std::vector<int>", &t_trackType);
786  tree2->Branch("t_hltbits", "std::vector<bool>", &t_hltbits);
787 }

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

◆ beginRun()

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

Definition at line 790 of file HcalIsoTrkAnalyzer.cc.

790  {
791  hdc_ = &iSetup.getData(tok_ddrec_);
792 
793  bool changed_(true);
794  bool flag = hltConfig_.init(iRun, iSetup, processName_, changed_);
795  edm::LogVerbatim("HcalIsoTrack") << "Run[" << nRun_ << "] " << iRun.run() << " process " << processName_
796  << " init flag " << flag << " change flag " << changed_;
797  // check if trigger names in (new) config
798  if (changed_) {
799 #ifdef EDM_ML_DEBUG
800  edm::LogVerbatim("HcalIsoTrack") << "New trigger menu found !!!";
801 #endif
802  const unsigned int n(hltConfig_.size());
803  for (unsigned itrig = 0; itrig < trigNames_.size(); itrig++) {
804  unsigned int triggerindx = hltConfig_.triggerIndex(trigNames_[itrig]);
805  if (triggerindx >= n) {
806  edm::LogWarning("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " does not exist in "
807  << "the current menu";
808 #ifdef EDM_ML_DEBUG
809  } else {
810  edm::LogVerbatim("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " exists";
811 #endif
812  }
813  }
814  }
815 }

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

◆ dR()

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

Definition at line 1282 of file HcalIsoTrkAnalyzer.cc.

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

References reco::deltaR().

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

◆ endRun()

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

Definition at line 818 of file HcalIsoTrkAnalyzer.cc.

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

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

◆ eThreshold()

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

Definition at line 1333 of file HcalIsoTrkAnalyzer.cc.

1333  {
1334  const GlobalPoint& pos = geo->getPosition(id);
1335  double eta = std::abs(pos.eta());
1336  double eThr(hitEthrEB_);
1337  if (id.subdetId() != EcalBarrel) {
1338  eThr = (((eta * hitEthrEE3_ + hitEthrEE2_) * eta + hitEthrEE1_) * eta + hitEthrEE0_);
1339  if (eThr < hitEthrEELo_)
1340  eThr = hitEthrEELo_;
1341  else if (eThr > hitEthrEEHi_)
1342  eThr = hitEthrEEHi_;
1343  }
1344  return eThr;
1345 }

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

Referenced by fillTree().

◆ fillDescriptions()

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

Definition at line 823 of file HcalIsoTrkAnalyzer.cc.

823  {
824  edm::ParameterSetDescription desc;
825  std::vector<std::string> trig = {"HLT_PFJet40",
826  "HLT_PFJet60",
827  "HLT_PFJet80",
828  "HLT_PFJet140",
829  "HLT_PFJet200",
830  "HLT_PFJet260",
831  "HLT_PFJet320",
832  "HLT_PFJet400",
833  "HLT_PFJet450",
834  "HLT_PFJet500"};
835  desc.add<std::vector<std::string> >("triggers", trig);
836  desc.add<std::string>("processName", "HLT");
837  desc.add<std::string>("l1Filter", "");
838  desc.add<std::string>("l2Filter", "L2Filter");
839  desc.add<std::string>("l3Filter", "Filter");
840  // following 10 parameters are parameters to select good tracks
841  desc.add<std::string>("trackQuality", "highPurity");
842  desc.add<double>("minTrackPt", 1.0);
843  desc.add<double>("maxDxyPV", 0.02);
844  desc.add<double>("maxDzPV", 0.02);
845  desc.add<double>("maxChi2", 5.0);
846  desc.add<double>("maxDpOverP", 0.1);
847  desc.add<int>("minOuterHit", 4);
848  desc.add<int>("minLayerCrossed", 8);
849  desc.add<int>("maxInMiss", 0);
850  desc.add<int>("maxOutMiss", 0);
851  // Minimum momentum of selected isolated track and signal zone
852  desc.add<double>("minimumTrackP", 10.0);
853  desc.add<double>("coneRadius", 34.98);
854  // signal zone in ECAL and MIP energy cutoff
855  desc.add<double>("coneRadiusMIP", 14.0);
856  desc.add<double>("coneRadiusMIP2", 18.0);
857  desc.add<double>("coneRadiusMIP3", 20.0);
858  desc.add<double>("coneRadiusMIP4", 22.0);
859  desc.add<double>("coneRadiusMIP5", 24.0);
860  desc.add<double>("maximumEcalEnergy", 2.0);
861  // following 4 parameters are for isolation cuts and described in the code
862  desc.add<double>("maxTrackP", 8.0);
863  desc.add<double>("slopeTrackP", 0.05090504066);
864  desc.add<double>("isolationEnergyTight", 2.0);
865  desc.add<double>("isolationEnergyLoose", 10.0);
866  // energy thershold for ECAL (from Egamma group)
867  desc.add<double>("EBHitEnergyThreshold", 0.08);
868  desc.add<double>("EEHitEnergyThreshold0", 0.30);
869  desc.add<double>("EEHitEnergyThreshold1", 0.00);
870  desc.add<double>("EEHitEnergyThreshold2", 0.00);
871  desc.add<double>("EEHitEnergyThreshold3", 0.00);
872  desc.add<double>("EEHitEnergyThresholdLow", 0.30);
873  desc.add<double>("EEHitEnergyThresholdHigh", 0.30);
874  // prescale factors
875  desc.add<double>("momentumLow", 40.0);
876  desc.add<double>("momentumHigh", 60.0);
877  desc.add<int>("prescaleLow", 1);
878  desc.add<int>("prescaleHigh", 1);
879  // various labels for collections used in the code
880  desc.add<edm::InputTag>("labelTriggerEvent", edm::InputTag("hltTriggerSummaryAOD", "", "HLT"));
881  desc.add<edm::InputTag>("labelTriggerResult", edm::InputTag("TriggerResults", "", "HLT"));
882  desc.add<std::string>("labelTrack", "generalTracks");
883  desc.add<std::string>("labelVertex", "offlinePrimaryVertices");
884  desc.add<std::string>("labelEBRecHit", "EcalRecHitsEB");
885  desc.add<std::string>("labelEERecHit", "EcalRecHitsEE");
886  desc.add<std::string>("labelHBHERecHit", "hbhereco");
887  desc.add<std::string>("labelBeamSpot", "offlineBeamSpot");
888  desc.add<std::string>("labelCaloTower", "towerMaker");
889  desc.add<std::string>("labelMuon", "muons");
890  desc.add<edm::InputTag>("algInputTag", edm::InputTag("gtStage2Digis"));
891  desc.add<edm::InputTag>("extInputTag", edm::InputTag("gtStage2Digis"));
892  desc.addUntracked<std::string>("moduleName", "");
893  desc.addUntracked<std::string>("producerName", "");
894  // Various flags used for selecting tracks, choice of energy Method2/0
895  // Data type 0/1 for single jet trigger or others
896  desc.addUntracked<int>("useRaw", 0);
897  desc.addUntracked<bool>("ignoreTriggers", false);
898  desc.addUntracked<bool>("useL1Trigger", false);
899  desc.addUntracked<double>("hcalScale", 1.0);
900  desc.addUntracked<int>("dataType", 0);
901  desc.addUntracked<bool>("unCorrect", false);
902  desc.addUntracked<bool>("collapseDepth", false);
903  desc.addUntracked<std::string>("l1TrigName", "L1_SingleJet60");
904  desc.addUntracked<int>("outMode", 11);
905  std::vector<int> dummy;
906  desc.addUntracked<std::vector<int> >("oldID", dummy);
907  desc.addUntracked<std::vector<int> >("newDepth", dummy);
908  desc.addUntracked<bool>("hep17", false);
909  descriptions.add("HcalIsoTrkAnalyzer", desc);
910 }

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

◆ fillTree()

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

Definition at line 912 of file HcalIsoTrkAnalyzer.cc.

928  {
929  int nSave(0), nLoose(0), nTight(0);
930  //Loop over tracks
931  std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
932  unsigned int nTracks(0), nselTracks(0);
933  t_nTrk = trkCaloDirections.size();
934  t_rhoh = (tower.isValid()) ? rhoh(tower) : 0;
935  for (trkDetItr = trkCaloDirections.begin(), nTracks = 0; trkDetItr != trkCaloDirections.end();
936  trkDetItr++, nTracks++) {
937  const reco::Track* pTrack = &(*(trkDetItr->trkItr));
938  math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
939 #ifdef EDM_ML_DEBUG
940  edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) :" << pTrack->pt() << "|"
941  << pTrack->eta() << "|" << pTrack->phi() << "|" << pTrack->p();
942 #endif
943  t_mindR2 = 999;
944  for (unsigned int k = 0; k < vecL3.size(); ++k) {
945  double dr = dR(vecL3[k], v4);
946  if (dr < t_mindR2) {
947  t_mindR2 = dr;
948  }
949  }
950  t_mindR1 = (!vecL1.empty()) ? dR(vecL1[0], v4) : 999;
951 #ifdef EDM_ML_DEBUG
952  edm::LogVerbatim("HcalIsoTrack") << "Closest L3 object at dr :" << t_mindR2 << " and from L1 " << t_mindR1;
953 #endif
954  t_ieta = t_iphi = 0;
955  if (trkDetItr->okHCAL) {
956  HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
957  t_ieta = detId.ieta();
958  t_iphi = detId.iphi();
959  if (t_p > 40.0 && t_p <= 60.0)
960  t_ietaAll->emplace_back(t_ieta);
961  }
962  //Selection of good track
963  t_selectTk = spr::goodTrack(pTrack, leadPV, selectionParameter_, false);
964  spr::trackSelectionParameters oneCutParameters = selectionParameter_;
965  oneCutParameters.maxDxyPV = 10;
966  oneCutParameters.maxDzPV = 100;
967  oneCutParameters.maxInMiss = 2;
968  oneCutParameters.maxOutMiss = 2;
969  bool qltyFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
970  oneCutParameters = selectionParameter_;
971  oneCutParameters.maxDxyPV = 10;
972  oneCutParameters.maxDzPV = 100;
973  t_qltyMissFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
974  oneCutParameters = selectionParameter_;
975  oneCutParameters.maxInMiss = 2;
976  oneCutParameters.maxOutMiss = 2;
977  t_qltyPVFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
978  double eIsolation = maxRestrictionP_ * exp(slopeRestrictionP_ * std::abs((double)t_ieta));
979  if (eIsolation < eIsolate1_)
980  eIsolation = eIsolate1_;
981  if (eIsolation < eIsolate2_)
982  eIsolation = eIsolate2_;
983 #ifdef EDM_ML_DEBUG
984  edm::LogVerbatim("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " << qltyFlag << "|" << trkDetItr->okECAL << "|"
985  << trkDetItr->okHCAL << " eIsolation " << eIsolation;
986 #endif
987  t_qltyFlag = (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL);
988  bool notMuon = (muonh.isValid()) ? notaMuon(pTrack, muonh) : true;
989  if (t_qltyFlag && notMuon) {
990  nselTracks++;
991  int nNearTRKs(0);
993  std::vector<DetId> eIds;
994  std::vector<double> eHit;
995  t_eMipDR = spr::eCone_ecal(geo,
996  barrelRecHitsHandle,
997  endcapRecHitsHandle,
998  trkDetItr->pointHCAL,
999  trkDetItr->pointECAL,
1000  a_mipR_,
1001  trkDetItr->directionECAL,
1002  eIds,
1003  eHit);
1004  double eEcal(0);
1005  for (unsigned int k = 0; k < eIds.size(); ++k) {
1006  if (eHit[k] > eThreshold(eIds[k], geo))
1007  eEcal += eHit[k];
1008  }
1009 #ifdef EDM_ML_DEBUG
1010  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR << ":" << eEcal;
1011 #endif
1012  t_eMipDR = eEcal;
1015  std::vector<DetId> eIds2;
1016  std::vector<double> eHit2;
1017  t_eMipDR2 = spr::eCone_ecal(geo,
1018  barrelRecHitsHandle,
1019  endcapRecHitsHandle,
1020  trkDetItr->pointHCAL,
1021  trkDetItr->pointECAL,
1022  a_mipR2_,
1023  trkDetItr->directionECAL,
1024  eIds2,
1025  eHit2);
1026  double eEcal2(0);
1027  for (unsigned int k = 0; k < eIds2.size(); ++k) {
1028  if (eHit2[k] > eThreshold(eIds2[k], geo))
1029  eEcal2 += eHit2[k];
1030  }
1031 #ifdef EDM_ML_DEBUG
1032  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR2 << ":" << eEcal2;
1033 #endif
1034  t_eMipDR2 = eEcal2;
1037  std::vector<DetId> eIds3;
1038  std::vector<double> eHit3;
1039  t_eMipDR3 = spr::eCone_ecal(geo,
1040  barrelRecHitsHandle,
1041  endcapRecHitsHandle,
1042  trkDetItr->pointHCAL,
1043  trkDetItr->pointECAL,
1044  a_mipR3_,
1045  trkDetItr->directionECAL,
1046  eIds3,
1047  eHit3);
1048  double eEcal3(0);
1049  for (unsigned int k = 0; k < eIds3.size(); ++k) {
1050  if (eHit3[k] > eThreshold(eIds3[k], geo))
1051  eEcal3 += eHit3[k];
1052  }
1053 #ifdef EDM_ML_DEBUG
1054  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR3 << ":" << eEcal3;
1055 #endif
1056  t_eMipDR3 = eEcal3;
1059  std::vector<DetId> eIds4;
1060  std::vector<double> eHit4;
1061  t_eMipDR4 = spr::eCone_ecal(geo,
1062  barrelRecHitsHandle,
1063  endcapRecHitsHandle,
1064  trkDetItr->pointHCAL,
1065  trkDetItr->pointECAL,
1066  a_mipR4_,
1067  trkDetItr->directionECAL,
1068  eIds4,
1069  eHit4);
1070  double eEcal4(0);
1071  for (unsigned int k = 0; k < eIds4.size(); ++k) {
1072  if (eHit4[k] > eThreshold(eIds4[k], geo))
1073  eEcal4 += eHit4[k];
1074  }
1075 #ifdef EDM_ML_DEBUG
1076  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR4 << ":" << eEcal4;
1077 #endif
1078  t_eMipDR4 = eEcal4;
1081  std::vector<DetId> eIds5;
1082  std::vector<double> eHit5;
1083  t_eMipDR5 = spr::eCone_ecal(geo,
1084  barrelRecHitsHandle,
1085  endcapRecHitsHandle,
1086  trkDetItr->pointHCAL,
1087  trkDetItr->pointECAL,
1088  a_mipR5_,
1089  trkDetItr->directionECAL,
1090  eIds5,
1091  eHit5);
1092  double eEcal5(0);
1093  for (unsigned int k = 0; k < eIds5.size(); ++k) {
1094  if (eHit5[k] > eThreshold(eIds5[k], geo))
1095  eEcal5 += eHit5[k];
1096  }
1097 #ifdef EDM_ML_DEBUG
1098  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR5 << ":" << eEcal5;
1099 #endif
1100  t_eMipDR5 = eEcal5;
1102 
1103  t_emaxNearP = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 15, 15);
1104  const DetId cellE(trkDetItr->detIdECAL);
1105  std::pair<double, bool> e11x11P = spr::eECALmatrix(cellE,
1106  barrelRecHitsHandle,
1107  endcapRecHitsHandle,
1108  *theEcalChStatus,
1109  geo,
1110  caloTopology,
1111  theEcalSevlv,
1112  5,
1113  5,
1114  -100.0,
1115  -100.0,
1116  -100.0,
1117  100.0);
1118  std::pair<double, bool> e15x15P = spr::eECALmatrix(cellE,
1119  barrelRecHitsHandle,
1120  endcapRecHitsHandle,
1121  *theEcalChStatus,
1122  geo,
1123  caloTopology,
1124  theEcalSevlv,
1125  7,
1126  7,
1127  -100.0,
1128  -100.0,
1129  -100.0,
1130  100.0);
1131  if (e11x11P.second && e15x15P.second) {
1132  t_eAnnular = (e15x15P.first - e11x11P.first);
1133  } else {
1134  t_eAnnular = -(e15x15P.first - e11x11P.first);
1135  }
1136  t_hmaxNearP = spr::chargeIsolationCone(nTracks, trkCaloDirections, a_charIsoR_, nNearTRKs, false);
1137  const DetId cellH(trkDetItr->detIdHCAL);
1138  double h5x5 = spr::eHCALmatrix(
1139  theHBHETopology, cellH, hbhe, 2, 2, false, true, -100.0, -100.0, -100.0, -100.0, -100.0, 100.0);
1140  double h7x7 = spr::eHCALmatrix(
1141  theHBHETopology, cellH, hbhe, 3, 3, false, true, -100.0, -100.0, -100.0, -100.0, -100.0, 100.0);
1142  t_hAnnular = h7x7 - h5x5;
1143 #ifdef EDM_ML_DEBUG
1144  edm::LogVerbatim("HcalIsoTrack") << "max p Near (Ecal) " << t_emaxNearP << " (Hcal) " << t_hmaxNearP
1145  << " Annular E (Ecal) " << e11x11P.first << ":" << e15x15P.first << ":"
1146  << t_eAnnular << " (Hcal) " << h5x5 << ":" << h7x7 << ":" << t_hAnnular;
1147 #endif
1148  t_gentrackP = trackP(pTrack, genParticles);
1149  if (t_eMipDR < eEcalMax_ && t_hmaxNearP < eIsolation) {
1150  t_DetIds->clear();
1151  t_HitEnergies->clear();
1152  t_DetIds1->clear();
1153  t_HitEnergies1->clear();
1154  t_DetIds3->clear();
1155  t_HitEnergies3->clear();
1156  int nRecHits(-999), nRecHits1(-999), nRecHits3(-999);
1157  std::vector<DetId> ids, ids1, ids3;
1158  std::vector<double> edet0, edet1, edet3;
1159  t_eHcal = spr::eCone_hcal(geo,
1160  hbhe,
1161  trkDetItr->pointHCAL,
1162  trkDetItr->pointECAL,
1163  a_coneR_,
1164  trkDetItr->directionHCAL,
1165  nRecHits,
1166  ids,
1167  edet0,
1168  useRaw_);
1169  if (!oldID_.empty()) {
1170  for (unsigned k = 0; k < ids.size(); ++k)
1171  ids[k] = newId(ids[k]);
1172  }
1173  storeEnergy(0, respCorrs, ids, edet0, t_eHcal, t_DetIds, t_HitEnergies);
1174 
1175  //----- hcal energy in the extended cone 1 (a_coneR+10) --------------
1176  t_eHcal10 = spr::eCone_hcal(geo,
1177  hbhe,
1178  trkDetItr->pointHCAL,
1179  trkDetItr->pointECAL,
1180  a_coneR1_,
1181  trkDetItr->directionHCAL,
1182  nRecHits1,
1183  ids1,
1184  edet1,
1185  useRaw_);
1186  if (!oldID_.empty()) {
1187  for (unsigned k = 0; k < ids1.size(); ++k)
1188  ids1[k] = newId(ids1[k]);
1189  }
1190  storeEnergy(1, respCorrs, ids1, edet1, t_eHcal10, t_DetIds1, t_HitEnergies1);
1191 
1192  //----- hcal energy in the extended cone 3 (a_coneR+30) --------------
1193  t_eHcal30 = spr::eCone_hcal(geo,
1194  hbhe,
1195  trkDetItr->pointHCAL,
1196  trkDetItr->pointECAL,
1197  a_coneR2_,
1198  trkDetItr->directionHCAL,
1199  nRecHits3,
1200  ids3,
1201  edet3,
1202  useRaw_);
1203  if (!oldID_.empty()) {
1204  for (unsigned k = 0; k < ids3.size(); ++k)
1205  ids3[k] = newId(ids3[k]);
1206  }
1207  storeEnergy(3, respCorrs, ids3, edet3, t_eHcal30, t_DetIds3, t_HitEnergies3);
1208 
1209  t_p = pTrack->p();
1210  t_pt = pTrack->pt();
1211  t_phi = pTrack->phi();
1212 
1213 #ifdef EDM_ML_DEBUG
1214  edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) :" << t_pt << "|"
1215  << pTrack->eta() << "|" << t_phi << "|" << t_p << " Generator Level p "
1216  << t_gentrackP;
1217  edm::LogVerbatim("HcalIsoTrack") << "e_MIP " << t_eMipDR << " Chg Isolation " << t_hmaxNearP << " eHcal"
1218  << t_eHcal << " ieta " << t_ieta << " Quality " << t_qltyMissFlag << ":"
1219  << t_qltyPVFlag << ":" << t_selectTk;
1220  for (unsigned int ll = 0; ll < t_DetIds->size(); ll++) {
1221  edm::LogVerbatim("HcalIsoTrack")
1222  << "det id is = " << t_DetIds->at(ll) << " hit enery is = " << t_HitEnergies->at(ll);
1223  }
1224  for (unsigned int ll = 0; ll < t_DetIds1->size(); ll++) {
1225  edm::LogVerbatim("HcalIsoTrack")
1226  << "det id is = " << t_DetIds1->at(ll) << " hit enery is = " << t_HitEnergies1->at(ll);
1227  }
1228  for (unsigned int ll = 0; ll < t_DetIds3->size(); ll++) {
1229  edm::LogVerbatim("HcalIsoTrack")
1230  << "det id is = " << t_DetIds3->at(ll) << " hit enery is = " << t_HitEnergies3->at(ll);
1231  }
1232 #endif
1233  bool accept(false);
1234  if (t_p > pTrackMin_) {
1235  if (t_p < pTrackLow_) {
1236  ++nLow_;
1237  if (prescaleLow_ <= 1)
1238  accept = true;
1239  else if (nLow_ % prescaleLow_ == 1)
1240  accept = true;
1241  } else if (t_p > pTrackHigh_) {
1242  ++nHigh_;
1243  if (prescaleHigh_ <= 1)
1244  accept = true;
1245  else if (nHigh_ % prescaleHigh_ == 1)
1246  accept = true;
1247  } else {
1248  accept = true;
1249  }
1250  }
1251  if (accept) {
1252  tree->Fill();
1253  nSave++;
1254  int type(0);
1255  if (t_eMipDR < 1.0) {
1256  if (t_hmaxNearP < eIsolate2_) {
1257  ++nLoose;
1258  type = 1;
1259  }
1260  if (t_hmaxNearP < eIsolate1_) {
1261  ++nTight;
1262  type = 2;
1263  }
1264  }
1265  if (t_p > 40.0 && t_p <= 60.0 && t_selectTk) {
1266  t_ietaGood->emplace_back(t_ieta);
1267  t_trackType->emplace_back(type);
1268  }
1269 #ifdef EDM_ML_DEBUG
1270  for (unsigned int k = 0; k < t_trgbits->size(); k++) {
1271  edm::LogVerbatim("HcalIsoTrack") << "trigger bit is = " << t_trgbits->at(k);
1272  }
1273 #endif
1274  }
1275  }
1276  }
1277  }
1278  std::array<int, 3> i3{{nSave, nLoose, nTight}};
1279  return i3;
1280 }

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

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

◆ newId()

DetId HcalIsoTrkAnalyzer::newId ( const DetId &  id)
private

Definition at line 1347 of file HcalIsoTrkAnalyzer.cc.

1347  {
1348  HcalDetId hid(id);
1349  if (hep17_ && ((hid.iphi() < 63) || (hid.iphi() > 66) || (hid.zside() < 0)))
1350  return id;
1351  for (unsigned int k = 0; k < oldID_.size(); ++k) {
1352  if ((hid.subdetId() == oldDet_[k]) && (hid.ietaAbs() == oldEta_[k]) && (hid.depth() == oldDepth_[k])) {
1353  return static_cast<DetId>(HcalDetId(hid.subdet(), hid.ieta(), hid.iphi(), newDepth_[k]));
1354  }
1355  }
1356  return id;
1357 }

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

Referenced by fillTree().

◆ notaMuon()

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

Definition at line 1418 of file HcalIsoTrkAnalyzer.cc.

1418  {
1419  bool flag(true);
1420  for (reco::MuonCollection::const_iterator recMuon = muonh->begin(); recMuon != muonh->end(); ++recMuon) {
1421  if (recMuon->innerTrack().isNonnull()) {
1422  const reco::Track* pTrack = (recMuon->innerTrack()).get();
1423  bool mediumMuon = (((recMuon->isPFMuon()) && (recMuon->isGlobalMuon() || recMuon->isTrackerMuon())) &&
1424  (recMuon->innerTrack()->validFraction() > 0.49));
1425  if (mediumMuon) {
1426  double chiGlobal = ((recMuon->globalTrack().isNonnull()) ? recMuon->globalTrack()->normalizedChi2() : 999);
1427  bool goodGlob = (recMuon->isGlobalMuon() && chiGlobal < 3 &&
1428  recMuon->combinedQuality().chi2LocalPosition < 12 && recMuon->combinedQuality().trkKink < 20);
1429  mediumMuon = muon::segmentCompatibility(*recMuon) > (goodGlob ? 0.303 : 0.451);
1430  }
1431  if (mediumMuon) {
1432  double dR = reco::deltaR(pTrack->eta(), pTrack->phi(), pTrack0->eta(), pTrack0->phi());
1433  if (dR < 0.1) {
1434  flag = false;
1435  break;
1436  }
1437  }
1438  }
1439  }
1440  return flag;
1441 }

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

Referenced by fillTree().

◆ rhoh()

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

Definition at line 1302 of file HcalIsoTrkAnalyzer.cc.

1302  {
1303  std::vector<double> sumPFNallSMDQH2;
1304  sumPFNallSMDQH2.reserve(phibins_.size() * etabins_.size());
1305 
1306  for (auto eta : etabins_) {
1307  for (auto phi : phibins_) {
1308  double hadder = 0;
1309  for (const auto& pf_it : (*tower)) {
1310  if (fabs(eta - pf_it.eta()) > etahalfdist_)
1311  continue;
1312  if (fabs(reco::deltaPhi(phi, pf_it.phi())) > phihalfdist_)
1313  continue;
1314  hadder += pf_it.hadEt();
1315  }
1316  sumPFNallSMDQH2.emplace_back(hadder);
1317  }
1318  }
1319 
1320  double evt_smdq(0);
1321  std::sort(sumPFNallSMDQH2.begin(), sumPFNallSMDQH2.end());
1322  if (sumPFNallSMDQH2.size() % 2)
1323  evt_smdq = sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 1) / 2];
1324  else
1325  evt_smdq = (sumPFNallSMDQH2[sumPFNallSMDQH2.size() / 2] + sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 2) / 2]) / 2.;
1326  double rhoh = evt_smdq / (etadist_ * phidist_);
1327 #ifdef EDM_ML_DEBUG
1328  edm::LogVerbatim("HcalIsoTrack") << "Rho " << evt_smdq << ":" << rhoh;
1329 #endif
1330  return rhoh;
1331 }

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

Referenced by fillTree().

◆ storeEnergy()

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

Definition at line 1359 of file HcalIsoTrkAnalyzer.cc.

1365  {
1366  double ehcal(0);
1367  if (unCorrect_) {
1368  for (unsigned int k = 0; k < ids.size(); ++k) {
1369  double corr = (respCorrs->getValues(ids[k]))->getValue();
1370  if (corr != 0)
1371  edet[k] /= corr;
1372  ehcal += edet[k];
1373  }
1374  } else {
1375  for (const auto& en : edet)
1376  ehcal += en;
1377  }
1378  if ((std::abs(ehcal - eHcal) > 0.001) && (!unCorrect_))
1379  edm::LogWarning("HcalIsoTrack") << "Check inconsistent energies: " << indx << " " << eHcal << ":" << ehcal
1380  << " from " << ids.size() << " cells";
1381  eHcal = hcalScale_ * ehcal;
1382 
1383  if (collapseDepth_) {
1384  std::map<HcalDetId, double> hitMap;
1385  for (unsigned int k = 0; k < ids.size(); ++k) {
1386  HcalDetId id = hdc_->mergedDepthDetId(HcalDetId(ids[k]));
1387  auto itr = hitMap.find(id);
1388  if (itr == hitMap.end()) {
1389  hitMap[id] = edet[k];
1390  } else {
1391  (itr->second) += edet[k];
1392  }
1393  }
1394  detIds->reserve(hitMap.size());
1395  hitEnergies->reserve(hitMap.size());
1396  for (const auto& hit : hitMap) {
1397  detIds->emplace_back(hit.first.rawId());
1398  hitEnergies->emplace_back(hit.second);
1399  }
1400  } else {
1401  detIds->reserve(ids.size());
1402  hitEnergies->reserve(ids.size());
1403  for (unsigned int k = 0; k < ids.size(); ++k) {
1404  detIds->emplace_back(ids[k].rawId());
1405  hitEnergies->emplace_back(edet[k]);
1406  }
1407  }
1408 #ifdef EDM_ML_DEBUG
1409  edm::LogVerbatim("HcalIsoTrack") << "Input to storeEnergy with " << ids.size() << " cells";
1410  for (unsigned int k = 0; k < ids.size(); ++k)
1411  edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId(ids[k]) << " E " << edet[k];
1412  edm::LogVerbatim("HcalIsoTrack") << "Output of storeEnergy with " << detIds->size() << " cells and Etot " << eHcal;
1413  for (unsigned int k = 0; k < detIds->size(); ++k)
1414  edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId((*detIds)[k]) << " E " << (*hitEnergies)[k];
1415 #endif
1416 }

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

Referenced by fillTree().

◆ trackP()

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

Definition at line 1286 of file HcalIsoTrkAnalyzer.cc.

1287  {
1288  double pmom = -1.0;
1289  if (genParticles.isValid()) {
1290  double mindR(999.9);
1291  for (const auto& p : (*genParticles)) {
1292  double dR = reco::deltaR(pTrack->eta(), pTrack->phi(), p.momentum().Eta(), p.momentum().Phi());
1293  if (dR < mindR) {
1294  mindR = dR;
1295  pmom = p.momentum().R();
1296  }
1297  }
1298  }
1299  return pmom;
1300 }

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

Referenced by fillTree().

Member Data Documentation

◆ a_charIsoR_

double HcalIsoTrkAnalyzer::a_charIsoR_
private

Definition at line 161 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_coneR1_

double HcalIsoTrkAnalyzer::a_coneR1_
private

Definition at line 161 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_coneR2_

double HcalIsoTrkAnalyzer::a_coneR2_
private

Definition at line 161 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_coneR_

const double HcalIsoTrkAnalyzer::a_coneR_
private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR2_

const double HcalIsoTrkAnalyzer::a_mipR2_
private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR3_

const double HcalIsoTrkAnalyzer::a_mipR3_
private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR4_

const double HcalIsoTrkAnalyzer::a_mipR4_
private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR5_

const double HcalIsoTrkAnalyzer::a_mipR5_
private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ a_mipR_

const double HcalIsoTrkAnalyzer::a_mipR_
private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ collapseDepth_

const bool HcalIsoTrkAnalyzer::collapseDepth_
private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

◆ dataType_

const int HcalIsoTrkAnalyzer::dataType_
private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ eEcalMax_

const double HcalIsoTrkAnalyzer::eEcalMax_
private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ eIsolate1_

const double HcalIsoTrkAnalyzer::eIsolate1_
private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ eIsolate2_

const double HcalIsoTrkAnalyzer::eIsolate2_
private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ etabins_

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

Definition at line 163 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ etadist_

double HcalIsoTrkAnalyzer::etadist_
private

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ etahalfdist_

double HcalIsoTrkAnalyzer::etahalfdist_
private

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ fs

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

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob().

◆ hcalScale_

const double HcalIsoTrkAnalyzer::hcalScale_
private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

◆ hdc_

const HcalDDDRecConstants* HcalIsoTrkAnalyzer::hdc_
private

Definition at line 162 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and storeEnergy().

◆ hep17_

const bool HcalIsoTrkAnalyzer::hep17_
private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ hitEthrEB_

const double HcalIsoTrkAnalyzer::hitEthrEB_
private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE0_

const double HcalIsoTrkAnalyzer::hitEthrEE0_
private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE1_

const double HcalIsoTrkAnalyzer::hitEthrEE1_
private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE2_

const double HcalIsoTrkAnalyzer::hitEthrEE2_
private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEE3_

const double HcalIsoTrkAnalyzer::hitEthrEE3_
private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEEHi_

const double HcalIsoTrkAnalyzer::hitEthrEEHi_
private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hitEthrEELo_

const double HcalIsoTrkAnalyzer::hitEthrEELo_
private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

◆ hltConfig_

HLTConfigProvider HcalIsoTrkAnalyzer::hltConfig_
private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginRun().

◆ ignoreTrigger_

const bool HcalIsoTrkAnalyzer::ignoreTrigger_
private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l1Filter_

const std::string HcalIsoTrkAnalyzer::l1Filter_
private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l1GtUtils_

l1t::L1TGlobalUtil* HcalIsoTrkAnalyzer::l1GtUtils_
private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l1TrigName_

const std::string HcalIsoTrkAnalyzer::l1TrigName_
private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l2Filter_

const std::string HcalIsoTrkAnalyzer::l2Filter_
private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ l3Filter_

const std::string HcalIsoTrkAnalyzer::l3Filter_
private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelEB_

const std::string HcalIsoTrkAnalyzer::labelEB_
private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelEE_

const std::string HcalIsoTrkAnalyzer::labelEE_
private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelGenTrack_

const std::string HcalIsoTrkAnalyzer::labelGenTrack_
private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelHBHE_

const std::string HcalIsoTrkAnalyzer::labelHBHE_
private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ labelRecVtx_

const std::string HcalIsoTrkAnalyzer::labelRecVtx_
private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ labelTower_

const std::string HcalIsoTrkAnalyzer::labelTower_
private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ maxRestrictionP_

const double HcalIsoTrkAnalyzer::maxRestrictionP_
private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ mode_

const int HcalIsoTrkAnalyzer::mode_
private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ newDepth_

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

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ nHigh_

unsigned int HcalIsoTrkAnalyzer::nHigh_
private

Definition at line 160 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree().

◆ nLow_

unsigned int HcalIsoTrkAnalyzer::nLow_
private

Definition at line 160 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree().

◆ nRun_

unsigned int HcalIsoTrkAnalyzer::nRun_
private

Definition at line 160 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and endRun().

◆ oldDepth_

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

Definition at line 164 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ oldDet_

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

Definition at line 164 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ oldEta_

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

Definition at line 164 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

◆ oldID_

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

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

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

◆ phibins_

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

Definition at line 163 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ phidist_

double HcalIsoTrkAnalyzer::phidist_
private

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ phihalfdist_

double HcalIsoTrkAnalyzer::phihalfdist_
private

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

◆ prescaleHigh_

const int HcalIsoTrkAnalyzer::prescaleHigh_
private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ prescaleLow_

const int HcalIsoTrkAnalyzer::prescaleLow_
private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ processName_

const std::string HcalIsoTrkAnalyzer::processName_
private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkAnalyzer().

◆ pTrackHigh_

const double HcalIsoTrkAnalyzer::pTrackHigh_
private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ pTrackLow_

const double HcalIsoTrkAnalyzer::pTrackLow_
private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ pTrackMin_

const double HcalIsoTrkAnalyzer::pTrackMin_
private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ selectionParameter_

spr::trackSelectionParameters HcalIsoTrkAnalyzer::selectionParameter_
private

Definition at line 139 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ slopeRestrictionP_

const double HcalIsoTrkAnalyzer::slopeRestrictionP_
private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

◆ t_allvertex

int HcalIsoTrkAnalyzer::t_allvertex
private

Definition at line 208 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_DataType

int HcalIsoTrkAnalyzer::t_DataType
private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_DetIds

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

Definition at line 203 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_DetIds1

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

Definition at line 203 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_DetIds3

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

Definition at line 203 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eAnnular

double HcalIsoTrkAnalyzer::t_eAnnular
private

Definition at line 199 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal

double HcalIsoTrkAnalyzer::t_eHcal
private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal10

double HcalIsoTrkAnalyzer::t_eHcal10
private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eHcal30

double HcalIsoTrkAnalyzer::t_eHcal30
private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_emaxNearP

double HcalIsoTrkAnalyzer::t_emaxNearP
private

Definition at line 199 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR

double HcalIsoTrkAnalyzer::t_eMipDR
private

Definition at line 197 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR2

double HcalIsoTrkAnalyzer::t_eMipDR2
private

Definition at line 197 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR3

double HcalIsoTrkAnalyzer::t_eMipDR3
private

Definition at line 197 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR4

double HcalIsoTrkAnalyzer::t_eMipDR4
private

Definition at line 197 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_eMipDR5

double HcalIsoTrkAnalyzer::t_eMipDR5
private

Definition at line 198 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Event

int HcalIsoTrkAnalyzer::t_Event
private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_EventNo

unsigned int HcalIsoTrkAnalyzer::t_EventNo
private

Definition at line 190 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_EventWeight

double HcalIsoTrkAnalyzer::t_EventWeight
private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_gentrackP

double HcalIsoTrkAnalyzer::t_gentrackP
private

Definition at line 198 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_goodPV

int HcalIsoTrkAnalyzer::t_goodPV
private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_hAnnular

double HcalIsoTrkAnalyzer::t_hAnnular
private

Definition at line 199 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies

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

Definition at line 204 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies1

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

Definition at line 204 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_HitEnergies3

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

Definition at line 204 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_hltbits

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

Definition at line 205 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_hmaxNearP

double HcalIsoTrkAnalyzer::t_hmaxNearP
private

Definition at line 198 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_ieta

int HcalIsoTrkAnalyzer::t_ieta
private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_ietaAll

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

Definition at line 209 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_ietaGood

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

Definition at line 209 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_iphi

int HcalIsoTrkAnalyzer::t_iphi
private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_L1Bit

bool HcalIsoTrkAnalyzer::t_L1Bit
private

Definition at line 206 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1eta

double HcalIsoTrkAnalyzer::t_l1eta
private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1phi

double HcalIsoTrkAnalyzer::t_l1phi
private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l1pt

double HcalIsoTrkAnalyzer::t_l1pt
private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3eta

double HcalIsoTrkAnalyzer::t_l3eta
private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3phi

double HcalIsoTrkAnalyzer::t_l3phi
private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_l3pt

double HcalIsoTrkAnalyzer::t_l3pt
private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_mindR1

double HcalIsoTrkAnalyzer::t_mindR1
private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_mindR2

double HcalIsoTrkAnalyzer::t_mindR2
private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_nTrk

int HcalIsoTrkAnalyzer::t_nTrk
private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_nVtx

int HcalIsoTrkAnalyzer::t_nVtx
private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_p

double HcalIsoTrkAnalyzer::t_p
private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_phi

double HcalIsoTrkAnalyzer::t_phi
private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_pt

double HcalIsoTrkAnalyzer::t_pt
private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyFlag

bool HcalIsoTrkAnalyzer::t_qltyFlag
private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyMissFlag

bool HcalIsoTrkAnalyzer::t_qltyMissFlag
private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_qltyPVFlag

bool HcalIsoTrkAnalyzer::t_qltyPVFlag
private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_rhoh

double HcalIsoTrkAnalyzer::t_rhoh
private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Run

int HcalIsoTrkAnalyzer::t_Run
private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_RunNo

unsigned int HcalIsoTrkAnalyzer::t_RunNo
private

Definition at line 190 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_selectTk

bool HcalIsoTrkAnalyzer::t_selectTk
private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

◆ t_Tracks

int HcalIsoTrkAnalyzer::t_Tracks
private

Definition at line 207 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksLoose

int HcalIsoTrkAnalyzer::t_TracksLoose
private

Definition at line 208 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksProp

int HcalIsoTrkAnalyzer::t_TracksProp
private

Definition at line 207 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksSaved

int HcalIsoTrkAnalyzer::t_TracksSaved
private

Definition at line 207 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TracksTight

int HcalIsoTrkAnalyzer::t_TracksTight
private

Definition at line 208 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_trackType

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

Definition at line 209 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_trgbits

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

Definition at line 205 of file HcalIsoTrkAnalyzer.cc.

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

◆ t_TrigPass

bool HcalIsoTrkAnalyzer::t_TrigPass
private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ t_TrigPassSel

bool HcalIsoTrkAnalyzer::t_TrigPassSel
private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ theTrackQuality_

const std::string HcalIsoTrkAnalyzer::theTrackQuality_
private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ theTriggerResultsLabel_

const edm::InputTag HcalIsoTrkAnalyzer::theTriggerResultsLabel_
private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

◆ tok_alg_

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

Definition at line 177 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_bFieldH_

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

Definition at line 181 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_bs_

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

Definition at line 171 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_cala_

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

Definition at line 175 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_caloTopology_

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

Definition at line 185 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_ddrec_

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

Definition at line 180 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkAnalyzer().

◆ tok_EB_

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

Definition at line 172 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_ecalChStatus_

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

Definition at line 182 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_EE_

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

Definition at line 173 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_ew_

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

Definition at line 176 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_genTrack_

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

Definition at line 169 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_geom_

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

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_hbhe_

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

Definition at line 174 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_htopo_

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

Definition at line 186 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_Muon_

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

Definition at line 178 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_parts_

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

Definition at line 168 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_recVtx_

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

Definition at line 170 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_resp_

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

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_sevlv_

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

Definition at line 183 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_trigEvt_

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

Definition at line 166 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tok_trigRes_

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

Definition at line 167 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ tree

TTree* HcalIsoTrkAnalyzer::tree
private

Definition at line 189 of file HcalIsoTrkAnalyzer.cc.

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

◆ tree2

TTree * HcalIsoTrkAnalyzer::tree2
private

Definition at line 189 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ triggerEvent_

const edm::InputTag HcalIsoTrkAnalyzer::triggerEvent_
private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ trigNames_

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

Definition at line 138 of file HcalIsoTrkAnalyzer.cc.

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

◆ unCorrect_

const bool HcalIsoTrkAnalyzer::unCorrect_
private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

◆ useL1Trigger_

const bool HcalIsoTrkAnalyzer::useL1Trigger_
private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

◆ useRaw_

const int HcalIsoTrkAnalyzer::useRaw_
private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

HcalIsoTrkAnalyzer::phidist_
double phidist_
Definition: HcalIsoTrkAnalyzer.cc:165
spr::trackSelectionParameters::maxDxyPV
double maxDxyPV
Definition: TrackSelection.h:22
HcalDDDRecConstants::mergedDepthDetId
HcalDetId mergedDepthDetId(const HcalDetId &id) const
Definition: HcalDDDRecConstants.cc:647
HcalIsoTrkAnalyzer::tok_cala_
edm::EDGetTokenT< CaloTowerCollection > tok_cala_
Definition: HcalIsoTrkAnalyzer.cc:175
EcalSeverityLevelAlgo
Definition: EcalSeverityLevelAlgo.h:33
jetSelector_cfi.nLoose
nLoose
Definition: jetSelector_cfi.py:11
DDAxes::y
HcalIsoTrkAnalyzer::tok_bFieldH_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > tok_bFieldH_
Definition: HcalIsoTrkAnalyzer.cc:181
HcalIsoTrkAnalyzer::t_p
double t_p
Definition: HcalIsoTrkAnalyzer.cc:193
testProducerWithPsetDescEmpty_cfi.i3
i3
Definition: testProducerWithPsetDescEmpty_cfi.py:47
HcalRespCorrs
Definition: HcalRespCorrs.h:17
HcalIsoTrkAnalyzer::tok_EB_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: HcalIsoTrkAnalyzer.cc:172
mps_fire.i
i
Definition: mps_fire.py:428
HcalIsoTrkAnalyzer::t_TracksTight
int t_TracksTight
Definition: HcalIsoTrkAnalyzer.cc:208
HcalIsoTrkAnalyzer::prescaleLow_
const int prescaleLow_
Definition: HcalIsoTrkAnalyzer.cc:148
HcalIsoTrkAnalyzer::l1TrigName_
const std::string l1TrigName_
Definition: HcalIsoTrkAnalyzer.cc:157
HcalIsoTrkAnalyzer::hitEthrEE3_
const double hitEthrEE3_
Definition: HcalIsoTrkAnalyzer.cc:153
HcalIsoTrkAnalyzer::labelRecVtx_
const std::string labelRecVtx_
Definition: HcalIsoTrkAnalyzer.cc:156
LogMessageMonitor_cff.modules
modules
Definition: LogMessageMonitor_cff.py:7
HcalIsoTrkAnalyzer::t_trgbits
std::vector< bool > * t_trgbits
Definition: HcalIsoTrkAnalyzer.cc:205
HcalIsoTrkAnalyzer::t_nVtx
int t_nVtx
Definition: HcalIsoTrkAnalyzer.cc:192
genParticles2HepMC_cfi.genParticles
genParticles
Definition: genParticles2HepMC_cfi.py:4
HcalIsoTrkAnalyzer::t_DetIds3
std::vector< unsigned int > * t_DetIds3
Definition: HcalIsoTrkAnalyzer.cc:203
HcalIsoTrkAnalyzer::fs
edm::Service< TFileService > fs
Definition: HcalIsoTrkAnalyzer.cc:136
HcalIsoTrkAnalyzer::t_eHcal10
double t_eHcal10
Definition: HcalIsoTrkAnalyzer.cc:200
HcalIsoTrkAnalyzer::a_mipR4_
const double a_mipR4_
Definition: HcalIsoTrkAnalyzer.cc:144
dqmiodumpmetadata.n
n
Definition: dqmiodumpmetadata.py:28
HcalCondObjectContainerBase::setTopo
void setTopo(const HcalTopology *topo)
Definition: HcalCondObjectContainerBase.cc:17
edm::Handle::product
T const * product() const
Definition: Handle.h:70
HcalIsoTrkAnalyzer::tok_ecalChStatus_
edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > tok_ecalChStatus_
Definition: HcalIsoTrkAnalyzer.cc:182
vec1
std::vector< double > vec1
Definition: HCALResponse.h:15
spr::trackSelectionParameters::maxDzPV
double maxDzPV
Definition: TrackSelection.h:22
HcalIsoTrkAnalyzer::t_nTrk
int t_nTrk
Definition: HcalIsoTrkAnalyzer.cc:192
HcalDetId::iphi
constexpr int iphi() const
get the cell iphi
Definition: HcalDetId.h:157
spr::trackSelectionParameters::maxChi2
double maxChi2
Definition: TrackSelection.h:22
muon::segmentCompatibility
float segmentCompatibility(const reco::Muon &muon, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
Definition: MuonSelectors.cc:61
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:631
HcalIsoTrkAnalyzer::selectionParameter_
spr::trackSelectionParameters selectionParameter_
Definition: HcalIsoTrkAnalyzer.cc:139
HcalIsoTrkAnalyzer::oldID_
const std::vector< int > oldID_
Definition: HcalIsoTrkAnalyzer.cc:158
HcalIsoTrkAnalyzer::tok_bs_
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
Definition: HcalIsoTrkAnalyzer.cc:171
CaloGeometry::getPosition
GlobalPoint getPosition(const DetId &id) const
Get the position of a given detector id.
Definition: CaloGeometry.cc:50
reco::deltaPhi
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:26
HcalIsoTrkAnalyzer::a_coneR1_
double a_coneR1_
Definition: HcalIsoTrkAnalyzer.cc:161
HcalIsoTrkAnalyzer::nLow_
unsigned int nLow_
Definition: HcalIsoTrkAnalyzer.cc:160
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
tree
Definition: tree.py:1
HcalIsoTrkAnalyzer::t_mindR1
double t_mindR1
Definition: HcalIsoTrkAnalyzer.cc:196
pos
Definition: PixelAliasList.h:18
HcalTopology
Definition: HcalTopology.h:26
HcalIsoTrkAnalyzer::nHigh_
unsigned int nHigh_
Definition: HcalIsoTrkAnalyzer.cc:160
HcalIsoTrkAnalyzer::tok_htopo_
edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tok_htopo_
Definition: HcalIsoTrkAnalyzer.cc:186
HcalIsoTrkAnalyzer::t_DetIds1
std::vector< unsigned int > * t_DetIds1
Definition: HcalIsoTrkAnalyzer.cc:203
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89285
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
hgcalTowerProducer_cfi.tower
tower
Definition: hgcalTowerProducer_cfi.py:4
HcalIsoTrkAnalyzer::theTriggerResultsLabel_
const edm::InputTag theTriggerResultsLabel_
Definition: HcalIsoTrkAnalyzer.cc:155
triggerResults
static const std::string triggerResults
Definition: EdmProvDump.cc:45
HcalIsoTrkAnalyzer::eIsolate2_
const double eIsolate2_
Definition: HcalIsoTrkAnalyzer.cc:146
HcalIsoTrkAnalyzer::a_coneR_
const double a_coneR_
Definition: HcalIsoTrkAnalyzer.cc:143
HcalIsoTrkAnalyzer::tok_trigEvt_
edm::EDGetTokenT< trigger::TriggerEvent > tok_trigEvt_
Definition: HcalIsoTrkAnalyzer.cc:166
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:150
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:47
HcalIsoTrkAnalyzer::t_qltyFlag
bool t_qltyFlag
Definition: HcalIsoTrkAnalyzer.cc:201
HcalIsoTrkAnalyzer::tok_genTrack_
edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
Definition: HcalIsoTrkAnalyzer.cc:169
HcalIsoTrkAnalyzer::a_mipR_
const double a_mipR_
Definition: HcalIsoTrkAnalyzer.cc:143
reco::TrackBase::px
double px() const
x coordinate of momentum vector
Definition: TrackBase.h:640
HcalIsoTrkAnalyzer::t_TracksLoose
int t_TracksLoose
Definition: HcalIsoTrkAnalyzer.cc:208
HcalIsoTrkAnalyzer::tok_ew_
edm::EDGetTokenT< GenEventInfoProduct > tok_ew_
Definition: HcalIsoTrkAnalyzer.cc:176
HcalIsoTrkAnalyzer::t_hltbits
std::vector< bool > * t_hltbits
Definition: HcalIsoTrkAnalyzer.cc:205
spr::goodTrack
bool goodTrack(const reco::Track *pTrack, math::XYZPoint leadPV, trackSelectionParameters parameters, bool debug=false)
Definition: TrackSelection.cc:10
l1t::UseEventSetupIn::Event
HcalIsoTrkAnalyzer::t_DataType
int t_DataType
Definition: HcalIsoTrkAnalyzer.cc:191
HcalIsoTrkAnalyzer::t_qltyPVFlag
bool t_qltyPVFlag
Definition: HcalIsoTrkAnalyzer.cc:202
DDAxes::x
HcalIsoTrkAnalyzer::t_l3phi
double t_l3phi
Definition: HcalIsoTrkAnalyzer.cc:195
HcalIsoTrkAnalyzer::etabins_
std::vector< double > etabins_
Definition: HcalIsoTrkAnalyzer.cc:163
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
HcalIsoTrkAnalyzer::tok_ddrec_
edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord > tok_ddrec_
Definition: HcalIsoTrkAnalyzer.cc:180
EcalCondObjectContainer
Definition: EcalCondObjectContainer.h:13
vec2
std::vector< vec1 > vec2
Definition: HCALResponse.h:16
BeamSpotPI::nTracks
Definition: BeamSpotPayloadInspectorHelper.h:42
HcalIsoTrkAnalyzer::useRaw_
const int useRaw_
Definition: HcalIsoTrkAnalyzer.cc:149
convertSQLiteXML.ok
bool ok
Definition: convertSQLiteXML.py:98
HcalIsoTrkAnalyzer::t_HitEnergies
std::vector< double > * t_HitEnergies
Definition: HcalIsoTrkAnalyzer.cc:204
edm::Handle< reco::GenParticleCollection >
HcalIsoTrkAnalyzer::trigNames_
const std::vector< std::string > trigNames_
Definition: HcalIsoTrkAnalyzer.cc:138
HcalIsoTrkAnalyzer::l3Filter_
const std::string l3Filter_
Definition: HcalIsoTrkAnalyzer.cc:142
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
CaloTopology
Definition: CaloTopology.h:19
HcalIsoTrkAnalyzer::t_hmaxNearP
double t_hmaxNearP
Definition: HcalIsoTrkAnalyzer.cc:198
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
EcalBarrel
Definition: EcalSubdetector.h:10
HcalIsoTrkAnalyzer::tok_caloTopology_
edm::ESGetToken< CaloTopology, CaloTopologyRecord > tok_caloTopology_
Definition: HcalIsoTrkAnalyzer.cc:185
HcalIsoTrkAnalyzer::mode_
const int mode_
Definition: HcalIsoTrkAnalyzer.cc:149
HcalIsoTrkAnalyzer::pTrackHigh_
const double pTrackHigh_
Definition: HcalIsoTrkAnalyzer.cc:147
HcalIsoTrkAnalyzer::rhoh
double rhoh(const edm::Handle< CaloTowerCollection > &)
Definition: HcalIsoTrkAnalyzer.cc:1302
HcalIsoTrkAnalyzer::oldDet_
std::vector< int > oldDet_
Definition: HcalIsoTrkAnalyzer.cc:164
HcalIsoTrkAnalyzer::t_Event
int t_Event
Definition: HcalIsoTrkAnalyzer.cc:191
HcalIsoTrkAnalyzer::t_l1eta
double t_l1eta
Definition: HcalIsoTrkAnalyzer.cc:194
HcalIsoTrkAnalyzer::t_ietaGood
std::vector< int > * t_ietaGood
Definition: HcalIsoTrkAnalyzer.cc:209
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:637
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
HLTConfigProvider::moduleLabels
const std::vector< std::string > & moduleLabels(unsigned int trigger) const
label(s) of module(s) on a trigger path
Definition: HLTConfigProvider.h:76
HcalIsoTrkAnalyzer::hitEthrEEHi_
const double hitEthrEEHi_
Definition: HcalIsoTrkAnalyzer.cc:154
DetId
Definition: DetId.h:17
HcalIsoTrkAnalyzer::t_hAnnular
double t_hAnnular
Definition: HcalIsoTrkAnalyzer.cc:199
HcalIsoTrkAnalyzer::tok_EE_
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: HcalIsoTrkAnalyzer.cc:173
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
CaloGeometry
Definition: CaloGeometry.h:21
HcalIsoTrkAnalyzer::t_eMipDR5
double t_eMipDR5
Definition: HcalIsoTrkAnalyzer.cc:198
HcalIsoTrkAnalyzer::t_phi
double t_phi
Definition: HcalIsoTrkAnalyzer.cc:193
HcalIsoTrkAnalyzer::dR
double dR(math::XYZTLorentzVector &, math::XYZTLorentzVector &)
Definition: HcalIsoTrkAnalyzer.cc:1282
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)
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
reco::TrackBase::py
double py() const
y coordinate of momentum vector
Definition: TrackBase.h:643
l1t::L1TGlobalUtil::decisionsFinal
const std::vector< std::pair< std::string, bool > > & decisionsFinal()
Definition: L1TGlobalUtil.h:156
names
const std::string names[nVars_]
Definition: PhotonIDValueMapProducer.cc:124
fileCollector.done
done
Definition: fileCollector.py:123
L1TEGammaOffline_cfi.triggerNames
triggerNames
Definition: L1TEGammaOffline_cfi.py:40
alignCSCRings.corr
dictionary corr
Definition: alignCSCRings.py:124
PVValHelper::eta
Definition: PVValidationHelpers.h:70
l1t::L1TGlobalUtil
Definition: L1TGlobalUtil.h:44
HcalIsoTrkAnalyzer::t_l1phi
double t_l1phi
Definition: HcalIsoTrkAnalyzer.cc:194
HcalIsoTrkAnalyzer::phibins_
std::vector< double > phibins_
Definition: HcalIsoTrkAnalyzer.cc:163
HcalIsoTrkAnalyzer::newId
DetId newId(const DetId &)
Definition: HcalIsoTrkAnalyzer.cc:1347
HcalIsoTrkAnalyzer::labelEE_
const std::string labelEE_
Definition: HcalIsoTrkAnalyzer.cc:157
HcalIsoTrkAnalyzer::t_TrigPass
bool t_TrigPass
Definition: HcalIsoTrkAnalyzer.cc:202
DDAxes::z
reco::Track
Definition: Track.h:27
HcalIsoTrkAnalyzer::t_L1Bit
bool t_L1Bit
Definition: HcalIsoTrkAnalyzer.cc:206
spr::trackSelectionParameters::minPt
double minPt
Definition: TrackSelection.h:20
HcalIsoTrkAnalyzer::phihalfdist_
double phihalfdist_
Definition: HcalIsoTrkAnalyzer.cc:165
spr::trackSelectionParameters
Definition: TrackSelection.h:12
spr::trackSelectionParameters::minOuterHit
int minOuterHit
Definition: TrackSelection.h:23
HcalIsoTrkAnalyzer::t_selectTk
bool t_selectTk
Definition: HcalIsoTrkAnalyzer.cc:201
HcalIsoTrkAnalyzer::tok_trigRes_
edm::EDGetTokenT< edm::TriggerResults > tok_trigRes_
Definition: HcalIsoTrkAnalyzer.cc:167
dqmdumpme.k
k
Definition: dqmdumpme.py:60
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:59
Point3DBase< float, GlobalTag >
HcalIsoTrkAnalyzer::t_eMipDR2
double t_eMipDR2
Definition: HcalIsoTrkAnalyzer.cc:197
HcalIsoTrkAnalyzer::t_goodPV
int t_goodPV
Definition: HcalIsoTrkAnalyzer.cc:192
HcalIsoTrkAnalyzer::hitEthrEELo_
const double hitEthrEELo_
Definition: HcalIsoTrkAnalyzer.cc:154
spr::chargeIsolationCone
double chargeIsolationCone(unsigned int trkIndex, std::vector< spr::propagatedTrackDirection > &trkDirs, double dR, int &nNearTRKs, bool debug=false)
Definition: ChargeIsolation.cc:288
HcalIsoTrkAnalyzer::useL1Trigger_
const bool useL1Trigger_
Definition: HcalIsoTrkAnalyzer.cc:150
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:649
trigger::TriggerObject
Single trigger physics object (e.g., an isolated muon)
Definition: TriggerObject.h:21
HcalIsoTrkAnalyzer::t_Run
int t_Run
Definition: HcalIsoTrkAnalyzer.cc:191
spr::trackSelectionParameters::maxDpOverP
double maxDpOverP
Definition: TrackSelection.h:22
spr::trackSelectionParameters::minLayerCrossed
int minLayerCrossed
Definition: TrackSelection.h:23
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HcalIsoTrkAnalyzer::t_eMipDR3
double t_eMipDR3
Definition: HcalIsoTrkAnalyzer.cc:197
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:63
HcalIsoTrkAnalyzer::t_RunNo
unsigned int t_RunNo
Definition: HcalIsoTrkAnalyzer.cc:190
HcalIsoTrkAnalyzer::tok_alg_
edm::EDGetTokenT< BXVector< GlobalAlgBlk > > tok_alg_
Definition: HcalIsoTrkAnalyzer.cc:177
HcalIsoTrkAnalyzer::a_mipR2_
const double a_mipR2_
Definition: HcalIsoTrkAnalyzer.cc:143
HcalDetId::ieta
constexpr int ieta() const
get the cell ieta
Definition: HcalDetId.h:155
HcalIsoTrkAnalyzer::a_charIsoR_
double a_charIsoR_
Definition: HcalIsoTrkAnalyzer.cc:161
HcalIsoTrkAnalyzer::t_eHcal
double t_eHcal
Definition: HcalIsoTrkAnalyzer.cc:200
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
HcalIsoTrkAnalyzer::t_qltyMissFlag
bool t_qltyMissFlag
Definition: HcalIsoTrkAnalyzer.cc:201
HcalIsoTrkAnalyzer::newDepth_
const std::vector< int > newDepth_
Definition: HcalIsoTrkAnalyzer.cc:158
HcalIsoTrkAnalyzer::tok_resp_
edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_resp_
Definition: HcalIsoTrkAnalyzer.cc:187
HcalIsoTrkAnalyzer::eEcalMax_
const double eEcalMax_
Definition: HcalIsoTrkAnalyzer.cc:144
type
type
Definition: SiPixelVCal_PayloadInspector.cc:37
jetUpdater_cfi.sort
sort
Definition: jetUpdater_cfi.py:29
HcalIsoTrkAnalyzer::t_trackType
std::vector< int > * t_trackType
Definition: HcalIsoTrkAnalyzer.cc:209
trigger::TriggerObjectCollection
std::vector< TriggerObject > TriggerObjectCollection
collection of trigger physics objects (e.g., all isolated muons)
Definition: TriggerObject.h:75
HcalIsoTrkAnalyzer::t_allvertex
int t_allvertex
Definition: HcalIsoTrkAnalyzer.cc:208
gainCalibHelper::gainCalibPI::type
type
Definition: SiPixelGainCalibHelper.h:40
spr::trackSelectionParameters::maxOutMiss
int maxOutMiss
Definition: TrackSelection.h:24
HcalIsoTrkAnalyzer::t_emaxNearP
double t_emaxNearP
Definition: HcalIsoTrkAnalyzer.cc:199
HcalIsoTrkAnalyzer::hdc_
const HcalDDDRecConstants * hdc_
Definition: HcalIsoTrkAnalyzer.cc:162
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:652
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
HcalIsoTrkAnalyzer::oldEta_
std::vector< int > oldEta_
Definition: HcalIsoTrkAnalyzer.cc:164
l1t::L1TGlobalUtil::retrieveL1
void retrieveL1(const edm::Event &iEvent, const edm::EventSetup &evSetup)
initialize the class (mainly reserve)
Definition: L1TGlobalUtil.cc:66
HcalIsoTrkAnalyzer::t_eMipDR4
double t_eMipDR4
Definition: HcalIsoTrkAnalyzer.cc:197
HcalDetId
Definition: HcalDetId.h:12
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)
iEvent
int iEvent
Definition: GenABIO.cc:224
HcalIsoTrkAnalyzer::labelGenTrack_
const std::string labelGenTrack_
Definition: HcalIsoTrkAnalyzer.cc:156
HcalIsoTrkAnalyzer::tok_Muon_
edm::EDGetTokenT< reco::MuonCollection > tok_Muon_
Definition: HcalIsoTrkAnalyzer.cc:178
HcalIsoTrkAnalyzer::tok_geom_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord > tok_geom_
Definition: HcalIsoTrkAnalyzer.cc:184
HcalIsoTrkAnalyzer::a_mipR5_
const double a_mipR5_
Definition: HcalIsoTrkAnalyzer.cc:144
HcalIsoTrkAnalyzer::hcalScale_
const double hcalScale_
Definition: HcalIsoTrkAnalyzer.cc:146
HcalIsoTrkAnalyzer::tok_sevlv_
edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd > tok_sevlv_
Definition: HcalIsoTrkAnalyzer.cc:183
trigger::Keys
std::vector< size_type > Keys
Definition: TriggerTypeDefs.h:19
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:49
HcalIsoTrkAnalyzer::eThreshold
double eThreshold(const DetId &id, const CaloGeometry *geo) const
Definition: HcalIsoTrkAnalyzer.cc:1333
HcalIsoTrkAnalyzer::t_l3eta
double t_l3eta
Definition: HcalIsoTrkAnalyzer.cc:195
HcalIsoTrkAnalyzer::slopeRestrictionP_
const double slopeRestrictionP_
Definition: HcalIsoTrkAnalyzer.cc:145
HcalIsoTrkAnalyzer::maxRestrictionP_
const double maxRestrictionP_
Definition: HcalIsoTrkAnalyzer.cc:145
trigger::TriggerEvent
The single EDProduct to be saved for each event (AOD case)
Definition: TriggerEvent.h:25
HcalIsoTrkAnalyzer::t_DetIds
std::vector< unsigned int > * t_DetIds
Definition: HcalIsoTrkAnalyzer.cc:203
get
#define get
PDWG_DiPhoton_SD_cff.triggerEvent
triggerEvent
Definition: PDWG_DiPhoton_SD_cff.py:39
spr::trackSelectionParameters::minQuality
reco::TrackBase::TrackQuality minQuality
Definition: TrackSelection.h:21
HcalIsoTrkAnalyzer::t_ieta
int t_ieta
Definition: HcalIsoTrkAnalyzer.cc:191
ClusterTask_cfi.trig
trig
Definition: ClusterTask_cfi.py:398
HcalIsoTrkAnalyzer::a_coneR2_
double a_coneR2_
Definition: HcalIsoTrkAnalyzer.cc:161
HcalIsoTrkAnalyzer::t_iphi
int t_iphi
Definition: HcalIsoTrkAnalyzer.cc:191
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:73
photonIsolationHIProducer_cfi.hbhe
hbhe
Definition: photonIsolationHIProducer_cfi.py:8
HcalIsoTrkAnalyzer::notaMuon
bool notaMuon(const reco::Track *pTrack0, const edm::Handle< reco::MuonCollection > &muonh)
Definition: HcalIsoTrkAnalyzer.cc:1418
HcalIsoTrkAnalyzer::nRun_
unsigned int nRun_
Definition: HcalIsoTrkAnalyzer.cc:160
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::ignoreTrigger_
const bool ignoreTrigger_
Definition: HcalIsoTrkAnalyzer.cc:150
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
HcalIsoTrkAnalyzer::hitEthrEE0_
const double hitEthrEE0_
Definition: HcalIsoTrkAnalyzer.cc:152
DDAxes::phi
HcalIsoTrkAnalyzer::hitEthrEE2_
const double hitEthrEE2_
Definition: HcalIsoTrkAnalyzer.cc:153
HcalIsoTrkAnalyzer::t_l1pt
double t_l1pt
Definition: HcalIsoTrkAnalyzer.cc:194
HcalIsoTrkAnalyzer::theTrackQuality_
const std::string theTrackQuality_
Definition: HcalIsoTrkAnalyzer.cc:140
HcalIsoTrkAnalyzer::t_HitEnergies1
std::vector< double > * t_HitEnergies1
Definition: HcalIsoTrkAnalyzer.cc:204
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, const edm::Handle< reco::MuonCollection > &muonh)
Definition: HcalIsoTrkAnalyzer.cc:912
HcalIsoTrkAnalyzer::t_TracksProp
int t_TracksProp
Definition: HcalIsoTrkAnalyzer.cc:207
submitPVResolutionJobs.desc
string desc
Definition: submitPVResolutionJobs.py:251
HcalIsoTrkAnalyzer::pTrackLow_
const double pTrackLow_
Definition: HcalIsoTrkAnalyzer.cc:147
HcalIsoTrkAnalyzer::tok_parts_
edm::EDGetTokenT< reco::GenParticleCollection > tok_parts_
Definition: HcalIsoTrkAnalyzer.cc:168
HcalIsoTrkAnalyzer::t_ietaAll
std::vector< int > * t_ietaAll
Definition: HcalIsoTrkAnalyzer.cc:209
genParticles2HepMC_cfi.genEventInfo
genEventInfo
Definition: genParticles2HepMC_cfi.py:6
HLTConfigProvider::init
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d'tor
Definition: HLTConfigProvider.cc:36
HcalIsoTrkAnalyzer::etahalfdist_
double etahalfdist_
Definition: HcalIsoTrkAnalyzer.cc:165
Calorimetry_cff.bField
bField
Definition: Calorimetry_cff.py:284
HcalIsoTrkAnalyzer::t_mindR2
double t_mindR2
Definition: HcalIsoTrkAnalyzer.cc:196
HcalIsoTrkAnalyzer::hitEthrEB_
const double hitEthrEB_
Definition: HcalIsoTrkAnalyzer.cc:152
HcalIsoTrkAnalyzer::labelTower_
const std::string labelTower_
Definition: HcalIsoTrkAnalyzer.cc:157
HcalIsoTrkAnalyzer::triggerEvent_
const edm::InputTag triggerEvent_
Definition: HcalIsoTrkAnalyzer.cc:155
HcalIsoTrkAnalyzer::hitEthrEE1_
const double hitEthrEE1_
Definition: HcalIsoTrkAnalyzer.cc:152
HcalIsoTrkAnalyzer::t_eAnnular
double t_eAnnular
Definition: HcalIsoTrkAnalyzer.cc:199
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:29
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
spr::trackSelectionParameters::maxInMiss
int maxInMiss
Definition: TrackSelection.h:24
HcalIsoTrkAnalyzer::l1GtUtils_
l1t::L1TGlobalUtil * l1GtUtils_
Definition: HcalIsoTrkAnalyzer.cc:135
HcalIsoTrkAnalyzer::trackP
double trackP(const reco::Track *, const edm::Handle< reco::GenParticleCollection > &)
Definition: HcalIsoTrkAnalyzer.cc:1286
HcalIsoTrkAnalyzer::t_EventWeight
double t_EventWeight
Definition: HcalIsoTrkAnalyzer.cc:193
flavorHistoryFilter_cfi.dr
dr
Definition: flavorHistoryFilter_cfi.py:37
TFileService::kSharedResource
static const std::string kSharedResource
Definition: TFileService.h:76
edm::TriggerNames
Definition: TriggerNames.h:55
HcalIsoTrkAnalyzer::a_mipR3_
const double a_mipR3_
Definition: HcalIsoTrkAnalyzer.cc:143
HcalIsoTrkAnalyzer::unCorrect_
const bool unCorrect_
Definition: HcalIsoTrkAnalyzer.cc:151
HcalIsoTrkAnalyzer::t_gentrackP
double t_gentrackP
Definition: HcalIsoTrkAnalyzer.cc:198
HcalIsoTrkAnalyzer::l1Filter_
const std::string l1Filter_
Definition: HcalIsoTrkAnalyzer.cc:141
reco::deltaR
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
HcalIsoTrkAnalyzer::collapseDepth_
const bool collapseDepth_
Definition: HcalIsoTrkAnalyzer.cc:151
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_TrigPassSel
bool t_TrigPassSel
Definition: HcalIsoTrkAnalyzer.cc:202
HcalIsoTrkAnalyzer::hltConfig_
HLTConfigProvider hltConfig_
Definition: HcalIsoTrkAnalyzer.cc:137
HcalIsoTrkAnalyzer::tree2
TTree * tree2
Definition: HcalIsoTrkAnalyzer.cc:189
HcalIsoTrkAnalyzer::t_EventNo
unsigned int t_EventNo
Definition: HcalIsoTrkAnalyzer.cc:190
HcalIsoTrkAnalyzer::oldDepth_
std::vector< int > oldDepth_
Definition: HcalIsoTrkAnalyzer.cc:164
HcalIsoTrkAnalyzer::pTrackMin_
const double pTrackMin_
Definition: HcalIsoTrkAnalyzer.cc:144
HcalIsoTrkAnalyzer::etadist_
double etadist_
Definition: HcalIsoTrkAnalyzer.cc:165
reco::TrackBase::pz
double pz() const
z coordinate of momentum vector
Definition: TrackBase.h:646
HcalIsoTrkAnalyzer::t_rhoh
double t_rhoh
Definition: HcalIsoTrkAnalyzer.cc:200
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
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_TracksSaved
int t_TracksSaved
Definition: HcalIsoTrkAnalyzer.cc:207
dummy
Definition: DummySelector.h:38
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
HcalIsoTrkAnalyzer::t_pt
double t_pt
Definition: HcalIsoTrkAnalyzer.cc:193
HcalIsoTrkAnalyzer::t_l3pt
double t_l3pt
Definition: HcalIsoTrkAnalyzer.cc:195
HcalIsoTrkAnalyzer::tok_recVtx_
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
Definition: HcalIsoTrkAnalyzer.cc:170
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
HcalIsoTrkAnalyzer::t_eHcal30
double t_eHcal30
Definition: HcalIsoTrkAnalyzer.cc:200
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)
HcalIsoTrkAnalyzer::eIsolate1_
const double eIsolate1_
Definition: HcalIsoTrkAnalyzer.cc:146
MagneticField
Definition: MagneticField.h:19
edm::Log
Definition: MessageLogger.h:70
HcalIsoTrkAnalyzer::t_eMipDR
double t_eMipDR
Definition: HcalIsoTrkAnalyzer.cc:197
HcalIsoTrkAnalyzer::t_HitEnergies3
std::vector< double > * t_HitEnergies3
Definition: HcalIsoTrkAnalyzer.cc:204
HcalIsoTrkAnalyzer::processName_
const std::string processName_
Definition: HcalIsoTrkAnalyzer.cc:141
edm::InputTag
Definition: InputTag.h:15
label
const char * label
Definition: PFTauDecayModeTools.cc:11
HcalIsoTrkAnalyzer::labelHBHE_
const std::string labelHBHE_
Definition: HcalIsoTrkAnalyzer.cc:157
HcalIsoTrkAnalyzer::labelEB_
const std::string labelEB_
Definition: HcalIsoTrkAnalyzer.cc:156
hit
Definition: SiStripHitEffFromCalibTree.cc:88
TFileService::make
T * make(const Args &... args) const
make new ROOT object
Definition: TFileService.h:64
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:1359
RemoveAddSevLevel.flag
flag
Definition: RemoveAddSevLevel.py:116
HcalIsoTrkAnalyzer::prescaleHigh_
const int prescaleHigh_
Definition: HcalIsoTrkAnalyzer.cc:148
HcalIsoTrkAnalyzer::t_Tracks
int t_Tracks
Definition: HcalIsoTrkAnalyzer.cc:207
ValidationMatrix.hlt
hlt
Definition: ValidationMatrix.py:459
HcalIsoTrkAnalyzer::dataType_
const int dataType_
Definition: HcalIsoTrkAnalyzer.cc:149
HcalIsoTrkAnalyzer::hep17_
const bool hep17_
Definition: HcalIsoTrkAnalyzer.cc:159
HcalIsoTrkAnalyzer::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: HcalIsoTrkAnalyzer.cc:174
HcalIsoTrkAnalyzer::l2Filter_
const std::string l2Filter_
Definition: HcalIsoTrkAnalyzer.cc:142
Generated for CMSSW Reference Manual by doxygen 1.8.16