CMS 3D CMS Logo
HcalIsoTrkAnalyzer.cc
Go to the documentation of this file.
1 // system include files
2 #include <atomic>
3 #include <cmath>
4 #include <memory>
5 #include <string>
6 #include <vector>
7 
8 // Root objects
9 #include "TROOT.h"
10 #include "TSystem.h"
11 #include "TFile.h"
12 #include "TProfile.h"
13 #include "TDirectory.h"
14 #include "TTree.h"
15 #include "TLorentzVector.h"
16 #include "TInterpreter.h"
17 
18 #include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
19 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
20 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
21 
22 //Tracks
23 #include "DataFormats/TrackReco/interface/Track.h"
24 #include "DataFormats/TrackReco/interface/TrackFwd.h"
25 #include "DataFormats/TrackReco/interface/HitPattern.h"
26 #include "DataFormats/TrackReco/interface/TrackBase.h"
27 // Vertices
28 #include "DataFormats/VertexReco/interface/VertexFwd.h"
29 #include "DataFormats/VertexReco/interface/Vertex.h"
30 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
31 
32 //Triggers
33 #include "DataFormats/Common/interface/TriggerResults.h"
34 #include "DataFormats/HLTReco/interface/TriggerEvent.h"
35 #include "DataFormats/HLTReco/interface/TriggerObject.h"
36 #include "DataFormats/Luminosity/interface/LumiDetails.h"
37 #include "DataFormats/Math/interface/LorentzVector.h"
38 #include "DataFormats/Math/interface/deltaR.h"
39 
40 #include "L1Trigger/L1TGlobal/interface/L1TGlobalUtil.h"
41 #include "DataFormats/L1TGlobal/interface/GlobalAlgBlk.h"
42 #include "DataFormats/L1TGlobal/interface/GlobalExtBlk.h"
43 #include "DataFormats/L1Trigger/interface/BXVector.h"
44 
45 #include "CondFormats/DataRecord/interface/EcalChannelStatusRcd.h"
46 #include "CondFormats/DataRecord/interface/HcalRespCorrsRcd.h"
47 #include "CondFormats/HcalObjects/interface/HcalRespCorrs.h"
48 
49 //Generator information
50 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
51 #include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"
52 #include "SimDataFormats/GeneratorProducts/interface/GenEventInfoProduct.h"
53 
54 #include "FWCore/Framework/interface/Frameworkfwd.h"
55 #include "FWCore/Framework/interface/one/EDAnalyzer.h"
56 
57 #include "FWCore/Framework/interface/Event.h"
58 #include "FWCore/Framework/interface/MakerMacros.h"
59 #include "FWCore/ServiceRegistry/interface/Service.h"
60 #include "FWCore/ParameterSet/interface/ParameterSet.h"
61 #include "FWCore/Common/interface/TriggerNames.h"
62 #include "CommonTools/UtilAlgos/interface/TFileService.h"
63 #include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
64 
65 #include "Calibration/IsolatedParticles/interface/CaloPropagateTrack.h"
66 #include "Calibration/IsolatedParticles/interface/ChargeIsolation.h"
67 #include "Calibration/IsolatedParticles/interface/eCone.h"
68 #include "Calibration/IsolatedParticles/interface/eECALMatrix.h"
69 #include "Calibration/IsolatedParticles/interface/eHCALMatrix.h"
70 #include "Calibration/IsolatedParticles/interface/TrackSelection.h"
71 
72 #include "MagneticField/Engine/interface/MagneticField.h"
73 #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
74 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
75 #include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
76 #include "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"
77 #include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
78 #include "Geometry/CaloTopology/interface/HcalTopology.h"
79 #include "Geometry/CaloTopology/interface/CaloTopology.h"
80 #include "Geometry/HcalCommonData/interface/HcalDDDRecConstants.h"
81 #include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"
82 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
83 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"
84 
85 //#define EDM_ML_DEBUG
86 
87 class HcalIsoTrkAnalyzer : public edm::one::EDAnalyzer<edm::one::WatchRuns, edm::one::SharedResources> {
88 public:
89  explicit HcalIsoTrkAnalyzer(edm::ParameterSet const&);
90  ~HcalIsoTrkAnalyzer() override {}
91 
92  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
93 
94 private:
95  void analyze(edm::Event const&, edm::EventSetup const&) override;
96  void beginJob() override;
97  void beginRun(edm::Run const&, edm::EventSetup const&) override;
98  void endRun(edm::Run const&, edm::EventSetup const&) override;
99 
100  std::array<int, 3> fillTree(std::vector<math::XYZTLorentzVector>& vecL1,
101  std::vector<math::XYZTLorentzVector>& vecL3,
102  math::XYZPoint& leadPV,
103  std::vector<spr::propagatedTrackDirection>& trkCaloDirections,
104  std::vector<spr::propagatedTrackID>& trkCaloDets,
105  const CaloGeometry* geo,
106  const CaloTopology* topo,
107  const HcalTopology* theHBHETopology,
108  const EcalChannelStatus* theEcalChStatus,
109  const EcalSeverityLevelAlgo* theEcalSevlv,
110  edm::Handle<EcalRecHitCollection>& barrelRecHitsHandle,
111  edm::Handle<EcalRecHitCollection>& endcapRecHitsHandle,
112  edm::Handle<HBHERecHitCollection>& hbhe,
113  edm::Handle<CaloTowerCollection>& towerHandle,
114  edm::Handle<reco::GenParticleCollection>& genParticles,
115  const HcalRespCorrs* respCorrs);
116  double dR(math::XYZTLorentzVector&, math::XYZTLorentzVector&);
117  double trackP(const reco::Track*, const edm::Handle<reco::GenParticleCollection>&);
118  double rhoh(const edm::Handle<CaloTowerCollection>&);
119  double eThreshold(const DetId& id, const CaloGeometry* geo) const;
120  DetId newId(const DetId&);
121  void storeEnergy(int indx,
122  const HcalRespCorrs* respCorrs,
123  const std::vector<DetId>& ids,
124  std::vector<double>& edet,
125  double& eHcal,
126  std::vector<unsigned int>* detIds,
127  std::vector<double>* hitEnergies);
128 
129  l1t::L1TGlobalUtil* l1GtUtils_;
130  edm::Service<TFileService> fs;
131  HLTConfigProvider hltConfig_;
132  const std::vector<std::string> trigNames_;
133  spr::trackSelectionParameters selectionParameter_;
134  const std::string theTrackQuality_;
135  const std::string processName_, l1Filter_;
136  const std::string l2Filter_, l3Filter_;
137  const double a_coneR_, a_mipR_, a_mipR2_, a_mipR3_;
138  const double a_mipR4_, a_mipR5_, pTrackMin_, eEcalMax_;
139  const double maxRestrictionP_, slopeRestrictionP_;
140  const double hcalScale_, eIsolate1_, eIsolate2_;
141  const double pTrackLow_, pTrackHigh_;
142  const int prescaleLow_, prescaleHigh_;
143  const int useRaw_, dataType_, mode_;
144  const bool ignoreTrigger_, useL1Trigger_;
145  const bool unCorrect_, collapseDepth_;
146  const double hitEthrEB_, hitEthrEE0_, hitEthrEE1_;
147  const double hitEthrEE2_, hitEthrEE3_;
148  const double hitEthrEELo_, hitEthrEEHi_;
149  const edm::InputTag triggerEvent_, theTriggerResultsLabel_;
150  const std::string labelGenTrack_, labelRecVtx_, labelEB_;
151  const std::string labelEE_, labelHBHE_, labelTower_, l1TrigName_;
152  const std::vector<int> oldID_, newDepth_;
153  const bool hep17_;
154  unsigned int nRun_, nLow_, nHigh_;
155  double a_charIsoR_, a_coneR1_, a_coneR2_;
156  const HcalDDDRecConstants* hdc_;
157  std::vector<double> etabins_, phibins_;
158  std::vector<int> oldDet_, oldEta_, oldDepth_;
159  double etadist_, phidist_, etahalfdist_, phihalfdist_;
160  edm::EDGetTokenT<trigger::TriggerEvent> tok_trigEvt_;
161  edm::EDGetTokenT<edm::TriggerResults> tok_trigRes_;
162  edm::EDGetTokenT<reco::GenParticleCollection> tok_parts_;
163  edm::EDGetTokenT<reco::TrackCollection> tok_genTrack_;
164  edm::EDGetTokenT<reco::VertexCollection> tok_recVtx_;
165  edm::EDGetTokenT<reco::BeamSpot> tok_bs_;
166  edm::EDGetTokenT<EcalRecHitCollection> tok_EB_;
167  edm::EDGetTokenT<EcalRecHitCollection> tok_EE_;
168  edm::EDGetTokenT<HBHERecHitCollection> tok_hbhe_;
169  edm::EDGetTokenT<CaloTowerCollection> tok_cala_;
170  edm::EDGetTokenT<GenEventInfoProduct> tok_ew_;
171  edm::EDGetTokenT<BXVector<GlobalAlgBlk>> tok_alg_;
172 
173  edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> tok_ddrec_;
174  edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> tok_bFieldH_;
175  edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> tok_ecalChStatus_;
176  edm::ESGetToken<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd> tok_sevlv_;
177  edm::ESGetToken<CaloGeometry, CaloGeometryRecord> tok_geom_;
178  edm::ESGetToken<CaloTopology, CaloTopologyRecord> tok_caloTopology_;
179  edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> tok_htopo_;
180  edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> tok_resp_;
181 
182  TTree *tree, *tree2;
183  unsigned int t_RunNo, t_EventNo;
184  int t_Run, t_Event, t_DataType, t_ieta, t_iphi;
185  int t_goodPV, t_nVtx, t_nTrk;
186  double t_EventWeight, t_p, t_pt, t_phi;
187  double t_l1pt, t_l1eta, t_l1phi;
188  double t_l3pt, t_l3eta, t_l3phi;
189  double t_mindR1, t_mindR2;
190  double t_eMipDR, t_eMipDR2, t_eMipDR3, t_eMipDR4;
191  double t_eMipDR5, t_hmaxNearP, t_gentrackP;
192  double t_emaxNearP, t_eAnnular, t_hAnnular;
193  double t_eHcal, t_eHcal10, t_eHcal30, t_rhoh;
194  bool t_selectTk, t_qltyFlag, t_qltyMissFlag;
195  bool t_qltyPVFlag, t_TrigPass, t_TrigPassSel;
196  std::vector<unsigned int>*t_DetIds, *t_DetIds1, *t_DetIds3;
197  std::vector<double>*t_HitEnergies, *t_HitEnergies1, *t_HitEnergies3;
198  std::vector<bool>*t_trgbits, *t_hltbits;
199  bool t_L1Bit;
200  int t_Tracks, t_TracksProp, t_TracksSaved;
201  int t_TracksLoose, t_TracksTight, t_allvertex;
202  std::vector<int>*t_ietaAll, *t_ietaGood, *t_trackType;
203 };
204 
205 HcalIsoTrkAnalyzer::HcalIsoTrkAnalyzer(const edm::ParameterSet& iConfig)
206  : trigNames_(iConfig.getParameter<std::vector<std::string>>("triggers")),
207  theTrackQuality_(iConfig.getParameter<std::string>("trackQuality")),
208  processName_(iConfig.getParameter<std::string>("processName")),
209  l1Filter_(iConfig.getParameter<std::string>("l1Filter")),
210  l2Filter_(iConfig.getParameter<std::string>("l2Filter")),
211  l3Filter_(iConfig.getParameter<std::string>("l3Filter")),
212  a_coneR_(iConfig.getParameter<double>("coneRadius")),
213  a_mipR_(iConfig.getParameter<double>("coneRadiusMIP")),
214  a_mipR2_(iConfig.getParameter<double>("coneRadiusMIP2")),
215  a_mipR3_(iConfig.getParameter<double>("coneRadiusMIP3")),
216  a_mipR4_(iConfig.getParameter<double>("coneRadiusMIP4")),
217  a_mipR5_(iConfig.getParameter<double>("coneRadiusMIP5")),
218  pTrackMin_(iConfig.getParameter<double>("minimumTrackP")),
219  eEcalMax_(iConfig.getParameter<double>("maximumEcalEnergy")),
220  maxRestrictionP_(iConfig.getParameter<double>("maxTrackP")),
221  slopeRestrictionP_(iConfig.getParameter<double>("slopeTrackP")),
222  hcalScale_(iConfig.getUntrackedParameter<double>("hHcalScale", 1.0)),
223  eIsolate1_(iConfig.getParameter<double>("isolationEnergyTight")),
224  eIsolate2_(iConfig.getParameter<double>("isolationEnergyLoose")),
225  pTrackLow_(iConfig.getParameter<double>("momentumLow")),
226  pTrackHigh_(iConfig.getParameter<double>("momentumHigh")),
227  prescaleLow_(iConfig.getParameter<int>("prescaleLow")),
228  prescaleHigh_(iConfig.getParameter<int>("prescaleHigh")),
229  useRaw_(iConfig.getUntrackedParameter<int>("useRaw", 0)),
230  dataType_(iConfig.getUntrackedParameter<int>("dataType", 0)),
231  mode_(iConfig.getUntrackedParameter<int>("outMode", 11)),
232  ignoreTrigger_(iConfig.getUntrackedParameter<bool>("ignoreTriggers", false)),
233  useL1Trigger_(iConfig.getUntrackedParameter<bool>("useL1Trigger", false)),
234  unCorrect_(iConfig.getUntrackedParameter<bool>("unCorrect", false)),
235  collapseDepth_(iConfig.getUntrackedParameter<bool>("collapseDepth", false)),
236  hitEthrEB_(iConfig.getParameter<double>("EBHitEnergyThreshold")),
237  hitEthrEE0_(iConfig.getParameter<double>("EEHitEnergyThreshold0")),
238  hitEthrEE1_(iConfig.getParameter<double>("EEHitEnergyThreshold1")),
239  hitEthrEE2_(iConfig.getParameter<double>("EEHitEnergyThreshold2")),
240  hitEthrEE3_(iConfig.getParameter<double>("EEHitEnergyThreshold3")),
241  hitEthrEELo_(iConfig.getParameter<double>("EEHitEnergyThresholdLow")),
242  hitEthrEEHi_(iConfig.getParameter<double>("EEHitEnergyThresholdHigh")),
243  triggerEvent_(iConfig.getParameter<edm::InputTag>("labelTriggerEvent")),
244  theTriggerResultsLabel_(iConfig.getParameter<edm::InputTag>("labelTriggerResult")),
245  labelGenTrack_(iConfig.getParameter<std::string>("labelTrack")),
246  labelRecVtx_(iConfig.getParameter<std::string>("labelVertex")),
247  labelEB_(iConfig.getParameter<std::string>("labelEBRecHit")),
248  labelEE_(iConfig.getParameter<std::string>("labelEERecHit")),
249  labelHBHE_(iConfig.getParameter<std::string>("labelHBHERecHit")),
250  labelTower_(iConfig.getParameter<std::string>("labelCaloTower")),
251  l1TrigName_(iConfig.getUntrackedParameter<std::string>("l1TrigName", "L1_SingleJet60")),
252  oldID_(iConfig.getUntrackedParameter<std::vector<int>>("oldID")),
253  newDepth_(iConfig.getUntrackedParameter<std::vector<int>>("newDepth")),
254  hep17_(iConfig.getUntrackedParameter<bool>("hep17")),
255  nRun_(0),
256  nLow_(0),
257  nHigh_(0),
258  hdc_(nullptr) {
259  usesResource(TFileService::kSharedResource);
260 
261  //now do whatever initialization is needed
262  const double isolationRadius(28.9), innerR(10.0), outerR(30.0);
263  reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality_);
264  selectionParameter_.minPt = iConfig.getParameter<double>("minTrackPt");
265  ;
266  selectionParameter_.minQuality = trackQuality_;
267  selectionParameter_.maxDxyPV = iConfig.getParameter<double>("maxDxyPV");
268  selectionParameter_.maxDzPV = iConfig.getParameter<double>("maxDzPV");
269  selectionParameter_.maxChi2 = iConfig.getParameter<double>("maxChi2");
270  selectionParameter_.maxDpOverP = iConfig.getParameter<double>("maxDpOverP");
271  selectionParameter_.minOuterHit = iConfig.getParameter<int>("minOuterHit");
272  selectionParameter_.minLayerCrossed = iConfig.getParameter<int>("minLayerCrossed");
273  selectionParameter_.maxInMiss = iConfig.getParameter<int>("maxInMiss");
274  selectionParameter_.maxOutMiss = iConfig.getParameter<int>("maxOutMiss");
275  a_charIsoR_ = a_coneR_ + isolationRadius;
276  a_coneR1_ = a_coneR_ + innerR;
277  a_coneR2_ = a_coneR_ + outerR;
278  // Different isolation cuts are described in DN-2016/029
279  // Tight cut uses 2 GeV; Loose cut uses 10 GeV
280  // Eta dependent cut uses (maxRestrictionP_ * exp(|ieta|*log(2.5)/18))
281  // with the factor for exponential slopeRestrictionP_ = log(2.5)/18
282  // maxRestrictionP_ = 8 GeV as came from a study
283  std::string labelBS = iConfig.getParameter<std::string>("labelBeamSpot");
284  std::string modnam = iConfig.getUntrackedParameter<std::string>("moduleName", "");
285  std::string prdnam = iConfig.getUntrackedParameter<std::string>("producerName", "");
286  edm::InputTag algTag = iConfig.getParameter<edm::InputTag>("algInputTag");
287  edm::InputTag extTag = iConfig.getParameter<edm::InputTag>("extInputTag");
288  for (unsigned int k = 0; k < oldID_.size(); ++k) {
289  oldDet_.emplace_back((oldID_[k] / 10000) % 10);
290  oldEta_.emplace_back((oldID_[k] / 100) % 100);
291  oldDepth_.emplace_back(oldID_[k] % 100);
292  }
293 
294  l1GtUtils_ = new l1t::L1TGlobalUtil(iConfig, consumesCollector(), *this, algTag, extTag, l1t::UseEventSetupIn::Event);
295  // define tokens for access
296  tok_trigEvt_ = consumes<trigger::TriggerEvent>(triggerEvent_);
297  tok_trigRes_ = consumes<edm::TriggerResults>(theTriggerResultsLabel_);
298  tok_bs_ = consumes<reco::BeamSpot>(labelBS);
299  tok_genTrack_ = consumes<reco::TrackCollection>(labelGenTrack_);
300  tok_ew_ = consumes<GenEventInfoProduct>(edm::InputTag("generator"));
301  tok_parts_ = consumes<reco::GenParticleCollection>(edm::InputTag("genParticles"));
302  tok_cala_ = consumes<CaloTowerCollection>(labelTower_);
303  tok_alg_ = consumes<BXVector<GlobalAlgBlk>>(algTag);
304 
305  if (modnam.empty()) {
306  tok_recVtx_ = consumes<reco::VertexCollection>(labelRecVtx_);
307  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEB_));
308  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEE_));
309  tok_hbhe_ = consumes<HBHERecHitCollection>(labelHBHE_);
310  edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
311  << labelBS << " " << labelRecVtx_ << " " << labelGenTrack_ << " "
312  << edm::InputTag("ecalRecHit", labelEB_) << " "
313  << edm::InputTag("ecalRecHit", labelEE_) << " " << labelHBHE_ << " "
314  << labelTower_;
315  } else {
316  tok_recVtx_ = consumes<reco::VertexCollection>(edm::InputTag(modnam, labelRecVtx_, prdnam));
317  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag(modnam, labelEB_, prdnam));
318  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag(modnam, labelEE_, prdnam));
319  tok_hbhe_ = consumes<HBHERecHitCollection>(edm::InputTag(modnam, labelHBHE_, prdnam));
320  edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
321  << labelBS << " " << edm::InputTag(modnam, labelRecVtx_, prdnam) << " "
322  << labelGenTrack_ << " " << edm::InputTag(modnam, labelEB_, prdnam) << " "
323  << edm::InputTag(modnam, labelEE_, prdnam) << " "
324  << edm::InputTag(modnam, labelHBHE_, prdnam) << " " << labelTower_;
325  }
326 
327  tok_ddrec_ = esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>();
328  tok_bFieldH_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
329  tok_ecalChStatus_ = esConsumes<EcalChannelStatus, EcalChannelStatusRcd>();
330  tok_sevlv_ = esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>();
331  tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
332  tok_caloTopology_ = esConsumes<CaloTopology, CaloTopologyRecord>();
333  tok_htopo_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
334  tok_resp_ = esConsumes<HcalRespCorrs, HcalRespCorrsRcd>();
335 
336  edm::LogVerbatim("HcalIsoTrack")
337  << "Parameters read from config file \n"
338  << "\t minPt " << selectionParameter_.minPt << "\t theTrackQuality " << theTrackQuality_ << "\t minQuality "
339  << selectionParameter_.minQuality << "\t maxDxyPV " << selectionParameter_.maxDxyPV << "\t maxDzPV "
340  << selectionParameter_.maxDzPV << "\t maxChi2 " << selectionParameter_.maxChi2 << "\t maxDpOverP "
341  << selectionParameter_.maxDpOverP << "\t minOuterHit " << selectionParameter_.minOuterHit << "\t minLayerCrossed "
342  << selectionParameter_.minLayerCrossed << "\t maxInMiss " << selectionParameter_.maxInMiss << "\t maxOutMiss "
343  << selectionParameter_.maxOutMiss << "\t a_coneR " << a_coneR_ << ":" << a_coneR1_ << ":" << a_coneR2_
344  << "\t a_charIsoR " << a_charIsoR_ << "\t a_mipR " << a_mipR_ << "\t a_mipR2 " << a_mipR2_ << "\t a_mipR3 "
345  << a_mipR3_ << "\t a_mipR4 " << a_mipR4_ << "\t a_mipR5 " << a_mipR5_ << "\n pTrackMin_ " << pTrackMin_
346  << "\t eEcalMax_ " << eEcalMax_ << "\t maxRestrictionP_ " << maxRestrictionP_ << "\t slopeRestrictionP_ "
347  << slopeRestrictionP_ << "\t eIsolateStrong_ " << eIsolate1_ << "\t eIsolateSoft_ " << eIsolate2_
348  << "\t hcalScale_ " << hcalScale_ << "\n\t momentumLow_ " << pTrackLow_ << "\t prescaleLow_ " << prescaleLow_
349  << "\t momentumHigh_ " << pTrackHigh_ << "\t prescaleHigh_ " << prescaleHigh_ << "\n\t useRaw_ " << useRaw_
350  << "\t ignoreTrigger_ " << ignoreTrigger_ << "\n\t useL1Trigegr_ " << useL1Trigger_ << "\t dataType_ "
351  << dataType_ << "\t mode_ " << mode_ << "\t unCorrect_ " << unCorrect_ << "\t collapseDepth_ "
352  << collapseDepth_ << "\t L1TrigName_ " << l1TrigName_ << "\nThreshold for EB " << hitEthrEB_ << " EE "
353  << hitEthrEE0_ << ":" << hitEthrEE1_ << ":" << hitEthrEE2_ << ":" << hitEthrEE3_ << ":" << hitEthrEELo_ << ":"
354  << hitEthrEEHi_;
355  edm::LogVerbatim("HcalIsoTrack") << "Process " << processName_ << " L1Filter:" << l1Filter_
356  << " L2Filter:" << l2Filter_ << " L3Filter:" << l3Filter_;
357  for (unsigned int k = 0; k < trigNames_.size(); ++k) {
358  edm::LogVerbatim("HcalIsoTrack") << "Trigger[" << k << "] " << trigNames_[k];
359  }
360  edm::LogVerbatim("HcalIsoTrack") << oldID_.size() << " DetIDs to be corrected with HEP17 flag:" << hep17_;
361  for (unsigned int k = 0; k < oldID_.size(); ++k)
362  edm::LogVerbatim("HcalIsoTrack") << "[" << k << "] Det " << oldDet_[k] << " EtaAbs " << oldEta_[k] << " Depth "
363  << oldDepth_[k] << ":" << newDepth_[k];
364 
365  for (int i = 0; i < 10; i++)
366  phibins_.push_back(-M_PI + 0.1 * (2 * i + 1) * M_PI);
367  for (int i = 0; i < 8; ++i)
368  etabins_.push_back(-2.1 + 0.6 * i);
369  etadist_ = etabins_[1] - etabins_[0];
370  phidist_ = phibins_[1] - phibins_[0];
371  etahalfdist_ = 0.5 * etadist_;
372  phihalfdist_ = 0.5 * phidist_;
373  edm::LogVerbatim("HcalIsoTrack") << "EtaDist " << etadist_ << " " << etahalfdist_ << " PhiDist " << phidist_ << " "
374  << phihalfdist_;
375  unsigned int k1(0), k2(0);
376  for (auto phi : phibins_) {
377  edm::LogVerbatim("HcalIsoTrack") << "phibin_[" << k1 << "] " << phi;
378  ++k1;
379  }
380  for (auto eta : etabins_) {
381  edm::LogVerbatim("HcalIsoTrack") << "etabin_[" << k2 << "] " << eta;
382  ++k2;
383  }
384 }
385 
386 void HcalIsoTrkAnalyzer::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) {
387  t_Run = iEvent.id().run();
388  t_Event = iEvent.id().event();
389  t_DataType = dataType_;
390 #ifdef EDM_ML_DEBUG
391  edm::LogVerbatim("HcalIsoTrack") << "Run " << t_Run << " Event " << t_Event << " type " << t_DataType
392  << " Luminosity " << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing();
393 #endif
394  //Get magnetic field and ECAL channel status
395  const MagneticField* bField = &iSetup.getData(tok_bFieldH_);
396  const EcalChannelStatus* theEcalChStatus = &iSetup.getData(tok_ecalChStatus_);
397  const EcalSeverityLevelAlgo* theEcalSevlv = &iSetup.getData(tok_sevlv_);
398 
399  // get handles to calogeometry and calotopology
400  const CaloGeometry* geo = &iSetup.getData(tok_geom_);
401  const CaloTopology* caloTopology = &iSetup.getData(tok_caloTopology_);
402  const HcalTopology* theHBHETopology = &iSetup.getData(tok_htopo_);
403  const HcalRespCorrs* resp = &iSetup.getData(tok_resp_);
404  HcalRespCorrs* respCorrs = new HcalRespCorrs(*resp);
405  respCorrs->setTopo(theHBHETopology);
406 
407  //=== genParticle information
408  edm::Handle<reco::GenParticleCollection> genParticles;
409  iEvent.getByToken(tok_parts_, genParticles);
410 
411  bool okC(true);
412  //Get track collection
413  edm::Handle<reco::TrackCollection> trkCollection;
414  iEvent.getByToken(tok_genTrack_, trkCollection);
415  if (!trkCollection.isValid()) {
416  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelGenTrack_;
417  okC = false;
418  }
419 
420  //event weight for FLAT sample
421  t_EventWeight = 1.0;
422  edm::Handle<GenEventInfoProduct> genEventInfo;
423  iEvent.getByToken(tok_ew_, genEventInfo);
424  if (genEventInfo.isValid())
425  t_EventWeight = genEventInfo->weight();
426 
427  //Define the best vertex and the beamspot
428  edm::Handle<reco::VertexCollection> recVtxs;
429  iEvent.getByToken(tok_recVtx_, recVtxs);
430  edm::Handle<reco::BeamSpot> beamSpotH;
431  iEvent.getByToken(tok_bs_, beamSpotH);
432  math::XYZPoint leadPV(0, 0, 0);
433  t_goodPV = t_nVtx = 0;
434  if (recVtxs.isValid() && !(recVtxs->empty())) {
435  t_nVtx = recVtxs->size();
436  for (unsigned int k = 0; k < recVtxs->size(); ++k) {
437  if (!((*recVtxs)[k].isFake()) && ((*recVtxs)[k].ndof() > 4)) {
438  if (t_goodPV == 0)
439  leadPV = math::XYZPoint((*recVtxs)[k].x(), (*recVtxs)[k].y(), (*recVtxs)[k].z());
440  t_goodPV++;
441  }
442  }
443  }
444  if (t_goodPV == 0 && beamSpotH.isValid()) {
445  leadPV = beamSpotH->position();
446  }
447  t_allvertex = t_goodPV;
448 #ifdef EDM_ML_DEBUG
449  edm::LogVerbatim("HcalIsoTrack") << "Primary Vertex " << leadPV << " out of " << t_goodPV << " vertex";
450  if (beamSpotH.isValid()) {
451  edm::LogVerbatim("HcalIsoTrack") << " Beam Spot " << beamSpotH->position();
452  }
453 #endif
454  // RecHits
455  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
456  iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
457  if (!barrelRecHitsHandle.isValid()) {
458  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEB_;
459  okC = false;
460  }
461  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
462  iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
463  if (!endcapRecHitsHandle.isValid()) {
464  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEE_;
465  okC = false;
466  }
467  edm::Handle<HBHERecHitCollection> hbhe;
468  iEvent.getByToken(tok_hbhe_, hbhe);
469  if (!hbhe.isValid()) {
470  edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelHBHE_;
471  okC = false;
472  }
473  edm::Handle<CaloTowerCollection> caloTower;
474  iEvent.getByToken(tok_cala_, caloTower);
475 
476  //Propagate tracks to calorimeter surface)
477  std::vector<spr::propagatedTrackDirection> trkCaloDirections;
478  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDirections, false);
479  std::vector<spr::propagatedTrackID> trkCaloDets;
480  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDets, false);
481  std::vector<math::XYZTLorentzVector> vecL1, vecL3;
482  t_RunNo = iEvent.id().run();
483  t_EventNo = iEvent.id().event();
484  t_Tracks = trkCollection->size();
485  t_TracksProp = trkCaloDirections.size();
486  t_ietaAll->clear();
487  t_ietaGood->clear();
488  t_trackType->clear();
489  t_trgbits->clear();
490  t_hltbits->clear();
491 #ifdef EDM_ML_DEBUG
492  edm::LogVerbatim("HcalIsoTrack") << "# of propagated tracks " << t_TracksProp << " out of " << t_Tracks
493  << " with Trigger " << ignoreTrigger_;
494 #endif
495 
496  //Trigger
497  t_trgbits->assign(trigNames_.size(), false);
498  t_hltbits->assign(trigNames_.size(), false);
499  t_TracksSaved = t_TracksLoose = t_TracksTight = 0;
500  t_L1Bit = true;
501  t_TrigPass = false;
502 
503  //L1
504  l1GtUtils_->retrieveL1(iEvent, iSetup, tok_alg_);
505  const std::vector<std::pair<std::string, bool>>& finalDecisions = l1GtUtils_->decisionsFinal();
506  for (const auto& decision : finalDecisions) {
507  if (decision.first.find(l1TrigName_) != std::string::npos) {
508  t_L1Bit = decision.second;
509  break;
510  }
511  }
512 #ifdef EDM_ML_DEBUG
513  edm::LogVerbatim("HcalIsoTrack") << "Trigger Information for " << l1TrigName_ << " is " << t_L1Bit
514  << " from a list of " << finalDecisions.size() << " decisions";
515 #endif
516 
517  //HLT
518  edm::Handle<edm::TriggerResults> triggerResults;
519  iEvent.getByToken(tok_trigRes_, triggerResults);
520  if (triggerResults.isValid()) {
521  const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
522  const std::vector<std::string>& names = triggerNames.triggerNames();
523  if (!trigNames_.empty()) {
524  for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
525  int hlt = triggerResults->accept(iHLT);
526  for (unsigned int i = 0; i < trigNames_.size(); ++i) {
527  if (names[iHLT].find(trigNames_[i]) != std::string::npos) {
528  t_trgbits->at(i) = (hlt > 0);
529  t_hltbits->at(i) = (hlt > 0);
530  if (hlt > 0)
531  t_TrigPass = true;
532 #ifdef EDM_ML_DEBUG
533  edm::LogVerbatim("HcalIsoTrack")
534  << "This trigger " << names[iHLT] << " Flag " << hlt << ":" << t_trgbits->at(i);
535 #endif
536  }
537  }
538  }
539  }
540  }
541 #ifdef EDM_ML_DEBUG
542  edm::LogVerbatim("HcalIsoTrack") << "HLT Information shows " << t_TrigPass << ":" << trigNames_.empty() << ":" << okC;
543 #endif
544 
545  std::array<int, 3> ntksave{{0, 0, 0}};
546  if (ignoreTrigger_ || useL1Trigger_) {
547  t_l1pt = t_l1eta = t_l1phi = 0;
548  t_l3pt = t_l3eta = t_l3phi = 0;
549  if (ignoreTrigger_ || t_L1Bit)
550  ntksave = fillTree(vecL1,
551  vecL3,
552  leadPV,
553  trkCaloDirections,
554  trkCaloDets,
555  geo,
556  caloTopology,
557  theHBHETopology,
558  theEcalChStatus,
559  theEcalSevlv,
560  barrelRecHitsHandle,
561  endcapRecHitsHandle,
562  hbhe,
563  caloTower,
564  genParticles,
565  respCorrs);
566  t_TracksSaved = ntksave[0];
567  t_TracksLoose = ntksave[1];
568  t_TracksTight = ntksave[2];
569  } else {
570  trigger::TriggerEvent triggerEvent;
571  edm::Handle<trigger::TriggerEvent> triggerEventHandle;
572  iEvent.getByToken(tok_trigEvt_, triggerEventHandle);
573  if (!triggerEventHandle.isValid()) {
574  edm::LogWarning("HcalIsoTrack") << "Error! Can't get the product " << triggerEvent_.label();
575  } else if (okC) {
576  triggerEvent = *(triggerEventHandle.product());
577  const trigger::TriggerObjectCollection& TOC(triggerEvent.getObjects());
578  bool done(false);
579  if (triggerResults.isValid()) {
580  std::vector<std::string> modules;
581  const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
582  const std::vector<std::string>& names = triggerNames.triggerNames();
583  for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
584  bool ok = (t_TrigPass) || (trigNames_.empty());
585  if (ok) {
586  unsigned int triggerindx = hltConfig_.triggerIndex(names[iHLT]);
587  const std::vector<std::string>& moduleLabels(hltConfig_.moduleLabels(triggerindx));
588  std::vector<math::XYZTLorentzVector> vecL2;
589  vecL1.clear();
590  vecL3.clear();
591  //loop over all trigger filters in event (i.e. filters passed)
592  for (unsigned int ifilter = 0; ifilter < triggerEvent.sizeFilters(); ++ifilter) {
593  std::vector<int> Keys;
594  std::string label = triggerEvent.filterTag(ifilter).label();
595  //loop over keys to objects passing this filter
596  for (unsigned int imodule = 0; imodule < moduleLabels.size(); imodule++) {
597  if (label.find(moduleLabels[imodule]) != std::string::npos) {
598 #ifdef EDM_ML_DEBUG
599  edm::LogVerbatim("HcalIsoTrack") << "FilterName " << label;
600 #endif
601  for (unsigned int ifiltrKey = 0; ifiltrKey < triggerEvent.filterKeys(ifilter).size(); ++ifiltrKey) {
602  Keys.push_back(triggerEvent.filterKeys(ifilter)[ifiltrKey]);
603  const trigger::TriggerObject& TO(TOC[Keys[ifiltrKey]]);
604  math::XYZTLorentzVector v4(TO.px(), TO.py(), TO.pz(), TO.energy());
605  if (label.find(l2Filter_) != std::string::npos) {
606  vecL2.push_back(v4);
607  } else if (label.find(l3Filter_) != std::string::npos) {
608  vecL3.push_back(v4);
609  } else if ((label.find(l1Filter_) != std::string::npos) || (l1Filter_.empty())) {
610  vecL1.push_back(v4);
611  }
612 #ifdef EDM_ML_DEBUG
613  edm::LogVerbatim("HcalIsoTrack")
614  << "key " << ifiltrKey << " : pt " << TO.pt() << " eta " << TO.eta() << " phi " << TO.phi()
615  << " mass " << TO.mass() << " Id " << TO.id();
616 #endif
617  }
618 #ifdef EDM_ML_DEBUG
619  edm::LogVerbatim("HcalIsoTrack")
620  << "sizes " << vecL1.size() << ":" << vecL2.size() << ":" << vecL3.size();
621 #endif
622  }
623  }
624  }
626  math::XYZTLorentzVector mindRvec1;
627  double mindR1(999);
628  for (unsigned int i = 0; i < vecL2.size(); i++) {
629  double dr = dR(vecL1[0], vecL2[i]);
630 #ifdef EDM_ML_DEBUG
631  edm::LogVerbatim("HcalIsoTrack") << "lvl2[" << i << "] dR " << dr;
632 #endif
633  if (dr < mindR1) {
634  mindR1 = dr;
635  mindRvec1 = vecL2[i];
636  }
637  }
638 #ifdef EDM_ML_DEBUG
639  edm::LogVerbatim("HcalIsoTrack") << "L2 object closest to L1 " << mindRvec1 << " at Dr " << mindR1;
640 #endif
641 
642  if (!vecL1.empty()) {
643  t_l1pt = vecL1[0].pt();
644  t_l1eta = vecL1[0].eta();
645  t_l1phi = vecL1[0].phi();
646  } else {
647  t_l1pt = t_l1eta = t_l1phi = 0;
648  }
649  if (!vecL3.empty()) {
650  t_l3pt = vecL3[0].pt();
651  t_l3eta = vecL3[0].eta();
652  t_l3phi = vecL3[0].phi();
653  } else {
654  t_l3pt = t_l3eta = t_l3phi = 0;
655  }
656  // Now fill in the tree for each selected track
657  if (!done) {
658  ntksave = fillTree(vecL1,
659  vecL3,
660  leadPV,
661  trkCaloDirections,
662  trkCaloDets,
663  geo,
664  caloTopology,
665  theHBHETopology,
666  theEcalChStatus,
667  theEcalSevlv,
668  barrelRecHitsHandle,
669  endcapRecHitsHandle,
670  hbhe,
671  caloTower,
672  genParticles,
673  respCorrs);
674  t_TracksSaved += ntksave[0];
675  t_TracksLoose += ntksave[1];
676  t_TracksTight += ntksave[2];
677  done = true;
678  }
679  }
680  }
681  }
682  }
683  }
684 #ifdef EDM_ML_DEBUG
685  edm::LogVerbatim("HcalIsoTrack") << "Final results on selected tracks " << t_TracksSaved << ":" << t_TracksLoose
686  << ":" << t_TracksTight;
687 #endif
688  t_TrigPassSel = (t_TracksSaved > 0);
689  tree2->Fill();
690 }
691 
692 void HcalIsoTrkAnalyzer::beginJob() {
693  tree = fs->make<TTree>("CalibTree", "CalibTree");
694 
695  tree->Branch("t_Run", &t_Run, "t_Run/I");
696  tree->Branch("t_Event", &t_Event, "t_Event/I");
697  tree->Branch("t_DataType", &t_DataType, "t_DataType/I");
698  tree->Branch("t_ieta", &t_ieta, "t_ieta/I");
699  tree->Branch("t_iphi", &t_iphi, "t_iphi/I");
700  tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
701  tree->Branch("t_nVtx", &t_nVtx, "t_nVtx/I");
702  tree->Branch("t_nTrk", &t_nTrk, "t_nTrk/I");
703  tree->Branch("t_goodPV", &t_goodPV, "t_goodPV/I");
704  tree->Branch("t_l1pt", &t_l1pt, "t_l1pt/D");
705  tree->Branch("t_l1eta", &t_l1eta, "t_l1eta/D");
706  tree->Branch("t_l1phi", &t_l1phi, "t_l1phi/D");
707  tree->Branch("t_l3pt", &t_l3pt, "t_l3pt/D");
708  tree->Branch("t_l3eta", &t_l3eta, "t_l3eta/D");
709  tree->Branch("t_l3phi", &t_l3phi, "t_l3phi/D");
710  tree->Branch("t_p", &t_p, "t_p/D");
711  tree->Branch("t_pt", &t_pt, "t_pt/D");
712  tree->Branch("t_phi", &t_phi, "t_phi/D");
713  tree->Branch("t_mindR1", &t_mindR1, "t_mindR1/D");
714  tree->Branch("t_mindR2", &t_mindR2, "t_mindR2/D");
715  tree->Branch("t_eMipDR", &t_eMipDR, "t_eMipDR/D");
716  tree->Branch("t_eMipDR2", &t_eMipDR2, "t_eMipDR2/D");
717  tree->Branch("t_eMipDR3", &t_eMipDR3, "t_eMipDR3/D");
718  tree->Branch("t_eMipDR4", &t_eMipDR4, "t_eMipDR4/D");
719  tree->Branch("t_eMipDR5", &t_eMipDR5, "t_eMipDR5/D");
720  tree->Branch("t_eHcal", &t_eHcal, "t_eHcal/D");
721  tree->Branch("t_eHcal10", &t_eHcal10, "t_eHcal10/D");
722  tree->Branch("t_eHcal30", &t_eHcal30, "t_eHcal30/D");
723  tree->Branch("t_hmaxNearP", &t_hmaxNearP, "t_hmaxNearP/D");
724  tree->Branch("t_emaxNearP", &t_emaxNearP, "t_emaxNearP/D");
725  tree->Branch("t_eAnnular", &t_eAnnular, "t_eAnnular/D");
726  tree->Branch("t_hAnnular", &t_hAnnular, "t_hAnnular/D");
727  tree->Branch("t_rhoh", &t_rhoh, "t_rhoh/D");
728  tree->Branch("t_selectTk", &t_selectTk, "t_selectTk/O");
729  tree->Branch("t_qltyFlag", &t_qltyFlag, "t_qltyFlag/O");
730  tree->Branch("t_qltyMissFlag", &t_qltyMissFlag, "t_qltyMissFlag/O");
731  tree->Branch("t_qltyPVFlag", &t_qltyPVFlag, "t_qltyPVFlag/O");
732  tree->Branch("t_gentrackP", &t_gentrackP, "t_gentrackP/D");
733 
734  t_DetIds = new std::vector<unsigned int>();
735  t_DetIds1 = new std::vector<unsigned int>();
736  t_DetIds3 = new std::vector<unsigned int>();
737  t_HitEnergies = new std::vector<double>();
738  t_HitEnergies1 = new std::vector<double>();
739  t_HitEnergies3 = new std::vector<double>();
740  t_trgbits = new std::vector<bool>();
741  tree->Branch("t_DetIds", "std::vector<unsigned int>", &t_DetIds);
742  tree->Branch("t_HitEnergies", "std::vector<double>", &t_HitEnergies);
743  tree->Branch("t_trgbits", "std::vector<bool>", &t_trgbits);
744  tree->Branch("t_DetIds1", "std::vector<unsigned int>", &t_DetIds1);
745  tree->Branch("t_DetIds3", "std::vector<unsigned int>", &t_DetIds3);
746  tree->Branch("t_HitEnergies1", "std::vector<double>", &t_HitEnergies1);
747  tree->Branch("t_HitEnergies3", "std::vector<double>", &t_HitEnergies3);
748 
749  tree2 = fs->make<TTree>("EventInfo", "Event Information");
750 
751  tree2->Branch("t_RunNo", &t_RunNo, "t_RunNo/i");
752  tree2->Branch("t_EventNo", &t_EventNo, "t_EventNo/i");
753  tree2->Branch("t_Tracks", &t_Tracks, "t_Tracks/I");
754  tree2->Branch("t_TracksProp", &t_TracksProp, "t_TracksProp/I");
755  tree2->Branch("t_TracksSaved", &t_TracksSaved, "t_TracksSaved/I");
756  tree2->Branch("t_TracksLoose", &t_TracksLoose, "t_TracksLoose/I");
757  tree2->Branch("t_TracksTight", &t_TracksTight, "t_TracksTight/I");
758  tree2->Branch("t_TrigPass", &t_TrigPass, "t_TrigPass/O");
759  tree2->Branch("t_TrigPassSel", &t_TrigPassSel, "t_TrigPassSel/O");
760  tree2->Branch("t_L1Bit", &t_L1Bit, "t_L1Bit/O");
761  tree2->Branch("t_allvertex", &t_allvertex, "t_allvertex/I");
762  t_hltbits = new std::vector<bool>();
763  t_ietaAll = new std::vector<int>();
764  t_ietaGood = new std::vector<int>();
765  t_trackType = new std::vector<int>();
766  tree2->Branch("t_ietaAll", "std::vector<int>", &t_ietaAll);
767  tree2->Branch("t_ietaGood", "std::vector<int>", &t_ietaGood);
768  tree2->Branch("t_trackType", "std::vector<int>", &t_trackType);
769  tree2->Branch("t_hltbits", "std::vector<bool>", &t_hltbits);
770 }
771 
772 // ------------ method called when starting to processes a run ------------
773 void HcalIsoTrkAnalyzer::beginRun(edm::Run const& iRun, edm::EventSetup const& iSetup) {
774  hdc_ = &iSetup.getData(tok_ddrec_);
775 
776  bool changed_(true);
777  bool flag = hltConfig_.init(iRun, iSetup, processName_, changed_);
778  edm::LogVerbatim("HcalIsoTrack") << "Run[" << nRun_ << "] " << iRun.run() << " process " << processName_
779  << " init flag " << flag << " change flag " << changed_;
780  // check if trigger names in (new) config
781  if (changed_) {
782 #ifdef EDM_ML_DEBUG
783  edm::LogVerbatim("HcalIsoTrack") << "New trigger menu found !!!";
784 #endif
785  const unsigned int n(hltConfig_.size());
786  for (unsigned itrig = 0; itrig < trigNames_.size(); itrig++) {
787  unsigned int triggerindx = hltConfig_.triggerIndex(trigNames_[itrig]);
788  if (triggerindx >= n) {
789  edm::LogWarning("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " does not exist in "
790  << "the current menu";
791 #ifdef EDM_ML_DEBUG
792  } else {
793  edm::LogVerbatim("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " exists";
794 #endif
795  }
796  }
797  }
798 }
799 
800 // ------------ method called when ending the processing of a run ------------
801 void HcalIsoTrkAnalyzer::endRun(edm::Run const& iRun, edm::EventSetup const&) {
802  nRun_++;
803  edm::LogVerbatim("HcalIsoTrack") << "endRun[" << nRun_ << "] " << iRun.run();
804 }
805 
806 void HcalIsoTrkAnalyzer::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
807  edm::ParameterSetDescription desc;
808  std::vector<std::string> trig = {"HLT_PFJet40",
809  "HLT_PFJet60",
810  "HLT_PFJet80",
811  "HLT_PFJet140",
812  "HLT_PFJet200",
813  "HLT_PFJet260",
814  "HLT_PFJet320",
815  "HLT_PFJet400",
816  "HLT_PFJet450",
817  "HLT_PFJet500"};
818  desc.add<std::vector<std::string>>("triggers", trig);
819  desc.add<std::string>("processName", "HLT");
820  desc.add<std::string>("l1Filter", "");
821  desc.add<std::string>("l2Filter", "L2Filter");
822  desc.add<std::string>("l3Filter", "Filter");
823  // following 10 parameters are parameters to select good tracks
824  desc.add<std::string>("trackQuality", "highPurity");
825  desc.add<double>("minTrackPt", 1.0);
826  desc.add<double>("maxDxyPV", 0.02);
827  desc.add<double>("maxDzPV", 0.02);
828  desc.add<double>("maxChi2", 5.0);
829  desc.add<double>("maxDpOverP", 0.1);
830  desc.add<int>("minOuterHit", 4);
831  desc.add<int>("minLayerCrossed", 8);
832  desc.add<int>("maxInMiss", 0);
833  desc.add<int>("maxOutMiss", 0);
834  // Minimum momentum of selected isolated track and signal zone
835  desc.add<double>("minimumTrackP", 20.0);
836  desc.add<double>("coneRadius", 34.98);
837  // signal zone in ECAL and MIP energy cutoff
838  desc.add<double>("coneRadiusMIP", 14.0);
839  desc.add<double>("coneRadiusMIP2", 18.0);
840  desc.add<double>("coneRadiusMIP3", 20.0);
841  desc.add<double>("coneRadiusMIP4", 22.0);
842  desc.add<double>("coneRadiusMIP5", 24.0);
843  desc.add<double>("maximumEcalEnergy", 2.0);
844  // following 4 parameters are for isolation cuts and described in the code
845  desc.add<double>("maxTrackP", 8.0);
846  desc.add<double>("slopeTrackP", 0.05090504066);
847  desc.add<double>("isolationEnergyTight", 2.0);
848  desc.add<double>("isolationEnergyLoose", 10.0);
849  // energy thershold for ECAL (from Egamma group)
850  desc.add<double>("EBHitEnergyThreshold", 0.08);
851  desc.add<double>("EEHitEnergyThreshold0", 0.30);
852  desc.add<double>("EEHitEnergyThreshold1", 0.00);
853  desc.add<double>("EEHitEnergyThreshold2", 0.00);
854  desc.add<double>("EEHitEnergyThreshold3", 0.00);
855  desc.add<double>("EEHitEnergyThresholdLow", 0.30);
856  desc.add<double>("EEHitEnergyThresholdHigh", 0.30);
857  // prescale factors
858  desc.add<double>("momentumLow", 40.0);
859  desc.add<double>("momentumHigh", 60.0);
860  desc.add<int>("prescaleLow", 1);
861  desc.add<int>("prescaleHigh", 1);
862  // various labels for collections used in the code
863  desc.add<edm::InputTag>("labelTriggerEvent", edm::InputTag("hltTriggerSummaryAOD", "", "HLT"));
864  desc.add<edm::InputTag>("labelTriggerResult", edm::InputTag("TriggerResults", "", "HLT"));
865  desc.add<std::string>("labelTrack", "generalTracks");
866  desc.add<std::string>("labelVertex", "offlinePrimaryVertices");
867  desc.add<std::string>("labelEBRecHit", "EcalRecHitsEB");
868  desc.add<std::string>("labelEERecHit", "EcalRecHitsEE");
869  desc.add<std::string>("labelHBHERecHit", "hbhereco");
870  desc.add<std::string>("labelBeamSpot", "offlineBeamSpot");
871  desc.add<std::string>("labelCaloTower", "towerMaker");
872  desc.add<edm::InputTag>("algInputTag", edm::InputTag("gtStage2Digis"));
873  desc.add<edm::InputTag>("extInputTag", edm::InputTag("gtStage2Digis"));
874  desc.addUntracked<std::string>("moduleName", "");
875  desc.addUntracked<std::string>("producerName", "");
876  // Various flags used for selecting tracks, choice of energy Method2/0
877  // Data type 0/1 for single jet trigger or others
878  desc.addUntracked<int>("useRaw", 0);
879  desc.addUntracked<bool>("ignoreTriggers", false);
880  desc.addUntracked<bool>("useL1Trigger", false);
881  desc.addUntracked<double>("hcalScale", 1.0);
882  desc.addUntracked<int>("dataType", 0);
883  desc.addUntracked<bool>("unCorrect", false);
884  desc.addUntracked<bool>("collapseDepth", false);
885  desc.addUntracked<std::string>("l1TrigName", "L1_SingleJet60");
886  desc.addUntracked<int>("outMode", 11);
887  std::vector<int> dummy;
888  desc.addUntracked<std::vector<int>>("oldID", dummy);
889  desc.addUntracked<std::vector<int>>("newDepth", dummy);
890  desc.addUntracked<bool>("hep17", false);
891  descriptions.add("HcalIsoTrkAnalyzer", desc);
892 }
893 
894 std::array<int, 3> HcalIsoTrkAnalyzer::fillTree(std::vector<math::XYZTLorentzVector>& vecL1,
895  std::vector<math::XYZTLorentzVector>& vecL3,
896  math::XYZPoint& leadPV,
897  std::vector<spr::propagatedTrackDirection>& trkCaloDirections,
898  std::vector<spr::propagatedTrackID>& trkCaloDets,
899  const CaloGeometry* geo,
900  const CaloTopology* caloTopology,
901  const HcalTopology* theHBHETopology,
902  const EcalChannelStatus* theEcalChStatus,
903  const EcalSeverityLevelAlgo* theEcalSevlv,
904  edm::Handle<EcalRecHitCollection>& barrelRecHitsHandle,
905  edm::Handle<EcalRecHitCollection>& endcapRecHitsHandle,
906  edm::Handle<HBHERecHitCollection>& hbhe,
907  edm::Handle<CaloTowerCollection>& tower,
908  edm::Handle<reco::GenParticleCollection>& genParticles,
909  const HcalRespCorrs* respCorrs) {
910  int nSave(0), nLoose(0), nTight(0);
911  //Loop over tracks
912  std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
913  unsigned int nTracks(0), nselTracks(0);
914  t_nTrk = trkCaloDirections.size();
915  t_rhoh = (tower.isValid()) ? rhoh(tower) : 0;
916  for (trkDetItr = trkCaloDirections.begin(), nTracks = 0; trkDetItr != trkCaloDirections.end();
917  trkDetItr++, nTracks++) {
918  const reco::Track* pTrack = &(*(trkDetItr->trkItr));
919  math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
920 #ifdef EDM_ML_DEBUG
921  edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) :" << pTrack->pt() << "|"
922  << pTrack->eta() << "|" << pTrack->phi() << "|" << pTrack->p();
923 #endif
924  t_mindR2 = 999;
925  for (unsigned int k = 0; k < vecL3.size(); ++k) {
926  double dr = dR(vecL3[k], v4);
927  if (dr < t_mindR2) {
928  t_mindR2 = dr;
929  }
930  }
931  t_mindR1 = (!vecL1.empty()) ? dR(vecL1[0], v4) : 999;
932 #ifdef EDM_ML_DEBUG
933  edm::LogVerbatim("HcalIsoTrack") << "Closest L3 object at dr :" << t_mindR2 << " and from L1 " << t_mindR1;
934 #endif
935  t_ieta = t_iphi = 0;
936  if (trkDetItr->okHCAL) {
937  HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
938  t_ieta = detId.ieta();
939  t_iphi = detId.iphi();
940  if (t_p > 40.0 && t_p <= 60.0)
941  t_ietaAll->emplace_back(t_ieta);
942  }
943  //Selection of good track
944  t_selectTk = spr::goodTrack(pTrack, leadPV, selectionParameter_, false);
945  spr::trackSelectionParameters oneCutParameters = selectionParameter_;
946  oneCutParameters.maxDxyPV = 10;
947  oneCutParameters.maxDzPV = 100;
948  oneCutParameters.maxInMiss = 2;
949  oneCutParameters.maxOutMiss = 2;
950  bool qltyFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
951  oneCutParameters = selectionParameter_;
952  oneCutParameters.maxDxyPV = 10;
953  oneCutParameters.maxDzPV = 100;
954  t_qltyMissFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
955  oneCutParameters = selectionParameter_;
956  oneCutParameters.maxInMiss = 2;
957  oneCutParameters.maxOutMiss = 2;
958  t_qltyPVFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
959  double eIsolation = maxRestrictionP_ * exp(slopeRestrictionP_ * std::abs((double)t_ieta));
960  if (eIsolation < eIsolate1_)
961  eIsolation = eIsolate1_;
962  if (eIsolation < eIsolate2_)
963  eIsolation = eIsolate2_;
964 #ifdef EDM_ML_DEBUG
965  edm::LogVerbatim("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " << qltyFlag << "|" << trkDetItr->okECAL << "|"
966  << trkDetItr->okHCAL << " eIsolation " << eIsolation;
967 #endif
968  t_qltyFlag = (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL);
969  if (t_qltyFlag) {
970  nselTracks++;
971  int nNearTRKs(0);
973  std::vector<DetId> eIds;
974  std::vector<double> eHit;
975  t_eMipDR = spr::eCone_ecal(geo,
976  barrelRecHitsHandle,
977  endcapRecHitsHandle,
978  trkDetItr->pointHCAL,
979  trkDetItr->pointECAL,
980  a_mipR_,
981  trkDetItr->directionECAL,
982  eIds,
983  eHit);
984  double eEcal(0);
985  for (unsigned int k = 0; k < eIds.size(); ++k) {
986  if (eHit[k] > eThreshold(eIds[k], geo))
987  eEcal += eHit[k];
988  }
989 #ifdef EDM_ML_DEBUG
990  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR << ":" << eEcal;
991 #endif
992  t_eMipDR = eEcal;
995  std::vector<DetId> eIds2;
996  std::vector<double> eHit2;
997  t_eMipDR2 = spr::eCone_ecal(geo,
998  barrelRecHitsHandle,
999  endcapRecHitsHandle,
1000  trkDetItr->pointHCAL,
1001  trkDetItr->pointECAL,
1002  a_mipR2_,
1003  trkDetItr->directionECAL,
1004  eIds2,
1005  eHit2);
1006  double eEcal2(0);
1007  for (unsigned int k = 0; k < eIds2.size(); ++k) {
1008  if (eHit2[k] > eThreshold(eIds2[k], geo))
1009  eEcal2 += eHit2[k];
1010  }
1011 #ifdef EDM_ML_DEBUG
1012  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR2 << ":" << eEcal2;
1013 #endif
1014  t_eMipDR2 = eEcal2;
1017  std::vector<DetId> eIds3;
1018  std::vector<double> eHit3;
1019  t_eMipDR3 = spr::eCone_ecal(geo,
1020  barrelRecHitsHandle,
1021  endcapRecHitsHandle,
1022  trkDetItr->pointHCAL,
1023  trkDetItr->pointECAL,
1024  a_mipR3_,
1025  trkDetItr->directionECAL,
1026  eIds3,
1027  eHit3);
1028  double eEcal3(0);
1029  for (unsigned int k = 0; k < eIds3.size(); ++k) {
1030  if (eHit3[k] > eThreshold(eIds3[k], geo))
1031  eEcal3 += eHit3[k];
1032  }
1033 #ifdef EDM_ML_DEBUG
1034  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR3 << ":" << eEcal3;
1035 #endif
1036  t_eMipDR3 = eEcal3;
1039  std::vector<DetId> eIds4;
1040  std::vector<double> eHit4;
1041  t_eMipDR4 = spr::eCone_ecal(geo,
1042  barrelRecHitsHandle,
1043  endcapRecHitsHandle,
1044  trkDetItr->pointHCAL,
1045  trkDetItr->pointECAL,
1046  a_mipR4_,
1047  trkDetItr->directionECAL,
1048  eIds4,
1049  eHit4);
1050  double eEcal4(0);
1051  for (unsigned int k = 0; k < eIds4.size(); ++k) {
1052  if (eHit4[k] > eThreshold(eIds4[k], geo))
1053  eEcal4 += eHit4[k];
1054  }
1055 #ifdef EDM_ML_DEBUG
1056  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR4 << ":" << eEcal4;
1057 #endif
1058  t_eMipDR4 = eEcal4;
1061  std::vector<DetId> eIds5;
1062  std::vector<double> eHit5;
1063  t_eMipDR5 = spr::eCone_ecal(geo,
1064  barrelRecHitsHandle,
1065  endcapRecHitsHandle,
1066  trkDetItr->pointHCAL,
1067  trkDetItr->pointECAL,
1068  a_mipR5_,
1069  trkDetItr->directionECAL,
1070  eIds5,
1071  eHit5);
1072  double eEcal5(0);
1073  for (unsigned int k = 0; k < eIds5.size(); ++k) {
1074  if (eHit5[k] > eThreshold(eIds5[k], geo))
1075  eEcal5 += eHit5[k];
1076  }
1077 #ifdef EDM_ML_DEBUG
1078  edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR5 << ":" << eEcal5;
1079 #endif
1080  t_eMipDR5 = eEcal5;
1082 
1083  t_emaxNearP = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 15, 15);
1084  const DetId cellE(trkDetItr->detIdECAL);
1085  std::pair<double, bool> e11x11P = spr::eECALmatrix(cellE,
1086  barrelRecHitsHandle,
1087  endcapRecHitsHandle,
1088  *theEcalChStatus,
1089  geo,
1090  caloTopology,
1091  theEcalSevlv,
1092  5,
1093  5,
1094  -100.0,
1095  -100.0,
1096  -100.0,
1097  100.0);
1098  std::pair<double, bool> e15x15P = spr::eECALmatrix(cellE,
1099  barrelRecHitsHandle,
1100  endcapRecHitsHandle,
1101  *theEcalChStatus,
1102  geo,
1103  caloTopology,
1104  theEcalSevlv,
1105  7,
1106  7,
1107  -100.0,
1108  -100.0,
1109  -100.0,
1110  100.0);
1111  if (e11x11P.second && e15x15P.second) {
1112  t_eAnnular = (e15x15P.first - e11x11P.first);
1113  } else {
1114  t_eAnnular = -(e15x15P.first - e11x11P.first);
1115  }
1116  t_hmaxNearP = spr::chargeIsolationCone(nTracks, trkCaloDirections, a_charIsoR_, nNearTRKs, false);
1117  const DetId cellH(trkDetItr->detIdHCAL);
1118  double h5x5 = spr::eHCALmatrix(
1119  theHBHETopology, cellH, hbhe, 2, 2, false, true, -100.0, -100.0, -100.0, -100.0, -100.0, 100.0);
1120  double h7x7 = spr::eHCALmatrix(
1121  theHBHETopology, cellH, hbhe, 3, 3, false, true, -100.0, -100.0, -100.0, -100.0, -100.0, 100.0);
1122  t_hAnnular = h7x7 - h5x5;
1123 #ifdef EDM_ML_DEBUG
1124  edm::LogVerbatim("HcalIsoTrack") << "max p Near (Ecal) " << t_emaxNearP << " (Hcal) " << t_hmaxNearP
1125  << " Annular E (Ecal) " << e11x11P.first << ":" << e15x15P.first << ":"
1126  << t_eAnnular << " (Hcal) " << h5x5 << ":" << h7x7 << ":" << t_hAnnular;
1127 #endif
1128  t_gentrackP = trackP(pTrack, genParticles);
1129  if (t_eMipDR < eEcalMax_ && t_hmaxNearP < eIsolation) {
1130  t_DetIds->clear();
1131  t_HitEnergies->clear();
1132  t_DetIds1->clear();
1133  t_HitEnergies1->clear();
1134  t_DetIds3->clear();
1135  t_HitEnergies3->clear();
1136  int nRecHits(-999), nRecHits1(-999), nRecHits3(-999);
1137  std::vector<DetId> ids, ids1, ids3;
1138  std::vector<double> edet0, edet1, edet3;
1139  t_eHcal = spr::eCone_hcal(geo,
1140  hbhe,
1141  trkDetItr->pointHCAL,
1142  trkDetItr->pointECAL,
1143  a_coneR_,
1144  trkDetItr->directionHCAL,
1145  nRecHits,
1146  ids,
1147  edet0,
1148  useRaw_);
1149  if (!oldID_.empty()) {
1150  for (unsigned k = 0; k < ids.size(); ++k)
1151  ids[k] = newId(ids[k]);
1152  }
1153  storeEnergy(0, respCorrs, ids, edet0, t_eHcal, t_DetIds, t_HitEnergies);
1154 
1155  //----- hcal energy in the extended cone 1 (a_coneR+10) --------------
1156  t_eHcal10 = spr::eCone_hcal(geo,
1157  hbhe,
1158  trkDetItr->pointHCAL,
1159  trkDetItr->pointECAL,
1160  a_coneR1_,
1161  trkDetItr->directionHCAL,
1162  nRecHits1,
1163  ids1,
1164  edet1,
1165  useRaw_);
1166  if (!oldID_.empty()) {
1167  for (unsigned k = 0; k < ids1.size(); ++k)
1168  ids1[k] = newId(ids1[k]);
1169  }
1170  storeEnergy(1, respCorrs, ids1, edet1, t_eHcal10, t_DetIds1, t_HitEnergies1);
1171 
1172  //----- hcal energy in the extended cone 3 (a_coneR+30) --------------
1173  t_eHcal30 = spr::eCone_hcal(geo,
1174  hbhe,
1175  trkDetItr->pointHCAL,
1176  trkDetItr->pointECAL,
1177  a_coneR2_,
1178  trkDetItr->directionHCAL,
1179  nRecHits3,
1180  ids3,
1181  edet3,
1182  useRaw_);
1183  if (!oldID_.empty()) {
1184  for (unsigned k = 0; k < ids3.size(); ++k)
1185  ids3[k] = newId(ids3[k]);
1186  }
1187  storeEnergy(3, respCorrs, ids3, edet3, t_eHcal30, t_DetIds3, t_HitEnergies3);
1188 
1189  t_p = pTrack->p();
1190  t_pt = pTrack->pt();
1191  t_phi = pTrack->phi();
1192 
1193 #ifdef EDM_ML_DEBUG
1194  edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) :" << t_pt << "|"
1195  << pTrack->eta() << "|" << t_phi << "|" << t_p << " Generator Level p "
1196  << t_gentrackP;
1197  edm::LogVerbatim("HcalIsoTrack") << "e_MIP " << t_eMipDR << " Chg Isolation " << t_hmaxNearP << " eHcal"
1198  << t_eHcal << " ieta " << t_ieta << " Quality " << t_qltyMissFlag << ":"
1199  << t_qltyPVFlag << ":" << t_selectTk;
1200  for (unsigned int ll = 0; ll < t_DetIds->size(); ll++) {
1201  edm::LogVerbatim("HcalIsoTrack")
1202  << "det id is = " << t_DetIds->at(ll) << " hit enery is = " << t_HitEnergies->at(ll);
1203  }
1204  for (unsigned int ll = 0; ll < t_DetIds1->size(); ll++) {
1205  edm::LogVerbatim("HcalIsoTrack")
1206  << "det id is = " << t_DetIds1->at(ll) << " hit enery is = " << t_HitEnergies1->at(ll);
1207  }
1208  for (unsigned int ll = 0; ll < t_DetIds3->size(); ll++) {
1209  edm::LogVerbatim("HcalIsoTrack")
1210  << "det id is = " << t_DetIds3->at(ll) << " hit enery is = " << t_HitEnergies3->at(ll);
1211  }
1212 #endif
1213  bool accept(false);
1214  if (t_p > pTrackMin_) {
1215  if (t_p < pTrackLow_) {
1216  ++nLow_;
1217  if (prescaleLow_ <= 1)
1218  accept = true;
1219  else if (nLow_ % prescaleLow_ == 1)
1220  accept = true;
1221  } else if (t_p > pTrackHigh_) {
1222  ++nHigh_;
1223  if (prescaleHigh_ <= 1)
1224  accept = true;
1225  else if (nHigh_ % prescaleHigh_ == 1)
1226  accept = true;
1227  } else {
1228  accept = true;
1229  }
1230  }
1231  if (accept) {
1232  tree->Fill();
1233  nSave++;
1234  int type(0);
1235  if (t_eMipDR < 1.0) {
1236  if (t_hmaxNearP < eIsolate2_) {
1237  ++nLoose;
1238  type = 1;
1239  }
1240  if (t_hmaxNearP < eIsolate1_) {
1241  ++nTight;
1242  type = 2;
1243  }
1244  }
1245  if (t_p > 40.0 && t_p <= 60.0 && t_selectTk) {
1246  t_ietaGood->emplace_back(t_ieta);
1247  t_trackType->emplace_back(type);
1248  }
1249 #ifdef EDM_ML_DEBUG
1250  for (unsigned int k = 0; k < t_trgbits->size(); k++) {
1251  edm::LogVerbatim("HcalIsoTrack") << "trigger bit is = " << t_trgbits->at(k);
1252  }
1253 #endif
1254  }
1255  }
1256  }
1257  }
1258  std::array<int, 3> i3{{nSave, nLoose, nTight}};
1259  return i3;
1260 }
1261 
1262 double HcalIsoTrkAnalyzer::dR(math::XYZTLorentzVector& vec1, math::XYZTLorentzVector& vec2) {
1263  return reco::deltaR(vec1.eta(), vec1.phi(), vec2.eta(), vec2.phi());
1264 }
1265 
1266 double HcalIsoTrkAnalyzer::trackP(const reco::Track* pTrack,
1267  const edm::Handle<reco::GenParticleCollection>& genParticles) {
1268  double pmom = -1.0;
1269  if (genParticles.isValid()) {
1270  double mindR(999.9);
1271  for (const auto& p : (*genParticles)) {
1272  double dR = reco::deltaR(pTrack->eta(), pTrack->phi(), p.momentum().Eta(), p.momentum().Phi());
1273  if (dR < mindR) {
1274  mindR = dR;
1275  pmom = p.momentum().R();
1276  }
1277  }
1278  }
1279  return pmom;
1280 }
1281 
1282 double HcalIsoTrkAnalyzer::rhoh(const edm::Handle<CaloTowerCollection>& tower) {
1283  std::vector<double> sumPFNallSMDQH2;
1284  sumPFNallSMDQH2.reserve(phibins_.size() * etabins_.size());
1285 
1286  for (auto eta : etabins_) {
1287  for (auto phi : phibins_) {
1288  double hadder = 0;
1289  for (const auto& pf_it : (*tower)) {
1290  if (fabs(eta - pf_it.eta()) > etahalfdist_)
1291  continue;
1292  if (fabs(reco::deltaPhi(phi, pf_it.phi())) > phihalfdist_)
1293  continue;
1294  hadder += pf_it.hadEt();
1295  }
1296  sumPFNallSMDQH2.emplace_back(hadder);
1297  }
1298  }
1299 
1300  double evt_smdq(0);
1301  std::sort(sumPFNallSMDQH2.begin(), sumPFNallSMDQH2.end());
1302  if (sumPFNallSMDQH2.size() % 2)
1303  evt_smdq = sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 1) / 2];
1304  else
1305  evt_smdq = (sumPFNallSMDQH2[sumPFNallSMDQH2.size() / 2] + sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 2) / 2]) / 2.;
1306  double rhoh = evt_smdq / (etadist_ * phidist_);
1307 #ifdef EDM_ML_DEBUG
1308  edm::LogVerbatim("HcalIsoTrack") << "Rho " << evt_smdq << ":" << rhoh;
1309 #endif
1310  return rhoh;
1311 }
1312 
1313 double HcalIsoTrkAnalyzer::eThreshold(const DetId& id, const CaloGeometry* geo) const {
1314  const GlobalPoint& pos = geo->getPosition(id);
1315  double eta = std::abs(pos.eta());
1316  double eThr(hitEthrEB_);
1317  if (id.subdetId() != EcalBarrel) {
1318  eThr = (((eta * hitEthrEE3_ + hitEthrEE2_) * eta + hitEthrEE1_) * eta + hitEthrEE0_);
1319  if (eThr < hitEthrEELo_)
1320  eThr = hitEthrEELo_;
1321  else if (eThr > hitEthrEEHi_)
1322  eThr = hitEthrEEHi_;
1323  }
1324  return eThr;
1325 }
1326 
1327 DetId HcalIsoTrkAnalyzer::newId(const DetId& id) {
1328  HcalDetId hid(id);
1329  if (hep17_ && ((hid.iphi() < 63) || (hid.iphi() > 66) || (hid.zside() < 0)))
1330  return id;
1331  for (unsigned int k = 0; k < oldID_.size(); ++k) {
1332  if ((hid.subdetId() == oldDet_[k]) && (hid.ietaAbs() == oldEta_[k]) && (hid.depth() == oldDepth_[k])) {
1333  return static_cast<DetId>(HcalDetId(hid.subdet(), hid.ieta(), hid.iphi(), newDepth_[k]));
1334  }
1335  }
1336  return id;
1337 }
1338 
1339 void HcalIsoTrkAnalyzer::storeEnergy(int indx,
1340  const HcalRespCorrs* respCorrs,
1341  const std::vector<DetId>& ids,
1342  std::vector<double>& edet,
1343  double& eHcal,
1344  std::vector<unsigned int>* detIds,
1345  std::vector<double>* hitEnergies) {
1346  double ehcal(0);
1347  if (unCorrect_) {
1348  for (unsigned int k = 0; k < ids.size(); ++k) {
1349  double corr = (respCorrs->getValues(ids[k]))->getValue();
1350  if (corr != 0)
1351  edet[k] /= corr;
1352  ehcal += edet[k];
1353  }
1354  } else {
1355  for (const auto& en : edet)
1356  ehcal += en;
1357  }
1358  if ((std::abs(ehcal - eHcal) > 0.001) && (!unCorrect_))
1359  edm::LogWarning("HcalIsoTrack") << "Check inconsistent energies: " << indx << " " << eHcal << ":" << ehcal
1360  << " from " << ids.size() << " cells";
1361  eHcal = hcalScale_ * ehcal;
1362 
1363  if (collapseDepth_) {
1364  std::map<HcalDetId, double> hitMap;
1365  for (unsigned int k = 0; k < ids.size(); ++k) {
1366  HcalDetId id = hdc_->mergedDepthDetId(HcalDetId(ids[k]));
1367  auto itr = hitMap.find(id);
1368  if (itr == hitMap.end()) {
1369  hitMap[id] = edet[k];
1370  } else {
1371  (itr->second) += edet[k];
1372  }
1373  }
1374  detIds->reserve(hitMap.size());
1375  hitEnergies->reserve(hitMap.size());
1376  for (const auto& hit : hitMap) {
1377  detIds->emplace_back(hit.first.rawId());
1378  hitEnergies->emplace_back(hit.second);
1379  }
1380  } else {
1381  detIds->reserve(ids.size());
1382  hitEnergies->reserve(ids.size());
1383  for (unsigned int k = 0; k < ids.size(); ++k) {
1384  detIds->emplace_back(ids[k].rawId());
1385  hitEnergies->emplace_back(edet[k]);
1386  }
1387  }
1388 #ifdef EDM_ML_DEBUG
1389  edm::LogVerbatim("HcalIsoTrack") << "Input to storeEnergy with " << ids.size() << " cells";
1390  for (unsigned int k = 0; k < ids.size(); ++k)
1391  edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId(ids[k]) << " E " << edet[k];
1392  edm::LogVerbatim("HcalIsoTrack") << "Output of storeEnergy with " << detIds->size() << " cells and Etot " << eHcal;
1393  for (unsigned int k = 0; k < detIds->size(); ++k)
1394  edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId((*detIds)[k]) << " E " << (*hitEnergies)[k];
1395 #endif
1396 }
1397 
1398 //define this as a plug-in
1399 DEFINE_FWK_MODULE(HcalIsoTrkAnalyzer);
trigger::TriggerObject::phi
float phi() const
Definition: TriggerObject.h:54
HcalIsoTrkAnalyzer::phidist_
double phidist_
Definition: HcalIsoTrkAnalyzer.cc:159
spr::trackSelectionParameters::maxDxyPV
double maxDxyPV
Definition: TrackSelection.h:22
HcalIsoTrkAnalyzer::beginRun
void beginRun(edm::Run const &, edm::EventSetup const &) override
Definition: HcalIsoTrkAnalyzer.cc:773
HcalDDDRecConstants::mergedDepthDetId
HcalDetId mergedDepthDetId(const HcalDetId &id) const
Definition: HcalDDDRecConstants.cc:647
HcalIsoTrkAnalyzer::tok_cala_
edm::EDGetTokenT< CaloTowerCollection > tok_cala_
Definition: HcalIsoTrkAnalyzer.cc:169
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:174
HcalIsoTrkAnalyzer::t_p
double t_p
Definition: HcalIsoTrkAnalyzer.cc:186
HcalRespCorrsRcd.h
testProducerWithPsetDescEmpty_cfi.i3
i3
Definition: testProducerWithPsetDescEmpty_cfi.py:47
electrons_cff.bool
bool
Definition: electrons_cff.py:393
EDAnalyzer.h
HcalRespCorrs
Definition: HcalRespCorrs.h:17
HcalIsoTrkAnalyzer::tok_EB_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: HcalIsoTrkAnalyzer.cc:166
mps_fire.i
i
Definition: mps_fire.py:428
HcalIsoTrkAnalyzer::t_TracksTight
int t_TracksTight
Definition: HcalIsoTrkAnalyzer.cc:201
HcalIsoTrkAnalyzer::prescaleLow_
const int prescaleLow_
Definition: HcalIsoTrkAnalyzer.cc:142
HcalIsoTrkAnalyzer::l1TrigName_
const std::string l1TrigName_
Definition: HcalIsoTrkAnalyzer.cc:151
HcalIsoTrkAnalyzer::hitEthrEE3_
const double hitEthrEE3_
Definition: HcalIsoTrkAnalyzer.cc:147
HcalIsoTrkAnalyzer::labelRecVtx_
const std::string labelRecVtx_
Definition: HcalIsoTrkAnalyzer.cc:150
LogMessageMonitor_cff.modules
modules
Definition: LogMessageMonitor_cff.py:7
HcalIsoTrkAnalyzer::t_trgbits
std::vector< bool > * t_trgbits
Definition: HcalIsoTrkAnalyzer.cc:198
HcalIsoTrkAnalyzer::t_nVtx
int t_nVtx
Definition: HcalIsoTrkAnalyzer.cc:185
genParticles2HepMC_cfi.genParticles
genParticles
Definition: genParticles2HepMC_cfi.py:4
HcalIsoTrkAnalyzer::t_DetIds3
std::vector< unsigned int > * t_DetIds3
Definition: HcalIsoTrkAnalyzer.cc:196
HcalIsoTrkAnalyzer::fs
edm::Service< TFileService > fs
Definition: HcalIsoTrkAnalyzer.cc:130
HcalIsoTrkAnalyzer::t_eHcal10
double t_eHcal10
Definition: HcalIsoTrkAnalyzer.cc:193
funct::false
false
Definition: Factorize.h:29
HcalIsoTrkAnalyzer::a_mipR4_
const double a_mipR4_
Definition: HcalIsoTrkAnalyzer.cc:138
BXVector.h
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:175
TriggerResults.h
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:185
HcalDetId::iphi
constexpr int iphi() const
get the cell iphi
Definition: HcalDetId.h:157
spr::trackSelectionParameters::maxChi2
double maxChi2
Definition: TrackSelection.h:22
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:631
HcalIsoTrkAnalyzer::selectionParameter_
spr::trackSelectionParameters selectionParameter_
Definition: HcalIsoTrkAnalyzer.cc:133
HcalIsoTrkAnalyzer::oldID_
const std::vector< int > oldID_
Definition: HcalIsoTrkAnalyzer.cc:152
HcalIsoTrkAnalyzer::tok_bs_
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
Definition: HcalIsoTrkAnalyzer.cc:165
EcalSeverityLevelAlgoRcd.h
edm::Run
Definition: Run.h:45
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
TriggerEvent.h
edm::EDGetTokenT< trigger::TriggerEvent >
TrackBase.h
edm
HLT enums.
Definition: AlignableModifier.h:19
HcalIsoTrkAnalyzer::a_coneR1_
double a_coneR1_
Definition: HcalIsoTrkAnalyzer.cc:155
GlobalExtBlk.h
HcalIsoTrkAnalyzer::nLow_
unsigned int nLow_
Definition: HcalIsoTrkAnalyzer.cc:154
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
tree
Definition: tree.py:1
HcalIsoTrkAnalyzer::t_mindR1
double t_mindR1
Definition: HcalIsoTrkAnalyzer.cc:189
pos
Definition: PixelAliasList.h:18
HcalIsoTrkAnalyzer::nHigh_
unsigned int nHigh_
Definition: HcalIsoTrkAnalyzer.cc:154
HcalTopology
Definition: HcalTopology.h:26
HcalIsoTrkAnalyzer::tok_htopo_
edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tok_htopo_
Definition: HcalIsoTrkAnalyzer.cc:179
HcalIsoTrkAnalyzer::t_DetIds1
std::vector< unsigned int > * t_DetIds1
Definition: HcalIsoTrkAnalyzer.cc:196
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89353
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
hgcalTowerProducer_cfi.tower
tower
Definition: hgcalTowerProducer_cfi.py:4
HcalIsoTrkAnalyzer::theTriggerResultsLabel_
const edm::InputTag theTriggerResultsLabel_
Definition: HcalIsoTrkAnalyzer.cc:149
triggerResults
static const std::string triggerResults
Definition: EdmProvDump.cc:45
HcalIsoTrkAnalyzer::eIsolate2_
const double eIsolate2_
Definition: HcalIsoTrkAnalyzer.cc:140
HcalIsoTrkAnalyzer::a_coneR_
const double a_coneR_
Definition: HcalIsoTrkAnalyzer.cc:137
HcalIsoTrkAnalyzer::tok_trigEvt_
edm::EDGetTokenT< trigger::TriggerEvent > tok_trigEvt_
Definition: HcalIsoTrkAnalyzer.cc:160
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:46
HcalIsoTrkAnalyzer::t_qltyFlag
bool t_qltyFlag
Definition: HcalIsoTrkAnalyzer.cc:194
HcalIsoTrkAnalyzer::tok_genTrack_
edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
Definition: HcalIsoTrkAnalyzer.cc:163
HcalIsoTrkAnalyzer::a_mipR_
const double a_mipR_
Definition: HcalIsoTrkAnalyzer.cc:137
HcalDetId::depth
constexpr int depth() const
get the tower depth
Definition: HcalDetId.h:164
reco::TrackBase::px
double px() const
x coordinate of momentum vector
Definition: TrackBase.h:640
HcalIsoTrkAnalyzer::t_TracksLoose
int t_TracksLoose
Definition: HcalIsoTrkAnalyzer.cc:201
HcalIsoTrkAnalyzer::tok_ew_
edm::EDGetTokenT< GenEventInfoProduct > tok_ew_
Definition: HcalIsoTrkAnalyzer.cc:170
HcalIsoTrkAnalyzer::t_hltbits
std::vector< bool > * t_hltbits
Definition: HcalIsoTrkAnalyzer.cc:198
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:184
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
HcalIsoTrkAnalyzer::t_qltyPVFlag
bool t_qltyPVFlag
Definition: HcalIsoTrkAnalyzer.cc:195
DDAxes::x
HcalIsoTrkAnalyzer::t_l3phi
double t_l3phi
Definition: HcalIsoTrkAnalyzer.cc:188
HcalIsoTrkAnalyzer::etabins_
std::vector< double > etabins_
Definition: HcalIsoTrkAnalyzer.cc:157
edm::one::EDAnalyzer
Definition: EDAnalyzer.h:30
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:173
EcalCondObjectContainer
Definition: EcalCondObjectContainer.h:13
vec2
std::vector< vec1 > vec2
Definition: HCALResponse.h:16
BeamSpotPI::nTracks
Definition: BeamSpotPayloadInspectorHelper.h:41
HcalIsoTrkAnalyzer::useRaw_
const int useRaw_
Definition: HcalIsoTrkAnalyzer.cc:143
convertSQLiteXML.ok
bool ok
Definition: convertSQLiteXML.py:98
HcalIsoTrkAnalyzer::t_HitEnergies
std::vector< double > * t_HitEnergies
Definition: HcalIsoTrkAnalyzer.cc:197
edm::Handle
Definition: AssociativeIterator.h:50
HcalIsoTrkAnalyzer::trigNames_
const std::vector< std::string > trigNames_
Definition: HcalIsoTrkAnalyzer.cc:132
edm::RunBase::run
RunNumber_t run() const
Definition: RunBase.h:40
HcalIsoTrkAnalyzer::l3Filter_
const std::string l3Filter_
Definition: HcalIsoTrkAnalyzer.cc:136
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:191
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
EcalBarrel
Definition: EcalSubdetector.h:10
EcalRecHitCollections.h
HcalIsoTrkAnalyzer::tok_caloTopology_
edm::ESGetToken< CaloTopology, CaloTopologyRecord > tok_caloTopology_
Definition: HcalIsoTrkAnalyzer.cc:178
HcalIsoTrkAnalyzer::mode_
const int mode_
Definition: HcalIsoTrkAnalyzer.cc:143
HcalIsoTrkAnalyzer::pTrackHigh_
const double pTrackHigh_
Definition: HcalIsoTrkAnalyzer.cc:141
GenParticle.h
HcalIsoTrkAnalyzer::rhoh
double rhoh(const edm::Handle< CaloTowerCollection > &)
Definition: HcalIsoTrkAnalyzer.cc:1282
HcalIsoTrkAnalyzer::oldDet_
std::vector< int > oldDet_
Definition: HcalIsoTrkAnalyzer.cc:158
HcalIsoTrkAnalyzer::t_Event
int t_Event
Definition: HcalIsoTrkAnalyzer.cc:184
HcalIsoTrkAnalyzer::t_l1eta
double t_l1eta
Definition: HcalIsoTrkAnalyzer.cc:187
HcalIsoTrkAnalyzer::t_ietaGood
std::vector< int > * t_ietaGood
Definition: HcalIsoTrkAnalyzer.cc:202
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:148
DetId
Definition: DetId.h:17
HcalIsoTrkAnalyzer::t_hAnnular
double t_hAnnular
Definition: HcalIsoTrkAnalyzer.cc:192
HcalIsoTrkAnalyzer::tok_EE_
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: HcalIsoTrkAnalyzer.cc:167
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:30
MakerMacros.h
trigger::TriggerObject::py
float py() const
Definition: TriggerObject.h:58
trigger::TriggerObject::pt
float pt() const
Definition: TriggerObject.h:52
CaloGeometry
Definition: CaloGeometry.h:21
Track.h
HcalIsoTrkAnalyzer::t_eMipDR5
double t_eMipDR5
Definition: HcalIsoTrkAnalyzer.cc:191
HcalIsoTrkAnalyzer::t_phi
double t_phi
Definition: HcalIsoTrkAnalyzer.cc:186
TrackFwd.h
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
HcalIsoTrkAnalyzer::dR
double dR(math::XYZTLorentzVector &, math::XYZTLorentzVector &)
Definition: HcalIsoTrkAnalyzer.cc:1262
BeamSpot.h
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
Service.h
PVValHelper::eta
Definition: PVValidationHelpers.h:69
l1t::L1TGlobalUtil
Definition: L1TGlobalUtil.h:44
ChargeIsolation.h
HcalIsoTrkAnalyzer::t_l1phi
double t_l1phi
Definition: HcalIsoTrkAnalyzer.cc:187
HcalIsoTrkAnalyzer::phibins_
std::vector< double > phibins_
Definition: HcalIsoTrkAnalyzer.cc:157
HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer
~HcalIsoTrkAnalyzer() override
Definition: HcalIsoTrkAnalyzer.cc:90
HcalIsoTrkAnalyzer::newId
DetId newId(const DetId &)
Definition: HcalIsoTrkAnalyzer.cc:1327
eHCALMatrix.h
HcalIsoTrkAnalyzer::labelEE_
const std::string labelEE_
Definition: HcalIsoTrkAnalyzer.cc:151
GenParticleFwd.h
HcalIsoTrkAnalyzer::t_TrigPass
bool t_TrigPass
Definition: HcalIsoTrkAnalyzer.cc:195
DDAxes::z
reco::Track
Definition: Track.h:27
IdealMagneticFieldRecord.h
HcalIsoTrkAnalyzer::t_L1Bit
bool t_L1Bit
Definition: HcalIsoTrkAnalyzer.cc:199
spr::trackSelectionParameters::minPt
double minPt
Definition: TrackSelection.h:20
HcalIsoTrkAnalyzer::phihalfdist_
double phihalfdist_
Definition: HcalIsoTrkAnalyzer.cc:159
spr::trackSelectionParameters
Definition: TrackSelection.h:12
spr::trackSelectionParameters::minOuterHit
int minOuterHit
Definition: TrackSelection.h:23
HcalIsoTrkAnalyzer::t_selectTk
bool t_selectTk
Definition: HcalIsoTrkAnalyzer.cc:194
HcalIsoTrkAnalyzer::tok_trigRes_
edm::EDGetTokenT< edm::TriggerResults > tok_trigRes_
Definition: HcalIsoTrkAnalyzer.cc:161
trigger::TriggerObject::px
float px() const
Definition: TriggerObject.h:57
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:190
HcalIsoTrkAnalyzer::t_goodPV
int t_goodPV
Definition: HcalIsoTrkAnalyzer.cc:185
HcalIsoTrkAnalyzer::hitEthrEELo_
const double hitEthrEELo_
Definition: HcalIsoTrkAnalyzer.cc:148
eECALMatrix.h
spr::chargeIsolationCone
double chargeIsolationCone(unsigned int trkIndex, std::vector< spr::propagatedTrackDirection > &trkDirs, double dR, int &nNearTRKs, bool debug=false)
Definition: ChargeIsolation.cc:288
eCone.h
EcalSeverityLevelAlgo.h
HcalIsoTrkAnalyzer::useL1Trigger_
const bool useL1Trigger_
Definition: HcalIsoTrkAnalyzer.cc:144
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:184
DetId::subdetId
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum)
Definition: DetId.h:48
spr::trackSelectionParameters::maxDpOverP
double maxDpOverP
Definition: TrackSelection.h:22
edm::ConfigurationDescriptions
Definition: ConfigurationDescriptions.h:28
spr::trackSelectionParameters::minLayerCrossed
int minLayerCrossed
Definition: TrackSelection.h:23
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
Vertex.h
HcalIsoTrkAnalyzer::t_eMipDR3
double t_eMipDR3
Definition: HcalIsoTrkAnalyzer.cc:190
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:183
HcalIsoTrkAnalyzer::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: HcalIsoTrkAnalyzer.cc:806
HcalIsoTrkAnalyzer::tok_alg_
edm::EDGetTokenT< BXVector< GlobalAlgBlk > > tok_alg_
Definition: HcalIsoTrkAnalyzer.cc:171
TFileService.h
trigger::TriggerObject::energy
float energy() const
Definition: TriggerObject.h:61
HcalIsoTrkAnalyzer::a_mipR2_
const double a_mipR2_
Definition: HcalIsoTrkAnalyzer.cc:137
CaloSubdetectorGeometry.h
HcalDetId::ieta
constexpr int ieta() const
get the cell ieta
Definition: HcalDetId.h:155
HcalIsoTrkAnalyzer::a_charIsoR_
double a_charIsoR_
Definition: HcalIsoTrkAnalyzer.cc:155
HcalIsoTrkAnalyzer::t_eHcal
double t_eHcal
Definition: HcalIsoTrkAnalyzer.cc:193
edm::ParameterSet
Definition: ParameterSet.h:47
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
HcalIsoTrkAnalyzer::t_qltyMissFlag
bool t_qltyMissFlag
Definition: HcalIsoTrkAnalyzer.cc:194
HcalIsoTrkAnalyzer::newDepth_
const std::vector< int > newDepth_
Definition: HcalIsoTrkAnalyzer.cc:152
GenEventInfoProduct.h
Event.h
HcalIsoTrkAnalyzer::tok_resp_
edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_resp_
Definition: HcalIsoTrkAnalyzer.cc:180
LumiDetails.h
HcalIsoTrkAnalyzer::eEcalMax_
const double eEcalMax_
Definition: HcalIsoTrkAnalyzer.cc:138
deltaR.h
HcalIsoTrkAnalyzer::fillTree
std::array< int, 3 > fillTree(std::vector< math::XYZTLorentzVector > &vecL1, std::vector< math::XYZTLorentzVector > &vecL3, math::XYZPoint &leadPV, std::vector< spr::propagatedTrackDirection > &trkCaloDirections, std::vector< spr::propagatedTrackID > &trkCaloDets, const CaloGeometry *geo, const CaloTopology *topo, const HcalTopology *theHBHETopology, const EcalChannelStatus *theEcalChStatus, const EcalSeverityLevelAlgo *theEcalSevlv, edm::Handle< EcalRecHitCollection > &barrelRecHitsHandle, edm::Handle< EcalRecHitCollection > &endcapRecHitsHandle, edm::Handle< HBHERecHitCollection > &hbhe, edm::Handle< CaloTowerCollection > &towerHandle, edm::Handle< reco::GenParticleCollection > &genParticles, const HcalRespCorrs *respCorrs)
Definition: HcalIsoTrkAnalyzer.cc:894
type
type
Definition: SiPixelVCal_PayloadInspector.cc:37
TrackSelection.h
HcalIsoTrkAnalyzer::t_trackType
std::vector< int > * t_trackType
Definition: HcalIsoTrkAnalyzer.cc:202
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:201
gainCalibHelper::gainCalibPI::type
type
Definition: SiPixelGainCalibHelper.h:39
spr::trackSelectionParameters::maxOutMiss
int maxOutMiss
Definition: TrackSelection.h:24
HcalIsoTrkAnalyzer::HcalIsoTrkAnalyzer
HcalIsoTrkAnalyzer(edm::ParameterSet const &)
Definition: HcalIsoTrkAnalyzer.cc:205
HcalIsoTrkAnalyzer::t_emaxNearP
double t_emaxNearP
Definition: HcalIsoTrkAnalyzer.cc:192
HcalIsoTrkAnalyzer::hdc_
const HcalDDDRecConstants * hdc_
Definition: HcalIsoTrkAnalyzer.cc:156
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:652
trigger::TriggerObject::id
int id() const
getters
Definition: TriggerObject.h:51
CaloTopologyRecord.h
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:158
HcalDetId::subdet
constexpr HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:138
LorentzVector.h
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:190
HcalDetId
Definition: HcalDetId.h:12
edm::Service< TFileService >
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)
createfilelist.int
int
Definition: createfilelist.py:10
TriggerNames.h
iEvent
int iEvent
Definition: GenABIO.cc:224
HcalIsoTrkAnalyzer::labelGenTrack_
const std::string labelGenTrack_
Definition: HcalIsoTrkAnalyzer.cc:150
HcalIsoTrkAnalyzer::tok_geom_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord > tok_geom_
Definition: HcalIsoTrkAnalyzer.cc:177
HcalIsoTrkAnalyzer::a_mipR5_
const double a_mipR5_
Definition: HcalIsoTrkAnalyzer.cc:138
HcalIsoTrkAnalyzer::hcalScale_
const double hcalScale_
Definition: HcalIsoTrkAnalyzer.cc:140
HcalIsoTrkAnalyzer::tok_sevlv_
edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd > tok_sevlv_
Definition: HcalIsoTrkAnalyzer.cc:176
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:1313
trackerHitRTTI::vector
Definition: trackerHitRTTI.h:21
GlobalAlgBlk.h
trigger::TriggerObject::eta
float eta() const
Definition: TriggerObject.h:53
HcalIsoTrkAnalyzer::t_l3eta
double t_l3eta
Definition: HcalIsoTrkAnalyzer.cc:188
HcalIsoTrkAnalyzer::slopeRestrictionP_
const double slopeRestrictionP_
Definition: HcalIsoTrkAnalyzer.cc:139
HcalIsoTrkAnalyzer::maxRestrictionP_
const double maxRestrictionP_
Definition: HcalIsoTrkAnalyzer.cc:139
MagneticField.h
edm::EventSetup
Definition: EventSetup.h:57
HcalRespCorrs.h
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:196
HLTConfigProvider.h
itr
std::vector< std::pair< float, float > >::iterator itr
Definition: HGCDigitizer.cc:29
HitPattern.h
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:184
ClusterTask_cfi.trig
trig
Definition: ClusterTask_cfi.py:398
HcalIsoTrkAnalyzer::a_coneR2_
double a_coneR2_
Definition: HcalIsoTrkAnalyzer.cc:155
HcalIsoTrkAnalyzer::t_iphi
int t_iphi
Definition: HcalIsoTrkAnalyzer.cc:184
edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord >
TriggerObject.h
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
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:120
CaloTopology.h
CaloTowerCollection.h
HcalIsoTrkAnalyzer::tree
TTree * tree
Definition: HcalIsoTrkAnalyzer.cc:182
photonIsolationHIProducer_cfi.hbhe
hbhe
Definition: photonIsolationHIProducer_cfi.py:8
HcalIsoTrkAnalyzer::nRun_
unsigned int nRun_
Definition: HcalIsoTrkAnalyzer.cc:154
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:144
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
HcalIsoTrkAnalyzer::hitEthrEE0_
const double hitEthrEE0_
Definition: HcalIsoTrkAnalyzer.cc:146
VertexFwd.h
DDAxes::phi
HcalIsoTrkAnalyzer::hitEthrEE2_
const double hitEthrEE2_
Definition: HcalIsoTrkAnalyzer.cc:147
CaloSubdetectorTopology.h
HcalIsoTrkAnalyzer::t_l1pt
double t_l1pt
Definition: HcalIsoTrkAnalyzer.cc:187
HcalIsoTrkAnalyzer::theTrackQuality_
const std::string theTrackQuality_
Definition: HcalIsoTrkAnalyzer.cc:134
HcalTopology.h
HcalIsoTrkAnalyzer::t_HitEnergies1
std::vector< double > * t_HitEnergies1
Definition: HcalIsoTrkAnalyzer.cc:197
EcalTrigTowerConstituentsMap.h
HcalIsoTrkAnalyzer::t_TracksProp
int t_TracksProp
Definition: HcalIsoTrkAnalyzer.cc:200
submitPVResolutionJobs.desc
string desc
Definition: submitPVResolutionJobs.py:251
HLTConfigProvider
Definition: HLTConfigProvider.h:29
std
Definition: JetResolutionObject.h:76
HcalIsoTrkAnalyzer::pTrackLow_
const double pTrackLow_
Definition: HcalIsoTrkAnalyzer.cc:141
trigger::TriggerObject::pz
float pz() const
Definition: TriggerObject.h:59
HcalIsoTrkAnalyzer::tok_parts_
edm::EDGetTokenT< reco::GenParticleCollection > tok_parts_
Definition: HcalIsoTrkAnalyzer.cc:162
HcalIsoTrkAnalyzer::t_ietaAll
std::vector< int > * t_ietaAll
Definition: HcalIsoTrkAnalyzer.cc:202
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:159
Calorimetry_cff.bField
bField
Definition: Calorimetry_cff.py:292
HcalIsoTrkAnalyzer::t_mindR2
double t_mindR2
Definition: HcalIsoTrkAnalyzer.cc:189
HcalIsoTrkAnalyzer::hitEthrEB_
const double hitEthrEB_
Definition: HcalIsoTrkAnalyzer.cc:146
HcalIsoTrkAnalyzer::labelTower_
const std::string labelTower_
Definition: HcalIsoTrkAnalyzer.cc:151
HcalIsoTrkAnalyzer::triggerEvent_
const edm::InputTag triggerEvent_
Definition: HcalIsoTrkAnalyzer.cc:149
HcalIsoTrkAnalyzer::hitEthrEE1_
const double hitEthrEE1_
Definition: HcalIsoTrkAnalyzer.cc:146
Frameworkfwd.h
HcalIsoTrkAnalyzer::t_eAnnular
double t_eAnnular
Definition: HcalIsoTrkAnalyzer.cc:192
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
CaloPropagateTrack.h
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:31
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:129
HcalIsoTrkAnalyzer::trackP
double trackP(const reco::Track *, const edm::Handle< reco::GenParticleCollection > &)
Definition: HcalIsoTrkAnalyzer.cc:1266
HcalIsoTrkAnalyzer::t_EventWeight
double t_EventWeight
Definition: HcalIsoTrkAnalyzer.cc:186
flavorHistoryFilter_cfi.dr
dr
Definition: flavorHistoryFilter_cfi.py:37
HcalIsoTrkAnalyzer
Definition: HcalIsoTrkAnalyzer.cc:87
CaloGeometry.h
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:137
HcalIsoTrkAnalyzer::unCorrect_
const bool unCorrect_
Definition: HcalIsoTrkAnalyzer.cc:145
HcalDDDRecConstants.h
HcalIsoTrkAnalyzer::t_gentrackP
double t_gentrackP
Definition: HcalIsoTrkAnalyzer.cc:191
HcalIsoTrkAnalyzer::l1Filter_
const std::string l1Filter_
Definition: HcalIsoTrkAnalyzer.cc:135
reco::deltaR
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
HcalRecHitCollections.h
HcalIsoTrkAnalyzer::collapseDepth_
const bool collapseDepth_
Definition: HcalIsoTrkAnalyzer.cc:145
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:195
L1TGlobalUtil.h
HcalIsoTrkAnalyzer::hltConfig_
HLTConfigProvider hltConfig_
Definition: HcalIsoTrkAnalyzer.cc:131
HcalIsoTrkAnalyzer::tree2
TTree * tree2
Definition: HcalIsoTrkAnalyzer.cc:182
HcalDDDRecConstants
Definition: HcalDDDRecConstants.h:23
HcalIsoTrkAnalyzer::analyze
void analyze(edm::Event const &, edm::EventSetup const &) override
Definition: HcalIsoTrkAnalyzer.cc:386
HcalIsoTrkAnalyzer::t_EventNo
unsigned int t_EventNo
Definition: HcalIsoTrkAnalyzer.cc:183
HcalIsoTrkAnalyzer::oldDepth_
std::vector< int > oldDepth_
Definition: HcalIsoTrkAnalyzer.cc:158
HcalIsoTrkAnalyzer::pTrackMin_
const double pTrackMin_
Definition: HcalIsoTrkAnalyzer.cc:138
HcalIsoTrkAnalyzer::etadist_
double etadist_
Definition: HcalIsoTrkAnalyzer.cc:159
reco::TrackBase::pz
double pz() const
z coordinate of momentum vector
Definition: TrackBase.h:646
HcalIsoTrkAnalyzer::t_rhoh
double t_rhoh
Definition: HcalIsoTrkAnalyzer.cc:193
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:200
HcalDetId::ietaAbs
constexpr int ietaAbs() const
get the absolute value of the cell ieta
Definition: HcalDetId.h:148
ParameterSet.h
dummy
Definition: DummySelector.h:38
HcalIsoTrkAnalyzer::beginJob
void beginJob() override
Definition: HcalIsoTrkAnalyzer.cc:692
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
HcalDetId::zside
constexpr int zside() const
get the z-side of the cell (1/-1)
Definition: HcalDetId.h:141
HcalIsoTrkAnalyzer::t_l3pt
double t_l3pt
Definition: HcalIsoTrkAnalyzer.cc:188
HcalIsoTrkAnalyzer::t_pt
double t_pt
Definition: HcalIsoTrkAnalyzer.cc:186
HcalIsoTrkAnalyzer::tok_recVtx_
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
Definition: HcalIsoTrkAnalyzer.cc:164
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
HcalIsoTrkAnalyzer::t_eHcal30
double t_eHcal30
Definition: HcalIsoTrkAnalyzer.cc:193
edm::Event
Definition: Event.h:73
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:140
MagneticField
Definition: MagneticField.h:19
edm::Log
Definition: MessageLogger.h:70
HcalIsoTrkAnalyzer::t_eMipDR
double t_eMipDR
Definition: HcalIsoTrkAnalyzer.cc:190
HcalIsoTrkAnalyzer::t_HitEnergies3
std::vector< double > * t_HitEnergies3
Definition: HcalIsoTrkAnalyzer.cc:197
HcalIsoTrkAnalyzer::endRun
void endRun(edm::Run const &, edm::EventSetup const &) override
Definition: HcalIsoTrkAnalyzer.cc:801
HcalIsoTrkAnalyzer::processName_
const std::string processName_
Definition: HcalIsoTrkAnalyzer.cc:135
trigger::TriggerObject::mass
float mass() const
Definition: TriggerObject.h:55
edm::InputTag
Definition: InputTag.h:15
label
const char * label
Definition: PFTauDecayModeTools.cc:11
HcalIsoTrkAnalyzer::labelHBHE_
const std::string labelHBHE_
Definition: HcalIsoTrkAnalyzer.cc:151
EcalChannelStatusRcd.h
HcalIsoTrkAnalyzer::labelEB_
const std::string labelEB_
Definition: HcalIsoTrkAnalyzer.cc:150
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:1339
RemoveAddSevLevel.flag
flag
Definition: RemoveAddSevLevel.py:116
HcalIsoTrkAnalyzer::prescaleHigh_
const int prescaleHigh_
Definition: HcalIsoTrkAnalyzer.cc:142
HcalIsoTrkAnalyzer::t_Tracks
int t_Tracks
Definition: HcalIsoTrkAnalyzer.cc:200
ValidationMatrix.hlt
hlt
Definition: ValidationMatrix.py:459
HcalIsoTrkAnalyzer::dataType_
const int dataType_
Definition: HcalIsoTrkAnalyzer.cc:143
HcalIsoTrkAnalyzer::hep17_
const bool hep17_
Definition: HcalIsoTrkAnalyzer.cc:153
HcalIsoTrkAnalyzer::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: HcalIsoTrkAnalyzer.cc:168
HcalIsoTrkAnalyzer::l2Filter_
const std::string l2Filter_
Definition: HcalIsoTrkAnalyzer.cc:136
Generated for CMSSW Reference Manual by doxygen 1.8.16