CMS 3D CMS Logo
EvtPlaneProducer.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: EvtPlaneProducer
4 // Class: EvtPlaneProducer
5 //
13 //
14 // Original Author: Sergey Petrushanko
15 // Created: Fri Jul 11 10:05:00 2008
16 //
17 //
18 
19 // system include files
20 #include <memory>
21 #include <iostream>
22 #include <ctime>
23 #include <cmath>
24 
25 // user include files
26 #include "FWCore/Framework/interface/Frameworkfwd.h"
27 #include "FWCore/Framework/interface/stream/EDProducer.h"
28 #include "FWCore/Framework/interface/Event.h"
29 #include "FWCore/Framework/interface/EventSetup.h"
30 
31 #include "FWCore/Framework/interface/MakerMacros.h"
32 #include "FWCore/Framework/interface/ESHandle.h"
33 #include "FWCore/Framework/interface/ESWatcher.h"
34 #include "FWCore/ParameterSet/interface/ParameterSet.h"
35 
36 #include "DataFormats/Candidate/interface/Candidate.h"
37 #include "DataFormats/HeavyIonEvent/interface/EvtPlane.h"
38 #include "DataFormats/RecoCandidate/interface/RecoChargedCandidate.h"
39 #include "DataFormats/RecoCandidate/interface/RecoChargedCandidateFwd.h"
40 
41 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
42 #include "DataFormats/Common/interface/Ref.h"
43 
44 #include "DataFormats/Common/interface/Handle.h"
45 #include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
46 #include "DataFormats/CastorReco/interface/CastorTower.h"
47 
48 #include "FWCore/ServiceRegistry/interface/Service.h"
49 #include "DataFormats/TrackReco/interface/TrackFwd.h"
50 #include "DataFormats/TrackReco/interface/Track.h"
51 #include "DataFormats/VertexReco/interface/Vertex.h"
52 #include "DataFormats/VertexReco/interface/VertexFwd.h"
53 #include <cstdlib>
54 #include "RecoHI/HiEvtPlaneAlgos/interface/HiEvtPlaneList.h"
55 #include "CondFormats/HIObjects/interface/RPFlatParams.h"
56 #include "CondFormats/DataRecord/interface/HeavyIonRPRcd.h"
57 #include "CondFormats/DataRecord/interface/HeavyIonRcd.h"
58 #include "CondFormats/HIObjects/interface/CentralityTable.h"
59 
60 #include "DataFormats/Math/interface/Point3D.h"
61 #include "DataFormats/Common/interface/RefProd.h"
62 #include "DataFormats/Common/interface/Ref.h"
63 #include "DataFormats/Common/interface/RefVector.h"
64 
65 #include "RecoHI/HiEvtPlaneAlgos/interface/HiEvtPlaneFlatten.h"
66 #include "RecoHI/HiEvtPlaneAlgos/interface/LoadEPDB.h"
67 
68 using namespace std;
69 using namespace hi;
70 
71 //
72 // class decleration
73 //
74 
75 namespace hi {
76  class GenPlane {
77  public:
78  GenPlane(string name, double etaminval1, double etamaxval1, double etaminval2, double etamaxval2, int orderval) {
79  epname = name;
80  etamin1 = etaminval1;
81  etamax1 = etamaxval1;
82  etamin2 = etaminval2;
83  etamax2 = etamaxval2;
84  sumsin = 0;
85  sumcos = 0;
86  sumsinNoWgt = 0;
87  sumcosNoWgt = 0;
88 
89  mult = 0;
90  order = (double)orderval;
91  }
92  ~GenPlane() { ; }
93  void addParticle(double w, double PtOrEt, double s, double c, double eta) {
94  if ((eta >= etamin1 && eta < etamax1) || (etamin2 != etamax2 && eta >= etamin2 && eta < etamax2)) {
95  sumsin += w * s;
96  sumcos += w * c;
97  sumsinNoWgt += s;
98  sumcosNoWgt += c;
99 
100  sumw += fabs(w);
101  sumw2 += w * w;
102  sumPtOrEt += PtOrEt;
103  sumPtOrEt2 += PtOrEt * PtOrEt;
104  ++mult;
105  }
106  }
107 
108  double getAngle(double &ang,
109  double &sv,
110  double &cv,
111  double &svNoWgt,
112  double &cvNoWgt,
113  double &w,
114  double &w2,
115  double &PtOrEt,
116  double &PtOrEt2,
117  uint &epmult) {
118  ang = -10;
119  sv = 0;
120  cv = 0;
121  sv = sumsin;
122  cv = sumcos;
123  svNoWgt = sumsinNoWgt;
124  cvNoWgt = sumcosNoWgt;
125  w = sumw;
126  w2 = sumw2;
127  PtOrEt = sumPtOrEt;
128  PtOrEt2 = sumPtOrEt2;
129  epmult = mult;
130  double q = sv * sv + cv * cv;
131  if (q > 0)
132  ang = atan2(sv, cv) / order;
133  return ang;
134  }
135  void reset() {
136  sumsin = 0;
137  sumcos = 0;
138  sumsinNoWgt = 0;
139  sumcosNoWgt = 0;
140  sumw = 0;
141  sumw2 = 0;
142  mult = 0;
143  sumPtOrEt = 0;
144  sumPtOrEt2 = 0;
145  }
146 
147  private:
148  string epname;
149  double etamin1;
150  double etamax1;
151 
152  double etamin2;
153  double etamax2;
154  double sumsin;
155  double sumcos;
156  double sumsinNoWgt;
157  double sumcosNoWgt;
158  uint mult;
159  double sumw;
160  double sumw2;
161  double sumPtOrEt;
162  double sumPtOrEt2;
163  double order;
164  };
165 } // namespace hi
166 
167 class EvtPlaneProducer : public edm::stream::EDProducer<> {
168 public:
169  explicit EvtPlaneProducer(const edm::ParameterSet &);
170  ~EvtPlaneProducer() override;
171 
172 private:
173  GenPlane *rp[NumEPNames];
174 
175  void produce(edm::Event &, const edm::EventSetup &) override;
176 
177  // ----------member data ---------------------------
178 
179  std::string centralityVariable_;
180  std::string centralityLabel_;
181  std::string centralityMC_;
182 
183  edm::InputTag centralityBinTag_;
184  edm::EDGetTokenT<int> centralityBinToken;
185 
186  edm::InputTag vertexTag_;
187  edm::EDGetTokenT<std::vector<reco::Vertex>> vertexToken;
188  edm::Handle<std::vector<reco::Vertex>> vertex_;
189 
190  edm::InputTag caloTag_;
191  edm::EDGetTokenT<CaloTowerCollection> caloToken;
192  edm::Handle<CaloTowerCollection> caloCollection_;
193 
194  edm::InputTag castorTag_;
195  edm::EDGetTokenT<std::vector<reco::CastorTower>> castorToken;
196  edm::Handle<std::vector<reco::CastorTower>> castorCollection_;
197 
198  edm::InputTag trackTag_;
199  edm::EDGetTokenT<reco::TrackCollection> trackToken;
200  edm::Handle<reco::TrackCollection> trackCollection_;
201 
202  edm::ESWatcher<HeavyIonRcd> hiWatcher;
203  edm::ESWatcher<HeavyIonRPRcd> hirpWatcher;
204 
205  bool loadDB_;
206  double minet_;
207  double maxet_;
208  double minpt_;
209  double maxpt_;
210  double minvtx_;
211  double maxvtx_;
212  double dzerr_;
213  double chi2_;
214  int FlatOrder_;
215  int NumFlatBins_;
216  double nCentBins_;
217  double caloCentRef_;
218  double caloCentRefWidth_;
219  int CentBinCompression_;
220  HiEvtPlaneFlatten *flat[NumEPNames];
221 };
222 
223 EvtPlaneProducer::EvtPlaneProducer(const edm::ParameterSet &iConfig)
224  : centralityVariable_(iConfig.getParameter<std::string>("centralityVariable")),
225  centralityBinTag_(iConfig.getParameter<edm::InputTag>("centralityBinTag")),
226  vertexTag_(iConfig.getParameter<edm::InputTag>("vertexTag")),
227  caloTag_(iConfig.getParameter<edm::InputTag>("caloTag")),
228  castorTag_(iConfig.getParameter<edm::InputTag>("castorTag")),
229  trackTag_(iConfig.getParameter<edm::InputTag>("trackTag")),
230  loadDB_(iConfig.getParameter<bool>("loadDB")),
231  minet_(iConfig.getParameter<double>("minet")),
232  maxet_(iConfig.getParameter<double>("maxet")),
233  minpt_(iConfig.getParameter<double>("minpt")),
234  maxpt_(iConfig.getParameter<double>("maxpt")),
235  minvtx_(iConfig.getParameter<double>("minvtx")),
236  maxvtx_(iConfig.getParameter<double>("maxvtx")),
237  dzerr_(iConfig.getParameter<double>("dzerr")),
238  chi2_(iConfig.getParameter<double>("chi2")),
239  FlatOrder_(iConfig.getParameter<int>("FlatOrder")),
240  NumFlatBins_(iConfig.getParameter<int>("NumFlatBins")),
241  caloCentRef_(iConfig.getParameter<double>("caloCentRef")),
242  caloCentRefWidth_(iConfig.getParameter<double>("caloCentRefWidth")),
243  CentBinCompression_(iConfig.getParameter<int>("CentBinCompression")) {
244  nCentBins_ = 200.;
245 
246  if (iConfig.exists("nonDefaultGlauberModel")) {
247  centralityMC_ = iConfig.getParameter<std::string>("nonDefaultGlauberModel");
248  }
249  centralityLabel_ = centralityVariable_ + centralityMC_;
250 
251  centralityBinToken = consumes<int>(centralityBinTag_);
252 
253  vertexToken = consumes<std::vector<reco::Vertex>>(vertexTag_);
254 
255  caloToken = consumes<CaloTowerCollection>(caloTag_);
256 
257  castorToken = consumes<std::vector<reco::CastorTower>>(castorTag_);
258 
259  trackToken = consumes<reco::TrackCollection>(trackTag_);
260 
261  produces<reco::EvtPlaneCollection>();
262  for (int i = 0; i < NumEPNames; i++) {
263  rp[i] = new GenPlane(EPNames[i], EPEtaMin1[i], EPEtaMax1[i], EPEtaMin2[i], EPEtaMax2[i], EPOrder[i]);
264  }
265  for (int i = 0; i < NumEPNames; i++) {
266  flat[i] = new HiEvtPlaneFlatten();
267  flat[i]->init(FlatOrder_, NumFlatBins_, EPNames[i], EPOrder[i]);
268  }
269 }
270 
271 EvtPlaneProducer::~EvtPlaneProducer() {
272  // do anything here that needs to be done at desctruction time
273  // (e.g. close files, deallocate resources etc.)
274  for (int i = 0; i < NumEPNames; i++) {
275  delete flat[i];
276  }
277 }
278 
279 //
280 // member functions
281 //
282 
283 // ------------ method called to produce the data ------------
284 void EvtPlaneProducer::produce(edm::Event &iEvent, const edm::EventSetup &iSetup) {
285  using namespace edm;
286  using namespace std;
287  using namespace reco;
288 
289  if (loadDB_ && (hiWatcher.check(iSetup) || hirpWatcher.check(iSetup))) {
290  //
291  //Get Size of Centrality Table
292  //
293  edm::ESHandle<CentralityTable> centDB_;
294  iSetup.get<HeavyIonRcd>().get(centralityLabel_, centDB_);
295  nCentBins_ = centDB_->m_table.size();
296  for (int i = 0; i < NumEPNames; i++) {
297  if (caloCentRef_ > 0) {
298  int minbin = (caloCentRef_ - caloCentRefWidth_ / 2.) * nCentBins_ / 100.;
299  int maxbin = (caloCentRef_ + caloCentRefWidth_ / 2.) * nCentBins_ / 100.;
300  minbin /= CentBinCompression_;
301  maxbin /= CentBinCompression_;
302  if (minbin > 0 && maxbin >= minbin) {
303  if (EPDet[i] == HF || EPDet[i] == Castor)
304  flat[i]->setCaloCentRefBins(minbin, maxbin);
305  }
306  }
307  }
308 
309  //
310  //Get flattening parameter file.
311  //
312  if (loadDB_) {
313  edm::ESHandle<RPFlatParams> flatparmsDB_;
314  iSetup.get<HeavyIonRPRcd>().get(flatparmsDB_);
315  LoadEPDB db(flatparmsDB_, flat);
316  if (!db.IsSuccess()) {
317  loadDB_ = kFALSE;
318  }
319  }
320 
321  } //rp record change
322 
323  //
324  //Get Centrality
325  //
326  int bin = 0;
327  if (loadDB_) {
328  edm::Handle<int> cbin_;
329  iEvent.getByToken(centralityBinToken, cbin_);
330  int cbin = *cbin_;
331  bin = cbin / CentBinCompression_;
332  }
333  //
334  //Get Vertex
335  //
336  int vs_sell = 0.;
337  float vzr_sell;
338  iEvent.getByToken(vertexToken, vertex_);
339  const reco::VertexCollection *vertices3 = nullptr;
340  if (vertex_.isValid()) {
341  vertices3 = vertex_.product();
342  vs_sell = vertices3->size();
343  }
344  if (vs_sell > 0) {
345  vzr_sell = vertices3->begin()->z();
346  } else
347  vzr_sell = -999.9;
348  //
349  for (int i = 0; i < NumEPNames; i++)
350  rp[i]->reset();
351  if (vzr_sell < minvtx_ or vzr_sell > maxvtx_)
352  return;
353 
354  //calorimetry part
355 
356  double tower_eta, tower_phi;
357  double tower_energyet, tower_energyet_e, tower_energyet_h;
358 
359  iEvent.getByToken(caloToken, caloCollection_);
360 
361  if (caloCollection_.isValid()) {
362  for (CaloTowerCollection::const_iterator j = caloCollection_->begin(); j != caloCollection_->end(); j++) {
363  tower_eta = j->eta();
364  tower_phi = j->phi();
365  tower_energyet_e = j->emEt();
366  tower_energyet_h = j->hadEt();
367  tower_energyet = tower_energyet_e + tower_energyet_h;
368  double minet = minet_;
369  double maxet = maxet_;
370  for (int i = 0; i < NumEPNames; i++) {
371  if (minet_ < 0)
372  minet = minTransverse[i];
373  if (maxet_ < 0)
374  maxet = maxTransverse[i];
375  if (tower_energyet < minet)
376  continue;
377  if (tower_energyet > maxet)
378  continue;
379  if (EPDet[i] == HF) {
380  double w = tower_energyet;
381  if (loadDB_)
382  w = tower_energyet * flat[i]->getEtScale(vzr_sell, bin);
383  if (EPOrder[i] == 1) {
384  if (MomConsWeight[i][0] == 'y' && loadDB_) {
385  w = flat[i]->getW(tower_energyet, vzr_sell, bin);
386  }
387  if (tower_eta < 0)
388  w = -w;
389  }
390  rp[i]->addParticle(w, tower_energyet, sin(EPOrder[i] * tower_phi), cos(EPOrder[i] * tower_phi), tower_eta);
391  }
392  }
393  }
394  }
395 
396  //Castor part
397 
398  iEvent.getByToken(castorToken, castorCollection_);
399 
400  if (castorCollection_.isValid()) {
401  for (std::vector<reco::CastorTower>::const_iterator j = castorCollection_->begin(); j != castorCollection_->end();
402  j++) {
403  tower_eta = j->eta();
404  tower_phi = j->phi();
405  tower_energyet = j->et();
406  double minet = minet_;
407  double maxet = maxet_;
408  for (int i = 0; i < NumEPNames; i++) {
409  if (EPDet[i] == Castor) {
410  if (minet_ < 0)
411  minet = minTransverse[i];
412  if (maxet_ < 0)
413  maxet = maxTransverse[i];
414  if (tower_energyet < minet)
415  continue;
416  if (tower_energyet > maxet)
417  continue;
418  double w = tower_energyet;
419  if (EPOrder[i] == 1) {
420  if (MomConsWeight[i][0] == 'y' && loadDB_) {
421  w = flat[i]->getW(tower_energyet, vzr_sell, bin);
422  }
423  if (tower_eta < 0)
424  w = -w;
425  }
426  rp[i]->addParticle(w, tower_energyet, sin(EPOrder[i] * tower_phi), cos(EPOrder[i] * tower_phi), tower_eta);
427  }
428  }
429  }
430  }
431 
432  //Tracking part
433 
434  double track_eta;
435  double track_phi;
436  double track_pt;
437 
438  double vzErr2 = 0.0, vxyErr = 0.0;
439  math::XYZPoint vtxPoint(0.0, 0.0, 0.0);
440  if (vertex_.isValid() && !vertex_->empty()) {
441  vtxPoint = vertex_->begin()->position();
442  vzErr2 = (vertex_->begin()->zError()) * (vertex_->begin()->zError());
443  vxyErr = vertex_->begin()->xError() * vertex_->begin()->yError();
444  }
445 
446  iEvent.getByToken(trackToken, trackCollection_);
447  if (trackCollection_.isValid()) {
448  for (reco::TrackCollection::const_iterator j = trackCollection_->begin(); j != trackCollection_->end(); j++) {
449  bool accepted = true;
450  bool isPixel = false;
451  // determine if the track is a pixel track
452  if (j->numberOfValidHits() < 7)
453  isPixel = true;
454 
455  // determine the vertex significance
456  double d0 = 0.0, dz = 0.0, d0sigma = 0.0, dzsigma = 0.0;
457  d0 = -1. * j->dxy(vtxPoint);
458  dz = j->dz(vtxPoint);
459  d0sigma = sqrt(j->d0Error() * j->d0Error() + vxyErr);
460  dzsigma = sqrt(j->dzError() * j->dzError() + vzErr2);
461 
462  // cuts for pixel tracks
463  if (isPixel) {
464  // dz significance cut
465  if (fabs(dz / dzsigma) > dzerr_)
466  accepted = false;
467  // chi2/ndof cut
468  if (j->normalizedChi2() > chi2_)
469  accepted = false;
470  }
471  // cuts for full tracks
472  if (!isPixel) {
473  // dz and d0 significance cuts
474  if (fabs(dz / dzsigma) > 3)
475  accepted = false;
476  if (fabs(d0 / d0sigma) > 3)
477  accepted = false;
478  // pt resolution cut
479  if (j->ptError() / j->pt() > 0.1)
480  accepted = false;
481  // number of valid hits cut
482  if (j->numberOfValidHits() < 12)
483  accepted = false;
484  }
485  if (accepted) {
486  track_eta = j->eta();
487  track_phi = j->phi();
488  track_pt = j->pt();
489  double minpt = minpt_;
490  double maxpt = maxpt_;
491  for (int i = 0; i < NumEPNames; i++) {
492  if (minpt_ < 0)
493  minpt = minTransverse[i];
494  if (maxpt_ < 0)
495  maxpt = maxTransverse[i];
496  if (track_pt < minpt)
497  continue;
498  if (track_pt > maxpt)
499  continue;
500  if (EPDet[i] == Tracker) {
501  double w = track_pt;
502  if (w > 2.5)
503  w = 2.0; //v2 starts decreasing above ~2.5 GeV/c
504  if (EPOrder[i] == 1) {
505  if (MomConsWeight[i][0] == 'y' && loadDB_) {
506  w = flat[i]->getW(track_pt, vzr_sell, bin);
507  }
508  if (track_eta < 0)
509  w = -w;
510  }
511  rp[i]->addParticle(w, track_pt, sin(EPOrder[i] * track_phi), cos(EPOrder[i] * track_phi), track_eta);
512  }
513  }
514  }
515  } //end for
516  }
517 
518  auto evtplaneOutput = std::make_unique<EvtPlaneCollection>();
519 
520  double ang = -10;
521  double sv = 0;
522  double cv = 0;
523  double svNoWgt = 0;
524  double cvNoWgt = 0;
525 
526  double wv = 0;
527  double wv2 = 0;
528  double pe = 0;
529  double pe2 = 0;
530  uint epmult = 0;
531 
532  for (int i = 0; i < NumEPNames; i++) {
533  rp[i]->getAngle(ang, sv, cv, svNoWgt, cvNoWgt, wv, wv2, pe, pe2, epmult);
534  evtplaneOutput->push_back(EvtPlane(i, 0, ang, sv, cv, wv, wv2, pe, pe2, epmult));
535  evtplaneOutput->back().addLevel(3, 0., svNoWgt, cvNoWgt);
536  }
537 
538  iEvent.put(std::move(evtplaneOutput));
539 }
540 
541 //define this as a plug-in
542 DEFINE_FWK_MODULE(EvtPlaneProducer);
EvtPlaneProducer::rp
GenPlane * rp[NumEPNames]
Definition: EvtPlaneProducer.cc:173
EvtPlaneProducer::caloToken
edm::EDGetTokenT< CaloTowerCollection > caloToken
Definition: EvtPlaneProducer.cc:191
edm::ESWatcher::check
bool check(const edm::EventSetup &iSetup)
Definition: ESWatcher.h:52
EvtPlaneProducer::~EvtPlaneProducer
~EvtPlaneProducer() override
Definition: EvtPlaneProducer.cc:271
RefProd.h
EvtPlane.h
Handle.h
w2
common ppss p3p6s2 common epss epspn46 common const1 w2
Definition: inclppp.h:1
hi::GenPlane::sumw
double sumw
Definition: EvtPlaneProducer.cc:159
electrons_cff.bool
bool
Definition: electrons_cff.py:372
mps_fire.i
i
Definition: mps_fire.py:355
hi::MomConsWeight
const std::string MomConsWeight[]
Definition: HiEvtPlaneList.h:113
edm::SortedCollection< CaloTower >::const_iterator
std::vector< CaloTower >::const_iterator const_iterator
Definition: SortedCollection.h:80
EvtPlaneProducer::centralityBinTag_
edm::InputTag centralityBinTag_
Definition: EvtPlaneProducer.cc:183
CastorTower.h
EvtPlaneProducer::FlatOrder_
int FlatOrder_
Definition: EvtPlaneProducer.cc:214
HiEvtPlaneFlatten::setCaloCentRefBins
void setCaloCentRefBins(const int caloCentRefMinBin, const int caloCentRefMaxBin)
Definition: HiEvtPlaneFlatten.h:128
EvtPlaneProducer::hiWatcher
edm::ESWatcher< HeavyIonRcd > hiWatcher
Definition: EvtPlaneProducer.cc:202
edm::Handle::product
T const * product() const
Definition: Handle.h:70
hi::GenPlane::GenPlane
GenPlane(string name, double etaminval1, double etamaxval1, double etaminval2, double etamaxval2, int orderval)
Definition: EvtPlaneProducer.cc:78
hi::EPNames
const std::string EPNames[]
Definition: HiEvtPlaneList.h:73
edm::ESWatcher< HeavyIonRcd >
hi::GenPlane
Definition: EvtPlaneProducer.cc:76
ESHandle.h
HiEvtPlaneFlatten::getW
double getW(double pt, double vtx, int centbin) const
Definition: HiEvtPlaneFlatten.h:149
cuy.cv
cv
Definition: cuy.py:364
hi::GenPlane::sumw2
double sumw2
Definition: EvtPlaneProducer.cc:160
hi::GenPlane::addParticle
void addParticle(double w, double PtOrEt, double s, double c, double eta)
Definition: EvtPlaneProducer.cc:93
edm::EDGetTokenT< int >
edm
HLT enums.
Definition: AlignableModifier.h:19
EvtPlaneProducer::centralityBinToken
edm::EDGetTokenT< int > centralityBinToken
Definition: EvtPlaneProducer.cc:184
CentralityTable.h
HiEvtPlaneFlatten
Definition: HiEvtPlaneFlatten.h:27
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
EvtPlaneProducer::maxpt_
double maxpt_
Definition: EvtPlaneProducer.cc:209
LoadEPDB
Definition: LoadEPDB.h:30
EvtPlaneProducer::dzerr_
double dzerr_
Definition: EvtPlaneProducer.cc:212
hi::GenPlane::epname
string epname
Definition: EvtPlaneProducer.cc:148
HeavyIonRPRcd
Definition: HeavyIonRPRcd.h:24
hi::GenPlane::sumsinNoWgt
double sumsinNoWgt
Definition: EvtPlaneProducer.cc:156
data-class-funcs.q
q
Definition: data-class-funcs.py:169
EvtPlaneProducer::loadDB_
bool loadDB_
Definition: EvtPlaneProducer.cc:205
EvtPlaneProducer::minet_
double minet_
Definition: EvtPlaneProducer.cc:206
hi::EPDet
const int EPDet[]
Definition: HiEvtPlaneList.h:81
EDProducer.h
hi::EPEtaMax2
const double EPEtaMax2[]
Definition: HiEvtPlaneList.h:98
EvtPlaneProducer::trackTag_
edm::InputTag trackTag_
Definition: EvtPlaneProducer.cc:198
HeavyIonRPRcd.h
EvtPlaneProducer::trackToken
edm::EDGetTokenT< reco::TrackCollection > trackToken
Definition: EvtPlaneProducer.cc:199
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
edm::Handle
Definition: AssociativeIterator.h:50
HiEvtPlane_cfi.minpt
minpt
Definition: HiEvtPlane_cfi.py:19
parallelization.uint
uint
Definition: parallelization.py:124
hi::EPOrder
const int EPOrder[]
Definition: HiEvtPlaneList.h:85
HiEvtPlane_cfi.minet
minet
Definition: HiEvtPlane_cfi.py:17
EvtPlaneProducer::caloTag_
edm::InputTag caloTag_
Definition: EvtPlaneProducer.cc:190
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
align::Tracker
Definition: StructureType.h:70
HeavyIonRcd.h
EvtPlaneProducer::flat
HiEvtPlaneFlatten * flat[NumEPNames]
Definition: EvtPlaneProducer.cc:220
dqmiodatasetharvest.db
db
Definition: dqmiodatasetharvest.py:119
MakerMacros.h
alignCSCRings.s
s
Definition: alignCSCRings.py:92
EvtPlaneProducer::chi2_
double chi2_
Definition: EvtPlaneProducer.cc:213
EvtPlaneProducer::caloCentRefWidth_
double caloCentRefWidth_
Definition: EvtPlaneProducer.cc:218
EvtPlaneProducer::hirpWatcher
edm::ESWatcher< HeavyIonRPRcd > hirpWatcher
Definition: EvtPlaneProducer.cc:203
EvtPlaneProducer::centralityLabel_
std::string centralityLabel_
Definition: EvtPlaneProducer.cc:180
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
Track.h
hi::GenPlane::reset
void reset()
Definition: EvtPlaneProducer.cc:135
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
TrackFwd.h
EvtPlaneProducer::centralityVariable_
std::string centralityVariable_
Definition: EvtPlaneProducer.cc:179
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
Service.h
PVValHelper::eta
Definition: PVValidationHelpers.h:69
w
const double w
Definition: UKUtility.cc:23
EvtPlaneProducer::centralityMC_
std::string centralityMC_
Definition: EvtPlaneProducer.cc:181
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
EvtPlaneProducer::produce
void produce(edm::Event &, const edm::EventSetup &) override
Definition: EvtPlaneProducer.cc:284
cms::dd::accepted
bool accepted(std::vector< std::regex > const &, std::string_view)
pfDeepBoostedJetPreprocessParams_cfi.sv
sv
Definition: pfDeepBoostedJetPreprocessParams_cfi.py:226
CentralityTable::m_table
std::vector< CBin > m_table
Definition: CentralityTable.h:39
edm::ESHandle< CentralityTable >
EvtPlaneProducer::nCentBins_
double nCentBins_
Definition: EvtPlaneProducer.cc:216
EvtPlaneProducer::CentBinCompression_
int CentBinCompression_
Definition: EvtPlaneProducer.cc:219
hi::GenPlane::order
double order
Definition: EvtPlaneProducer.cc:163
EvtPlaneProducer::vertexTag_
edm::InputTag vertexTag_
Definition: EvtPlaneProducer.cc:186
EvtPlaneProducer::castorToken
edm::EDGetTokenT< std::vector< reco::CastorTower > > castorToken
Definition: EvtPlaneProducer.cc:195
hi::GenPlane::~GenPlane
~GenPlane()
Definition: EvtPlaneProducer.cc:92
EvtPlaneProducer::vertexToken
edm::EDGetTokenT< std::vector< reco::Vertex > > vertexToken
Definition: EvtPlaneProducer.cc:187
EvtPlaneProducer::maxvtx_
double maxvtx_
Definition: EvtPlaneProducer.cc:211
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
Vertex.h
Castor
hi::maxTransverse
const double maxTransverse[]
Definition: HiEvtPlaneList.h:105
EvtPlaneProducer::minpt_
double minpt_
Definition: EvtPlaneProducer.cc:208
hi::GenPlane::mult
uint mult
Definition: EvtPlaneProducer.cc:158
edm::ParameterSet::exists
bool exists(std::string const &parameterName) const
checks if a parameter exists
Definition: ParameterSet.cc:674
EvtPlaneProducer::castorCollection_
edm::Handle< std::vector< reco::CastorTower > > castorCollection_
Definition: EvtPlaneProducer.cc:196
HLT_2018_cff.InputTag
InputTag
Definition: HLT_2018_cff.py:79016
hi::GenPlane::sumsin
double sumsin
Definition: EvtPlaneProducer.cc:154
edm::ParameterSet
Definition: ParameterSet.h:36
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
RefVector.h
Event.h
fastTrackerRecHitType::isPixel
bool isPixel(HitType hitType)
Definition: FastTrackerRecHit.h:37
EvtPlaneProducer
Definition: EvtPlaneProducer.cc:167
hi::GenPlane::sumcos
double sumcos
Definition: EvtPlaneProducer.cc:155
hi::GenPlane::etamax1
double etamax1
Definition: EvtPlaneProducer.cc:150
hi::GenPlane::etamax2
double etamax2
Definition: EvtPlaneProducer.cc:153
reco::EvtPlane
Definition: EvtPlane.h:20
hi::GenPlane::sumPtOrEt
double sumPtOrEt
Definition: EvtPlaneProducer.cc:161
eventshapeDQM_cfi.order
order
Definition: eventshapeDQM_cfi.py:8
LoadEPDB.h
createfilelist.int
int
Definition: createfilelist.py:10
iEvent
int iEvent
Definition: GenABIO.cc:224
hi::GenPlane::etamin2
double etamin2
Definition: EvtPlaneProducer.cc:152
hi::GenPlane::sumcosNoWgt
double sumcosNoWgt
Definition: EvtPlaneProducer.cc:157
edm::stream::EDProducer
Definition: EDProducer.h:38
edm::EventSetup
Definition: EventSetup.h:57
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
HiEvtPlaneList.h
get
#define get
hi::minTransverse
const double minTransverse[]
Definition: HiEvtPlaneList.h:101
EvtPlaneProducer::caloCollection_
edm::Handle< CaloTowerCollection > caloCollection_
Definition: EvtPlaneProducer.cc:192
HiEvtPlaneFlatten::getEtScale
double getEtScale(double vtx, int centbin) const
Definition: HiEvtPlaneFlatten.h:133
CaloTowerCollection.h
newFWLiteAna.bin
bin
Definition: newFWLiteAna.py:161
VertexFwd.h
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
hi
Definition: HiEvtPlaneList.h:38
eostools.move
def move(src, dest)
Definition: eostools.py:511
std
Definition: JetResolutionObject.h:76
Ref.h
EvtPlaneProducer::vertex_
edm::Handle< std::vector< reco::Vertex > > vertex_
Definition: EvtPlaneProducer.cc:188
EvtPlaneProducer::caloCentRef_
double caloCentRef_
Definition: EvtPlaneProducer.cc:217
HiEvtPlane_cfi.maxpt
maxpt
Definition: HiEvtPlane_cfi.py:20
RecoChargedCandidate.h
PVValHelper::dz
Definition: PVValidationHelpers.h:50
Frameworkfwd.h
ESWatcher.h
hi::EPEtaMax1
const double EPEtaMax1[]
Definition: HiEvtPlaneList.h:91
hi::GenPlane::getAngle
double getAngle(double &ang, double &sv, double &cv, double &svNoWgt, double &cvNoWgt, double &w, double &w2, double &PtOrEt, double &PtOrEt2, uint &epmult)
Definition: EvtPlaneProducer.cc:108
EvtPlaneProducer::castorTag_
edm::InputTag castorTag_
Definition: EvtPlaneProducer.cc:194
CaloTPGTranscoder_cfi.HF
HF
Definition: CaloTPGTranscoder_cfi.py:4
EvtPlaneProducer::NumFlatBins_
int NumFlatBins_
Definition: EvtPlaneProducer.cc:215
Point3D.h
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
EventSetup.h
HcalRecHitCollections.h
VarParsing.mult
mult
Definition: VarParsing.py:659
HiEvtPlaneFlatten::init
void init(int order, int nbins, std::string tag, int vord)
Definition: HiEvtPlaneFlatten.h:35
hi::EPEtaMin2
const double EPEtaMin2[]
Definition: HiEvtPlaneList.h:95
reset
void reset(double vett[256])
Definition: TPedValues.cc:11
hi::NumEPNames
static const int NumEPNames
Definition: HiEvtPlaneList.h:126
EvtPlaneProducer::maxet_
double maxet_
Definition: EvtPlaneProducer.cc:207
Candidate.h
EvtPlaneProducer::trackCollection_
edm::Handle< reco::TrackCollection > trackCollection_
Definition: EvtPlaneProducer.cc:200
HeavyIonRcd
Definition: HeavyIonRcd.h:24
ParameterSet.h
hi::EPEtaMin1
const double EPEtaMin1[]
Definition: HiEvtPlaneList.h:87
hi::GenPlane::etamin1
double etamin1
Definition: EvtPlaneProducer.cc:149
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
hi::GenPlane::sumPtOrEt2
double sumPtOrEt2
Definition: EvtPlaneProducer.cc:162
edm::Event
Definition: Event.h:73
EvtPlaneProducer::EvtPlaneProducer
EvtPlaneProducer(const edm::ParameterSet &)
Definition: EvtPlaneProducer.cc:223
RPFlatParams.h
RecoChargedCandidateFwd.h
d0
static constexpr float d0
Definition: L1EGammaCrystalsEmulatorProducer.cc:85
edm::InputTag
Definition: InputTag.h:15
HiEvtPlane_cfi.maxet
maxet
Definition: HiEvtPlane_cfi.py:18
EvtPlaneProducer::minvtx_
double minvtx_
Definition: EvtPlaneProducer.cc:210
HiEvtPlaneFlatten.h
Generated for CMSSW Reference Manual by doxygen 1.8.16