CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
IsolatedTracksHcalScale.cc
Go to the documentation of this file.
1 
2 #include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
3 #include "DataFormats/TrackReco/interface/TrackBase.h"
4 #include "Calibration/IsolatedParticles/plugins/IsolatedTracksHcalScale.h"
5 
6 #include "Calibration/IsolatedParticles/interface/CaloPropagateTrack.h"
7 #include "Calibration/IsolatedParticles/interface/ChargeIsolation.h"
8 #include "Calibration/IsolatedParticles/interface/eCone.h"
9 #include "Calibration/IsolatedParticles/interface/eHCALMatrixExtra.h"
10 #include "MagneticField/Engine/interface/MagneticField.h"
11 #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
12 #include "CondFormats/DataRecord/interface/EcalChannelStatusRcd.h"
13 #include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"
14 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
15 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"
16 
17 IsolatedTracksHcalScale::IsolatedTracksHcalScale(const edm::ParameterSet& iConfig) {
18 
19  //now do what ever initialization is needed
20  doMC = iConfig.getUntrackedParameter<bool>("DoMC", false);
21  myverbose = iConfig.getUntrackedParameter<int>("Verbosity", 5 );
22  theTrackQuality = iConfig.getUntrackedParameter<std::string>("TrackQuality","highPurity");
23  reco::TrackBase::TrackQuality trackQuality_=reco::TrackBase::qualityByName(theTrackQuality);
24  selectionParameters.minPt = iConfig.getUntrackedParameter<double>("MinTrackPt", 10.0);
25  selectionParameters.minQuality = trackQuality_;
26  selectionParameters.maxDxyPV = iConfig.getUntrackedParameter<double>("MaxDxyPV", 0.2);
27  selectionParameters.maxDzPV = iConfig.getUntrackedParameter<double>("MaxDzPV", 5.0);
28  selectionParameters.maxChi2 = iConfig.getUntrackedParameter<double>("MaxChi2", 5.0);
29  selectionParameters.maxDpOverP = iConfig.getUntrackedParameter<double>("MaxDpOverP", 0.1);
30  selectionParameters.minOuterHit = iConfig.getUntrackedParameter<int>("MinOuterHit", 4);
31  selectionParameters.minLayerCrossed = iConfig.getUntrackedParameter<int>("MinLayerCrossed", 8);
32  selectionParameters.maxInMiss = iConfig.getUntrackedParameter<int>("MaxInMiss", 0);
33  selectionParameters.maxOutMiss = iConfig.getUntrackedParameter<int>("MaxOutMiss", 0);
34  a_coneR = iConfig.getUntrackedParameter<double>("ConeRadius",34.98);
35  a_charIsoR = a_coneR + 28.9;
36  a_neutIsoR = a_charIsoR*0.726;
37  a_mipR = iConfig.getUntrackedParameter<double>("ConeRadiusMIP",14.0);
38  tMinE_ = iConfig.getUntrackedParameter<double>("TimeMinCutECAL", -500.);
39  tMaxE_ = iConfig.getUntrackedParameter<double>("TimeMaxCutECAL", 500.);
40 
41  if (myverbose>=0) {
42  std::cout <<"Parameters read from config file \n"
43  <<" doMC " << doMC
44  <<"\t myverbose " << myverbose
45  <<"\t minPt " << selectionParameters.minPt
46  <<"\t theTrackQuality " << theTrackQuality
47  <<"\t minQuality " << selectionParameters.minQuality
48  <<"\t maxDxyPV " << selectionParameters.maxDxyPV
49  <<"\t maxDzPV " << selectionParameters.maxDzPV
50  <<"\t maxChi2 " << selectionParameters.maxChi2
51  <<"\t maxDpOverP " << selectionParameters.maxDpOverP
52  <<"\t minOuterHit " << selectionParameters.minOuterHit
53  <<"\t minLayerCrossed " << selectionParameters.minLayerCrossed
54  <<"\t maxInMiss " << selectionParameters.maxInMiss
55  <<"\t maxOutMiss " << selectionParameters.maxOutMiss
56  <<"\t a_coneR " << a_coneR
57  <<"\t a_charIsoR " << a_charIsoR
58  <<"\t a_neutIsoR " << a_neutIsoR
59  <<"\t a_mipR " << a_mipR
60  <<"\t time Range (" << tMinE_ << ":" << tMaxE_ << ")"
61  << std::endl;
62  }
63  initL1 = false;
64 
65 }
66 
67 IsolatedTracksHcalScale::~IsolatedTracksHcalScale() {}
68 
69 void IsolatedTracksHcalScale::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
70 
71  edm::ESHandle<MagneticField> bFieldH;
72  iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
73  bField = bFieldH.product();
74 
75  // get handles to calogeometry and calotopology
76  edm::ESHandle<CaloGeometry> pG;
77  iSetup.get<CaloGeometryRecord>().get(pG);
78  const CaloGeometry* geo = pG.product();
79 
80  edm::ESHandle<CaloTopology> theCaloTopology;
81  iSetup.get<CaloTopologyRecord>().get(theCaloTopology);
82  const CaloTopology *caloTopology = theCaloTopology.product();
83 
84  /*
85  edm::ESHandle<HcalTopology> htopo;
86  iSetup.get<IdealGeometryRecord>().get(htopo);
87  const HcalTopology* theHBHETopology = htopo.product();
88  */
89 
90  // Retrieve the good/bad ECAL channels from the DB
91  edm::ESHandle<EcalChannelStatus> ecalChStatus;
92  iSetup.get<EcalChannelStatusRcd>().get(ecalChStatus);
93  const EcalChannelStatus* theEcalChStatus = ecalChStatus.product();
94 
95  /*
96  // Retrieve trigger tower map
97  edm::ESHandle<EcalTrigTowerConstituentsMap> hTtmap;
98  iSetup.get<IdealGeometryRecord>().get(hTtmap);
99  const EcalTrigTowerConstituentsMap& ttMap = *hTtmap;
100  */
101 
102  clearTreeVectors();
103 
104  nEventProc++;
105 
106  t_RunNo = iEvent.id().run();
107  t_EvtNo = iEvent.id().event();
108  t_Lumi = iEvent.luminosityBlock();
109  t_Bunch = iEvent.bunchCrossing();
110  if (myverbose>0) std::cout << nEventProc << " Run " << t_RunNo << " Event " << t_EvtNo << " Lumi " << t_Lumi << " Bunch " << t_Bunch << std::endl;
111 
112  edm::Handle<reco::TrackCollection> trkCollection;
113  iEvent.getByLabel("generalTracks", trkCollection);
114 
115  edm::Handle<reco::VertexCollection> recVtxs;
116  iEvent.getByLabel("offlinePrimaryVertices",recVtxs);
117 
118  // Get the beamspot
119  edm::Handle<reco::BeamSpot> beamSpotH;
120  iEvent.getByLabel("offlineBeamSpot", beamSpotH);
121 
122  math::XYZPoint leadPV(0,0,0);
123  if (recVtxs->size()>0 && !((*recVtxs)[0].isFake())) {
124  leadPV = math::XYZPoint( (*recVtxs)[0].x(),(*recVtxs)[0].y(), (*recVtxs)[0].z() );
125  } else if (beamSpotH.isValid()) {
126  leadPV = beamSpotH->position();
127  }
128 
129  if (myverbose>0) {
130  std::cout << "Primary Vertex " << leadPV;
131  if (beamSpotH.isValid()) std::cout << " Beam Spot " << beamSpotH->position();
132  std::cout << std::endl;
133  }
134 
135  std::vector<spr::propagatedTrackDirection> trkCaloDirections;
136  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality, trkCaloDirections, (myverbose>2));
137  std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
138 
139  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
140  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
141  iEvent.getByLabel("ecalRecHit","EcalRecHitsEB",barrelRecHitsHandle);
142  iEvent.getByLabel("ecalRecHit","EcalRecHitsEE",endcapRecHitsHandle);
143 
144  edm::Handle<HBHERecHitCollection> hbhe;
145  iEvent.getByLabel("hbhereco",hbhe);
146  const HBHERecHitCollection Hithbhe = *(hbhe.product());
147 
148  //get Handles to SimTracks and SimHits
149  edm::Handle<edm::SimTrackContainer> SimTk;
150  edm::SimTrackContainer::const_iterator simTrkItr;
151  edm::Handle<edm::SimVertexContainer> SimVtx;
152  edm::SimVertexContainer::const_iterator vtxItr = SimVtx->begin();
153 
154  //get Handles to PCaloHitContainers of eb/ee/hbhe
155  edm::Handle<edm::PCaloHitContainer> pcaloeb;
156  edm::Handle<edm::PCaloHitContainer> pcaloee;
157  edm::Handle<edm::PCaloHitContainer> pcalohh;
158 
159  //associates tracker rechits/simhits to a track
160  TrackerHitAssociator* associate=0;
161 
162  if (doMC) {
163  iEvent.getByLabel("g4SimHits",SimTk);
164  iEvent.getByLabel("g4SimHits",SimVtx);
165  iEvent.getByLabel("g4SimHits", "EcalHitsEB", pcaloeb);
166  iEvent.getByLabel("g4SimHits", "EcalHitsEE", pcaloee);
167  iEvent.getByLabel("g4SimHits", "HcalHits", pcalohh);
168  associate = new TrackerHitAssociator(iEvent);
169  }
170 
171  unsigned int nTracks=0;
172  for (trkDetItr = trkCaloDirections.begin(),nTracks=0; trkDetItr != trkCaloDirections.end(); trkDetItr++,nTracks++){
173  const reco::Track* pTrack = &(*(trkDetItr->trkItr));
174  if (spr::goodTrack(pTrack,leadPV,selectionParameters,(myverbose>2)) && trkDetItr->okECAL && trkDetItr->okHCAL) {
175  int nRH_eMipDR=0, nRH_eDR=0, nNearTRKs=0, nRecHitsCone=-99;
176  double distFromHotCell=-99.0, distFromHotCell2=-99.0;
177  int ietaHotCell=-99, iphiHotCell=-99;
178  int ietaHotCell2=-99, iphiHotCell2=-99;
179  GlobalPoint gposHotCell(0.,0.,0.), gposHotCell2(0.,0.,0.);
180  std::vector<DetId> coneRecHitDetIds, coneRecHitDetIds2;
181  std::pair<double, bool> e11x11_20SigP, e15x15_20SigP;
182  double hCone = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL,
183  trkDetItr->pointECAL,
184  a_coneR, trkDetItr->directionHCAL,
185  nRecHitsCone, coneRecHitDetIds,
186  distFromHotCell, ietaHotCell, iphiHotCell,
187  gposHotCell, -1);
188  double hConeHB = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL,
189  trkDetItr->pointECAL,
190  a_coneR, trkDetItr->directionHCAL,
191  nRecHitsCone, coneRecHitDetIds,
192  distFromHotCell, ietaHotCell,
193  iphiHotCell, gposHotCell,
194  (int)(HcalBarrel));
195  double eHCALDR = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL,
196  trkDetItr->pointECAL, a_charIsoR,
197  trkDetItr->directionHCAL, nRecHitsCone,
198  coneRecHitDetIds2, distFromHotCell2,
199  ietaHotCell2, iphiHotCell2, gposHotCell2,
200  -1);
201  double eHCALDRHB = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL,
202  trkDetItr->pointECAL, a_charIsoR,
203  trkDetItr->directionHCAL, nRecHitsCone,
204  coneRecHitDetIds2, distFromHotCell2,
205  ietaHotCell2, iphiHotCell2,
206  gposHotCell2, (int)(HcalBarrel));
207 
208  double conehmaxNearP = spr::chargeIsolationCone(nTracks, trkCaloDirections, a_charIsoR, nNearTRKs, (myverbose>3));
209 
210  double eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle,
211  endcapRecHitsHandle,trkDetItr->pointHCAL,
212  trkDetItr->pointECAL, a_mipR,
213  trkDetItr->directionECAL, nRH_eMipDR);
214  double eECALDR = spr::eCone_ecal(geo, barrelRecHitsHandle,
215  endcapRecHitsHandle,trkDetItr->pointHCAL,
216  trkDetItr->pointECAL, a_neutIsoR,
217  trkDetItr->directionECAL, nRH_eDR);
218  double eMipDR_1= spr::eCone_ecal(geo, barrelRecHitsHandle,
219  endcapRecHitsHandle,trkDetItr->pointHCAL,
220  trkDetItr->pointECAL, a_mipR,
221  trkDetItr->directionECAL, nRH_eMipDR,
222  0.030, 0.150);
223  double eECALDR_1=spr::eCone_ecal(geo, barrelRecHitsHandle,
224  endcapRecHitsHandle,trkDetItr->pointHCAL,
225  trkDetItr->pointECAL, a_neutIsoR,
226  trkDetItr->directionECAL, nRH_eDR,
227  0.030, 0.150);
228  double eMipDR_2= spr::eCone_ecal(geo, barrelRecHitsHandle,
229  endcapRecHitsHandle,trkDetItr->pointHCAL,
230  trkDetItr->pointECAL, a_mipR,
231  trkDetItr->directionECAL, nRH_eMipDR,
232  0.060, 0.300);
233  double eECALDR_2=spr::eCone_ecal(geo, barrelRecHitsHandle,
234  endcapRecHitsHandle,trkDetItr->pointHCAL,
235  trkDetItr->pointECAL, a_neutIsoR,
236  trkDetItr->directionECAL, nRH_eDR,
237  0.060, 0.300);
238 
239  HcalDetId closestCell = (HcalDetId)(trkDetItr->detIdHCAL);
240 
241  edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
242  iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
243 
244  e11x11_20SigP = spr::eECALmatrix(trkDetItr->detIdECAL,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),5,5, 0.060, 0.300, tMinE_,tMaxE_);
245  e15x15_20SigP = spr::eECALmatrix(trkDetItr->detIdECAL,barrelRecHitsHandle,endcapRecHitsHandle, *theEcalChStatus, geo, caloTopology,sevlv.product(),7,7, 0.060, 0.300, tMinE_,tMaxE_);
246 
247  // Fill the tree Branches here
248  t_trackP ->push_back( pTrack->p() );
249  t_trackPt ->push_back( pTrack->pt() );
250  t_trackEta ->push_back( pTrack->momentum().eta() );
251  t_trackPhi ->push_back( pTrack->momentum().phi() );
252  t_trackHcalEta ->push_back( closestCell.ieta() );
253  t_trackHcalPhi ->push_back( closestCell.iphi() );
254  t_hCone ->push_back( hCone);
255  t_conehmaxNearP ->push_back( conehmaxNearP);
256  t_eMipDR ->push_back( eMipDR);
257  t_eECALDR ->push_back( eECALDR);
258  t_eHCALDR ->push_back( eHCALDR);
259  t_e11x11_20Sig ->push_back( e11x11_20SigP.first );
260  t_e15x15_20Sig ->push_back( e15x15_20SigP.first );
261  t_eMipDR_1 ->push_back( eMipDR_1);
262  t_eECALDR_1 ->push_back( eECALDR_1);
263  t_eMipDR_2 ->push_back( eMipDR_2);
264  t_eECALDR_2 ->push_back( eECALDR_2);
265  t_hConeHB ->push_back( hConeHB);
266  t_eHCALDRHB ->push_back( eHCALDRHB);
267 
268  if (myverbose > 0) {
269  std::cout << "Track p " << pTrack->p() << " pt " << pTrack->pt()
270  << " eta " << pTrack->momentum().eta() << " phi "
271  << pTrack->momentum().phi() << " ieta/iphi ("
272  << closestCell.ieta() << ", " << closestCell.iphi()
273  << ") Energy in cone " << hCone << " Charge Isolation "
274  << conehmaxNearP << " eMIP (" << eMipDR << ", "
275  << eMipDR_1 << ", " << eMipDR_2 << ")"
276  << " Neutral isolation (ECAL) (" << eECALDR-eMipDR << ", "
277  << eECALDR_1-eMipDR_1 << ", " << eECALDR_2-eMipDR_2 << ")"
278  << " (ECAL NxN) " << e15x15_20SigP.first-e11x11_20SigP.first
279  << " (HCAL) " << eHCALDR-hCone << std::endl;
280  }
281 
282  if (doMC) {
283  int nSimHits = -999;
284  double hsim;
285  std::map<std::string, double> hsimInfo;
286  std::vector<int> multiplicity;
287  hsim = spr::eCone_hcal(geo, pcalohh, trkDetItr->pointHCAL,
288  trkDetItr->pointECAL, a_coneR,
289  trkDetItr->directionHCAL, nSimHits);
290  hsimInfo = spr::eHCALSimInfoCone(iEvent, pcalohh, SimTk, SimVtx,
291  pTrack, *associate, geo,
292  trkDetItr->pointHCAL,
293  trkDetItr->pointECAL, a_coneR,
294  trkDetItr->directionHCAL,
295  multiplicity);
296 
297  t_hsimInfoMatched ->push_back(hsimInfo["eMatched" ]);
298  t_hsimInfoRest ->push_back(hsimInfo["eRest" ]);
299  t_hsimInfoPhoton ->push_back(hsimInfo["eGamma" ]);
300  t_hsimInfoNeutHad ->push_back(hsimInfo["eNeutralHad"]);
301  t_hsimInfoCharHad ->push_back(hsimInfo["eChargedHad"]);
302  t_hsimInfoPdgMatched->push_back(hsimInfo["pdgMatched" ]);
303  t_hsimInfoTotal ->push_back(hsimInfo["eTotal" ]);
304 
305  t_hsimInfoNMatched ->push_back(multiplicity.at(0));
306  t_hsimInfoNTotal ->push_back(multiplicity.at(1));
307  t_hsimInfoNNeutHad ->push_back(multiplicity.at(2));
308  t_hsimInfoNCharHad ->push_back(multiplicity.at(3));
309  t_hsimInfoNPhoton ->push_back(multiplicity.at(4));
310  t_hsimInfoNRest ->push_back(multiplicity.at(5));
311 
312  t_hsim ->push_back(hsim );
313  t_nSimHits ->push_back(nSimHits );
314 
315  if (myverbose > 0) {
316  std::cout << "Matched (E) " << hsimInfo["eMatched"] << " (N) "
317  << multiplicity.at(0) << " Rest (E) " << hsimInfo["eRest"]
318  << " (N) " << multiplicity.at(1) << " Gamma (E) "
319  << hsimInfo["eGamma"] << " (N) " << multiplicity.at(2)
320  << " Neutral Had (E) " << hsimInfo["eNeutralHad"]
321  << " (N) " << multiplicity.at(3) << " Charged Had (E) "
322  << hsimInfo["eChargedHad"] << " (N) " << multiplicity.at(4)
323  << " Total (E) " << hsimInfo["eTotal"] << " (N) "
324  << multiplicity.at(5) << " PDG " << hsimInfo["pdgMatched"]
325  << " Total E " << hsim << " NHit " << nSimHits <<std::endl;
326  }
327  }
328  }
329  }
330 
331  // delete associate;
332  if (associate) delete associate;
333  tree->Fill();
334 }
335 
336 void IsolatedTracksHcalScale::beginJob() {
337 
338  nEventProc=0;
339 
340 
342  tree = fs->make<TTree>("tree", "tree");
343  tree->SetAutoSave(10000);
344 
345  tree->Branch("t_RunNo" ,&t_RunNo ,"t_RunNo/I");
346  tree->Branch("t_Lumi" ,&t_Lumi ,"t_Lumi/I");
347  tree->Branch("t_Bunch" ,&t_Bunch ,"t_Bunch/I");
348 
349  t_trackP = new std::vector<double>();
350  t_trackPt = new std::vector<double>();
351  t_trackEta = new std::vector<double>();
352  t_trackPhi = new std::vector<double>();
353  t_trackHcalEta = new std::vector<double>();
354  t_trackHcalPhi = new std::vector<double>();
355  t_hCone = new std::vector<double>();
356  t_conehmaxNearP = new std::vector<double>();
357  t_eMipDR = new std::vector<double>();
358  t_eECALDR = new std::vector<double>();
359  t_eHCALDR = new std::vector<double>();
360  t_e11x11_20Sig = new std::vector<double>();
361  t_e15x15_20Sig = new std::vector<double>();
362  t_eMipDR_1 = new std::vector<double>();
363  t_eECALDR_1 = new std::vector<double>();
364  t_eMipDR_2 = new std::vector<double>();
365  t_eECALDR_2 = new std::vector<double>();
366  t_hConeHB = new std::vector<double>();
367  t_eHCALDRHB = new std::vector<double>();
368 
369  tree->Branch("t_trackP", "vector<double>", &t_trackP );
370  tree->Branch("t_trackPt", "vector<double>", &t_trackPt );
371  tree->Branch("t_trackEta", "vector<double>", &t_trackEta );
372  tree->Branch("t_trackPhi", "vector<double>", &t_trackPhi );
373  tree->Branch("t_trackHcalEta", "vector<double>", &t_trackHcalEta );
374  tree->Branch("t_trackHcalPhi", "vector<double>", &t_trackHcalPhi );
375  tree->Branch("t_hCone", "vector<double>", &t_hCone );
376  tree->Branch("t_conehmaxNearP", "vector<double>", &t_conehmaxNearP );
377  tree->Branch("t_eMipDR", "vector<double>", &t_eMipDR );
378  tree->Branch("t_eECALDR", "vector<double>", &t_eECALDR );
379  tree->Branch("t_eHCALDR", "vector<double>", &t_eHCALDR );
380  tree->Branch("t_e11x11_20Sig", "vector<double>", &t_e11x11_20Sig );
381  tree->Branch("t_e15x15_20Sig", "vector<double>", &t_e15x15_20Sig );
382  tree->Branch("t_eMipDR_1", "vector<double>", &t_eMipDR_1 );
383  tree->Branch("t_eECALDR_1", "vector<double>", &t_eECALDR_1 );
384  tree->Branch("t_eMipDR_2", "vector<double>", &t_eMipDR_2 );
385  tree->Branch("t_eECALDR_2", "vector<double>", &t_eECALDR_2 );
386  tree->Branch("t_hConeHB", "vector<double>", &t_hConeHB );
387  tree->Branch("t_eHCALDRHB", "vector<double>", &t_eHCALDRHB );
388 
389  if (doMC) {
390  t_hsimInfoMatched = new std::vector<double>();
391  t_hsimInfoRest = new std::vector<double>();
392  t_hsimInfoPhoton = new std::vector<double>();
393  t_hsimInfoNeutHad = new std::vector<double>();
394  t_hsimInfoCharHad = new std::vector<double>();
395  t_hsimInfoPdgMatched = new std::vector<double>();
396  t_hsimInfoTotal = new std::vector<double>();
397  t_hsimInfoNMatched = new std::vector<int>();
398  t_hsimInfoNTotal = new std::vector<int>();
399  t_hsimInfoNNeutHad = new std::vector<int>();
400  t_hsimInfoNCharHad = new std::vector<int>();
401  t_hsimInfoNPhoton = new std::vector<int>();
402  t_hsimInfoNRest = new std::vector<int>();
403  t_hsim = new std::vector<double>();
404  t_nSimHits = new std::vector<int>();
405 
406  tree->Branch("t_hsimInfoMatched", "vector<double>", &t_hsimInfoMatched );
407  tree->Branch("t_hsimInfoRest", "vector<double>", &t_hsimInfoRest );
408  tree->Branch("t_hsimInfoPhoton", "vector<double>", &t_hsimInfoPhoton );
409  tree->Branch("t_hsimInfoNeutHad", "vector<double>", &t_hsimInfoNeutHad );
410  tree->Branch("t_hsimInfoCharHad", "vector<double>", &t_hsimInfoCharHad );
411  tree->Branch("t_hsimInfoPdgMatched", "vector<double>", &t_hsimInfoPdgMatched );
412  tree->Branch("t_hsimInfoTotal", "vector<double>", &t_hsimInfoTotal );
413  tree->Branch("t_hsimInfoNMatched", "vector<int>", &t_hsimInfoNMatched );
414  tree->Branch("t_hsimInfoNTotal", "vector<int>", &t_hsimInfoNTotal );
415  tree->Branch("t_hsimInfoNNeutHad", "vector<int>", &t_hsimInfoNNeutHad );
416  tree->Branch("t_hsimInfoNCharHad", "vector<int>", &t_hsimInfoNCharHad );
417  tree->Branch("t_hsimInfoNPhoton", "vector<int>", &t_hsimInfoNPhoton );
418  tree->Branch("t_hsimInfoNRest", "vector<int>", &t_hsimInfoNRest );
419  tree->Branch("t_hsim", "vector<double>", &t_hsim );
420  tree->Branch("t_nSimHits", "vector<int>", &t_nSimHits );
421  }
422 }
423 
424 void IsolatedTracksHcalScale::endJob() {
425 
426  std::cout << "Number of Events Processed " << nEventProc << std::endl;
427 }
428 
429 void IsolatedTracksHcalScale::clearTreeVectors() {
430 
431  t_trackP ->clear();
432  t_trackPt ->clear();
433  t_trackEta ->clear();
434  t_trackPhi ->clear();
435  t_trackHcalEta ->clear();
436  t_trackHcalPhi ->clear();
437  t_hCone ->clear();
438  t_conehmaxNearP ->clear();
439  t_eMipDR ->clear();
440  t_eECALDR ->clear();
441  t_eHCALDR ->clear();
442  t_e11x11_20Sig ->clear();
443  t_e15x15_20Sig ->clear();
444  t_eMipDR_1 ->clear();
445  t_eECALDR_1 ->clear();
446  t_eMipDR_2 ->clear();
447  t_eECALDR_2 ->clear();
448  t_hConeHB ->clear();
449  t_eHCALDRHB ->clear();
450 
451  if (doMC) {
452  t_hsimInfoMatched ->clear();
453  t_hsimInfoRest ->clear();
454  t_hsimInfoPhoton ->clear();
455  t_hsimInfoNeutHad ->clear();
456  t_hsimInfoCharHad ->clear();
457  t_hsimInfoPdgMatched ->clear();
458  t_hsimInfoTotal ->clear();
459  t_hsimInfoNMatched ->clear();
460  t_hsimInfoNTotal ->clear();
461  t_hsimInfoNNeutHad ->clear();
462  t_hsimInfoNCharHad ->clear();
463  t_hsimInfoNPhoton ->clear();
464  t_hsimInfoNRest ->clear();
465  t_hsim ->clear();
466  t_nSimHits ->clear();
467  }
468 }
469 
470 
471 //define this as a plug-in
472 DEFINE_FWK_MODULE(IsolatedTracksHcalScale);
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:39
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:591
IsolatedTracksHcalScale::t_hsimInfoMatched
std::vector< double > * t_hsimInfoMatched
Definition: IsolatedTracksHcalScale.h:170
IsolatedTracksHcalScale::fs
edm::Service< TFileService > fs
Definition: IsolatedTracksHcalScale.h:176
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:41
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
IsolatedTracksHcalScale::t_eHCALDR
std::vector< double > * t_eHCALDR
Definition: IsolatedTracksHcalScale.h:165
IsolatedTracksHcalScale::t_hsimInfoNRest
std::vector< int > * t_hsimInfoNRest
Definition: IsolatedTracksHcalScale.h:174
IsolatedTracksHcalScale::t_hCone
std::vector< double > * t_hCone
Definition: IsolatedTracksHcalScale.h:166
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:18
IsolatedTracksHcalScale::IsolatedTracksHcalScale
IsolatedTracksHcalScale(const edm::ParameterSet &)
Definition: IsolatedTracksHcalScale.cc:17
IsolatedTracksHcalScale::t_hsimInfoNTotal
std::vector< int > * t_hsimInfoNTotal
Definition: IsolatedTracksHcalScale.h:173
IsolatedTracksHcalScale::t_hsimInfoNCharHad
std::vector< int > * t_hsimInfoNCharHad
Definition: IsolatedTracksHcalScale.h:174
IsolatedTracksHcalScale::t_hsimInfoNMatched
std::vector< int > * t_hsimInfoNMatched
Definition: IsolatedTracksHcalScale.h:173
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:113
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
IsolatedTracksHcalScale::t_hsim
std::vector< double > * t_hsim
Definition: IsolatedTracksHcalScale.h:172
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:62
IsolatedTracksHcalScale::t_hsimInfoNeutHad
std::vector< double > * t_hsimInfoNeutHad
Definition: IsolatedTracksHcalScale.h:171
edm::EventBase::luminosityBlock
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:59
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
IsolatedTracksHcalScale::t_eMipDR_2
std::vector< double > * t_eMipDR_2
Definition: IsolatedTracksHcalScale.h:168
IsolatedTracksHcalScale::t_hsimInfoNPhoton
std::vector< int > * t_hsimInfoNPhoton
Definition: IsolatedTracksHcalScale.h:174
spr::chargeIsolationCone
double chargeIsolationCone(unsigned int trkIndex, std::vector< spr::propagatedTrackDirection > &trkDirs, double dR, int &nNearTRKs, bool debug=false)
Definition: ChargeIsolation.cc:232
reco::TrackBase::momentum
const Vector & momentum() const
track momentum vector
Definition: TrackBase.h:651
detailsBasic3DVector::z
float float float z
Definition: extBasic3DVector.h:15
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
spr::eCone_hcal
double eCone_hcal(const CaloGeometry *geo, edm::Handle< T > &hits, const GlobalPoint &hpoint1, const GlobalPoint &point1, double dR, const GlobalVector &trackMom, int &nRecHits, double hbThr=-100, double heThr=-100, double hfThr=-100, double hoThr=-100, double tMin=-500, double tMax=500, int detOnly=-1)
iEvent
int iEvent
Definition: GenABIO.cc:230
IsolatedTracksHcalScale::t_e11x11_20Sig
std::vector< double > * t_e11x11_20Sig
Definition: IsolatedTracksHcalScale.h:167
IsolatedTracksHcalScale::t_hsimInfoNNeutHad
std::vector< int > * t_hsimInfoNNeutHad
Definition: IsolatedTracksHcalScale.h:173
spr::goodTrack
bool goodTrack(const reco::Track *pTrack, math::XYZPoint leadPV, trackSelectionParameters parameters, bool debug=false)
Definition: TrackSelection.cc:8
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:597
IsolatedTracksHcalScale::analyze
void analyze(const edm::Event &, const edm::EventSetup &)
Definition: IsolatedTracksHcalScale.cc:69
IsolatedTracksHcalScale::t_trackEta
std::vector< double > * t_trackEta
Definition: IsolatedTracksHcalScale.h:164
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:36
IsolatedTracksHcalScale::t_eECALDR_2
std::vector< double > * t_eECALDR_2
Definition: IsolatedTracksHcalScale.h:168
IsolatedTracksHcalScale::t_eECALDR
std::vector< double > * t_eECALDR
Definition: IsolatedTracksHcalScale.h:166
IsolatedTracksHcalScale::selectionParameters
spr::trackSelectionParameters selectionParameters
Definition: IsolatedTracksHcalScale.h:131
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:402
IsolatedTracksHcalScale::t_eMipDR_1
std::vector< double > * t_eMipDR_1
Definition: IsolatedTracksHcalScale.h:168
HcalDetId::iphi
int iphi() const
get the cell iphi
Definition: HcalDetId.h:38
IsolatedTracksHcalScale::t_trackPt
std::vector< double > * t_trackPt
Definition: IsolatedTracksHcalScale.h:164
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:102
IsolatedTracksHcalScale::bField
const MagneticField * bField
Definition: IsolatedTracksHcalScale.h:155
edm::Handle::product
T const * product() const
Definition: Handle.h:81
IsolatedTracksHcalScale::t_trackHcalEta
std::vector< double > * t_trackHcalEta
Definition: IsolatedTracksHcalScale.h:165
IsolatedTracksHcalScale::t_trackP
std::vector< double > * t_trackP
Definition: IsolatedTracksHcalScale.h:164
IsolatedTracksHcalScale::t_eECALDR_1
std::vector< double > * t_eECALDR_1
Definition: IsolatedTracksHcalScale.h:168
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
IsolatedTracksHcalScale::t_trackHcalPhi
std::vector< double > * t_trackHcalPhi
Definition: IsolatedTracksHcalScale.h:165
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
IsolatedTracksHcalScale::t_hsimInfoPhoton
std::vector< double > * t_hsimInfoPhoton
Definition: IsolatedTracksHcalScale.h:170
IsolatedTracksHcalScale::t_hsimInfoTotal
std::vector< double > * t_hsimInfoTotal
Definition: IsolatedTracksHcalScale.h:172
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
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)
IsolatedTracksHcalScale::t_hsimInfoCharHad
std::vector< double > * t_hsimInfoCharHad
Definition: IsolatedTracksHcalScale.h:171
spr::trackSelectionParameters::minQuality
reco::TrackBase::TrackQuality minQuality
Definition: TrackSelection.h:22
IsolatedTracksHcalScale::t_trackPhi
std::vector< double > * t_trackPhi
Definition: IsolatedTracksHcalScale.h:164
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
IsolatedTracksHcalScale::t_eHCALDRHB
std::vector< double > * t_eHCALDRHB
Definition: IsolatedTracksHcalScale.h:169
IsolatedTracksHcalScale::t_hsimInfoPdgMatched
std::vector< double > * t_hsimInfoPdgMatched
Definition: IsolatedTracksHcalScale.h:171
x
Definition: DDAxes.h:10
IsolatedTracksHcalScale::t_nSimHits
std::vector< int > * t_nSimHits
Definition: IsolatedTracksHcalScale.h:175
IsolatedTracksHcalScale::t_hConeHB
std::vector< double > * t_hConeHB
Definition: IsolatedTracksHcalScale.h:169
IsolatedTracksHcalScale::t_hsimInfoRest
std::vector< double > * t_hsimInfoRest
Definition: IsolatedTracksHcalScale.h:170
IsolatedTracksHcalScale::t_conehmaxNearP
std::vector< double > * t_conehmaxNearP
Definition: IsolatedTracksHcalScale.h:166
IsolatedTracksHcalScale::t_e15x15_20Sig
std::vector< double > * t_e15x15_20Sig
Definition: IsolatedTracksHcalScale.h:167
IsolatedTracksHcalScale::t_eMipDR
std::vector< double > * t_eMipDR
Definition: IsolatedTracksHcalScale.h:166
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