CMS 3D CMS Logo
CentralityProducer.cc
Go to the documentation of this file.
1 //
2 // Original Author: Yetkin Yilmaz, Young Soo Park
3 // Created: Wed Jun 11 15:31:41 CEST 2008
4 //
5 //
6 
7 // system include files
8 #include <memory>
9 #include <iostream>
10 
11 // user include files
12 #include "FWCore/Framework/interface/Frameworkfwd.h"
13 #include "FWCore/Framework/interface/global/EDProducer.h"
14 #include "FWCore/Utilities/interface/StreamID.h"
15 #include "FWCore/Framework/interface/Event.h"
16 #include "FWCore/Framework/interface/EventSetup.h"
17 
18 #include "FWCore/Framework/interface/MakerMacros.h"
19 #include "FWCore/Framework/interface/ESHandle.h"
20 #include "FWCore/ParameterSet/interface/ParameterSet.h"
21 #include "FWCore/Utilities/interface/InputTag.h"
22 #include "FWCore/Utilities/interface/ESGetToken.h"
23 
24 #include "DataFormats/Candidate/interface/Candidate.h"
25 
26 #include "DataFormats/HeavyIonEvent/interface/Centrality.h"
27 #include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
28 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
29 
30 #include "DataFormats/Common/interface/Ref.h"
31 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHitCollection.h"
32 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHit.h"
33 #include "DataFormats/TrackReco/interface/Track.h"
34 #include "DataFormats/VertexReco/interface/Vertex.h"
35 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
36 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
37 
38 #include "Geometry/Records/interface/CaloGeometryRecord.h"
39 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
40 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
41 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
42 #include "Geometry/CommonDetUnit/interface/PixelGeomDetUnit.h"
43 
44 using namespace std;
45 using namespace edm;
46 using namespace reco;
47 
48 //
49 // class declaration
50 //
51 
52 namespace reco {
53  class CentralityProducer : public edm::global::EDProducer<> {
54  public:
55  explicit CentralityProducer(const edm::ParameterSet&);
56  ~CentralityProducer() override = default;
57  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
58 
59  private:
60  void beginJob() override;
61  void produce(edm::StreamID, edm::Event&, const edm::EventSetup&) const override;
62  void endJob() override;
63 
64  // ----------member data ---------------------------
65 
66  const bool produceHFhits_;
67  const bool produceHFtowers_;
68  const bool produceEcalhits_;
69  const bool produceZDChits_;
70  const bool lowGainZDC_;
71  const bool produceETmidRap_;
72  const bool producePixelhits_;
73  const bool doPixelCut_;
74  const bool produceTracks_;
75  const bool producePixelTracks_;
76 
77  const double midRapidityRange_;
78  const double trackPtCut_;
79  const double trackEtaCut_;
80  const double hfEtaCut_;
81 
82  const bool reuseAny_;
83  const bool useQuality_;
84  const reco::TrackBase::TrackQuality trackQuality_;
85 
86  const edm::EDGetTokenT<HFRecHitCollection> srcHFhits_;
87  const edm::EDGetTokenT<CaloTowerCollection> srcTowers_;
88  const edm::EDGetTokenT<EcalRecHitCollection> srcEBhits_;
89  const edm::EDGetTokenT<EcalRecHitCollection> srcEEhits_;
90  const edm::EDGetTokenT<ZDCRecHitCollection> srcZDChits_;
91  const edm::EDGetTokenT<SiPixelRecHitCollection> srcPixelhits_;
92  const edm::EDGetTokenT<TrackCollection> srcTracks_;
93  const edm::EDGetTokenT<TrackCollection> srcPixelTracks_;
94  const edm::EDGetTokenT<VertexCollection> srcVertex_;
95  const edm::EDGetTokenT<Centrality> reuseTag_;
96 
97  const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> caloGeom_;
98  const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> trackerGeom_;
99  const edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> trackerTopo_;
100  };
101 
102  //
103  // constants, enums and typedefs
104  //
105 
106  //
107  // static data member definitions
108  //
109 
110  //
111  // constructors and destructor
112  //
113  CentralityProducer::CentralityProducer(const edm::ParameterSet& iConfig)
114  : produceHFhits_(iConfig.getParameter<bool>("produceHFhits")),
115  produceHFtowers_(iConfig.getParameter<bool>("produceHFtowers")),
116  produceEcalhits_(iConfig.getParameter<bool>("produceEcalhits")),
117  produceZDChits_(iConfig.getParameter<bool>("produceZDChits")),
118  lowGainZDC_(iConfig.getParameter<bool>("lowGainZDC")),
119  produceETmidRap_(iConfig.getParameter<bool>("produceETmidRapidity")),
120  producePixelhits_(iConfig.getParameter<bool>("producePixelhits")),
121  doPixelCut_(iConfig.getParameter<bool>("doPixelCut")),
122  produceTracks_(iConfig.getParameter<bool>("produceTracks")),
123  producePixelTracks_(iConfig.getParameter<bool>("producePixelTracks")),
124  midRapidityRange_(iConfig.getParameter<double>("midRapidityRange")),
125  trackPtCut_(iConfig.getParameter<double>("trackPtCut")),
126  trackEtaCut_(iConfig.getParameter<double>("trackEtaCut")),
127  hfEtaCut_(iConfig.getParameter<double>("hfEtaCut")),
128  reuseAny_(iConfig.getParameter<bool>("reUseCentrality")),
129  useQuality_(iConfig.getParameter<bool>("useQuality")),
130  trackQuality_(TrackBase::qualityByName(iConfig.getParameter<std::string>("trackQuality"))),
131  srcHFhits_(produceHFhits_ ? consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcHFhits"))
132  : edm::EDGetTokenT<HFRecHitCollection>()),
133  srcTowers_((produceHFtowers_ || produceETmidRap_)
134  ? consumes<CaloTowerCollection>(iConfig.getParameter<edm::InputTag>("srcTowers"))
135  : edm::EDGetTokenT<CaloTowerCollection>()),
136  srcEBhits_(produceEcalhits_ ? consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcEBhits"))
137  : edm::EDGetTokenT<EcalRecHitCollection>()),
138  srcEEhits_(produceEcalhits_ ? consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcEEhits"))
139  : edm::EDGetTokenT<EcalRecHitCollection>()),
140  srcZDChits_(produceZDChits_ ? consumes<ZDCRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcZDChits"))
141  : edm::EDGetTokenT<ZDCRecHitCollection>()),
142  srcPixelhits_(producePixelhits_
143  ? consumes<SiPixelRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcPixelhits"))
144  : edm::EDGetTokenT<SiPixelRecHitCollection>()),
145  srcTracks_(produceTracks_ ? consumes<TrackCollection>(iConfig.getParameter<edm::InputTag>("srcTracks"))
146  : edm::EDGetTokenT<TrackCollection>()),
147  srcPixelTracks_(producePixelTracks_
148  ? consumes<TrackCollection>(iConfig.getParameter<edm::InputTag>("srcPixelTracks"))
149  : edm::EDGetTokenT<TrackCollection>()),
150  srcVertex_((produceTracks_ || producePixelTracks_)
151  ? consumes<VertexCollection>(iConfig.getParameter<edm::InputTag>("srcVertex"))
152  : edm::EDGetTokenT<VertexCollection>()),
153  reuseTag_(reuseAny_ ? consumes<Centrality>(iConfig.getParameter<edm::InputTag>("srcReUse"))
154  : edm::EDGetTokenT<Centrality>()),
155  caloGeom_(produceEcalhits_ ? esConsumes<CaloGeometry, CaloGeometryRecord>()
156  : edm::ESGetToken<CaloGeometry, CaloGeometryRecord>()),
157  trackerGeom_(producePixelhits_ ? esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>()
158  : edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord>()),
159  trackerTopo_(producePixelhits_ ? esConsumes<TrackerTopology, TrackerTopologyRcd>()
160  : edm::ESGetToken<TrackerTopology, TrackerTopologyRcd>()) {
161  produces<Centrality>();
162  }
163 
164  //
165  // member functions
166  //
167 
168  // ------------ method called to produce the data ------------
169  void CentralityProducer::produce(edm::StreamID sid, edm::Event& iEvent, const edm::EventSetup& iSetup) const {
170  auto creco = std::make_unique<Centrality>();
171  Handle<Centrality> inputCentrality;
172  if (reuseAny_)
173  iEvent.getByToken(reuseTag_, inputCentrality);
174 
175  if (produceHFhits_) {
176  creco->etHFhitSumPlus_ = 0;
177  creco->etHFhitSumMinus_ = 0;
178  Handle<HFRecHitCollection> hits;
179  iEvent.getByToken(srcHFhits_, hits);
180  for (size_t ihit = 0; ihit < hits->size(); ++ihit) {
181  const HFRecHit& rechit = (*hits)[ihit];
182  if (rechit.id().ieta() > 0)
183  creco->etHFhitSumPlus_ += rechit.energy();
184  if (rechit.id().ieta() < 0)
185  creco->etHFhitSumMinus_ += rechit.energy();
186  }
187  } else {
188  if (reuseAny_) {
189  creco->etHFhitSumMinus_ = inputCentrality->EtHFhitSumMinus();
190  creco->etHFhitSumPlus_ = inputCentrality->EtHFhitSumPlus();
191  }
192  }
193 
194  if (produceHFtowers_ || produceETmidRap_) {
195  creco->etHFtowerSumPlus_ = 0;
196  creco->etHFtowerSumMinus_ = 0;
197  creco->etHFtowerSumECutPlus_ = 0;
198  creco->etHFtowerSumECutMinus_ = 0;
199  creco->etMidRapiditySum_ = 0;
200 
201  Handle<CaloTowerCollection> towers;
202  iEvent.getByToken(srcTowers_, towers);
203 
204  for (size_t i = 0; i < towers->size(); ++i) {
205  const CaloTower& tower = (*towers)[i];
206  double eta = tower.eta();
207  if (produceHFtowers_) {
208  bool isHF = tower.ietaAbs() > 29;
209  if (isHF && eta > 0) {
210  creco->etHFtowerSumPlus_ += tower.pt();
211  if (tower.energy() > 1.5)
212  creco->etHFtowerSumECutPlus_ += tower.pt();
213  if (eta > hfEtaCut_)
214  creco->etHFtruncatedPlus_ += tower.pt();
215  } else if (isHF && eta < 0) {
216  creco->etHFtowerSumMinus_ += tower.pt();
217  if (tower.energy() > 1.5)
218  creco->etHFtowerSumECutMinus_ += tower.pt();
219  if (eta < -hfEtaCut_)
220  creco->etHFtruncatedMinus_ += tower.pt();
221  }
222  } else {
223  if (reuseAny_) {
224  creco->etHFtowerSumMinus_ = inputCentrality->EtHFtowerSumMinus();
225  creco->etHFtowerSumPlus_ = inputCentrality->EtHFtowerSumPlus();
226  creco->etHFtowerSumECutMinus_ = inputCentrality->EtHFtowerSumECutMinus();
227  creco->etHFtowerSumECutPlus_ = inputCentrality->EtHFtowerSumECutPlus();
228  creco->etHFtruncatedMinus_ = inputCentrality->EtHFtruncatedMinus();
229  creco->etHFtruncatedPlus_ = inputCentrality->EtHFtruncatedPlus();
230  }
231  }
232  if (produceETmidRap_) {
233  if (std::abs(eta) < midRapidityRange_)
234  creco->etMidRapiditySum_ += tower.pt() / (midRapidityRange_ * 2.);
235  } else if (reuseAny_)
236  creco->etMidRapiditySum_ = inputCentrality->EtMidRapiditySum();
237  }
238  } else {
239  if (reuseAny_) {
240  creco->etHFtowerSumMinus_ = inputCentrality->EtHFtowerSumMinus();
241  creco->etHFtowerSumPlus_ = inputCentrality->EtHFtowerSumPlus();
242  creco->etHFtowerSumECutMinus_ = inputCentrality->EtHFtowerSumECutMinus();
243  creco->etHFtowerSumECutPlus_ = inputCentrality->EtHFtowerSumECutPlus();
244  creco->etMidRapiditySum_ = inputCentrality->EtMidRapiditySum();
245  }
246  }
247 
248  if (produceEcalhits_) {
249  creco->etEESumPlus_ = 0;
250  creco->etEESumMinus_ = 0;
251  creco->etEBSum_ = 0;
252 
253  Handle<EcalRecHitCollection> ebHits;
254  Handle<EcalRecHitCollection> eeHits;
255  iEvent.getByToken(srcEBhits_, ebHits);
256  iEvent.getByToken(srcEEhits_, eeHits);
257 
258  edm::ESHandle<CaloGeometry> cGeo = iSetup.getHandle(caloGeom_);
259  for (unsigned int i = 0; i < ebHits->size(); ++i) {
260  const EcalRecHit& hit = (*ebHits)[i];
261  const GlobalPoint& pos = cGeo->getPosition(hit.id());
262  double et = hit.energy() * (pos.perp() / pos.mag());
263  creco->etEBSum_ += et;
264  }
265 
266  for (unsigned int i = 0; i < eeHits->size(); ++i) {
267  const EcalRecHit& hit = (*eeHits)[i];
268  const GlobalPoint& pos = cGeo->getPosition(hit.id());
269  double et = hit.energy() * (pos.perp() / pos.mag());
270  if (pos.z() > 0) {
271  creco->etEESumPlus_ += et;
272  } else {
273  creco->etEESumMinus_ += et;
274  }
275  }
276  } else {
277  if (reuseAny_) {
278  creco->etEESumMinus_ = inputCentrality->EtEESumMinus();
279  creco->etEESumPlus_ = inputCentrality->EtEESumPlus();
280  creco->etEBSum_ = inputCentrality->EtEBSum();
281  }
282  }
283 
284  if (producePixelhits_) {
285  edm::ESHandle<TrackerGeometry> tGeo = iSetup.getHandle(trackerGeom_);
286  edm::ESHandle<TrackerTopology> topo = iSetup.getHandle(trackerTopo_);
287  creco->pixelMultiplicity_ = 0;
288  const SiPixelRecHitCollection* rechits;
289  Handle<SiPixelRecHitCollection> rchts;
290  iEvent.getByToken(srcPixelhits_, rchts);
291  rechits = rchts.product();
292  int nPixel = 0;
293  int nPixel_plus = 0;
294  int nPixel_minus = 0;
295  for (SiPixelRecHitCollection::const_iterator it = rechits->begin(); it != rechits->end(); it++) {
296  SiPixelRecHitCollection::DetSet hits = *it;
297  DetId detId = DetId(hits.detId());
298  SiPixelRecHitCollection::const_iterator recHitMatch = rechits->find(detId);
299  const SiPixelRecHitCollection::DetSet recHitRange = *recHitMatch;
300  for (SiPixelRecHitCollection::DetSet::const_iterator recHitIterator = recHitRange.begin();
301  recHitIterator != recHitRange.end();
302  ++recHitIterator) {
303  // add selection if needed, now all hits.
304  const SiPixelRecHit* recHit = &(*recHitIterator);
305  const PixelGeomDetUnit* pixelLayer =
306  dynamic_cast<const PixelGeomDetUnit*>(tGeo->idToDet(recHit->geographicalId()));
307  GlobalPoint gpos = pixelLayer->toGlobal(recHit->localPosition());
308  math::XYZVector rechitPos(gpos.x(), gpos.y(), gpos.z());
309  double eta = rechitPos.eta();
310  int clusterSize = recHit->cluster()->size();
311  unsigned layer = topo->layer(detId);
312  if (doPixelCut_) {
313  if (detId.det() == DetId::Tracker && detId.subdetId() == PixelSubdetector::PixelBarrel) {
314  double abeta = std::abs(eta);
315  if (layer == 1) {
316  if (18 * abeta - 40 > clusterSize)
317  continue;
318  } else if (layer == 2) {
319  if (6 * abeta - 7.2 > clusterSize)
320  continue;
321  } else if (layer == 3 || layer == 4) {
322  if (4 * abeta - 2.4 > clusterSize)
323  continue;
324  }
325  }
326  }
327  nPixel++;
328  if (eta >= 0)
329  nPixel_plus++;
330  else if (eta < 0)
331  nPixel_minus++;
332  }
333  }
334  creco->pixelMultiplicity_ = nPixel;
335  creco->pixelMultiplicityPlus_ = nPixel_plus;
336  creco->pixelMultiplicityMinus_ = nPixel_minus;
337  } else {
338  if (reuseAny_) {
339  creco->pixelMultiplicity_ = inputCentrality->multiplicityPixel();
340  creco->pixelMultiplicityPlus_ = inputCentrality->multiplicityPixelPlus();
341  creco->pixelMultiplicityMinus_ = inputCentrality->multiplicityPixelMinus();
342  }
343  }
344 
345  if (produceTracks_) {
346  double vx = -999.;
347  double vy = -999.;
348  double vz = -999.;
349  double vxError = -999.;
350  double vyError = -999.;
351  double vzError = -999.;
352  math::XYZVector vtxPos(0, 0, 0);
353 
354  Handle<VertexCollection> recoVertices;
355  iEvent.getByToken(srcVertex_, recoVertices);
356  unsigned int daughter = 0;
357  int greatestvtx = 0;
358 
359  for (unsigned int i = 0; i < recoVertices->size(); ++i) {
360  daughter = (*recoVertices)[i].tracksSize();
361  if (daughter > (*recoVertices)[greatestvtx].tracksSize())
362  greatestvtx = i;
363  }
364 
365  if (!recoVertices->empty()) {
366  vx = (*recoVertices)[greatestvtx].position().x();
367  vy = (*recoVertices)[greatestvtx].position().y();
368  vz = (*recoVertices)[greatestvtx].position().z();
369  vxError = (*recoVertices)[greatestvtx].xError();
370  vyError = (*recoVertices)[greatestvtx].yError();
371  vzError = (*recoVertices)[greatestvtx].zError();
372  }
373 
374  vtxPos = math::XYZVector(vx, vy, vz);
375 
376  Handle<TrackCollection> tracks;
377  iEvent.getByToken(srcTracks_, tracks);
378  int nTracks = 0;
379 
380  double trackCounter = 0;
381  double trackCounterEta = 0;
382  double trackCounterEtaPt = 0;
383 
384  for (unsigned int i = 0; i < tracks->size(); ++i) {
385  const Track& track = (*tracks)[i];
386  if (useQuality_ && !track.quality(trackQuality_))
387  continue;
388 
389  if (track.pt() > trackPtCut_)
390  trackCounter++;
391  if (std::abs(track.eta()) < trackEtaCut_) {
392  trackCounterEta++;
393  if (track.pt() > trackPtCut_)
394  trackCounterEtaPt++;
395  }
396 
397  math::XYZPoint v1(vx, vy, vz);
398  double dz = track.dz(v1);
399  double dzsigma2 = track.dzError() * track.dzError() + vzError * vzError;
400  double dxy = track.dxy(v1);
401  double dxysigma2 = track.dxyError() * track.dxyError() + vxError * vyError;
402 
403  const double pterrcut = 0.1;
404  const double dzrelcut = 3.0;
405  const double dxyrelcut = 3.0;
406 
407  if (track.quality(trackQuality_) && track.pt() > 0.4 && std::abs(track.eta()) < 2.4 &&
408  track.ptError() / track.pt() < pterrcut && dz * dz < dzrelcut * dzrelcut * dzsigma2 &&
409  dxy * dxy < dxyrelcut * dxyrelcut * dxysigma2) {
410  nTracks++;
411  }
412  }
413 
414  creco->trackMultiplicity_ = nTracks;
415  creco->ntracksPtCut_ = trackCounter;
416  creco->ntracksEtaCut_ = trackCounterEta;
417  creco->ntracksEtaPtCut_ = trackCounterEtaPt;
418 
419  } else {
420  if (reuseAny_) {
421  creco->trackMultiplicity_ = inputCentrality->Ntracks();
422  creco->ntracksPtCut_ = inputCentrality->NtracksPtCut();
423  creco->ntracksEtaCut_ = inputCentrality->NtracksEtaCut();
424  creco->ntracksEtaPtCut_ = inputCentrality->NtracksEtaPtCut();
425  }
426  }
427 
428  if (producePixelTracks_) {
429  Handle<TrackCollection> pixeltracks;
430  iEvent.getByToken(srcPixelTracks_, pixeltracks);
431  int nPixelTracks = pixeltracks->size();
432  int nPixelTracksPlus = 0;
433  int nPixelTracksMinus = 0;
434 
435  for (auto const& track : *pixeltracks) {
436  if (track.eta() < 0)
437  nPixelTracksMinus++;
438  else
439  nPixelTracksPlus++;
440  }
441  creco->nPixelTracks_ = nPixelTracks;
442  creco->nPixelTracksPlus_ = nPixelTracksPlus;
443  creco->nPixelTracksMinus_ = nPixelTracksMinus;
444  } else {
445  if (reuseAny_) {
446  creco->nPixelTracks_ = inputCentrality->NpixelTracks();
447  creco->nPixelTracksPlus_ = inputCentrality->NpixelTracksPlus();
448  creco->nPixelTracksMinus_ = inputCentrality->NpixelTracksMinus();
449  }
450  }
451 
452  if (produceZDChits_) {
453  creco->zdcSumPlus_ = 0;
454  creco->zdcSumMinus_ = 0;
455 
456  Handle<ZDCRecHitCollection> hits;
457  bool zdcAvailable = iEvent.getByToken(srcZDChits_, hits);
458  if (zdcAvailable) {
459  for (size_t ihit = 0; ihit < hits->size(); ++ihit) {
460  const ZDCRecHit& rechit = (*hits)[ihit];
461  if (rechit.id().zside() > 0) {
462  if (lowGainZDC_) {
463  creco->zdcSumPlus_ += rechit.lowGainEnergy();
464  } else {
465  creco->zdcSumPlus_ += rechit.energy();
466  }
467  }
468  if (rechit.id().zside() < 0) {
469  if (lowGainZDC_) {
470  creco->zdcSumMinus_ += rechit.lowGainEnergy();
471  } else {
472  creco->zdcSumMinus_ += rechit.energy();
473  }
474  }
475  }
476  } else {
477  creco->zdcSumPlus_ = -9;
478  creco->zdcSumMinus_ = -9;
479  }
480  } else {
481  if (reuseAny_) {
482  creco->zdcSumMinus_ = inputCentrality->zdcSumMinus();
483  creco->zdcSumPlus_ = inputCentrality->zdcSumPlus();
484  }
485  }
486 
487  iEvent.put(std::move(creco));
488  }
489 
490  void CentralityProducer::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
491  edm::ParameterSetDescription desc;
492  desc.add<bool>("produceHFhits", true);
493  desc.add<bool>("produceHFtowers", true);
494  desc.add<bool>("produceEcalhits", true);
495  desc.add<bool>("produceZDChits", true);
496  desc.add<bool>("produceETmidRapidity", true);
497  desc.add<bool>("producePixelhits", true);
498  desc.add<bool>("produceTracks", true);
499  desc.add<bool>("producePixelTracks", true);
500  desc.add<bool>("reUseCentrality", false);
501  desc.add<edm::InputTag>("srcHFhits", edm::InputTag("hfreco"));
502  desc.add<edm::InputTag>("srcTowers", edm::InputTag("towerMaker"));
503  desc.add<edm::InputTag>("srcEBhits", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
504  desc.add<edm::InputTag>("srcEEhits", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
505  desc.add<edm::InputTag>("srcZDChits", edm::InputTag("zdcreco"));
506  desc.add<edm::InputTag>("srcPixelhits", edm::InputTag("siPixelRecHits"));
507  desc.add<edm::InputTag>("srcTracks", edm::InputTag("hiGeneralTracks"));
508  desc.add<edm::InputTag>("srcVertex", edm::InputTag("hiSelectedVertex"));
509  desc.add<edm::InputTag>("srcReUse", edm::InputTag("hiCentrality"));
510  desc.add<edm::InputTag>("srcPixelTracks", edm::InputTag("hiPixel3PrimTracks"));
511  desc.add<bool>("doPixelCut", true);
512  desc.add<bool>("useQuality", true);
513  desc.add<string>("trackQuality", "highPurity");
514  desc.add<double>("trackEtaCut", 2);
515  desc.add<double>("trackPtCut", 1);
516  desc.add<double>("hfEtaCut", 4)->setComment("hf above the absolute value of this cut is used");
517  desc.add<double>("midRapidityRange", 1);
518  desc.add<bool>("lowGainZDC", true);
519 
520  descriptions.addDefault(desc);
521  }
522 
523  // ------------ method called once each job just before starting event loop ------------
524  void CentralityProducer::beginJob() {}
525 
526  // ------------ method called once each job just after ending the event loop ------------
527  void CentralityProducer::endJob() {}
528 
529  //define this as a plug-in
530  DEFINE_FWK_MODULE(CentralityProducer);
531 
532 } // namespace reco
TrackerGeometry::idToDet
const TrackerGeomDet * idToDet(DetId) const override
Definition: TrackerGeometry.cc:193
edm::StreamID
Definition: StreamID.h:30
HLT_FULL_cff.towers
towers
Definition: HLT_FULL_cff.py:36362
reco::CentralityProducer::srcPixelhits_
const edm::EDGetTokenT< SiPixelRecHitCollection > srcPixelhits_
Definition: CentralityProducer.cc:91
bk::beginJob
void beginJob()
Definition: Breakpoints.cc:14
EcalRecHit
Definition: EcalRecHit.h:15
electrons_cff.bool
bool
Definition: electrons_cff.py:366
mps_fire.i
i
Definition: mps_fire.py:428
HLT_FULL_cff.track
track
Definition: HLT_FULL_cff.py:11713
reco::Centrality::EtHFhitSumPlus
double EtHFhitSumPlus() const
Definition: Centrality.h:22
TrackerGeometry.h
hit::id
unsigned int id
Definition: SiStripHitEffFromCalibTree.cc:92
CaloRecHit::energy
constexpr float energy() const
Definition: CaloRecHit.h:29
PixelSubdetector::PixelBarrel
Definition: PixelSubdetector.h:11
edm::Handle::product
T const * product() const
Definition: Handle.h:70
reco::CentralityProducer::srcPixelTracks_
const edm::EDGetTokenT< TrackCollection > srcPixelTracks_
Definition: CentralityProducer.cc:93
hfClusterShapes_cfi.hits
hits
Definition: hfClusterShapes_cfi.py:5
reco::CentralityProducer::producePixelhits_
const bool producePixelhits_
Definition: CentralityProducer.cc:72
ESHandle.h
reco::CentralityProducer::caloGeom_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord > caloGeom_
Definition: CentralityProducer.cc:97
reco::CentralityProducer::doPixelCut_
const bool doPixelCut_
Definition: CentralityProducer.cc:73
CaloGeometry::getPosition
GlobalPoint getPosition(const DetId &id) const
Get the position of a given detector id.
Definition: CaloGeometry.cc:50
reco::CentralityProducer
Definition: CentralityProducer.cc:53
edm::EDGetTokenT
Definition: EDGetToken.h:33
CaloGeometryRecord
Definition: CaloGeometryRecord.h:30
edm
HLT enums.
Definition: AlignableModifier.h:19
reco::Centrality::NpixelTracksPlus
double NpixelTracksPlus() const
Definition: Centrality.h:51
TrackerTopology
Definition: TrackerTopology.h:16
DetId::det
constexpr Detector det() const
get the detector field from this detid
Definition: DetId.h:46
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
pos
Definition: PixelAliasList.h:18
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89285
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
hgcalTowerProducer_cfi.tower
tower
Definition: hgcalTowerProducer_cfi.py:4
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:197
reco::CentralityProducer::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: CentralityProducer.cc:490
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:150
HFRecHit::id
constexpr HcalDetId id() const
Definition: HFRecHit.h:26
reco::Centrality::EtEESumPlus
double EtEESumPlus() const
Definition: Centrality.h:34
edm::SortedCollection
Definition: SortedCollection.h:49
reco::CentralityProducer::beginJob
void beginJob() override
Definition: CentralityProducer.cc:524
TrackerTopology::layer
unsigned int layer(const DetId &id) const
Definition: TrackerTopology.cc:47
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:215
reco::CentralityProducer::trackerGeom_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > trackerGeom_
Definition: CentralityProducer.cc:98
reco::CentralityProducer::srcTowers_
const edm::EDGetTokenT< CaloTowerCollection > srcTowers_
Definition: CentralityProducer.cc:87
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
reco::Centrality::EtMidRapiditySum
double EtMidRapiditySum() const
Definition: Centrality.h:57
reco::Centrality::NpixelTracksMinus
double NpixelTracksMinus() const
Definition: Centrality.h:52
BeamSpotPI::nTracks
Definition: BeamSpotPayloadInspectorHelper.h:42
edm::Handle
Definition: AssociativeIterator.h:50
reco::CentralityProducer::srcEBhits_
const edm::EDGetTokenT< EcalRecHitCollection > srcEBhits_
Definition: CentralityProducer.cc:88
ESGetToken.h
rpcPointValidation_cfi.recHit
recHit
Definition: rpcPointValidation_cfi.py:7
EcalRecHitCollections.h
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23
reco::Centrality::NtracksPtCut
double NtracksPtCut() const
Definition: Centrality.h:47
reco::CentralityProducer::lowGainZDC_
const bool lowGainZDC_
Definition: CentralityProducer.cc:70
DetId
Definition: DetId.h:17
reco::CentralityProducer::srcVertex_
const edm::EDGetTokenT< VertexCollection > srcVertex_
Definition: CentralityProducer.cc:94
edmNew::DetSet::end
iterator end()
Definition: DetSetNew.h:56
MakerMacros.h
TrackerTopology.h
CaloGeometry
Definition: CaloGeometry.h:21
reco::CentralityProducer::trackEtaCut_
const double trackEtaCut_
Definition: CentralityProducer.cc:79
Track.h
HFRecHit
Definition: HFRecHit.h:11
reco::CentralityProducer::hfEtaCut_
const double hfEtaCut_
Definition: CentralityProducer.cc:80
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
PVValHelper::eta
Definition: PVValidationHelpers.h:70
reco::Centrality::EtEBSum
double EtEBSum() const
Definition: Centrality.h:39
Centrality.h
reco::Centrality::EtHFtowerSumPlus
double EtHFtowerSumPlus() const
Definition: Centrality.h:25
HcalObjRepresent::isHF
bool isHF(int etabin, int depth)
Definition: HcalObjRepresent.h:880
PixelGeomDetUnit
Definition: PixelGeomDetUnit.h:15
TrackerDigiGeometryRecord
Definition: TrackerDigiGeometryRecord.h:15
fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
reco::Centrality::EtHFtruncatedMinus
double EtHFtruncatedMinus() const
Definition: Centrality.h:32
reco::Track
Definition: Track.h:27
edm::ESHandle< CaloGeometry >
edmNew::DetSet
Definition: DetSetNew.h:22
HI_PhotonSkim_cff.rechits
rechits
Definition: HI_PhotonSkim_cff.py:76
reco::Centrality::zdcSumPlus
double zdcSumPlus() const
Definition: Centrality.h:55
reco::CentralityProducer::srcZDChits_
const edm::EDGetTokenT< ZDCRecHitCollection > srcZDChits_
Definition: CentralityProducer.cc:90
SiPixelRecHit.h
reco::CentralityProducer::trackQuality_
const reco::TrackBase::TrackQuality trackQuality_
Definition: CentralityProducer.cc:84
reco::Centrality::EtHFhitSumMinus
double EtHFhitSumMinus() const
Definition: Centrality.h:23
reco::Centrality::multiplicityPixelMinus
double multiplicityPixelMinus() const
Definition: Centrality.h:45
DetId::subdetId
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum)
Definition: DetId.h:48
reco::CentralityProducer::produce
void produce(edm::StreamID, edm::Event &, const edm::EventSetup &) const override
Definition: CentralityProducer.cc:169
edm::global::EDProducer
Definition: EDProducer.h:32
CaloGeometryRecord.h
edm::ConfigurationDescriptions
Definition: ConfigurationDescriptions.h:28
phase1PixelTopology::layer
constexpr std::array< uint8_t, layerIndexSize > layer
Definition: phase1PixelTopology.h:99
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
Vertex.h
reco::Centrality::EtHFtowerSumECutPlus
double EtHFtowerSumECutPlus() const
Definition: Centrality.h:28
HcalZDCDetId::zside
int zside() const
get the z-side of the cell (1/-1)
Definition: HcalZDCDetId.h:39
TrackerDigiGeometryRecord.h
SiPixelRecHitCollection.h
reco::CentralityProducer::produceZDChits_
const bool produceZDChits_
Definition: CentralityProducer.cc:69
HcalDetId::ieta
constexpr int ieta() const
get the cell ieta
Definition: HcalDetId.h:155
GeomDet::toGlobal
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:49
edm::ParameterSet
Definition: ParameterSet.h:47
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
DetId::Tracker
Definition: DetId.h:25
reco::CentralityProducer::producePixelTracks_
const bool producePixelTracks_
Definition: CentralityProducer.cc:75
Event.h
tracks
const uint32_t *__restrict__ const HitContainer *__restrict__ TkSoA *__restrict__ tracks
Definition: CAHitNtupletGeneratorKernelsImpl.h:159
reco::Centrality::NtracksEtaPtCut
double NtracksEtaPtCut() const
Definition: Centrality.h:49
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:31
reco::Centrality::multiplicityPixelPlus
double multiplicityPixelPlus() const
Definition: Centrality.h:44
iEvent
int iEvent
Definition: GenABIO.cc:224
CaloTower
Definition: CaloTower.h:26
reco::TrackBase
Definition: TrackBase.h:62
edm::EventSetup::getHandle
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:155
ZDCRecHit::id
HcalZDCDetId id() const
get the id
Definition: ZDCRecHit.h:18
reco::Centrality
Definition: Centrality.h:11
reco::Centrality::EtHFtruncatedPlus
double EtHFtruncatedPlus() const
Definition: Centrality.h:31
edmNew::DetSet::begin
iterator begin()
Definition: DetSetNew.h:54
EgHLTOffHistBins_cfi.et
et
Definition: EgHLTOffHistBins_cfi.py:8
edm::EventSetup
Definition: EventSetup.h:58
reco::CentralityProducer::endJob
void endJob() override
Definition: CentralityProducer.cc:527
reco::Centrality::NtracksEtaCut
double NtracksEtaCut() const
Definition: Centrality.h:48
reco::CentralityProducer::midRapidityRange_
const double midRapidityRange_
Definition: CentralityProducer.cc:77
reco::CentralityProducer::produceETmidRap_
const bool produceETmidRap_
Definition: CentralityProducer.cc:71
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >
reco::Centrality::EtEESumMinus
double EtEESumMinus() const
Definition: Centrality.h:35
InputTag.h
reco::CentralityProducer::trackPtCut_
const double trackPtCut_
Definition: CentralityProducer.cc:78
CaloTowerCollection.h
reco::CentralityProducer::srcTracks_
const edm::EDGetTokenT< TrackCollection > srcTracks_
Definition: CentralityProducer.cc:92
reco::CentralityProducer::srcEEhits_
const edm::EDGetTokenT< EcalRecHitCollection > srcEEhits_
Definition: CentralityProducer.cc:89
edmNew::DetSetVector
Definition: DetSetNew.h:13
submitPVResolutionJobs.desc
string desc
Definition: submitPVResolutionJobs.py:251
reco::Centrality::NpixelTracks
double NpixelTracks() const
Definition: Centrality.h:50
eostools.move
def move(src, dest)
Definition: eostools.py:511
std
Definition: JetResolutionObject.h:76
Ref.h
PVValHelper::dxy
Definition: PVValidationHelpers.h:48
reco::CentralityProducer::produceTracks_
const bool produceTracks_
Definition: CentralityProducer.cc:74
PVValHelper::dz
Definition: PVValidationHelpers.h:51
Frameworkfwd.h
reco::CentralityProducer::trackerTopo_
const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd > trackerTopo_
Definition: CentralityProducer.cc:99
reco::CentralityProducer::reuseAny_
const bool reuseAny_
Definition: CentralityProducer.cc:82
CaloGeometry.h
PixelGeomDetUnit.h
reco::Centrality::zdcSumMinus
double zdcSumMinus() const
Definition: Centrality.h:56
reco::Centrality::multiplicityPixel
double multiplicityPixel() const
Definition: Centrality.h:43
EventSetup.h
ZDCRecHit
Definition: ZDCRecHit.h:11
reco::CentralityProducer::useQuality_
const bool useQuality_
Definition: CentralityProducer.cc:83
HcalRecHitCollections.h
ZDCRecHit::lowGainEnergy
float lowGainEnergy() const
Definition: ZDCRecHit.h:20
math::GlobalPoint
ROOT::Math::PositionVector3D< ROOT::Math::Cartesian3D< float >, ROOT::Math::GlobalCoordinateSystemTag > GlobalPoint
point in global coordinate system
Definition: Point3D.h:18
reco::CentralityProducer::srcHFhits_
const edm::EDGetTokenT< HFRecHitCollection > srcHFhits_
Definition: CentralityProducer.cc:86
reco::Centrality::EtHFtowerSumECutMinus
double EtHFtowerSumECutMinus() const
Definition: Centrality.h:29
Candidate.h
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
TrackerTopologyRcd
Definition: TrackerTopologyRcd.h:10
ParameterSet.h
EDProducer.h
reco::CentralityProducer::produceHFtowers_
const bool produceHFtowers_
Definition: CentralityProducer.cc:67
reco::CentralityProducer::reuseTag_
const edm::EDGetTokenT< Centrality > reuseTag_
Definition: CentralityProducer.cc:95
DeDxTools::esConsumes
ESGetTokenH3DDVariant esConsumes(std::string const &Reccord, edm::ConsumesCollector &)
Definition: DeDxTools.cc:283
reco::CentralityProducer::produceEcalhits_
const bool produceEcalhits_
Definition: CentralityProducer.cc:68
edm::Event
Definition: Event.h:73
edm::ConfigurationDescriptions::addDefault
void addDefault(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:99
StreamID.h
reco::Centrality::Ntracks
double Ntracks() const
Definition: Centrality.h:46
edm::InputTag
Definition: InputTag.h:15
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:14
hit
Definition: SiStripHitEffFromCalibTree.cc:88
reco::Centrality::EtHFtowerSumMinus
double EtHFtowerSumMinus() const
Definition: Centrality.h:26
TrackerGeometry
Definition: TrackerGeometry.h:14
reco::CentralityProducer::produceHFhits_
const bool produceHFhits_
Definition: CentralityProducer.cc:66
Generated for CMSSW Reference Manual by doxygen 1.8.16