CMS 3D CMS Logo
IsolatedTracksNxN.cc
Go to the documentation of this file.
1 // -*- C++ -*
2 //
3 // Package: IsolatedTracksNxN
4 // Class: IsolatedTracksNxN
5 //
13 //
14 // Original Author: Seema Sharma
15 // Created: Mon Aug 10 15:30:40 CST 2009
16 //
17 //
18 
19 #include "DataFormats/TrackReco/interface/TrackBase.h"
20 #include "Calibration/IsolatedParticles/plugins/IsolatedTracksNxN.h"
21 
22 #include "Calibration/IsolatedParticles/interface/CaloPropagateTrack.h"
23 #include "Calibration/IsolatedParticles/interface/ChargeIsolation.h"
24 #include "Calibration/IsolatedParticles/interface/eHCALMatrixExtra.h"
25 #include "MagneticField/Engine/interface/MagneticField.h"
26 #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
27 #include "CondFormats/DataRecord/interface/EcalChannelStatusRcd.h"
28 #include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"
29 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
30 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"
31 #include "CondFormats/L1TObjects/interface/L1GtTriggerMask.h"
32 #include "CondFormats/DataRecord/interface/L1GtTriggerMaskAlgoTrigRcd.h"
33 #include "CondFormats/DataRecord/interface/L1GtTriggerMaskTechTrigRcd.h"
34 
35 static const bool useL1EventSetup(true);
36 static const bool useL1GtTriggerMenuLite(true);
37 
38 IsolatedTracksNxN::IsolatedTracksNxN(const edm::ParameterSet& iConfig) :
39  m_l1GtUtils(iConfig, consumesCollector(), useL1GtTriggerMenuLite, *this),
40  trackerHitAssociatorConfig_(consumesCollector()) {
41 
42  //now do what ever initialization is needed
43  doMC = iConfig.getUntrackedParameter<bool> ("DoMC", false);
44  writeAllTracks = iConfig.getUntrackedParameter<bool> ("WriteAllTracks", false );
45  myverbose_ = iConfig.getUntrackedParameter<int> ("Verbosity", 5 );
46  pvTracksPtMin_ = iConfig.getUntrackedParameter<double>("PVTracksPtMin", 1.0 );
47  debugTrks_ = iConfig.getUntrackedParameter<int> ("DebugTracks" );
48  printTrkHitPattern_ = iConfig.getUntrackedParameter<bool> ("PrintTrkHitPattern" );
49  minTrackP_ = iConfig.getUntrackedParameter<double>("minTrackP", 1.0 );
50  maxTrackEta_ = iConfig.getUntrackedParameter<double>("maxTrackEta", 5.0 );
51  debugL1Info_ = iConfig.getUntrackedParameter<bool> ("DebugL1Info", false );
52  L1TriggerAlgoInfo_ = iConfig.getUntrackedParameter<bool> ("L1TriggerAlgoInfo");
53  tMinE_ = iConfig.getUntrackedParameter<double>("TimeMinCutECAL", -500.);
54  tMaxE_ = iConfig.getUntrackedParameter<double>("TimeMaxCutECAL", 500.);
55  tMinH_ = iConfig.getUntrackedParameter<double>("TimeMinCutHCAL", -500.);
56  tMaxH_ = iConfig.getUntrackedParameter<double>("TimeMaxCutHCAL", 500.);
57 
58  edm::InputTag L1extraTauJetSource_ = iConfig.getParameter<edm::InputTag> ("L1extraTauJetSource" );
59  edm::InputTag L1extraCenJetSource_ = iConfig.getParameter<edm::InputTag> ("L1extraCenJetSource" );
60  edm::InputTag L1extraFwdJetSource_ = iConfig.getParameter<edm::InputTag> ("L1extraFwdJetSource" );
61  edm::InputTag L1extraMuonSource_ = iConfig.getParameter<edm::InputTag> ("L1extraMuonSource" );
62  edm::InputTag L1extraIsoEmSource_ = iConfig.getParameter<edm::InputTag> ("L1extraIsoEmSource" );
63  edm::InputTag L1extraNonIsoEmSource_ = iConfig.getParameter<edm::InputTag> ("L1extraNonIsoEmSource" );
64  edm::InputTag L1GTReadoutRcdSource_ = iConfig.getParameter<edm::InputTag> ("L1GTReadoutRcdSource" );
65  edm::InputTag L1GTObjectMapRcdSource_= iConfig.getParameter<edm::InputTag> ("L1GTObjectMapRcdSource");
66  edm::InputTag JetSrc_ = iConfig.getParameter<edm::InputTag> ("JetSource");
67  edm::InputTag JetExtender_ = iConfig.getParameter<edm::InputTag> ("JetExtender");
68  edm::InputTag HBHERecHitSource_ = iConfig.getParameter<edm::InputTag> ("HBHERecHitSource");
69 
70  // define tokens for access
71  tok_L1extTauJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraTauJetSource_);
72  tok_L1extCenJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraCenJetSource_);
73  tok_L1extFwdJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraFwdJetSource_);
74  tok_L1extMu_ = consumes<l1extra::L1MuonParticleCollection>(L1extraMuonSource_);
75  tok_L1extIsoEm_ = consumes<l1extra::L1EmParticleCollection>(L1extraIsoEmSource_);
76  tok_L1extNoIsoEm_ = consumes<l1extra::L1EmParticleCollection>(L1extraNonIsoEmSource_);
77  tok_jets_ = consumes<reco::CaloJetCollection>(JetSrc_);
78  tok_hbhe_ = consumes<HBHERecHitCollection>(HBHERecHitSource_);
79 
80  tok_genTrack_ = consumes<reco::TrackCollection>(edm::InputTag("generalTracks"));
81  tok_recVtx_ = consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"));
82  tok_bs_ = consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"));
83  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEB"));
84  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEE"));
85 
86  tok_simTk_ = consumes<edm::SimTrackContainer>(edm::InputTag("g4SimHits"));
87  tok_simVtx_ = consumes<edm::SimVertexContainer>(edm::InputTag("g4SimHits"));
88  tok_caloEB_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "EcalHitsEB"));
89  tok_caloEE_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "EcalHitsEE"));
90  tok_caloHH_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "HcalHits"));
91  nbad = 0;
92 
93  if(myverbose_>=0) {
94  std::cout <<"Parameters read from config file \n"
95  <<" doMC " << doMC
96  <<"\t myverbose_ " << myverbose_
97  <<"\t minTrackP_ " << minTrackP_
98  << "\t maxTrackEta_ " << maxTrackEta_
99  << "\t tMinE_ " << tMinE_
100  << "\t tMaxE_ " << tMaxE_
101  << "\t tMinH_ " << tMinH_
102  << "\t tMaxH_ " << tMaxH_
103  << "\n debugL1Info_ " << debugL1Info_
104  << "\t L1TriggerAlgoInfo_ " << L1TriggerAlgoInfo_
105  << "\n" << std::endl;
106  }
107 
108  initL1 = false;
109 
110 }
111 
112 IsolatedTracksNxN::~IsolatedTracksNxN() {
113 }
114 
115 void IsolatedTracksNxN::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
116 
117  bool haveIsoTrack=false;
118 
119  edm::ESHandle<MagneticField> bFieldH;
120  iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
121  bField = bFieldH.product();
122 
123  clearTreeVectors();
124 
125  t_RunNo = iEvent.id().run();
126  t_EvtNo = iEvent.id().event();
127  t_Lumi = iEvent.luminosityBlock();
128  t_Bunch = iEvent.bunchCrossing();
129 
130  nEventProc++;
131 
132  edm::Handle<reco::TrackCollection> trkCollection;
133  iEvent.getByToken(tok_genTrack_, trkCollection);
134  reco::TrackCollection::const_iterator trkItr;
135  if(debugTrks_>1){
136  std::cout << "Track Collection: " << std::endl;
137  std::cout << "Number of Tracks " << trkCollection->size() << std::endl;
138  }
139  std::string theTrackQuality = "highPurity";
140  reco::TrackBase::TrackQuality trackQuality_=reco::TrackBase::qualityByName(theTrackQuality);
141 
142  //===================== save L1 Trigger information =======================
143  if( L1TriggerAlgoInfo_ ) {
144 
145  m_l1GtUtils.getL1GtRunCache(iEvent, iSetup, useL1EventSetup, useL1GtTriggerMenuLite);
146 
147  int iErrorCode = -1;
148  int l1ConfCode = -1;
149  const bool l1Conf = m_l1GtUtils.availableL1Configuration(iErrorCode, l1ConfCode);
150  if( !l1Conf) {
151  std::cout << "\nL1 configuration code:" << l1ConfCode
152  << "\nNo valid L1 trigger configuration available."
153  << "\nSee text above for error code interpretation"
154  << "\nNo return here, in order to test each method, protected against configuration error."
155  << std::endl;
156  }
157 
158  const L1GtTriggerMenu* m_l1GtMenu = m_l1GtUtils.ptrL1TriggerMenuEventSetup(iErrorCode);
159  const AlgorithmMap& algorithmMap = m_l1GtMenu->gtAlgorithmMap();
160  const std::string& menuName = m_l1GtMenu->gtTriggerMenuName();
161 
162  if (!initL1) {
163  initL1=true;
164  std::cout << "menuName " << menuName << std::endl;
165  for (CItAlgo itAlgo = algorithmMap.begin(); itAlgo != algorithmMap.end(); itAlgo++) {
166  std::string algName = itAlgo->first;
167  int algBitNumber = ( itAlgo->second ).algoBitNumber();
168  l1AlgoMap.insert( std::pair<std::pair<unsigned int,std::string>,int>( std::pair<unsigned int,std::string>(algBitNumber, algName) , 0) ) ;
169  }
170  std::map< std::pair<unsigned int,std::string>, int>::iterator itr;
171  for(itr=l1AlgoMap.begin(); itr!=l1AlgoMap.end(); itr++) {
172  std::cout << " ********** " << (itr->first).first <<" "<<(itr->first).second <<" "<<itr->second << std::endl;
173  }
174  }
175 
176  std::vector<int> algbits;
177  for (CItAlgo itAlgo = algorithmMap.begin(); itAlgo != algorithmMap.end(); itAlgo++) {
178  std::string algName = itAlgo->first;
179  int algBitNumber = ( itAlgo->second ).algoBitNumber();
180  bool decision = m_l1GtUtils.decision (iEvent, itAlgo->first, iErrorCode);
181  int preScale = m_l1GtUtils.prescaleFactor (iEvent, itAlgo->first, iErrorCode);
182 
183  // save the algo names which fired
184  if( decision ){
185  l1AlgoMap[std::pair<unsigned int,std::string>(algBitNumber, algName)] += 1;
186  h_L1AlgoNames->Fill( algBitNumber );
187  t_L1AlgoNames->push_back(itAlgo->first);
188  t_L1PreScale->push_back(preScale);
189  t_L1Decision[algBitNumber] = 1;
190  algbits.push_back(algBitNumber);
191  }
192  }
193 
194  if(debugL1Info_) {
195  for(unsigned int ii=0; ii<t_L1AlgoNames->size(); ii++){
196  std::cout<<ii<<" "<<(*t_L1AlgoNames)[ii]<<" "<<(*t_L1PreScale)[ii]<<" "<<algbits[ii]<<std::endl;
197  }
198  std::cout << "L1Decision: ";
199  for (int i=0; i<128; ++i) {
200  std::cout << " " << i << ":" << t_L1Decision[i];
201  if ((i+1)%32 == 0) std::cout << std::endl;
202  }
203  }
204 
205  // L1Taus
206  edm::Handle<l1extra::L1JetParticleCollection> l1TauHandle;
207  iEvent.getByToken(tok_L1extTauJet_,l1TauHandle);
208  l1extra::L1JetParticleCollection::const_iterator itr;
209  int iL1Obj=0;
210  for(itr = l1TauHandle->begin(),iL1Obj=0; itr != l1TauHandle->end(); ++itr,iL1Obj++) {
211  if(iL1Obj<1) {
212  t_L1TauJetPt ->push_back( itr->pt() );
213  t_L1TauJetEta ->push_back( itr->eta() );
214  t_L1TauJetPhi ->push_back( itr->phi() );
215  }
216  if(debugL1Info_) {
217  std::cout << "tauJ p/pt " << itr->momentum() << " " << itr->pt()
218  << " eta/phi " << itr->eta() << " " << itr->phi()
219  << std::endl;
220  }
221  }
222 
223  // L1 Central Jets
224  edm::Handle<l1extra::L1JetParticleCollection> l1CenJetHandle;
225  iEvent.getByToken(tok_L1extCenJet_,l1CenJetHandle);
226  for( itr = l1CenJetHandle->begin(),iL1Obj=0; itr != l1CenJetHandle->end(); ++itr,iL1Obj++ ) {
227  if(iL1Obj<1) {
228  t_L1CenJetPt ->push_back( itr->pt() );
229  t_L1CenJetEta ->push_back( itr->eta() );
230  t_L1CenJetPhi ->push_back( itr->phi() );
231  }
232  if(debugL1Info_) {
233  std::cout << "cenJ p/pt " << itr->momentum() << " " << itr->pt()
234  << " eta/phi " << itr->eta() << " " << itr->phi()
235  << std::endl;
236  }
237  }
238 
239  // L1 Forward Jets
240  edm::Handle<l1extra::L1JetParticleCollection> l1FwdJetHandle;
241  iEvent.getByToken(tok_L1extFwdJet_,l1FwdJetHandle);
242  for( itr = l1FwdJetHandle->begin(),iL1Obj=0; itr != l1FwdJetHandle->end(); ++itr,iL1Obj++ ) {
243  if(iL1Obj<1) {
244  t_L1FwdJetPt ->push_back( itr->pt() );
245  t_L1FwdJetEta ->push_back( itr->eta() );
246  t_L1FwdJetPhi ->push_back( itr->phi() );
247  }
248  if(debugL1Info_) {
249  std::cout << "fwdJ p/pt " << itr->momentum() << " " << itr->pt()
250  << " eta/phi " << itr->eta() << " " << itr->phi()
251  << std::endl;
252  }
253  }
254 
255  // L1 Isolated EM onjects
256  l1extra::L1EmParticleCollection::const_iterator itrEm;
257  edm::Handle<l1extra::L1EmParticleCollection> l1IsoEmHandle ;
258  iEvent.getByToken(tok_L1extIsoEm_, l1IsoEmHandle);
259  for( itrEm = l1IsoEmHandle->begin(),iL1Obj=0; itrEm != l1IsoEmHandle->end(); ++itrEm,iL1Obj++ ) {
260  if(iL1Obj<1) {
261  t_L1IsoEMPt ->push_back( itrEm->pt() );
262  t_L1IsoEMEta ->push_back( itrEm->eta() );
263  t_L1IsoEMPhi ->push_back( itrEm->phi() );
264  }
265  if(debugL1Info_) {
266  std::cout << "isoEm p/pt " << itrEm->momentum() << " " << itrEm->pt()
267  << " eta/phi " << itrEm->eta() << " " << itrEm->phi()
268  << std::endl;
269  }
270  }
271 
272  // L1 Non-Isolated EM onjects
273  edm::Handle<l1extra::L1EmParticleCollection> l1NonIsoEmHandle ;
274  iEvent.getByToken(tok_L1extNoIsoEm_, l1NonIsoEmHandle);
275  for( itrEm = l1NonIsoEmHandle->begin(),iL1Obj=0; itrEm != l1NonIsoEmHandle->end(); ++itrEm,iL1Obj++ ) {
276  if(iL1Obj<1) {
277  t_L1NonIsoEMPt ->push_back( itrEm->pt() );
278  t_L1NonIsoEMEta ->push_back( itrEm->eta() );
279  t_L1NonIsoEMPhi ->push_back( itrEm->phi() );
280  }
281  if(debugL1Info_) {
282  std::cout << "nonIsoEm p/pt " << itrEm->momentum() << " " << itrEm->pt()
283  << " eta/phi " << itrEm->eta() << " " << itrEm->phi()
284  << std::endl;
285  }
286  }
287 
288  // L1 Muons
289  l1extra::L1MuonParticleCollection::const_iterator itrMu;
290  edm::Handle<l1extra::L1MuonParticleCollection> l1MuHandle ;
291  iEvent.getByToken(tok_L1extMu_, l1MuHandle);
292  for( itrMu = l1MuHandle->begin(),iL1Obj=0; itrMu != l1MuHandle->end(); ++itrMu,iL1Obj++ ) {
293  if(iL1Obj<1) {
294  t_L1MuonPt ->push_back( itrMu->pt() );
295  t_L1MuonEta ->push_back( itrMu->eta() );
296  t_L1MuonPhi ->push_back( itrMu->phi() );
297  }
298  if(debugL1Info_) {
299  std::cout << "l1muon p/pt " << itrMu->momentum() << " " << itrMu->pt()
300  << " eta/phi " << itrMu->eta() << " " << itrMu->phi()
301  << std::endl;
302  }
303  }
304  }
305 
306  //============== store the information about all the Non-Fake vertices ===============
307 
308  edm::Handle<reco::VertexCollection> recVtxs;
309  iEvent.getByToken(tok_recVtx_,recVtxs);
310 
311  std::vector<reco::Track> svTracks;
312  math::XYZPoint leadPV(0,0,0);
313  double sumPtMax = -1.0;
314  for(unsigned int ind=0;ind<recVtxs->size();ind++) {
315 
316  if ( !((*recVtxs)[ind].isFake()) ) {
317  double vtxTrkSumPt=0.0, vtxTrkSumPtWt=0.0;
318  int vtxTrkNWt =0;
319  double vtxTrkSumPtHP=0.0, vtxTrkSumPtHPWt=0.0;
320  int vtxTrkNHP =0, vtxTrkNHPWt =0;
321 
322  reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
323 
324  for (vtxTrack = (*recVtxs)[ind].tracks_begin(); vtxTrack!=(*recVtxs)[ind].tracks_end(); vtxTrack++) {
325 
326  if((*vtxTrack)->pt()<pvTracksPtMin_) continue;
327 
328  bool trkQuality = (*vtxTrack)->quality(trackQuality_);
329 
330  vtxTrkSumPt += (*vtxTrack)->pt();
331 
332  vtxTrkSumPt += (*vtxTrack)->pt();
333  if( trkQuality ) {
334  vtxTrkSumPtHP += (*vtxTrack)->pt();
335  vtxTrkNHP++;
336  }
337 
338  double weight = (*recVtxs)[ind].trackWeight(*vtxTrack);
339  h_PVTracksWt ->Fill(weight);
340  if( weight>0.5) {
341  vtxTrkSumPtWt += (*vtxTrack)->pt();
342  vtxTrkNWt++;
343  if( trkQuality ) {
344  vtxTrkSumPtHPWt += (*vtxTrack)->pt();
345  vtxTrkNHPWt++;
346  }
347  }
348  }
349 
350  if(vtxTrkSumPt>sumPtMax) {
351  sumPtMax = vtxTrkSumPt;
352  leadPV = math::XYZPoint( (*recVtxs)[ind].x(),(*recVtxs)[ind].y(), (*recVtxs)[ind].z() );
353  }
354 
355  t_PVx ->push_back( (*recVtxs)[ind].x() );
356  t_PVy ->push_back( (*recVtxs)[ind].y() );
357  t_PVz ->push_back( (*recVtxs)[ind].z() );
358  t_PVisValid ->push_back( (*recVtxs)[ind].isValid() );
359  t_PVNTracks ->push_back( (*recVtxs)[ind].tracksSize() );
360  t_PVndof ->push_back( (*recVtxs)[ind].ndof() );
361  t_PVNTracksWt ->push_back( vtxTrkNWt );
362  t_PVTracksSumPt ->push_back( vtxTrkSumPt );
363  t_PVTracksSumPtWt->push_back( vtxTrkSumPtWt );
364 
365  t_PVNTracksHP ->push_back( vtxTrkNHP );
366  t_PVNTracksHPWt ->push_back( vtxTrkNHPWt );
367  t_PVTracksSumPtHP ->push_back( vtxTrkSumPtHP );
368  t_PVTracksSumPtHPWt->push_back( vtxTrkSumPtHPWt );
369 
370  if(myverbose_==4) {
371  std::cout<<"PV "<<ind<<" isValid "<<(*recVtxs)[ind].isValid()<<" isFake "<<(*recVtxs)[ind].isFake()
372  <<" hasRefittedTracks() "<<ind<<" "<<(*recVtxs)[ind].hasRefittedTracks()
373  <<" refittedTrksSize "<<(*recVtxs)[ind].refittedTracks().size()
374  <<" tracksSize() "<<(*recVtxs)[ind].tracksSize()<<" sumPt "<<vtxTrkSumPt
375  <<std::endl;
376  }
377  } // if vtx is not Fake
378  } // loop over PVs
379  //===================================================================================
380 
381  // Get the beamspot
382  edm::Handle<reco::BeamSpot> beamSpotH;
383  iEvent.getByToken(tok_bs_, beamSpotH);
384  math::XYZPoint bspot;
385  bspot = ( beamSpotH.isValid() ) ? beamSpotH->position() : math::XYZPoint(0, 0, 0);
386 
387  //=====================================================================
388 
389  edm::Handle<reco::CaloJetCollection> jets;
390  iEvent.getByToken(tok_jets_,jets);
391  // edm::Handle<reco::JetExtendedAssociation::Container> jetExtender;
392  // iEvent.getByLabel(JetExtender_,jetExtender);
393 
394  for(unsigned int ijet=0;ijet<(*jets).size();ijet++) {
395  t_jetPt ->push_back( (*jets)[ijet].pt() );
396  t_jetEta ->push_back( (*jets)[ijet].eta() );
397  t_jetPhi ->push_back( (*jets)[ijet].phi() );
398  t_nTrksJetVtx ->push_back( -1.0);
399  t_nTrksJetCalo->push_back( -1.0);
400  }
401 
402  //===================================================================================
403 
404  // get handles to calogeometry and calotopology
405  edm::ESHandle<CaloGeometry> pG;
406  iSetup.get<CaloGeometryRecord>().get(pG);
407  const CaloGeometry* geo = pG.product();
408 
409  edm::ESHandle<CaloTopology> theCaloTopology;
410  iSetup.get<CaloTopologyRecord>().get(theCaloTopology);
411  const CaloTopology *caloTopology = theCaloTopology.product();
412 
413  edm::ESHandle<HcalTopology> htopo;
414  iSetup.get<HcalRecNumberingRecord>().get(htopo);
415  const HcalTopology* theHBHETopology = htopo.product();
416 
417  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
418  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
419  iEvent.getByToken(tok_EB_,barrelRecHitsHandle);
420  iEvent.getByToken(tok_EE_,endcapRecHitsHandle);
421 
422  // Retrieve the good/bad ECAL channels from the DB
423  edm::ESHandle<EcalChannelStatus> ecalChStatus;
424  iSetup.get<EcalChannelStatusRcd>().get(ecalChStatus);
425  const EcalChannelStatus* theEcalChStatus = ecalChStatus.product();
426 
427  // Retrieve trigger tower map
428  edm::ESHandle<EcalTrigTowerConstituentsMap> hTtmap;
429  iSetup.get<IdealGeometryRecord>().get(hTtmap);
430  const EcalTrigTowerConstituentsMap& ttMap = *hTtmap;
431 
432  edm::Handle<HBHERecHitCollection> hbhe;
433  iEvent.getByToken(tok_hbhe_, hbhe);
434  if (!hbhe.isValid()) {
435  nbad++;
436  if (nbad < 10) std::cout << "No HBHE rechit collection\n";
437  return;
438  }
439  const HBHERecHitCollection Hithbhe = *(hbhe.product());
440 
441  //get Handles to SimTracks and SimHits
442  edm::Handle<edm::SimTrackContainer> SimTk;
443  if (doMC) iEvent.getByToken(tok_simTk_,SimTk);
444  edm::SimTrackContainer::const_iterator simTrkItr;
445 
446  edm::Handle<edm::SimVertexContainer> SimVtx;
447  if (doMC) iEvent.getByToken(tok_simVtx_,SimVtx);
448 
449  //get Handles to PCaloHitContainers of eb/ee/hbhe
450  edm::Handle<edm::PCaloHitContainer> pcaloeb;
451  if (doMC) iEvent.getByToken(tok_caloEB_, pcaloeb);
452 
453  edm::Handle<edm::PCaloHitContainer> pcaloee;
454  if (doMC) iEvent.getByToken(tok_caloEE_, pcaloee);
455 
456  edm::Handle<edm::PCaloHitContainer> pcalohh;
457  if (doMC) iEvent.getByToken(tok_caloHH_, pcalohh);
458 
459  //associates tracker rechits/simhits to a track
460  std::unique_ptr<TrackerHitAssociator> associate;
461  if (doMC) associate.reset(new TrackerHitAssociator(iEvent, trackerHitAssociatorConfig_));
462 
463  //===================================================================================
464 
465  h_nTracks->Fill(trkCollection->size());
466 
467  int nTracks = 0;
468 
469  t_nTracks = trkCollection->size();
470 
471  // get the list of DetIds closest to the impact point of track on surface calorimeters
472  std::vector<spr::propagatedTrackID> trkCaloDets;
473  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality, trkCaloDets, false);
474  std::vector<spr::propagatedTrackID>::const_iterator trkDetItr;
475 
476  if(myverbose_>2) {
477  for(trkDetItr = trkCaloDets.begin(); trkDetItr != trkCaloDets.end(); trkDetItr++){
478  std::cout<<trkDetItr->trkItr->p()<<" "<<trkDetItr->trkItr->eta()<<" "<<trkDetItr->okECAL<<" ";
479  if(trkDetItr->detIdECAL.subdetId() == EcalBarrel) std::cout << (EBDetId)trkDetItr->detIdECAL <<" ";
480  else std::cout << (EEDetId)trkDetItr->detIdECAL <<" ";
481  std::cout<<trkDetItr->okHCAL<<" ";
482  if(trkDetItr->okHCAL) std::cout<<(HcalDetId)trkDetItr->detIdHCAL;
483  std::cout << std::endl;
484  }
485  }
486 
487  int nvtxTracks=0;
488  for(trkDetItr = trkCaloDets.begin(),nTracks=0; trkDetItr != trkCaloDets.end(); trkDetItr++,nTracks++){
489 
490  const reco::Track* pTrack = &(*(trkDetItr->trkItr));
491 
492  // find vertex index the track is associated with
493  int pVtxTkId = -1;
494  for(unsigned int ind=0; ind<recVtxs->size(); ind++) {
495  if (!((*recVtxs)[ind].isFake())) {
496  reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
497  for (vtxTrack = (*recVtxs)[ind].tracks_begin(); vtxTrack!=(*recVtxs)[ind].tracks_end(); vtxTrack++) {
498 
499  const edm::RefToBase<reco::Track> pvtxTrack = (*vtxTrack);
500  if ( pTrack == pvtxTrack.get() ) {
501  pVtxTkId = ind;
502  break;
503  if(myverbose_>2) {
504  if( pTrack->pt()>1.0) {
505  std::cout<<"Debug the track association with vertex "<<std::endl;
506  std::cout<< pTrack << " "<< pvtxTrack.get() << std::endl;
507  std::cout<<" trkVtxIndex "<<nvtxTracks<<" vtx "<<ind<<" pt "<<pTrack->pt()
508  <<" eta "<<pTrack->eta()<<" "<<pTrack->pt()-pvtxTrack->pt()
509  <<" "<< pTrack->eta()-pvtxTrack->eta()
510  << std::endl;
511  nvtxTracks++;
512  }
513  }
514  }
515  }
516  }
517  }
518 
519  const reco::HitPattern& hitp = pTrack->hitPattern();
520  int nLayersCrossed = hitp.trackerLayersWithMeasurement() ;
521  int nOuterHits = hitp.stripTOBLayersWithMeasurement()+hitp.stripTECLayersWithMeasurement() ;
522 
523  bool ifGood = pTrack->quality(trackQuality_);
524  double pt1 = pTrack->pt();
525  double p1 = pTrack->p();
526  double eta1 = pTrack->momentum().eta();
527  double phi1 = pTrack->momentum().phi();
528  double etaEcal1 = trkDetItr->etaECAL;
529  double phiEcal1 = trkDetItr->phiECAL;
530  double etaHcal1 = trkDetItr->etaHCAL;
531  double phiHcal1 = trkDetItr->phiHCAL;
532  double dxy1 = pTrack->dxy();
533  double dz1 = pTrack->dz();
534  double chisq1 = pTrack->normalizedChi2();
535  double dxybs1 = beamSpotH.isValid() ? pTrack->dxy(bspot) : pTrack->dxy();
536  double dzbs1 = beamSpotH.isValid() ? pTrack->dz(bspot) : pTrack->dz();
537  double dxypv1 = pTrack->dxy();
538  double dzpv1 = pTrack->dz();
539  if(pVtxTkId>=0) {
540  math::XYZPoint thisTkPV = math::XYZPoint( (*recVtxs)[pVtxTkId].x(),(*recVtxs)[pVtxTkId].y(), (*recVtxs)[pVtxTkId].z() );
541  dxypv1 = pTrack->dxy(thisTkPV);
542  dzpv1 = pTrack->dz (thisTkPV);
543  }
544 
545  h_recEtaPt_0->Fill(eta1, pt1);
546  h_recEtaP_0 ->Fill(eta1, p1);
547  h_recPt_0 ->Fill(pt1);
548  h_recP_0 ->Fill(p1);
549  h_recEta_0 ->Fill(eta1);
550  h_recPhi_0 ->Fill(phi1);
551 
552  if(ifGood && nLayersCrossed>7 ) {
553  h_recEtaPt_1->Fill(eta1, pt1);
554  h_recEtaP_1 ->Fill(eta1, p1);
555  h_recPt_1 ->Fill(pt1);
556  h_recP_1 ->Fill(p1);
557  h_recEta_1 ->Fill(eta1);
558  h_recPhi_1 ->Fill(phi1);
559  }
560 
561  if( ! ifGood ) continue;
562 
563  if( writeAllTracks && p1>2.0 && nLayersCrossed>7) {
564  t_trackPAll ->push_back( p1 );
565  t_trackEtaAll ->push_back( eta1 );
566  t_trackPhiAll ->push_back( phi1 );
567  t_trackPtAll ->push_back( pt1 );
568  t_trackDxyAll ->push_back( dxy1 );
569  t_trackDzAll ->push_back( dz1 );
570  t_trackDxyPVAll ->push_back( dxypv1 );
571  t_trackDzPVAll ->push_back( dzpv1 );
572  t_trackChiSqAll ->push_back( chisq1 );
573  }
574  if (doMC) {
575  edm::SimTrackContainer::const_iterator matchedSimTrkAll = spr::matchedSimTrack(iEvent, SimTk, SimVtx, pTrack, *associate, false);
576  if( writeAllTracks && matchedSimTrkAll != SimTk->end()) t_trackPdgIdAll->push_back( matchedSimTrkAll->type() );
577  }
578 
579  if( pt1>minTrackP_ && std::abs(eta1)<maxTrackEta_ && trkDetItr->okECAL) {
580 
581  double maxNearP31x31=999.0, maxNearP25x25=999.0, maxNearP21x21=999.0, maxNearP15x15=999.0;
582  maxNearP31x31 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 15,15);
583  maxNearP25x25 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 12,12);
584  maxNearP21x21 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 10,10);
585  maxNearP15x15 = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 7, 7);
586 
587  int iTrkEtaBin=-1, iTrkMomBin=-1;
588  for(unsigned int ieta=0; ieta<NEtaBins; ieta++) {
589  if(std::abs(eta1)>genPartEtaBins[ieta] && std::abs(eta1)<genPartEtaBins[ieta+1] ) iTrkEtaBin = ieta;
590  }
591  for(unsigned int ipt=0; ipt<NPBins; ipt++) {
592  if( p1>genPartPBins[ipt] && p1<genPartPBins[ipt+1] ) iTrkMomBin = ipt;
593  }
594  if( iTrkMomBin>=0 && iTrkEtaBin>=0 ) {
595  h_maxNearP31x31[iTrkMomBin][iTrkEtaBin]->Fill( maxNearP31x31 );
596  h_maxNearP25x25[iTrkMomBin][iTrkEtaBin]->Fill( maxNearP25x25 );
597  h_maxNearP21x21[iTrkMomBin][iTrkEtaBin]->Fill( maxNearP21x21 );
598  h_maxNearP15x15[iTrkMomBin][iTrkEtaBin]->Fill( maxNearP15x15 );
599  }
600  if( maxNearP31x31<0.0 && nLayersCrossed>7 && nOuterHits>4) {
601  h_recEtaPt_2->Fill(eta1, pt1);
602  h_recEtaP_2 ->Fill(eta1, p1);
603  h_recPt_2 ->Fill(pt1);
604  h_recP_2 ->Fill(p1);
605  h_recEta_2 ->Fill(eta1);
606  h_recPhi_2 ->Fill(phi1);
607  }
608 
609  // if charge isolated in ECAL, store the further quantities
610  if( maxNearP31x31<1.0) {
611 
612  haveIsoTrack = true;
613 
614  // get the matching simTrack
615  double simTrackP = -1;
616  if (doMC) {
617  edm::SimTrackContainer::const_iterator matchedSimTrk = spr::matchedSimTrack(iEvent, SimTk, SimVtx, pTrack, *associate, false);
618  if( matchedSimTrk != SimTk->end() )simTrackP = matchedSimTrk->momentum().P();
619  }
620 
621  edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
622  iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
623 
624  // get ECal Tranverse Profile
625  std::pair<double, bool> e7x7P, e9x9P, e11x11P, e15x15P;
626  std::pair<double, bool> e7x7_10SigP, e9x9_10SigP, e11x11_10SigP, e15x15_10SigP;
627  std::pair<double, bool> e7x7_15SigP, e9x9_15SigP, e11x11_15SigP, e15x15_15SigP;
628  std::pair<double, bool> e7x7_20SigP, e9x9_20SigP, e11x11_20SigP, e15x15_20SigP;
629  std::pair<double, bool> e7x7_25SigP, e9x9_25SigP, e11x11_25SigP, e15x15_25SigP;
630  std::pair<double, bool> e7x7_30SigP, e9x9_30SigP, e11x11_30SigP, e15x15_30SigP;
631 
632  spr::caloSimInfo simInfo3x3, simInfo5x5, simInfo7x7, simInfo9x9;
633  spr::caloSimInfo simInfo11x11, simInfo13x13, simInfo15x15, simInfo21x21, simInfo25x25, simInfo31x31;
634  double trkEcalEne=0;
635 
636  const DetId isoCell = trkDetItr->detIdECAL;
637  e7x7P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),3,3, -100.0, -100.0, tMinE_,tMaxE_);
638  e9x9P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),4,4, -100.0, -100.0, tMinE_,tMaxE_);
639  e11x11P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5, -100.0, -100.0, tMinE_,tMaxE_);
640  e15x15P = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7, -100.0, -100.0, tMinE_,tMaxE_);
641 
642  e7x7_10SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),3,3, 0.030, 0.150, tMinE_,tMaxE_);
643  e9x9_10SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),4,4, 0.030, 0.150, tMinE_,tMaxE_);
644  e11x11_10SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5, 0.030, 0.150, tMinE_,tMaxE_);
645  e15x15_10SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7, 0.030, 0.150, tMinE_,tMaxE_);
646 
647  e7x7_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, sevlv.product(),ttMap, 3,3, 0.20,0.45, tMinE_,tMaxE_);
648  e9x9_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, sevlv.product(),ttMap, 4,4, 0.20,0.45, tMinE_,tMaxE_);
649  e11x11_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(), ttMap, 5,5, 0.20,0.45, tMinE_,tMaxE_);
650  e15x15_15SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, sevlv.product(),ttMap, 7,7, 0.20,0.45, tMinE_,tMaxE_, false);
651 
652  e7x7_20SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),3,3, 0.060, 0.300, tMinE_,tMaxE_);
653  e9x9_20SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),4,4, 0.060, 0.300, tMinE_,tMaxE_);
654  e11x11_20SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5, 0.060, 0.300, tMinE_,tMaxE_);
655  e15x15_20SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7, 0.060, 0.300, tMinE_,tMaxE_);
656 
657  e7x7_25SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),3,3, 0.075, 0.375, tMinE_,tMaxE_);
658  e9x9_25SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),4,4, 0.075, 0.375, tMinE_,tMaxE_);
659  e11x11_25SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5, 0.075, 0.375, tMinE_,tMaxE_);
660  e15x15_25SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7, 0.075, 0.375, tMinE_,tMaxE_);
661 
662  e7x7_30SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),3,3, 0.090, 0.450, tMinE_,tMaxE_);
663  e9x9_30SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),4,4, 0.090, 0.450, tMinE_,tMaxE_);
664  e11x11_30SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5, 0.090, 0.450, tMinE_,tMaxE_);
665  e15x15_30SigP = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7, 0.090, 0.450, tMinE_,tMaxE_);
666  if(myverbose_ == 2) {
667  std::cout << "clean ecal rechit " << std::endl;
668  std::cout<<"e7x7 "<<e7x7P.first<<" e9x9 "<<e9x9P.first<<" e11x11 " << e11x11P.first << " e15x15 "<<e15x15P.first<<std::endl;
669  std::cout<<"e7x7_10Sig "<<e7x7_10SigP.first<<" e11x11_10Sig "<<e11x11_10SigP.first<<" e15x15_10Sig "<<e15x15_10SigP.first<<std::endl;
670  }
671 
672  if (doMC) {
673  // check the energy from SimHits
674  spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 1,1, simInfo3x3);
675  spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 2,2, simInfo5x5);
676  spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 3,3, simInfo7x7);
677  spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 4,4, simInfo9x9);
678  spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 5,5, simInfo11x11);
679  spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 6,6, simInfo13x13);
680  spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 7,7, simInfo15x15, 150.0,false);
681  spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 10,10, simInfo21x21);
682  spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 12,12, simInfo25x25);
683  spr::eECALSimInfo(iEvent,isoCell,geo,caloTopology,pcaloeb,pcaloee,SimTk,SimVtx,pTrack, *associate, 15,15, simInfo31x31);
684 
685  trkEcalEne = spr::eCaloSimInfo(iEvent, geo, pcaloeb,pcaloee, SimTk, SimVtx, pTrack, *associate);
686  if(myverbose_ == 1) {
687  std::cout << "Track momentum " << pt1 << std::endl;
688 
689  std::cout << "ecal siminfo " << std::endl;
690  std::cout << "simInfo3x3: " << "eTotal " << simInfo3x3.eTotal << " eMatched " << simInfo3x3.eMatched << " eRest " << simInfo3x3.eRest << " eGamma "<<simInfo3x3.eGamma << " eNeutralHad " << simInfo3x3.eNeutralHad << " eChargedHad " << simInfo3x3.eChargedHad << std::endl;
691  std::cout << "simInfo5x5: " << "eTotal " << simInfo5x5.eTotal << " eMatched " << simInfo5x5.eMatched << " eRest " << simInfo5x5.eRest << " eGamma "<<simInfo5x5.eGamma << " eNeutralHad " << simInfo5x5.eNeutralHad << " eChargedHad " << simInfo5x5.eChargedHad << std::endl;
692  std::cout << "simInfo7x7: " << "eTotal " << simInfo7x7.eTotal << " eMatched " << simInfo7x7.eMatched << " eRest " << simInfo7x7.eRest << " eGamma "<<simInfo7x7.eGamma << " eNeutralHad " << simInfo7x7.eNeutralHad << " eChargedHad " << simInfo7x7.eChargedHad << std::endl;
693  std::cout << "simInfo9x9: " << "eTotal " << simInfo9x9.eTotal << " eMatched " << simInfo9x9.eMatched << " eRest " << simInfo9x9.eRest << " eGamma "<<simInfo9x9.eGamma << " eNeutralHad " << simInfo9x9.eNeutralHad << " eChargedHad " << simInfo9x9.eChargedHad << std::endl;
694  std::cout << "simInfo11x11: " << "eTotal " << simInfo11x11.eTotal << " eMatched " << simInfo11x11.eMatched << " eRest " << simInfo11x11.eRest << " eGamma "<<simInfo11x11.eGamma << " eNeutralHad " << simInfo11x11.eNeutralHad << " eChargedHad " << simInfo11x11.eChargedHad << std::endl;
695  std::cout << "simInfo15x15: " << "eTotal " << simInfo15x15.eTotal << " eMatched " << simInfo15x15.eMatched << " eRest " << simInfo15x15.eRest << " eGamma "<<simInfo15x15.eGamma << " eNeutralHad " << simInfo15x15.eNeutralHad << " eChargedHad " << simInfo15x15.eChargedHad << std::endl;
696  std::cout << "simInfo31x31: " << "eTotal " << simInfo31x31.eTotal << " eMatched " << simInfo31x31.eMatched << " eRest " << simInfo31x31.eRest << " eGamma "<<simInfo31x31.eGamma << " eNeutralHad " << simInfo31x31.eNeutralHad << " eChargedHad " << simInfo31x31.eChargedHad << std::endl;
697  std::cout << "trkEcalEne" << trkEcalEne << std::endl;
698  }
699  }
700 
701  // ======= Get HCAL information
702  double hcalScale=1.0;
703  if( std::abs(pTrack->eta())<1.4 ) {
704  hcalScale=120.0;
705  } else {
706  hcalScale=135.0;
707  }
708 
709  double maxNearHcalP3x3=-1, maxNearHcalP5x5=-1, maxNearHcalP7x7=-1;
710  maxNearHcalP3x3 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 1,1);
711  maxNearHcalP5x5 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 2,2);
712  maxNearHcalP7x7 = spr::chargeIsolationHcal(nTracks, trkCaloDets, theHBHETopology, 3,3);
713 
714  double h3x3=0, h5x5=0, h7x7=0;
715  double h3x3Sig=0, h5x5Sig=0, h7x7Sig=0;
716  double trkHcalEne = 0;
717  spr::caloSimInfo hsimInfo3x3, hsimInfo5x5, hsimInfo7x7;
718 
719  if(trkDetItr->okHCAL) {
720  const DetId ClosestCell(trkDetItr->detIdHCAL);
721  // bool includeHO=false, bool algoNew=true, bool debug=false
722  h3x3 = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,1,1, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_,tMaxH_);
723  h5x5 = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,2,2, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_,tMaxH_);
724  h7x7 = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,3,3, false, true, -100.0, -100.0, -100.0, -100.0, tMinH_,tMaxH_);
725  h3x3Sig = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,1,1, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_,tMaxH_);
726  h5x5Sig = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,2,2, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_,tMaxH_);
727  h7x7Sig = spr::eHCALmatrix(theHBHETopology, ClosestCell, hbhe,3,3, false, true, 0.7, 0.8, -100.0, -100.0, tMinH_,tMaxH_);
728 
729  if(myverbose_==2) {
730  std::cout << "HCAL 3x3 " << h3x3 << " " << h3x3Sig << " 5x5 " << h5x5 << " " << h5x5Sig << " 7x7 " << h7x7 << " " << h7x7Sig << std::endl;
731  }
732 
733  if (doMC) {
734  spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, geo,pcalohh, SimTk, SimVtx, pTrack, *associate, 1,1, hsimInfo3x3);
735  spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, geo,pcalohh, SimTk, SimVtx, pTrack, *associate, 2,2, hsimInfo5x5);
736  spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, geo,pcalohh, SimTk, SimVtx, pTrack, *associate, 3,3, hsimInfo7x7, 150.0, false,false);
737  trkHcalEne = spr::eCaloSimInfo(iEvent, geo,pcalohh, SimTk, SimVtx, pTrack, *associate);
738  if(myverbose_ == 1) {
739  std::cout << "Hcal siminfo " << std::endl;
740  std::cout << "hsimInfo3x3: " << "eTotal " << hsimInfo3x3.eTotal << " eMatched " << hsimInfo3x3.eMatched << " eRest " << hsimInfo3x3.eRest << " eGamma "<<hsimInfo3x3.eGamma << " eNeutralHad " << hsimInfo3x3.eNeutralHad << " eChargedHad " << hsimInfo3x3.eChargedHad << std::endl;
741  std::cout << "hsimInfo5x5: " << "eTotal " << hsimInfo5x5.eTotal << " eMatched " << hsimInfo5x5.eMatched << " eRest " << hsimInfo5x5.eRest << " eGamma "<<hsimInfo5x5.eGamma << " eNeutralHad " << hsimInfo5x5.eNeutralHad << " eChargedHad " << hsimInfo5x5.eChargedHad << std::endl;
742  std::cout << "hsimInfo7x7: " << "eTotal " << hsimInfo7x7.eTotal << " eMatched " << hsimInfo7x7.eMatched << " eRest " << hsimInfo7x7.eRest << " eGamma "<<hsimInfo7x7.eGamma << " eNeutralHad " << hsimInfo7x7.eNeutralHad << " eChargedHad " << hsimInfo7x7.eChargedHad << std::endl;
743  std::cout << "trkHcalEne " << trkHcalEne << std::endl;
744  }
745  }
746 
747  // debug the ecal and hcal matrix
748  if(myverbose_==4) {
749  std::cout<<"Run "<<iEvent.id().run()<<" Event "<<iEvent.id().event()<<std::endl;
750  std::vector<std::pair<DetId,double> > v7x7 = spr::eHCALmatrixCell(theHBHETopology, ClosestCell, hbhe,3,3, false, false);
751  double sumv=0.0;
752 
753  for(unsigned int iv=0; iv<v7x7.size(); iv++) {
754  sumv += v7x7[iv].second;
755  }
756  std::cout<<"h7x7 "<<h7x7<<" v7x7 "<<sumv << " in " << v7x7.size() <<std::endl;
757  for(unsigned int iv=0; iv<v7x7.size(); iv++) {
758  HcalDetId id = v7x7[iv].first;
759  std::cout << " Cell " << iv << " 0x" << std::hex << v7x7[iv].first() << std::dec << " " << id << " Energy " << v7x7[iv].second << std::endl;
760  }
761  }
762 
763  }
764  if (doMC) {
765  trkHcalEne = spr::eCaloSimInfo(iEvent, geo,pcalohh, SimTk, SimVtx, pTrack, *associate);
766  }
767 
768  // ====================================================================================================
769  // get diff between track outermost hit position and the propagation point at outermost surface of tracker
770  std::pair<math::XYZPoint,double> point2_TK0 = spr::propagateTrackerEnd( pTrack, bField, false);
771  math::XYZPoint diff(pTrack->outerPosition().X()-point2_TK0.first.X(),
772  pTrack->outerPosition().Y()-point2_TK0.first.Y(),
773  pTrack->outerPosition().Z()-point2_TK0.first.Z() );
774  double trackOutPosOutHitDr = diff.R();
775  double trackL = point2_TK0.second;
776  if (myverbose_==5) {
777  std::cout<<" propagted "<<point2_TK0.first<<" "<< point2_TK0.first.eta()<<" "<<point2_TK0.first.phi()<<std::endl;
778  std::cout<<" outerPosition() "<< pTrack->outerPosition() << " "<< pTrack->outerPosition().eta()<< " " << pTrack->outerPosition().phi()<< std::endl;
779  std::cout<<"diff " << diff << " diffR " <<diff.R()<<" diffR/L "<<diff.R()/point2_TK0.second <<std::endl;
780  }
781 
782  for(unsigned int ind=0;ind<recVtxs->size();ind++) {
783  if (!((*recVtxs)[ind].isFake())) {
784  reco::Vertex::trackRef_iterator vtxTrack = (*recVtxs)[ind].tracks_begin();
785  if( DeltaR(eta1,phi1, (*vtxTrack)->eta(),(*vtxTrack)->phi()) < 0.01 ) t_trackPVIdx ->push_back( ind );
786  else t_trackPVIdx ->push_back( -1 );
787  }
788  }
789 
790  // Fill the tree Branches here
791  t_trackPVIdx ->push_back( pVtxTkId );
792  t_trackP ->push_back( p1 );
793  t_trackPt ->push_back( pt1 );
794  t_trackEta ->push_back( eta1 );
795  t_trackPhi ->push_back( phi1 );
796  t_trackEcalEta ->push_back( etaEcal1 );
797  t_trackEcalPhi ->push_back( phiEcal1 );
798  t_trackHcalEta ->push_back( etaHcal1 );
799  t_trackHcalPhi ->push_back( phiHcal1 );
800  t_trackDxy ->push_back( dxy1 );
801  t_trackDz ->push_back( dz1 );
802  t_trackDxyBS ->push_back( dxybs1 );
803  t_trackDzBS ->push_back( dzbs1 );
804  t_trackDxyPV ->push_back( dxypv1 );
805  t_trackDzPV ->push_back( dzpv1 );
806  t_trackChiSq ->push_back( chisq1 );
807  t_trackNOuterHits ->push_back( nOuterHits );
808  t_NLayersCrossed ->push_back( nLayersCrossed );
809 
810  t_trackHitsTOB ->push_back( hitp.stripTOBLayersWithMeasurement() );
811  t_trackHitsTEC ->push_back( hitp.stripTECLayersWithMeasurement() );
812  t_trackHitInMissTOB ->push_back( hitp.stripTOBLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS) );
813  t_trackHitInMissTEC ->push_back( hitp.stripTECLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS) );
814  t_trackHitInMissTIB ->push_back( hitp.stripTIBLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS) );
815  t_trackHitInMissTID ->push_back( hitp.stripTIDLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS) );
816  t_trackHitInMissTIBTID ->push_back( hitp.stripTIBLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS) + hitp.stripTIDLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS) );
817 
818  t_trackHitOutMissTOB ->push_back( hitp.stripTOBLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS) );
819  t_trackHitOutMissTEC ->push_back( hitp.stripTECLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS) );
820  t_trackHitOutMissTIB ->push_back( hitp.stripTIBLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS) );
821  t_trackHitOutMissTID ->push_back( hitp.stripTIDLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS) );
822  t_trackHitOutMissTOBTEC ->push_back( hitp.stripTOBLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS) + hitp.stripTECLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS) );
823 
824  t_trackHitInMeasTOB ->push_back( hitp.stripTOBLayersWithMeasurement() );
825  t_trackHitInMeasTEC ->push_back( hitp.stripTECLayersWithMeasurement() );
826  t_trackHitInMeasTIB ->push_back( hitp.stripTIBLayersWithMeasurement() );
827  t_trackHitInMeasTID ->push_back( hitp.stripTIDLayersWithMeasurement() );
828  t_trackHitOutMeasTOB ->push_back( hitp.stripTOBLayersWithMeasurement() );
829  t_trackHitOutMeasTEC ->push_back( hitp.stripTECLayersWithMeasurement() );
830  t_trackHitOutMeasTIB ->push_back( hitp.stripTIBLayersWithMeasurement() );
831  t_trackHitOutMeasTID ->push_back( hitp.stripTIDLayersWithMeasurement() );
832  t_trackOutPosOutHitDr ->push_back( trackOutPosOutHitDr );
833  t_trackL ->push_back( trackL );
834 
835  t_maxNearP31x31 ->push_back( maxNearP31x31 );
836  t_maxNearP21x21 ->push_back( maxNearP21x21 );
837 
838  t_ecalSpike11x11 ->push_back( e11x11P.second );
839  t_e7x7 ->push_back( e7x7P.first );
840  t_e9x9 ->push_back( e9x9P.first );
841  t_e11x11 ->push_back( e11x11P.first );
842  t_e15x15 ->push_back( e15x15P.first );
843 
844  t_e7x7_10Sig ->push_back( e7x7_10SigP.first );
845  t_e9x9_10Sig ->push_back( e9x9_10SigP.first );
846  t_e11x11_10Sig ->push_back( e11x11_10SigP.first );
847  t_e15x15_10Sig ->push_back( e15x15_10SigP.first );
848  t_e7x7_15Sig ->push_back( e7x7_15SigP.first );
849  t_e9x9_15Sig ->push_back( e9x9_15SigP.first );
850  t_e11x11_15Sig ->push_back( e11x11_15SigP.first );
851  t_e15x15_15Sig ->push_back( e15x15_15SigP.first );
852  t_e7x7_20Sig ->push_back( e7x7_20SigP.first );
853  t_e9x9_20Sig ->push_back( e9x9_20SigP.first );
854  t_e11x11_20Sig ->push_back( e11x11_20SigP.first );
855  t_e15x15_20Sig ->push_back( e15x15_20SigP.first );
856  t_e7x7_25Sig ->push_back( e7x7_25SigP.first );
857  t_e9x9_25Sig ->push_back( e9x9_25SigP.first );
858  t_e11x11_25Sig ->push_back( e11x11_25SigP.first );
859  t_e15x15_25Sig ->push_back( e15x15_25SigP.first );
860  t_e7x7_30Sig ->push_back( e7x7_30SigP.first );
861  t_e9x9_30Sig ->push_back( e9x9_30SigP.first );
862  t_e11x11_30Sig ->push_back( e11x11_30SigP.first );
863  t_e15x15_30Sig ->push_back( e15x15_30SigP.first );
864 
865  if (doMC) {
866  t_esim7x7 ->push_back( simInfo7x7.eTotal );
867  t_esim9x9 ->push_back( simInfo9x9.eTotal );
868  t_esim11x11 ->push_back( simInfo11x11.eTotal );
869  t_esim15x15 ->push_back( simInfo15x15.eTotal );
870 
871  t_esim7x7Matched ->push_back( simInfo7x7.eMatched );
872  t_esim9x9Matched ->push_back( simInfo9x9.eMatched );
873  t_esim11x11Matched ->push_back( simInfo11x11.eMatched );
874  t_esim15x15Matched ->push_back( simInfo15x15.eMatched );
875 
876  t_esim7x7Rest ->push_back( simInfo7x7.eRest );
877  t_esim9x9Rest ->push_back( simInfo9x9.eRest );
878  t_esim11x11Rest ->push_back( simInfo11x11.eRest );
879  t_esim15x15Rest ->push_back( simInfo15x15.eRest );
880 
881  t_esim7x7Photon ->push_back( simInfo7x7.eGamma );
882  t_esim9x9Photon ->push_back( simInfo9x9.eGamma );
883  t_esim11x11Photon ->push_back( simInfo11x11.eGamma );
884  t_esim15x15Photon ->push_back( simInfo15x15.eGamma );
885 
886  t_esim7x7NeutHad ->push_back( simInfo7x7.eNeutralHad );
887  t_esim9x9NeutHad ->push_back( simInfo9x9.eNeutralHad );
888  t_esim11x11NeutHad ->push_back( simInfo11x11.eNeutralHad );
889  t_esim15x15NeutHad ->push_back( simInfo15x15.eNeutralHad );
890 
891  t_esim7x7CharHad ->push_back( simInfo7x7.eChargedHad );
892  t_esim9x9CharHad ->push_back( simInfo9x9.eChargedHad );
893  t_esim11x11CharHad ->push_back( simInfo11x11.eChargedHad );
894  t_esim15x15CharHad ->push_back( simInfo15x15.eChargedHad );
895 
896  t_trkEcalEne ->push_back( trkEcalEne );
897  t_simTrackP ->push_back( simTrackP );
898  t_esimPdgId ->push_back( simInfo11x11.pdgMatched );
899  }
900 
901  t_maxNearHcalP3x3 ->push_back( maxNearHcalP3x3 );
902  t_maxNearHcalP5x5 ->push_back( maxNearHcalP5x5 );
903  t_maxNearHcalP7x7 ->push_back( maxNearHcalP7x7 );
904 
905  t_h3x3 ->push_back( h3x3 );
906  t_h5x5 ->push_back( h5x5 );
907  t_h7x7 ->push_back( h7x7 );
908  t_h3x3Sig ->push_back( h3x3Sig );
909  t_h5x5Sig ->push_back( h5x5Sig );
910  t_h7x7Sig ->push_back( h7x7Sig );
911 
912  t_infoHcal ->push_back( trkDetItr->okHCAL );
913  if (doMC) {
914  t_trkHcalEne ->push_back( hcalScale*trkHcalEne );
915 
916  t_hsim3x3 ->push_back( hcalScale*hsimInfo3x3.eTotal );
917  t_hsim5x5 ->push_back( hcalScale*hsimInfo5x5.eTotal );
918  t_hsim7x7 ->push_back( hcalScale*hsimInfo7x7.eTotal );
919 
920  t_hsim3x3Matched ->push_back( hcalScale*hsimInfo3x3.eMatched );
921  t_hsim5x5Matched ->push_back( hcalScale*hsimInfo5x5.eMatched );
922  t_hsim7x7Matched ->push_back( hcalScale*hsimInfo7x7.eMatched );
923 
924  t_hsim3x3Rest ->push_back( hcalScale*hsimInfo3x3.eRest );
925  t_hsim5x5Rest ->push_back( hcalScale*hsimInfo5x5.eRest );
926  t_hsim7x7Rest ->push_back( hcalScale*hsimInfo7x7.eRest );
927 
928  t_hsim3x3Photon ->push_back( hcalScale*hsimInfo3x3.eGamma );
929  t_hsim5x5Photon ->push_back( hcalScale*hsimInfo5x5.eGamma );
930  t_hsim7x7Photon ->push_back( hcalScale*hsimInfo7x7.eGamma );
931 
932  t_hsim3x3NeutHad ->push_back( hcalScale*hsimInfo3x3.eNeutralHad );
933  t_hsim5x5NeutHad ->push_back( hcalScale*hsimInfo5x5.eNeutralHad );
934  t_hsim7x7NeutHad ->push_back( hcalScale*hsimInfo7x7.eNeutralHad );
935 
936  t_hsim3x3CharHad ->push_back( hcalScale*hsimInfo3x3.eChargedHad );
937  t_hsim5x5CharHad ->push_back( hcalScale*hsimInfo5x5.eChargedHad );
938  t_hsim7x7CharHad ->push_back( hcalScale*hsimInfo7x7.eChargedHad );
939  }
940 
941  } // if loosely isolated track
942  } // check p1/eta
943  } // loop over track collection
944 
945  if(haveIsoTrack) tree->Fill();
946 }
947 
948 // ----- method called once each job just before starting event loop ----
949 void IsolatedTracksNxN::beginJob() {
950 
951  nEventProc=0;
952 
953  // double tempgen_TH[21] = { 1.0, 2.0, 3.0, 4.0, 5.0,
954  double tempgen_TH[16] = { 0.0, 1.0, 2.0, 3.0, 4.0,
955  5.0, 6.0, 7.0, 9.0, 11.0,
956  15.0, 20.0, 30.0, 50.0, 75.0, 100.0};
957 
958  for(int i=0; i<16; i++) genPartPBins[i] = tempgen_TH[i];
959 
960  double tempgen_Eta[4] = {0.0, 1.131, 1.653, 2.172};
961 
962  for(int i=0; i<4; i++) genPartEtaBins[i] = tempgen_Eta[i];
963 
964  BookHistograms();
965 }
966 
967 // ----- method called once each job just after ending the event loop ----
968 void IsolatedTracksNxN::endJob() {
969 
970  if(L1TriggerAlgoInfo_) {
971  std::map< std::pair<unsigned int,std::string>, int>::iterator itr;
972  for(itr=l1AlgoMap.begin(); itr!=l1AlgoMap.end(); itr++) {
973  std::cout << " ****endjob**** " << (itr->first).first <<" "
974  <<(itr->first).second <<" "<<itr->second
975  << std::endl;
976  int ibin = (itr->first).first;
977  TString name( (itr->first).second );
978  h_L1AlgoNames->GetXaxis()->SetBinLabel(ibin+1,name);
979  }
980  std::cout << "Number of Events Processed " << nEventProc << std::endl;
981  }
982 
983 }
984 
985 //========================================================================================================
986 
987 void IsolatedTracksNxN::clearTreeVectors() {
988 
989  t_PVx ->clear();
990  t_PVy ->clear();
991  t_PVz ->clear();
992  t_PVisValid ->clear();
993  t_PVndof ->clear();
994  t_PVNTracks ->clear();
995  t_PVNTracksWt ->clear();
996  t_PVTracksSumPt ->clear();
997  t_PVTracksSumPtWt ->clear();
998  t_PVNTracksHP ->clear();
999  t_PVNTracksHPWt ->clear();
1000  t_PVTracksSumPtHP ->clear();
1001  t_PVTracksSumPtHPWt ->clear();
1002 
1003  for(int i=0; i<128; i++) t_L1Decision[i]=0;
1004  t_L1AlgoNames ->clear();
1005  t_L1PreScale ->clear();
1006 
1007  t_L1CenJetPt ->clear();
1008  t_L1CenJetEta ->clear();
1009  t_L1CenJetPhi ->clear();
1010  t_L1FwdJetPt ->clear();
1011  t_L1FwdJetEta ->clear();
1012  t_L1FwdJetPhi ->clear();
1013  t_L1TauJetPt ->clear();
1014  t_L1TauJetEta ->clear();
1015  t_L1TauJetPhi ->clear();
1016  t_L1MuonPt ->clear();
1017  t_L1MuonEta ->clear();
1018  t_L1MuonPhi ->clear();
1019  t_L1IsoEMPt ->clear();
1020  t_L1IsoEMEta ->clear();
1021  t_L1IsoEMPhi ->clear();
1022  t_L1NonIsoEMPt ->clear();
1023  t_L1NonIsoEMEta ->clear();
1024  t_L1NonIsoEMPhi ->clear();
1025  t_L1METPt ->clear();
1026  t_L1METEta ->clear();
1027  t_L1METPhi ->clear();
1028 
1029  t_jetPt ->clear();
1030  t_jetEta ->clear();
1031  t_jetPhi ->clear();
1032  t_nTrksJetCalo ->clear();
1033  t_nTrksJetVtx ->clear();
1034 
1035  t_trackPAll ->clear();
1036  t_trackEtaAll ->clear();
1037  t_trackPhiAll ->clear();
1038  t_trackPdgIdAll ->clear();
1039  t_trackPtAll ->clear();
1040  t_trackDxyAll ->clear();
1041  t_trackDzAll ->clear();
1042  t_trackDxyPVAll ->clear();
1043  t_trackDzPVAll ->clear();
1044  t_trackChiSqAll ->clear();
1045 
1046  t_trackP ->clear();
1047  t_trackPt ->clear();
1048  t_trackEta ->clear();
1049  t_trackPhi ->clear();
1050  t_trackEcalEta ->clear();
1051  t_trackEcalPhi ->clear();
1052  t_trackHcalEta ->clear();
1053  t_trackHcalPhi ->clear();
1054  t_NLayersCrossed ->clear();
1055  t_trackNOuterHits ->clear();
1056  t_trackDxy ->clear();
1057  t_trackDxyBS ->clear();
1058  t_trackDz ->clear();
1059  t_trackDzBS ->clear();
1060  t_trackDxyPV ->clear();
1061  t_trackDzPV ->clear();
1062  t_trackChiSq ->clear();
1063  t_trackPVIdx ->clear();
1064  t_trackHitsTOB ->clear();
1065  t_trackHitsTEC ->clear();
1066  t_trackHitInMissTOB ->clear();
1067  t_trackHitInMissTEC ->clear();
1068  t_trackHitInMissTIB ->clear();
1069  t_trackHitInMissTID ->clear();
1070  t_trackHitInMissTIBTID ->clear();
1071  t_trackHitOutMissTOB ->clear();
1072  t_trackHitOutMissTEC ->clear();
1073  t_trackHitOutMissTIB ->clear();
1074  t_trackHitOutMissTID ->clear();
1075  t_trackHitOutMissTOBTEC ->clear();
1076  t_trackHitInMeasTOB ->clear();
1077  t_trackHitInMeasTEC ->clear();
1078  t_trackHitInMeasTIB ->clear();
1079  t_trackHitInMeasTID ->clear();
1080  t_trackHitOutMeasTOB ->clear();
1081  t_trackHitOutMeasTEC ->clear();
1082  t_trackHitOutMeasTIB ->clear();
1083  t_trackHitOutMeasTID ->clear();
1084  t_trackOutPosOutHitDr ->clear();
1085  t_trackL ->clear();
1086 
1087  t_maxNearP31x31 ->clear();
1088  t_maxNearP21x21 ->clear();
1089 
1090  t_ecalSpike11x11 ->clear();
1091  t_e7x7 ->clear();
1092  t_e9x9 ->clear();
1093  t_e11x11 ->clear();
1094  t_e15x15 ->clear();
1095 
1096  t_e7x7_10Sig ->clear();
1097  t_e9x9_10Sig ->clear();
1098  t_e11x11_10Sig ->clear();
1099  t_e15x15_10Sig ->clear();
1100  t_e7x7_15Sig ->clear();
1101  t_e9x9_15Sig ->clear();
1102  t_e11x11_15Sig ->clear();
1103  t_e15x15_15Sig ->clear();
1104  t_e7x7_20Sig ->clear();
1105  t_e9x9_20Sig ->clear();
1106  t_e11x11_20Sig ->clear();
1107  t_e15x15_20Sig ->clear();
1108  t_e7x7_25Sig ->clear();
1109  t_e9x9_25Sig ->clear();
1110  t_e11x11_25Sig ->clear();
1111  t_e15x15_25Sig ->clear();
1112  t_e7x7_30Sig ->clear();
1113  t_e9x9_30Sig ->clear();
1114  t_e11x11_30Sig ->clear();
1115  t_e15x15_30Sig ->clear();
1116 
1117  if (doMC) {
1118  t_simTrackP ->clear();
1119  t_esimPdgId ->clear();
1120  t_trkEcalEne ->clear();
1121 
1122  t_esim7x7 ->clear();
1123  t_esim9x9 ->clear();
1124  t_esim11x11 ->clear();
1125  t_esim15x15 ->clear();
1126 
1127  t_esim7x7Matched ->clear();
1128  t_esim9x9Matched ->clear();
1129  t_esim11x11Matched ->clear();
1130  t_esim15x15Matched ->clear();
1131 
1132  t_esim7x7Rest ->clear();
1133  t_esim9x9Rest ->clear();
1134  t_esim11x11Rest ->clear();
1135  t_esim15x15Rest ->clear();
1136 
1137  t_esim7x7Photon ->clear();
1138  t_esim9x9Photon ->clear();
1139  t_esim11x11Photon ->clear();
1140  t_esim15x15Photon ->clear();
1141 
1142  t_esim7x7NeutHad ->clear();
1143  t_esim9x9NeutHad ->clear();
1144  t_esim11x11NeutHad ->clear();
1145  t_esim15x15NeutHad ->clear();
1146 
1147  t_esim7x7CharHad ->clear();
1148  t_esim9x9CharHad ->clear();
1149  t_esim11x11CharHad ->clear();
1150  t_esim15x15CharHad ->clear();
1151  }
1152 
1153  t_maxNearHcalP3x3 ->clear();
1154  t_maxNearHcalP5x5 ->clear();
1155  t_maxNearHcalP7x7 ->clear();
1156 
1157  t_h3x3 ->clear();
1158  t_h5x5 ->clear();
1159  t_h7x7 ->clear();
1160  t_h3x3Sig ->clear();
1161  t_h5x5Sig ->clear();
1162  t_h7x7Sig ->clear();
1163 
1164  t_infoHcal ->clear();
1165 
1166  if (doMC) {
1167  t_trkHcalEne ->clear();
1168 
1169  t_hsim3x3 ->clear();
1170  t_hsim5x5 ->clear();
1171  t_hsim7x7 ->clear();
1172  t_hsim3x3Matched ->clear();
1173  t_hsim5x5Matched ->clear();
1174  t_hsim7x7Matched ->clear();
1175  t_hsim3x3Rest ->clear();
1176  t_hsim5x5Rest ->clear();
1177  t_hsim7x7Rest ->clear();
1178  t_hsim3x3Photon ->clear();
1179  t_hsim5x5Photon ->clear();
1180  t_hsim7x7Photon ->clear();
1181  t_hsim3x3NeutHad ->clear();
1182  t_hsim5x5NeutHad ->clear();
1183  t_hsim7x7NeutHad ->clear();
1184  t_hsim3x3CharHad ->clear();
1185  t_hsim5x5CharHad ->clear();
1186  t_hsim7x7CharHad ->clear();
1187  }
1188 }
1189 
1190 void IsolatedTracksNxN::BookHistograms(){
1191 
1192  char hname[100], htit[100];
1193 
1194  TFileDirectory dir = fs->mkdir("nearMaxTrackP");
1195 
1196  for(unsigned int ieta=0; ieta<NEtaBins; ieta++) {
1197  double lowEta=-5.0, highEta= 5.0;
1198  lowEta = genPartEtaBins[ieta];
1199  highEta = genPartEtaBins[ieta+1];
1200 
1201  for(unsigned int ipt=0; ipt<NPBins; ipt++) {
1202  double lowP=0.0, highP=300.0;
1203  lowP = genPartPBins[ipt];
1204  highP = genPartPBins[ipt+1];
1205  sprintf(hname, "h_maxNearP31x31_ptBin%i_etaBin%i",ipt, ieta);
1206  sprintf(htit, "maxNearP in 31x31 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP );
1207  h_maxNearP31x31[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
1208  h_maxNearP31x31[ipt][ieta] ->Sumw2();
1209  sprintf(hname, "h_maxNearP25x25_ptBin%i_etaBin%i",ipt, ieta);
1210  sprintf(htit, "maxNearP in 25x25 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP );
1211  h_maxNearP25x25[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
1212  h_maxNearP25x25[ipt][ieta] ->Sumw2();
1213  sprintf(hname, "h_maxNearP21x21_ptBin%i_etaBin%i",ipt, ieta);
1214  sprintf(htit, "maxNearP in 21x21 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP );
1215  h_maxNearP21x21[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
1216  h_maxNearP21x21[ipt][ieta] ->Sumw2();
1217  sprintf(hname, "h_maxNearP15x15_ptBin%i_etaBin%i",ipt, ieta);
1218  sprintf(htit, "maxNearP in 15x15 (%3.2f<|#eta|<%3.2f), (%2.0f<trkP<%3.0f)", lowEta, highEta, lowP, highP );
1219  h_maxNearP15x15[ipt][ieta] = dir.make<TH1F>(hname, htit, 220, -2.0, 20.0);
1220  h_maxNearP15x15[ipt][ieta] ->Sumw2();
1221  }
1222  }
1223 
1224  h_L1AlgoNames = fs->make<TH1I>("h_L1AlgoNames", "h_L1AlgoNames:Bin Labels", 128, -0.5, 127.5);
1225 
1226  // Reconstructed Tracks
1227 
1228  h_PVTracksWt = fs->make<TH1F>("h_PVTracksWt", "h_PVTracksWt", 600, -0.1, 1.1);
1229 
1230  h_nTracks = fs->make<TH1F>("h_nTracks", "h_nTracks", 1000, -0.5, 999.5);
1231 
1232  sprintf(hname, "h_recEtaPt_0");
1233  sprintf(htit, "h_recEtaPt (all tracks Eta vs pT)");
1234  h_recEtaPt_0 = fs->make<TH2F>(hname, htit, 30, -3.0,3.0, 15, genPartPBins);
1235 
1236  sprintf(hname, "h_recEtaP_0");
1237  sprintf(htit, "h_recEtaP (all tracks Eta vs pT)");
1238  h_recEtaP_0 = fs->make<TH2F>(hname, htit, 30, -3.0,3.0, 15, genPartPBins);
1239 
1240  h_recPt_0 = fs->make<TH1F>("h_recPt_0", "Pt (all tracks)", 15, genPartPBins);
1241  h_recP_0 = fs->make<TH1F>("h_recP_0", "P (all tracks)", 15, genPartPBins);
1242  h_recEta_0 = fs->make<TH1F>("h_recEta_0", "Eta (all tracks)", 60, -3.0, 3.0);
1243  h_recPhi_0 = fs->make<TH1F>("h_recPhi_0", "Phi (all tracks)", 100, -3.2, 3.2);
1244  //-------------------------
1245  sprintf(hname, "h_recEtaPt_1");
1246  sprintf(htit, "h_recEtaPt (all good tracks Eta vs pT)");
1247  h_recEtaPt_1 = fs->make<TH2F>(hname, htit, 30, -3.0,3.0, 15, genPartPBins);
1248 
1249  sprintf(hname, "h_recEtaP_1");
1250  sprintf(htit, "h_recEtaP (all good tracks Eta vs pT)");
1251  h_recEtaP_1 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
1252 
1253  h_recPt_1 = fs->make<TH1F>("h_recPt_1", "Pt (all good tracks)", 15, genPartPBins);
1254  h_recP_1 = fs->make<TH1F>("h_recP_1", "P (all good tracks)", 15, genPartPBins);
1255  h_recEta_1 = fs->make<TH1F>("h_recEta_1", "Eta (all good tracks)", 60, -3.0, 3.0);
1256  h_recPhi_1 = fs->make<TH1F>("h_recPhi_1", "Phi (all good tracks)", 100, -3.2, 3.2);
1257  //-------------------------
1258  sprintf(hname, "h_recEtaPt_2");
1259  sprintf(htit, "h_recEtaPt (charge isolation Eta vs pT)");
1260  h_recEtaPt_2 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
1261 
1262  sprintf(hname, "h_recEtaP_2");
1263  sprintf(htit, "h_recEtaP (charge isolation Eta vs pT)");
1264  h_recEtaP_2 = fs->make<TH2F>(hname, htit, 30, -3.0, 3.0, 15, genPartPBins);
1265 
1266  h_recPt_2 = fs->make<TH1F>("h_recPt_2", "Pt (charge isolation)", 15, genPartPBins);
1267  h_recP_2 = fs->make<TH1F>("h_recP_2", "P (charge isolation)", 15, genPartPBins);
1268  h_recEta_2 = fs->make<TH1F>("h_recEta_2", "Eta (charge isolation)", 60, -3.0, 3.0);
1269  h_recPhi_2 = fs->make<TH1F>("h_recPhi_2", "Phi (charge isolation)", 100, -3.2, 3.2);
1270 
1271 
1272  tree = fs->make<TTree>("tree", "tree");
1273  tree->SetAutoSave(10000);
1274 
1275 
1276  tree->Branch("t_EvtNo" ,&t_EvtNo ,"t_EvtNo/I");
1277  tree->Branch("t_RunNo" ,&t_RunNo ,"t_RunNo/I");
1278  tree->Branch("t_Lumi" ,&t_Lumi ,"t_Lumi/I");
1279  tree->Branch("t_Bunch" ,&t_Bunch ,"t_Bunch/I");
1280 
1281 
1282  t_PVx = new std::vector<double>();
1283  t_PVy = new std::vector<double>();
1284  t_PVz = new std::vector<double>();
1285  t_PVisValid = new std::vector<int>();
1286  t_PVndof = new std::vector<int>();
1287  t_PVNTracks = new std::vector<int>();
1288  t_PVNTracksWt = new std::vector<int>();
1289  t_PVTracksSumPt = new std::vector<double>();
1290  t_PVTracksSumPtWt = new std::vector<double>();
1291  t_PVNTracksHP = new std::vector<int>();
1292  t_PVNTracksHPWt = new std::vector<int>();
1293  t_PVTracksSumPtHP = new std::vector<double>();
1294  t_PVTracksSumPtHPWt = new std::vector<double>();
1295 
1296  tree->Branch("PVx" ,"vector<double>" ,&t_PVx);
1297  tree->Branch("PVy" ,"vector<double>" ,&t_PVy);
1298  tree->Branch("PVz" ,"vector<double>" ,&t_PVz);
1299  tree->Branch("PVisValid" ,"vector<int>" ,&t_PVisValid);
1300  tree->Branch("PVndof" ,"vector<int>" ,&t_PVndof);
1301  tree->Branch("PVNTracks" ,"vector<int>" ,&t_PVNTracks);
1302  tree->Branch("PVNTracksWt" ,"vector<int>" ,&t_PVNTracksWt);
1303  tree->Branch("t_PVTracksSumPt" ,"vector<double>" ,&t_PVTracksSumPt);
1304  tree->Branch("t_PVTracksSumPtWt" ,"vector<double>" ,&t_PVTracksSumPtWt);
1305  tree->Branch("PVNTracksHP" ,"vector<int>" ,&t_PVNTracksHP);
1306  tree->Branch("PVNTracksHPWt" ,"vector<int>" ,&t_PVNTracksHPWt);
1307  tree->Branch("t_PVTracksSumPtHP" ,"vector<double>" ,&t_PVTracksSumPtHP);
1308  tree->Branch("t_PVTracksSumPtHPWt" ,"vector<double>" ,&t_PVTracksSumPtHPWt);
1309 
1310  //----- L1Trigger information
1311  for(int i=0; i<128; i++) t_L1Decision[i]=0;
1312  t_L1AlgoNames = new std::vector<std::string>();
1313  t_L1PreScale = new std::vector<int>();
1314  t_L1CenJetPt = new std::vector<double>();
1315  t_L1CenJetEta = new std::vector<double>();
1316  t_L1CenJetPhi = new std::vector<double>();
1317  t_L1FwdJetPt = new std::vector<double>();
1318  t_L1FwdJetEta = new std::vector<double>();
1319  t_L1FwdJetPhi = new std::vector<double>();
1320  t_L1TauJetPt = new std::vector<double>();
1321  t_L1TauJetEta = new std::vector<double>();
1322  t_L1TauJetPhi = new std::vector<double>();
1323  t_L1MuonPt = new std::vector<double>();
1324  t_L1MuonEta = new std::vector<double>();
1325  t_L1MuonPhi = new std::vector<double>();
1326  t_L1IsoEMPt = new std::vector<double>();
1327  t_L1IsoEMEta = new std::vector<double>();
1328  t_L1IsoEMPhi = new std::vector<double>();
1329  t_L1NonIsoEMPt = new std::vector<double>();
1330  t_L1NonIsoEMEta = new std::vector<double>();
1331  t_L1NonIsoEMPhi = new std::vector<double>();
1332  t_L1METPt = new std::vector<double>();
1333  t_L1METEta = new std::vector<double>();
1334  t_L1METPhi = new std::vector<double>();
1335 
1336  tree->Branch("t_L1Decision", t_L1Decision, "t_L1Decision[128]/I");
1337  tree->Branch("t_L1AlgoNames", "vector<string>", &t_L1AlgoNames);
1338  tree->Branch("t_L1PreScale", "vector<int>", &t_L1PreScale);
1339  tree->Branch("t_L1CenJetPt", "vector<double>", &t_L1CenJetPt);
1340  tree->Branch("t_L1CenJetEta", "vector<double>", &t_L1CenJetEta);
1341  tree->Branch("t_L1CenJetPhi", "vector<double>", &t_L1CenJetPhi);
1342  tree->Branch("t_L1FwdJetPt", "vector<double>", &t_L1FwdJetPt);
1343  tree->Branch("t_L1FwdJetEta", "vector<double>", &t_L1FwdJetEta);
1344  tree->Branch("t_L1FwdJetPhi", "vector<double>", &t_L1FwdJetPhi);
1345  tree->Branch("t_L1TauJetPt", "vector<double>", &t_L1TauJetPt);
1346  tree->Branch("t_L1TauJetEta", "vector<double>", &t_L1TauJetEta);
1347  tree->Branch("t_L1TauJetPhi", "vector<double>", &t_L1TauJetPhi);
1348  tree->Branch("t_L1MuonPt", "vector<double>", &t_L1MuonPt);
1349  tree->Branch("t_L1MuonEta", "vector<double>", &t_L1MuonEta);
1350  tree->Branch("t_L1MuonPhi", "vector<double>", &t_L1MuonPhi);
1351  tree->Branch("t_L1IsoEMPt", "vector<double>", &t_L1IsoEMPt);
1352  tree->Branch("t_L1IsoEMEta", "vector<double>", &t_L1IsoEMEta);
1353  tree->Branch("t_L1IsoEMPhi", "vector<double>", &t_L1IsoEMPhi);
1354  tree->Branch("t_L1NonIsoEMPt", "vector<double>", &t_L1NonIsoEMPt);
1355  tree->Branch("t_L1NonIsoEMEta", "vector<double>", &t_L1NonIsoEMEta);
1356  tree->Branch("t_L1NonIsoEMPhi", "vector<double>", &t_L1NonIsoEMPhi);
1357  tree->Branch("t_L1METPt", "vector<double>", &t_L1METPt);
1358  tree->Branch("t_L1METEta", "vector<double>", &t_L1METEta);
1359  tree->Branch("t_L1METPhi", "vector<double>", &t_L1METPhi);
1360 
1361  t_jetPt = new std::vector<double>();
1362  t_jetEta = new std::vector<double>();
1363  t_jetPhi = new std::vector<double>();
1364  t_nTrksJetCalo = new std::vector<double>();
1365  t_nTrksJetVtx = new std::vector<double>();
1366  tree->Branch("t_jetPt", "vector<double>",&t_jetPt);
1367  tree->Branch("t_jetEta", "vector<double>",&t_jetEta);
1368  tree->Branch("t_jetPhi", "vector<double>",&t_jetPhi);
1369  tree->Branch("t_nTrksJetCalo", "vector<double>",&t_nTrksJetCalo);
1370  tree->Branch("t_nTrksJetVtx", "vector<double>",&t_nTrksJetVtx);
1371 
1372  t_trackPAll = new std::vector<double>();
1373  t_trackEtaAll = new std::vector<double>();
1374  t_trackPhiAll = new std::vector<double>();
1375  t_trackPdgIdAll = new std::vector<double>();
1376  t_trackPtAll = new std::vector<double>();
1377  t_trackDxyAll = new std::vector<double>();
1378  t_trackDzAll = new std::vector<double>();
1379  t_trackDxyPVAll = new std::vector<double>();
1380  t_trackDzPVAll = new std::vector<double>();
1381  t_trackChiSqAll = new std::vector<double>();
1382  tree->Branch("t_trackPAll", "vector<double>", &t_trackPAll );
1383  tree->Branch("t_trackPhiAll", "vector<double>", &t_trackPhiAll );
1384  tree->Branch("t_trackEtaAll", "vector<double>", &t_trackEtaAll );
1385  tree->Branch("t_trackPtAll", "vector<double>", &t_trackPtAll );
1386  tree->Branch("t_trackDxyAll", "vector<double>", &t_trackDxyAll );
1387  tree->Branch("t_trackDzAll", "vector<double>", &t_trackDzAll );
1388  tree->Branch("t_trackDxyPVAll", "vector<double>", &t_trackDxyPVAll );
1389  tree->Branch("t_trackDzPVAll", "vector<double>", &t_trackDzPVAll );
1390  tree->Branch("t_trackChiSqAll", "vector<double>", &t_trackChiSqAll );
1391  //tree->Branch("t_trackPdgIdAll", "vector<double>", &t_trackPdgIdAll);
1392 
1393  t_trackP = new std::vector<double>();
1394  t_trackPt = new std::vector<double>();
1395  t_trackEta = new std::vector<double>();
1396  t_trackPhi = new std::vector<double>();
1397  t_trackEcalEta = new std::vector<double>();
1398  t_trackEcalPhi = new std::vector<double>();
1399  t_trackHcalEta = new std::vector<double>();
1400  t_trackHcalPhi = new std::vector<double>();
1401  t_trackNOuterHits = new std::vector<int>();
1402  t_NLayersCrossed = new std::vector<int>();
1403  t_trackDxy = new std::vector<double>();
1404  t_trackDxyBS = new std::vector<double>();
1405  t_trackDz = new std::vector<double>();
1406  t_trackDzBS = new std::vector<double>();
1407  t_trackDxyPV = new std::vector<double>();
1408  t_trackDzPV = new std::vector<double>();
1409  t_trackPVIdx = new std::vector<int>();
1410  t_trackChiSq = new std::vector<double>();
1411  t_trackHitsTOB = new std::vector<int>();
1412  t_trackHitsTEC = new std::vector<int>();
1413  t_trackHitInMissTOB = new std::vector<int>();
1414  t_trackHitInMissTEC = new std::vector<int>();
1415  t_trackHitInMissTIB = new std::vector<int>();
1416  t_trackHitInMissTID = new std::vector<int>();
1417  t_trackHitOutMissTOB = new std::vector<int>();
1418  t_trackHitOutMissTEC = new std::vector<int>();
1419  t_trackHitOutMissTIB = new std::vector<int>();
1420  t_trackHitOutMissTID = new std::vector<int>();
1421  t_trackHitInMissTIBTID = new std::vector<int>();
1422  t_trackHitOutMissTOB = new std::vector<int>();
1423  t_trackHitOutMissTEC = new std::vector<int>();
1424  t_trackHitOutMissTIB = new std::vector<int>();
1425  t_trackHitOutMissTID = new std::vector<int>();
1426  t_trackHitOutMissTOBTEC= new std::vector<int>();
1427  t_trackHitInMeasTOB = new std::vector<int>();
1428  t_trackHitInMeasTEC = new std::vector<int>();
1429  t_trackHitInMeasTIB = new std::vector<int>();
1430  t_trackHitInMeasTID = new std::vector<int>();
1431  t_trackHitOutMeasTOB = new std::vector<int>();
1432  t_trackHitOutMeasTEC = new std::vector<int>();
1433  t_trackHitOutMeasTIB = new std::vector<int>();
1434  t_trackHitOutMeasTID = new std::vector<int>();
1435  t_trackOutPosOutHitDr = new std::vector<double>();
1436  t_trackL = new std::vector<double>();
1437 
1438  tree->Branch("t_trackP", "vector<double>", &t_trackP );
1439  tree->Branch("t_trackPt", "vector<double>", &t_trackPt );
1440  tree->Branch("t_trackEta", "vector<double>", &t_trackEta );
1441  tree->Branch("t_trackPhi", "vector<double>", &t_trackPhi );
1442  tree->Branch("t_trackEcalEta", "vector<double>", &t_trackEcalEta );
1443  tree->Branch("t_trackEcalPhi", "vector<double>", &t_trackEcalPhi );
1444  tree->Branch("t_trackHcalEta", "vector<double>", &t_trackHcalEta );
1445  tree->Branch("t_trackHcalPhi", "vector<double>", &t_trackHcalPhi );
1446 
1447  tree->Branch("t_trackNOuterHits", "vector<int>", &t_trackNOuterHits );
1448  tree->Branch("t_NLayersCrossed", "vector<int>", &t_NLayersCrossed );
1449  tree->Branch("t_trackHitsTOB", "vector<int>", &t_trackHitsTOB );
1450  tree->Branch("t_trackHitsTEC", "vector<int>", &t_trackHitsTEC );
1451  tree->Branch("t_trackHitInMissTOB", "vector<int>", &t_trackHitInMissTOB );
1452  tree->Branch("t_trackHitInMissTEC", "vector<int>", &t_trackHitInMissTEC );
1453  tree->Branch("t_trackHitInMissTIB", "vector<int>", &t_trackHitInMissTIB );
1454  tree->Branch("t_trackHitInMissTID", "vector<int>", &t_trackHitInMissTID );
1455  tree->Branch("t_trackHitInMissTIBTID", "vector<int>", &t_trackHitInMissTIBTID );
1456  tree->Branch("t_trackHitOutMissTOB", "vector<int>", &t_trackHitOutMissTOB);
1457  tree->Branch("t_trackHitOutMissTEC", "vector<int>", &t_trackHitOutMissTEC);
1458  tree->Branch("t_trackHitOutMissTIB", "vector<int>", &t_trackHitOutMissTIB);
1459  tree->Branch("t_trackHitOutMissTID", "vector<int>", &t_trackHitOutMissTID);
1460  tree->Branch("t_trackHitOutMissTOBTEC","vector<int>", &t_trackHitOutMissTOBTEC);
1461  tree->Branch("t_trackHitInMeasTOB", "vector<int>", &t_trackHitInMeasTOB );
1462  tree->Branch("t_trackHitInMeasTEC", "vector<int>", &t_trackHitInMeasTEC );
1463  tree->Branch("t_trackHitInMeasTIB", "vector<int>", &t_trackHitInMeasTIB );
1464  tree->Branch("t_trackHitInMeasTID", "vector<int>", &t_trackHitInMeasTID );
1465  tree->Branch("t_trackHitOutMeasTOB", "vector<int>", &t_trackHitOutMeasTOB);
1466  tree->Branch("t_trackHitOutMeasTEC", "vector<int>", &t_trackHitOutMeasTEC);
1467  tree->Branch("t_trackHitOutMeasTIB", "vector<int>", &t_trackHitOutMeasTIB);
1468  tree->Branch("t_trackHitOutMeasTID", "vector<int>", &t_trackHitOutMeasTID);
1469  tree->Branch("t_trackOutPosOutHitDr", "vector<double>", &t_trackOutPosOutHitDr);
1470  tree->Branch("t_trackL", "vector<double>", &t_trackL);
1471 
1472  tree->Branch("t_trackDxy", "vector<double>", &t_trackDxy );
1473  tree->Branch("t_trackDxyBS", "vector<double>", &t_trackDxyBS );
1474  tree->Branch("t_trackDz", "vector<double>", &t_trackDz );
1475  tree->Branch("t_trackDzBS", "vector<double>", &t_trackDzBS );
1476  tree->Branch("t_trackDxyPV", "vector<double>", &t_trackDxyPV );
1477  tree->Branch("t_trackDzPV", "vector<double>", &t_trackDzPV );
1478  tree->Branch("t_trackChiSq", "vector<double>", &t_trackChiSq );
1479  tree->Branch("t_trackPVIdx", "vector<int>", &t_trackPVIdx );
1480 
1481  t_maxNearP31x31 = new std::vector<double>();
1482  t_maxNearP21x21 = new std::vector<double>();
1483 
1484  tree->Branch("t_maxNearP31x31", "vector<double>", &t_maxNearP31x31);
1485  tree->Branch("t_maxNearP21x21", "vector<double>", &t_maxNearP21x21);
1486 
1487  t_ecalSpike11x11 = new std::vector<int>();
1488  t_e7x7 = new std::vector<double>();
1489  t_e9x9 = new std::vector<double>();
1490  t_e11x11 = new std::vector<double>();
1491  t_e15x15 = new std::vector<double>();
1492 
1493  tree->Branch("t_ecalSpike11x11", "vector<int>", &t_ecalSpike11x11);
1494  tree->Branch("t_e7x7", "vector<double>", &t_e7x7);
1495  tree->Branch("t_e9x9", "vector<double>", &t_e9x9);
1496  tree->Branch("t_e11x11", "vector<double>", &t_e11x11);
1497  tree->Branch("t_e15x15", "vector<double>", &t_e15x15);
1498 
1499  t_e7x7_10Sig = new std::vector<double>();
1500  t_e9x9_10Sig = new std::vector<double>();
1501  t_e11x11_10Sig = new std::vector<double>();
1502  t_e15x15_10Sig = new std::vector<double>();
1503  t_e7x7_15Sig = new std::vector<double>();
1504  t_e9x9_15Sig = new std::vector<double>();
1505  t_e11x11_15Sig = new std::vector<double>();
1506  t_e15x15_15Sig = new std::vector<double>();
1507  t_e7x7_20Sig = new std::vector<double>();
1508  t_e9x9_20Sig = new std::vector<double>();
1509  t_e11x11_20Sig = new std::vector<double>();
1510  t_e15x15_20Sig = new std::vector<double>();
1511  t_e7x7_25Sig = new std::vector<double>();
1512  t_e9x9_25Sig = new std::vector<double>();
1513  t_e11x11_25Sig = new std::vector<double>();
1514  t_e15x15_25Sig = new std::vector<double>();
1515  t_e7x7_30Sig = new std::vector<double>();
1516  t_e9x9_30Sig = new std::vector<double>();
1517  t_e11x11_30Sig = new std::vector<double>();
1518  t_e15x15_30Sig = new std::vector<double>();
1519 
1520  tree->Branch("t_e7x7_10Sig" ,"vector<double>", &t_e7x7_10Sig);
1521  tree->Branch("t_e9x9_10Sig" ,"vector<double>", &t_e9x9_10Sig);
1522  tree->Branch("t_e11x11_10Sig" ,"vector<double>", &t_e11x11_10Sig);
1523  tree->Branch("t_e15x15_10Sig" ,"vector<double>", &t_e15x15_10Sig);
1524  tree->Branch("t_e7x7_15Sig" ,"vector<double>", &t_e7x7_15Sig);
1525  tree->Branch("t_e9x9_15Sig" ,"vector<double>", &t_e9x9_15Sig);
1526  tree->Branch("t_e11x11_15Sig" ,"vector<double>", &t_e11x11_15Sig);
1527  tree->Branch("t_e15x15_15Sig" ,"vector<double>", &t_e15x15_15Sig);
1528  tree->Branch("t_e7x7_20Sig" ,"vector<double>", &t_e7x7_20Sig);
1529  tree->Branch("t_e9x9_20Sig" ,"vector<double>", &t_e9x9_20Sig);
1530  tree->Branch("t_e11x11_20Sig" ,"vector<double>", &t_e11x11_20Sig);
1531  tree->Branch("t_e15x15_20Sig" ,"vector<double>", &t_e15x15_20Sig);
1532  tree->Branch("t_e7x7_25Sig" ,"vector<double>", &t_e7x7_25Sig);
1533  tree->Branch("t_e9x9_25Sig" ,"vector<double>", &t_e9x9_25Sig);
1534  tree->Branch("t_e11x11_25Sig" ,"vector<double>", &t_e11x11_25Sig);
1535  tree->Branch("t_e15x15_25Sig" ,"vector<double>", &t_e15x15_25Sig);
1536  tree->Branch("t_e7x7_30Sig" ,"vector<double>", &t_e7x7_30Sig);
1537  tree->Branch("t_e9x9_30Sig" ,"vector<double>", &t_e9x9_30Sig);
1538  tree->Branch("t_e11x11_30Sig" ,"vector<double>", &t_e11x11_30Sig);
1539  tree->Branch("t_e15x15_30Sig" ,"vector<double>", &t_e15x15_30Sig);
1540 
1541  if (doMC) {
1542  t_esim7x7 = new std::vector<double>();
1543  t_esim9x9 = new std::vector<double>();
1544  t_esim11x11 = new std::vector<double>();
1545  t_esim15x15 = new std::vector<double>();
1546 
1547  t_esim7x7Matched = new std::vector<double>();
1548  t_esim9x9Matched = new std::vector<double>();
1549  t_esim11x11Matched = new std::vector<double>();
1550  t_esim15x15Matched = new std::vector<double>();
1551 
1552  t_esim7x7Rest = new std::vector<double>();
1553  t_esim9x9Rest = new std::vector<double>();
1554  t_esim11x11Rest = new std::vector<double>();
1555  t_esim15x15Rest = new std::vector<double>();
1556 
1557  t_esim7x7Photon = new std::vector<double>();
1558  t_esim9x9Photon = new std::vector<double>();
1559  t_esim11x11Photon = new std::vector<double>();
1560  t_esim15x15Photon = new std::vector<double>();
1561 
1562  t_esim7x7NeutHad = new std::vector<double>();
1563  t_esim9x9NeutHad = new std::vector<double>();
1564  t_esim11x11NeutHad = new std::vector<double>();
1565  t_esim15x15NeutHad = new std::vector<double>();
1566 
1567  t_esim7x7CharHad = new std::vector<double>();
1568  t_esim9x9CharHad = new std::vector<double>();
1569  t_esim11x11CharHad = new std::vector<double>();
1570  t_esim15x15CharHad = new std::vector<double>();
1571 
1572  t_trkEcalEne = new std::vector<double>();
1573  t_simTrackP = new std::vector<double>();
1574  t_esimPdgId = new std::vector<double>();
1575 
1576  tree->Branch("t_esim7x7", "vector<double>", &t_esim7x7);
1577  tree->Branch("t_esim9x9", "vector<double>", &t_esim9x9);
1578  tree->Branch("t_esim11x11", "vector<double>", &t_esim11x11);
1579  tree->Branch("t_esim15x15", "vector<double>", &t_esim15x15);
1580 
1581  tree->Branch("t_esim7x7Matched", "vector<double>", &t_esim7x7Matched);
1582  tree->Branch("t_esim9x9Matched", "vector<double>", &t_esim9x9Matched);
1583  tree->Branch("t_esim11x11Matched", "vector<double>", &t_esim11x11Matched);
1584  tree->Branch("t_esim15x15Matched", "vector<double>", &t_esim15x15Matched);
1585 
1586  tree->Branch("t_esim7x7Rest", "vector<double>", &t_esim7x7Rest);
1587  tree->Branch("t_esim9x9Rest", "vector<double>", &t_esim9x9Rest);
1588  tree->Branch("t_esim11x11Rest", "vector<double>", &t_esim11x11Rest);
1589  tree->Branch("t_esim15x15Rest", "vector<double>", &t_esim15x15Rest);
1590 
1591  tree->Branch("t_esim7x7Photon", "vector<double>", &t_esim7x7Photon);
1592  tree->Branch("t_esim9x9Photon", "vector<double>", &t_esim9x9Photon);
1593  tree->Branch("t_esim11x11Photon", "vector<double>", &t_esim11x11Photon);
1594  tree->Branch("t_esim15x15Photon", "vector<double>", &t_esim15x15Photon);
1595 
1596 ;
1597  tree->Branch("t_esim7x7NeutHad", "vector<double>", &t_esim7x7NeutHad);
1598  tree->Branch("t_esim9x9NeutHad", "vector<double>", &t_esim9x9NeutHad);
1599  tree->Branch("t_esim11x11NeutHad", "vector<double>", &t_esim11x11NeutHad);
1600  tree->Branch("t_esim15x15NeutHad", "vector<double>", &t_esim15x15NeutHad);
1601 
1602  tree->Branch("t_esim7x7CharHad", "vector<double>", &t_esim7x7CharHad);
1603  tree->Branch("t_esim9x9CharHad", "vector<double>", &t_esim9x9CharHad);
1604  tree->Branch("t_esim11x11CharHad", "vector<double>", &t_esim11x11CharHad);
1605  tree->Branch("t_esim15x15CharHad", "vector<double>", &t_esim15x15CharHad);
1606 
1607  tree->Branch("t_trkEcalEne", "vector<double>", &t_trkEcalEne);
1608  tree->Branch("t_simTrackP", "vector<double>", &t_simTrackP);
1609  tree->Branch("t_esimPdgId", "vector<double>", &t_esimPdgId);
1610  }
1611 
1612  t_maxNearHcalP3x3 = new std::vector<double>();
1613  t_maxNearHcalP5x5 = new std::vector<double>();
1614  t_maxNearHcalP7x7 = new std::vector<double>();
1615  t_h3x3 = new std::vector<double>();
1616  t_h5x5 = new std::vector<double>();
1617  t_h7x7 = new std::vector<double>();
1618  t_h3x3Sig = new std::vector<double>();
1619  t_h5x5Sig = new std::vector<double>();
1620  t_h7x7Sig = new std::vector<double>();
1621  t_infoHcal = new std::vector<int>();
1622 
1623  if (doMC) {
1624  t_trkHcalEne = new std::vector<double>();
1625  t_hsim3x3 = new std::vector<double>();
1626  t_hsim5x5 = new std::vector<double>();
1627  t_hsim7x7 = new std::vector<double>();
1628  t_hsim3x3Matched = new std::vector<double>();
1629  t_hsim5x5Matched = new std::vector<double>();
1630  t_hsim7x7Matched = new std::vector<double>();
1631  t_hsim3x3Rest = new std::vector<double>();
1632  t_hsim5x5Rest = new std::vector<double>();
1633  t_hsim7x7Rest = new std::vector<double>();
1634  t_hsim3x3Photon = new std::vector<double>();
1635  t_hsim5x5Photon = new std::vector<double>();
1636  t_hsim7x7Photon = new std::vector<double>();
1637  t_hsim3x3NeutHad = new std::vector<double>();
1638  t_hsim5x5NeutHad = new std::vector<double>();
1639  t_hsim7x7NeutHad = new std::vector<double>();
1640  t_hsim3x3CharHad = new std::vector<double>();
1641  t_hsim5x5CharHad = new std::vector<double>();
1642  t_hsim7x7CharHad = new std::vector<double>();
1643  }
1644 
1645  tree->Branch("t_maxNearHcalP3x3", "vector<double>", &t_maxNearHcalP3x3);
1646  tree->Branch("t_maxNearHcalP5x5", "vector<double>", &t_maxNearHcalP5x5);
1647  tree->Branch("t_maxNearHcalP7x7", "vector<double>", &t_maxNearHcalP7x7);
1648  tree->Branch("t_h3x3", "vector<double>", &t_h3x3);
1649  tree->Branch("t_h5x5", "vector<double>", &t_h5x5);
1650  tree->Branch("t_h7x7", "vector<double>", &t_h7x7);
1651  tree->Branch("t_h3x3Sig", "vector<double>", &t_h3x3Sig);
1652  tree->Branch("t_h5x5Sig", "vector<double>", &t_h5x5Sig);
1653  tree->Branch("t_h7x7Sig", "vector<double>", &t_h7x7Sig);
1654  tree->Branch("t_infoHcal", "vector<int>", &t_infoHcal);
1655 
1656  if (doMC) {
1657  tree->Branch("t_trkHcalEne", "vector<double>", &t_trkHcalEne);
1658  tree->Branch("t_hsim3x3", "vector<double>", &t_hsim3x3);
1659  tree->Branch("t_hsim5x5", "vector<double>", &t_hsim5x5);
1660  tree->Branch("t_hsim7x7", "vector<double>", &t_hsim7x7);
1661  tree->Branch("t_hsim3x3Matched", "vector<double>", &t_hsim3x3Matched);
1662  tree->Branch("t_hsim5x5Matched", "vector<double>", &t_hsim5x5Matched);
1663  tree->Branch("t_hsim7x7Matched", "vector<double>", &t_hsim7x7Matched);
1664  tree->Branch("t_hsim3x3Rest", "vector<double>", &t_hsim3x3Rest);
1665  tree->Branch("t_hsim5x5Rest", "vector<double>", &t_hsim5x5Rest);
1666  tree->Branch("t_hsim7x7Rest", "vector<double>", &t_hsim7x7Rest);
1667  tree->Branch("t_hsim3x3Photon", "vector<double>", &t_hsim3x3Photon);
1668  tree->Branch("t_hsim5x5Photon", "vector<double>", &t_hsim5x5Photon);
1669  tree->Branch("t_hsim7x7Photon", "vector<double>", &t_hsim7x7Photon);
1670  tree->Branch("t_hsim3x3NeutHad", "vector<double>", &t_hsim3x3NeutHad);
1671  tree->Branch("t_hsim5x5NeutHad", "vector<double>", &t_hsim5x5NeutHad);
1672  tree->Branch("t_hsim7x7NeutHad", "vector<double>", &t_hsim7x7NeutHad);
1673  tree->Branch("t_hsim3x3CharHad", "vector<double>", &t_hsim3x3CharHad);
1674  tree->Branch("t_hsim5x5CharHad", "vector<double>", &t_hsim5x5CharHad);
1675  tree->Branch("t_hsim7x7CharHad", "vector<double>", &t_hsim7x7CharHad);
1676  }
1677  tree->Branch("t_nTracks", &t_nTracks, "t_nTracks/I");
1678 
1679 }
1680 
1681 
1682 double IsolatedTracksNxN::DeltaPhi(double v1, double v2) {
1683  // Computes the correctly normalized phi difference
1684  // v1, v2 = phi of object 1 and 2
1685 
1686  double pi = 3.141592654;
1687  double twopi = 6.283185307;
1688 
1689  double diff = std::abs(v2 - v1);
1690  double corr = twopi - diff;
1691  if (diff < pi){ return diff;} else { return corr;}
1692 }
1693 
1694 double IsolatedTracksNxN::DeltaR(double eta1, double phi1, double eta2, double phi2) {
1695  double deta = eta1 - eta2;
1696  double dphi = DeltaPhi(phi1, phi2);
1697  return std::sqrt(deta*deta + dphi*dphi);
1698 }
1699 
1700 void IsolatedTracksNxN::printTrack(const reco::Track* pTrack) {
1701 
1702  std::string theTrackQuality = "highPurity";
1703  reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality);
1704 
1705  std::cout << " Reference Point " << pTrack->referencePoint() <<"\n"
1706  << " TrackMmentum " << pTrack->momentum()
1707  << " (pt,eta,phi)(" << pTrack->pt()<<","<<pTrack->eta()<<","<<pTrack->phi()<<")"
1708  << " p " << pTrack->p() << "\n"
1709  << " Normalized chi2 " << pTrack->normalizedChi2() <<" charge " << pTrack->charge()
1710  << " qoverp() " << pTrack->qoverp() <<"+-" << pTrack->qoverpError()
1711  << " d0 " << pTrack->d0() << "\n"
1712  << " NValidHits " << pTrack->numberOfValidHits() << " NLostHits " << pTrack->numberOfLostHits()
1713  << " TrackQuality " << pTrack->qualityName(trackQuality_) << " " << pTrack->quality(trackQuality_)
1714  << std::endl;
1715 
1716  if( printTrkHitPattern_ ) {
1717  const reco::HitPattern& p = pTrack->hitPattern();
1718 
1719  std::cout<<"default " << std::endl;
1720  for (int i=0; i<p.numberOfAllHits(reco::HitPattern::TRACK_HITS); i++) {
1721  p.printHitPattern(reco::HitPattern::TRACK_HITS, i, std::cout);
1722  }
1723  std::cout<<"trackerMissingInnerHits() " << std::endl;
1724  for (int i=0; i<p.numberOfAllHits(reco::HitPattern::MISSING_INNER_HITS); i++) {
1725  p.printHitPattern(reco::HitPattern::MISSING_INNER_HITS, i, std::cout);
1726  }
1727  std::cout<<"trackerMissingOuterHits() " << std::endl;
1728  for (int i=0; i<p.numberOfAllHits(reco::HitPattern::MISSING_OUTER_HITS); i++) {
1729  p.printHitPattern(reco::HitPattern::MISSING_OUTER_HITS, i, std::cout);
1730  }
1731 
1732 
1733  std::cout << "\n \t trackerLayersWithMeasurement() " << p.trackerLayersWithMeasurement()
1734  << "\n \t pixelLayersWithMeasurement() " << p.pixelLayersWithMeasurement()
1735  << "\n \t stripLayersWithMeasurement() " << p.stripLayersWithMeasurement()
1736  << "\n \t pixelBarrelLayersWithMeasurement() " << p.pixelBarrelLayersWithMeasurement()
1737  << "\n \t pixelEndcapLayersWithMeasurement() " << p.pixelEndcapLayersWithMeasurement()
1738  << "\n \t stripTIBLayersWithMeasurement() " << p.stripTIBLayersWithMeasurement()
1739  << "\n \t stripTIDLayersWithMeasurement() " << p.stripTIDLayersWithMeasurement()
1740  << "\n \t stripTOBLayersWithMeasurement() " << p.stripTOBLayersWithMeasurement()
1741  << "\n \t stripTECLayersWithMeasurement() " << p.stripTECLayersWithMeasurement()
1742  << std::endl;
1743 
1744  }
1745 }
1746 
1747 //define this as a plug-in
1748 DEFINE_FWK_MODULE(IsolatedTracksNxN);
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:39
IsolatedTracksNxN::t_hsim7x7CharHad
std::vector< double > * t_hsim7x7CharHad
Definition: IsolatedTracksNxN.h:291
reco::TrackBase::qoverp
double qoverp() const
q / p
Definition: TrackBase.h:573
IsolatedTracksNxN::DeltaPhi
double DeltaPhi(double v1, double v2)
Definition: IsolatedTracksNxN.cc:1682
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:615
IsolatedTracksNxN::t_hsim3x3Matched
std::vector< double > * t_hsim3x3Matched
Definition: IsolatedTracksNxN.h:287
reco::HitPattern::stripTOBLayersWithMeasurement
int stripTOBLayersWithMeasurement() const
Definition: HitPattern.cc:604
IsolatedTracksNxN::t_esim7x7CharHad
std::vector< double > * t_esim7x7CharHad
Definition: IsolatedTracksNxN.h:278
IsolatedTracksNxN::tok_L1extCenJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extCenJet_
Definition: IsolatedTracksNxN.h:160
L1GtTriggerMenu::gtTriggerMenuName
const std::string & gtTriggerMenuName() const
Definition: L1GtTriggerMenu.h:99
IsolatedTracksNxN::t_e15x15
std::vector< double > * t_e15x15
Definition: IsolatedTracksNxN.h:262
edm::RefToBase::get
value_type const * get() const
Definition: RefToBase.h:234
IsolatedTracksNxN::t_nTrksJetCalo
std::vector< double > * t_nTrksJetCalo
Definition: IsolatedTracksNxN.h:237
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:41
IsolatedTracksNxN::t_maxNearHcalP7x7
std::vector< double > * t_maxNearHcalP7x7
Definition: IsolatedTracksNxN.h:280
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
reco::TrackBase::referencePoint
const Point & referencePoint() const
Reference point on the track.
Definition: TrackBase.h:681
IsolatedTracksNxN::t_trackPt
std::vector< double > * t_trackPt
Definition: IsolatedTracksNxN.h:243
IsolatedTracksNxN::t_trackPAll
std::vector< double > * t_trackPAll
Definition: IsolatedTracksNxN.h:239
IsolatedTracksNxN::t_e7x7
std::vector< double > * t_e7x7
Definition: IsolatedTracksNxN.h:262
IsolatedTracksNxN::t_trackHcalPhi
std::vector< double > * t_trackHcalPhi
Definition: IsolatedTracksNxN.h:244
IsolatedTracksNxN::t_trackDxyPVAll
std::vector< double > * t_trackDxyPVAll
Definition: IsolatedTracksNxN.h:241
nTracks
const unsigned int nTracks(const reco::Vertex &sv)
Definition: TemplatedVertexArbitrator.h:48
IsolatedTracksNxN::~IsolatedTracksNxN
~IsolatedTracksNxN() override
Definition: IsolatedTracksNxN.cc:112
reco::TrackBase::qualityName
static std::string qualityName(TrackQuality)
Definition: TrackBase.h:532
IsolatedTracksNxN::t_trackHitOutMeasTEC
std::vector< int > * t_trackHitOutMeasTEC
Definition: IsolatedTracksNxN.h:255
IsolatedTracksNxN::t_trackPdgIdAll
std::vector< double > * t_trackPdgIdAll
Definition: IsolatedTracksNxN.h:239
mps_fire.i
i
Definition: mps_fire.py:269
spr::propagateCALO
std::vector< spr::propagatedTrackID > propagateCALO(edm::Handle< reco::TrackCollection > &trkCollection, const CaloGeometry *geo, const MagneticField *bField, std::string &theTrackQuality, bool debug=false)
Definition: CaloPropagateTrack.cc:117
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::t_trackOutPosOutHitDr
std::vector< double > * t_trackOutPosOutHitDr
Definition: IsolatedTracksNxN.h:256
IsolatedTracksNxN::t_maxNearHcalP5x5
std::vector< double > * t_maxNearHcalP5x5
Definition: IsolatedTracksNxN.h:280
IsolatedTracksNxN::t_L1NonIsoEMPt
std::vector< double > * t_L1NonIsoEMPt
Definition: IsolatedTracksNxN.h:233
IsolatedTracksNxN::t_trackHitOutMeasTIB
std::vector< int > * t_trackHitOutMeasTIB
Definition: IsolatedTracksNxN.h:255
IsolatedTracksNxN::t_esim7x7
std::vector< double > * t_esim7x7
Definition: IsolatedTracksNxN.h:273
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
IsolatedTracksNxN::t_trackHitInMissTOB
std::vector< int > * t_trackHitInMissTOB
Definition: IsolatedTracksNxN.h:252
IsolatedTracksNxN::t_trackDz
std::vector< double > * t_trackDz
Definition: IsolatedTracksNxN.h:245
reco::TrackBase::d0
double d0() const
dxy parameter in perigee convention (d0 = -dxy)
Definition: TrackBase.h:597
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::t_trackEcalEta
std::vector< double > * t_trackEcalEta
Definition: IsolatedTracksNxN.h:244
reco::HitPattern::stripTIBLayersWithMeasurement
int stripTIBLayersWithMeasurement() const
Definition: HitPattern.cc:582
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:561
IsolatedTracksNxN::t_trackEtaAll
std::vector< double > * t_trackEtaAll
Definition: IsolatedTracksNxN.h:239
IsolatedTracksNxN::t_esim11x11Matched
std::vector< double > * t_esim11x11Matched
Definition: IsolatedTracksNxN.h:274
IsolatedTracksNxN::t_L1IsoEMEta
std::vector< double > * t_L1IsoEMEta
Definition: IsolatedTracksNxN.h:232
IsolatedTracksNxN::t_PVTracksSumPt
std::vector< double > * t_PVTracksSumPt
Definition: IsolatedTracksNxN.h:223
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:508
IsolatedTracksNxN::t_e9x9_10Sig
std::vector< double > * t_e9x9_10Sig
Definition: IsolatedTracksNxN.h:263
IsolatedTracksNxN::h_maxNearP25x25
TH1F * h_maxNearP25x25[NPBins][NEtaBins]
Definition: IsolatedTracksNxN.h:193
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:151
IsolatedTracksNxN::t_hsim3x3CharHad
std::vector< double > * t_hsim3x3CharHad
Definition: IsolatedTracksNxN.h:291
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
reco::HitPattern::stripTOBLayersWithoutMeasurement
int stripTOBLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:672
IsolatedTracksNxN::t_PVndof
std::vector< int > * t_PVndof
Definition: IsolatedTracksNxN.h:225
AlgorithmMap
std::map< std::string, L1GtAlgorithm > AlgorithmMap
map containing the algorithms
Definition: L1GtTriggerMenuFwd.h:31
IsolatedTracksNxN::t_hsim7x7Rest
std::vector< double > * t_hsim7x7Rest
Definition: IsolatedTracksNxN.h:288
IsolatedTracksNxN::t_L1TauJetPhi
std::vector< double > * t_L1TauJetPhi
Definition: IsolatedTracksNxN.h:230
IsolatedTracksNxN::t_PVTracksSumPtWt
std::vector< double > * t_PVTracksSumPtWt
Definition: IsolatedTracksNxN.h:224
L1GtUtils::availableL1Configuration
const bool availableL1Configuration(int &errorCode, int &l1ConfCode) const
Definition: L1GtUtils.cc:1948
IsolatedTracksNxN::t_infoHcal
std::vector< int > * t_infoHcal
Definition: IsolatedTracksNxN.h:283
IsolatedTracksNxN::t_e11x11_15Sig
std::vector< double > * t_e11x11_15Sig
Definition: IsolatedTracksNxN.h:264
IsolatedTracksNxN::NEtaBins
static const size_t NEtaBins
Definition: IsolatedTracksNxN.h:191
IsolatedTracksNxN::t_L1TauJetPt
std::vector< double > * t_L1TauJetPt
Definition: IsolatedTracksNxN.h:230
IsolatedTracksNxN::t_PVNTracksWt
std::vector< int > * t_PVNTracksWt
Definition: IsolatedTracksNxN.h:225
IsolatedTracksNxN::t_trackEcalPhi
std::vector< double > * t_trackEcalPhi
Definition: IsolatedTracksNxN.h:244
IsolatedTracksNxN::t_esim11x11NeutHad
std::vector< double > * t_esim11x11NeutHad
Definition: IsolatedTracksNxN.h:277
IsolatedTracksNxN::t_L1TauJetEta
std::vector< double > * t_L1TauJetEta
Definition: IsolatedTracksNxN.h:230
IsolatedTracksNxN::t_e15x15_15Sig
std::vector< double > * t_e15x15_15Sig
Definition: IsolatedTracksNxN.h:264
IsolatedTracksNxN::t_L1MuonEta
std::vector< double > * t_L1MuonEta
Definition: IsolatedTracksNxN.h:231
IsolatedTracksNxN::t_L1MuonPt
std::vector< double > * t_L1MuonPt
Definition: IsolatedTracksNxN.h:231
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, bool useRaw=false, bool debug=false)
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:66
IsolatedTracksNxN::t_h3x3
std::vector< double > * t_h3x3
Definition: IsolatedTracksNxN.h:281
edm::EventBase::luminosityBlock
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:63
IsolatedTracksNxN::t_L1MuonPhi
std::vector< double > * t_L1MuonPhi
Definition: IsolatedTracksNxN.h:231
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:645
reco::TrackBase::numberOfLostHits
unsigned short numberOfLostHits() const
number of cases where track crossed a layer without getting a hit.
Definition: TrackBase.h:826
L1GtUtils::decision
const bool decision(const edm::Event &iEvent, const std::string &nameAlgoTechTrig, int &errorCode) const
Definition: L1GtUtils.cc:1208
weight
Definition: weight.py:1
IsolatedTracksNxN::t_trackDxyAll
std::vector< double > * t_trackDxyAll
Definition: IsolatedTracksNxN.h:241
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
IsolatedTracksNxN::t_maxNearHcalP3x3
std::vector< double > * t_maxNearHcalP3x3
Definition: IsolatedTracksNxN.h:280
IsolatedTracksNxN::IsolatedTracksNxN
IsolatedTracksNxN(const edm::ParameterSet &)
Definition: IsolatedTracksNxN.cc:38
IsolatedTracksNxN::tok_simVtx_
edm::EDGetTokenT< edm::SimVertexContainer > tok_simVtx_
Definition: IsolatedTracksNxN.h:177
IsolatedTracksNxN::t_hsim3x3Photon
std::vector< double > * t_hsim3x3Photon
Definition: IsolatedTracksNxN.h:289
IsolatedTracksNxN::tok_bs_
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
Definition: IsolatedTracksNxN.h:172
reco::HitPattern::pixelLayersWithMeasurement
int pixelLayersWithMeasurement() const
Definition: HitPattern.cc:499
reco::TrackBase::momentum
const Vector & momentum() const
track momentum vector
Definition: TrackBase.h:675
IsolatedTracksNxN::t_PVTracksSumPtHP
std::vector< double > * t_PVTracksSumPtHP
Definition: IsolatedTracksNxN.h:224
reco::Track::outerPosition
const math::XYZPoint & outerPosition() const
position of the outermost hit
Definition: Track.h:65
IsolatedTracksNxN::t_esim11x11Rest
std::vector< double > * t_esim11x11Rest
Definition: IsolatedTracksNxN.h:275
reco::HitPattern::trackerLayersWithMeasurement
int trackerLayersWithMeasurement() const
Definition: HitPattern.cc:518
reco::HitPattern::pixelEndcapLayersWithMeasurement
int pixelEndcapLayersWithMeasurement() const
Definition: HitPattern.cc:570
IsolatedTracksNxN::t_trackNOuterHits
std::vector< int > * t_trackNOuterHits
Definition: IsolatedTracksNxN.h:250
IsolatedTracksNxN::t_h7x7Sig
std::vector< double > * t_h7x7Sig
Definition: IsolatedTracksNxN.h:282
IsolatedTracksNxN::t_e15x15_30Sig
std::vector< double > * t_e15x15_30Sig
Definition: IsolatedTracksNxN.h:267
IsolatedTracksNxN::t_trackPhiAll
std::vector< double > * t_trackPhiAll
Definition: IsolatedTracksNxN.h:239
IsolatedTracksNxN::t_e7x7_25Sig
std::vector< double > * t_e7x7_25Sig
Definition: IsolatedTracksNxN.h:266
IsolatedTracksNxN::t_trkHcalEne
std::vector< double > * t_trkHcalEne
Definition: IsolatedTracksNxN.h:285
IsolatedTracksNxN::t_e7x7_20Sig
std::vector< double > * t_e7x7_20Sig
Definition: IsolatedTracksNxN.h:265
IsolatedTracksNxN::t_e11x11_10Sig
std::vector< double > * t_e11x11_10Sig
Definition: IsolatedTracksNxN.h:263
pi
const Double_t pi
Definition: trackSplitPlot.h:36
CItAlgo
AlgorithmMap::const_iterator CItAlgo
iterators through map containing the algorithms
Definition: L1GtTriggerMenuFwd.h:38
IsolatedTracksNxN::t_trackPtAll
std::vector< double > * t_trackPtAll
Definition: IsolatedTracksNxN.h:240
IsolatedTracksNxN::t_h3x3Sig
std::vector< double > * t_h3x3Sig
Definition: IsolatedTracksNxN.h:282
photonIsolationHIProducer_cfi.hbhe
hbhe
Definition: photonIsolationHIProducer_cfi.py:8
IsolatedTracksNxN::fs
edm::Service< TFileService > fs
Definition: IsolatedTracksNxN.h:293
IsolatedTracksNxN::t_esimPdgId
std::vector< double > * t_esimPdgId
Definition: IsolatedTracksNxN.h:269
spr::caloSimInfo::eNeutralHad
double eNeutralHad
Definition: CaloSimInfo.h:47
reco::HitPattern::printHitPattern
void printHitPattern(HitCategory category, int position, std::ostream &stream) const
Definition: HitPattern.cc:832
spr::caloSimInfo::eChargedHad
double eChargedHad
Definition: CaloSimInfo.h:48
IsolatedTracksNxN::t_e11x11_20Sig
std::vector< double > * t_e11x11_20Sig
Definition: IsolatedTracksNxN.h:265
IsolatedTracksNxN::t_jetEta
std::vector< double > * t_jetEta
Definition: IsolatedTracksNxN.h:236
IsolatedTracksNxN::l1AlgoMap
std::map< std::pair< unsigned int, std::string >, int > l1AlgoMap
Definition: IsolatedTracksNxN.h:148
IsolatedTracksNxN::t_maxNearP31x31
std::vector< double > * t_maxNearP31x31
Definition: IsolatedTracksNxN.h:258
IsolatedTracksNxN::t_trackHitInMissTEC
std::vector< int > * t_trackHitInMissTEC
Definition: IsolatedTracksNxN.h:252
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:16
iEvent
int iEvent
Definition: GenABIO.cc:230
IsolatedTracksNxN::t_trackEta
std::vector< double > * t_trackEta
Definition: IsolatedTracksNxN.h:243
IsolatedTracksNxN::t_esim9x9Photon
std::vector< double > * t_esim9x9Photon
Definition: IsolatedTracksNxN.h:276
IsolatedTracksNxN::t_trackHitInMeasTID
std::vector< int > * t_trackHitInMeasTID
Definition: IsolatedTracksNxN.h:254
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:651
IsolatedTracksNxN::t_hsim5x5Rest
std::vector< double > * t_hsim5x5Rest
Definition: IsolatedTracksNxN.h:288
IsolatedTracksNxN::t_e11x11
std::vector< double > * t_e11x11
Definition: IsolatedTracksNxN.h:262
IsolatedTracksNxN::t_esim15x15
std::vector< double > * t_esim15x15
Definition: IsolatedTracksNxN.h:273
IsolatedTracksNxN::tok_simTk_
edm::EDGetTokenT< edm::SimTrackContainer > tok_simTk_
Definition: IsolatedTracksNxN.h:176
IsolatedTracksNxN::t_L1CenJetPt
std::vector< double > * t_L1CenJetPt
Definition: IsolatedTracksNxN.h:228
IsolatedTracksNxN::t_trackChiSqAll
std::vector< double > * t_trackChiSqAll
Definition: IsolatedTracksNxN.h:241
IsolatedTracksNxN::t_e15x15_25Sig
std::vector< double > * t_e15x15_25Sig
Definition: IsolatedTracksNxN.h:266
IsolatedTracksNxN::t_trackHitOutMissTIB
std::vector< int > * t_trackHitOutMissTIB
Definition: IsolatedTracksNxN.h:253
IsolatedTracksNxN::tok_caloHH_
edm::EDGetTokenT< edm::PCaloHitContainer > tok_caloHH_
Definition: IsolatedTracksNxN.h:180
IsolatedTracksNxN::t_esim15x15Rest
std::vector< double > * t_esim15x15Rest
Definition: IsolatedTracksNxN.h:275
IsolatedTracksNxN::t_esim11x11
std::vector< double > * t_esim11x11
Definition: IsolatedTracksNxN.h:273
IsolatedTracksNxN::t_hsim7x7
std::vector< double > * t_hsim7x7
Definition: IsolatedTracksNxN.h:286
IsolatedTracksNxN::t_PVNTracksHPWt
std::vector< int > * t_PVNTracksHPWt
Definition: IsolatedTracksNxN.h:226
reco::HitPattern::stripTIDLayersWithMeasurement
int stripTIDLayersWithMeasurement() const
Definition: HitPattern.cc:593
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
IsolatedTracksNxN::t_trackHitOutMissTID
std::vector< int > * t_trackHitOutMissTID
Definition: IsolatedTracksNxN.h:253
IsolatedTracksNxN::t_trackHitOutMissTOBTEC
std::vector< int > * t_trackHitOutMissTOBTEC
Definition: IsolatedTracksNxN.h:253
fwrapper::jets
vector< PseudoJet > jets
Definition: fastjetfortran_madfks.cc:42
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:621
IsolatedTracksNxN::t_trackHitsTEC
std::vector< int > * t_trackHitsTEC
Definition: IsolatedTracksNxN.h:251
IsolatedTracksNxN::t_e9x9_30Sig
std::vector< double > * t_e9x9_30Sig
Definition: IsolatedTracksNxN.h:267
IsolatedTracksNxN::bField
const MagneticField * bField
Definition: IsolatedTracksNxN.h:186
IsolatedTracksNxN::tok_caloEB_
edm::EDGetTokenT< edm::PCaloHitContainer > tok_caloEB_
Definition: IsolatedTracksNxN.h:178
IsolatedTracksNxN::t_L1CenJetPhi
std::vector< double > * t_L1CenJetPhi
Definition: IsolatedTracksNxN.h:228
reco::HitPattern::numberOfAllHits
int numberOfAllHits(HitCategory category) const
Definition: HitPattern.h:807
IsolatedTracksNxN::t_L1FwdJetEta
std::vector< double > * t_L1FwdJetEta
Definition: IsolatedTracksNxN.h:229
IsolatedTracksNxN::t_h7x7
std::vector< double > * t_h7x7
Definition: IsolatedTracksNxN.h:281
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
IsolatedTracksNxN::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: IsolatedTracksNxN.cc:115
IsolatedTracksNxN::t_ecalSpike11x11
std::vector< int > * t_ecalSpike11x11
Definition: IsolatedTracksNxN.h:261
IsolatedTracksNxN::t_esim15x15NeutHad
std::vector< double > * t_esim15x15NeutHad
Definition: IsolatedTracksNxN.h:277
IsolatedTracksNxN::t_e9x9_20Sig
std::vector< double > * t_e9x9_20Sig
Definition: IsolatedTracksNxN.h:265
IsolatedTracksNxN::t_trackHitOutMeasTOB
std::vector< int > * t_trackHitOutMeasTOB
Definition: IsolatedTracksNxN.h:255
IsolatedTracksNxN::t_trackHitInMeasTIB
std::vector< int > * t_trackHitInMeasTIB
Definition: IsolatedTracksNxN.h:254
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
reco::TrackBase::numberOfValidHits
unsigned short numberOfValidHits() const
number of valid hits found
Definition: TrackBase.h:820
IsolatedTracksNxN::tok_genTrack_
edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
Definition: IsolatedTracksNxN.h:170
TFileDirectory::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileDirectory.h:59
IsolatedTracksNxN::t_trackDzPV
std::vector< double > * t_trackDzPV
Definition: IsolatedTracksNxN.h:246
IsolatedTracksNxN::t_PVNTracksHP
std::vector< int > * t_PVNTracksHP
Definition: IsolatedTracksNxN.h:226
IsolatedTracksNxN::t_L1FwdJetPhi
std::vector< double > * t_L1FwdJetPhi
Definition: IsolatedTracksNxN.h:229
IsolatedTracksNxN::tok_L1extFwdJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extFwdJet_
Definition: IsolatedTracksNxN.h:161
IsolatedTracksNxN::tok_L1extIsoEm_
edm::EDGetTokenT< l1extra::L1EmParticleCollection > tok_L1extIsoEm_
Definition: IsolatedTracksNxN.h:164
IsolatedTracksNxN::t_nTrksJetVtx
std::vector< double > * t_nTrksJetVtx
Definition: IsolatedTracksNxN.h:237
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
IsolatedTracksNxN::t_esim7x7NeutHad
std::vector< double > * t_esim7x7NeutHad
Definition: IsolatedTracksNxN.h:277
IsolatedTracksNxN::t_hsim3x3NeutHad
std::vector< double > * t_hsim3x3NeutHad
Definition: IsolatedTracksNxN.h:290
IsolatedTracksNxN::t_L1NonIsoEMPhi
std::vector< double > * t_L1NonIsoEMPhi
Definition: IsolatedTracksNxN.h:233
IsolatedTracksNxN::t_trackDzBS
std::vector< double > * t_trackDzBS
Definition: IsolatedTracksNxN.h:245
IsolatedTracksNxN::t_trackHitOutMissTOB
std::vector< int > * t_trackHitOutMissTOB
Definition: IsolatedTracksNxN.h:253
IsolatedTracksNxN::t_L1CenJetEta
std::vector< double > * t_L1CenJetEta
Definition: IsolatedTracksNxN.h:228
IsolatedTracksNxN::t_trackDzPVAll
std::vector< double > * t_trackDzPVAll
Definition: IsolatedTracksNxN.h:241
corr
JetCorrectorParameters corr
Definition: classes.h:5
IsolatedTracksNxN::tok_L1extMu_
edm::EDGetTokenT< l1extra::L1MuonParticleCollection > tok_L1extMu_
Definition: IsolatedTracksNxN.h:163
cuy.ii
ii
Definition: cuy.py:588
IsolatedTracksNxN::t_esim7x7Matched
std::vector< double > * t_esim7x7Matched
Definition: IsolatedTracksNxN.h:274
IsolatedTracksNxN::t_esim9x9
std::vector< double > * t_esim9x9
Definition: IsolatedTracksNxN.h:273
IsolatedTracksNxN::t_hsim5x5
std::vector< double > * t_hsim5x5
Definition: IsolatedTracksNxN.h:286
IsolatedTracksNxN::t_h5x5
std::vector< double > * t_h5x5
Definition: IsolatedTracksNxN.h:281
reco::TrackBase::qoverpError
double qoverpError() const
error on signed transverse curvature
Definition: TrackBase.h:757
IsolatedTracksNxN::t_trackDxyPV
std::vector< double > * t_trackDxyPV
Definition: IsolatedTracksNxN.h:246
IsolatedTracksNxN::t_hsim5x5NeutHad
std::vector< double > * t_hsim5x5NeutHad
Definition: IsolatedTracksNxN.h:290
IsolatedTracksNxN::t_trackChiSq
std::vector< double > * t_trackChiSq
Definition: IsolatedTracksNxN.h:247
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:609
hiDetachedQuadStep_cff.pt1
pt1
Definition: hiDetachedQuadStep_cff.py:70
IsolatedTracksNxN::t_hsim5x5Photon
std::vector< double > * t_hsim5x5Photon
Definition: IsolatedTracksNxN.h:289
IsolatedTracksNxN::t_PVy
std::vector< double > * t_PVy
Definition: IsolatedTracksNxN.h:223
IsolatedTracksNxN::t_L1PreScale
std::vector< int > * t_L1PreScale
Definition: IsolatedTracksNxN.h:220
IsolatedTracksNxN::t_hsim3x3Rest
std::vector< double > * t_hsim3x3Rest
Definition: IsolatedTracksNxN.h:288
reco::HitPattern::stripTIDLayersWithoutMeasurement
int stripTIDLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:661
DetId
Definition: DetId.h:18
IsolatedTracksNxN::t_PVisValid
std::vector< int > * t_PVisValid
Definition: IsolatedTracksNxN.h:225
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:125
IsolatedTracksNxN::t_L1METPhi
std::vector< double > * t_L1METPhi
Definition: IsolatedTracksNxN.h:234
IsolatedTracksNxN::t_esim9x9Rest
std::vector< double > * t_esim9x9Rest
Definition: IsolatedTracksNxN.h:275
IsolatedTracksNxN::t_trackHitInMissTIB
std::vector< int > * t_trackHitInMissTIB
Definition: IsolatedTracksNxN.h:252
TFileService::mkdir
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
Definition: TFileService.h:69
IsolatedTracksNxN::t_e9x9_15Sig
std::vector< double > * t_e9x9_15Sig
Definition: IsolatedTracksNxN.h:264
IsolatedTracksNxN::t_e7x7_15Sig
std::vector< double > * t_e7x7_15Sig
Definition: IsolatedTracksNxN.h:264
IsolatedTracksNxN::t_trackPhi
std::vector< double > * t_trackPhi
Definition: IsolatedTracksNxN.h:243
edm::Handle::product
T const * product() const
Definition: Handle.h:81
IsolatedTracksNxN::t_e7x7_10Sig
std::vector< double > * t_e7x7_10Sig
Definition: IsolatedTracksNxN.h:263
reco::TrackBase::hitPattern
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:446
IsolatedTracksNxN::endJob
void endJob() override
Definition: IsolatedTracksNxN.cc:968
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
IsolatedTracksNxN::t_trkEcalEne
std::vector< double > * t_trkEcalEne
Definition: IsolatedTracksNxN.h:271
IsolatedTracksNxN::tok_L1extNoIsoEm_
edm::EDGetTokenT< l1extra::L1EmParticleCollection > tok_L1extNoIsoEm_
Definition: IsolatedTracksNxN.h:165
IsolatedTracksNxN::t_trackDzAll
std::vector< double > * t_trackDzAll
Definition: IsolatedTracksNxN.h:241
IsolatedTracksNxN::tok_recVtx_
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
Definition: IsolatedTracksNxN.h:171
IsolatedTracksNxN::h_maxNearP31x31
TH1F * h_maxNearP31x31[NPBins][NEtaBins]
Definition: IsolatedTracksNxN.h:193
IsolatedTracksNxN::t_esim11x11CharHad
std::vector< double > * t_esim11x11CharHad
Definition: IsolatedTracksNxN.h:278
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
reco::HitPattern::stripLayersWithMeasurement
int stripLayersWithMeasurement() const
Definition: HitPattern.h:1035
reco::HitPattern::pixelBarrelLayersWithMeasurement
int pixelBarrelLayersWithMeasurement() const
Definition: HitPattern.cc:558
IsolatedTracksNxN::tok_caloEE_
edm::EDGetTokenT< edm::PCaloHitContainer > tok_caloEE_
Definition: IsolatedTracksNxN.h:179
IsolatedTracksNxN::tok_EB_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: IsolatedTracksNxN.h:174
IsolatedTracksNxN::h_maxNearP15x15
TH1F * h_maxNearP15x15[NPBins][NEtaBins]
Definition: IsolatedTracksNxN.h:193
IsolatedTracksNxN::t_e9x9_25Sig
std::vector< double > * t_e9x9_25Sig
Definition: IsolatedTracksNxN.h:266
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_esim7x7Photon
std::vector< double > * t_esim7x7Photon
Definition: IsolatedTracksNxN.h:276
useL1GtTriggerMenuLite
static const bool useL1GtTriggerMenuLite(true)
IsolatedTracksNxN::t_esim9x9NeutHad
std::vector< double > * t_esim9x9NeutHad
Definition: IsolatedTracksNxN.h:277
IsolatedTracksNxN::t_L1FwdJetPt
std::vector< double > * t_L1FwdJetPt
Definition: IsolatedTracksNxN.h:229
IsolatedTracksNxN::t_simTrackP
std::vector< double > * t_simTrackP
Definition: IsolatedTracksNxN.h:269
IsolatedTracksNxN::t_hsim7x7Matched
std::vector< double > * t_hsim7x7Matched
Definition: IsolatedTracksNxN.h:287
IsolatedTracksNxN::t_trackPVIdx
std::vector< int > * t_trackPVIdx
Definition: IsolatedTracksNxN.h:248
reco::TrackBase::quality
bool quality(const TrackQuality) const
Track quality.
Definition: TrackBase.h:510
IsolatedTracksNxN::t_trackHitsTOB
std::vector< int > * t_trackHitsTOB
Definition: IsolatedTracksNxN.h:251
IsolatedTracksNxN::t_esim9x9Matched
std::vector< double > * t_esim9x9Matched
Definition: IsolatedTracksNxN.h:274
IsolatedTracksNxN::t_hsim7x7NeutHad
std::vector< double > * t_hsim7x7NeutHad
Definition: IsolatedTracksNxN.h:290
IsolatedTracksNxN::t_L1IsoEMPhi
std::vector< double > * t_L1IsoEMPhi
Definition: IsolatedTracksNxN.h:232
IsolatedTracksNxN::t_PVNTracks
std::vector< int > * t_PVNTracks
Definition: IsolatedTracksNxN.h:225
IsolatedTracksNxN::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: IsolatedTracksNxN.h:168
IsolatedTracksNxN::t_maxNearP21x21
std::vector< double > * t_maxNearP21x21
Definition: IsolatedTracksNxN.h:259
IsolatedTracksNxN::trackerHitAssociatorConfig_
TrackerHitAssociator::Config trackerHitAssociatorConfig_
Definition: IsolatedTracksNxN.h:157
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:60
IsolatedTracksNxN::t_L1NonIsoEMEta
std::vector< double > * t_L1NonIsoEMEta
Definition: IsolatedTracksNxN.h:233
p1
double p1[4]
Definition: TauolaWrapper.h:89
IsolatedTracksNxN::t_PVz
std::vector< double > * t_PVz
Definition: IsolatedTracksNxN.h:223
spr::propagateTrackerEnd
std::pair< math::XYZPoint, double > propagateTrackerEnd(const reco::Track *, const MagneticField *, bool debug=false)
Definition: CaloPropagateTrack.cc:664
IsolatedTracksNxN::printTrack
void printTrack(const reco::Track *pTrack)
Definition: IsolatedTracksNxN.cc:1700
IsolatedTracksNxN::t_e15x15_20Sig
std::vector< double > * t_e15x15_20Sig
Definition: IsolatedTracksNxN.h:265
L1GtUtils::getL1GtRunCache
void getL1GtRunCache(const edm::Run &, const edm::EventSetup &, const bool, const bool)
get all the run-constant quantities for L1 trigger and cache them
Definition: L1GtUtils.cc:319
IsolatedTracksNxN::t_PVTracksSumPtHPWt
std::vector< double > * t_PVTracksSumPtHPWt
Definition: IsolatedTracksNxN.h:224
IsolatedTracksNxN::t_trackDxy
std::vector< double > * t_trackDxy
Definition: IsolatedTracksNxN.h:245
IsolatedTracksNxN::t_L1AlgoNames
std::vector< std::string > * t_L1AlgoNames
Definition: IsolatedTracksNxN.h:219
IsolatedTracksNxN::t_hsim5x5CharHad
std::vector< double > * t_hsim5x5CharHad
Definition: IsolatedTracksNxN.h:291
IsolatedTracksNxN::t_e15x15_10Sig
std::vector< double > * t_e15x15_10Sig
Definition: IsolatedTracksNxN.h:263
L1GtTriggerMenu::gtAlgorithmMap
const AlgorithmMap & gtAlgorithmMap() const
get / set the algorithm map (by name)
Definition: L1GtTriggerMenu.h:241
reco::HitPattern::stripTECLayersWithMeasurement
int stripTECLayersWithMeasurement() const
Definition: HitPattern.cc:615
reco::Vertex::trackRef_iterator
std::vector< TrackBaseRef >::const_iterator trackRef_iterator
The iteratator for the vector<TrackRef>
Definition: Vertex.h:37
IsolatedTracksNxN::h_maxNearP21x21
TH1F * h_maxNearP21x21[NPBins][NEtaBins]
Definition: IsolatedTracksNxN.h:193
reco::HitPattern::stripTECLayersWithoutMeasurement
int stripTECLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:683
IsolatedTracksNxN::t_jetPhi
std::vector< double > * t_jetPhi
Definition: IsolatedTracksNxN.h:236
IsolatedTracksNxN::t_trackL
std::vector< double > * t_trackL
Definition: IsolatedTracksNxN.h:256
IsolatedTracksNxN::t_esim7x7Rest
std::vector< double > * t_esim7x7Rest
Definition: IsolatedTracksNxN.h:275
IsolatedTracksNxN::t_trackDxyBS
std::vector< double > * t_trackDxyBS
Definition: IsolatedTracksNxN.h:245
reco::TrackBase::charge
int charge() const
track electric charge
Definition: TrackBase.h:567
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:62
IsolatedTracksNxN::t_trackHcalEta
std::vector< double > * t_trackHcalEta
Definition: IsolatedTracksNxN.h:244
IsolatedTracksNxN::t_L1METPt
std::vector< double > * t_L1METPt
Definition: IsolatedTracksNxN.h:234
gather_cfg.cout
cout
Definition: gather_cfg.py:145
dir
dbl *** dir
Definition: mlp_gen.cc:35
tree
Definition: tree.py:1
IsolatedTracksNxN::t_NLayersCrossed
std::vector< int > * t_NLayersCrossed
Definition: IsolatedTracksNxN.h:250
IsolatedTracksNxN::tok_jets_
edm::EDGetTokenT< reco::CaloJetCollection > tok_jets_
Definition: IsolatedTracksNxN.h:167
IsolatedTracksNxN::DeltaR
double DeltaR(double eta1, double phi1, double eta2, double phi2)
Definition: IsolatedTracksNxN.cc:1694
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)
IsolatedTracksNxN::t_trackHitInMeasTOB
std::vector< int > * t_trackHitInMeasTOB
Definition: IsolatedTracksNxN.h:254
IsolatedTracksNxN::t_PVx
std::vector< double > * t_PVx
Definition: IsolatedTracksNxN.h:223
IsolatedTracksNxN::t_trackHitInMeasTEC
std::vector< int > * t_trackHitInMeasTEC
Definition: IsolatedTracksNxN.h:254
IsolatedTracksNxN::beginJob
void beginJob() override
Definition: IsolatedTracksNxN.cc:949
IsolatedTracksNxN::t_esim15x15Photon
std::vector< double > * t_esim15x15Photon
Definition: IsolatedTracksNxN.h:276
IsolatedTracksNxN::t_e9x9
std::vector< double > * t_e9x9
Definition: IsolatedTracksNxN.h:262
IsolatedTracksNxN::t_hsim5x5Matched
std::vector< double > * t_hsim5x5Matched
Definition: IsolatedTracksNxN.h:287
L1GtUtils::ptrL1TriggerMenuEventSetup
const L1GtTriggerMenu * ptrL1TriggerMenuEventSetup(int &errorCode)
return a pointer to the L1 trigger menu from event setup
Definition: L1GtUtils.cc:1883
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:591
IsolatedTracksNxN::t_esim15x15Matched
std::vector< double > * t_esim15x15Matched
Definition: IsolatedTracksNxN.h:274
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:99
IsolatedTracksNxN::t_trackHitOutMeasTID
std::vector< int > * t_trackHitOutMeasTID
Definition: IsolatedTracksNxN.h:255
IsolatedTracksNxN::t_L1METEta
std::vector< double > * t_L1METEta
Definition: IsolatedTracksNxN.h:234
IsolatedTracksNxN::t_e7x7_30Sig
std::vector< double > * t_e7x7_30Sig
Definition: IsolatedTracksNxN.h:267
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
IsolatedTracksNxN::t_trackHitOutMissTEC
std::vector< int > * t_trackHitOutMissTEC
Definition: IsolatedTracksNxN.h:253
IsolatedTracksNxN::t_e11x11_25Sig
std::vector< double > * t_e11x11_25Sig
Definition: IsolatedTracksNxN.h:266
IsolatedTracksNxN::t_hsim3x3
std::vector< double > * t_hsim3x3
Definition: IsolatedTracksNxN.h:286
IsolatedTracksNxN::t_trackP
std::vector< double > * t_trackP
Definition: IsolatedTracksNxN.h:243
IsolatedTracksNxN::NPBins
static const size_t NPBins
Definition: IsolatedTracksNxN.h:190
IsolatedTracksNxN::t_hsim7x7Photon
std::vector< double > * t_hsim7x7Photon
Definition: IsolatedTracksNxN.h:289
IsolatedTracksNxN::t_jetPt
std::vector< double > * t_jetPt
Definition: IsolatedTracksNxN.h:236
useL1EventSetup
static const bool useL1EventSetup(true)
IsolatedTracksNxN::t_L1IsoEMPt
std::vector< double > * t_L1IsoEMPt
Definition: IsolatedTracksNxN.h:232
IsolatedTracksNxN::t_esim11x11Photon
std::vector< double > * t_esim11x11Photon
Definition: IsolatedTracksNxN.h:276
reco::HitPattern::stripTIBLayersWithoutMeasurement
int stripTIBLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:650
IsolatedTracksNxN::t_h5x5Sig
std::vector< double > * t_h5x5Sig
Definition: IsolatedTracksNxN.h:282
IsolatedTracksNxN::t_trackHitInMissTID
std::vector< int > * t_trackHitInMissTID
Definition: IsolatedTracksNxN.h:252
IsolatedTracksNxN::t_e11x11_30Sig
std::vector< double > * t_e11x11_30Sig
Definition: IsolatedTracksNxN.h:267
IsolatedTracksNxN::t_trackHitInMissTIBTID
std::vector< int > * t_trackHitInMissTIBTID
Definition: IsolatedTracksNxN.h:252
IsolatedTracksNxN::t_esim9x9CharHad
std::vector< double > * t_esim9x9CharHad
Definition: IsolatedTracksNxN.h:278
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::tok_EE_
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: IsolatedTracksNxN.h:175
IsolatedTracksNxN::tok_L1extTauJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extTauJet_
Definition: IsolatedTracksNxN.h:159
IsolatedTracksNxN::t_esim15x15CharHad
std::vector< double > * t_esim15x15CharHad
Definition: IsolatedTracksNxN.h:278
L1GtUtils::prescaleFactor
const int prescaleFactor(const edm::Event &iEvent, const std::string &nameAlgoTechTrig, int &errorCode) const
return prescale factor for a given algorithm or technical trigger
Definition: L1GtUtils.cc:1223
Generated for CMSSW Reference Manual by doxygen 1.8.11