CMS 3D CMS Logo
IsolatedTracksNxN.cc
Go to the documentation of this file.
1 // -*- C++ -*
2 //
3 // Package: IsolatedParticles
4 // Class: IsolatedTracksNxN
5 //
13 //
14 // Original Author: Seema Sharma
15 // Created: Mon Aug 10 15:30:40 CST 2009
16 //
17 //
18 
19 // system include files
20 #include <cmath>
21 #include <map>
22 #include <memory>
23 #include <string>
24 #include <vector>
25 
26 // user include files
27 #include <Math/GenVector/VectorUtil.h>
28 
29 // root objects
30 #include "TROOT.h"
31 #include "TSystem.h"
32 #include "TFile.h"
33 #include "TH1F.h"
34 #include "TH2F.h"
35 #include "TProfile.h"
36 #include "TDirectory.h"
37 #include "TTree.h"
38 
39 #include "Calibration/IsolatedParticles/interface/CaloSimInfo.h"
40 #include "Calibration/IsolatedParticles/interface/CaloPropagateTrack.h"
41 #include "Calibration/IsolatedParticles/interface/ChargeIsolation.h"
42 #include "Calibration/IsolatedParticles/interface/eECALMatrix.h"
43 #include "Calibration/IsolatedParticles/interface/eHCALMatrix.h"
44 #include "Calibration/IsolatedParticles/interface/eHCALMatrixExtra.h"
45 #include "Calibration/IsolatedParticles/interface/FindCaloHit.h"
46 #include "Calibration/IsolatedParticles/interface/MatchingSimTrack.h"
47 
48 #include "CondFormats/DataRecord/interface/EcalChannelStatusRcd.h"
49 //L1 trigger Menus etc
50 #include "CondFormats/DataRecord/interface/L1GtTriggerMenuRcd.h"
51 #include "CondFormats/DataRecord/interface/L1GtTriggerMenuRcd.h"
52 #include "CondFormats/DataRecord/interface/L1GtTriggerMaskAlgoTrigRcd.h"
53 #include "CondFormats/DataRecord/interface/L1GtTriggerMaskTechTrigRcd.h"
54 #include "CondFormats/L1TObjects/interface/L1GtTriggerMask.h"
55 #include "CondFormats/L1TObjects/interface/L1GtTriggerMenu.h"
56 #include "CondFormats/L1TObjects/interface/L1GtTriggerMenuFwd.h"
57 #include "CondFormats/L1TObjects/interface/L1GtTriggerMenu.h"
58 
59 #include "DataFormats/Common/interface/Ref.h"
60 #include "DataFormats/Common/interface/Handle.h"
61 #include "DataFormats/Math/interface/Point3D.h"
62 #include "DataFormats/Candidate/interface/Candidate.h"
63 // muons and tracks
64 #include "DataFormats/TrackReco/interface/Track.h"
65 #include "DataFormats/TrackReco/interface/TrackFwd.h"
66 #include "DataFormats/TrackReco/interface/TrackBase.h"
67 #include "DataFormats/TrackReco/interface/HitPattern.h"
68 // Vertices
69 #include "DataFormats/VertexReco/interface/VertexFwd.h"
70 #include "DataFormats/VertexReco/interface/Vertex.h"
71 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
72 // Calorimeters
73 #include "DataFormats/DetId/interface/DetId.h"
74 #include "DataFormats/EcalDetId/interface/EBDetId.h"
75 #include "DataFormats/EcalDetId/interface/EEDetId.h"
76 #include "DataFormats/HcalDetId/interface/HcalDetId.h"
77 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
78 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
79 // Trigger
80 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutSetupFwd.h"
81 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutSetup.h"
82 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
83 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapRecord.h"
84 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapFwd.h"
85 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMap.h"
86 #include "L1Trigger/GlobalTriggerAnalyzer/interface/L1GtUtils.h"
87 //L1 objects
88 #include "DataFormats/L1Trigger/interface/L1JetParticle.h"
89 #include "DataFormats/L1Trigger/interface/L1JetParticleFwd.h"
90 #include "DataFormats/L1Trigger/interface/L1EmParticle.h"
91 #include "DataFormats/L1Trigger/interface/L1EmParticleFwd.h"
92 #include "DataFormats/L1Trigger/interface/L1MuonParticle.h"
93 #include "DataFormats/L1Trigger/interface/L1MuonParticleFwd.h"
94 // Jets in the event
95 #include "DataFormats/JetReco/interface/CaloJet.h"
96 #include "DataFormats/JetReco/interface/CaloJetCollection.h"
97 #include "DataFormats/JetReco/interface/JetExtendedAssociation.h"
98 
99 #include "DataFormats/Math/interface/deltaPhi.h"
100 #include "DataFormats/Math/interface/deltaR.h"
101 
102 #include "FWCore/Framework/interface/Frameworkfwd.h"
103 #include "FWCore/Framework/interface/one/EDAnalyzer.h"
104 #include "FWCore/Framework/interface/Event.h"
105 #include "FWCore/Framework/interface/EventSetup.h"
106 #include "FWCore/Framework/interface/MakerMacros.h"
107 #include "FWCore/MessageLogger/interface/MessageLogger.h"
108 #include "FWCore/ParameterSet/interface/ParameterSet.h"
109 // TFile Service
110 #include "FWCore/ServiceRegistry/interface/Service.h"
111 #include "CommonTools/UtilAlgos/interface/TFileService.h"
112 
113 // ecal / hcal
114 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
115 #include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
116 #include "Geometry/Records/interface/CaloTopologyRecord.h"
117 #include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
118 #include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"
119 #include "Geometry/CaloTopology/interface/HcalTopology.h"
120 #include "Geometry/CaloTopology/interface/CaloTopology.h"
121 
122 #include "MagneticField/Engine/interface/MagneticField.h"
123 #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
124 
125 #include "RecoCaloTools/Navigation/interface/CaloNavigator.h"
126 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
127 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"
128 
129 // SimHit + SimTrack
130 #include "SimDataFormats/CaloHit/interface/PCaloHitContainer.h"
131 #include "SimDataFormats/Track/interface/SimTrack.h"
132 #include "SimDataFormats/Track/interface/SimTrackContainer.h"
133 #include "SimDataFormats/Vertex/interface/SimVertex.h"
134 #include "SimDataFormats/Vertex/interface/SimVertexContainer.h"
135 
136 // track associator
137 #include "SimTracker/Records/interface/TrackAssociatorRecord.h"
138 #include "SimTracker/TrackerHitAssociation/interface/TrackerHitAssociator.h"
139 #include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
140 #include "TrackingTools/TrackAssociator/interface/TrackDetectorAssociator.h"
141 
142 class IsolatedTracksNxN : public edm::one::EDAnalyzer<edm::one::SharedResources> {
143 public:
144  explicit IsolatedTracksNxN(const edm::ParameterSet &);
145  ~IsolatedTracksNxN() override;
146 
147  static void fillDescriptions(edm::ConfigurationDescriptions &descriptions);
148 
149 private:
150  void beginJob() override;
151  void analyze(const edm::Event &, const edm::EventSetup &) override;
152  void endJob() override;
153 
154  void printTrack(const reco::Track *pTrack);
155 
156  void bookHistograms();
157 
158  void clearTreeVectors();
159 
160 private:
161  std::unique_ptr<L1GtUtils> m_l1GtUtils;
162  static constexpr size_t nL1BitsMax = 128;
163  TrackerHitAssociator::Config trackerHitAssociatorConfig_;
164 
165  // map of trig bit, algo name and num events passed
166  std::map<std::pair<unsigned int, std::string>, int> l1AlgoMap_;
167  std::vector<unsigned int> m_triggerMaskAlgoTrig;
168 
169  const bool doMC_, writeAllTracks_;
170  const int myverbose_;
171  const double pvTracksPtMin_;
172  const int debugTrks_;
173  const bool printTrkHitPattern_;
174  const double minTrackP_, maxTrackEta_;
175  const bool debugL1Info_, L1TriggerAlgoInfo_;
176  const double tMinE_, tMaxE_, tMinH_, tMaxH_;
177  bool initL1_;
178  int nEventProc_, nbad_;
179 
180  edm::EDGetTokenT<l1extra::L1JetParticleCollection> tok_L1extTauJet_;
181  edm::EDGetTokenT<l1extra::L1JetParticleCollection> tok_L1extCenJet_;
182  edm::EDGetTokenT<l1extra::L1JetParticleCollection> tok_L1extFwdJet_;
183 
184  edm::EDGetTokenT<l1extra::L1MuonParticleCollection> tok_L1extMu_;
185  edm::EDGetTokenT<l1extra::L1EmParticleCollection> tok_L1extIsoEm_;
186  edm::EDGetTokenT<l1extra::L1EmParticleCollection> tok_L1extNoIsoEm_;
187 
188  edm::EDGetTokenT<reco::CaloJetCollection> tok_jets_;
189  edm::EDGetTokenT<HBHERecHitCollection> tok_hbhe_;
190 
191  edm::EDGetTokenT<reco::TrackCollection> tok_genTrack_;
192  edm::EDGetTokenT<reco::VertexCollection> tok_recVtx_;
193  edm::EDGetTokenT<reco::BeamSpot> tok_bs_;
194 
195  edm::EDGetTokenT<EcalRecHitCollection> tok_EB_;
196  edm::EDGetTokenT<EcalRecHitCollection> tok_EE_;
197  edm::EDGetTokenT<edm::SimTrackContainer> tok_simTk_;
198  edm::EDGetTokenT<edm::SimVertexContainer> tok_simVtx_;
199  edm::EDGetTokenT<edm::PCaloHitContainer> tok_caloEB_;
200  edm::EDGetTokenT<edm::PCaloHitContainer> tok_caloEE_;
201  edm::EDGetTokenT<edm::PCaloHitContainer> tok_caloHH_;
202 
203  edm::ESGetToken<CaloGeometry, CaloGeometryRecord> tok_geom_;
204  edm::ESGetToken<CaloTopology, CaloTopologyRecord> tok_caloTopology_;
205  edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> tok_topo_;
206  edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> tok_magField_;
207  edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> tok_ecalChStatus_;
208  edm::ESGetToken<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd> tok_sevlv_;
209  edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> tok_htmap_;
210  static constexpr size_t NPBins = 15;
211  static constexpr size_t NEtaBins = 3;
212  double genPartPBins[NPBins + 1], genPartEtaBins[NEtaBins + 1];
213 
214  TH1F *h_maxNearP15x15[NPBins][NEtaBins], *h_maxNearP21x21[NPBins][NEtaBins], *h_maxNearP25x25[NPBins][NEtaBins],
215  *h_maxNearP31x31[NPBins][NEtaBins];
216  TH1I *h_L1AlgoNames;
217  TH1F *h_PVTracksWt;
218  TH1F *h_nTracks;
219  TH1F *h_recPt_0, *h_recP_0, *h_recEta_0, *h_recPhi_0;
220  TH2F *h_recEtaPt_0, *h_recEtaP_0;
221  TH1F *h_recPt_1, *h_recP_1, *h_recEta_1, *h_recPhi_1;
222  TH2F *h_recEtaPt_1, *h_recEtaP_1;
223  TH1F *h_recPt_2, *h_recP_2, *h_recEta_2, *h_recPhi_2;
224  TH2F *h_recEtaPt_2, *h_recEtaP_2;
225 
226  TTree *tree_;
227 
228  int t_nTracks;
229  int t_RunNo, t_EvtNo, t_Lumi, t_Bunch;
230  std::vector<std::string> *t_L1AlgoNames;
231  std::vector<int> *t_L1PreScale;
232  int t_L1Decision[128];
233 
234  std::vector<double> *t_L1CenJetPt, *t_L1CenJetEta, *t_L1CenJetPhi;
235  std::vector<double> *t_L1FwdJetPt, *t_L1FwdJetEta, *t_L1FwdJetPhi;
236  std::vector<double> *t_L1TauJetPt, *t_L1TauJetEta, *t_L1TauJetPhi;
237  std::vector<double> *t_L1MuonPt, *t_L1MuonEta, *t_L1MuonPhi;
238  std::vector<double> *t_L1IsoEMPt, *t_L1IsoEMEta, *t_L1IsoEMPhi;
239  std::vector<double> *t_L1NonIsoEMPt, *t_L1NonIsoEMEta, *t_L1NonIsoEMPhi;
240  std::vector<double> *t_L1METPt, *t_L1METEta, *t_L1METPhi;
241 
242  std::vector<double> *t_PVx, *t_PVy, *t_PVz, *t_PVTracksSumPt;
243  std::vector<double> *t_PVTracksSumPtWt, *t_PVTracksSumPtHP, *t_PVTracksSumPtHPWt;
244  std::vector<int> *t_PVisValid, *t_PVNTracks, *t_PVNTracksWt, *t_PVndof;
245  std::vector<int> *t_PVNTracksHP, *t_PVNTracksHPWt;
246 
247  std::vector<double> *t_jetPt, *t_jetEta, *t_jetPhi;
248  std::vector<double> *t_nTrksJetCalo, *t_nTrksJetVtx;
249 
250  std::vector<double> *t_trackPAll, *t_trackEtaAll, *t_trackPhiAll, *t_trackPdgIdAll;
251  std::vector<double> *t_trackPtAll;
252  std::vector<double> *t_trackDxyAll, *t_trackDzAll, *t_trackDxyPVAll, *t_trackDzPVAll, *t_trackChiSqAll;
253 
254  std::vector<double> *t_trackP, *t_trackPt, *t_trackEta, *t_trackPhi;
255  std::vector<double> *t_trackEcalEta, *t_trackEcalPhi, *t_trackHcalEta, *t_trackHcalPhi;
256  std::vector<double> *t_trackDxy, *t_trackDxyBS, *t_trackDz, *t_trackDzBS;
257  std::vector<double> *t_trackDxyPV, *t_trackDzPV;
258  std::vector<double> *t_trackChiSq;
259  std::vector<int> *t_trackPVIdx;
260 
261  std::vector<int> *t_NLayersCrossed, *t_trackNOuterHits;
262  std::vector<int> *t_trackHitsTOB, *t_trackHitsTEC;
263  std::vector<int> *t_trackHitInMissTOB, *t_trackHitInMissTEC, *t_trackHitInMissTIB, *t_trackHitInMissTID,
264  *t_trackHitInMissTIBTID;
265  std::vector<int> *t_trackHitOutMissTOB, *t_trackHitOutMissTEC, *t_trackHitOutMissTIB, *t_trackHitOutMissTID,
266  *t_trackHitOutMissTOBTEC;
267  std::vector<int> *t_trackHitInMeasTOB, *t_trackHitInMeasTEC, *t_trackHitInMeasTIB, *t_trackHitInMeasTID;
268  std::vector<int> *t_trackHitOutMeasTOB, *t_trackHitOutMeasTEC, *t_trackHitOutMeasTIB, *t_trackHitOutMeasTID;
269  std::vector<double> *t_trackOutPosOutHitDr, *t_trackL;
270 
271  std::vector<double> *t_maxNearP31x31;
272  std::vector<double> *t_maxNearP21x21;
273  std::vector<int> *t_ecalSpike11x11;
274 
275  std::vector<double> *t_e7x7, *t_e9x9, *t_e11x11, *t_e15x15;
276  std::vector<double> *t_e7x7_10Sig, *t_e9x9_10Sig, *t_e11x11_10Sig, *t_e15x15_10Sig;
277  std::vector<double> *t_e7x7_15Sig, *t_e9x9_15Sig, *t_e11x11_15Sig, *t_e15x15_15Sig;
278  std::vector<double> *t_e7x7_20Sig, *t_e9x9_20Sig, *t_e11x11_20Sig, *t_e15x15_20Sig;
279  std::vector<double> *t_e7x7_25Sig, *t_e9x9_25Sig, *t_e11x11_25Sig, *t_e15x15_25Sig;
280  std::vector<double> *t_e7x7_30Sig, *t_e9x9_30Sig, *t_e11x11_30Sig, *t_e15x15_30Sig;
281 
282  std::vector<double> *t_esimPdgId, *t_simTrackP, *t_trkEcalEne;
283  std::vector<double> *t_esim7x7, *t_esim9x9, *t_esim11x11, *t_esim15x15;
284  std::vector<double> *t_esim7x7Matched, *t_esim9x9Matched, *t_esim11x11Matched, *t_esim15x15Matched;
285  std::vector<double> *t_esim7x7Rest, *t_esim9x9Rest, *t_esim11x11Rest, *t_esim15x15Rest;
286  std::vector<double> *t_esim7x7Photon, *t_esim9x9Photon, *t_esim11x11Photon, *t_esim15x15Photon;
287  std::vector<double> *t_esim7x7NeutHad, *t_esim9x9NeutHad, *t_esim11x11NeutHad, *t_esim15x15NeutHad;
288  std::vector<double> *t_esim7x7CharHad, *t_esim9x9CharHad, *t_esim11x11CharHad, *t_esim15x15CharHad;
289 
290  std::vector<double> *t_maxNearHcalP3x3, *t_maxNearHcalP5x5, *t_maxNearHcalP7x7;
291  std::vector<double> *t_h3x3, *t_h5x5, *t_h7x7;
292  std::vector<double> *t_h3x3Sig, *t_h5x5Sig, *t_h7x7Sig;
293  std::vector<int> *t_infoHcal;
294 
295  std::vector<double> *t_trkHcalEne;
296  std::vector<double> *t_hsim3x3, *t_hsim5x5, *t_hsim7x7;
297  std::vector<double> *t_hsim3x3Matched, *t_hsim5x5Matched, *t_hsim7x7Matched;
298  std::vector<double> *t_hsim3x3Rest, *t_hsim5x5Rest, *t_hsim7x7Rest;
299  std::vector<double> *t_hsim3x3Photon, *t_hsim5x5Photon, *t_hsim7x7Photon;
300  std::vector<double> *t_hsim3x3NeutHad, *t_hsim5x5NeutHad, *t_hsim7x7NeutHad;
301  std::vector<double> *t_hsim3x3CharHad, *t_hsim5x5CharHad, *t_hsim7x7CharHad;
302 };
303 
304 static const bool useL1EventSetup(true);
305 static const bool useL1GtTriggerMenuLite(true);
306 
307 IsolatedTracksNxN::IsolatedTracksNxN(const edm::ParameterSet &iConfig)
308  : trackerHitAssociatorConfig_(consumesCollector()),
309  doMC_(iConfig.getUntrackedParameter<bool>("doMC", false)),
310  writeAllTracks_(iConfig.getUntrackedParameter<bool>("writeAllTracks", false)),
311  myverbose_(iConfig.getUntrackedParameter<int>("verbosity", 5)),
312  pvTracksPtMin_(iConfig.getUntrackedParameter<double>("pvTracksPtMin", 1.0)),
313  debugTrks_(iConfig.getUntrackedParameter<int>("debugTracks", 0)),
314  printTrkHitPattern_(iConfig.getUntrackedParameter<bool>("printTrkHitPattern", false)),
315  minTrackP_(iConfig.getUntrackedParameter<double>("minTrackP", 1.0)),
316  maxTrackEta_(iConfig.getUntrackedParameter<double>("maxTrackEta", 5.0)),
317  debugL1Info_(iConfig.getUntrackedParameter<bool>("debugL1Info", false)),
318  L1TriggerAlgoInfo_(iConfig.getUntrackedParameter<bool>("l1TriggerAlgoInfo", false)),
319  tMinE_(iConfig.getUntrackedParameter<double>("timeMinCutECAL", -500.)),
320  tMaxE_(iConfig.getUntrackedParameter<double>("timeMaxCutECAL", 500.)),
321  tMinH_(iConfig.getUntrackedParameter<double>("timeMinCutHCAL", -500.)),
322  tMaxH_(iConfig.getUntrackedParameter<double>("timeMaxCutHCAL", 500.)),
323  t_L1AlgoNames(nullptr),
324  t_L1PreScale(nullptr),
325  t_L1CenJetPt(nullptr),
326  t_L1CenJetEta(nullptr),
327  t_L1CenJetPhi(nullptr),
328  t_L1FwdJetPt(nullptr),
329  t_L1FwdJetEta(nullptr),
330  t_L1FwdJetPhi(nullptr),
331  t_L1TauJetPt(nullptr),
332  t_L1TauJetEta(nullptr),
333  t_L1TauJetPhi(nullptr),
334  t_L1MuonPt(nullptr),
335  t_L1MuonEta(nullptr),
336  t_L1MuonPhi(nullptr),
337  t_L1IsoEMPt(nullptr),
338  t_L1IsoEMEta(nullptr),
339  t_L1IsoEMPhi(nullptr),
340  t_L1NonIsoEMPt(nullptr),
341  t_L1NonIsoEMEta(nullptr),
342  t_L1NonIsoEMPhi(nullptr),
343  t_L1METPt(nullptr),
344  t_L1METEta(nullptr),
345  t_L1METPhi(nullptr),
346  t_PVx(nullptr),
347  t_PVy(nullptr),
348  t_PVz(nullptr),
349  t_PVTracksSumPt(nullptr),
350  t_PVTracksSumPtWt(nullptr),
351  t_PVTracksSumPtHP(nullptr),
352  t_PVTracksSumPtHPWt(nullptr),
353  t_PVisValid(nullptr),
354  t_PVNTracks(nullptr),
355  t_PVNTracksWt(nullptr),
356  t_PVndof(nullptr),
357  t_PVNTracksHP(nullptr),
358  t_PVNTracksHPWt(nullptr),
359  t_jetPt(nullptr),
360  t_jetEta(nullptr),
361  t_jetPhi(nullptr),
362  t_nTrksJetCalo(nullptr),
363  t_nTrksJetVtx(nullptr),
364  t_trackPAll(nullptr),
365  t_trackEtaAll(nullptr),
366  t_trackPhiAll(nullptr),
367  t_trackPdgIdAll(nullptr),
368  t_trackPtAll(nullptr),
369  t_trackDxyAll(nullptr),
370  t_trackDzAll(nullptr),
371  t_trackDxyPVAll(nullptr),
372  t_trackDzPVAll(nullptr),
373  t_trackChiSqAll(nullptr),
374  t_trackP(nullptr),
375  t_trackPt(nullptr),
376  t_trackEta(nullptr),
377  t_trackPhi(nullptr),
378  t_trackEcalEta(nullptr),
379  t_trackEcalPhi(nullptr),
380  t_trackHcalEta(nullptr),
381  t_trackHcalPhi(nullptr),
382  t_trackDxy(nullptr),
383  t_trackDxyBS(nullptr),
384  t_trackDz(nullptr),
385  t_trackDzBS(nullptr),
386  t_trackDxyPV(nullptr),
387  t_trackDzPV(nullptr),
388  t_trackChiSq(nullptr),
389  t_trackPVIdx(nullptr),
390  t_NLayersCrossed(nullptr),
391  t_trackNOuterHits(nullptr),
392  t_trackHitsTOB(nullptr),
393  t_trackHitsTEC(nullptr),
394  t_trackHitInMissTOB(nullptr),
395  t_trackHitInMissTEC(nullptr),
396  t_trackHitInMissTIB(nullptr),
397  t_trackHitInMissTID(nullptr),
398  t_trackHitInMissTIBTID(nullptr),
399  t_trackHitOutMissTOB(nullptr),
400  t_trackHitOutMissTEC(nullptr),
401  t_trackHitOutMissTIB(nullptr),
402  t_trackHitOutMissTID(nullptr),
403  t_trackHitOutMissTOBTEC(nullptr),
404  t_trackHitInMeasTOB(nullptr),
405  t_trackHitInMeasTEC(nullptr),
406  t_trackHitInMeasTIB(nullptr),
407  t_trackHitInMeasTID(nullptr),
408  t_trackHitOutMeasTOB(nullptr),
409  t_trackHitOutMeasTEC(nullptr),
410  t_trackHitOutMeasTIB(nullptr),
411  t_trackHitOutMeasTID(nullptr),
412  t_trackOutPosOutHitDr(nullptr),
413  t_trackL(nullptr),
414  t_maxNearP31x31(nullptr),
415  t_maxNearP21x21(nullptr),
416  t_ecalSpike11x11(nullptr),
417  t_e7x7(nullptr),
418  t_e9x9(nullptr),
419  t_e11x11(nullptr),
420  t_e15x15(nullptr),
421  t_e7x7_10Sig(nullptr),
422  t_e9x9_10Sig(nullptr),
423  t_e11x11_10Sig(nullptr),
424  t_e15x15_10Sig(nullptr),
425  t_e7x7_15Sig(nullptr),
426  t_e9x9_15Sig(nullptr),
427  t_e11x11_15Sig(nullptr),
428  t_e15x15_15Sig(nullptr),
429  t_e7x7_20Sig(nullptr),
430  t_e9x9_20Sig(nullptr),
431  t_e11x11_20Sig(nullptr),
432  t_e15x15_20Sig(nullptr),
433  t_e7x7_25Sig(nullptr),
434  t_e9x9_25Sig(nullptr),
435  t_e11x11_25Sig(nullptr),
436  t_e15x15_25Sig(nullptr),
437  t_e7x7_30Sig(nullptr),
438  t_e9x9_30Sig(nullptr),
439  t_e11x11_30Sig(nullptr),
440  t_e15x15_30Sig(nullptr),
441  t_esimPdgId(nullptr),
442  t_simTrackP(nullptr),
443  t_trkEcalEne(nullptr),
444  t_esim7x7(nullptr),
445  t_esim9x9(nullptr),
446  t_esim11x11(nullptr),
447  t_esim15x15(nullptr),
448  t_esim7x7Matched(nullptr),
449  t_esim9x9Matched(nullptr),
450  t_esim11x11Matched(nullptr),
451  t_esim15x15Matched(nullptr),
452  t_esim7x7Rest(nullptr),
453  t_esim9x9Rest(nullptr),
454  t_esim11x11Rest(nullptr),
455  t_esim15x15Rest(nullptr),
456  t_esim7x7Photon(nullptr),
457  t_esim9x9Photon(nullptr),
458  t_esim11x11Photon(nullptr),
459  t_esim15x15Photon(nullptr),
460  t_esim7x7NeutHad(nullptr),
461  t_esim9x9NeutHad(nullptr),
462  t_esim11x11NeutHad(nullptr),
463  t_esim15x15NeutHad(nullptr),
464  t_esim7x7CharHad(nullptr),
465  t_esim9x9CharHad(nullptr),
466  t_esim11x11CharHad(nullptr),
467  t_esim15x15CharHad(nullptr),
468  t_maxNearHcalP3x3(nullptr),
469  t_maxNearHcalP5x5(nullptr),
470  t_maxNearHcalP7x7(nullptr),
471  t_h3x3(nullptr),
472  t_h5x5(nullptr),
473  t_h7x7(nullptr),
474  t_h3x3Sig(nullptr),
475  t_h5x5Sig(nullptr),
476  t_h7x7Sig(nullptr),
477  t_infoHcal(nullptr),
478  t_trkHcalEne(nullptr),
479  t_hsim3x3(nullptr),
480  t_hsim5x5(nullptr),
481  t_hsim7x7(nullptr),
482  t_hsim3x3Matched(nullptr),
483  t_hsim5x5Matched(nullptr),
484  t_hsim7x7Matched(nullptr),
485  t_hsim3x3Rest(nullptr),
486  t_hsim5x5Rest(nullptr),
487  t_hsim7x7Rest(nullptr),
488  t_hsim3x3Photon(nullptr),
489  t_hsim5x5Photon(nullptr),
490  t_hsim7x7Photon(nullptr),
491  t_hsim3x3NeutHad(nullptr),
492  t_hsim5x5NeutHad(nullptr),
493  t_hsim7x7NeutHad(nullptr),
494  t_hsim3x3CharHad(nullptr),
495  t_hsim5x5CharHad(nullptr),
496  t_hsim7x7CharHad(nullptr) {
497  if (L1TriggerAlgoInfo_) {
498  m_l1GtUtils = std::make_unique<L1GtUtils>(
499  iConfig, consumesCollector(), useL1GtTriggerMenuLite, *this, L1GtUtils::UseEventSetupIn::Event);
500  }
501 
502  usesResource(TFileService::kSharedResource);
503 
504  //now do what ever initialization is needed
505 
506  edm::InputTag L1extraTauJetSource_ = iConfig.getParameter<edm::InputTag>("l1extraTauJetSource");
507  edm::InputTag L1extraCenJetSource_ = iConfig.getParameter<edm::InputTag>("l1extraCenJetSource");
508  edm::InputTag L1extraFwdJetSource_ = iConfig.getParameter<edm::InputTag>("l1extraFwdJetSource");
509  edm::InputTag L1extraMuonSource_ = iConfig.getParameter<edm::InputTag>("l1extraMuonSource");
510  edm::InputTag L1extraIsoEmSource_ = iConfig.getParameter<edm::InputTag>("l1extraIsoEmSource");
511  edm::InputTag L1extraNonIsoEmSource_ = iConfig.getParameter<edm::InputTag>("l1extraNonIsoEmSource");
512  edm::InputTag L1GTReadoutRcdSource_ = iConfig.getParameter<edm::InputTag>("l1GTReadoutRcdSource");
513  edm::InputTag L1GTObjectMapRcdSource_ = iConfig.getParameter<edm::InputTag>("l1GTObjectMapRcdSource");
514  edm::InputTag JetSrc_ = iConfig.getParameter<edm::InputTag>("jetSource");
515  edm::InputTag JetExtender_ = iConfig.getParameter<edm::InputTag>("jetExtender");
516  edm::InputTag HBHERecHitSource_ = iConfig.getParameter<edm::InputTag>("hbheRecHitSource");
517 
518  // define tokens for access
519  tok_L1extTauJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraTauJetSource_);
520  tok_L1extCenJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraCenJetSource_);
521  tok_L1extFwdJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraFwdJetSource_);
522  tok_L1extMu_ = consumes<l1extra::L1MuonParticleCollection>(L1extraMuonSource_);
523  tok_L1extIsoEm_ = consumes<l1extra::L1EmParticleCollection>(L1extraIsoEmSource_);
524  tok_L1extNoIsoEm_ = consumes<l1extra::L1EmParticleCollection>(L1extraNonIsoEmSource_);
525  tok_jets_ = consumes<reco::CaloJetCollection>(JetSrc_);
526  tok_hbhe_ = consumes<HBHERecHitCollection>(HBHERecHitSource_);
527 
528  tok_genTrack_ = consumes<reco::TrackCollection>(edm::InputTag("generalTracks"));
529  tok_recVtx_ = consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"));
530  tok_bs_ = consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"));
531  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
532  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
533 
534  tok_simTk_ = consumes<edm::SimTrackContainer>(edm::InputTag("g4SimHits"));
535  tok_simVtx_ = consumes<edm::SimVertexContainer>(edm::InputTag("g4SimHits"));
536  tok_caloEB_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "EcalHitsEB"));
537  tok_caloEE_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "EcalHitsEE"));
538  tok_caloHH_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "HcalHits"));
539 
540  if (myverbose_ >= 0) {
541  edm::LogVerbatim("IsoTrack") << "Parameters read from config file \n"
542  << " doMC " << doMC_ << "\t myverbose " << myverbose_ << "\t minTrackP "
543  << minTrackP_ << "\t maxTrackEta " << maxTrackEta_ << "\t tMinE " << tMinE_
544  << "\t tMaxE " << tMaxE_ << "\t tMinH " << tMinH_ << "\t tMaxH " << tMaxH_
545  << "\n debugL1Info " << debugL1Info_ << "\t L1TriggerAlgoInfo " << L1TriggerAlgoInfo_
546  << "\n";
547  }
548 
549  tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
550  tok_caloTopology_ = esConsumes<CaloTopology, CaloTopologyRecord>();
551  tok_topo_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
552  tok_magField_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
553  tok_ecalChStatus_ = esConsumes<EcalChannelStatus, EcalChannelStatusRcd>();
554  tok_sevlv_ = esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>();
555  tok_htmap_ = esConsumes<EcalTrigTowerConstituentsMap, IdealGeometryRecord>();
556 }
557 
558 IsolatedTracksNxN::~IsolatedTracksNxN() {
559  delete t_PVx;
560  delete t_PVy;
561  delete t_PVz;
562  delete t_PVisValid;
563  delete t_PVndof;
564  delete t_PVNTracks;
565  delete t_PVNTracksWt;
566  delete t_PVTracksSumPt;
567  delete t_PVTracksSumPtWt;
568  delete t_PVNTracksHP;
569  delete t_PVNTracksHPWt;
570  delete t_PVTracksSumPtHP;
571  delete t_PVTracksSumPtHPWt;
572  delete t_L1AlgoNames;
573  delete t_L1PreScale;
574  delete t_L1CenJetPt;
575  delete t_L1CenJetEta;
576  delete t_L1CenJetPhi;
577  delete t_L1FwdJetPt;
578  delete t_L1FwdJetEta;
579  delete t_L1FwdJetPhi;
580  delete t_L1TauJetPt;
581  delete t_L1TauJetEta;
582  delete t_L1TauJetPhi;
583  delete t_L1MuonPt;
584  delete t_L1MuonEta;
585  delete t_L1MuonPhi;
586  delete t_L1IsoEMPt;
587  delete t_L1IsoEMEta;
588  delete t_L1IsoEMPhi;
589  delete t_L1NonIsoEMPt;
590  delete t_L1NonIsoEMEta;
591  delete t_L1NonIsoEMPhi;
592  delete t_L1METPt;
593  delete t_L1METEta;
594  delete t_L1METPhi;
595  delete t_jetPt;
596  delete t_jetEta;
597  delete t_jetPhi;
598  delete t_nTrksJetCalo;
599  delete t_nTrksJetVtx;
600  delete t_trackPAll;
601  delete t_trackEtaAll;
602  delete t_trackPhiAll;
603  delete t_trackPdgIdAll;
604  delete t_trackPtAll;
605  delete t_trackDxyAll;
606  delete t_trackDzAll;
607  delete t_trackDxyPVAll;
608  delete t_trackDzPVAll;
609  delete t_trackChiSqAll;
610  delete t_trackP;
611  delete t_trackPt;
612  delete t_trackEta;
613  delete t_trackPhi;
614  delete t_trackEcalEta;
615  delete t_trackEcalPhi;
616  delete t_trackHcalEta;
617  delete t_trackHcalPhi;
618  delete t_trackNOuterHits;
619  delete t_NLayersCrossed;
620  delete t_trackDxy;
621  delete t_trackDxyBS;
622  delete t_trackDz;
623  delete t_trackDzBS;
624  delete t_trackDxyPV;
625  delete t_trackDzPV;
626  delete t_trackPVIdx;
627  delete t_trackChiSq;
628  delete t_trackHitsTOB;
629  delete t_trackHitsTEC;
630  delete t_trackHitInMissTOB;
631  delete t_trackHitInMissTEC;
632  delete t_trackHitInMissTIB;
633  delete t_trackHitInMissTID;
634  delete t_trackHitInMissTIBTID;
635  delete t_trackHitOutMissTOB;
636  delete t_trackHitOutMissTEC;
637  delete t_trackHitOutMissTIB;
638  delete t_trackHitOutMissTID;
639  delete t_trackHitOutMissTOBTEC;
640  delete t_trackHitInMeasTOB;
641  delete t_trackHitInMeasTEC;
642  delete t_trackHitInMeasTIB;
643  delete t_trackHitInMeasTID;
644  delete t_trackHitOutMeasTOB;
645  delete t_trackHitOutMeasTEC;
646  delete t_trackHitOutMeasTIB;
647  delete t_trackHitOutMeasTID;
648  delete t_trackOutPosOutHitDr;
649  delete t_trackL;
650  delete t_maxNearP31x31;
651  delete t_maxNearP21x21;
652  delete t_ecalSpike11x11;
653  delete t_e7x7;
654  delete t_e9x9;
655  delete t_e11x11;
656  delete t_e15x15;
657  delete t_e7x7_10Sig;
658  delete t_e9x9_10Sig;
659  delete t_e11x11_10Sig;
660  delete t_e15x15_10Sig;
661  delete t_e7x7_15Sig;
662  delete t_e9x9_15Sig;
663  delete t_e11x11_15Sig;
664  delete t_e15x15_15Sig;
665  delete t_e7x7_20Sig;
666  delete t_e9x9_20Sig;
667  delete t_e11x11_20Sig;
668  delete t_e15x15_20Sig;
669  delete t_e7x7_25Sig;
670  delete t_e9x9_25Sig;
671  delete t_e11x11_25Sig;
672  delete t_e15x15_25Sig;
673  delete t_e7x7_30Sig;
674  delete t_e9x9_30Sig;
675  delete t_e11x11_30Sig;
676  delete t_e15x15_30Sig;
677  delete t_esim7x7;
678  delete t_esim9x9;
679  delete t_esim11x11;
680  delete t_esim15x15;
681  delete t_esim7x7Matched;
682  delete t_esim9x9Matched;
683  delete t_esim11x11Matched;
684  delete t_esim15x15Matched;
685  delete t_esim7x7Rest;
686  delete t_esim9x9Rest;
687  delete t_esim11x11Rest;
688  delete t_esim15x15Rest;
689  delete t_esim7x7Photon;
690  delete t_esim9x9Photon;
691  delete t_esim11x11Photon;
692  delete t_esim15x15Photon;
693  delete t_esim7x7NeutHad;
694  delete t_esim9x9NeutHad;
695  delete t_esim11x11NeutHad;
696  delete t_esim15x15NeutHad;
697  delete t_esim7x7CharHad;
698  delete t_esim9x9CharHad;
699  delete t_esim11x11CharHad;
700  delete t_esim15x15CharHad;
701  delete t_trkEcalEne;
702  delete t_simTrackP;
703  delete t_esimPdgId;
704  delete t_maxNearHcalP3x3;
705  delete t_maxNearHcalP5x5;
706  delete t_maxNearHcalP7x7;
707  delete t_h3x3;
708  delete t_h5x5;
709  delete t_h7x7;
710  delete t_h3x3Sig;
711  delete t_h5x5Sig;
712  delete t_h7x7Sig;
713  delete t_infoHcal;
714  delete t_trkHcalEne;
715  delete t_hsim3x3;
716  delete t_hsim5x5;
717  delete t_hsim7x7;
718  delete t_hsim3x3Matched;
719  delete t_hsim5x5Matched;
720  delete t_hsim7x7Matched;
721  delete t_hsim3x3Rest;
722  delete t_hsim5x5Rest;
723  delete t_hsim7x7Rest;
724  delete t_hsim3x3Photon;
725  delete t_hsim5x5Photon;
726  delete t_hsim7x7Photon;
727  delete t_hsim3x3NeutHad;
728  delete t_hsim5x5NeutHad;
729  delete t_hsim7x7NeutHad;
730  delete t_hsim3x3CharHad;
731  delete t_hsim5x5CharHad;
732  delete t_hsim7x7CharHad;
733 }
734 
735 void IsolatedTracksNxN::fillDescriptions(edm::ConfigurationDescriptions &descriptions) {
736  edm::ParameterSetDescription desc;
737  desc.addUntracked<bool>("doMC", false);
738  desc.addUntracked<bool>("writeAllTracks", false);
739  desc.addUntracked<int>("verbosity", 1);
740  desc.addUntracked<double>("pvTracksPtMin", 0.200);
741  desc.addUntracked<int>("debugTracks", 0);
742  desc.addUntracked<bool>("printTrkHitPattern", true);
743  desc.addUntracked<double>("minTrackP", 1.0);
744  desc.addUntracked<double>("maxTrackEta", 2.6);
745  desc.addUntracked<bool>("debugL1Info", false);
746  desc.addUntracked<bool>("l1TriggerAlgoInfo", false);
747  desc.add<edm::InputTag>("l1extraTauJetSource", edm::InputTag("l1extraParticles", "Tau"));
748  desc.add<edm::InputTag>("l1extraCenJetSource", edm::InputTag("l1extraParticles", "Central"));
749  desc.add<edm::InputTag>("l1extraFwdJetSource", edm::InputTag("l1extraParticles", "Forward"));
750  desc.add<edm::InputTag>("l1extraMuonSource", edm::InputTag("l1extraParticles"));
751  desc.add<edm::InputTag>("l1extraIsoEmSource", edm::InputTag("l1extraParticles", "Isolated"));
752  desc.add<edm::InputTag>("l1extraNonIsoEmSource", edm::InputTag("l1extraParticles", "NonIsolated"));
753  desc.add<edm::InputTag>("l1GTReadoutRcdSource", edm::InputTag("gtDigis"));
754  desc.add<edm::InputTag>("l1GTObjectMapRcdSource", edm::InputTag("hltL1GtObjectMap"));
755  desc.add<edm::InputTag>("jetSource", edm::InputTag("iterativeCone5CaloJets"));
756  desc.add<edm::InputTag>("jetExtender", edm::InputTag("iterativeCone5JetExtender"));
757  desc.add<edm::InputTag>("hbheRecHitSource", edm::InputTag("hbhereco"));
758  desc.addUntracked<double>("maxNearTrackPT", 1.0);
759  desc.addUntracked<double>("timeMinCutECAL", -500.0);
760  desc.addUntracked<double>("timeMaxCutECAL", 500.0);
761  desc.addUntracked<double>("timeMinCutHCAL", -500.0);
762  desc.addUntracked<double>("timeMaxCutHCAL", 500.0);
763  descriptions.add("isolatedTracksNxN", desc);
764 }
765 
766 void IsolatedTracksNxN::analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup) {
767  bool haveIsoTrack = false;
768 
769  const MagneticField *bField = &iSetup.getData(tok_magField_);
770 
771  clearTreeVectors();
772 
773  t_RunNo = iEvent.id().run();
774  t_EvtNo = iEvent.id().event();
775  t_Lumi = iEvent.luminosityBlock();
776  t_Bunch = iEvent.bunchCrossing();
777 
778  ++nEventProc_;
779 
780  edm::Handle<reco::TrackCollection> trkCollection;
781  iEvent.getByToken(tok_genTrack_, trkCollection);
782  if (debugTrks_ > 1) {
783  edm::LogVerbatim("IsoTrack") << "Track Collection: ";
784  edm::LogVerbatim("IsoTrack") << "Number of Tracks " << trkCollection->size();
785  }
786  std::string theTrackQuality = "highPurity";
787  reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality);
788 
789  //===================== save L1 Trigger information =======================
790  if (L1TriggerAlgoInfo_) {
791  m_l1GtUtils->getL1GtRunCache(iEvent, iSetup, useL1EventSetup, useL1GtTriggerMenuLite);
792 
793  int iErrorCode = -1;
794  int l1ConfCode = -1;
795  const bool l1Conf = m_l1GtUtils->availableL1Configuration(iErrorCode, l1ConfCode);
796  if (!l1Conf) {
797  edm::LogVerbatim("IsoTrack")
798  << "\nL1 configuration code:" << l1ConfCode << "\nNo valid L1 trigger configuration available."
799  << "\nSee text above for error code interpretation"
800  << "\nNo return here, in order to test each method, protected against configuration error.";
801  }
802 
803  const L1GtTriggerMenu *m_l1GtMenu = m_l1GtUtils->ptrL1TriggerMenuEventSetup(iErrorCode);
804  const AlgorithmMap &algorithmMap = m_l1GtMenu->gtAlgorithmMap();
805  const std::string &menuName = m_l1GtMenu->gtTriggerMenuName();
806 
807  if (!initL1_) {
808  initL1_ = true;
809  edm::LogVerbatim("IsoTrack") << "menuName " << menuName;
810  for (CItAlgo itAlgo = algorithmMap.begin(); itAlgo != algorithmMap.end(); itAlgo++) {
811  std::string algName = itAlgo->first;
812  int algBitNumber = (itAlgo->second).algoBitNumber();
813  l1AlgoMap_.insert(std::pair<std::pair<unsigned int, std::string>, int>(
814  std::pair<unsigned int, std::string>(algBitNumber, algName), 0));
815  }
816  std::map<std::pair<unsigned int, std::string>, int>::iterator itr;
817  for (itr = l1AlgoMap_.begin(); itr != l1AlgoMap_.end(); itr++) {
818  edm::LogVerbatim("IsoTrack") << " ********** " << (itr->first).first << " " << (itr->first).second << " "
819  << itr->second;
820  }
821  }
822 
823  std::vector<int> algbits;
824  for (CItAlgo itAlgo = algorithmMap.begin(); itAlgo != algorithmMap.end(); itAlgo++) {
825  std::string algName = itAlgo->first;
826  int algBitNumber = (itAlgo->second).algoBitNumber();
827  bool decision = m_l1GtUtils->decision(iEvent, itAlgo->first, iErrorCode);
828  int preScale = m_l1GtUtils->prescaleFactor(iEvent, itAlgo->first, iErrorCode);
829 
830  // save the algo names which fired
831  if (decision) {
832  l1AlgoMap_[std::pair<unsigned int, std::string>(algBitNumber, algName)] += 1;
833  h_L1AlgoNames->Fill(algBitNumber);
834  t_L1AlgoNames->push_back(itAlgo->first);
835  t_L1PreScale->push_back(preScale);
836  t_L1Decision[algBitNumber] = 1;
837  algbits.push_back(algBitNumber);
838  }
839  }
840 
841  if (debugL1Info_) {
842  for (unsigned int ii = 0; ii < t_L1AlgoNames->size(); ii++) {
843  edm::LogVerbatim("IsoTrack") << ii << " " << (*t_L1AlgoNames)[ii] << " " << (*t_L1PreScale)[ii] << " "
844  << algbits[ii];
845  }
846  for (int i = 0; i < 128; ++i) {
847  edm::LogVerbatim("IsoTrack") << "L1Decision: " << i << ":" << t_L1Decision[i];
848  }
849  }
850 
851  // L1Taus
852  edm::Handle<l1extra::L1JetParticleCollection> l1TauHandle;
853  iEvent.getByToken(tok_L1extTauJet_, l1TauHandle);
854  l1extra::L1JetParticleCollection::const_iterator itr;
855  int iL1Obj = 0;
856  for (itr = l1TauHandle->begin(), iL1Obj = 0; itr != l1TauHandle->end(); ++itr, iL1Obj++) {
857  if (iL1Obj < 1) {
858  t_L1TauJetPt->push_back(itr->pt());
859  t_L1TauJetEta->push_back(itr->eta());
860  t_L1TauJetPhi->push_back(itr->phi());
861  }
862  if (debugL1Info_) {
863  edm::LogVerbatim("IsoTrack") << "tauJ p/pt " << itr->momentum() << " " << itr->pt() << " eta/phi "
864  << itr->eta() << " " << itr->phi();
865  }
866  }
867 
868  // L1 Central Jets
869  edm::Handle<l1extra::L1JetParticleCollection> l1CenJetHandle;
870  iEvent.getByToken(tok_L1extCenJet_, l1CenJetHandle);
871  for (itr = l1CenJetHandle->begin(), iL1Obj = 0; itr != l1CenJetHandle->end(); ++itr, iL1Obj++) {
872  if (iL1Obj < 1) {
873  t_L1CenJetPt->push_back(itr->pt());
874  t_L1CenJetEta->push_back(itr->eta());
875  t_L1CenJetPhi->push_back(itr->phi());
876  }
877  if (debugL1Info_) {
878  edm::LogVerbatim("IsoTrack") << "cenJ p/pt " << itr->momentum() << " " << itr->pt() << " eta/phi "
879  << itr->eta() << " " << itr->phi();
880  }
881  }
882 
883  // L1 Forward Jets
884  edm::Handle<l1extra::L1JetParticleCollection> l1FwdJetHandle;
885  iEvent.getByToken(tok_L1extFwdJet_, l1FwdJetHandle);
886  for (itr = l1FwdJetHandle->begin(), iL1Obj = 0; itr != l1FwdJetHandle->end(); ++itr, iL1Obj++) {
887  if (iL1Obj < 1) {
888  t_L1FwdJetPt->push_back(itr->pt());
889  t_L1FwdJetEta->push_back(itr->eta());
890  t_L1FwdJetPhi->push_back(itr->phi());
891  }
892  if (debugL1Info_) {
893  edm::LogVerbatim("IsoTrack") << "fwdJ p/pt " << itr->momentum() << " " << itr->pt() << " eta/phi "
894  << itr->eta() << " " << itr->phi();
895  }
896  }
897 
898  // L1 Isolated EM onjects
899  l1extra::L1EmParticleCollection::const_iterator itrEm;
900  edm::Handle<l1extra::L1EmParticleCollection> l1IsoEmHandle;
901  iEvent.getByToken(tok_L1extIsoEm_, l1IsoEmHandle);
902  for (itrEm = l1IsoEmHandle->begin(), iL1Obj = 0; itrEm != l1IsoEmHandle->end(); ++itrEm, iL1Obj++) {
903  if (iL1Obj < 1) {
904  t_L1IsoEMPt->push_back(itrEm->pt());
905  t_L1IsoEMEta->push_back(itrEm->eta());
906  t_L1IsoEMPhi->push_back(itrEm->phi());
907  }
908  if (debugL1Info_) {
909  edm::LogVerbatim("IsoTrack") << "isoEm p/pt " << itrEm->momentum() << " " << itrEm->pt() << " eta/phi "
910  << itrEm->eta() << " " << itrEm->phi();
911  }
912  }
913 
914  // L1 Non-Isolated EM onjects
915  edm::Handle<l1extra::L1EmParticleCollection> l1NonIsoEmHandle;
916  iEvent.getByToken(tok_L1extNoIsoEm_, l1NonIsoEmHandle);
917  for (itrEm = l1NonIsoEmHandle->begin(), iL1Obj = 0; itrEm != l1NonIsoEmHandle->end(); ++itrEm, iL1Obj++) {
918  if (iL1Obj < 1) {
919  t_L1NonIsoEMPt->push_back(itrEm->pt());
920  t_L1NonIsoEMEta->push_back(itrEm->eta());
921  t_L1NonIsoEMPhi->push_back(itrEm->phi());
922  }
923  if (debugL1Info_) {
924  edm::LogVerbatim("IsoTrack") << "nonIsoEm p/pt " << itrEm->momentum() << " " << itrEm->pt() << " eta/phi "
925  << itrEm->eta() << " " << itrEm->phi();
926  }
927  }
928 
929  // L1 Muons
930  l1extra::L1MuonParticleCollection::const_iterator itrMu;
931  edm::Handle<l1extra::L1MuonParticleCollection> l1MuHandle;
932  iEvent.getByToken(tok_L1extMu_, l1MuHandle);
933  for (itrMu = l1MuHandle->begin(), iL1Obj = 0; itrMu != l1MuHandle->end(); ++itrMu, iL1Obj++) {
934  if (iL1Obj < 1) {
935  t_L1MuonPt->push_back(itrMu->pt());
936  t_L1MuonEta->push_back(itrMu->eta());
937  t_L1MuonPhi->push_back(itrMu->phi());
938  }
939  if (debugL1Info_) {
940  edm::LogVerbatim("IsoTrack") << "l1muon p/pt " << itrMu->momentum() << " " << itrMu->pt() << " eta/phi "
941  << itrMu->eta() << " " << itrMu->phi();
942  }
943  }
944  }
945 
946  //============== store the information about all the Non-Fake vertices ===============
947 
948  edm::Handle<reco::VertexCollection> recVtxs;
949  iEvent.getByToken(tok_recVtx_, recVtxs);
950 
951  std::vector<reco::Track> svTracks;
952  math::XYZPoint leadPV(0, 0, 0);
953  double sumPtMax = -1.0;
954  for (unsigned int ind = 0; ind < recVtxs->size(); ind++) {
955  if (!((*recVtxs)[ind].isFake())) {
956  double vtxTrkSumPt = 0.0, vtxTrkSumPtWt = 0.0;
957  int vtxTrkNWt = 0;
958  double vtxTrkSumPtHP = 0.0, vtxTrkSumPtHPWt = 0.0;
959  int vtxTrkNHP = 0, vtxTrkNHPWt = 0;
960 
961  reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
962 
963  for (vtxTrack = (*recVtxs)[ind].tracks_begin(); vtxTrack != (*recVtxs)[ind].tracks_end(); vtxTrack++) {
964  if ((*vtxTrack)->pt() < pvTracksPtMin_)
965  continue;
966 
967  bool trkQuality = (*vtxTrack)->quality(trackQuality_);
968 
969  vtxTrkSumPt += (*vtxTrack)->pt();
970 
971  vtxTrkSumPt += (*vtxTrack)->pt();
972  if (trkQuality) {
973  vtxTrkSumPtHP += (*vtxTrack)->pt();
974  vtxTrkNHP++;
975  }
976 
977  double weight = (*recVtxs)[ind].trackWeight(*vtxTrack);
978  h_PVTracksWt->Fill(weight);
979  if (weight > 0.5) {
980  vtxTrkSumPtWt += (*vtxTrack)->pt();
981  vtxTrkNWt++;
982  if (trkQuality) {
983  vtxTrkSumPtHPWt += (*vtxTrack)->pt();
984  vtxTrkNHPWt++;
985  }
986  }
987  }
988 
989  if (vtxTrkSumPt > sumPtMax) {
990  sumPtMax = vtxTrkSumPt;
991  leadPV = math::XYZPoint((*recVtxs)[ind].x(), (*recVtxs)[ind].y(), (*recVtxs)[ind].z());
992  }
993 
994  t_PVx->push_back((*recVtxs)[ind].x());
995  t_PVy->push_back((*recVtxs)[ind].y());
996  t_PVz->push_back((*recVtxs)[ind].z());
997  t_PVisValid->push_back((*recVtxs)[ind].isValid());
998  t_PVNTracks->push_back((*recVtxs)[ind].tracksSize());
999  t_PVndof->push_back((*recVtxs)[ind].ndof());
1000  t_PVNTracksWt->push_back(vtxTrkNWt);
1001  t_PVTracksSumPt->push_back(vtxTrkSumPt);
1002  t_PVTracksSumPtWt->push_back(vtxTrkSumPtWt);
1003 
1004  t_PVNTracksHP->push_back(vtxTrkNHP);
1005  t_PVNTracksHPWt->push_back(vtxTrkNHPWt);
1006  t_PVTracksSumPtHP->push_back(vtxTrkSumPtHP);
1007  t_PVTracksSumPtHPWt->push_back(vtxTrkSumPtHPWt);
1008 
1009  if (myverbose_ == 4) {
1010  edm::LogVerbatim("IsoTrack") << "PV " << ind << " isValid " << (*recVtxs)[ind].isValid() << " isFake "
1011  << (*recVtxs)[ind].isFake() << " hasRefittedTracks() " << ind << " "
1012  << (*recVtxs)[ind].hasRefittedTracks() << " refittedTrksSize "
1013  << (*recVtxs)[ind].refittedTracks().size() << " tracksSize() "
1014  << (*recVtxs)[ind].tracksSize() << " sumPt " << vtxTrkSumPt;
1015  }
1016  } // if vtx is not Fake
1017  } // loop over PVs
1018  //===================================================================================
1019 
1020  // Get the beamspot
1021  edm::Handle<reco::BeamSpot> beamSpotH;
1022  iEvent.getByToken(tok_bs_, beamSpotH);
1023  math::XYZPoint bspot;
1024  bspot = (beamSpotH.isValid()) ? beamSpotH->position() : math::XYZPoint(0, 0, 0);
1025 
1026  //=====================================================================
1027 
1028  edm::Handle<reco::CaloJetCollection> jets;
1029  iEvent.getByToken(tok_jets_, jets);
1030  // edm::Handle<reco::JetExtendedAssociation::Container> jetExtender;
1031  // iEvent.getByLabel(JetExtender_,jetExtender);
1032 
1033  for (unsigned int ijet = 0; ijet < (*jets).size(); ijet++) {
1034  t_jetPt->push_back((*jets)[ijet].pt());
1035  t_jetEta->push_back((*jets)[ijet].eta());
1036  t_jetPhi->push_back((*jets)[ijet].phi());
1037  t_nTrksJetVtx->push_back(-1.0);
1038  t_nTrksJetCalo->push_back(-1.0);
1039  }
1040 
1041  //===================================================================================
1042 
1043  // get handles to calogeometry and calotopology
1044  const CaloGeometry *geo = &iSetup.getData(tok_geom_);
1045  const CaloTopology *caloTopology = &iSetup.getData(tok_caloTopology_);
1046  const HcalTopology *theHBHETopology = &iSetup.getData(tok_topo_);
1047 
1048  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
1049  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
1050  iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
1051  iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
1052 
1053  // Retrieve the good/bad ECAL channels from the DB
1054  const EcalChannelStatus *theEcalChStatus = &iSetup.getData(tok_ecalChStatus_);
1055  const EcalSeverityLevelAlgo *sevlv = &iSetup.getData(tok_sevlv_);
1056 
1057  // Retrieve trigger tower map
1058  const EcalTrigTowerConstituentsMap &ttMap = iSetup.getData(tok_htmap_);
1059 
1060  edm::Handle<HBHERecHitCollection> hbhe;
1061  iEvent.getByToken(tok_hbhe_, hbhe);
1062  if (!hbhe.isValid()) {
1063  ++nbad_;
1064  if (nbad_ < 10)
1065  edm::LogVerbatim("IsoTrack") << "No HBHE rechit collection";
1066  return;
1067  }
1068  const HBHERecHitCollection Hithbhe = *(hbhe.product());
1069 
1070  //get Handles to SimTracks and SimHits
1071  edm::Handle<edm::SimTrackContainer> SimTk;
1072  if (doMC_)
1073  iEvent.getByToken(tok_simTk_, SimTk);
1074 
1075  edm::Handle<edm::SimVertexContainer> SimVtx;
1076  if (doMC_)
1077  iEvent.getByToken(tok_simVtx_, SimVtx);
1078 
1079  //get Handles to PCaloHitContainers of eb/ee/hbhe
1080  edm::Handle<edm::PCaloHitContainer> pcaloeb;
1081  if (doMC_)
1082  iEvent.getByToken(tok_caloEB_, pcaloeb);
1083 
1084  edm::Handle<edm::PCaloHitContainer> pcaloee;
1085  if (doMC_)
1086  iEvent.getByToken(tok_caloEE_, pcaloee);
1087 
1088  edm::Handle<edm::PCaloHitContainer> pcalohh;
1089  if (doMC_)
1090  iEvent.getByToken(tok_caloHH_, pcalohh);
1091 
1092  //associates tracker rechits/simhits to a track
1093  std::unique_ptr<TrackerHitAssociator> associate;
1094  if (doMC_)
1095  associate = std::make_unique<TrackerHitAssociator>(iEvent, trackerHitAssociatorConfig_);
1096 
1097  //===================================================================================
1098 
1099  h_nTracks->Fill(trkCollection->size());
1100 
1101  int nTracks = 0;
1102 
1103  t_nTracks = trkCollection->size();
1104 
1105  // get the list of DetIds closest to the impact point of track on surface calorimeters
1106  std::vector<spr::propagatedTrackID> trkCaloDets;
1107  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality, trkCaloDets, false);
1108  std::vector<spr::propagatedTrackID>::const_iterator trkDetItr;
1109 
1110  if (myverbose_ > 2) {
1111  for (trkDetItr = trkCaloDets.begin(); trkDetItr != trkCaloDets.end(); trkDetItr++) {
1112  edm::LogVerbatim("IsoTrack") << trkDetItr->trkItr->p() << " " << trkDetItr->trkItr->eta() << " "
1113  << trkDetItr->okECAL << " " << trkDetItr->okHCAL;
1114  if (trkDetItr->okECAL) {
1115  if (trkDetItr->detIdECAL.subdetId() == EcalBarrel)
1116  edm::LogVerbatim("IsoTrack") << (EBDetId)trkDetItr->detIdECAL;
1117  else
1118  edm::LogVerbatim("IsoTrack") << (EEDetId)trkDetItr->detIdECAL;
1119  }
1120  if (trkDetItr->okHCAL)
1121  edm::LogVerbatim("IsoTrack") << (HcalDetId)trkDetItr->detIdHCAL;
1122  }
1123  }
1124 
1125  int nvtxTracks = 0;
1126  for (trkDetItr = trkCaloDets.begin(), nTracks = 0; trkDetItr != trkCaloDets.end(); trkDetItr++, nTracks++) {
1127  const reco::Track *pTrack = &(*(trkDetItr->trkItr));
1128 
1129  // find vertex index the track is associated with
1130  int pVtxTkId = -1;
1131  for (unsigned int ind = 0; ind < recVtxs->size(); ind++) {
1132  if (!((*recVtxs)[ind].isFake())) {
1133  reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
1134  for (vtxTrack = (*recVtxs)[ind].tracks_begin(); vtxTrack != (*recVtxs)[ind].tracks_end(); vtxTrack++) {
1135  const edm::RefToBase<reco::Track> pvtxTrack = (*vtxTrack);
1136  if (pTrack == pvtxTrack.get()) {
1137  pVtxTkId = ind;
1138  break;
1139  if (myverbose_ > 2) {
1140  if (pTrack->pt() > 1.0) {
1141  edm::LogVerbatim("IsoTrack") << "Debug the track association with vertex ";
1142  edm::LogVerbatim("IsoTrack") << pTrack << " " << pvtxTrack.get();
1143  edm::LogVerbatim("IsoTrack")
1144  << " trkVtxIndex " << nvtxTracks << " vtx " << ind << " pt " << pTrack->pt() << " eta "
1145  << pTrack->eta() << " " << pTrack->pt() - pvtxTrack->pt() << " "
1146  << pTrack->eta() - pvtxTrack->eta();
1147  nvtxTracks++;
1148  }
1149  }
1150  }
1151  }
1152  }
1153  }
1154 
1155  const reco::HitPattern &hitp = pTrack->hitPattern();
1156  int nLayersCrossed = hitp.trackerLayersWithMeasurement();
1157  int nOuterHits = hitp.stripTOBLayersWithMeasurement() + hitp.stripTECLayersWithMeasurement();
1158 
1159  bool ifGood = pTrack->quality(trackQuality_);
1160  double pt1 = pTrack->pt();
1161  double p1 = pTrack->p();
1162  double eta1 = pTrack->momentum().eta();
1163  double phi1 = pTrack->momentum().phi();
1164  double etaEcal1 = trkDetItr->etaECAL;
1165  double phiEcal1 = trkDetItr->phiECAL;
1166  double etaHcal1 = trkDetItr->etaHCAL;
1167  double phiHcal1 = trkDetItr->phiHCAL;
1168  double dxy1 = pTrack->dxy();
1169  double dz1 = pTrack->dz();
1170  double chisq1 = pTrack->normalizedChi2();
1171  double dxybs1 = beamSpotH.isValid() ? pTrack->dxy(bspot) : pTrack->dxy();
1172  double dzbs1 = beamSpotH.isValid() ? pTrack->dz(bspot) : pTrack->dz();
1173  double dxypv1 = pTrack->dxy();
1174  double dzpv1 = pTrack->dz();
1175  if (pVtxTkId >= 0) {
1176  math::XYZPoint thisTkPV =
1177  math::XYZPoint((*recVtxs)[pVtxTkId].x(), (*recVtxs)[pVtxTkId].y(), (*recVtxs)[pVtxTkId].z());
1178  dxypv1 = pTrack->dxy(thisTkPV);
1179  dzpv1 = pTrack->dz(thisTkPV);
1180  }
1181 
1182  h_recEtaPt_0->Fill(eta1, pt1);
1183  h_recEtaP_0->Fill(eta1, p1);
1184  h_recPt_0->Fill(pt1);
1185  h_recP_0->Fill(p1);
1186  h_recEta_0->Fill(eta1);
1187  h_recPhi_0->Fill(phi1);
1188 
1189  if (ifGood && nLayersCrossed > 7) {
1190  h_recEtaPt_1->Fill(eta1, pt1);
1191  h_recEtaP_1->Fill(eta1, p1);
1192  h_recPt_1->Fill(pt1);
1193  h_recP_1->Fill(p1);
1194  h_recEta_1->Fill(eta1);
1195  h_recPhi_1->Fill(phi1);
1196  }
1197 
1198  if (!ifGood)
1199  continue;
1200 
1201  if (writeAllTracks_ && p1 > 2.0 && nLayersCrossed > 7) {
1202  t_trackPAll->push_back(p1);
1203  t_trackEtaAll->push_back(eta1);
1204  t_trackPhiAll->push_back(phi1);
1205  t_trackPtAll->push_back(pt1);
1206  t_trackDxyAll->push_back(dxy1);
1207  t_trackDzAll->push_back(dz1);
1208  t_trackDxyPVAll->push_back(dxypv1);
1209  t_trackDzPVAll->push_back(dzpv1);
1210  t_trackChiSqAll->push_back(chisq1);
1211  }
1212  if (doMC_) {
1213  edm::SimTrackContainer::const_iterator matchedSimTrkAll =
1214  spr::matchedSimTrack(iEvent, SimTk, SimVtx, pTrack, *associate, false);
1215  if (writeAllTracks_ && matchedSimTrkAll != SimTk->end())
1216  t_trackPdgIdAll->push_back(matchedSimTrkAll->type());
1217  }
1218 
1219  if (pt1 > minTrackP_ && std::abs(eta1) < maxTrackEta_ && trkDetItr->okECAL) {
1220  double maxNearP31x31 = 999.0, maxNearP25x25 = 999.0, maxNearP21x21 = 999.0, maxNearP15x15 = 999.0;
1221  maxNearP31x31 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 15, 15);
1222  maxNearP25x25 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 12, 12);
1223  maxNearP21x21 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 10, 10);
1224  maxNearP15x15 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 7, 7);
1225 
1226  int iTrkEtaBin = -1, iTrkMomBin = -1;
1227  for (unsigned int ieta = 0; ieta < NEtaBins; ieta++) {
1228  if (std::abs(eta1) > genPartEtaBins[ieta] && std::abs(eta1) < genPartEtaBins[ieta + 1])
1229  iTrkEtaBin = ieta;
1230  }
1231  for (unsigned int ipt = 0; ipt < NPBins; ipt++) {
1232  if (p1 > genPartPBins[ipt] && p1 < genPartPBins[ipt + 1])
1233  iTrkMomBin = ipt;
1234  }
1235  if (iTrkMomBin >= 0 && iTrkEtaBin >= 0) {
1236  h_maxNearP31x31[iTrkMomBin][iTrkEtaBin]->Fill(maxNearP31x31);
1237  h_maxNearP25x25[iTrkMomBin][iTrkEtaBin]->Fill(maxNearP25x25);
1238  h_maxNearP21x21[iTrkMomBin][iTrkEtaBin]->Fill(maxNearP21x21);
1239  h_maxNearP15x15[iTrkMomBin][iTrkEtaBin]->Fill(maxNearP15x15);
1240  }
1241  if (maxNearP31x31 < 0.0 && nLayersCrossed > 7 && nOuterHits > 4) {
1242  h_recEtaPt_2->Fill(eta1, pt1);
1243  h_recEtaP_2->Fill(eta1, p1);
1244  h_recPt_2->Fill(pt1);
1245  h_recP_2->Fill(p1);
1246  h_recEta_2->Fill(eta1);
1247  h_recPhi_2->Fill(phi1);
1248  }
1249 
1250  // if charge isolated in ECAL, store the further quantities
1251  if (maxNearP31x31 < 1.0) {
1252  haveIsoTrack = true;
1253 
1254  // get the matching simTrack
1255  double simTrackP = -1;
1256  if (doMC_) {
1257  edm::SimTrackContainer::const_iterator matchedSimTrk =
1258  spr::matchedSimTrack(iEvent, SimTk, SimVtx, pTrack, *associate, false);
1259  if (matchedSimTrk != SimTk->end())
1260  simTrackP = matchedSimTrk->momentum().P();
1261  }
1262 
1263  // get ECal Tranverse Profile
1264  std::pair<double, bool> e7x7P, e9x9P, e11x11P, e15x15P;
1265  std::pair<double, bool> e7x7_10SigP, e9x9_10SigP, e11x11_10SigP, e15x15_10SigP;
1266  std::pair<double, bool> e7x7_15SigP, e9x9_15SigP, e11x11_15SigP, e15x15_15SigP;
1267  std::pair<double, bool> e7x7_20SigP, e9x9_20SigP, e11x11_20SigP, e15x15_20SigP;
1268  std::pair<double, bool> e7x7_25SigP, e9x9_25SigP, e11x11_25SigP, e15x15_25SigP;
1269  std::pair<double, bool> e7x7_30SigP, e9x9_30SigP, e11x11_30SigP, e15x15_30SigP;
1270 
1271  spr::caloSimInfo simInfo3x3, simInfo5x5, simInfo7x7, simInfo9x9;
1272  spr::caloSimInfo simInfo11x11, simInfo13x13, simInfo15x15, simInfo21x21, simInfo25x25, simInfo31x31;
1273  double trkEcalEne = 0;
1274 
1275  const DetId isoCell = trkDetItr->detIdECAL;
1276  e7x7P = spr::eECALmatrix(isoCell,
1277  barrelRecHitsHandle,
1278  endcapRecHitsHandle,
1279  *theEcalChStatus,
1280  geo,
1281  caloTopology,
1282  sevlv,
1283  3,
1284  3,
1285  -100.0,
1286  -100.0,
1287  tMinE_,
1288  tMaxE_);
1289  e9x9P = spr::eECALmatrix(isoCell,
1290  barrelRecHitsHandle,
1291  endcapRecHitsHandle,
1292  *theEcalChStatus,
1293  geo,
1294  caloTopology,
1295  sevlv,
1296  4,
1297  4,
1298  -100.0,
1299  -100.0,
1300  tMinE_,
1301  tMaxE_);
1302  e11x11P = spr::eECALmatrix(isoCell,
1303  barrelRecHitsHandle,
1304  endcapRecHitsHandle,
1305  *theEcalChStatus,
1306  geo,
1307  caloTopology,
1308  sevlv,
1309  5,
1310  5,
1311  -100.0,
1312  -100.0,
1313  tMinE_,
1314  tMaxE_);
1315  e15x15P = spr::eECALmatrix(isoCell,
1316  barrelRecHitsHandle,
1317  endcapRecHitsHandle,
1318  *theEcalChStatus,
1319  geo,
1320  caloTopology,
1321  sevlv,
1322  7,
1323  7,
1324  -100.0,
1325  -100.0,
1326  tMinE_,
1327  tMaxE_);
1328 
1329  e7x7_10SigP = spr::eECALmatrix(isoCell,
1330  barrelRecHitsHandle,
1331  endcapRecHitsHandle,
1332  *theEcalChStatus,
1333  geo,
1334  caloTopology,
1335  sevlv,
1336  3,
1337  3,
1338  0.030,
1339  0.150,
1340  tMinE_,
1341  tMaxE_);
1342  e9x9_10SigP = spr::eECALmatrix(isoCell,
1343  barrelRecHitsHandle,
1344  endcapRecHitsHandle,
1345  *theEcalChStatus,
1346  geo,
1347  caloTopology,
1348  sevlv,
1349  4,
1350  4,
1351  0.030,
1352  0.150,
1353  tMinE_,
1354  tMaxE_);
1355  e11x11_10SigP = spr::eECALmatrix(isoCell,
1356  barrelRecHitsHandle,
1357  endcapRecHitsHandle,
1358  *theEcalChStatus,
1359  geo,
1360  caloTopology,
1361  sevlv,
1362  5,
1363  5,
1364  0.030,
1365  0.150,
1366  tMinE_,
1367  tMaxE_);
1368  e15x15_10SigP = spr::eECALmatrix(isoCell,
1369  barrelRecHitsHandle,
1370  endcapRecHitsHandle,
1371  *theEcalChStatus,
1372  geo,
1373  caloTopology,
1374  sevlv,
1375  7,
1376  7,
1377  0.030,
1378  0.150,
1379  tMinE_,
1380  tMaxE_);
1381 
1382  e7x7_15SigP = spr::eECALmatrix(isoCell,
1383  barrelRecHitsHandle,
1384  endcapRecHitsHandle,
1385  *theEcalChStatus,
1386  geo,
1387  caloTopology,
1388  sevlv,
1389  ttMap,
1390  3,
1391  3,
1392  0.20,
1393  0.45,
1394  tMinE_,
1395  tMaxE_);
1396  e9x9_15SigP = spr::eECALmatrix(isoCell,
1397  barrelRecHitsHandle,
1398  endcapRecHitsHandle,
1399  *theEcalChStatus,
1400  geo,
1401  caloTopology,
1402  sevlv,
1403  ttMap,
1404  4,
1405  4,
1406  0.20,
1407  0.45,
1408  tMinE_,
1409  tMaxE_);
1410  e11x11_15SigP = spr::eECALmatrix(isoCell,
1411  barrelRecHitsHandle,
1412  endcapRecHitsHandle,
1413  *theEcalChStatus,
1414  geo,
1415  caloTopology,
1416  sevlv,
1417  ttMap,
1418  5,
1419  5,
1420  0.20,
1421  0.45,
1422  tMinE_,
1423  tMaxE_);
1424  e15x15_15SigP = spr::eECALmatrix(isoCell,
1425  barrelRecHitsHandle,
1426  endcapRecHitsHandle,
1427  *theEcalChStatus,
1428  geo,
1429  caloTopology,
1430  sevlv,
1431  ttMap,
1432  7,
1433  7,
1434  0.20,
1435  0.45,
1436  tMinE_,
1437  tMaxE_,
1438  false);
1439 
1440  e7x7_20SigP = spr::eECALmatrix(isoCell,
1441  barrelRecHitsHandle,
1442  endcapRecHitsHandle,
1443  *theEcalChStatus,
1444  geo,
1445  caloTopology,
1446  sevlv,
1447  3,
1448  3,
1449  0.060,
1450  0.300,
1451  tMinE_,
1452  tMaxE_);
1453  e9x9_20SigP = spr::eECALmatrix(isoCell,
1454  barrelRecHitsHandle,
1455  endcapRecHitsHandle,
1456  *theEcalChStatus,
1457  geo,
1458  caloTopology,
1459  sevlv,
1460  4,
1461  4,
1462  0.060,
1463  0.300,
1464  tMinE_,
1465  tMaxE_);
1466  e11x11_20SigP = spr::eECALmatrix(isoCell,
1467  barrelRecHitsHandle,
1468  endcapRecHitsHandle,
1469  *theEcalChStatus,
1470  geo,
1471  caloTopology,
1472  sevlv,
1473  5,
1474  5,
1475  0.060,
1476  0.300,
1477  tMinE_,
1478  tMaxE_);
1479  e15x15_20SigP = spr::eECALmatrix(isoCell,
1480  barrelRecHitsHandle,
1481  endcapRecHitsHandle,
1482  *theEcalChStatus,
1483  geo,
1484  caloTopology,
1485  sevlv,
1486  7,
1487  7,
1488  0.060,
1489  0.300,
1490  tMinE_,
1491  tMaxE_);
1492 
1493  e7x7_25SigP = spr::eECALmatrix(isoCell,
1494  barrelRecHitsHandle,
1495  endcapRecHitsHandle,
1496  *theEcalChStatus,
1497  geo,
1498  caloTopology,
1499  sevlv,
1500  3,
1501  3,
1502  0.075,
1503  0.375,
1504  tMinE_,
1505  tMaxE_);
1506  e9x9_25SigP = spr::eECALmatrix(isoCell,
1507  barrelRecHitsHandle,
1508  endcapRecHitsHandle,
1509  *theEcalChStatus,
1510  geo,
1511  caloTopology,
1512  sevlv,
1513  4,
1514  4,
1515  0.075,
1516  0.375,
1517  tMinE_,
1518  tMaxE_);
1519  e11x11_25SigP = spr::eECALmatrix(isoCell,
1520  barrelRecHitsHandle,
1521  endcapRecHitsHandle,
1522  *theEcalChStatus,
1523  geo,
1524  caloTopology,
1525  sevlv,
1526  5,
1527  5,
1528  0.075,
1529  0.375,
1530  tMinE_,
1531  tMaxE_);
1532  e15x15_25SigP = spr::eECALmatrix(isoCell,
1533  barrelRecHitsHandle,
1534  endcapRecHitsHandle,
1535  *theEcalChStatus,
1536  geo,
1537  caloTopology,
1538  sevlv,
1539  7,
1540  7,
1541  0.075,
1542  0.375,
1543  tMinE_,
1544  tMaxE_);
1545 
1546  e7x7_30SigP = spr::eECALmatrix(isoCell,
1547  barrelRecHitsHandle,
1548  endcapRecHitsHandle,
1549  *theEcalChStatus,
1550  geo,
1551  caloTopology,
1552  sevlv,
1553  3,
1554  3,
1555  0.090,
1556  0.450,
1557  tMinE_,
1558  tMaxE_);
1559  e9x9_30SigP = spr::eECALmatrix(isoCell,
1560  barrelRecHitsHandle,
1561  endcapRecHitsHandle,
1562  *theEcalChStatus,
1563  geo,
1564  caloTopology,
1565  sevlv,
1566  4,
1567  4,
1568  0.090,
1569  0.450,
1570  tMinE_,
1571  tMaxE_);
1572  e11x11_30SigP = spr::eECALmatrix(isoCell,
1573  barrelRecHitsHandle,
1574  endcapRecHitsHandle,
1575  *theEcalChStatus,
1576  geo,
1577  caloTopology,
1578  sevlv,
1579  5,
1580  5,
1581  0.090,
1582  0.450,
1583  tMinE_,
1584  tMaxE_);
1585  e15x15_30SigP = spr::eECALmatrix(isoCell,
1586  barrelRecHitsHandle,
1587  endcapRecHitsHandle,
1588  *theEcalChStatus,
1589  geo,
1590  caloTopology,
1591  sevlv,
1592  7,
1593  7,
1594  0.090,
1595  0.450,
1596  tMinE_,
1597  tMaxE_);
1598  if (myverbose_ == 2) {
1599  edm::LogVerbatim("IsoTrack") << "clean ecal rechit ";
1600  edm::LogVerbatim("IsoTrack") << "e7x7 " << e7x7P.first << " e9x9 " << e9x9P.first << " e11x11 "
1601  << e11x11P.first << " e15x15 " << e15x15P.first;
1602  edm::LogVerbatim("IsoTrack") << "e7x7_10Sig " << e7x7_10SigP.first << " e11x11_10Sig " << e11x11_10SigP.first
1603  << " e15x15_10Sig " << e15x15_10SigP.first;
1604  }
1605 
1606  if (doMC_) {
1607  // check the energy from SimHits
1608  spr::eECALSimInfo(
1609  iEvent, isoCell, geo, caloTopology, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate, 1, 1, simInfo3x3);
1610  spr::eECALSimInfo(
1611  iEvent, isoCell, geo, caloTopology, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate, 2, 2, simInfo5x5);
1612  spr::eECALSimInfo(
1613  iEvent, isoCell, geo, caloTopology, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate, 3, 3, simInfo7x7);
1614  spr::eECALSimInfo(
1615  iEvent, isoCell, geo, caloTopology, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate, 4, 4, simInfo9x9);
1616  spr::eECALSimInfo(iEvent,
1617  isoCell,
1618  geo,
1619  caloTopology,
1620  pcaloeb,
1621  pcaloee,
1622  SimTk,
1623  SimVtx,
1624  pTrack,
1625  *associate,
1626  5,
1627  5,
1628  simInfo11x11);
1629  spr::eECALSimInfo(iEvent,
1630  isoCell,
1631  geo,
1632  caloTopology,
1633  pcaloeb,
1634  pcaloee,
1635  SimTk,
1636  SimVtx,
1637  pTrack,
1638  *associate,
1639  6,
1640  6,
1641  simInfo13x13);
1642  spr::eECALSimInfo(iEvent,
1643  isoCell,
1644  geo,
1645  caloTopology,
1646  pcaloeb,
1647  pcaloee,
1648  SimTk,
1649  SimVtx,
1650  pTrack,
1651  *associate,
1652  7,
1653  7,
1654  simInfo15x15,
1655  150.0,
1656  false);
1657  spr::eECALSimInfo(iEvent,
1658  isoCell,
1659  geo,
1660  caloTopology,
1661  pcaloeb,
1662  pcaloee,
1663  SimTk,
1664  SimVtx,
1665  pTrack,
1666  *associate,
1667  10,
1668  10,
1669  simInfo21x21);
1670  spr::eECALSimInfo(iEvent,
1671  isoCell,
1672  geo,
1673  caloTopology,
1674  pcaloeb,
1675  pcaloee,
1676  SimTk,
1677  SimVtx,
1678  pTrack,
1679  *associate,
1680  12,
1681  12,
1682  simInfo25x25);
1683  spr::eECALSimInfo(iEvent,
1684  isoCell,
1685  geo,
1686  caloTopology,
1687  pcaloeb,
1688  pcaloee,
1689  SimTk,
1690  SimVtx,
1691  pTrack,
1692  *associate,
1693  15,
1694  15,
1695  simInfo31x31);
1696 
1697  trkEcalEne = spr::eCaloSimInfo(iEvent, geo, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate);
1698  if (myverbose_ == 1) {
1699  edm::LogVerbatim("IsoTrack") << "Track momentum " << pt1;
1700  edm::LogVerbatim("IsoTrack") << "ecal siminfo ";
1701  edm::LogVerbatim("IsoTrack") << "simInfo3x3: eTotal " << simInfo3x3.eTotal << " eMatched "
1702  << simInfo3x3.eMatched << " eRest " << simInfo3x3.eRest << " eGamma "
1703  << simInfo3x3.eGamma << " eNeutralHad " << simInfo3x3.eNeutralHad
1704  << " eChargedHad " << simInfo3x3.eChargedHad;
1705  edm::LogVerbatim("IsoTrack") << "simInfo5x5: eTotal " << simInfo5x5.eTotal << " eMatched "
1706  << simInfo5x5.eMatched << " eRest " << simInfo5x5.eRest << " eGamma "
1707  << simInfo5x5.eGamma << " eNeutralHad " << simInfo5x5.eNeutralHad
1708  << " eChargedHad " << simInfo5x5.eChargedHad;
1709  edm::LogVerbatim("IsoTrack") << "simInfo7x7: eTotal " << simInfo7x7.eTotal << " eMatched "
1710  << simInfo7x7.eMatched << " eRest " << simInfo7x7.eRest << " eGamma "
1711  << simInfo7x7.eGamma << " eNeutralHad " << simInfo7x7.eNeutralHad
1712  << " eChargedHad " << simInfo7x7.eChargedHad;
1713  edm::LogVerbatim("IsoTrack") << "simInfo9x9: eTotal " << simInfo9x9.eTotal << " eMatched "
1714  << simInfo9x9.eMatched << " eRest " << simInfo9x9.eRest << " eGamma "
1715  << simInfo9x9.eGamma << " eNeutralHad " << simInfo9x9.eNeutralHad
1716  << " eChargedHad " << simInfo9x9.eChargedHad;
1717  edm::LogVerbatim("IsoTrack") << "simInfo11x11: eTotal " << simInfo11x11.eTotal << " eMatched "
1718  << simInfo11x11.eMatched << " eRest " << simInfo11x11.eRest << " eGamma "
1719  << simInfo11x11.eGamma << " eNeutralHad " << simInfo11x11.eNeutralHad
1720  << " eChargedHad " << simInfo11x11.eChargedHad;
1721  edm::LogVerbatim("IsoTrack") << "simInfo15x15: eTotal " << simInfo15x15.eTotal << " eMatched "
1722  << simInfo15x15.eMatched << " eRest " << simInfo15x15.eRest << " eGamma "
1723  << simInfo15x15.eGamma << " eNeutralHad " << simInfo15x15.eNeutralHad
1724  << " eChargedHad " << simInfo15x15.eChargedHad;
1725  edm::LogVerbatim("IsoTrack") << "simInfo31x31: eTotal " << simInfo31x31.eTotal << " eMatched "
1726  << simInfo31x31.eMatched << " eRest " << simInfo31x31.eRest << " eGamma "
1727  << simInfo31x31.eGamma << " eNeutralHad " << simInfo31x31.eNeutralHad
1728  << " eChargedHad " << simInfo31x31.eChargedHad;
1729  edm::LogVerbatim("IsoTrack") << "trkEcalEne" << trkEcalEne;
1730  }
1731  }
1732 
1733  // ======= Get HCAL information
1734  double hcalScale = 1.0;
1735  if (std::abs(pTrack->eta()) < 1.4) {
1736  hcalScale = 120.0;
1737  } else {
1738  hcalScale = 135.0;
1739  }
1740 
1741  double maxNearHcalP3x3 = -1, maxNearHcalP5x5 = -1, maxNearHcalP7x7 = -1;
1742  maxNearHcalP3x3 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 1, 1);
1743  maxNearHcalP5x5 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 2, 2);
1744  maxNearHcalP7x7 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 3, 3);
1745 
1746  double h3x3 = 0, h5x5 = 0, h7x7 = 0;
1747  double h3x3Sig = 0, h5x5Sig = 0, h7x7Sig = 0;
1748  double trkHcalEne = 0;
1749  spr::caloSimInfo hsimInfo3x3, hsimInfo5x5, hsimInfo7x7;
1750 
1751  if (trkDetItr->okHCAL) {
1752  const DetId ClosestCell(trkDetItr->detIdHCAL);
1753  // bool includeHO=false, bool algoNew=true, bool debug=false
1754  h3x3 = spr::eHCALmatrix(
1755  theHBHETopology, ClosestCell, hbhe, 1, 1, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_, tMaxH_);
1756  h5x5 = spr::eHCALmatrix(
1757  theHBHETopology, ClosestCell, hbhe, 2, 2, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_, tMaxH_);
1758  h7x7 = spr::eHCALmatrix(
1759  theHBHETopology, ClosestCell, hbhe, 3, 3, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_, tMaxH_);
1760  h3x3Sig = spr::eHCALmatrix(
1761  theHBHETopology, ClosestCell, hbhe, 1, 1, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_, tMaxH_);
1762  h5x5Sig = spr::eHCALmatrix(
1763  theHBHETopology, ClosestCell, hbhe, 2, 2, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_, tMaxH_);
1764  h7x7Sig = spr::eHCALmatrix(
1765  theHBHETopology, ClosestCell, hbhe, 3, 3, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_, tMaxH_);
1766 
1767  if (myverbose_ == 2) {
1768  edm::LogVerbatim("IsoTrack") << "HCAL 3x3 " << h3x3 << " " << h3x3Sig << " 5x5 " << h5x5 << " " << h5x5Sig
1769  << " 7x7 " << h7x7 << " " << h7x7Sig;
1770  }
1771 
1772  if (doMC_) {
1773  spr::eHCALSimInfo(
1774  iEvent, theHBHETopology, ClosestCell, geo, pcalohh, SimTk, SimVtx, pTrack, *associate, 1, 1, hsimInfo3x3);
1775  spr::eHCALSimInfo(
1776  iEvent, theHBHETopology, ClosestCell, geo, pcalohh, SimTk, SimVtx, pTrack, *associate, 2, 2, hsimInfo5x5);
1777  spr::eHCALSimInfo(iEvent,
1778  theHBHETopology,
1779  ClosestCell,
1780  geo,
1781  pcalohh,
1782  SimTk,
1783  SimVtx,
1784  pTrack,
1785  *associate,
1786  3,
1787  3,
1788  hsimInfo7x7,
1789  150.0,
1790  false,
1791  false);
1792  trkHcalEne = spr::eCaloSimInfo(iEvent, geo, pcalohh, SimTk, SimVtx, pTrack, *associate);
1793  if (myverbose_ == 1) {
1794  edm::LogVerbatim("IsoTrack") << "Hcal siminfo ";
1795  edm::LogVerbatim("IsoTrack")
1796  << "hsimInfo3x3: eTotal " << hsimInfo3x3.eTotal << " eMatched " << hsimInfo3x3.eMatched << " eRest "
1797  << hsimInfo3x3.eRest << " eGamma " << hsimInfo3x3.eGamma << " eNeutralHad " << hsimInfo3x3.eNeutralHad
1798  << " eChargedHad " << hsimInfo3x3.eChargedHad;
1799  edm::LogVerbatim("IsoTrack")
1800  << "hsimInfo5x5: eTotal " << hsimInfo5x5.eTotal << " eMatched " << hsimInfo5x5.eMatched << " eRest "
1801  << hsimInfo5x5.eRest << " eGamma " << hsimInfo5x5.eGamma << " eNeutralHad " << hsimInfo5x5.eNeutralHad
1802  << " eChargedHad " << hsimInfo5x5.eChargedHad;
1803  edm::LogVerbatim("IsoTrack")
1804  << "hsimInfo7x7: eTotal " << hsimInfo7x7.eTotal << " eMatched " << hsimInfo7x7.eMatched << " eRest "
1805  << hsimInfo7x7.eRest << " eGamma " << hsimInfo7x7.eGamma << " eNeutralHad " << hsimInfo7x7.eNeutralHad
1806  << " eChargedHad " << hsimInfo7x7.eChargedHad;
1807  edm::LogVerbatim("IsoTrack") << "trkHcalEne " << trkHcalEne;
1808  }
1809  }
1810 
1811  // debug the ecal and hcal matrix
1812  if (myverbose_ == 4) {
1813  edm::LogVerbatim("IsoTrack") << "Run " << iEvent.id().run() << " Event " << iEvent.id().event();
1814  std::vector<std::pair<DetId, double> > v7x7 =
1815  spr::eHCALmatrixCell(theHBHETopology, ClosestCell, hbhe, 3, 3, false, false);
1816  double sumv = 0.0;
1817 
1818  for (unsigned int iv = 0; iv < v7x7.size(); iv++) {
1819  sumv += v7x7[iv].second;
1820  }
1821  edm::LogVerbatim("IsoTrack") << "h7x7 " << h7x7 << " v7x7 " << sumv << " in " << v7x7.size();
1822  for (unsigned int iv = 0; iv < v7x7.size(); iv++) {
1823  HcalDetId id = v7x7[iv].first;
1824  edm::LogVerbatim("IsoTrack") << " Cell " << iv << " 0x" << std::hex << v7x7[iv].first() << std::dec << " "
1825  << id << " Energy " << v7x7[iv].second;
1826  }
1827  }
1828  }
1829  if (doMC_) {
1830  trkHcalEne = spr::eCaloSimInfo(iEvent, geo, pcalohh, SimTk, SimVtx, pTrack, *associate);
1831  }
1832 
1833  // ====================================================================================================
1834  // get diff between track outermost hit position and the propagation point at outermost surface of tracker
1835  std::pair<math::XYZPoint, double> point2_TK0 = spr::propagateTrackerEnd(pTrack, bField, false);
1836  math::XYZPoint diff(pTrack->outerPosition().X() - point2_TK0.first.X(),
1837  pTrack->outerPosition().Y() - point2_TK0.first.Y(),
1838  pTrack->outerPosition().Z() - point2_TK0.first.Z());
1839  double trackOutPosOutHitDr = diff.R();
1840  double trackL = point2_TK0.second;
1841  if (myverbose_ == 5) {
1842  edm::LogVerbatim("IsoTrack") << " propagted " << point2_TK0.first << " " << point2_TK0.first.eta() << " "
1843  << point2_TK0.first.phi();
1844  edm::LogVerbatim("IsoTrack") << " outerPosition() " << pTrack->outerPosition() << " "
1845  << pTrack->outerPosition().eta() << " " << pTrack->outerPosition().phi();
1846  edm::LogVerbatim("IsoTrack") << "diff " << diff << " diffR " << diff.R() << " diffR/L "
1847  << diff.R() / point2_TK0.second;
1848  }
1849 
1850  for (unsigned int ind = 0; ind < recVtxs->size(); ind++) {
1851  if (!((*recVtxs)[ind].isFake())) {
1852  reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
1853  if (reco::deltaR(eta1, phi1, (*vtxTrack)->eta(), (*vtxTrack)->phi()) < 0.01)
1854  t_trackPVIdx->push_back(ind);
1855  else
1856  t_trackPVIdx->push_back(-1);
1857  }
1858  }
1859 
1860  // Fill the tree Branches here
1861  t_trackPVIdx->push_back(pVtxTkId);
1862  t_trackP->push_back(p1);
1863  t_trackPt->push_back(pt1);
1864  t_trackEta->push_back(eta1);
1865  t_trackPhi->push_back(phi1);
1866  t_trackEcalEta->push_back(etaEcal1);
1867  t_trackEcalPhi->push_back(phiEcal1);
1868  t_trackHcalEta->push_back(etaHcal1);
1869  t_trackHcalPhi->push_back(phiHcal1);
1870  t_trackDxy->push_back(dxy1);
1871  t_trackDz->push_back(dz1);
1872  t_trackDxyBS->push_back(dxybs1);
1873  t_trackDzBS->push_back(dzbs1);
1874  t_trackDxyPV->push_back(dxypv1);
1875  t_trackDzPV->push_back(dzpv1);
1876  t_trackChiSq->push_back(chisq1);
1877  t_trackNOuterHits->push_back(nOuterHits);
1878  t_NLayersCrossed->push_back(nLayersCrossed);
1879 
1880  t_trackHitsTOB->push_back(hitp.stripTOBLayersWithMeasurement());
1881  t_trackHitsTEC->push_back(hitp.stripTECLayersWithMeasurement());
1882  t_trackHitInMissTOB->push_back(hitp.stripTOBLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
1883  t_trackHitInMissTEC->push_back(hitp.stripTECLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
1884  t_trackHitInMissTIB->push_back(hitp.stripTIBLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
1885  t_trackHitInMissTID->push_back(hitp.stripTIDLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
1886  t_trackHitInMissTIBTID->push_back(hitp.stripTIBLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS) +
1887  hitp.stripTIDLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
1888 
1889  t_trackHitOutMissTOB->push_back(hitp.stripTOBLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
1890  t_trackHitOutMissTEC->push_back(hitp.stripTECLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
1891  t_trackHitOutMissTIB->push_back(hitp.stripTIBLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
1892  t_trackHitOutMissTID->push_back(hitp.stripTIDLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
1893  t_trackHitOutMissTOBTEC->push_back(hitp.stripTOBLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS) +
1894  hitp.stripTECLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
1895 
1896  t_trackHitInMeasTOB->push_back(hitp.stripTOBLayersWithMeasurement());
1897  t_trackHitInMeasTEC->push_back(hitp.stripTECLayersWithMeasurement());
1898  t_trackHitInMeasTIB->push_back(hitp.stripTIBLayersWithMeasurement());
1899  t_trackHitInMeasTID->push_back(hitp.stripTIDLayersWithMeasurement());
1900  t_trackHitOutMeasTOB->push_back(hitp.stripTOBLayersWithMeasurement());
1901  t_trackHitOutMeasTEC->push_back(hitp.stripTECLayersWithMeasurement());
1902  t_trackHitOutMeasTIB->push_back(hitp.stripTIBLayersWithMeasurement());
1903  t_trackHitOutMeasTID->push_back(hitp.stripTIDLayersWithMeasurement());
1904  t_trackOutPosOutHitDr->push_back(trackOutPosOutHitDr);
1905  t_trackL->push_back(trackL);
1906 
1907  t_maxNearP31x31->push_back(maxNearP31x31);
1908  t_maxNearP21x21->push_back(maxNearP21x21);
1909 
1910  t_ecalSpike11x11->push_back(e11x11P.second);
1911  t_e7x7->push_back(e7x7P.first);
1912  t_e9x9->push_back(e9x9P.first);
1913  t_e11x11->push_back(e11x11P.first);
1914  t_e15x15->push_back(e15x15P.first);
1915 
1916  t_e7x7_10Sig->push_back(e7x7_10SigP.first);
1917  t_e9x9_10Sig->push_back(e9x9_10SigP.first);
1918  t_e11x11_10Sig->push_back(e11x11_10SigP.first);
1919  t_e15x15_10Sig->push_back(e15x15_10SigP.first);
1920  t_e7x7_15Sig->push_back(e7x7_15SigP.first);
1921  t_e9x9_15Sig->push_back(e9x9_15SigP.first);
1922  t_e11x11_15Sig->push_back(e11x11_15SigP.first);
1923  t_e15x15_15Sig->push_back(e15x15_15SigP.first);
1924  t_e7x7_20Sig->push_back(e7x7_20SigP.first);
1925  t_e9x9_20Sig->push_back(e9x9_20SigP.first);
1926  t_e11x11_20Sig->push_back(e11x11_20SigP.first);
1927  t_e15x15_20Sig->push_back(e15x15_20SigP.first);
1928  t_e7x7_25Sig->push_back(e7x7_25SigP.first);
1929  t_e9x9_25Sig->push_back(e9x9_25SigP.first);
1930  t_e11x11_25Sig->push_back(e11x11_25SigP.first);
1931  t_e15x15_25Sig->push_back(e15x15_25SigP.first);
1932  t_e7x7_30Sig->push_back(e7x7_30SigP.first);
1933  t_e9x9_30Sig->push_back(e9x9_30SigP.first);
1934  t_e11x11_30Sig->push_back(e11x11_30SigP.first);
1935  t_e15x15_30Sig->push_back(e15x15_30SigP.first);
1936 
1937  if (doMC_) {
1938  t_esim7x7->push_back(simInfo7x7.eTotal);
1939  t_esim9x9->push_back(simInfo9x9.eTotal);
1940  t_esim11x11->push_back(simInfo11x11.eTotal);
1941  t_esim15x15->push_back(simInfo15x15.eTotal);
1942 
1943  t_esim7x7Matched->push_back(simInfo7x7.eMatched);
1944  t_esim9x9Matched->push_back(simInfo9x9.eMatched);
1945  t_esim11x11Matched->push_back(simInfo11x11.eMatched);
1946  t_esim15x15Matched->push_back(simInfo15x15.eMatched);
1947 
1948  t_esim7x7Rest->push_back(simInfo7x7.eRest);
1949  t_esim9x9Rest->push_back(simInfo9x9.eRest);
1950  t_esim11x11Rest->push_back(simInfo11x11.eRest);
1951  t_esim15x15Rest->push_back(simInfo15x15.eRest);
1952 
1953  t_esim7x7Photon->push_back(simInfo7x7.eGamma);
1954  t_esim9x9Photon->push_back(simInfo9x9.eGamma);
1955  t_esim11x11Photon->push_back(simInfo11x11.eGamma);
1956  t_esim15x15Photon->push_back(simInfo15x15.eGamma);
1957 
1958  t_esim7x7NeutHad->push_back(simInfo7x7.eNeutralHad);
1959  t_esim9x9NeutHad->push_back(simInfo9x9.eNeutralHad);
1960  t_esim11x11NeutHad->push_back(simInfo11x11.eNeutralHad);
1961  t_esim15x15NeutHad->push_back(simInfo15x15.eNeutralHad);
1962 
1963  t_esim7x7CharHad->push_back(simInfo7x7.eChargedHad);
1964  t_esim9x9CharHad->push_back(simInfo9x9.eChargedHad);
1965  t_esim11x11CharHad->push_back(simInfo11x11.eChargedHad);
1966  t_esim15x15CharHad->push_back(simInfo15x15.eChargedHad);
1967 
1968  t_trkEcalEne->push_back(trkEcalEne);
1969  t_simTrackP->push_back(simTrackP);
1970  t_esimPdgId->push_back(simInfo11x11.pdgMatched);
1971  }
1972 
1973  t_maxNearHcalP3x3->push_back(maxNearHcalP3x3);
1974  t_maxNearHcalP5x5->push_back(maxNearHcalP5x5);
1975  t_maxNearHcalP7x7->push_back(maxNearHcalP7x7);
1976 
1977  t_h3x3->push_back(h3x3);
1978  t_h5x5->push_back(h5x5);
1979  t_h7x7->push_back(h7x7);
1980  t_h3x3Sig->push_back(h3x3Sig);
1981  t_h5x5Sig->push_back(h5x5Sig);
1982  t_h7x7Sig->push_back(h7x7Sig);
1983 
1984  t_infoHcal->push_back(trkDetItr->okHCAL);
1985  if (doMC_) {
1986  t_trkHcalEne->push_back(hcalScale * trkHcalEne);
1987 
1988  t_hsim3x3->push_back(hcalScale * hsimInfo3x3.eTotal);
1989  t_hsim5x5->push_back(hcalScale * hsimInfo5x5.eTotal);
1990  t_hsim7x7->push_back(hcalScale * hsimInfo7x7.eTotal);
1991 
1992  t_hsim3x3Matched->push_back(hcalScale * hsimInfo3x3.eMatched);
1993  t_hsim5x5Matched->push_back(hcalScale * hsimInfo5x5.eMatched);
1994  t_hsim7x7Matched->push_back(hcalScale * hsimInfo7x7.eMatched);
1995 
1996  t_hsim3x3Rest->push_back(hcalScale * hsimInfo3x3.eRest);
1997  t_hsim5x5Rest->push_back(hcalScale * hsimInfo5x5.eRest);
1998  t_hsim7x7Rest->push_back(hcalScale * hsimInfo7x7.eRest);
1999 
2000  t_hsim3x3Photon->push_back(hcalScale * hsimInfo3x3.eGamma);
2001  t_hsim5x5Photon->push_back(hcalScale * hsimInfo5x5.eGamma);
2002  t_hsim7x7Photon->push_back(hcalScale * hsimInfo7x7.eGamma);
2003 
2004  t_hsim3x3NeutHad->push_back(hcalScale * hsimInfo3x3.eNeutralHad);
2005  t_hsim5x5NeutHad->push_back(hcalScale * hsimInfo5x5.eNeutralHad);
2006  t_hsim7x7NeutHad->push_back(hcalScale * hsimInfo7x7.eNeutralHad);
2007 
2008  t_hsim3x3CharHad->push_back(hcalScale * hsimInfo3x3.eChargedHad);
2009  t_hsim5x5CharHad->push_back(hcalScale * hsimInfo5x5.eChargedHad);
2010  t_hsim7x7CharHad->push_back(hcalScale * hsimInfo7x7.eChargedHad);
2011  }
2012 
2013  } // if loosely isolated track
2014  } // check p1/eta
2015  } // loop over track collection
2016 
2017  if (haveIsoTrack)
2018  tree_->Fill();
2019 }
2020 
2021 // ----- method called once each job just before starting event loop ----
2022 void IsolatedTracksNxN::beginJob() {
2023  nEventProc_ = 0;
2024  nbad_ = 0;
2025  initL1_ = false;
2026  double tempgen_TH[NPBins + 1] = {
2027  0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 9.0, 11.0, 15.0, 20.0, 30.0, 50.0, 75.0, 100.0};
2028 
2029  for (unsigned int i = 0; i < NPBins + 1; i++)
2030  genPartPBins[i] = tempgen_TH[i];
2031 
2032  double tempgen_Eta[NEtaBins + 1] = {0.0, 1.131, 1.653, 2.172};
2033 
2034  for (unsigned int i = 0; i < NEtaBins + 1; i++)
2035  genPartEtaBins[i] = tempgen_Eta[i];
2036 
2037  bookHistograms();
2038 }
2039 
2040 // ----- method called once each job just after ending the event loop ----
2041 void IsolatedTracksNxN::endJob() {
2042  if (L1TriggerAlgoInfo_) {
2043  std::map<std::pair<unsigned int, std::string>, int>::iterator itr;
2044  for (itr = l1AlgoMap_.begin(); itr != l1AlgoMap_.end(); itr++) {
2045  edm::LogVerbatim("IsoTrack") << " ****endjob**** " << (itr->first).first << " " << (itr->first).second << " "
2046  << itr->second;
2047  int ibin = (itr->first).first;
2048  TString name((itr->first).second);
2049  h_L1AlgoNames->GetXaxis()->SetBinLabel(ibin + 1, name);
2050  }
2051  edm::LogVerbatim("IsoTrack") << "Number of Events Processed " << nEventProc_;
2052  }
2053 }
2054 
2055 //========================================================================================================
2056 
2057 void IsolatedTracksNxN::clearTreeVectors() {
2058  t_PVx->clear();
2059  t_PVy->clear();
2060  t_PVz->clear();
2061  t_PVisValid->clear();
2062  t_PVndof->clear();
2063  t_PVNTracks->clear();
2064  t_PVNTracksWt->clear();
2065  t_PVTracksSumPt->clear();
2066  t_PVTracksSumPtWt->clear();
2067  t_PVNTracksHP->clear();
2068  t_PVNTracksHPWt->clear();
2069  t_PVTracksSumPtHP->clear();
2070  t_PVTracksSumPtHPWt->clear();
2071 
2072  for (int i = 0; i < 128; i++)
2073  t_L1Decision[i] = 0;
2074  t_L1AlgoNames->clear();
2075  t_L1PreScale->clear();
2076 
2077  t_L1CenJetPt->clear();
2078  t_L1CenJetEta->clear();
2079  t_L1CenJetPhi->clear();
2080  t_L1FwdJetPt->clear();
2081  t_L1FwdJetEta->clear();
2082  t_L1FwdJetPhi->clear();
2083  t_L1TauJetPt->clear();
2084  t_L1TauJetEta->clear();
2085  t_L1TauJetPhi->clear();
2086  t_L1MuonPt->clear();
2087  t_L1MuonEta->clear();
2088  t_L1MuonPhi->clear();
2089  t_L1IsoEMPt->clear();
2090  t_L1IsoEMEta->clear();
2091  t_L1IsoEMPhi->clear();
2092  t_L1NonIsoEMPt->clear();
2093  t_L1NonIsoEMEta->clear();
2094  t_L1NonIsoEMPhi->clear();
2095  t_L1METPt->clear();
2096  t_L1METEta->clear();
2097  t_L1METPhi->clear();
2098 
2099  t_jetPt->clear();
2100  t_jetEta->clear();
2101  t_jetPhi->clear();
2102  t_nTrksJetCalo->clear();
2103  t_nTrksJetVtx->clear();
2104 
2105  t_trackPAll->clear();
2106  t_trackEtaAll->clear();
2107  t_trackPhiAll->clear();
2108  t_trackPdgIdAll->clear();
2109  t_trackPtAll->clear();
2110  t_trackDxyAll->clear();
2111  t_trackDzAll->clear();
2112  t_trackDxyPVAll->clear();
2113  t_trackDzPVAll->clear();
2114  t_trackChiSqAll->clear();
2115 
2116  t_trackP->clear();
2117  t_trackPt->clear();
2118  t_trackEta->clear();
2119  t_trackPhi->clear();
2120  t_trackEcalEta->clear();
2121  t_trackEcalPhi->clear();
2122  t_trackHcalEta->clear();
2123  t_trackHcalPhi->clear();
2124  t_NLayersCrossed->clear();
2125  t_trackNOuterHits->clear();
2126  t_trackDxy->clear();
2127  t_trackDxyBS->clear();
2128  t_trackDz->clear();
2129  t_trackDzBS->clear();
2130  t_trackDxyPV->clear();
2131  t_trackDzPV->clear();
2132  t_trackChiSq->clear();
2133  t_trackPVIdx->clear();
2134  t_trackHitsTOB->clear();
2135  t_trackHitsTEC->clear();
2136  t_trackHitInMissTOB->clear();
2137  t_trackHitInMissTEC->clear();
2138  t_trackHitInMissTIB->clear();
2139  t_trackHitInMissTID->clear();
2140  t_trackHitInMissTIBTID->clear();
2141  t_trackHitOutMissTOB->clear();
2142  t_trackHitOutMissTEC->clear();
2143  t_trackHitOutMissTIB->clear();
2144  t_trackHitOutMissTID->clear();
2145  t_trackHitOutMissTOBTEC->clear();
2146  t_trackHitInMeasTOB->clear();
2147  t_trackHitInMeasTEC->clear();
2148  t_trackHitInMeasTIB->clear();
2149  t_trackHitInMeasTID->clear();
2150  t_trackHitOutMeasTOB->clear();
2151  t_trackHitOutMeasTEC->clear();
2152  t_trackHitOutMeasTIB->clear();
2153  t_trackHitOutMeasTID->clear();
2154  t_trackOutPosOutHitDr->clear();
2155  t_trackL->clear();
2156 
2157  t_maxNearP31x31->clear();
2158  t_maxNearP21x21->clear();
2159 
2160  t_ecalSpike11x11->clear();
2161  t_e7x7->clear();
2162  t_e9x9->clear();
2163  t_e11x11->clear();
2164  t_e15x15->clear();
2165 
2166  t_e7x7_10Sig->clear();
2167  t_e9x9_10Sig->clear();
2168  t_e11x11_10Sig->clear();
2169  t_e15x15_10Sig->clear();
2170  t_e7x7_15Sig->clear();
2171  t_e9x9_15Sig->clear();
2172  t_e11x11_15Sig->clear();
2173  t_e15x15_15Sig->clear();
2174  t_e7x7_20Sig->clear();
2175  t_e9x9_20Sig->clear();
2176  t_e11x11_20Sig->clear();
2177  t_e15x15_20Sig->clear();
2178  t_e7x7_25Sig->clear();
2179  t_e9x9_25Sig->clear();
2180  t_e11x11_25Sig->clear();
2181  t_e15x15_25Sig->clear();
2182  t_e7x7_30Sig->clear();
2183  t_e9x9_30Sig->clear();
2184  t_e11x11_30Sig->clear();
2185  t_e15x15_30Sig->clear();
2186 
2187  if (doMC_) {
2188  t_simTrackP->clear();
2189  t_esimPdgId->clear();
2190  t_trkEcalEne->clear();
2191 
2192  t_esim7x7->clear();
2193  t_esim9x9->clear();
2194  t_esim11x11->clear();
2195  t_esim15x15->clear();
2196 
2197  t_esim7x7Matched->clear();
2198  t_esim9x9Matched->clear();
2199  t_esim11x11Matched->clear();
2200  t_esim15x15Matched->clear();
2201 
2202  t_esim7x7Rest->clear();
2203  t_esim9x9Rest->clear();
2204  t_esim11x11Rest->clear();
2205  t_esim15x15Rest->clear();
2206 
2207  t_esim7x7Photon->clear();
2208  t_esim9x9Photon->clear();
2209  t_esim11x11Photon->clear();
2210  t_esim15x15Photon->clear();
2211 
2212  t_esim7x7NeutHad->clear();
2213  t_esim9x9NeutHad->clear();
2214  t_esim11x11NeutHad->clear();
2215  t_esim15x15NeutHad->clear();
2216 
2217  t_esim7x7CharHad->clear();
2218  t_esim9x9CharHad->clear();
2219  t_esim11x11CharHad->clear();
2220  t_esim15x15CharHad->clear();
2221  }
2222 
2223  t_maxNearHcalP3x3->clear();
2224  t_maxNearHcalP5x5->clear();
2225  t_maxNearHcalP7x7->clear();
2226 
2227  t_h3x3->clear();
2228  t_h5x5->clear();
2229  t_h7x7->clear();
2230  t_h3x3Sig->clear();
2231  t_h5x5Sig->clear();
2232  t_h7x7Sig->clear();
2233 
2234  t_infoHcal->clear();
2235 
2236  if (doMC_) {
2237  t_trkHcalEne->clear();
2238 
2239  t_hsim3x3->clear();
2240  t_hsim5x5->clear();
2241  t_hsim7x7->clear();
2242  t_hsim3x3Matched->clear();
2243  t_hsim5x5Matched->clear();
2244  t_hsim7x7Matched->clear();
2245  t_hsim3x3Rest->clear();
2246  t_hsim5x5Rest->clear();
2247  t_hsim7x7Rest->clear();
2248  t_hsim3x3Photon->clear();
2249  t_hsim5x5Photon->clear();
2250  t_hsim7x7Photon->clear();
2251  t_hsim3x3NeutHad->clear();
2252  t_hsim5x5NeutHad->clear();
2253  t_hsim7x7NeutHad->clear();
2254  t_hsim3x3CharHad->clear();
2255  t_hsim5x5CharHad->clear();
2256  t_hsim7x7CharHad->clear();
2257  }
2258 }
2259 
2260 void IsolatedTracksNxN::bookHistograms() {
2261  char hname[100], htit[100];
2262 
2263  edm::Service<TFileService> fs;
2264  TFileDirectory dir = fs->mkdir("nearMaxTrackP");
2265 
2266  for (unsigned int ieta = 0; ieta < NEtaBins; ieta++) {
2267  double lowEta = -5.0, highEta = 5.0;
2268  lowEta = genPartEtaBins[ieta];
2269  highEta = genPartEtaBins[ieta + 1];
2270 
2271  for (unsigned int ipt = 0; ipt < NPBins; ipt++) {
2272  double lowP = 0.0, highP = 300.0;
2273  lowP = genPartPBins[ipt];
2274  highP = genPartPBins[ipt + 1];
2275  sprintf(hname, "h_maxNearP31x31_ptBin%i_etaBin%i", ipt, ieta);
2276  sprintf(htit, "maxNearP in 31x31 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP);
2277  h_maxNearP31x31[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
2278  h_maxNearP31x31[ipt][ieta]->Sumw2();
2279  sprintf(hname, "h_maxNearP25x25_ptBin%i_etaBin%i", ipt, ieta);
2280  sprintf(htit, "maxNearP in 25x25 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP);
2281  h_maxNearP25x25[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
2282  h_maxNearP25x25[ipt][ieta]->Sumw2();
2283  sprintf(hname, "h_maxNearP21x21_ptBin%i_etaBin%i", ipt, ieta);
2284  sprintf(htit, "maxNearP in 21x21 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP);
2285  h_maxNearP21x21[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
2286  h_maxNearP21x21[ipt][ieta]->Sumw2();
2287  sprintf(hname, "h_maxNearP15x15_ptBin%i_etaBin%i", ipt, ieta);
2288  sprintf(htit, "maxNearP in 15x15 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP);
2289  h_maxNearP15x15[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
2290  h_maxNearP15x15[ipt][ieta]->Sumw2();
2291  }
2292  }
2293 
2294  h_L1AlgoNames = fs->make<TH1I>("h_L1AlgoNames", "h_L1AlgoNames:Bin Labels", 128, -0.5, 127.5);
2295 
2296  // Reconstructed Tracks
2297 
2298  h_PVTracksWt = fs->make<TH1F>("h_PVTracksWt", "h_PVTracksWt", 600, -0.1, 1.1);
2299 
2300  h_nTracks = fs->make<TH1F>("h_nTracks", "h_nTracks", 1000, -0.5, 999.5);
2301 
2302  sprintf(hname, "h_recEtaPt_0");
2303  sprintf(htit, "h_recEtaPt (all tracks Eta vs pT)");
2304  h_recEtaPt_0 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
2305 
2306  sprintf(hname, "h_recEtaP_0");
2307  sprintf(htit, "h_recEtaP (all tracks Eta vs pT)");
2308  h_recEtaP_0 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
2309 
2310  h_recPt_0 = fs->make<TH1F>("h_recPt_0", "Pt (all tracks)", 15, genPartPBins);
2311  h_recP_0 = fs->make<TH1F>("h_recP_0", "P (all tracks)", 15, genPartPBins);
2312  h_recEta_0 = fs->make<TH1F>("h_recEta_0", "Eta (all tracks)", 60, -3.0, 3.0);
2313  h_recPhi_0 = fs->make<TH1F>("h_recPhi_0", "Phi (all tracks)", 100, -3.2, 3.2);
2314  //-------------------------
2315  sprintf(hname, "h_recEtaPt_1");
2316  sprintf(htit, "h_recEtaPt (all good tracks Eta vs pT)");
2317  h_recEtaPt_1 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
2318 
2319  sprintf(hname, "h_recEtaP_1");
2320  sprintf(htit, "h_recEtaP (all good tracks Eta vs pT)");
2321  h_recEtaP_1 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
2322 
2323  h_recPt_1 = fs->make<TH1F>("h_recPt_1", "Pt (all good tracks)", 15, genPartPBins);
2324  h_recP_1 = fs->make<TH1F>("h_recP_1", "P (all good tracks)", 15, genPartPBins);
2325  h_recEta_1 = fs->make<TH1F>("h_recEta_1", "Eta (all good tracks)", 60, -3.0, 3.0);
2326  h_recPhi_1 = fs->make<TH1F>("h_recPhi_1", "Phi (all good tracks)", 100, -3.2, 3.2);
2327  //-------------------------
2328  sprintf(hname, "h_recEtaPt_2");
2329  sprintf(htit, "h_recEtaPt (charge isolation Eta vs pT)");
2330  h_recEtaPt_2 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
2331 
2332  sprintf(hname, "h_recEtaP_2");
2333  sprintf(htit, "h_recEtaP (charge isolation Eta vs pT)");
2334  h_recEtaP_2 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
2335 
2336  h_recPt_2 = fs->make<TH1F>("h_recPt_2", "Pt (charge isolation)", 15, genPartPBins);
2337  h_recP_2 = fs->make<TH1F>("h_recP_2", "P (charge isolation)", 15, genPartPBins);
2338  h_recEta_2 = fs->make<TH1F>("h_recEta_2", "Eta (charge isolation)", 60, -3.0, 3.0);
2339  h_recPhi_2 = fs->make<TH1F>("h_recPhi_2", "Phi (charge isolation)", 100, -3.2, 3.2);
2340 
2341  tree_ = fs->make<TTree>("tree", "tree");
2342  tree_->SetAutoSave(10000);
2343 
2344  tree_->Branch("t_EvtNo", &t_EvtNo, "t_EvtNo/I");
2345  tree_->Branch("t_RunNo", &t_RunNo, "t_RunNo/I");
2346  tree_->Branch("t_Lumi", &t_Lumi, "t_Lumi/I");
2347  tree_->Branch("t_Bunch", &t_Bunch, "t_Bunch/I");
2348 
2349  t_PVx = new std::vector<double>();
2350  t_PVy = new std::vector<double>();
2351  t_PVz = new std::vector<double>();
2352  t_PVisValid = new std::vector<int>();
2353  t_PVndof = new std::vector<int>();
2354  t_PVNTracks = new std::vector<int>();
2355  t_PVNTracksWt = new std::vector<int>();
2356  t_PVTracksSumPt = new std::vector<double>();
2357  t_PVTracksSumPtWt = new std::vector<double>();
2358  t_PVNTracksHP = new std::vector<int>();
2359  t_PVNTracksHPWt = new std::vector<int>();
2360  t_PVTracksSumPtHP = new std::vector<double>();
2361  t_PVTracksSumPtHPWt = new std::vector<double>();
2362 
2363  tree_->Branch("PVx", "std::vector<double>", &t_PVx);
2364  tree_->Branch("PVy", "std::vector<double>", &t_PVy);
2365  tree_->Branch("PVz", "std::vector<double>", &t_PVz);
2366  tree_->Branch("PVisValid", "std::vector<int>", &t_PVisValid);
2367  tree_->Branch("PVndof", "std::vector<int>", &t_PVndof);
2368  tree_->Branch("PVNTracks", "std::vector<int>", &t_PVNTracks);
2369  tree_->Branch("PVNTracksWt", "std::vector<int>", &t_PVNTracksWt);
2370  tree_->Branch("t_PVTracksSumPt", "std::vector<double>", &t_PVTracksSumPt);
2371  tree_->Branch("t_PVTracksSumPtWt", "std::vector<double>", &t_PVTracksSumPtWt);
2372  tree_->Branch("PVNTracksHP", "std::vector<int>", &t_PVNTracksHP);
2373  tree_->Branch("PVNTracksHPWt", "std::vector<int>", &t_PVNTracksHPWt);
2374  tree_->Branch("t_PVTracksSumPtHP", "std::vector<double>", &t_PVTracksSumPtHP);
2375  tree_->Branch("t_PVTracksSumPtHPWt", "std::vector<double>", &t_PVTracksSumPtHPWt);
2376 
2377  //----- L1Trigger information
2378  for (int i = 0; i < 128; i++)
2379  t_L1Decision[i] = 0;
2380  t_L1AlgoNames = new std::vector<std::string>();
2381  t_L1PreScale = new std::vector<int>();
2382  t_L1CenJetPt = new std::vector<double>();
2383  t_L1CenJetEta = new std::vector<double>();
2384  t_L1CenJetPhi = new std::vector<double>();
2385  t_L1FwdJetPt = new std::vector<double>();
2386  t_L1FwdJetEta = new std::vector<double>();
2387  t_L1FwdJetPhi = new std::vector<double>();
2388  t_L1TauJetPt = new std::vector<double>();
2389  t_L1TauJetEta = new std::vector<double>();
2390  t_L1TauJetPhi = new std::vector<double>();
2391  t_L1MuonPt = new std::vector<double>();
2392  t_L1MuonEta = new std::vector<double>();
2393  t_L1MuonPhi = new std::vector<double>();
2394  t_L1IsoEMPt = new std::vector<double>();
2395  t_L1IsoEMEta = new std::vector<double>();
2396  t_L1IsoEMPhi = new std::vector<double>();
2397  t_L1NonIsoEMPt = new std::vector<double>();
2398  t_L1NonIsoEMEta = new std::vector<double>();
2399  t_L1NonIsoEMPhi = new std::vector<double>();
2400  t_L1METPt = new std::vector<double>();
2401  t_L1METEta = new std::vector<double>();
2402  t_L1METPhi = new std::vector<double>();
2403 
2404  tree_->Branch("t_L1Decision", t_L1Decision, "t_L1Decision[128]/I");
2405  tree_->Branch("t_L1AlgoNames", "std::vector<string>", &t_L1AlgoNames);
2406  tree_->Branch("t_L1PreScale", "std::vector<int>", &t_L1PreScale);
2407  tree_->Branch("t_L1CenJetPt", "std::vector<double>", &t_L1CenJetPt);
2408  tree_->Branch("t_L1CenJetEta", "std::vector<double>", &t_L1CenJetEta);
2409  tree_->Branch("t_L1CenJetPhi", "std::vector<double>", &t_L1CenJetPhi);
2410  tree_->Branch("t_L1FwdJetPt", "std::vector<double>", &t_L1FwdJetPt);
2411  tree_->Branch("t_L1FwdJetEta", "std::vector<double>", &t_L1FwdJetEta);
2412  tree_->Branch("t_L1FwdJetPhi", "std::vector<double>", &t_L1FwdJetPhi);
2413  tree_->Branch("t_L1TauJetPt", "std::vector<double>", &t_L1TauJetPt);
2414  tree_->Branch("t_L1TauJetEta", "std::vector<double>", &t_L1TauJetEta);
2415  tree_->Branch("t_L1TauJetPhi", "std::vector<double>", &t_L1TauJetPhi);
2416  tree_->Branch("t_L1MuonPt", "std::vector<double>", &t_L1MuonPt);
2417  tree_->Branch("t_L1MuonEta", "std::vector<double>", &t_L1MuonEta);
2418  tree_->Branch("t_L1MuonPhi", "std::vector<double>", &t_L1MuonPhi);
2419  tree_->Branch("t_L1IsoEMPt", "std::vector<double>", &t_L1IsoEMPt);
2420  tree_->Branch("t_L1IsoEMEta", "std::vector<double>", &t_L1IsoEMEta);
2421  tree_->Branch("t_L1IsoEMPhi", "std::vector<double>", &t_L1IsoEMPhi);
2422  tree_->Branch("t_L1NonIsoEMPt", "std::vector<double>", &t_L1NonIsoEMPt);
2423  tree_->Branch("t_L1NonIsoEMEta", "std::vector<double>", &t_L1NonIsoEMEta);
2424  tree_->Branch("t_L1NonIsoEMPhi", "std::vector<double>", &t_L1NonIsoEMPhi);
2425  tree_->Branch("t_L1METPt", "std::vector<double>", &t_L1METPt);
2426  tree_->Branch("t_L1METEta", "std::vector<double>", &t_L1METEta);
2427  tree_->Branch("t_L1METPhi", "std::vector<double>", &t_L1METPhi);
2428 
2429  t_jetPt = new std::vector<double>();
2430  t_jetEta = new std::vector<double>();
2431  t_jetPhi = new std::vector<double>();
2432  t_nTrksJetCalo = new std::vector<double>();
2433  t_nTrksJetVtx = new std::vector<double>();
2434  tree_->Branch("t_jetPt", "std::vector<double>", &t_jetPt);
2435  tree_->Branch("t_jetEta", "std::vector<double>", &t_jetEta);
2436  tree_->Branch("t_jetPhi", "std::vector<double>", &t_jetPhi);
2437  tree_->Branch("t_nTrksJetCalo", "std::vector<double>", &t_nTrksJetCalo);
2438  tree_->Branch("t_nTrksJetVtx", "std::vector<double>", &t_nTrksJetVtx);
2439 
2440  t_trackPAll = new std::vector<double>();
2441  t_trackEtaAll = new std::vector<double>();
2442  t_trackPhiAll = new std::vector<double>();
2443  t_trackPdgIdAll = new std::vector<double>();
2444  t_trackPtAll = new std::vector<double>();
2445  t_trackDxyAll = new std::vector<double>();
2446  t_trackDzAll = new std::vector<double>();
2447  t_trackDxyPVAll = new std::vector<double>();
2448  t_trackDzPVAll = new std::vector<double>();
2449  t_trackChiSqAll = new std::vector<double>();
2450  tree_->Branch("t_trackPAll", "std::vector<double>", &t_trackPAll);
2451  tree_->Branch("t_trackPhiAll", "std::vector<double>", &t_trackPhiAll);
2452  tree_->Branch("t_trackEtaAll", "std::vector<double>", &t_trackEtaAll);
2453  tree_->Branch("t_trackPtAll", "std::vector<double>", &t_trackPtAll);
2454  tree_->Branch("t_trackDxyAll", "std::vector<double>", &t_trackDxyAll);
2455  tree_->Branch("t_trackDzAll", "std::vector<double>", &t_trackDzAll);
2456  tree_->Branch("t_trackDxyPVAll", "std::vector<double>", &t_trackDxyPVAll);
2457  tree_->Branch("t_trackDzPVAll", "std::vector<double>", &t_trackDzPVAll);
2458  tree_->Branch("t_trackChiSqAll", "std::vector<double>", &t_trackChiSqAll);
2459  //tree_->Branch("t_trackPdgIdAll", "std::vector<double>", &t_trackPdgIdAll);
2460 
2461  t_trackP = new std::vector<double>();
2462  t_trackPt = new std::vector<double>();
2463  t_trackEta = new std::vector<double>();
2464  t_trackPhi = new std::vector<double>();
2465  t_trackEcalEta = new std::vector<double>();
2466  t_trackEcalPhi = new std::vector<double>();
2467  t_trackHcalEta = new std::vector<double>();
2468  t_trackHcalPhi = new std::vector<double>();
2469  t_trackNOuterHits = new std::vector<int>();
2470  t_NLayersCrossed = new std::vector<int>();
2471  t_trackDxy = new std::vector<double>();
2472  t_trackDxyBS = new std::vector<double>();
2473  t_trackDz = new std::vector<double>();
2474  t_trackDzBS = new std::vector<double>();
2475  t_trackDxyPV = new std::vector<double>();
2476  t_trackDzPV = new std::vector<double>();
2477  t_trackPVIdx = new std::vector<int>();
2478  t_trackChiSq = new std::vector<double>();
2479  t_trackHitsTOB = new std::vector<int>();
2480  t_trackHitsTEC = new std::vector<int>();
2481  t_trackHitInMissTOB = new std::vector<int>();
2482  t_trackHitInMissTEC = new std::vector<int>();
2483  t_trackHitInMissTIB = new std::vector<int>();
2484  t_trackHitInMissTID = new std::vector<int>();
2485  t_trackHitOutMissTOB = new std::vector<int>();
2486  t_trackHitOutMissTEC = new std::vector<int>();
2487  t_trackHitOutMissTIB = new std::vector<int>();
2488  t_trackHitOutMissTID = new std::vector<int>();
2489  t_trackHitInMissTIBTID = new std::vector<int>();
2490  t_trackHitOutMissTOB = new std::vector<int>();
2491  t_trackHitOutMissTEC = new std::vector<int>();
2492  t_trackHitOutMissTIB = new std::vector<int>();
2493  t_trackHitOutMissTID = new std::vector<int>();
2494  t_trackHitOutMissTOBTEC = new std::vector<int>();
2495  t_trackHitInMeasTOB = new std::vector<int>();
2496  t_trackHitInMeasTEC = new std::vector<int>();
2497  t_trackHitInMeasTIB = new std::vector<int>();
2498  t_trackHitInMeasTID = new std::vector<int>();
2499  t_trackHitOutMeasTOB = new std::vector<int>();
2500  t_trackHitOutMeasTEC = new std::vector<int>();
2501  t_trackHitOutMeasTIB = new std::vector<int>();
2502  t_trackHitOutMeasTID = new std::vector<int>();
2503  t_trackOutPosOutHitDr = new std::vector<double>();
2504  t_trackL = new std::vector<double>();
2505 
2506  tree_->Branch("t_trackP", "std::vector<double>", &t_trackP);
2507  tree_->Branch("t_trackPt", "std::vector<double>", &t_trackPt);
2508  tree_->Branch("t_trackEta", "std::vector<double>", &t_trackEta);
2509  tree_->Branch("t_trackPhi", "std::vector<double>", &t_trackPhi);
2510  tree_->Branch("t_trackEcalEta", "std::vector<double>", &t_trackEcalEta);
2511  tree_->Branch("t_trackEcalPhi", "std::vector<double>", &t_trackEcalPhi);
2512  tree_->Branch("t_trackHcalEta", "std::vector<double>", &t_trackHcalEta);
2513  tree_->Branch("t_trackHcalPhi", "std::vector<double>", &t_trackHcalPhi);
2514 
2515  tree_->Branch("t_trackNOuterHits", "std::vector<int>", &t_trackNOuterHits);
2516  tree_->Branch("t_NLayersCrossed", "std::vector<int>", &t_NLayersCrossed);
2517  tree_->Branch("t_trackHitsTOB", "std::vector<int>", &t_trackHitsTOB);
2518  tree_->Branch("t_trackHitsTEC", "std::vector<int>", &t_trackHitsTEC);
2519  tree_->Branch("t_trackHitInMissTOB", "std::vector<int>", &t_trackHitInMissTOB);
2520  tree_->Branch("t_trackHitInMissTEC", "std::vector<int>", &t_trackHitInMissTEC);
2521  tree_->Branch("t_trackHitInMissTIB", "std::vector<int>", &t_trackHitInMissTIB);
2522  tree_->Branch("t_trackHitInMissTID", "std::vector<int>", &t_trackHitInMissTID);
2523  tree_->Branch("t_trackHitInMissTIBTID", "std::vector<int>", &t_trackHitInMissTIBTID);
2524  tree_->Branch("t_trackHitOutMissTOB", "std::vector<int>", &t_trackHitOutMissTOB);
2525  tree_->Branch("t_trackHitOutMissTEC", "std::vector<int>", &t_trackHitOutMissTEC);
2526  tree_->Branch("t_trackHitOutMissTIB", "std::vector<int>", &t_trackHitOutMissTIB);
2527  tree_->Branch("t_trackHitOutMissTID", "std::vector<int>", &t_trackHitOutMissTID);
2528  tree_->Branch("t_trackHitOutMissTOBTEC", "std::vector<int>", &t_trackHitOutMissTOBTEC);
2529  tree_->Branch("t_trackHitInMeasTOB", "std::vector<int>", &t_trackHitInMeasTOB);
2530  tree_->Branch("t_trackHitInMeasTEC", "std::vector<int>", &t_trackHitInMeasTEC);
2531  tree_->Branch("t_trackHitInMeasTIB", "std::vector<int>", &t_trackHitInMeasTIB);
2532  tree_->Branch("t_trackHitInMeasTID", "std::vector<int>", &t_trackHitInMeasTID);
2533  tree_->Branch("t_trackHitOutMeasTOB", "std::vector<int>", &t_trackHitOutMeasTOB);
2534  tree_->Branch("t_trackHitOutMeasTEC", "std::vector<int>", &t_trackHitOutMeasTEC);
2535  tree_->Branch("t_trackHitOutMeasTIB", "std::vector<int>", &t_trackHitOutMeasTIB);
2536  tree_->Branch("t_trackHitOutMeasTID", "std::vector<int>", &t_trackHitOutMeasTID);
2537  tree_->Branch("t_trackOutPosOutHitDr", "std::vector<double>", &t_trackOutPosOutHitDr);
2538  tree_->Branch("t_trackL", "std::vector<double>", &t_trackL);
2539 
2540  tree_->Branch("t_trackDxy", "std::vector<double>", &t_trackDxy);
2541  tree_->Branch("t_trackDxyBS", "std::vector<double>", &t_trackDxyBS);
2542  tree_->Branch("t_trackDz", "std::vector<double>", &t_trackDz);
2543  tree_->Branch("t_trackDzBS", "std::vector<double>", &t_trackDzBS);
2544  tree_->Branch("t_trackDxyPV", "std::vector<double>", &t_trackDxyPV);
2545  tree_->Branch("t_trackDzPV", "std::vector<double>", &t_trackDzPV);
2546  tree_->Branch("t_trackChiSq", "std::vector<double>", &t_trackChiSq);
2547  tree_->Branch("t_trackPVIdx", "std::vector<int>", &t_trackPVIdx);
2548 
2549  t_maxNearP31x31 = new std::vector<double>();
2550  t_maxNearP21x21 = new std::vector<double>();
2551 
2552  tree_->Branch("t_maxNearP31x31", "std::vector<double>", &t_maxNearP31x31);
2553  tree_->Branch("t_maxNearP21x21", "std::vector<double>", &t_maxNearP21x21);
2554 
2555  t_ecalSpike11x11 = new std::vector<int>();
2556  t_e7x7 = new std::vector<double>();
2557  t_e9x9 = new std::vector<double>();
2558  t_e11x11 = new std::vector<double>();
2559  t_e15x15 = new std::vector<double>();
2560 
2561  tree_->Branch("t_ecalSpike11x11", "std::vector<int>", &t_ecalSpike11x11);
2562  tree_->Branch("t_e7x7", "std::vector<double>", &t_e7x7);
2563  tree_->Branch("t_e9x9", "std::vector<double>", &t_e9x9);
2564  tree_->Branch("t_e11x11", "std::vector<double>", &t_e11x11);
2565  tree_->Branch("t_e15x15", "std::vector<double>", &t_e15x15);
2566 
2567  t_e7x7_10Sig = new std::vector<double>();
2568  t_e9x9_10Sig = new std::vector<double>();
2569  t_e11x11_10Sig = new std::vector<double>();
2570  t_e15x15_10Sig = new std::vector<double>();
2571  t_e7x7_15Sig = new std::vector<double>();
2572  t_e9x9_15Sig = new std::vector<double>();
2573  t_e11x11_15Sig = new std::vector<double>();
2574  t_e15x15_15Sig = new std::vector<double>();
2575  t_e7x7_20Sig = new std::vector<double>();
2576  t_e9x9_20Sig = new std::vector<double>();
2577  t_e11x11_20Sig = new std::vector<double>();
2578  t_e15x15_20Sig = new std::vector<double>();
2579  t_e7x7_25Sig = new std::vector<double>();
2580  t_e9x9_25Sig = new std::vector<double>();
2581  t_e11x11_25Sig = new std::vector<double>();
2582  t_e15x15_25Sig = new std::vector<double>();
2583  t_e7x7_30Sig = new std::vector<double>();
2584  t_e9x9_30Sig = new std::vector<double>();
2585  t_e11x11_30Sig = new std::vector<double>();
2586  t_e15x15_30Sig = new std::vector<double>();
2587 
2588  tree_->Branch("t_e7x7_10Sig", "std::vector<double>", &t_e7x7_10Sig);
2589  tree_->Branch("t_e9x9_10Sig", "std::vector<double>", &t_e9x9_10Sig);
2590  tree_->Branch("t_e11x11_10Sig", "std::vector<double>", &t_e11x11_10Sig);
2591  tree_->Branch("t_e15x15_10Sig", "std::vector<double>", &t_e15x15_10Sig);
2592  tree_->Branch("t_e7x7_15Sig", "std::vector<double>", &t_e7x7_15Sig);
2593  tree_->Branch("t_e9x9_15Sig", "std::vector<double>", &t_e9x9_15Sig);
2594  tree_->Branch("t_e11x11_15Sig", "std::vector<double>", &t_e11x11_15Sig);
2595  tree_->Branch("t_e15x15_15Sig", "std::vector<double>", &t_e15x15_15Sig);
2596  tree_->Branch("t_e7x7_20Sig", "std::vector<double>", &t_e7x7_20Sig);
2597  tree_->Branch("t_e9x9_20Sig", "std::vector<double>", &t_e9x9_20Sig);
2598  tree_->Branch("t_e11x11_20Sig", "std::vector<double>", &t_e11x11_20Sig);
2599  tree_->Branch("t_e15x15_20Sig", "std::vector<double>", &t_e15x15_20Sig);
2600  tree_->Branch("t_e7x7_25Sig", "std::vector<double>", &t_e7x7_25Sig);
2601  tree_->Branch("t_e9x9_25Sig", "std::vector<double>", &t_e9x9_25Sig);
2602  tree_->Branch("t_e11x11_25Sig", "std::vector<double>", &t_e11x11_25Sig);
2603  tree_->Branch("t_e15x15_25Sig", "std::vector<double>", &t_e15x15_25Sig);
2604  tree_->Branch("t_e7x7_30Sig", "std::vector<double>", &t_e7x7_30Sig);
2605  tree_->Branch("t_e9x9_30Sig", "std::vector<double>", &t_e9x9_30Sig);
2606  tree_->Branch("t_e11x11_30Sig", "std::vector<double>", &t_e11x11_30Sig);
2607  tree_->Branch("t_e15x15_30Sig", "std::vector<double>", &t_e15x15_30Sig);
2608 
2609  if (doMC_) {
2610  t_esim7x7 = new std::vector<double>();
2611  t_esim9x9 = new std::vector<double>();
2612  t_esim11x11 = new std::vector<double>();
2613  t_esim15x15 = new std::vector<double>();
2614 
2615  t_esim7x7Matched = new std::vector<double>();
2616  t_esim9x9Matched = new std::vector<double>();
2617  t_esim11x11Matched = new std::vector<double>();
2618  t_esim15x15Matched = new std::vector<double>();
2619 
2620  t_esim7x7Rest = new std::vector<double>();
2621  t_esim9x9Rest = new std::vector<double>();
2622  t_esim11x11Rest = new std::vector<double>();
2623  t_esim15x15Rest = new std::vector<double>();
2624 
2625  t_esim7x7Photon = new std::vector<double>();
2626  t_esim9x9Photon = new std::vector<double>();
2627  t_esim11x11Photon = new std::vector<double>();
2628  t_esim15x15Photon = new std::vector<double>();
2629 
2630  t_esim7x7NeutHad = new std::vector<double>();
2631  t_esim9x9NeutHad = new std::vector<double>();
2632  t_esim11x11NeutHad = new std::vector<double>();
2633  t_esim15x15NeutHad = new std::vector<double>();
2634 
2635  t_esim7x7CharHad = new std::vector<double>();
2636  t_esim9x9CharHad = new std::vector<double>();
2637  t_esim11x11CharHad = new std::vector<double>();
2638  t_esim15x15CharHad = new std::vector<double>();
2639 
2640  t_trkEcalEne = new std::vector<double>();
2641  t_simTrackP = new std::vector<double>();
2642  t_esimPdgId = new std::vector<double>();
2643 
2644  tree_->Branch("t_esim7x7", "std::vector<double>", &t_esim7x7);
2645  tree_->Branch("t_esim9x9", "std::vector<double>", &t_esim9x9);
2646  tree_->Branch("t_esim11x11", "std::vector<double>", &t_esim11x11);
2647  tree_->Branch("t_esim15x15", "std::vector<double>", &t_esim15x15);
2648 
2649  tree_->Branch("t_esim7x7Matched", "std::vector<double>", &t_esim7x7Matched);
2650  tree_->Branch("t_esim9x9Matched", "std::vector<double>", &t_esim9x9Matched);
2651  tree_->Branch("t_esim11x11Matched", "std::vector<double>", &t_esim11x11Matched);
2652  tree_->Branch("t_esim15x15Matched", "std::vector<double>", &t_esim15x15Matched);
2653 
2654  tree_->Branch("t_esim7x7Rest", "std::vector<double>", &t_esim7x7Rest);
2655  tree_->Branch("t_esim9x9Rest", "std::vector<double>", &t_esim9x9Rest);
2656  tree_->Branch("t_esim11x11Rest", "std::vector<double>", &t_esim11x11Rest);
2657  tree_->Branch("t_esim15x15Rest", "std::vector<double>", &t_esim15x15Rest);
2658 
2659  tree_->Branch("t_esim7x7Photon", "std::vector<double>", &t_esim7x7Photon);
2660  tree_->Branch("t_esim9x9Photon", "std::vector<double>", &t_esim9x9Photon);
2661  tree_->Branch("t_esim11x11Photon", "std::vector<double>", &t_esim11x11Photon);
2662  tree_->Branch("t_esim15x15Photon", "std::vector<double>", &t_esim15x15Photon);
2663 
2664  tree_->Branch("t_esim7x7NeutHad", "std::vector<double>", &t_esim7x7NeutHad);
2665  tree_->Branch("t_esim9x9NeutHad", "std::vector<double>", &t_esim9x9NeutHad);
2666  tree_->Branch("t_esim11x11NeutHad", "std::vector<double>", &t_esim11x11NeutHad);
2667  tree_->Branch("t_esim15x15NeutHad", "std::vector<double>", &t_esim15x15NeutHad);
2668 
2669  tree_->Branch("t_esim7x7CharHad", "std::vector<double>", &t_esim7x7CharHad);
2670  tree_->Branch("t_esim9x9CharHad", "std::vector<double>", &t_esim9x9CharHad);
2671  tree_->Branch("t_esim11x11CharHad", "std::vector<double>", &t_esim11x11CharHad);
2672  tree_->Branch("t_esim15x15CharHad", "std::vector<double>", &t_esim15x15CharHad);
2673 
2674  tree_->Branch("t_trkEcalEne", "std::vector<double>", &t_trkEcalEne);
2675  tree_->Branch("t_simTrackP", "std::vector<double>", &t_simTrackP);
2676  tree_->Branch("t_esimPdgId", "std::vector<double>", &t_esimPdgId);
2677  }
2678 
2679  t_maxNearHcalP3x3 = new std::vector<double>();
2680  t_maxNearHcalP5x5 = new std::vector<double>();
2681  t_maxNearHcalP7x7 = new std::vector<double>();
2682  t_h3x3 = new std::vector<double>();
2683  t_h5x5 = new std::vector<double>();
2684  t_h7x7 = new std::vector<double>();
2685  t_h3x3Sig = new std::vector<double>();
2686  t_h5x5Sig = new std::vector<double>();
2687  t_h7x7Sig = new std::vector<double>();
2688  t_infoHcal = new std::vector<int>();
2689 
2690  if (doMC_) {
2691  t_trkHcalEne = new std::vector<double>();
2692  t_hsim3x3 = new std::vector<double>();
2693  t_hsim5x5 = new std::vector<double>();
2694  t_hsim7x7 = new std::vector<double>();
2695  t_hsim3x3Matched = new std::vector<double>();
2696  t_hsim5x5Matched = new std::vector<double>();
2697  t_hsim7x7Matched = new std::vector<double>();
2698  t_hsim3x3Rest = new std::vector<double>();
2699  t_hsim5x5Rest = new std::vector<double>();
2700  t_hsim7x7Rest = new std::vector<double>();
2701  t_hsim3x3Photon = new std::vector<double>();
2702  t_hsim5x5Photon = new std::vector<double>();
2703  t_hsim7x7Photon = new std::vector<double>();
2704  t_hsim3x3NeutHad = new std::vector<double>();
2705  t_hsim5x5NeutHad = new std::vector<double>();
2706  t_hsim7x7NeutHad = new std::vector<double>();
2707  t_hsim3x3CharHad = new std::vector<double>();
2708  t_hsim5x5CharHad = new std::vector<double>();
2709  t_hsim7x7CharHad = new std::vector<double>();
2710  }
2711 
2712  tree_->Branch("t_maxNearHcalP3x3", "std::vector<double>", &t_maxNearHcalP3x3);
2713  tree_->Branch("t_maxNearHcalP5x5", "std::vector<double>", &t_maxNearHcalP5x5);
2714  tree_->Branch("t_maxNearHcalP7x7", "std::vector<double>", &t_maxNearHcalP7x7);
2715  tree_->Branch("t_h3x3", "std::vector<double>", &t_h3x3);
2716  tree_->Branch("t_h5x5", "std::vector<double>", &t_h5x5);
2717  tree_->Branch("t_h7x7", "std::vector<double>", &t_h7x7);
2718  tree_->Branch("t_h3x3Sig", "std::vector<double>", &t_h3x3Sig);
2719  tree_->Branch("t_h5x5Sig", "std::vector<double>", &t_h5x5Sig);
2720  tree_->Branch("t_h7x7Sig", "std::vector<double>", &t_h7x7Sig);
2721  tree_->Branch("t_infoHcal", "std::vector<int>", &t_infoHcal);
2722 
2723  if (doMC_) {
2724  tree_->Branch("t_trkHcalEne", "std::vector<double>", &t_trkHcalEne);
2725  tree_->Branch("t_hsim3x3", "std::vector<double>", &t_hsim3x3);
2726  tree_->Branch("t_hsim5x5", "std::vector<double>", &t_hsim5x5);
2727  tree_->Branch("t_hsim7x7", "std::vector<double>", &t_hsim7x7);
2728  tree_->Branch("t_hsim3x3Matched", "std::vector<double>", &t_hsim3x3Matched);
2729  tree_->Branch("t_hsim5x5Matched", "std::vector<double>", &t_hsim5x5Matched);
2730  tree_->Branch("t_hsim7x7Matched", "std::vector<double>", &t_hsim7x7Matched);
2731  tree_->Branch("t_hsim3x3Rest", "std::vector<double>", &t_hsim3x3Rest);
2732  tree_->Branch("t_hsim5x5Rest", "std::vector<double>", &t_hsim5x5Rest);
2733  tree_->Branch("t_hsim7x7Rest", "std::vector<double>", &t_hsim7x7Rest);
2734  tree_->Branch("t_hsim3x3Photon", "std::vector<double>", &t_hsim3x3Photon);
2735  tree_->Branch("t_hsim5x5Photon", "std::vector<double>", &t_hsim5x5Photon);
2736  tree_->Branch("t_hsim7x7Photon", "std::vector<double>", &t_hsim7x7Photon);
2737  tree_->Branch("t_hsim3x3NeutHad", "std::vector<double>", &t_hsim3x3NeutHad);
2738  tree_->Branch("t_hsim5x5NeutHad", "std::vector<double>", &t_hsim5x5NeutHad);
2739  tree_->Branch("t_hsim7x7NeutHad", "std::vector<double>", &t_hsim7x7NeutHad);
2740  tree_->Branch("t_hsim3x3CharHad", "std::vector<double>", &t_hsim3x3CharHad);
2741  tree_->Branch("t_hsim5x5CharHad", "std::vector<double>", &t_hsim5x5CharHad);
2742  tree_->Branch("t_hsim7x7CharHad", "std::vector<double>", &t_hsim7x7CharHad);
2743  }
2744  tree_->Branch("t_nTracks", &t_nTracks, "t_nTracks/I");
2745 }
2746 
2747 void IsolatedTracksNxN::printTrack(const reco::Track *pTrack) {
2748  std::string theTrackQuality = "highPurity";
2749  reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality);
2750 
2751  edm::LogVerbatim("IsoTrack") << " Reference Point " << pTrack->referencePoint() << "\n TrackMmentum "
2752  << pTrack->momentum() << " (pt,eta,phi)(" << pTrack->pt() << "," << pTrack->eta() << ","
2753  << pTrack->phi() << ")"
2754  << " p " << pTrack->p() << "\n Normalized chi2 " << pTrack->normalizedChi2()
2755  << " charge " << pTrack->charge() << " qoverp() " << pTrack->qoverp() << "+-"
2756  << pTrack->qoverpError() << " d0 " << pTrack->d0() << "\n NValidHits "
2757  << pTrack->numberOfValidHits() << " NLostHits " << pTrack->numberOfLostHits()
2758  << " TrackQuality " << pTrack->qualityName(trackQuality_) << " "
2759  << pTrack->quality(trackQuality_);
2760 
2761  if (printTrkHitPattern_) {
2762  const reco::HitPattern &p = pTrack->hitPattern();
2763 
2764  edm::LogVerbatim("IsoTrack") << "default ";
2765  for (int i = 0; i < p.numberOfAllHits(reco::HitPattern::TRACK_HITS); i++) {
2766  p.printHitPattern(reco::HitPattern::TRACK_HITS, i, std::cout);
2767  }
2768  edm::LogVerbatim("IsoTrack") << "trackerMissingInnerHits() ";
2769  for (int i = 0; i < p.numberOfAllHits(reco::HitPattern::MISSING_INNER_HITS); i++) {
2770  p.printHitPattern(reco::HitPattern::MISSING_INNER_HITS, i, std::cout);
2771  }
2772  edm::LogVerbatim("IsoTrack") << "trackerMissingOuterHits() ";
2773  for (int i = 0; i < p.numberOfAllHits(reco::HitPattern::MISSING_OUTER_HITS); i++) {
2774  p.printHitPattern(reco::HitPattern::MISSING_OUTER_HITS, i, std::cout);
2775  }
2776 
2777  edm::LogVerbatim("IsoTrack") << "\n \t trackerLayersWithMeasurement() " << p.trackerLayersWithMeasurement()
2778  << "\n \t pixelLayersWithMeasurement() " << p.pixelLayersWithMeasurement()
2779  << "\n \t stripLayersWithMeasurement() " << p.stripLayersWithMeasurement()
2780  << "\n \t pixelBarrelLayersWithMeasurement() " << p.pixelBarrelLayersWithMeasurement()
2781  << "\n \t pixelEndcapLayersWithMeasurement() " << p.pixelEndcapLayersWithMeasurement()
2782  << "\n \t stripTIBLayersWithMeasurement() " << p.stripTIBLayersWithMeasurement()
2783  << "\n \t stripTIDLayersWithMeasurement() " << p.stripTIDLayersWithMeasurement()
2784  << "\n \t stripTOBLayersWithMeasurement() " << p.stripTOBLayersWithMeasurement()
2785  << "\n \t stripTECLayersWithMeasurement() " << p.stripTECLayersWithMeasurement();
2786  }
2787 }
2788 
2789 //define this as a plug-in
2790 DEFINE_FWK_MODULE(IsolatedTracksNxN);
IsolatedTracksNxN::t_trackPVIdx
std::vector< int > * t_trackPVIdx
Definition: IsolatedTracksNxN.cc:259
reco::Track::outerPosition
const math::XYZPoint & outerPosition() const
position of the outermost hit
Definition: Track.h:62
IsolatedTracksNxN::IsolatedTracksNxN
IsolatedTracksNxN(const edm::ParameterSet &)
Definition: IsolatedTracksNxN.cc:307
spr::caloSimInfo::eRest
double eRest
Definition: CaloSimInfo.h:49
IsolatedTracksNxN::t_esimPdgId
std::vector< double > * t_esimPdgId
Definition: IsolatedTracksNxN.cc:282
IsolatedTracksNxN::t_trackHcalEta
std::vector< double > * t_trackHcalEta
Definition: IsolatedTracksNxN.cc:255
reco::Vertex::trackRef_iterator
std::vector< TrackBaseRef >::const_iterator trackRef_iterator
The iteratator for the vector<TrackRef>
Definition: Vertex.h:38
change_name.diff
diff
Definition: change_name.py:13
EcalSeverityLevelAlgo
Definition: EcalSeverityLevelAlgo.h:33
IsolatedTracksNxN::h_recEtaP_0
TH2F * h_recEtaP_0
Definition: IsolatedTracksNxN.cc:220
reco::HitPattern::MISSING_OUTER_HITS
Definition: HitPattern.h:155
IsolatedTracksNxN
Definition: IsolatedTracksNxN.cc:142
IsolatedTracksNxN::t_esim9x9Rest
std::vector< double > * t_esim9x9Rest
Definition: IsolatedTracksNxN.cc:285
L1GtTriggerMenu.h
CaloJetCollection.h
DDAxes::y
IsolatedTracksNxN::t_e11x11_15Sig
std::vector< double > * t_e11x11_15Sig
Definition: IsolatedTracksNxN.cc:277
Handle.h
IsolatedTracksNxN::t_trackPhi
std::vector< double > * t_trackPhi
Definition: IsolatedTracksNxN.cc:254
IsolatedTracksNxN::t_NLayersCrossed
std::vector< int > * t_NLayersCrossed
Definition: IsolatedTracksNxN.cc:261
IsolatedTracksNxN::t_esim11x11Rest
std::vector< double > * t_esim11x11Rest
Definition: IsolatedTracksNxN.cc:285
IsolatedTracksNxN::t_trackHitInMissTIBTID
std::vector< int > * t_trackHitInMissTIBTID
Definition: IsolatedTracksNxN.cc:263
IsolatedTracksNxN::t_esim9x9CharHad
std::vector< double > * t_esim9x9CharHad
Definition: IsolatedTracksNxN.cc:288
CaloJet.h
IsolatedTracksNxN::t_L1NonIsoEMEta
std::vector< double > * t_L1NonIsoEMEta
Definition: IsolatedTracksNxN.cc:239
L1GtTriggerMenu::gtTriggerMenuName
const std::string & gtTriggerMenuName() const
Definition: L1GtTriggerMenu.h:91
IsolatedTracksNxN::t_trkEcalEne
std::vector< double > * t_trkEcalEne
Definition: IsolatedTracksNxN.cc:282
useL1EventSetup
static const bool useL1EventSetup(true)
IsolatedTracksNxN::t_e7x7_10Sig
std::vector< double > * t_e7x7_10Sig
Definition: IsolatedTracksNxN.cc:276
IsolatedTracksNxN::tok_sevlv_
edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd > tok_sevlv_
Definition: IsolatedTracksNxN.cc:208
electrons_cff.bool
bool
Definition: electrons_cff.py:366
EDAnalyzer.h
IsolatedTracksNxN::t_esim15x15CharHad
std::vector< double > * t_esim15x15CharHad
Definition: IsolatedTracksNxN.cc:288
mps_fire.i
i
Definition: mps_fire.py:428
reco::TrackBase::qoverpError
double qoverpError() const
error on signed transverse curvature
Definition: TrackBase.h:732
FreeTrajectoryState.h
IsolatedTracksNxN::tok_magField_
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > tok_magField_
Definition: IsolatedTracksNxN.cc:206
IsolatedTracksNxN::t_hsim7x7CharHad
std::vector< double > * t_hsim7x7CharHad
Definition: IsolatedTracksNxN.cc:301
CaloNavigator.h
MessageLogger.h
IsolatedTracksNxN::t_L1METEta
std::vector< double > * t_L1METEta
Definition: IsolatedTracksNxN.cc:240
IsolatedTracksNxN::t_esim7x7CharHad
std::vector< double > * t_esim7x7CharHad
Definition: IsolatedTracksNxN.cc:288
IsolatedTracksNxN::h_L1AlgoNames
TH1I * h_L1AlgoNames
Definition: IsolatedTracksNxN.cc:216
L1GtTriggerMaskAlgoTrigRcd.h
funct::false
false
Definition: Factorize.h:29
IsolatedTracksNxN::minTrackP_
const double minTrackP_
Definition: IsolatedTracksNxN.cc:174
IsolatedTracksNxN::h_recPhi_0
TH1F * h_recPhi_0
Definition: IsolatedTracksNxN.cc:219
IsolatedTracksNxN::t_L1MuonPhi
std::vector< double > * t_L1MuonPhi
Definition: IsolatedTracksNxN.cc:237
IsolatedTracksNxN::t_trackChiSqAll
std::vector< double > * t_trackChiSqAll
Definition: IsolatedTracksNxN.cc:252
IsolatedTracksNxN::t_e15x15_10Sig
std::vector< double > * t_e15x15_10Sig
Definition: IsolatedTracksNxN.cc:276
L1GtTriggerMask.h
L1MuonParticleFwd.h
IsolatedTracksNxN::h_maxNearP21x21
TH1F * h_maxNearP21x21[NPBins][NEtaBins]
Definition: IsolatedTracksNxN.cc:214
reco::HitPattern::stripTOBLayersWithMeasurement
int stripTOBLayersWithMeasurement() const
Definition: HitPattern.cc:594
IsolatedTracksNxN::t_e7x7_20Sig
std::vector< double > * t_e7x7_20Sig
Definition: IsolatedTracksNxN.cc:278
IsolatedTracksNxN::t_e9x9_25Sig
std::vector< double > * t_e9x9_25Sig
Definition: IsolatedTracksNxN.cc:279
IsolatedTracksNxN::t_e11x11_25Sig
std::vector< double > * t_e11x11_25Sig
Definition: IsolatedTracksNxN.cc:279
IsolatedTracksNxN::h_PVTracksWt
TH1F * h_PVTracksWt
Definition: IsolatedTracksNxN.cc:217
spr::caloSimInfo::eMatched
double eMatched
Definition: CaloSimInfo.h:45
L1GtUtils::UseEventSetupIn::Event
IsolatedTracksNxN::t_trackDxyPVAll
std::vector< double > * t_trackDxyPVAll
Definition: IsolatedTracksNxN.cc:252
IsolatedTracksNxN::tok_caloHH_
edm::EDGetTokenT< edm::PCaloHitContainer > tok_caloHH_
Definition: IsolatedTracksNxN.cc:201
IsolatedTracksNxN::t_esim15x15Rest
std::vector< double > * t_esim15x15Rest
Definition: IsolatedTracksNxN.cc:285
reco::TrackBase::referencePoint
const Point & referencePoint() const
Reference point on the track.
Definition: TrackBase.h:667
IsolatedTracksNxN::myverbose_
const int myverbose_
Definition: IsolatedTracksNxN.cc:170
reco::HitPattern::stripTIBLayersWithoutMeasurement
int stripTIBLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:636
IsolatedTracksNxN::t_esim9x9Matched
std::vector< double > * t_esim9x9Matched
Definition: IsolatedTracksNxN.cc:284
IsolatedTracksNxN::t_hsim7x7Rest
std::vector< double > * t_hsim7x7Rest
Definition: IsolatedTracksNxN.cc:298
DiDispStaMuonMonitor_cfi.pt
pt
Definition: DiDispStaMuonMonitor_cfi.py:39
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:631
TrackDetectorAssociator.h
IsolatedTracksNxN::t_L1IsoEMPt
std::vector< double > * t_L1IsoEMPt
Definition: IsolatedTracksNxN.cc:238
EcalSeverityLevelAlgoRcd.h
L1GlobalTriggerReadoutSetupFwd.h
IsolatedTracksNxN::h_recP_2
TH1F * h_recP_2
Definition: IsolatedTracksNxN.cc:223
deltaPhi.h
IsolatedTracksNxN::h_nTracks
TH1F * h_nTracks
Definition: IsolatedTracksNxN.cc:218
edm::EDGetTokenT< l1extra::L1JetParticleCollection >
EBDetId
Definition: EBDetId.h:17
TrackBase.h
IsolatedTracksNxN::t_maxNearP21x21
std::vector< double > * t_maxNearP21x21
Definition: IsolatedTracksNxN.cc:272
IsolatedTracksNxN::t_trackEcalPhi
std::vector< double > * t_trackEcalPhi
Definition: IsolatedTracksNxN.cc:255
IsolatedTracksNxN::t_esim7x7
std::vector< double > * t_esim7x7
Definition: IsolatedTracksNxN.cc:283
IsolatedTracksNxN::t_trackHitOutMeasTIB
std::vector< int > * t_trackHitOutMeasTIB
Definition: IsolatedTracksNxN.cc:268
L1GtTriggerMenuRcd.h
gpuVertexFinder::iv
int32_t *__restrict__ iv
Definition: gpuClusterTracksDBSCAN.h:42
IsolatedTracksNxN::t_esim9x9Photon
std::vector< double > * t_esim9x9Photon
Definition: IsolatedTracksNxN.cc:286
IsolatedTracksNxN::t_L1MuonPt
std::vector< double > * t_L1MuonPt
Definition: IsolatedTracksNxN.cc:237
gather_cfg.cout
cout
Definition: gather_cfg.py:144
IsolatedTracksNxN::writeAllTracks_
const bool writeAllTracks_
Definition: IsolatedTracksNxN.cc:169
HcalTopology
Definition: HcalTopology.h:26
EBDetId.h
EEDetId.h
IsolatedTracksNxN::tMinE_
const double tMinE_
Definition: IsolatedTracksNxN.cc:176
IsolatedTracksNxN::t_trackNOuterHits
std::vector< int > * t_trackNOuterHits
Definition: IsolatedTracksNxN.cc:261
IsolatedTracksNxN::m_l1GtUtils
std::unique_ptr< L1GtUtils > m_l1GtUtils
Definition: IsolatedTracksNxN.cc:161
IsolatedTracksNxN::t_RunNo
int t_RunNo
Definition: IsolatedTracksNxN.cc:229
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89301
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
IsolatedTracksNxN::t_trackDxy
std::vector< double > * t_trackDxy
Definition: IsolatedTracksNxN.cc:256
spr::eHCALSimInfo
std::map< std::string, double > eHCALSimInfo(const edm::Event &, const HcalTopology *topology, const DetId &det, const CaloGeometry *geo, edm::Handle< T > &hits, edm::Handle< edm::SimTrackContainer > &SimTk, edm::Handle< edm::SimVertexContainer > &SimVtx, const reco::Track *pTrack, TrackerHitAssociator &associate, int ieta, int iphi, double timeCut=150, bool includeHO=false, bool debug=false)
IsolatedTracksNxN::h_recEtaPt_0
TH2F * h_recEtaPt_0
Definition: IsolatedTracksNxN.cc:220
IsolatedTracksNxN::genPartPBins
double genPartPBins[NPBins+1]
Definition: IsolatedTracksNxN.cc:212
IsolatedTracksNxN::t_hsim5x5CharHad
std::vector< double > * t_hsim5x5CharHad
Definition: IsolatedTracksNxN.cc:301
IsolatedTracksNxN::t_L1CenJetPt
std::vector< double > * t_L1CenJetPt
Definition: IsolatedTracksNxN.cc:234
IsolatedTracksNxN::t_trackHitOutMissTEC
std::vector< int > * t_trackHitOutMissTEC
Definition: IsolatedTracksNxN.cc:265
IsolatedTracksNxN::t_esim11x11Matched
std::vector< double > * t_esim11x11Matched
Definition: IsolatedTracksNxN.cc:284
IsolatedTracksNxN::t_trkHcalEne
std::vector< double > * t_trkHcalEne
Definition: IsolatedTracksNxN.cc:295
IsolatedTracksNxN::t_trackDxyAll
std::vector< double > * t_trackDxyAll
Definition: IsolatedTracksNxN.cc:252
IsolatedTracksNxN::t_PVTracksSumPt
std::vector< double > * t_PVTracksSumPt
Definition: IsolatedTracksNxN.cc:242
IsolatedTracksNxN::t_trackHcalPhi
std::vector< double > * t_trackHcalPhi
Definition: IsolatedTracksNxN.cc:255
IsolatedTracksNxN::t_trackPtAll
std::vector< double > * t_trackPtAll
Definition: IsolatedTracksNxN.cc:251
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:150
L1GtTriggerMaskTechTrigRcd.h
IsolatedTracksNxN::h_recEtaPt_1
TH2F * h_recEtaPt_1
Definition: IsolatedTracksNxN.cc:222
edm::SortedCollection
Definition: SortedCollection.h:49
IsolatedTracksNxN::t_trackDxyBS
std::vector< double > * t_trackDxyBS
Definition: IsolatedTracksNxN.cc:256
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:47
reco::TrackBase::qualityName
static std::string qualityName(TrackQuality)
Definition: TrackBase.h:572
IsolatedTracksNxN::t_hsim3x3Photon
std::vector< double > * t_hsim3x3Photon
Definition: IsolatedTracksNxN.cc:299
TrackerHitAssociator.h
IsolatedTracksNxN::clearTreeVectors
void clearTreeVectors()
Definition: IsolatedTracksNxN.cc:2057
IsolatedTracksNxN::t_infoHcal
std::vector< int > * t_infoHcal
Definition: IsolatedTracksNxN.cc:293
IsolatedTracksNxN::nbad_
int nbad_
Definition: IsolatedTracksNxN.cc:178
IsolatedTracksNxN::t_trackHitsTEC
std::vector< int > * t_trackHitsTEC
Definition: IsolatedTracksNxN.cc:262
IsolatedTracksNxN::t_e9x9_30Sig
std::vector< double > * t_e9x9_30Sig
Definition: IsolatedTracksNxN.cc:280
singleTopDQM_cfi.jets
jets
Definition: singleTopDQM_cfi.py:42
DDAxes::x
IsolatedTracksNxN::tok_geom_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord > tok_geom_
Definition: IsolatedTracksNxN.cc:203
IsolatedTracksNxN::t_trackOutPosOutHitDr
std::vector< double > * t_trackOutPosOutHitDr
Definition: IsolatedTracksNxN.cc:269
TFileDirectory
Definition: TFileDirectory.h:24
IsolatedTracksNxN::t_maxNearHcalP5x5
std::vector< double > * t_maxNearHcalP5x5
Definition: IsolatedTracksNxN.cc:290
IsolatedTracksNxN::t_trackHitInMissTEC
std::vector< int > * t_trackHitInMissTEC
Definition: IsolatedTracksNxN.cc:263
edm::one::EDAnalyzer
Definition: EDAnalyzer.h:30
IsolatedTracksNxN::t_trackHitInMissTOB
std::vector< int > * t_trackHitInMissTOB
Definition: IsolatedTracksNxN.cc:263
TrackerHitAssociator::Config
Definition: TrackerHitAssociator.h:57
IsolatedTracksNxN::t_esim7x7Matched
std::vector< double > * t_esim7x7Matched
Definition: IsolatedTracksNxN.cc:284
IsolatedTracksNxN::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: IsolatedTracksNxN.cc:766
EcalCondObjectContainer
Definition: EcalCondObjectContainer.h:13
BeamSpotPI::nTracks
Definition: BeamSpotPayloadInspectorHelper.h:42
IsolatedTracksNxN::t_trackChiSq
std::vector< double > * t_trackChiSq
Definition: IsolatedTracksNxN.cc:258
edm::Handle< reco::TrackCollection >
reco::TrackBase::numberOfValidHits
unsigned short numberOfValidHits() const
number of valid hits found
Definition: TrackBase.h:798
IsolatedTracksNxN::t_e11x11
std::vector< double > * t_e11x11
Definition: IsolatedTracksNxN.cc:275
IsolatedTracksNxN::t_L1PreScale
std::vector< int > * t_L1PreScale
Definition: IsolatedTracksNxN.cc:231
CaloTopology
Definition: CaloTopology.h:19
IsolatedTracksNxN::t_PVy
std::vector< double > * t_PVy
Definition: IsolatedTracksNxN.cc:242
IsolatedTracksNxN::m_triggerMaskAlgoTrig
std::vector< unsigned int > m_triggerMaskAlgoTrig
Definition: IsolatedTracksNxN.cc:167
EcalBarrel
Definition: EcalSubdetector.h:10
EcalRecHitCollections.h
spr::caloSimInfo::eGamma
double eGamma
Definition: CaloSimInfo.h:46
IsolatedTracksNxN::t_e15x15_25Sig
std::vector< double > * t_e15x15_25Sig
Definition: IsolatedTracksNxN.cc:279
reco::HitPattern::stripTOBLayersWithoutMeasurement
int stripTOBLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:656
IsolatedTracksNxN::t_trackHitOutMissTIB
std::vector< int > * t_trackHitOutMissTIB
Definition: IsolatedTracksNxN.cc:265
IsolatedTracksNxN::t_e9x9_10Sig
std::vector< double > * t_e9x9_10Sig
Definition: IsolatedTracksNxN.cc:276
IsolatedTracksNxN::h_recPt_1
TH1F * h_recPt_1
Definition: IsolatedTracksNxN.cc:221
IsolatedTracksNxN::t_esim11x11
std::vector< double > * t_esim11x11
Definition: IsolatedTracksNxN.cc:283
IsolatedTracksNxN::t_trackP
std::vector< double > * t_trackP
Definition: IsolatedTracksNxN.cc:254
IsolatedTracksNxN::t_L1METPhi
std::vector< double > * t_L1METPhi
Definition: IsolatedTracksNxN.cc:240
ndof
Definition: HIMultiTrackSelector.h:49
spr::caloSimInfo::eTotal
double eTotal
Definition: CaloSimInfo.h:50
IsolatedTracksNxN::NPBins
static constexpr size_t NPBins
Definition: IsolatedTracksNxN.cc:210
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:637
DetId
Definition: DetId.h:17
IsolatedTracksNxN::t_e11x11_20Sig
std::vector< double > * t_e11x11_20Sig
Definition: IsolatedTracksNxN.cc:278
IsolatedTracksNxN::tok_L1extNoIsoEm_
edm::EDGetTokenT< l1extra::L1EmParticleCollection > tok_L1extNoIsoEm_
Definition: IsolatedTracksNxN.cc:186
MakerMacros.h
IsolatedTracksNxN::t_L1TauJetPt
std::vector< double > * t_L1TauJetPt
Definition: IsolatedTracksNxN.cc:236
CaloGeometry
Definition: CaloGeometry.h:21
reco::HitPattern
Definition: HitPattern.h:147
FindCaloHit.h
IsolatedTracksNxN::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: IsolatedTracksNxN.cc:735
IsolatedTracksNxN::t_L1FwdJetEta
std::vector< double > * t_L1FwdJetEta
Definition: IsolatedTracksNxN.cc:235
IsolatedTracksNxN::t_trackDxyPV
std::vector< double > * t_trackDxyPV
Definition: IsolatedTracksNxN.cc:257
Track.h
TrackFwd.h
IsolatedTracksNxN::h_recPhi_2
TH1F * h_recPhi_2
Definition: IsolatedTracksNxN.cc:223
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
IsolatedTracksNxN::t_h3x3
std::vector< double > * t_h3x3
Definition: IsolatedTracksNxN.cc:291
BeamSpot.h
CItAlgo
AlgorithmMap::const_iterator CItAlgo
iterators through map containing the algorithms
Definition: L1GtTriggerMenuFwd.h:38
IsolatedTracksNxN::t_nTracks
int t_nTracks
Definition: IsolatedTracksNxN.cc:228
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)
IsolatedTracksNxN::t_esim7x7Photon
std::vector< double > * t_esim7x7Photon
Definition: IsolatedTracksNxN.cc:286
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
IsolatedTracksNxN::t_hsim5x5Photon
std::vector< double > * t_hsim5x5Photon
Definition: IsolatedTracksNxN.cc:299
reco::HitPattern::trackerLayersWithMeasurement
int trackerLayersWithMeasurement() const
Definition: HitPattern.cc:510
IsolatedTracksNxN::t_esim9x9NeutHad
std::vector< double > * t_esim9x9NeutHad
Definition: IsolatedTracksNxN.cc:287
IsolatedTracksNxN::t_trackHitInMeasTIB
std::vector< int > * t_trackHitInMeasTIB
Definition: IsolatedTracksNxN.cc:267
IsolatedTracksNxN::h_recEtaP_2
TH2F * h_recEtaP_2
Definition: IsolatedTracksNxN.cc:224
IsolatedTracksNxN::t_trackPt
std::vector< double > * t_trackPt
Definition: IsolatedTracksNxN.cc:254
spr::chargeIsolationHcal
double chargeIsolationHcal(unsigned int trkIndex, std::vector< spr::propagatedTrackID > &vdetIds, const HcalTopology *topology, int ieta, int iphi, bool debug=false)
Definition: ChargeIsolation.cc:114
IsolatedTracksNxN::printTrkHitPattern_
const bool printTrkHitPattern_
Definition: IsolatedTracksNxN.cc:173
IsolatedTracksNxN::t_e11x11_10Sig
std::vector< double > * t_e11x11_10Sig
Definition: IsolatedTracksNxN.cc:276
IsolatedTracksNxN::t_L1IsoEMEta
std::vector< double > * t_L1IsoEMEta
Definition: IsolatedTracksNxN.cc:238
Service.h
IsolatedTracksNxN::NEtaBins
static constexpr size_t NEtaBins
Definition: IsolatedTracksNxN.cc:211
PVValHelper::eta
Definition: PVValidationHelpers.h:70
IsolatedTracksNxN::tok_ecalChStatus_
edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > tok_ecalChStatus_
Definition: IsolatedTracksNxN.cc:207
IsolatedTracksNxN::t_maxNearHcalP3x3
std::vector< double > * t_maxNearHcalP3x3
Definition: IsolatedTracksNxN.cc:290
ChargeIsolation.h
IsolatedTracksNxN::t_PVNTracksHPWt
std::vector< int > * t_PVNTracksHPWt
Definition: IsolatedTracksNxN.cc:245
IsolatedTracksNxN::h_maxNearP25x25
TH1F * h_maxNearP25x25[NPBins][NEtaBins]
Definition: IsolatedTracksNxN.cc:214
SimVertex.h
IsolatedTracksNxN::~IsolatedTracksNxN
~IsolatedTracksNxN() override
Definition: IsolatedTracksNxN.cc:558
eHCALMatrix.h
DDAxes::z
reco::TrackBase::numberOfLostHits
unsigned short numberOfLostHits() const
number of cases where track crossed a layer without getting a hit.
Definition: TrackBase.h:801
reco::Track
Definition: Track.h:27
AlgorithmMap
std::map< std::string, L1GtAlgorithm > AlgorithmMap
map containing the algorithms
Definition: L1GtTriggerMenuFwd.h:31
IdealMagneticFieldRecord.h
IsolatedTracksNxN::t_trackPdgIdAll
std::vector< double > * t_trackPdgIdAll
Definition: IsolatedTracksNxN.cc:250
IsolatedTracksNxN::t_hsim3x3NeutHad
std::vector< double > * t_hsim3x3NeutHad
Definition: IsolatedTracksNxN.cc:300
IsolatedTracksNxN::t_L1NonIsoEMPhi
std::vector< double > * t_L1NonIsoEMPhi
Definition: IsolatedTracksNxN.cc:239
IsolatedTracksNxN::t_PVz
std::vector< double > * t_PVz
Definition: IsolatedTracksNxN.cc:242
IsolatedTracksNxN::t_PVNTracksWt
std::vector< int > * t_PVNTracksWt
Definition: IsolatedTracksNxN.cc:244
IsolatedTracksNxN::t_L1NonIsoEMPt
std::vector< double > * t_L1NonIsoEMPt
Definition: IsolatedTracksNxN.cc:239
IsolatedTracksNxN::t_L1Decision
int t_L1Decision[128]
Definition: IsolatedTracksNxN.cc:232
IsolatedTracksNxN::t_Lumi
int t_Lumi
Definition: IsolatedTracksNxN.cc:229
HLT_FULL_cff.pt1
pt1
Definition: HLT_FULL_cff.py:9893
IsolatedTracksNxN::t_jetEta
std::vector< double > * t_jetEta
Definition: IsolatedTracksNxN.cc:247
IsolatedTracksNxN::tok_caloEB_
edm::EDGetTokenT< edm::PCaloHitContainer > tok_caloEB_
Definition: IsolatedTracksNxN.cc:199
reco::TrackBase::charge
int charge() const
track electric charge
Definition: TrackBase.h:596
reco::TrackBase::dz
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:622
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:59
IsolatedTracksNxN::t_esim15x15NeutHad
std::vector< double > * t_esim15x15NeutHad
Definition: IsolatedTracksNxN.cc:287
IsolatedTracksNxN::t_L1AlgoNames
std::vector< std::string > * t_L1AlgoNames
Definition: IsolatedTracksNxN.cc:230
sistrip::SpyUtilities::isValid
const bool isValid(const Frame &aFrame, const FrameQuality &aQuality, const uint16_t aExpectedPos)
Definition: SiStripSpyUtilities.cc:124
eECALMatrix.h
IsolatedTracksNxN::t_esim15x15
std::vector< double > * t_esim15x15
Definition: IsolatedTracksNxN.cc:283
IsolatedTracksNxN::t_h7x7Sig
std::vector< double > * t_h7x7Sig
Definition: IsolatedTracksNxN.cc:292
IsolatedTracksNxN::bookHistograms
void bookHistograms()
Definition: IsolatedTracksNxN.cc:2260
edm::RefToBase::get
value_type const * get() const
Definition: RefToBase.h:209
EcalSeverityLevelAlgo.h
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:649
IsolatedTracksNxN::t_trackPhiAll
std::vector< double > * t_trackPhiAll
Definition: IsolatedTracksNxN.cc:250
first
auto first
Definition: CAHitNtupletGeneratorKernelsImpl.h:125
useL1GtTriggerMenuLite
static const bool useL1GtTriggerMenuLite(true)
IsolatedTracksNxN::t_trackEtaAll
std::vector< double > * t_trackEtaAll
Definition: IsolatedTracksNxN.cc:250
EEDetId
Definition: EEDetId.h:14
IsolatedTracksNxN::tok_htmap_
edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord > tok_htmap_
Definition: IsolatedTracksNxN.cc:209
edm::ConfigurationDescriptions
Definition: ConfigurationDescriptions.h:28
IsolatedTracksNxN::h_recEtaP_1
TH2F * h_recEtaP_1
Definition: IsolatedTracksNxN.cc:222
L1EmParticleFwd.h
IsolatedTracksNxN::t_L1FwdJetPhi
std::vector< double > * t_L1FwdJetPhi
Definition: IsolatedTracksNxN.cc:235
TFileService.h
EcalTrigTowerConstituentsMap
Definition: EcalTrigTowerConstituentsMap.h:19
HLT_FULL_cff.eta1
eta1
Definition: HLT_FULL_cff.py:9550
IsolatedTracksNxN::t_hsim7x7
std::vector< double > * t_hsim7x7
Definition: IsolatedTracksNxN.cc:296
reco::HitPattern::stripTIBLayersWithMeasurement
int stripTIBLayersWithMeasurement() const
Definition: HitPattern.cc:574
IsolatedTracksNxN::tok_simVtx_
edm::EDGetTokenT< edm::SimVertexContainer > tok_simVtx_
Definition: IsolatedTracksNxN.cc:198
IsolatedTracksNxN::t_PVNTracks
std::vector< int > * t_PVNTracks
Definition: IsolatedTracksNxN.cc:244
TFileService::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
IsolatedTracksNxN::tMaxH_
const double tMaxH_
Definition: IsolatedTracksNxN.cc:176
IsolatedTracksNxN::t_nTrksJetVtx
std::vector< double > * t_nTrksJetVtx
Definition: IsolatedTracksNxN.cc:248
IsolatedTracksNxN::t_PVndof
std::vector< int > * t_PVndof
Definition: IsolatedTracksNxN.cc:244
LEDCalibrationChannels.ieta
ieta
Definition: LEDCalibrationChannels.py:63
IsolatedTracksNxN::t_PVx
std::vector< double > * t_PVx
Definition: IsolatedTracksNxN.cc:242
CaloSubdetectorGeometry.h
L1JetParticleFwd.h
edm::ParameterSet
Definition: ParameterSet.h:47
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
IsolatedTracksNxN::nL1BitsMax
static constexpr size_t nL1BitsMax
Definition: IsolatedTracksNxN.cc:162
IsolatedTracksNxN::t_L1CenJetEta
std::vector< double > * t_L1CenJetEta
Definition: IsolatedTracksNxN.cc:234
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
IsolatedTracksNxN::tok_recVtx_
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
Definition: IsolatedTracksNxN.cc:192
IsolatedTracksNxN::t_maxNearP31x31
std::vector< double > * t_maxNearP31x31
Definition: IsolatedTracksNxN.cc:271
L1GtTriggerMenu
Definition: L1GtTriggerMenu.h:48
Event.h
L1GlobalTriggerReadoutSetup.h
IsolatedTracksNxN::tok_L1extMu_
edm::EDGetTokenT< l1extra::L1MuonParticleCollection > tok_L1extMu_
Definition: IsolatedTracksNxN.cc:184
IsolatedTracksNxN::t_e15x15_15Sig
std::vector< double > * t_e15x15_15Sig
Definition: IsolatedTracksNxN.cc:277
deltaR.h
spr::caloSimInfo::eNeutralHad
double eNeutralHad
Definition: CaloSimInfo.h:47
IsolatedTracksNxN::tok_topo_
edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tok_topo_
Definition: IsolatedTracksNxN.cc:205
IsolatedTracksNxN::t_hsim5x5Rest
std::vector< double > * t_hsim5x5Rest
Definition: IsolatedTracksNxN.cc:298
IsolatedTracksNxN::tok_caloEE_
edm::EDGetTokenT< edm::PCaloHitContainer > tok_caloEE_
Definition: IsolatedTracksNxN.cc:200
IsolatedTracksNxN::t_hsim5x5
std::vector< double > * t_hsim5x5
Definition: IsolatedTracksNxN.cc:296
IsolatedTracksNxN::t_e15x15_30Sig
std::vector< double > * t_e15x15_30Sig
Definition: IsolatedTracksNxN.cc:280
HcalDetId.h
IsolatedTracksNxN::t_trackL
std::vector< double > * t_trackL
Definition: IsolatedTracksNxN.cc:269
IsolatedTracksNxN::t_e9x9_20Sig
std::vector< double > * t_e9x9_20Sig
Definition: IsolatedTracksNxN.cc:278
IsolatedTracksNxN::tMaxE_
const double tMaxE_
Definition: IsolatedTracksNxN.cc:176
IsolatedTracksNxN::t_hsim3x3Rest
std::vector< double > * t_hsim3x3Rest
Definition: IsolatedTracksNxN.cc:298
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:652
MatchingSimTrack.h
IsolatedTracksNxN::t_L1IsoEMPhi
std::vector< double > * t_L1IsoEMPhi
Definition: IsolatedTracksNxN.cc:238
IsolatedTracksNxN::h_recEta_0
TH1F * h_recEta_0
Definition: IsolatedTracksNxN.cc:219
HcalDetId
Definition: HcalDetId.h:12
edm::Service< TFileService >
createfilelist.int
int
Definition: createfilelist.py:10
IsolatedTracksNxN::t_L1METPt
std::vector< double > * t_L1METPt
Definition: IsolatedTracksNxN.cc:240
iEvent
int iEvent
Definition: GenABIO.cc:224
IsolatedTracksNxN::t_e7x7_15Sig
std::vector< double > * t_e7x7_15Sig
Definition: IsolatedTracksNxN.cc:277
IsolatedTracksNxN::t_e9x9_15Sig
std::vector< double > * t_e9x9_15Sig
Definition: IsolatedTracksNxN.cc:277
IsolatedTracksNxN::t_esim11x11Photon
std::vector< double > * t_esim11x11Photon
Definition: IsolatedTracksNxN.cc:286
reco::TrackBase::d0
double d0() const
dxy parameter in perigee convention (d0 = -dxy)
Definition: TrackBase.h:611
IsolatedTracksNxN::t_trackHitInMeasTOB
std::vector< int > * t_trackHitInMeasTOB
Definition: IsolatedTracksNxN.cc:267
reco::HitPattern::stripTIDLayersWithoutMeasurement
int stripTIDLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:646
p1
double p1[4]
Definition: TauolaWrapper.h:89
IsolatedTracksNxN::t_PVTracksSumPtHP
std::vector< double > * t_PVTracksSumPtHP
Definition: IsolatedTracksNxN.cc:243
IsolatedTracksNxN::t_trackDzAll
std::vector< double > * t_trackDzAll
Definition: IsolatedTracksNxN.cc:252
MagneticField.h
IsolatedTracksNxN::t_trackHitOutMeasTID
std::vector< int > * t_trackHitOutMeasTID
Definition: IsolatedTracksNxN.cc:268
edm::EventSetup
Definition: EventSetup.h:58
IsolatedTracksNxN::t_PVTracksSumPtHPWt
std::vector< double > * t_PVTracksSumPtHPWt
Definition: IsolatedTracksNxN.cc:243
IsolatedTracksNxN::t_nTrksJetCalo
std::vector< double > * t_nTrksJetCalo
Definition: IsolatedTracksNxN.cc:248
TrackAssociatorRecord.h
L1EmParticle.h
IsolatedTracksNxN::t_trackEta
std::vector< double > * t_trackEta
Definition: IsolatedTracksNxN.cc:254
reco::TrackBase::normalizedChi2
double normalizedChi2() const
chi-squared divided by n.d.o.f. (or chi-squared * 1e6 if n.d.o.f. is zero)
Definition: TrackBase.h:593
spr::eECALSimInfo
void eECALSimInfo(const edm::Event &, const DetId &det, const CaloGeometry *geo, const CaloTopology *caloTopology, edm::Handle< T > &hitsEB, edm::Handle< T > &hitsEE, edm::Handle< edm::SimTrackContainer > &SimTk, edm::Handle< edm::SimVertexContainer > &SimVtx, const reco::Track *pTrack, TrackerHitAssociator &associate, int ieta, int iphi, caloSimInfo &info, double timeCut=150, bool debug=false)
IsolatedTracksNxN::t_trackHitInMeasTID
std::vector< int > * t_trackHitInMeasTID
Definition: IsolatedTracksNxN.cc:267
HitPattern.h
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >
IsolatedTracksNxN::t_hsim7x7Matched
std::vector< double > * t_hsim7x7Matched
Definition: IsolatedTracksNxN.cc:297
reco::HitPattern::TRACK_HITS
Definition: HitPattern.h:155
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:127
IsolatedTracksNxN::t_PVisValid
std::vector< int > * t_PVisValid
Definition: IsolatedTracksNxN.cc:244
IsolatedTracksNxN::t_hsim7x7NeutHad
std::vector< double > * t_hsim7x7NeutHad
Definition: IsolatedTracksNxN.cc:300
IsolatedTracksNxN::h_recP_1
TH1F * h_recP_1
Definition: IsolatedTracksNxN.cc:221
CaloTopology.h
photonIsolationHIProducer_cfi.hbhe
hbhe
Definition: photonIsolationHIProducer_cfi.py:8
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
IsolatedTracksNxN::nEventProc_
int nEventProc_
Definition: IsolatedTracksNxN.cc:178
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
VertexFwd.h
DDAxes::phi
IsolatedTracksNxN::t_trackHitInMissTIB
std::vector< int > * t_trackHitInMissTIB
Definition: IsolatedTracksNxN.cc:263
CaloSubdetectorTopology.h
IsolatedTracksNxN::l1AlgoMap_
std::map< std::pair< unsigned int, std::string >, int > l1AlgoMap_
Definition: IsolatedTracksNxN.cc:166
IsolatedTracksNxN::t_h3x3Sig
std::vector< double > * t_h3x3Sig
Definition: IsolatedTracksNxN.cc:292
IsolatedTracksNxN::t_trackHitOutMissTID
std::vector< int > * t_trackHitOutMissTID
Definition: IsolatedTracksNxN.cc:265
IsolatedTracksNxN::t_trackHitOutMissTOBTEC
std::vector< int > * t_trackHitOutMissTOBTEC
Definition: IsolatedTracksNxN.cc:265
reco::TrackBase::hitPattern
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:504
IsolatedTracksNxN::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: IsolatedTracksNxN.cc:189
CaloTopologyRecord.h
IsolatedTracksNxN::t_h5x5Sig
std::vector< double > * t_h5x5Sig
Definition: IsolatedTracksNxN.cc:292
IsolatedTracksNxN::tok_EE_
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: IsolatedTracksNxN.cc:196
HcalTopology.h
EcalTrigTowerConstituentsMap.h
submitPVResolutionJobs.desc
string desc
Definition: submitPVResolutionJobs.py:251
IsolatedTracksNxN::t_esim15x15Photon
std::vector< double > * t_esim15x15Photon
Definition: IsolatedTracksNxN.cc:286
CaloSimInfo.h
IsolatedTracksNxN::endJob
void endJob() override
Definition: IsolatedTracksNxN.cc:2041
IsolatedTracksNxN::t_hsim5x5Matched
std::vector< double > * t_hsim5x5Matched
Definition: IsolatedTracksNxN.cc:297
IsolatedTracksNxN::h_recP_0
TH1F * h_recP_0
Definition: IsolatedTracksNxN.cc:219
reco::HitPattern::stripTIDLayersWithMeasurement
int stripTIDLayersWithMeasurement() const
Definition: HitPattern.cc:584
Ref.h
IsolatedTracksNxN::genPartEtaBins
double genPartEtaBins[NEtaBins+1]
Definition: IsolatedTracksNxN.cc:212
L1GlobalTriggerObjectMapFwd.h
IsolatedTracksNxN::t_L1CenJetPhi
std::vector< double > * t_L1CenJetPhi
Definition: IsolatedTracksNxN.cc:234
IsolatedTracksNxN::t_trackEcalEta
std::vector< double > * t_trackEcalEta
Definition: IsolatedTracksNxN.cc:255
L1MuonParticle.h
IsolatedTracksNxN::L1TriggerAlgoInfo_
const bool L1TriggerAlgoInfo_
Definition: IsolatedTracksNxN.cc:175
JetExtendedAssociation.h
Calorimetry_cff.bField
bField
Definition: Calorimetry_cff.py:284
IsolatedTracksNxN::h_maxNearP15x15
TH1F * h_maxNearP15x15[NPBins][NEtaBins]
Definition: IsolatedTracksNxN.cc:214
Vertex.h
spr::propagateTrackerEnd
std::pair< math::XYZPoint, double > propagateTrackerEnd(const reco::Track *, const MagneticField *, bool debug=false)
Definition: CaloPropagateTrack.cc:782
DetId.h
Frameworkfwd.h
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
CaloPropagateTrack.h
IsolatedTracksNxN::tok_simTk_
edm::EDGetTokenT< edm::SimTrackContainer > tok_simTk_
Definition: IsolatedTracksNxN.cc:197
spr::matchedSimTrack
edm::SimTrackContainer::const_iterator matchedSimTrack(const edm::Event &iEvent, edm::Handle< edm::SimTrackContainer > &SimTk, edm::Handle< edm::SimVertexContainer > &SimVtx, const reco::Track *pTrack, TrackerHitAssociator &associate, bool debug=false)
Definition: MatchingSimTrack.cc:12
IsolatedTracksNxN::t_trackHitOutMeasTOB
std::vector< int > * t_trackHitOutMeasTOB
Definition: IsolatedTracksNxN.cc:268
L1GtUtils.h
IsolatedTracksNxN::t_hsim7x7Photon
std::vector< double > * t_hsim7x7Photon
Definition: IsolatedTracksNxN.cc:299
IsolatedTracksNxN::tok_caloTopology_
edm::ESGetToken< CaloTopology, CaloTopologyRecord > tok_caloTopology_
Definition: IsolatedTracksNxN.cc:204
CaloGeometry.h
TFileService::kSharedResource
static const std::string kSharedResource
Definition: TFileService.h:76
reco::HitPattern::MISSING_INNER_HITS
Definition: HitPattern.h:155
IsolatedTracksNxN::tok_genTrack_
edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
Definition: IsolatedTracksNxN.cc:191
IsolatedTracksNxN::h_recEta_1
TH1F * h_recEta_1
Definition: IsolatedTracksNxN.cc:221
IsolatedTracksNxN::t_trackHitOutMeasTEC
std::vector< int > * t_trackHitOutMeasTEC
Definition: IsolatedTracksNxN.cc:268
IsolatedTracksNxN::h_recPhi_1
TH1F * h_recPhi_1
Definition: IsolatedTracksNxN.cc:221
IsolatedTracksNxN::h_recEtaPt_2
TH2F * h_recEtaPt_2
Definition: IsolatedTracksNxN.cc:224
IsolatedTracksNxN::h_recEta_2
TH1F * h_recEta_2
Definition: IsolatedTracksNxN.cc:223
IsolatedTracksNxN::doMC_
const bool doMC_
Definition: IsolatedTracksNxN.cc:169
IsolatedTracksNxN::tok_L1extFwdJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extFwdJet_
Definition: IsolatedTracksNxN.cc:182
IsolatedTracksNxN::t_EvtNo
int t_EvtNo
Definition: IsolatedTracksNxN.cc:229
Point3D.h
spr::eHCALmatrixCell
std::vector< std::pair< DetId, double > > eHCALmatrixCell(const HcalTopology *topology, const DetId &det, edm::Handle< T > &hits, int ieta, int iphi, bool includeHO=false, double hbThr=-100, double heThr=-100, double hfThr=-100, double hoThr=-100, bool debug=false)
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
EventSetup.h
IsolatedTracksNxN::t_trackHitInMissTID
std::vector< int > * t_trackHitInMissTID
Definition: IsolatedTracksNxN.cc:263
IsolatedTracksNxN::trackerHitAssociatorConfig_
TrackerHitAssociator::Config trackerHitAssociatorConfig_
Definition: IsolatedTracksNxN.cc:163
spr::caloSimInfo::eChargedHad
double eChargedHad
Definition: CaloSimInfo.h:48
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
edm::RefToBase< reco::Track >
IsolatedTracksNxN::tMinH_
const double tMinH_
Definition: IsolatedTracksNxN.cc:176
IsolatedTracksNxN::t_e11x11_30Sig
std::vector< double > * t_e11x11_30Sig
Definition: IsolatedTracksNxN.cc:280
HcalRecHitCollections.h
IsolatedTracksNxN::t_L1TauJetEta
std::vector< double > * t_L1TauJetEta
Definition: IsolatedTracksNxN.cc:236
IsolatedTracksNxN::t_esim7x7NeutHad
std::vector< double > * t_esim7x7NeutHad
Definition: IsolatedTracksNxN.cc:287
PCaloHitContainer.h
IsolatedTracksNxN::tree_
TTree * tree_
Definition: IsolatedTracksNxN.cc:226
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)
IsolatedTracksNxN::t_trackDz
std::vector< double > * t_trackDz
Definition: IsolatedTracksNxN.cc:256
IsolatedTracksNxN::t_trackDzBS
std::vector< double > * t_trackDzBS
Definition: IsolatedTracksNxN.cc:256
IsolatedTracksNxN::t_e9x9
std::vector< double > * t_e9x9
Definition: IsolatedTracksNxN.cc:275
IsolatedTracksNxN::debugTrks_
const int debugTrks_
Definition: IsolatedTracksNxN.cc:172
IsolatedTracksNxN::pvTracksPtMin_
const double pvTracksPtMin_
Definition: IsolatedTracksNxN.cc:171
spr::caloSimInfo
Definition: CaloSimInfo.h:42
IsolatedTracksNxN::t_trackHitOutMissTOB
std::vector< int > * t_trackHitOutMissTOB
Definition: IsolatedTracksNxN.cc:265
IsolatedTracksNxN::t_esim11x11CharHad
std::vector< double > * t_esim11x11CharHad
Definition: IsolatedTracksNxN.cc:288
IsolatedTracksNxN::t_hsim3x3Matched
std::vector< double > * t_hsim3x3Matched
Definition: IsolatedTracksNxN.cc:297
IsolatedTracksNxN::printTrack
void printTrack(const reco::Track *pTrack)
Definition: IsolatedTracksNxN.cc:2747
IsolatedTracksNxN::t_trackDzPVAll
std::vector< double > * t_trackDzPVAll
Definition: IsolatedTracksNxN.cc:252
IsolatedTracksNxN::t_e15x15
std::vector< double > * t_e15x15
Definition: IsolatedTracksNxN.cc:275
IsolatedTracksNxN::t_h7x7
std::vector< double > * t_h7x7
Definition: IsolatedTracksNxN.cc:291
IsolatedTracksNxN::t_ecalSpike11x11
std::vector< int > * t_ecalSpike11x11
Definition: IsolatedTracksNxN.cc:273
Candidate.h
IsolatedTracksNxN::tok_EB_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: IsolatedTracksNxN.cc:195
IsolatedTracksNxN::t_trackPAll
std::vector< double > * t_trackPAll
Definition: IsolatedTracksNxN.cc:250
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
IsolatedTracksNxN::h_recPt_0
TH1F * h_recPt_0
Definition: IsolatedTracksNxN.cc:219
L1GtTriggerMenuFwd.h
IsolatedTracksNxN::t_e7x7
std::vector< double > * t_e7x7
Definition: IsolatedTracksNxN.cc:275
IsolatedTracksNxN::debugL1Info_
const bool debugL1Info_
Definition: IsolatedTracksNxN.cc:175
SimTrack.h
L1GlobalTriggerObjectMapRecord.h
ParameterSet.h
reco::HitPattern::stripTECLayersWithMeasurement
int stripTECLayersWithMeasurement() const
Definition: HitPattern.cc:604
IsolatedTracksNxN::t_Bunch
int t_Bunch
Definition: IsolatedTracksNxN.cc:229
IsolatedTracksNxN::t_hsim3x3CharHad
std::vector< double > * t_hsim3x3CharHad
Definition: IsolatedTracksNxN.cc:301
reco::TrackBase::momentum
const Vector & momentum() const
track momentum vector
Definition: TrackBase.h:664
IsolatedTracksNxN::t_trackHitsTOB
std::vector< int > * t_trackHitsTOB
Definition: IsolatedTracksNxN.cc:262
eHCALMatrixExtra.h
IsolatedTracksNxN::t_trackDzPV
std::vector< double > * t_trackDzPV
Definition: IsolatedTracksNxN.cc:257
IsolatedTracksNxN::t_jetPt
std::vector< double > * t_jetPt
Definition: IsolatedTracksNxN.cc:247
IsolatedTracksNxN::h_recPt_2
TH1F * h_recPt_2
Definition: IsolatedTracksNxN.cc:223
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
IsolatedTracksNxN::t_PVNTracksHP
std::vector< int > * t_PVNTracksHP
Definition: IsolatedTracksNxN.cc:245
IsolatedTracksNxN::tok_bs_
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
Definition: IsolatedTracksNxN.cc:193
IsolatedTracksNxN::t_L1TauJetPhi
std::vector< double > * t_L1TauJetPhi
Definition: IsolatedTracksNxN.cc:236
edm::Event
Definition: Event.h:73
IsolatedTracksNxN::t_PVTracksSumPtWt
std::vector< double > * t_PVTracksSumPtWt
Definition: IsolatedTracksNxN.cc:243
IsolatedTracksNxN::tok_L1extIsoEm_
edm::EDGetTokenT< l1extra::L1EmParticleCollection > tok_L1extIsoEm_
Definition: IsolatedTracksNxN.cc:185
MagneticField
Definition: MagneticField.h:19
L1GlobalTriggerObjectMap.h
edm::Log
Definition: MessageLogger.h:70
spr::eCaloSimInfo
void eCaloSimInfo(std::vector< DetId > vdets, const CaloGeometry *geo, edm::Handle< T > &hitsEB, edm::Handle< T > &hitsEE, edm::Handle< edm::SimTrackContainer > &SimTk, edm::Handle< edm::SimVertexContainer > &SimVtx, edm::SimTrackContainer::const_iterator trkInfo, caloSimInfo &info, double timeCut=150, bool debug=false)
IsolatedTracksNxN::tok_jets_
edm::EDGetTokenT< reco::CaloJetCollection > tok_jets_
Definition: IsolatedTracksNxN.cc:188
L1JetParticle.h
reco::TrackBase::dxy
double dxy() const
dxy parameter. (This is the transverse impact parameter w.r.t. to (0,0,0) ONLY if refPoint is close t...
Definition: TrackBase.h:608
IsolatedTracksNxN::t_esim11x11NeutHad
std::vector< double > * t_esim11x11NeutHad
Definition: IsolatedTracksNxN.cc:287
IsolatedTracksNxN::t_trackHitInMeasTEC
std::vector< int > * t_trackHitInMeasTEC
Definition: IsolatedTracksNxN.cc:267
IsolatedTracksNxN::tok_L1extTauJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extTauJet_
Definition: IsolatedTracksNxN.cc:180
TauDecayModes.dec
dec
Definition: TauDecayModes.py:142
SimTrackContainer.h
reco::TrackBase::quality
bool quality(const TrackQuality) const
Track quality.
Definition: TrackBase.h:552
cuy.ii
ii
Definition: cuy.py:589
IsolatedTracksNxN::beginJob
void beginJob() override
Definition: IsolatedTracksNxN.cc:2022
edm::InputTag
Definition: InputTag.h:15
L1GlobalTriggerReadoutRecord.h
IsolatedTracksNxN::t_L1MuonEta
std::vector< double > * t_L1MuonEta
Definition: IsolatedTracksNxN.cc:237
SimVertexContainer.h
IsolatedTracksNxN::maxTrackEta_
const double maxTrackEta_
Definition: IsolatedTracksNxN.cc:174
spr::caloSimInfo::pdgMatched
double pdgMatched
Definition: CaloSimInfo.h:44
IsolatedTracksNxN::h_maxNearP31x31
TH1F * h_maxNearP31x31[NPBins][NEtaBins]
Definition: IsolatedTracksNxN.cc:214
weight
Definition: weight.py:1
EcalChannelStatusRcd.h
IsolatedTracksNxN::t_esim15x15Matched
std::vector< double > * t_esim15x15Matched
Definition: IsolatedTracksNxN.cc:284
IsolatedTracksNxN::initL1_
bool initL1_
Definition: IsolatedTracksNxN.cc:177
IsolatedTracksNxN::t_e15x15_20Sig
std::vector< double > * t_e15x15_20Sig
Definition: IsolatedTracksNxN.cc:278
IsolatedTracksNxN::tok_L1extCenJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extCenJet_
Definition: IsolatedTracksNxN.cc:181
IsolatedTracksNxN::t_e7x7_30Sig
std::vector< double > * t_e7x7_30Sig
Definition: IsolatedTracksNxN.cc:280
L1GtTriggerMenu::gtAlgorithmMap
const AlgorithmMap & gtAlgorithmMap() const
get / set the algorithm map (by name)
Definition: L1GtTriggerMenu.h:187
reco::TrackBase::qoverp
double qoverp() const
q / p
Definition: TrackBase.h:599
IsolatedTracksNxN::t_e7x7_25Sig
std::vector< double > * t_e7x7_25Sig
Definition: IsolatedTracksNxN.cc:279
IsolatedTracksNxN::t_maxNearHcalP7x7
std::vector< double > * t_maxNearHcalP7x7
Definition: IsolatedTracksNxN.cc:290
reco::HitPattern::stripTECLayersWithoutMeasurement
int stripTECLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:666
IsolatedTracksNxN::t_esim9x9
std::vector< double > * t_esim9x9
Definition: IsolatedTracksNxN.cc:283
IsolatedTracksNxN::t_hsim5x5NeutHad
std::vector< double > * t_hsim5x5NeutHad
Definition: IsolatedTracksNxN.cc:300
IsolatedTracksNxN::t_h5x5
std::vector< double > * t_h5x5
Definition: IsolatedTracksNxN.cc:291
DeadROC_duringRun.dir
dir
Definition: DeadROC_duringRun.py:23
IsolatedTracksNxN::t_jetPhi
std::vector< double > * t_jetPhi
Definition: IsolatedTracksNxN.cc:247
IsolatedTracksNxN::t_L1FwdJetPt
std::vector< double > * t_L1FwdJetPt
Definition: IsolatedTracksNxN.cc:235
IsolatedTracksNxN::t_simTrackP
std::vector< double > * t_simTrackP
Definition: IsolatedTracksNxN.cc:282
IsolatedTracksNxN::t_hsim3x3
std::vector< double > * t_hsim3x3
Definition: IsolatedTracksNxN.cc:296
IsolatedTracksNxN::t_esim7x7Rest
std::vector< double > * t_esim7x7Rest
Definition: IsolatedTracksNxN.cc:285
Generated for CMSSW Reference Manual by doxygen 1.8.16