CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
IsolatedTracksCone.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: IsolatedTracksCone
4 // Class: IsolatedTracksCone
5 //
6 
15 //
16 // Original Author: Jim Hirschauer (adaptation of Seema Sharma's
17 // IsolatedTracksNew)
18 // Created: Thu Nov 6 15:30:40 CST 2008
19 //
20 //
21 #define _DEBUG_QUIET
22 
23 #include "Calibration/IsolatedParticles/plugins/IsolatedTracksCone.h"
24 #include "Calibration/IsolatedParticles/interface/ChargeIsolation.h"
25 
26 #include "DataFormats/HLTReco/interface/TriggerObject.h"
27 #include "FWCore/Common/interface/TriggerNames.h"
28 #include "FWCore/Framework/interface/LuminosityBlock.h"
29 #include "DataFormats/Common/interface/TriggerResults.h"
30 #include "DataFormats/HLTReco/interface/TriggerEvent.h"
31 #include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
32 
33 #include "CondFormats/DataRecord/interface/EcalChannelStatusRcd.h"
34 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
35 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"
36 
37 IsolatedTracksCone::IsolatedTracksCone(const edm::ParameterSet& iConfig) {
38 
39  //now do what ever initialization is needed
40  doMC = iConfig.getUntrackedParameter<bool> ("DoMC", false);
41  myverbose_ =
42  iConfig.getUntrackedParameter<int>( "Verbosity", 5 );
43  useJetTrigger_ =
44  iConfig.getUntrackedParameter<bool>( "useJetTrigger", false);
45  drLeadJetVeto_ =
46  iConfig.getUntrackedParameter<double>( "drLeadJetVeto", 1.2 );
47  ptMinLeadJet_ =
48  iConfig.getUntrackedParameter<double>( "ptMinLeadJet", 15.0 );
49 
50  debugTrks_ =
51  iConfig.getUntrackedParameter<int>("DebugTracks");
52  printTrkHitPattern_ =
53  iConfig.getUntrackedParameter<bool>("PrintTrkHitPattern");
54 
55  minTrackP_ =
56  iConfig.getUntrackedParameter<double>( "minTrackP", 10.0);
57  maxTrackEta_ =
58  iConfig.getUntrackedParameter<double>( "maxTrackEta", 5.0);
59  maxNearTrackP_ =
60  iConfig.getUntrackedParameter<double>( "maxNearTrackP", 1.0);
61 
62  debugEcalSimInfo_ =
63  iConfig.getUntrackedParameter<int>("DebugEcalSimInfo");
64 
65  applyEcalIsolation_ =
66  iConfig.getUntrackedParameter<bool>("ApplyEcalIsolation");
67 
68  _L1extraTauJetSource =
69  iConfig.getParameter<edm::InputTag>("L1extraTauJetSource");
70  _L1extraCenJetSource =
71  iConfig.getParameter<edm::InputTag>("L1extraCenJetSource");
72  _L1extraFwdJetSource =
73  iConfig.getParameter<edm::InputTag>("L1extraFwdJetSource");
74 
75  edm::ParameterSet parameters =
76  iConfig.getParameter<edm::ParameterSet>("TrackAssociatorParameters");
77  parameters_.loadParameters( parameters );
78  trackAssociator_ = new TrackDetectorAssociator();
79  trackAssociator_->useDefaultPropagator();
80 
81  if(myverbose_>=0) {
82  std::cout <<"Parameters read from config file \n"
83  << "myverbose_ " << myverbose_ << "\n"
84  << "useJetTrigger_ " << useJetTrigger_ << "\n"
85  << "drLeadJetVeto_ " << drLeadJetVeto_ << "\n"
86  << "minTrackP_ " << minTrackP_ << "\n"
87  << "maxTrackEta_ " << maxTrackEta_ << "\n"
88  << "maxNearTrackP_ " << maxNearTrackP_
89  << std::endl;
90  }
91 }
92 
93 IsolatedTracksCone::~IsolatedTracksCone() {
94  delete trackAssociator_;
95 }
96 
97 void IsolatedTracksCone::analyze(const edm::Event& iEvent,
98  const edm::EventSetup& iSetup) {
99 
100  unsigned int irun = (unsigned int)iEvent.id().run();
101  unsigned int ilum = (unsigned int)iEvent.getLuminosityBlock().luminosityBlock();
102  unsigned int ievt = (unsigned int)iEvent.id().event();
103 
104 
105 
106  clearTrackVectors();
107 
108  // check the L1 objects
109  bool L1Pass = false;
110  leadL1JetPT=-999, leadL1JetEta=-999, leadL1JetPhi=-999;
111  if( !useJetTrigger_ ) {
112  L1Pass = true;
113  } else {
114  edm::Handle<l1extra::L1JetParticleCollection> l1TauHandle;
115  iEvent.getByLabel(_L1extraTauJetSource,l1TauHandle);
116  l1extra::L1JetParticleCollection::const_iterator itr;
117  for(itr = l1TauHandle->begin(); itr != l1TauHandle->end(); ++itr )
118  {
119  if( itr->pt()>leadL1JetPT ) {
120  leadL1JetPT = itr->pt();
121  leadL1JetEta = itr->eta();
122  leadL1JetPhi = itr->phi();
123  }
124  }
125  edm::Handle<l1extra::L1JetParticleCollection> l1CenJetHandle;
126  iEvent.getByLabel(_L1extraCenJetSource,l1CenJetHandle);
127  for( itr = l1CenJetHandle->begin(); itr != l1CenJetHandle->end(); ++itr )
128  {
129  if( itr->pt()>leadL1JetPT ) {
130  leadL1JetPT = itr->pt();
131  leadL1JetEta = itr->eta();
132  leadL1JetPhi = itr->phi();
133  }
134  }
135  edm::Handle<l1extra::L1JetParticleCollection> l1FwdJetHandle;
136  iEvent.getByLabel(_L1extraFwdJetSource,l1FwdJetHandle);
137  for( itr = l1FwdJetHandle->begin(); itr != l1FwdJetHandle->end(); ++itr )
138  {
139  if( itr->pt()>leadL1JetPT ) {
140  leadL1JetPT = itr->pt();
141  leadL1JetEta = itr->eta();
142  leadL1JetPhi = itr->phi();
143  }
144  }
145  if(leadL1JetPT>ptMinLeadJet_)
146  {
147  L1Pass = true;
148  }
149 
150  }
151 
152 
154  // Break now if L1Pass is false
156  if (!L1Pass) {
157  nEVT_failL1++;
158  // // std::cout << "L1Pass is false : " << L1Pass << std::endl;
159  // return;
160  }
161 
163  // Get the collection handles
165 
166 
167  edm::ESHandle<CaloGeometry> pG;
168  iSetup.get<CaloGeometryRecord>().get(pG);
169  const CaloGeometry* geo = pG.product();
170  const CaloSubdetectorGeometry* gEB =
171  geo->getSubdetectorGeometry(DetId::Ecal,EcalBarrel);
172  const CaloSubdetectorGeometry* gEE =
173  geo->getSubdetectorGeometry(DetId::Ecal,EcalEndcap);
174  const CaloSubdetectorGeometry* gHB =
175  geo->getSubdetectorGeometry(DetId::Hcal,HcalBarrel);
176  const CaloSubdetectorGeometry* gHE =
177  geo->getSubdetectorGeometry(DetId::Hcal,HcalEndcap);
178 
179  edm::ESHandle<CaloTopology> theCaloTopology;
180  iSetup.get<CaloTopologyRecord>().get(theCaloTopology);
181  const CaloTopology *caloTopology = theCaloTopology.product();
182  // const CaloSubdetectorTopology* theEBTopology = theCaloTopology->getSubdetectorTopology(DetId::Ecal,EcalBarrel);
183  // const CaloSubdetectorTopology* theEETopology = theCaloTopology->getSubdetectorTopology(DetId::Ecal,EcalEndcap);
184 
185  edm::ESHandle<HcalTopology> htopo;
186  iSetup.get<IdealGeometryRecord>().get(htopo);
187  const HcalTopology* theHBHETopology = htopo.product();
188 
189  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
190  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
191  iEvent.getByLabel("ecalRecHit","EcalRecHitsEB",barrelRecHitsHandle);
192  iEvent.getByLabel("ecalRecHit","EcalRecHitsEE",endcapRecHitsHandle);
193 
194  // Retrieve the good/bad ECAL channels from the DB
195  edm::ESHandle<EcalChannelStatus> ecalChStatus;
196  iSetup.get<EcalChannelStatusRcd>().get(ecalChStatus);
197  const EcalChannelStatus* theEcalChStatus = ecalChStatus.product();
198 
199  edm::Handle<HBHERecHitCollection> hbhe;
200  iEvent.getByLabel("hbhereco",hbhe);
201  const HBHERecHitCollection Hithbhe = *(hbhe.product());
202 
203  edm::Handle<reco::TrackCollection> trkCollection;
204  iEvent.getByLabel("generalTracks", trkCollection);
205  reco::TrackCollection::const_iterator trkItr;
206  if(debugTrks_>1){
207  std::cout << "Track Collection: " << std::endl;
208  std::cout << "Number of Tracks " << trkCollection->size() << std::endl;
209  }
210  std::string theTrackQuality = "highPurity";
211  reco::TrackBase::TrackQuality trackQuality_=
212  reco::TrackBase::qualityByName(theTrackQuality);
213 
214  //get Handles to SimTracks and SimHits
215  edm::Handle<edm::SimTrackContainer> SimTk;
216  if (doMC) iEvent.getByLabel("g4SimHits",SimTk);
217  edm::SimTrackContainer::const_iterator simTrkItr;
218 
219  edm::Handle<edm::SimVertexContainer> SimVtx;
220  if (doMC) iEvent.getByLabel("g4SimHits",SimVtx);
221 
222  //get Handles to PCaloHitContainers of eb/ee/hbhe
223  edm::Handle<edm::PCaloHitContainer> pcaloeb;
224  if (doMC) iEvent.getByLabel("g4SimHits", "EcalHitsEB", pcaloeb);
225 
226  edm::Handle<edm::PCaloHitContainer> pcaloee;
227  if (doMC) iEvent.getByLabel("g4SimHits", "EcalHitsEE", pcaloee);
228 
229  edm::Handle<edm::PCaloHitContainer> pcalohh;
230  if (doMC) iEvent.getByLabel("g4SimHits", "HcalHits", pcalohh);
231 
232 
233 
235  // Get HLT_IsoTrackHB/HE Information
237 
238  edm::InputTag theTriggerResultsLabel ("TriggerResults","","HLT");
239  edm::Handle<edm::TriggerResults> triggerResults;
240  iEvent.getByLabel( theTriggerResultsLabel, triggerResults);
241 
242 
243 
244  std::vector<int> v_hlTriggers;
245  int hltHB(-99);
246  int hltHE(-99);
247  int hltL1Jet15 (-99);
248  int hltJet30 (-99);
249  int hltJet50 (-99);
250  int hltJet80 (-99);
251  int hltJet110 (-99);
252  int hltJet140 (-99);
253  int hltJet180 (-99);
254  int hltL1SingleEG5 (-99);
255  int hltZeroBias (-99);
256  int hltMinBiasHcal (-99);
257  int hltMinBiasEcal (-99);
258  int hltMinBiasPixel (-99);
259  int hltSingleIsoTau30_Trk5 (-99);
260  int hltDoubleLooseIsoTau15_Trk5(-99);
261 
262  if (triggerResults.isValid()) {
263 
264  const edm::TriggerNames & triggerNames = iEvent.triggerNames(*triggerResults);
265  // TriggerNames class triggerNames.init(*triggerResults);
266 
267 
268  for (unsigned int i=0; i<triggerResults->size(); i++){
269  // std::cout << "triggerNames.triggerName(" << i << ") = " << triggerNames.triggerName(i) << std::endl;
270  if (triggerNames.triggerName(i) == "HLT_IsoTrackHE_1E31") hltHE = triggerResults->accept(i);
271  if (triggerNames.triggerName(i) == "HLT_IsoTrackHB_1E31") hltHB = triggerResults->accept(i);
272  if (triggerNames.triggerName(i) == "HLT_L1Jet15" ) hltL1Jet15 = triggerResults->accept(i);
273  if (triggerNames.triggerName(i) == "HLT_Jet30" ) hltJet30 = triggerResults->accept(i);
274  if (triggerNames.triggerName(i) == "HLT_Jet50" ) hltJet50 = triggerResults->accept(i);
275  if (triggerNames.triggerName(i) == "HLT_Jet80" ) hltJet80 = triggerResults->accept(i);
276  if (triggerNames.triggerName(i) == "HLT_Jet110" ) hltJet110 = triggerResults->accept(i);
277  if (triggerNames.triggerName(i) == "HLT_Jet140" ) hltJet140 = triggerResults->accept(i);
278  if (triggerNames.triggerName(i) == "HLT_Jet180" ) hltJet180 = triggerResults->accept(i);
279  if (triggerNames.triggerName(i) == "HLT_L1SingleEG5" ) hltL1SingleEG5 = triggerResults->accept(i);
280  if (triggerNames.triggerName(i) == "HLT_ZeroBias" ) hltZeroBias = triggerResults->accept(i);
281  if (triggerNames.triggerName(i) == "HLT_MinBiasHcal" ) hltMinBiasHcal = triggerResults->accept(i);
282  if (triggerNames.triggerName(i) == "HLT_MinBiasEcal" ) hltMinBiasEcal = triggerResults->accept(i);
283  if (triggerNames.triggerName(i) == "HLT_MinBiasPixel" ) hltMinBiasPixel = triggerResults->accept(i);
284  if (triggerNames.triggerName(i) == "HLT_SingleIsoTau30_Trk5" ) hltSingleIsoTau30_Trk5 = triggerResults->accept(i);
285  if (triggerNames.triggerName(i) == "HLT_DoubleLooseIsoTau15_Trk5" ) hltDoubleLooseIsoTau15_Trk5 = triggerResults->accept(i);
286  }
287  }
288 
289 
290 
291 
293  // Primary loop over tracks
295  TrackerHitAssociator* associate=0;
296  if (doMC) associate = new TrackerHitAssociator(iEvent);
297 
298 
299  nTRK = 0;
300  nRawTRK = 0;
301  nFailPt = 0;
302  nFailEta = 0;
303  nFailHighPurityQaul = 0;
304  nMissEcal = 0;
305  nMissHcal = 0;
306 
307  for( trkItr = trkCollection->begin();
308  trkItr != trkCollection->end(); ++trkItr)
309  {
310 
311  nRawTRK++;
312 
313  const reco::Track* pTrack = &(*trkItr);
314 
316  // Check for min Pt and max Eta P
318 
319  bool trkQual = pTrack->quality(trackQuality_);
320  bool goodPt = pTrack->p()>minTrackP_;
321  bool goodEta = std::abs(pTrack->momentum().eta())<maxTrackEta_;
322 
323  double eta1 = pTrack->momentum().eta();
324  double phi1 = pTrack->momentum().phi();
325  double pt1 = pTrack->pt();
326  double p1 = pTrack->p();
327 
328 
329  if (!goodEta){
330  nFailEta++;
331  }
332  if (!goodPt){
333  nFailPt++;
334  }
335  if (!trkQual){
336  nFailHighPurityQaul++;
337  }
338 
339  hRawPt ->Fill(pt1 );
340  hRawP ->Fill(p1 );
341  hRawEta->Fill(eta1);
342  hRawPhi->Fill(phi1);
343 
344  if( !goodEta || !goodPt || !trkQual ) continue; // Skip to next track
345 
347  // Find track trajectory
349 
350 
351  const FreeTrajectoryState fts1 =
352  trackAssociator_->getFreeTrajectoryState(iSetup, *pTrack);
353 
354 
355  TrackDetMatchInfo info1 =
356  trackAssociator_->associate(iEvent, iSetup, fts1, parameters_);
357 
358 
359 
361  // First confirm track makes it to Hcal
363 
364  if (info1.trkGlobPosAtHcal.x()==0 &&
365  info1.trkGlobPosAtHcal.y()==0 &&
366  info1.trkGlobPosAtHcal.z()==0)
367  {
368  nMissHcal++;
369  continue;
370  }
371 
372  const GlobalPoint hpoint1(info1.trkGlobPosAtHcal.x(),
373  info1.trkGlobPosAtHcal.y(),
374  info1.trkGlobPosAtHcal.z());
375 
376 
377 
379  // Get basic quantities
381 
382  const reco::HitPattern& hitp = pTrack->hitPattern();
383  int nLayersCrossed = hitp.trackerLayersWithMeasurement();
384  int nOuterHits = hitp.stripTOBLayersWithMeasurement()
385  +hitp.stripTECLayersWithMeasurement() ;
386 
387 
388  double simP = 0;
389  if (doMC) {
390  edm::SimTrackContainer::const_iterator matchedSimTrk =
391  spr::matchedSimTrack(iEvent, SimTk, SimVtx, pTrack, *associate, false);
392  simP = matchedSimTrk->momentum().P();
393  }
395  // Get Ecal Point
397 
398  const GlobalPoint point1(info1.trkGlobPosAtEcal.x(),
399  info1.trkGlobPosAtEcal.y(),
400  info1.trkGlobPosAtEcal.z());
401 
402 
403  // Sanity check that track hits Ecal
404 
405  if (info1.trkGlobPosAtEcal.x()==0 &&
406  info1.trkGlobPosAtEcal.y()==0 &&
407  info1.trkGlobPosAtEcal.z()==0)
408  {
409  std::cout << "Track doesn't reach Ecal." << std::endl;
410  nMissEcal++;
411  continue;
412  }
413 
414  // Get Track Momentum - make sure you have latest version of TrackDetMatchInfo
415 
416  GlobalVector trackMomAtEcal = info1.trkMomAtEcal;
417  GlobalVector trackMomAtHcal = info1.trkMomAtHcal;
418 
420  // If using Jet trigger, get distance from leading jet
422 
423  double drFromLeadJet = 999.0;
424  if( useJetTrigger_ ) {
425  double dphi = DeltaPhi(phi1, leadL1JetPhi);
426  double deta = eta1 - leadL1JetEta;
427  drFromLeadJet = sqrt(dphi*dphi + deta*deta);
428  }
429 
430 
432  // Define Arrays for sizes of Charge, Neut Iso Radii and
433  // Clustering Cone Radius.
435 
436  const int a_size = 7;
437  double a_coneR[a_size];
438  double a_charIsoR[a_size];
439  double a_neutIsoR[a_size];
440 
441  a_coneR[0] = 17.49; // = area of 2x2
442  a_coneR[1] = 26.23; // = area of 3x3
443  a_coneR[2] = 30.61;
444  a_coneR[3] = 34.98; // = area of 4x4
445  a_coneR[4] = 39.35;
446  a_coneR[5] = 43.72; // = area of 5x5
447  a_coneR[6] = 52.46; // = area of 6x6
448 
449  for (int i=0; i<a_size; i++){
450  a_charIsoR[i] = a_coneR[i]+28.9; // 28.9 gives 55.1 for 3x3 benchmark
451  a_neutIsoR[i] = a_charIsoR[i]*0.726; // Ecal radius = 0.726*Hcal radius
452  }
453 
455  // Do Neutral Iso in radius on Ecal surface.
457 
458  // NxN cluster
459  double e3x3=-999.0;
460  double trkEcalEne =-999.0;
461  edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
462  iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
463 
464  if(std::abs(point1.eta())<1.479) {
465  const DetId isoCell = gEB->getClosestCell(point1);
466  e3x3 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, sevlv.product(),1,1).first;
467  trkEcalEne = spr::eCaloSimInfo(iEvent, geo, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate);
468  } else {
469  const DetId isoCell = gEE->getClosestCell(point1);
470  e3x3 = spr::eECALmatrix(isoCell,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology, sevlv.product(),1,1).first;
471  trkEcalEne = spr::eCaloSimInfo(iEvent, geo, pcaloeb, pcaloee, SimTk, SimVtx, pTrack, *associate);
472  }
473 
474  // Cone cluster
475 
476  // Set up array of cone sizes for MIP cut
477  const int a_mip_size = 5;
478  double a_mipR[a_mip_size];
479  a_mipR[0] = 3.84; // = area of 3x3 ecal
480  a_mipR[1] = 14.0;
481  a_mipR[2] = 19.0;
482  a_mipR[3] = 24.0;
483  a_mipR[4] = 9.0; // = from standard analyzer
484 
485  std::vector<double> v_eDR;
486  for (int i = 0 ; i < a_size ; i++){
487  int nRH_eDR = 0;
488 
489  // Cone in ecal
490  double eDR = spr::eCone_ecal(geo,
491  barrelRecHitsHandle,
492  endcapRecHitsHandle,
493  hpoint1, point1,
494  a_neutIsoR[i],
495  trackMomAtEcal, nRH_eDR);
496  v_eDR.push_back(eDR);
497 
498  }
499 
500  std::vector<double> v_eMipDR;
501  for (int i = 0 ; i < a_mip_size ; i++){
502  int nRH_eMipDR = 0;
503  double eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle,
504  endcapRecHitsHandle,
505  hpoint1, point1,
506  a_mipR[i], trackMomAtEcal, nRH_eMipDR);
507 
508  v_eMipDR.push_back(eMipDR);
509  }
510 
511 
513  // Do charge isolation in radius at Hcal surface for 5 different
514  // radii defined above in a_charIso
516 
517  std::vector<double> v_hmaxNearP_goodTrk;
518  std::vector<double> v_hmaxNearP ;
519  std::vector<int> v_hnNearTRKs ;
520  std::vector<int> v_hnLayers_maxNearP;
521  std::vector<int> v_htrkQual_maxNearP;
522 
523  std::vector<double> v_cone_hmaxNearP_goodTrk;
524  std::vector<double> v_cone_hmaxNearP ;
525  std::vector<int> v_cone_hnNearTRKs ;
526  std::vector<int> v_cone_hnLayers_maxNearP;
527  std::vector<int> v_cone_htrkQual_maxNearP;
528 
529  for (int i = 0 ; i < a_size ; i++){
530 
531  double hmaxNearP = -999.0;
532  int hnNearTRKs = 0;
533  int hnLayers_maxNearP = 0;
534  int htrkQual_maxNearP = -1;
535  double hmaxNearP_goodTrk = -999.0;
536 
537  double conehmaxNearP = -999.0;
538  int conehnNearTRKs = 0;
539  int conehnLayers_maxNearP = 0;
540  int conehtrkQual_maxNearP = -1;
541  double conehmaxNearP_goodTrk = -999.0;
542 
543  conehmaxNearP = spr::coneChargeIsolation(iEvent, iSetup,
544  trkItr, trkCollection,
545  *trackAssociator_, parameters_,
546  theTrackQuality, conehnNearTRKs,
547  conehnLayers_maxNearP,
548  conehtrkQual_maxNearP,
549  conehmaxNearP_goodTrk,
550  hpoint1, trackMomAtHcal,
551  a_charIsoR[i]);
552 
553  v_hmaxNearP_goodTrk.push_back(hmaxNearP_goodTrk);
554  v_hmaxNearP .push_back(hmaxNearP );
555  v_hnNearTRKs .push_back(hnNearTRKs );
556  v_hnLayers_maxNearP.push_back(hnLayers_maxNearP);
557  v_htrkQual_maxNearP.push_back(htrkQual_maxNearP);
558 
559  v_cone_hmaxNearP_goodTrk.push_back(conehmaxNearP_goodTrk);
560  v_cone_hmaxNearP .push_back(conehmaxNearP );
561  v_cone_hnNearTRKs .push_back(conehnNearTRKs );
562  v_cone_hnLayers_maxNearP.push_back(conehnLayers_maxNearP);
563  v_cone_htrkQual_maxNearP.push_back(conehtrkQual_maxNearP);
564 
565  }
566 
567  double h3x3=-999.0, h5x5=-999.0;
568  double hsim3x3=-999.0, hsim5x5=-999.0, trkHcalEne=-999.0;
569  std::map<std::string, double> hsimInfo3x3, hsimInfo5x5;
570  double distFromHotCell_h3x3 = -99.;
571  int ietaFromHotCell_h3x3 = -99;
572  int iphiFromHotCell_h3x3 = -99;
573  double distFromHotCell_h5x5 = -99.;
574  int ietaFromHotCell_h5x5 = -99;
575  int iphiFromHotCell_h5x5 = -99;
576 
577  GlobalPoint gPosHotCell_h3x3(0.,0.,0.);
578  GlobalPoint gPosHotCell_h5x5(0.,0.,0.);
579 
580  int nRH_h3x3(0), nRH_h5x5(0);
581 
582  // Hcal Energy Clustering
583 
584  // Get closetcell for ietaFromHotCell and iphiFromHotCell
585  DetId ClosestCell;
586  if( std::abs(pTrack->eta())<1.392 ) {
587  ClosestCell = gHB->getClosestCell(hpoint1);
588  } else {
589  ClosestCell = gHE->getClosestCell(hpoint1);
590  }
591  // Transform into HcalDetId so that I can get ieta, iphi later.
592  HcalDetId ClosestCell_HcalDetId(ClosestCell.rawId());
593 
594  // Using NxN Cluster
595  std::vector<int> v_RH_h3x3_ieta;
596  std::vector<int> v_RH_h3x3_iphi;
597  std::vector<double> v_RH_h3x3_ene;
598  std::vector<int> v_RH_h5x5_ieta;
599  std::vector<int> v_RH_h5x5_iphi;
600  std::vector<double> v_RH_h5x5_ene;
601 
602 
603  h3x3 = spr::eHCALmatrix(geo, theHBHETopology, ClosestCell, hbhe,1,1,
604  nRH_h3x3, v_RH_h3x3_ieta, v_RH_h3x3_iphi, v_RH_h3x3_ene,
605  gPosHotCell_h3x3);
606  distFromHotCell_h3x3 = spr::getDistInPlaneTrackDir(hpoint1, trackMomAtHcal, gPosHotCell_h3x3);
607 
608  h5x5 = spr::eHCALmatrix(geo, theHBHETopology, ClosestCell, hbhe,2,2,
609  nRH_h5x5, v_RH_h5x5_ieta, v_RH_h5x5_iphi, v_RH_h5x5_ene,
610  gPosHotCell_h5x5);
611  distFromHotCell_h5x5 = spr::getDistInPlaneTrackDir(hpoint1, trackMomAtHcal, gPosHotCell_h5x5);
612 
613 
614  // double heta = (double)hpoint1.eta();
615  // double hphi = (double)hpoint1.phi();
616  std::vector<int> multiplicity3x3;
617  std::vector<int> multiplicity5x5;
618  if (doMC) {
619  hsim3x3 = spr::eHCALmatrix(theHBHETopology, ClosestCell,
620  pcalohh,1,1);
621  hsim5x5 = spr::eHCALmatrix(theHBHETopology, ClosestCell,
622  pcalohh,2,2);
623 
624  hsimInfo3x3 = spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, pcalohh, SimTk, SimVtx, pTrack, *associate, 1,1, multiplicity3x3);
625  hsimInfo5x5 = spr::eHCALSimInfo(iEvent, theHBHETopology, ClosestCell, pcalohh, SimTk, SimVtx, pTrack, *associate, 2,2, multiplicity5x5);
626 
627  // Get energy from all simhits in hcal associated with iso track
628  trkHcalEne = spr::eCaloSimInfo(iEvent, geo, pcalohh, SimTk, SimVtx, pTrack, *associate);
629  }
630 
631  // Finally for cones of varying radii.
632  std::vector<double> v_hsimInfoConeMatched;
633  std::vector<double> v_hsimInfoConeRest ;
634  std::vector<double> v_hsimInfoConePhoton ;
635  std::vector<double> v_hsimInfoConeNeutHad;
636  std::vector<double> v_hsimInfoConeCharHad;
637  std::vector<double> v_hsimInfoConePdgMatched;
638  std::vector<double> v_hsimInfoConeTotal ;
639 
640  std::vector<int> v_hsimInfoConeNMatched;
641  std::vector<int> v_hsimInfoConeNTotal ;
642  std::vector<int> v_hsimInfoConeNNeutHad;
643  std::vector<int> v_hsimInfoConeNCharHad;
644  std::vector<int> v_hsimInfoConeNPhoton ;
645  std::vector<int> v_hsimInfoConeNRest ;
646 
647  std::vector<double> v_hsimCone ;
648  std::vector<double> v_hCone ;
649 
650  std::vector<int> v_nRecHitsCone ;
651  std::vector<int> v_nSimHitsCone ;
652 
653  std::vector<double> v_distFromHotCell;
654  std::vector<int> v_ietaFromHotCell;
655  std::vector<int> v_iphiFromHotCell;
656  GlobalPoint gposHotCell(0.,0.,0.);
657 
658 
659  std::vector<int> v_RH_r26_ieta;
660  std::vector<int> v_RH_r26_iphi;
661  std::vector<double> v_RH_r26_ene;
662  std::vector<int> v_RH_r44_ieta;
663  std::vector<int> v_RH_r44_iphi;
664  std::vector<double> v_RH_r44_ene;
665 
666 
667 
668  for (int i = 0 ; i < a_size ; i++){
669 
670 
671  std::map<std::string, double> hsimInfoCone;
672  double hsimCone = -999.0, hCone = -999.0;
673  double distFromHotCell = -99.0;
674  int ietaFromHotCell = -99;
675  int iphiFromHotCell = -99;
676  int ietaHotCell = -99;
677  int iphiHotCell = -99;
678  int nRecHitsCone = -999;
679  int nSimHitsCone = -999;
680 
681  std::vector<int> multiplicityCone;
682  std::vector<DetId> coneRecHitDetIds;
683  if (doMC)
684  hsimCone = spr::eCone_hcal(geo, pcalohh, hpoint1, point1,
685  a_coneR[i], trackMomAtHcal, nSimHitsCone);
686 
687  // If needed, get ieta and iphi of rechits for cones of 23.25
688  // and for hitmap for debugging
689  bool makeHitmaps = false;
690  if (a_coneR[i] == 26.23 && makeHitmaps)
691  {
692 
693  hCone = spr::eCone_hcal(geo, hbhe, hpoint1, point1,
694  a_coneR[i], trackMomAtHcal,nRecHitsCone,
695  v_RH_r26_ieta, v_RH_r26_iphi, v_RH_r26_ene,
696  coneRecHitDetIds, distFromHotCell,
697  ietaHotCell, iphiHotCell, gposHotCell);
698  }
699  else if (a_coneR[i] == 43.72 && makeHitmaps)
700  {
701 
702  hCone = spr::eCone_hcal(geo, hbhe, hpoint1, point1,
703  a_coneR[i], trackMomAtHcal,nRecHitsCone,
704  v_RH_r44_ieta, v_RH_r44_iphi, v_RH_r44_ene,
705  coneRecHitDetIds, distFromHotCell,
706  ietaHotCell, iphiHotCell, gposHotCell);
707  }
708  else
709  {
710 
711  hCone = spr::eCone_hcal(geo, hbhe, hpoint1, point1,
712  a_coneR[i], trackMomAtHcal, nRecHitsCone,
713  coneRecHitDetIds, distFromHotCell,
714  ietaHotCell, iphiHotCell, gposHotCell);
715  }
716 
717 
718 
719  if (ietaHotCell != 99){
720  ietaFromHotCell = ietaHotCell-ClosestCell_HcalDetId.ieta();
721  iphiFromHotCell = iphiHotCell-ClosestCell_HcalDetId.iphi();
722  }
723 
724  // SimHits NOT matched to RecHits
725  if (doMC) {
726  hsimInfoCone = spr::eHCALSimInfoCone(iEvent,pcalohh, SimTk, SimVtx, pTrack, *associate, geo, hpoint1, point1, a_coneR[i], trackMomAtHcal, multiplicityCone);
727 
728 
729 
730  // SimHits matched to RecHits
731  // hsimInfoCone = spr::eHCALSimInfoCone(iEvent,pcalohh, SimTk, SimVtx,
732  // pTrack, *associate,
733  // geo, hpoint1, point1,
734  // a_coneR[i], trackMomAtHcal,
735  // multiplicityCone,
736  // coneRecHitDetIds);
737 
738  v_hsimInfoConeMatched .push_back(hsimInfoCone["eMatched" ]);
739  v_hsimInfoConeRest .push_back(hsimInfoCone["eRest" ]);
740  v_hsimInfoConePhoton .push_back(hsimInfoCone["eGamma" ]);
741  v_hsimInfoConeNeutHad .push_back(hsimInfoCone["eNeutralHad"]);
742  v_hsimInfoConeCharHad .push_back(hsimInfoCone["eChargedHad"]);
743  v_hsimInfoConePdgMatched.push_back(hsimInfoCone["pdgMatched" ]);
744  v_hsimInfoConeTotal .push_back(hsimInfoCone["eTotal" ]);
745 
746  v_hsimInfoConeNMatched .push_back(multiplicityCone.at(0));
747 
748  v_hsimInfoConeNTotal .push_back(multiplicityCone.at(1));
749  v_hsimInfoConeNNeutHad .push_back(multiplicityCone.at(2));
750  v_hsimInfoConeNCharHad .push_back(multiplicityCone.at(3));
751  v_hsimInfoConeNPhoton .push_back(multiplicityCone.at(4));
752  v_hsimInfoConeNRest .push_back(multiplicityCone.at(5));
753 
754  v_hsimCone .push_back(hsimCone );
755  v_nSimHitsCone .push_back(nSimHitsCone );
756  }
757  v_hCone .push_back(hCone );
758  v_nRecHitsCone .push_back(nRecHitsCone );
759 
760  v_distFromHotCell .push_back(distFromHotCell );
761  v_ietaFromHotCell .push_back(ietaFromHotCell );
762  v_iphiFromHotCell .push_back(iphiFromHotCell );
763 
764 
765  }
766 
767 
769  // Fill Vectors that go into root file
771 
772  t_v_hnNearTRKs ->push_back(v_hnNearTRKs );
773  t_v_hnLayers_maxNearP ->push_back(v_hnLayers_maxNearP );
774  t_v_htrkQual_maxNearP ->push_back(v_htrkQual_maxNearP );
775  t_v_hmaxNearP_goodTrk ->push_back(v_hmaxNearP_goodTrk );
776  t_v_hmaxNearP ->push_back(v_hmaxNearP );
777 
778  t_v_cone_hnNearTRKs ->push_back(v_cone_hnNearTRKs );
779  t_v_cone_hnLayers_maxNearP ->push_back(v_cone_hnLayers_maxNearP);
780  t_v_cone_htrkQual_maxNearP ->push_back(v_cone_htrkQual_maxNearP);
781  t_v_cone_hmaxNearP_goodTrk ->push_back(v_cone_hmaxNearP_goodTrk);
782  t_v_cone_hmaxNearP ->push_back(v_cone_hmaxNearP );
783 
784  // t_hScale ->push_back(hScale );
785  t_trkNOuterHits ->push_back(nOuterHits );
786  t_trkNLayersCrossed ->push_back(nLayersCrossed );
787  t_dtFromLeadJet ->push_back(drFromLeadJet );
788  t_trkP ->push_back(p1 );
789  t_trkPt ->push_back(pt1 );
790  t_trkEta ->push_back(eta1 );
791  t_trkPhi ->push_back(phi1 );
792 
793  t_e3x3 ->push_back(e3x3 );
794  t_v_eDR ->push_back(v_eDR );
795  t_v_eMipDR ->push_back(v_eMipDR );
796 
797  t_h3x3 ->push_back(h3x3 );
798  t_h5x5 ->push_back(h5x5 );
799  t_nRH_h3x3 ->push_back(nRH_h3x3 );
800  t_nRH_h5x5 ->push_back(nRH_h5x5 );
801 
802  t_v_RH_h3x3_ieta ->push_back(v_RH_h3x3_ieta);
803  t_v_RH_h3x3_iphi ->push_back(v_RH_h3x3_iphi);
804  t_v_RH_h3x3_ene ->push_back(v_RH_h3x3_ene);
805  t_v_RH_h5x5_ieta ->push_back(v_RH_h5x5_ieta);
806  t_v_RH_h5x5_iphi ->push_back(v_RH_h5x5_iphi);
807  t_v_RH_h5x5_ene ->push_back(v_RH_h5x5_ene);
808 
809  if (doMC) {
810  t_simP ->push_back(simP );
811  t_hsim3x3 ->push_back(hsim3x3 );
812  t_hsim5x5 ->push_back(hsim5x5 );
813 
814  t_hsim3x3Matched ->push_back(hsimInfo3x3["eMatched"] );
815  t_hsim5x5Matched ->push_back(hsimInfo5x5["eMatched"] );
816  t_hsim3x3Rest ->push_back(hsimInfo3x3["eRest"] );
817  t_hsim5x5Rest ->push_back(hsimInfo5x5["eRest"] );
818  t_hsim3x3Photon ->push_back(hsimInfo3x3["eGamma"] );
819  t_hsim5x5Photon ->push_back(hsimInfo5x5["eGamma"] );
820  t_hsim3x3NeutHad ->push_back(hsimInfo3x3["eNeutralHad"] );
821  t_hsim5x5NeutHad ->push_back(hsimInfo5x5["eNeutralHad"] );
822  t_hsim3x3CharHad ->push_back(hsimInfo3x3["eChargedHad"] );
823  t_hsim5x5CharHad ->push_back(hsimInfo5x5["eChargedHad"] );
824  t_hsim3x3Total ->push_back(hsimInfo3x3["eTotal"] );
825  t_hsim5x5Total ->push_back(hsimInfo5x5["eTotal"] );
826  t_hsim3x3PdgMatched ->push_back(hsimInfo3x3["pdgMatched"] );
827  t_hsim5x5PdgMatched ->push_back(hsimInfo5x5["pdgMatched"] );
828 
829  t_hsim3x3NMatched ->push_back(multiplicity3x3.at(0));
830  t_hsim3x3NTotal ->push_back(multiplicity3x3.at(1));
831  t_hsim3x3NNeutHad ->push_back(multiplicity3x3.at(2));
832  t_hsim3x3NCharHad ->push_back(multiplicity3x3.at(3));
833  t_hsim3x3NPhoton ->push_back(multiplicity3x3.at(4));
834  t_hsim3x3NRest ->push_back(multiplicity3x3.at(5));
835 
836  t_hsim5x5NMatched ->push_back(multiplicity5x5.at(0));
837  t_hsim5x5NTotal ->push_back(multiplicity5x5.at(1));
838  t_hsim5x5NNeutHad ->push_back(multiplicity5x5.at(2));
839  t_hsim5x5NCharHad ->push_back(multiplicity5x5.at(3));
840  t_hsim5x5NPhoton ->push_back(multiplicity5x5.at(4));
841  t_hsim5x5NRest ->push_back(multiplicity5x5.at(5));
842  }
843 
844  t_distFromHotCell_h3x3->push_back(distFromHotCell_h3x3);
845  t_ietaFromHotCell_h3x3->push_back(ietaFromHotCell_h3x3);
846  t_iphiFromHotCell_h3x3->push_back(iphiFromHotCell_h3x3);
847  t_distFromHotCell_h5x5->push_back(distFromHotCell_h5x5);
848  t_ietaFromHotCell_h5x5->push_back(ietaFromHotCell_h5x5);
849  t_iphiFromHotCell_h5x5->push_back(iphiFromHotCell_h5x5);
850 
851  if (doMC) {
852  t_trkHcalEne ->push_back(trkHcalEne );
853  t_trkEcalEne ->push_back(trkEcalEne );
854 
855  t_v_hsimInfoConeMatched ->push_back(v_hsimInfoConeMatched );
856  t_v_hsimInfoConeRest ->push_back(v_hsimInfoConeRest );
857  t_v_hsimInfoConePhoton ->push_back(v_hsimInfoConePhoton );
858  t_v_hsimInfoConeNeutHad ->push_back(v_hsimInfoConeNeutHad );
859  t_v_hsimInfoConeCharHad ->push_back(v_hsimInfoConeCharHad );
860  t_v_hsimInfoConePdgMatched ->push_back(v_hsimInfoConePdgMatched);
861  t_v_hsimInfoConeTotal ->push_back(v_hsimInfoConeTotal );
862 
863  t_v_hsimInfoConeNMatched ->push_back(v_hsimInfoConeNMatched );
864  t_v_hsimInfoConeNTotal ->push_back(v_hsimInfoConeNTotal );
865  t_v_hsimInfoConeNNeutHad ->push_back(v_hsimInfoConeNNeutHad );
866  t_v_hsimInfoConeNCharHad ->push_back(v_hsimInfoConeNCharHad );
867  t_v_hsimInfoConeNPhoton ->push_back(v_hsimInfoConeNPhoton );
868  t_v_hsimInfoConeNRest ->push_back(v_hsimInfoConeNRest );
869 
870  t_v_hsimCone ->push_back(v_hsimCone );
871  t_v_hCone ->push_back(v_hCone );
872  t_v_nRecHitsCone->push_back(v_nRecHitsCone);
873  t_v_nSimHitsCone->push_back(v_nSimHitsCone);
874  }
875 
876 
877  t_v_distFromHotCell->push_back(v_distFromHotCell);
878  t_v_ietaFromHotCell->push_back(v_ietaFromHotCell);
879  t_v_iphiFromHotCell->push_back(v_iphiFromHotCell);
880 
881  t_v_RH_r26_ieta ->push_back(v_RH_r26_ieta);
882  t_v_RH_r26_iphi ->push_back(v_RH_r26_iphi);
883  t_v_RH_r26_ene ->push_back(v_RH_r26_ene);
884  t_v_RH_r44_ieta ->push_back(v_RH_r44_ieta);
885  t_v_RH_r44_iphi ->push_back(v_RH_r44_iphi);
886  t_v_RH_r44_ene ->push_back(v_RH_r44_ene);
887 
888 
889 
890  t_v_hlTriggers ->push_back(v_hlTriggers);
891  t_hltHB ->push_back(hltHB);
892  t_hltHE ->push_back(hltHE);
893  t_hltL1Jet15 ->push_back(hltL1Jet15 );
894  t_hltJet30 ->push_back(hltJet30 );
895  t_hltJet50 ->push_back(hltJet50 );
896  t_hltJet80 ->push_back(hltJet80 );
897  t_hltJet110 ->push_back(hltJet110 );
898  t_hltJet140 ->push_back(hltJet140 );
899  t_hltJet180 ->push_back(hltJet180 );
900  t_hltL1SingleEG5 ->push_back(hltL1SingleEG5 );
901  t_hltZeroBias ->push_back(hltZeroBias );
902  t_hltMinBiasHcal ->push_back(hltMinBiasHcal );
903  t_hltMinBiasEcal ->push_back(hltMinBiasEcal );
904  t_hltMinBiasPixel ->push_back(hltMinBiasPixel );
905  t_hltSingleIsoTau30_Trk5 ->push_back(hltSingleIsoTau30_Trk5 );
906  t_hltDoubleLooseIsoTau15_Trk5->push_back(hltDoubleLooseIsoTau15_Trk5);
907 
908  t_irun->push_back(irun);
909  t_ievt->push_back(ievt);
910  t_ilum->push_back(ilum);
911 
912  nTRK++;
913 
914 
915  } // Loop over track collection
916 
917  // std::cout << "nEVT= " << nEVT << std::endl;
918 
919  ntp->Fill();
920  nEVT++;
921 
922 
923  delete associate;
924 }
925 
926 
927 
928 void IsolatedTracksCone::beginJob(const edm::EventSetup&) {
929 
930  // hbScale = 120.0;
931  // heScale = 135.0;
932  nEVT=0;
933  nEVT_failL1=0;
934  nTRK=0;
935 
936  double tempgen_TH[22] = { 0.0, 1.0, 2.0, 3.0, 4.0,
937  5.0, 6.0, 7.0, 8.0, 9.0,
938  10.0, 12.0, 15.0, 20.0, 25.0,
939  30.0, 40.0, 50.0, 60.0, 70.0, 80.0, 100};
940  for(int i=0; i<22; i++) genPartPBins[i] = tempgen_TH[i];
941 
942 
943  double tempgen_Eta[5] = {0.0, 0.5, 1.1, 1.7, 2.0};
944  for(int i=0; i<5; i++) genPartEtaBins[i] = tempgen_Eta[i];
945 
946  t_v_hnNearTRKs = new std::vector<std::vector<int> > ();
947  t_v_hnLayers_maxNearP = new std::vector<std::vector<int> > ();
948  t_v_htrkQual_maxNearP = new std::vector<std::vector<int> > ();
949  t_v_hmaxNearP_goodTrk = new std::vector<std::vector<double> >();
950  t_v_hmaxNearP = new std::vector<std::vector<double> >();
951 
952  t_v_cone_hnNearTRKs = new std::vector<std::vector<int> > ();
953  t_v_cone_hnLayers_maxNearP = new std::vector<std::vector<int> > ();
954  t_v_cone_htrkQual_maxNearP = new std::vector<std::vector<int> > ();
955  t_v_cone_hmaxNearP_goodTrk = new std::vector<std::vector<double> >();
956  t_v_cone_hmaxNearP = new std::vector<std::vector<double> >();
957 
958  // t_hScale = new std::vector<double>();
959  t_trkNOuterHits = new std::vector<double>();
960  t_trkNLayersCrossed = new std::vector<double>();
961  t_dtFromLeadJet = new std::vector<double>();
962  t_trkP = new std::vector<double>();
963  t_trkPt = new std::vector<double>();
964  t_trkEta = new std::vector<double>();
965  t_trkPhi = new std::vector<double>();
966 
967  t_e3x3 = new std::vector<double>();
968  t_v_eDR = new std::vector<std::vector<double> >();
969  t_v_eMipDR = new std::vector<std::vector<double> >();
970 
971  t_h3x3 = new std::vector<double>();
972  t_h5x5 = new std::vector<double>();
973 
974  t_nRH_h3x3 = new std::vector<double>();
975  t_nRH_h5x5 = new std::vector<double>();
976 
977  if (doMC) {
978  t_simP = new std::vector<double>();
979  t_hsim3x3 = new std::vector<double>();
980  t_hsim5x5 = new std::vector<double>();
981 
982  t_hsim3x3Matched = new std::vector<double>();
983  t_hsim5x5Matched = new std::vector<double>();
984  t_hsim3x3Rest = new std::vector<double>();
985  t_hsim5x5Rest = new std::vector<double>();
986  t_hsim3x3Photon = new std::vector<double>();
987  t_hsim5x5Photon = new std::vector<double>();
988  t_hsim3x3NeutHad = new std::vector<double>();
989  t_hsim5x5NeutHad = new std::vector<double>();
990  t_hsim3x3CharHad = new std::vector<double>();
991  t_hsim5x5CharHad = new std::vector<double>();
992  t_hsim3x3PdgMatched = new std::vector<double>();
993  t_hsim5x5PdgMatched = new std::vector<double>();
994  t_hsim3x3Total = new std::vector<double>();
995  t_hsim5x5Total = new std::vector<double>();
996 
997  t_hsim3x3NMatched = new std::vector<int>();
998  t_hsim3x3NTotal = new std::vector<int>();
999  t_hsim3x3NNeutHad = new std::vector<int>();
1000  t_hsim3x3NCharHad = new std::vector<int>();
1001  t_hsim3x3NPhoton = new std::vector<int>();
1002  t_hsim3x3NRest = new std::vector<int>();
1003 
1004  t_hsim5x5NMatched = new std::vector<int>();
1005  t_hsim5x5NTotal = new std::vector<int>();
1006  t_hsim5x5NNeutHad = new std::vector<int>();
1007  t_hsim5x5NCharHad = new std::vector<int>();
1008  t_hsim5x5NPhoton = new std::vector<int>();
1009  t_hsim5x5NRest = new std::vector<int>();
1010 
1011  t_trkHcalEne = new std::vector<double>();
1012  t_trkEcalEne = new std::vector<double>();
1013  }
1014 
1015  t_distFromHotCell_h3x3 = new std::vector<double>();
1016  t_ietaFromHotCell_h3x3 = new std::vector<int>();
1017  t_iphiFromHotCell_h3x3 = new std::vector<int>();
1018  t_distFromHotCell_h5x5 = new std::vector<double>();
1019  t_ietaFromHotCell_h5x5 = new std::vector<int>();
1020  t_iphiFromHotCell_h5x5 = new std::vector<int>();
1021 
1022  if (doMC) {
1023  t_v_hsimInfoConeMatched = new std::vector<std::vector<double> >();
1024  t_v_hsimInfoConeRest = new std::vector<std::vector<double> >();
1025  t_v_hsimInfoConePhoton = new std::vector<std::vector<double> >();
1026  t_v_hsimInfoConeNeutHad = new std::vector<std::vector<double> >();
1027  t_v_hsimInfoConeCharHad = new std::vector<std::vector<double> >();
1028  t_v_hsimInfoConePdgMatched= new std::vector<std::vector<double> >();
1029  t_v_hsimInfoConeTotal = new std::vector<std::vector<double> >();
1030 
1031  t_v_hsimInfoConeNMatched = new std::vector<std::vector<int> >();
1032  t_v_hsimInfoConeNTotal = new std::vector<std::vector<int> >();
1033  t_v_hsimInfoConeNNeutHad = new std::vector<std::vector<int> >();
1034  t_v_hsimInfoConeNCharHad = new std::vector<std::vector<int> >();
1035  t_v_hsimInfoConeNPhoton = new std::vector<std::vector<int> >();
1036  t_v_hsimInfoConeNRest = new std::vector<std::vector<int> >();
1037 
1038  t_v_hsimCone = new std::vector<std::vector<double> >();
1039  }
1040 
1041  t_v_hCone = new std::vector<std::vector<double> >();
1042  t_v_nRecHitsCone= new std::vector<std::vector<int> >();
1043  t_v_nSimHitsCone= new std::vector<std::vector<int> >();
1044 
1045  t_v_distFromHotCell= new std::vector<std::vector<double> >();
1046  t_v_ietaFromHotCell= new std::vector<std::vector<int> >();
1047  t_v_iphiFromHotCell= new std::vector<std::vector<int> >();
1048 
1049  t_v_RH_h3x3_ieta = new std::vector<std::vector<int> >();
1050  t_v_RH_h3x3_iphi = new std::vector<std::vector<int> >();
1051  t_v_RH_h3x3_ene = new std::vector<std::vector<double> >();
1052  t_v_RH_h5x5_ieta = new std::vector<std::vector<int> >();
1053  t_v_RH_h5x5_iphi = new std::vector<std::vector<int> >();
1054  t_v_RH_h5x5_ene = new std::vector<std::vector<double> >();
1055  t_v_RH_r26_ieta = new std::vector<std::vector<int> >();
1056  t_v_RH_r26_iphi = new std::vector<std::vector<int> >();
1057  t_v_RH_r26_ene = new std::vector<std::vector<double> >();
1058  t_v_RH_r44_ieta = new std::vector<std::vector<int> >();
1059  t_v_RH_r44_iphi = new std::vector<std::vector<int> >();
1060  t_v_RH_r44_ene = new std::vector<std::vector<double> >();
1061 
1062 
1063  t_v_hlTriggers = new std::vector<std::vector<int> >();
1064 
1065  t_hltHE = new std::vector<int>();
1066  t_hltHB = new std::vector<int>();
1067  t_hltL1Jet15 = new std::vector<int>();
1068  t_hltJet30 = new std::vector<int>();
1069  t_hltJet50 = new std::vector<int>();
1070  t_hltJet80 = new std::vector<int>();
1071  t_hltJet110 = new std::vector<int>();
1072  t_hltJet140 = new std::vector<int>();
1073  t_hltJet180 = new std::vector<int>();
1074  t_hltL1SingleEG5 = new std::vector<int>();
1075  t_hltZeroBias = new std::vector<int>();
1076  t_hltMinBiasHcal = new std::vector<int>();
1077  t_hltMinBiasEcal = new std::vector<int>();
1078  t_hltMinBiasPixel = new std::vector<int>();
1079  t_hltSingleIsoTau30_Trk5 = new std::vector<int>();
1080  t_hltDoubleLooseIsoTau15_Trk5 = new std::vector<int>();
1081 
1082 
1083  t_irun = new std::vector<unsigned int>();
1084  t_ievt = new std::vector<unsigned int>();
1085  t_ilum = new std::vector<unsigned int>();
1086 
1087  BuildTree();
1088 }
1089 
1090 void IsolatedTracksCone::clearTrackVectors() {
1091 
1092  t_v_hnNearTRKs ->clear();
1093  t_v_hnLayers_maxNearP ->clear();
1094  t_v_htrkQual_maxNearP ->clear();
1095  t_v_hmaxNearP_goodTrk ->clear();
1096  t_v_hmaxNearP ->clear();
1097 
1098  t_v_cone_hnNearTRKs ->clear();
1099  t_v_cone_hnLayers_maxNearP ->clear();
1100  t_v_cone_htrkQual_maxNearP ->clear();
1101  t_v_cone_hmaxNearP_goodTrk ->clear();
1102  t_v_cone_hmaxNearP ->clear();
1103 
1104  // t_hScale ->clear();
1105  t_trkNOuterHits ->clear();
1106  t_trkNLayersCrossed ->clear();
1107  t_dtFromLeadJet ->clear();
1108  t_trkP ->clear();
1109  t_trkPt ->clear();
1110  t_trkEta ->clear();
1111  t_trkPhi ->clear();
1112  t_e3x3 ->clear();
1113  t_v_eDR ->clear();
1114  t_v_eMipDR ->clear();
1115  t_h3x3 ->clear();
1116  t_h5x5 ->clear();
1117  t_nRH_h3x3 ->clear();
1118  t_nRH_h5x5 ->clear();
1119 
1120  if (doMC) {
1121  t_simP ->clear();
1122  t_hsim3x3 ->clear();
1123  t_hsim5x5 ->clear();
1124  t_hsim3x3Matched ->clear();
1125  t_hsim5x5Matched ->clear();
1126  t_hsim3x3Rest ->clear();
1127  t_hsim5x5Rest ->clear();
1128  t_hsim3x3Photon ->clear();
1129  t_hsim5x5Photon ->clear();
1130  t_hsim3x3NeutHad ->clear();
1131  t_hsim5x5NeutHad ->clear();
1132  t_hsim3x3CharHad ->clear();
1133  t_hsim5x5CharHad ->clear();
1134  t_hsim3x3PdgMatched ->clear();
1135  t_hsim5x5PdgMatched ->clear();
1136  t_hsim3x3Total ->clear();
1137  t_hsim5x5Total ->clear();
1138 
1139  t_hsim3x3NMatched ->clear();
1140  t_hsim3x3NTotal ->clear();
1141  t_hsim3x3NNeutHad ->clear();
1142  t_hsim3x3NCharHad ->clear();
1143  t_hsim3x3NPhoton ->clear();
1144  t_hsim3x3NRest ->clear();
1145 
1146  t_hsim5x5NMatched ->clear();
1147  t_hsim5x5NTotal ->clear();
1148  t_hsim5x5NNeutHad ->clear();
1149  t_hsim5x5NCharHad ->clear();
1150  t_hsim5x5NPhoton ->clear();
1151  t_hsim5x5NRest ->clear();
1152 
1153  t_trkHcalEne ->clear();
1154  t_trkEcalEne ->clear();
1155  }
1156 
1157  t_distFromHotCell_h3x3 ->clear();
1158  t_ietaFromHotCell_h3x3 ->clear();
1159  t_iphiFromHotCell_h3x3 ->clear();
1160  t_distFromHotCell_h5x5 ->clear();
1161  t_ietaFromHotCell_h5x5 ->clear();
1162  t_iphiFromHotCell_h5x5 ->clear();
1163 
1164  if (doMC) {
1165  t_v_hsimInfoConeMatched ->clear();
1166  t_v_hsimInfoConeRest ->clear();
1167  t_v_hsimInfoConePhoton ->clear();
1168  t_v_hsimInfoConeNeutHad ->clear();
1169  t_v_hsimInfoConeCharHad ->clear();
1170  t_v_hsimInfoConePdgMatched->clear();
1171  t_v_hsimInfoConeTotal ->clear();
1172 
1173  t_v_hsimInfoConeNMatched ->clear();
1174  t_v_hsimInfoConeNRest ->clear();
1175  t_v_hsimInfoConeNPhoton ->clear();
1176  t_v_hsimInfoConeNNeutHad ->clear();
1177  t_v_hsimInfoConeNCharHad ->clear();
1178  t_v_hsimInfoConeNTotal ->clear();
1179 
1180  t_v_hsimCone ->clear();
1181  }
1182 
1183  t_v_hCone ->clear();
1184  t_v_nRecHitsCone->clear();
1185  t_v_nSimHitsCone->clear();
1186 
1187  t_v_distFromHotCell->clear();
1188  t_v_ietaFromHotCell->clear();
1189  t_v_iphiFromHotCell->clear();
1190 
1191  t_v_RH_h3x3_ieta ->clear();
1192  t_v_RH_h3x3_iphi ->clear();
1193  t_v_RH_h3x3_ene ->clear();
1194  t_v_RH_h5x5_ieta ->clear();
1195  t_v_RH_h5x5_iphi ->clear();
1196  t_v_RH_h5x5_ene ->clear();
1197  t_v_RH_r26_ieta ->clear();
1198  t_v_RH_r26_iphi ->clear();
1199  t_v_RH_r26_ene ->clear();
1200  t_v_RH_r44_ieta ->clear();
1201  t_v_RH_r44_iphi ->clear();
1202  t_v_RH_r44_ene ->clear();
1203 
1204  t_v_hlTriggers ->clear();
1205  t_hltHB->clear();
1206  t_hltHE->clear();
1207  t_hltL1Jet15 ->clear();
1208  t_hltJet30 ->clear();
1209  t_hltJet50 ->clear();
1210  t_hltJet80 ->clear();
1211  t_hltJet110 ->clear();
1212  t_hltJet140 ->clear();
1213  t_hltJet180 ->clear();
1214  t_hltL1SingleEG5 ->clear();
1215  t_hltZeroBias ->clear();
1216  t_hltMinBiasHcal ->clear();
1217  t_hltMinBiasEcal ->clear();
1218  t_hltMinBiasPixel ->clear();
1219  t_hltSingleIsoTau30_Trk5 ->clear();
1220  t_hltDoubleLooseIsoTau15_Trk5->clear();
1221 
1222  t_irun->clear();
1223  t_ievt->clear();
1224  t_ilum->clear();
1225 
1226 
1227 
1228 }
1229 
1230 // ----- method called once each job just after ending the event loop ----
1231 void IsolatedTracksCone::endJob() {
1232 
1233  std::cout << "Number of Events Processed " << nEVT << " failed L1 "
1234  << nEVT_failL1 << std::endl;
1235 
1236 }
1237 
1238 
1239 void IsolatedTracksCone::BuildTree(){
1240 
1241 
1242  hRawPt = fs->make<TH1F>("hRawPt ", "hRawPt ", 100, 0.0, 100.0);
1243  hRawP = fs->make<TH1F>("hRawP ", "hRawP ", 100, 0.0, 100.0);
1244  hRawEta = fs->make<TH1F>("hRawEta", "hRawEta", 15, 0.0, 3.0);
1245  hRawPhi = fs->make<TH1F>("hRawPhi", "hRawPhi", 100, -3.2, 3.2);
1246 
1247  ntp = fs->make<TTree>("ntp", "ntp");
1248 
1249 
1250  // Counters
1251  ntp->Branch("nEVT" , &nEVT , "nEVT/I" );
1252  ntp->Branch("leadL1JetPT" , &leadL1JetPT , "leadL1JetPT/D" );
1253  ntp->Branch("leadL1JetEta", &leadL1JetEta, "leadL1JetEta/D");
1254  ntp->Branch("leadL1JetPhi", &leadL1JetPhi, "leadL1JetPhi/D");
1255  ntp->Branch("nTRK", &nTRK, "nTRK/I");
1256  ntp->Branch("nRawTRK" , &nRawTRK ,"nRawTRK/I" );
1257  ntp->Branch("nFailHighPurityQaul", &nFailHighPurityQaul,"nFailHighPurityQaul/I");
1258  ntp->Branch("nFailPt" , &nFailPt ,"nFailPt/I" );
1259  ntp->Branch("nFailEta" , &nFailEta ,"nFailEta/I" );
1260  ntp->Branch("nMissEcal" , &nMissEcal ,"nMissEcal/I" );
1261  ntp->Branch("nMissHcal" , &nMissHcal ,"nMissHcal/I" );
1262 
1263  ntp->Branch("hnNearTRKs" ,"vector<vector<int> > ",&t_v_hnNearTRKs );
1264  ntp->Branch("hnLayers_maxNearP" ,"vector<vector<int> > ",&t_v_hnLayers_maxNearP );
1265  ntp->Branch("htrkQual_maxNearP" ,"vector<vector<int> > ",&t_v_htrkQual_maxNearP );
1266  ntp->Branch("hmaxNearP_goodTrk" ,"vector<vector<double> >",&t_v_hmaxNearP_goodTrk );
1267  ntp->Branch("hmaxNearP" ,"vector<vector<double> >",&t_v_hmaxNearP );
1268 
1269  ntp->Branch("cone_hnNearTRKs" ,"vector<vector<int> > ",&t_v_cone_hnNearTRKs );
1270  ntp->Branch("cone_hnLayers_maxNearP","vector<vector<int> > ",&t_v_cone_hnLayers_maxNearP);
1271  ntp->Branch("cone_htrkQual_maxNearP","vector<vector<int> > ",&t_v_cone_htrkQual_maxNearP);
1272  ntp->Branch("cone_hmaxNearP_goodTrk","vector<vector<double> >",&t_v_cone_hmaxNearP_goodTrk);
1273  ntp->Branch("cone_hmaxNearP" ,"vector<vector<double> >",&t_v_cone_hmaxNearP );
1274 
1275  // ntp->Branch("hScale" , "vector<double>", &t_hScale );
1276  ntp->Branch("trkNOuterHits" , "vector<double>", &t_trkNOuterHits );
1277  ntp->Branch("trkNLayersCrossed", "vector<double>", &t_trkNLayersCrossed);
1278  ntp->Branch("drFromLeadJet" , "vector<double>", &t_dtFromLeadJet );
1279  ntp->Branch("trkP" , "vector<double>", &t_trkP );
1280  ntp->Branch("trkPt" , "vector<double>", &t_trkPt );
1281  ntp->Branch("trkEta" , "vector<double>", &t_trkEta );
1282  ntp->Branch("trkPhi" , "vector<double>", &t_trkPhi );
1283  ntp->Branch("e3x3" , "vector<double>", &t_e3x3 );
1284 
1285  ntp->Branch("e3x3" , "vector<double>" , &t_e3x3 );
1286  ntp->Branch("v_eDR" , "vector<vector<double> >", &t_v_eDR);
1287  ntp->Branch("v_eMipDR" , "vector<vector<double> >", &t_v_eMipDR);
1288 
1289  ntp->Branch("h3x3" , "vector<double>", &t_h3x3 );
1290  ntp->Branch("h5x5" , "vector<double>", &t_h5x5 );
1291  ntp->Branch("nRH_h3x3" , "vector<double>", &t_nRH_h3x3 );
1292  ntp->Branch("nRH_h5x5" , "vector<double>", &t_nRH_h5x5 );
1293 
1294  if (doMC) {
1295  ntp->Branch("simP" , "vector<double>", &t_simP );
1296  ntp->Branch("hsim3x3" , "vector<double>", &t_hsim3x3 );
1297  ntp->Branch("hsim5x5" , "vector<double>", &t_hsim5x5 );
1298 
1299  ntp->Branch("hsim3x3Matched" , "vector<double>", &t_hsim3x3Matched );
1300  ntp->Branch("hsim5x5Matched" , "vector<double>", &t_hsim5x5Matched );
1301  ntp->Branch("hsim3x3Rest" , "vector<double>", &t_hsim3x3Rest );
1302  ntp->Branch("hsim5x5Rest" , "vector<double>", &t_hsim5x5Rest );
1303  ntp->Branch("hsim3x3Photon" , "vector<double>", &t_hsim3x3Photon );
1304  ntp->Branch("hsim5x5Photon" , "vector<double>", &t_hsim5x5Photon );
1305  ntp->Branch("hsim3x3NeutHad" , "vector<double>", &t_hsim3x3NeutHad );
1306  ntp->Branch("hsim5x5NeutHad" , "vector<double>", &t_hsim5x5NeutHad );
1307  ntp->Branch("hsim3x3CharHad" , "vector<double>", &t_hsim3x3CharHad );
1308  ntp->Branch("hsim5x5CharHad" , "vector<double>", &t_hsim5x5CharHad );
1309  ntp->Branch("hsim3x3PdgMatched", "vector<double>", &t_hsim3x3PdgMatched);
1310  ntp->Branch("hsim5x5PdgMatched", "vector<double>", &t_hsim5x5PdgMatched);
1311  ntp->Branch("hsim3x3Total" , "vector<double>", &t_hsim3x3Total );
1312  ntp->Branch("hsim5x5Total" , "vector<double>", &t_hsim5x5Total );
1313 
1314  ntp->Branch("hsim3x3NMatched" , "vector<int>", &t_hsim3x3NMatched );
1315  ntp->Branch("hsim3x3NRest" , "vector<int>", &t_hsim3x3NRest );
1316  ntp->Branch("hsim3x3NPhoton" , "vector<int>", &t_hsim3x3NPhoton );
1317  ntp->Branch("hsim3x3NNeutHad" , "vector<int>", &t_hsim3x3NNeutHad );
1318  ntp->Branch("hsim3x3NCharHad" , "vector<int>", &t_hsim3x3NCharHad );
1319  ntp->Branch("hsim3x3NTotal" , "vector<int>", &t_hsim3x3NTotal );
1320 
1321  ntp->Branch("hsim5x5NMatched" , "vector<int>", &t_hsim5x5NMatched );
1322  ntp->Branch("hsim5x5NRest" , "vector<int>", &t_hsim5x5NRest );
1323  ntp->Branch("hsim5x5NPhoton" , "vector<int>", &t_hsim5x5NPhoton );
1324  ntp->Branch("hsim5x5NNeutHad" , "vector<int>", &t_hsim5x5NNeutHad );
1325  ntp->Branch("hsim5x5NCharHad" , "vector<int>", &t_hsim5x5NCharHad );
1326  ntp->Branch("hsim5x5NTotal" , "vector<int>", &t_hsim5x5NTotal );
1327 
1328  ntp->Branch("trkHcalEne" , "vector<double>", &t_trkHcalEne );
1329  ntp->Branch("trkEcalEne" , "vector<double>", &t_trkEcalEne );
1330  }
1331 
1332  ntp->Branch("distFromHotCell_h3x3", "vector<double>", &t_distFromHotCell_h3x3);
1333  ntp->Branch("ietaFromHotCell_h3x3", "vector<int>", &t_ietaFromHotCell_h3x3);
1334  ntp->Branch("iphiFromHotCell_h3x3", "vector<int>", &t_iphiFromHotCell_h3x3);
1335  ntp->Branch("distFromHotCell_h5x5", "vector<double>", &t_distFromHotCell_h5x5);
1336  ntp->Branch("ietaFromHotCell_h5x5", "vector<int>", &t_ietaFromHotCell_h5x5);
1337  ntp->Branch("iphiFromHotCell_h5x5", "vector<int>", &t_iphiFromHotCell_h5x5);
1338 
1339  if (doMC) {
1340  ntp->Branch("v_hsimInfoConeMatched" ,"vector<vector<double> >",&t_v_hsimInfoConeMatched );
1341  ntp->Branch("v_hsimInfoConeRest" ,"vector<vector<double> >",&t_v_hsimInfoConeRest );
1342  ntp->Branch("v_hsimInfoConePhoton" ,"vector<vector<double> >",&t_v_hsimInfoConePhoton );
1343  ntp->Branch("v_hsimInfoConeNeutHad" ,"vector<vector<double> >",&t_v_hsimInfoConeNeutHad );
1344  ntp->Branch("v_hsimInfoConeCharHad" ,"vector<vector<double> >",&t_v_hsimInfoConeCharHad );
1345  ntp->Branch("v_hsimInfoConePdgMatched","vector<vector<double> >",&t_v_hsimInfoConePdgMatched);
1346  ntp->Branch("v_hsimInfoConeTotal" ,"vector<vector<double> >",&t_v_hsimInfoConeTotal );
1347 
1348  ntp->Branch("v_hsimInfoConeNMatched" ,"vector<vector<int> >" ,&t_v_hsimInfoConeNMatched );
1349  ntp->Branch("v_hsimInfoConeNRest" ,"vector<vector<int> >" ,&t_v_hsimInfoConeNRest );
1350  ntp->Branch("v_hsimInfoConeNPhoton" ,"vector<vector<int> >" ,&t_v_hsimInfoConeNPhoton );
1351  ntp->Branch("v_hsimInfoConeNNeutHad" ,"vector<vector<int> >" ,&t_v_hsimInfoConeNNeutHad );
1352  ntp->Branch("v_hsimInfoConeNCharHad" ,"vector<vector<int> >" ,&t_v_hsimInfoConeNCharHad );
1353  ntp->Branch("v_hsimInfoConeNTotal" ,"vector<vector<int> >" ,&t_v_hsimInfoConeNTotal );
1354 
1355  ntp->Branch("v_hsimCone" ,"vector<vector<double> >",&t_v_hsimCone );
1356  }
1357 
1358  ntp->Branch("v_hCone" ,"vector<vector<double> >",&t_v_hCone );
1359  ntp->Branch("v_nRecHitsCone" ,"vector<vector<int> >" ,&t_v_nRecHitsCone);
1360  ntp->Branch("v_nSimHitsCone" ,"vector<vector<int> >" ,&t_v_nSimHitsCone);
1361 
1362  ntp->Branch("v_distFromHotCell" ,"vector<vector<double> >",&t_v_distFromHotCell );
1363  ntp->Branch("v_ietaFromHotCell" ,"vector<vector<int> >",&t_v_ietaFromHotCell );
1364  ntp->Branch("v_iphiFromHotCell" ,"vector<vector<int> >",&t_v_iphiFromHotCell );
1365 
1366  ntp->Branch("v_RH_h3x3_ieta" ,"vector<vector<int> >" ,&t_v_RH_h3x3_ieta);
1367  ntp->Branch("v_RH_h3x3_iphi" ,"vector<vector<int> >" ,&t_v_RH_h3x3_iphi);
1368  ntp->Branch("v_RH_h3x3_ene" ,"vector<vector<double> >",&t_v_RH_h3x3_ene );
1369  ntp->Branch("v_RH_h5x5_ieta" ,"vector<vector<int> >" ,&t_v_RH_h5x5_ieta);
1370  ntp->Branch("v_RH_h5x5_iphi" ,"vector<vector<int> >" ,&t_v_RH_h5x5_iphi);
1371  ntp->Branch("v_RH_h5x5_ene" ,"vector<vector<double> >",&t_v_RH_h5x5_ene );
1372  ntp->Branch("v_RH_r26_ieta" ,"vector<vector<int> >" ,&t_v_RH_r26_ieta );
1373  ntp->Branch("v_RH_r26_iphi" ,"vector<vector<int> >" ,&t_v_RH_r26_iphi );
1374  ntp->Branch("v_RH_r26_ene" ,"vector<vector<double> >",&t_v_RH_r26_ene );
1375  ntp->Branch("v_RH_r44_ieta" ,"vector<vector<int> >" ,&t_v_RH_r44_ieta );
1376  ntp->Branch("v_RH_r44_iphi" ,"vector<vector<int> >" ,&t_v_RH_r44_iphi );
1377  ntp->Branch("v_RH_r44_ene" ,"vector<vector<double> >",&t_v_RH_r44_ene );
1378 
1379  ntp->Branch("v_hlTriggers" ,"vector<vector<int> >",&t_v_hlTriggers);
1380  ntp->Branch("v_hltHB" ,"vector<int>",&t_hltHB);
1381  ntp->Branch("v_hltHE" ,"vector<int>",&t_hltHE);
1382  ntp->Branch("v_hltL1Jet15" ,"vector<int>",&t_hltL1Jet15 );
1383  ntp->Branch("v_hltJet30" ,"vector<int>",&t_hltJet30 );
1384  ntp->Branch("v_hltJet50" ,"vector<int>",&t_hltJet50 );
1385  ntp->Branch("v_hltJet80" ,"vector<int>",&t_hltJet80 );
1386  ntp->Branch("v_hltJet110" ,"vector<int>",&t_hltJet110 );
1387  ntp->Branch("v_hltJet140" ,"vector<int>",&t_hltJet140 );
1388  ntp->Branch("v_hltJet180" ,"vector<int>",&t_hltJet180 );
1389  ntp->Branch("v_hltL1SingleEG5" ,"vector<int>",&t_hltL1SingleEG5 );
1390  ntp->Branch("v_hltZeroBias" ,"vector<int>",&t_hltZeroBias );
1391  ntp->Branch("v_hltMinBiasHcal" ,"vector<int>",&t_hltMinBiasHcal );
1392  ntp->Branch("v_hltMinBiasEcal" ,"vector<int>",&t_hltMinBiasEcal );
1393  ntp->Branch("v_hltMinBiasPixel" ,"vector<int>",&t_hltMinBiasPixel );
1394  ntp->Branch("v_hltSingleIsoTau30_Trk5" ,"vector<int>",&t_hltSingleIsoTau30_Trk5 );
1395  ntp->Branch("v_hltDoubleLooseIsoTau15_Trk5" ,"vector<int>",&t_hltDoubleLooseIsoTau15_Trk5);
1396 
1397  ntp->Branch("irun" ,"vector<unsigned int>", &t_irun);
1398  ntp->Branch("ievt" ,"vector<unsigned int>", &t_ievt);
1399  ntp->Branch("ilum" ,"vector<unsigned int>", &t_ilum);
1400 
1401 }
1402 
1403 
1404 void IsolatedTracksCone::printTrack(const reco::Track* pTrack) {
1405 
1406  std::string theTrackQuality = "highPurity";
1407  reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality);
1408 
1409  std::cout << " Reference Point " << pTrack->referencePoint() <<"\n"
1410  << " TrackMmentum " << pTrack->momentum()
1411  << " (pt,eta,phi)(" << pTrack->pt()<<","<<pTrack->eta()<<","<<pTrack->phi()<<")"
1412  << " p " << pTrack->p() << "\n"
1413  << " Normalized chi2 " << pTrack->normalizedChi2() <<" charge " << pTrack->charge()
1414  << " qoverp() " << pTrack->qoverp() <<"+-" << pTrack->qoverpError()
1415  << " d0 " << pTrack->d0() << "\n"
1416  << " NValidHits " << pTrack->numberOfValidHits() << " NLostHits " << pTrack->numberOfLostHits()
1417  << " TrackQuality " << pTrack->qualityName(trackQuality_) << " " << pTrack->quality(trackQuality_)
1418  << std::endl;
1419 
1420  if( printTrkHitPattern_ ) {
1421  const reco::HitPattern& p = pTrack->hitPattern();
1422 
1423  for (int i=0; i<p.numberOfHits(); i++) {
1424  p.printHitPattern(i, std::cout);
1425  }
1426  }
1427 
1428 }
1429 
1430 double IsolatedTracksCone::DeltaPhi(double v1, double v2) {
1431  // Computes the correctly normalized phi difference
1432  // v1, v2 = phi of object 1 and 2
1433 
1434  double pi = 3.141592654;
1435  double twopi = 6.283185307;
1436 
1437  double diff = std::abs(v2 - v1);
1438  double corr = twopi - diff;
1439  if (diff < pi){ return diff;} else { return corr;}
1440 }
1441 
1442 
1443 double IsolatedTracksCone::DeltaR(double eta1, double phi1,
1444  double eta2, double phi2) {
1445  double deta = eta1 - eta2;
1446  double dphi = DeltaPhi(phi1, phi2);
1447  return std::sqrt(deta*deta + dphi*dphi);
1448 }
1449 
1450 
1451 //define this as a plug-in
1452 DEFINE_FWK_MODULE(IsolatedTracksCone);
1453 
1454 
1455 
1456 
1457 
1458 
1459 
1460 
1461 
1462 
1463 
1464 
1465 
1466 
1467 
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:42
IsolatedTracksCone::t_hsim5x5NTotal
std::vector< int > * t_hsim5x5NTotal
Definition: IsolatedTracksCone.h:229
IsolatedTracksCone::t_v_distFromHotCell
std::vector< std::vector< double > > * t_v_distFromHotCell
Definition: IsolatedTracksCone.h:261
reco::TrackBase::qoverp
double qoverp() const
q/p
Definition: TrackBase.h:115
IsolatedTracksCone::t_hsim5x5CharHad
std::vector< double > * t_hsim5x5CharHad
Definition: IsolatedTracksCone.h:211
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:129
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
CaloGeometry::getSubdetectorGeometry
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:43
edm::EDAnalyzer::beginJob
virtual void beginJob()
Definition: EDAnalyzer.h:58
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:44
IsolatedTracksCone::t_v_hmaxNearP_goodTrk
std::vector< std::vector< double > > * t_v_hmaxNearP_goodTrk
Definition: IsolatedTracksCone.h:160
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
IsolatedTracksCone::t_v_hlTriggers
std::vector< std::vector< int > > * t_v_hlTriggers
Definition: IsolatedTracksCone.h:266
i
int i
Definition: DBlmapReader.cc:9
Parameters.parameters
dictionary parameters
Definition: Parameters.py:2
IsolatedTracksCone::t_trkEcalEne
std::vector< double > * t_trkEcalEne
Definition: IsolatedTracksCone.h:239
edm::Event::triggerNames
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
Definition: Event.cc:199
IsolatedTracksCone::t_v_hsimCone
std::vector< std::vector< double > > * t_v_hsimCone
Definition: IsolatedTracksCone.h:256
IsolatedTracksCone::parameters_
TrackAssociatorParameters parameters_
Definition: IsolatedTracksCone.h:132
reco::TrackBase::qualityName
static std::string qualityName(TrackQuality)
Definition: TrackBase.h:406
IsolatedTracksCone::t_hltJet30
std::vector< int > * t_hltJet30
Definition: IsolatedTracksCone.h:270
IsolatedTracksCone::t_trkPt
std::vector< double > * t_trkPt
Definition: IsolatedTracksCone.h:175
IsolatedTracksCone::t_ietaFromHotCell_h5x5
std::vector< int > * t_ietaFromHotCell_h5x5
Definition: IsolatedTracksCone.h:235
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)
IsolatedTracksCone::t_nRH_h3x3
std::vector< double > * t_nRH_h3x3
Definition: IsolatedTracksCone.h:188
IsolatedTracksCone::t_distFromHotCell_h3x3
std::vector< double > * t_distFromHotCell_h3x3
Definition: IsolatedTracksCone.h:231
reco::TrackBase::momentum
const Vector & momentum() const
track momentum vector
Definition: TrackBase.h:150
IsolatedTracksCone::t_v_hsimInfoConeNCharHad
std::vector< std::vector< int > > * t_v_hsimInfoConeNCharHad
Definition: IsolatedTracksCone.h:253
IsolatedTracksCone::t_hsim5x5NNeutHad
std::vector< int > * t_hsim5x5NNeutHad
Definition: IsolatedTracksCone.h:227
IsolatedTracksCone::t_hsim5x5Photon
std::vector< double > * t_hsim5x5Photon
Definition: IsolatedTracksCone.h:207
reco::TrackBase::d0
double d0() const
dxy parameter in perigee convention (d0 = - dxy)
Definition: TrackBase.h:123
spr::getDistInPlaneTrackDir
double getDistInPlaneTrackDir(const GlobalPoint &caloPoint, const GlobalVector &caloVector, const GlobalPoint &rechitPoint)
Definition: FindDistCone.cc:8
IsolatedTracksCone::t_hltMinBiasHcal
std::vector< int > * t_hltMinBiasHcal
Definition: IsolatedTracksCone.h:278
IsolatedTracksCone::t_v_cone_hmaxNearP_goodTrk
std::vector< std::vector< double > > * t_v_cone_hmaxNearP_goodTrk
Definition: IsolatedTracksCone.h:166
IsolatedTracksCone::t_hsim3x3CharHad
std::vector< double > * t_hsim3x3CharHad
Definition: IsolatedTracksCone.h:210
IsolatedTracksCone::t_hsim5x5NPhoton
std::vector< int > * t_hsim5x5NPhoton
Definition: IsolatedTracksCone.h:226
IsolatedTracksCone::t_v_RH_r44_ene
std::vector< std::vector< double > > * t_v_RH_r44_ene
Definition: IsolatedTracksCone.h:296
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)
IsolatedTracksCone::t_iphiFromHotCell_h3x3
std::vector< int > * t_iphiFromHotCell_h3x3
Definition: IsolatedTracksCone.h:233
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:111
IsolatedTracksCone::t_v_cone_hnLayers_maxNearP
std::vector< std::vector< int > > * t_v_cone_hnLayers_maxNearP
Definition: IsolatedTracksCone.h:164
IsolatedTracksCone::t_hsim5x5PdgMatched
std::vector< double > * t_hsim5x5PdgMatched
Definition: IsolatedTracksCone.h:213
IsolatedTracksCone::t_hsim3x3NeutHad
std::vector< double > * t_hsim3x3NeutHad
Definition: IsolatedTracksCone.h:208
IsolatedTracksCone::t_nRH_h5x5
std::vector< double > * t_nRH_h5x5
Definition: IsolatedTracksCone.h:189
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:95
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
TrackDetectorAssociator::getFreeTrajectoryState
static FreeTrajectoryState getFreeTrajectoryState(const edm::EventSetup &, const reco::Track &)
get FreeTrajectoryState from different track representations
Definition: TrackDetectorAssociator.cc:615
IsolatedTracksCone::t_v_hnLayers_maxNearP
std::vector< std::vector< int > > * t_v_hnLayers_maxNearP
Definition: IsolatedTracksCone.h:158
IsolatedTracksCone::t_hltHB
std::vector< int > * t_hltHB
Definition: IsolatedTracksCone.h:267
IsolatedTracksCone::t_v_hmaxNearP
std::vector< std::vector< double > > * t_v_hmaxNearP
Definition: IsolatedTracksCone.h:161
IsolatedTracksCone::t_v_RH_h3x3_ene
std::vector< std::vector< double > > * t_v_RH_h3x3_ene
Definition: IsolatedTracksCone.h:287
TrackDetectorAssociator::useDefaultPropagator
void useDefaultPropagator()
use the default propagator
Definition: TrackDetectorAssociator.cc:120
IsolatedTracksCone::t_hsim5x5NCharHad
std::vector< int > * t_hsim5x5NCharHad
Definition: IsolatedTracksCone.h:228
abs
#define abs(x)
Definition: mlp_lapack.h:159
IsolatedTracksCone::t_v_hsimInfoConeNNeutHad
std::vector< std::vector< int > > * t_v_hsimInfoConeNNeutHad
Definition: IsolatedTracksCone.h:252
IsolatedTracksCone::t_v_nRecHitsCone
std::vector< std::vector< int > > * t_v_nRecHitsCone
Definition: IsolatedTracksCone.h:258
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:139
IsolatedTracksCone::t_v_hsimInfoConeNMatched
std::vector< std::vector< int > > * t_v_hsimInfoConeNMatched
Definition: IsolatedTracksCone.h:249
reco::TrackBase::numberOfLostHits
unsigned short numberOfLostHits() const
number of cases where track crossed a layer without getting a hit.
Definition: TrackBase.h:234
IsolatedTracksCone::t_trkP
std::vector< double > * t_trkP
Definition: IsolatedTracksCone.h:173
IsolatedTracksCone::t_hltL1SingleEG5
std::vector< int > * t_hltL1SingleEG5
Definition: IsolatedTracksCone.h:276
IsolatedTracksCone::t_hsim5x5NRest
std::vector< int > * t_hsim5x5NRest
Definition: IsolatedTracksCone.h:225
IsolatedTracksCone::t_hsim5x5NeutHad
std::vector< double > * t_hsim5x5NeutHad
Definition: IsolatedTracksCone.h:209
IsolatedTracksCone::t_hltJet50
std::vector< int > * t_hltJet50
Definition: IsolatedTracksCone.h:271
IsolatedTracksCone::t_v_cone_hmaxNearP
std::vector< std::vector< double > > * t_v_cone_hmaxNearP
Definition: IsolatedTracksCone.h:167
IsolatedTracksCone::t_simP
std::vector< double > * t_simP
Definition: IsolatedTracksCone.h:174
reco::TrackBase::referencePoint
const Point & referencePoint() const
Reference point on the track.
Definition: TrackBase.h:153
IsolatedTracksCone::t_hsim3x3NCharHad
std::vector< int > * t_hsim3x3NCharHad
Definition: IsolatedTracksCone.h:221
IsolatedTracksCone::t_v_iphiFromHotCell
std::vector< std::vector< int > > * t_v_iphiFromHotCell
Definition: IsolatedTracksCone.h:263
IsolatedTracksCone::t_v_hnNearTRKs
std::vector< std::vector< int > > * t_v_hnNearTRKs
Definition: IsolatedTracksCone.h:157
IsolatedTracksCone::t_v_cone_htrkQual_maxNearP
std::vector< std::vector< int > > * t_v_cone_htrkQual_maxNearP
Definition: IsolatedTracksCone.h:165
TrackDetMatchInfo::trkGlobPosAtHcal
math::XYZPoint trkGlobPosAtHcal
Definition: TrackDetMatchInfo.h:40
IsolatedTracksCone::t_v_RH_r26_ene
std::vector< std::vector< double > > * t_v_RH_r26_ene
Definition: IsolatedTracksCone.h:293
edm::LuminosityBlockBase::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: LuminosityBlockBase.h:45
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:45
IsolatedTracksCone::t_hltDoubleLooseIsoTau15_Trk5
std::vector< int > * t_hltDoubleLooseIsoTau15_Trk5
Definition: IsolatedTracksCone.h:282
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 debug=false)
IsolatedTracksCone::t_hltJet140
std::vector< int > * t_hltJet140
Definition: IsolatedTracksCone.h:274
iEvent
int iEvent
Definition: GenABIO.cc:243
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:141
IsolatedTracksCone::t_hltMinBiasPixel
std::vector< int > * t_hltMinBiasPixel
Definition: IsolatedTracksCone.h:280
IsolatedTracksCone::t_trkPhi
std::vector< double > * t_trkPhi
Definition: IsolatedTracksCone.h:177
spr::eCone_hcal
double eCone_hcal(const CaloGeometry *geo, edm::Handle< T > &hits, const GlobalPoint &hpoint1, const GlobalPoint &point1, double dR, const GlobalVector &trackMom, int &nRecHits, double hbThr=-100, double heThr=-100, double hfThr=-100, double hoThr=-100, double tMin=-500, double tMax=500)
IsolatedTracksCone::t_hsim3x3Total
std::vector< double > * t_hsim3x3Total
Definition: IsolatedTracksCone.h:214
IsolatedTracksCone::t_hltHE
std::vector< int > * t_hltHE
Definition: IsolatedTracksCone.h:268
reco::HitPattern::trackerLayersWithMeasurement
int trackerLayersWithMeasurement() const
Definition: HitPattern.h:705
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:46
IsolatedTracksCone::t_v_RH_h3x3_iphi
std::vector< std::vector< int > > * t_v_RH_h3x3_iphi
Definition: IsolatedTracksCone.h:286
IsolatedTracksCone::t_iphiFromHotCell_h5x5
std::vector< int > * t_iphiFromHotCell_h5x5
Definition: IsolatedTracksCone.h:236
IsolatedTracksCone::t_v_hsimInfoConePdgMatched
std::vector< std::vector< double > > * t_v_hsimInfoConePdgMatched
Definition: IsolatedTracksCone.h:246
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:131
IsolatedTracksCone::t_v_htrkQual_maxNearP
std::vector< std::vector< int > > * t_v_htrkQual_maxNearP
Definition: IsolatedTracksCone.h:159
IsolatedTracksCone::t_hltMinBiasEcal
std::vector< int > * t_hltMinBiasEcal
Definition: IsolatedTracksCone.h:279
IsolatedTracksCone::t_v_RH_r44_iphi
std::vector< std::vector< int > > * t_v_RH_r44_iphi
Definition: IsolatedTracksCone.h:295
IsolatedTracksCone::t_hsim5x5Matched
std::vector< double > * t_hsim5x5Matched
Definition: IsolatedTracksCone.h:203
IsolatedTracksCone::t_v_hsimInfoConeTotal
std::vector< std::vector< double > > * t_v_hsimInfoConeTotal
Definition: IsolatedTracksCone.h:247
reco::HitPattern::numberOfHits
int numberOfHits() const
Definition: HitPattern.cc:213
IsolatedTracksCone::t_v_RH_h5x5_ene
std::vector< std::vector< double > > * t_v_RH_h5x5_ene
Definition: IsolatedTracksCone.h:290
IsolatedTracksCone::DeltaPhi
double DeltaPhi(double v1, double v2)
Definition: IsolatedTracksCone.cc:1430
IsolatedTracksCone::t_v_hsimInfoConeNRest
std::vector< std::vector< int > > * t_v_hsimInfoConeNRest
Definition: IsolatedTracksCone.h:250
edm::Event::getLuminosityBlock
LuminosityBlock const & getLuminosityBlock() const
Definition: Event.h:59
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:10
reco::TrackBase::numberOfValidHits
unsigned short numberOfValidHits() const
number of valid hits found
Definition: TrackBase.h:232
IsolatedTracksCone::t_v_hsimInfoConeCharHad
std::vector< std::vector< double > > * t_v_hsimInfoConeCharHad
Definition: IsolatedTracksCone.h:245
IsolatedTracksCone::_L1extraFwdJetSource
edm::InputTag _L1extraFwdJetSource
Definition: IsolatedTracksCone.h:123
reco::TrackBase::hitPattern
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:223
first
bool first
Definition: L1TdeRCT.cc:94
IsolatedTracksCone::t_hsim5x5Total
std::vector< double > * t_hsim5x5Total
Definition: IsolatedTracksCone.h:215
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
IsolatedTracksCone::t_hsim3x3Photon
std::vector< double > * t_hsim3x3Photon
Definition: IsolatedTracksCone.h:206
IsolatedTracksCone::t_hsim3x3PdgMatched
std::vector< double > * t_hsim3x3PdgMatched
Definition: IsolatedTracksCone.h:212
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
IsolatedTracksCone::t_v_ietaFromHotCell
std::vector< std::vector< int > > * t_v_ietaFromHotCell
Definition: IsolatedTracksCone.h:262
IsolatedTracksCone::t_hsim5x5Rest
std::vector< double > * t_hsim5x5Rest
Definition: IsolatedTracksCone.h:205
IsolatedTracksCone::t_hltJet110
std::vector< int > * t_hltJet110
Definition: IsolatedTracksCone.h:273
corr
JetCorrectorParameters corr
Definition: classes.h:9
IsolatedTracksCone::t_v_eMipDR
std::vector< std::vector< double > > * t_v_eMipDR
Definition: IsolatedTracksCone.h:181
IsolatedTracksCone::t_hltL1Jet15
std::vector< int > * t_hltL1Jet15
Definition: IsolatedTracksCone.h:269
reco::TrackBase::qoverpError
double qoverpError() const
error on signed transverse curvature
Definition: TrackBase.h:192
CaloSubdetectorGeometry::getClosestCell
virtual DetId getClosestCell(const GlobalPoint &r) const
Definition: CaloSubdetectorGeometry.cc:55
IsolatedTracksCone::t_hsim3x3NRest
std::vector< int > * t_hsim3x3NRest
Definition: IsolatedTracksCone.h:218
IsolatedTracksCone::t_v_hsimInfoConeRest
std::vector< std::vector< double > > * t_v_hsimInfoConeRest
Definition: IsolatedTracksCone.h:242
IsolatedTracksCone::t_v_hsimInfoConeNPhoton
std::vector< std::vector< int > > * t_v_hsimInfoConeNPhoton
Definition: IsolatedTracksCone.h:251
IsolatedTracksCone::t_v_RH_r26_ieta
std::vector< std::vector< int > > * t_v_RH_r26_ieta
Definition: IsolatedTracksCone.h:291
IsolatedTracksCone::t_dtFromLeadJet
std::vector< double > * t_dtFromLeadJet
Definition: IsolatedTracksCone.h:172
IsolatedTracksCone::t_v_RH_h5x5_ieta
std::vector< std::vector< int > > * t_v_RH_h5x5_ieta
Definition: IsolatedTracksCone.h:288
IsolatedTracksCone::t_v_RH_h3x3_ieta
std::vector< std::vector< int > > * t_v_RH_h3x3_ieta
Definition: IsolatedTracksCone.h:285
DetId
Definition: DetId.h:20
IsolatedTracksCone::t_ievt
std::vector< unsigned int > * t_ievt
Definition: IsolatedTracksCone.h:299
IsolatedTracksCone::t_trkNLayersCrossed
std::vector< double > * t_trkNLayersCrossed
Definition: IsolatedTracksCone.h:171
IsolatedTracksCone::_L1extraCenJetSource
edm::InputTag _L1extraCenJetSource
Definition: IsolatedTracksCone.h:123
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:46
IsolatedTracksCone::t_v_hsimInfoConeNeutHad
std::vector< std::vector< double > > * t_v_hsimInfoConeNeutHad
Definition: IsolatedTracksCone.h:244
IsolatedTracksCone::t_v_nSimHitsCone
std::vector< std::vector< int > > * t_v_nSimHitsCone
Definition: IsolatedTracksCone.h:259
IsolatedTracksCone::t_hltSingleIsoTau30_Trk5
std::vector< int > * t_hltSingleIsoTau30_Trk5
Definition: IsolatedTracksCone.h:281
IsolatedTracksCone::t_hsim5x5
std::vector< double > * t_hsim5x5
Definition: IsolatedTracksCone.h:186
TrackDetMatchInfo::trkMomAtEcal
GlobalVector trkMomAtEcal
Definition: TrackDetMatchInfo.h:43
edm::TriggerNames::triggerName
std::string const & triggerName(unsigned int index) const
Definition: TriggerNames.cc:27
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
IsolatedTracksCone::t_v_RH_h5x5_iphi
std::vector< std::vector< int > > * t_v_RH_h5x5_iphi
Definition: IsolatedTracksCone.h:289
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
IsolatedTracksCone::IsolatedTracksCone
IsolatedTracksCone(const edm::ParameterSet &)
Definition: IsolatedTracksCone.cc:37
IsolatedTracksCone::t_v_hsimInfoConeNTotal
std::vector< std::vector< int > > * t_v_hsimInfoConeNTotal
Definition: IsolatedTracksCone.h:254
IsolatedTracksCone::t_hsim3x3Matched
std::vector< double > * t_hsim3x3Matched
Definition: IsolatedTracksCone.h:202
edm::Handle::product
T const * product() const
Definition: Handle.h:74
reco::TrackBase::quality
bool quality(const TrackQuality) const
Track quality.
Definition: TrackBase.h:377
IsolatedTracksCone::t_hsim3x3NNeutHad
std::vector< int > * t_hsim3x3NNeutHad
Definition: IsolatedTracksCone.h:220
IsolatedTracksCone::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: IsolatedTracksCone.cc:97
IsolatedTracksCone::fs
edm::Service< TFileService > fs
Definition: IsolatedTracksCone.h:303
TrackDetMatchInfo::trkMomAtHcal
GlobalVector trkMomAtHcal
Definition: TrackDetMatchInfo.h:44
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
IsolatedTracksCone::t_hltZeroBias
std::vector< int > * t_hltZeroBias
Definition: IsolatedTracksCone.h:277
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
IsolatedTracksCone::t_h5x5
std::vector< double > * t_h5x5
Definition: IsolatedTracksCone.h:184
IsolatedTracksCone::t_trkNOuterHits
std::vector< double > * t_trkNOuterHits
Definition: IsolatedTracksCone.h:170
IsolatedTracksCone::t_v_hCone
std::vector< std::vector< double > > * t_v_hCone
Definition: IsolatedTracksCone.h:257
p1
double p1[4]
Definition: TauolaWrapper.h:89
spr::eCone_ecal
double eCone_ecal(const CaloGeometry *geo, edm::Handle< T > &barrelhits, edm::Handle< T > &endcaphits, const GlobalPoint &hpoint1, const GlobalPoint &point1, double dR, const GlobalVector &trackMom, int &nRecHits, double ebThr=-100, double eeThr=-100, double tMin=-500, double tMax=500)
IsolatedTracksCone::t_hsim3x3NPhoton
std::vector< int > * t_hsim3x3NPhoton
Definition: IsolatedTracksCone.h:219
IsolatedTracksCone::t_hsim3x3Rest
std::vector< double > * t_hsim3x3Rest
Definition: IsolatedTracksCone.h:204
IsolatedTracksCone::t_hsim3x3
std::vector< double > * t_hsim3x3
Definition: IsolatedTracksCone.h:185
IsolatedTracksCone::t_hsim3x3NMatched
std::vector< int > * t_hsim3x3NMatched
Definition: IsolatedTracksCone.h:217
IsolatedTracksCone::t_v_RH_r44_ieta
std::vector< std::vector< int > > * t_v_RH_r44_ieta
Definition: IsolatedTracksCone.h:294
IsolatedTracksCone::t_trkHcalEne
std::vector< double > * t_trkHcalEne
Definition: IsolatedTracksCone.h:238
IsolatedTracksCone::t_irun
std::vector< unsigned int > * t_irun
Definition: IsolatedTracksCone.h:298
IsolatedTracksCone::t_hltJet180
std::vector< int > * t_hltJet180
Definition: IsolatedTracksCone.h:275
IsolatedTracksCone::t_hsim3x3NTotal
std::vector< int > * t_hsim3x3NTotal
Definition: IsolatedTracksCone.h:222
TrackDetectorAssociator::associate
TrackDetMatchInfo associate(const edm::Event &, const edm::EventSetup &, const FreeTrajectoryState &, const AssociatorParameters &)
Definition: TrackDetectorAssociator.cc:159
TFileDirectory::make
T * make() const
make new ROOT object
Definition: TFileDirectory.h:55
IsolatedTracksCone::printTrack
void printTrack(const reco::Track *pTrack)
Definition: IsolatedTracksCone.cc:1404
patRefSel_triggerSelection_cff.triggerResults
tuple triggerResults
Definition: patRefSel_triggerSelection_cff.py:10
IsolatedTracksCone::t_ietaFromHotCell_h3x3
std::vector< int > * t_ietaFromHotCell_h3x3
Definition: IsolatedTracksCone.h:232
TrackDetMatchInfo::trkGlobPosAtEcal
math::XYZPoint trkGlobPosAtEcal
Track position at different parts of the calorimeter.
Definition: TrackDetMatchInfo.h:39
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
reco::TrackBase::charge
int charge() const
track electric charge
Definition: TrackBase.h:113
IsolatedTracksCone::t_distFromHotCell_h5x5
std::vector< double > * t_distFromHotCell_h5x5
Definition: IsolatedTracksCone.h:234
IsolatedTracksCone::_L1extraTauJetSource
edm::InputTag _L1extraTauJetSource
Definition: IsolatedTracksCone.h:123
IsolatedTracksCone::t_e3x3
std::vector< double > * t_e3x3
Definition: IsolatedTracksCone.h:178
reco::HitPattern::printHitPattern
void printHitPattern(int position, std::ostream &stream) const
Definition: HitPattern.cc:556
IsolatedTracksCone::DeltaR
double DeltaR(double eta1, double phi1, double eta2, double phi2)
Definition: IsolatedTracksCone.cc:1443
pi
double pi
Definition: Vx3DHLTAnalyzer.h:63
IsolatedTracksCone::t_trkEta
std::vector< double > * t_trkEta
Definition: IsolatedTracksCone.h:176
IsolatedTracksCone::t_h3x3
std::vector< double > * t_h3x3
Definition: IsolatedTracksCone.h:183
spr::coneChargeIsolation
double coneChargeIsolation(const edm::Event &iEvent, const edm::EventSetup &iSetup, reco::TrackCollection::const_iterator trkItr, edm::Handle< reco::TrackCollection > trkCollection, TrackDetectorAssociator &associator, TrackAssociatorParameters &parameters_, std::string theTrackQuality, int &nNearTRKs, int &nLayers_maxNearP, int &trkQual_maxNearP, double &maxNearP_goodTrk, const GlobalPoint &hpoint1, const GlobalVector &trackMom, double dR)
Definition: ChargeIsolation.cc:177
IsolatedTracksCone::t_v_hsimInfoConePhoton
std::vector< std::vector< double > > * t_v_hsimInfoConePhoton
Definition: IsolatedTracksCone.h:243
IsolatedTracksCone::t_ilum
std::vector< unsigned int > * t_ilum
Definition: IsolatedTracksCone.h:300
IsolatedTracksCone::trackAssociator_
TrackDetectorAssociator * trackAssociator_
Definition: IsolatedTracksCone.h:133
IsolatedTracksCone::t_hltJet80
std::vector< int > * t_hltJet80
Definition: IsolatedTracksCone.h:272
IsolatedTracksCone::t_v_eDR
std::vector< std::vector< double > > * t_v_eDR
Definition: IsolatedTracksCone.h:180
TrackAssociatorParameters::loadParameters
void loadParameters(const edm::ParameterSet &)
Definition: TrackAssociatorParameters.cc:20
IsolatedTracksCone::t_v_hsimInfoConeMatched
std::vector< std::vector< double > > * t_v_hsimInfoConeMatched
Definition: IsolatedTracksCone.h:241
IsolatedTracksCone::t_v_RH_r26_iphi
std::vector< std::vector< int > > * t_v_RH_r26_iphi
Definition: IsolatedTracksCone.h:292
IsolatedTracksCone::t_hsim5x5NMatched
std::vector< int > * t_hsim5x5NMatched
Definition: IsolatedTracksCone.h:224
IsolatedTracksCone::t_v_cone_hnNearTRKs
std::vector< std::vector< int > > * t_v_cone_hnNearTRKs
Definition: IsolatedTracksCone.h:163
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)
Generated for CMSSW Reference Manual by doxygen 1.8.5