CMS 3D CMS Logo
SimG4HGCalValidation.cc
Go to the documentation of this file.
1 // File: SimG4HGCalValidation.cc
3 // Description: Main analysis class for HGCal Validation of G4 Hits
5 
6 #include "SimG4Core/Notification/interface/BeginOfJob.h"
7 #include "SimG4Core/Notification/interface/BeginOfRun.h"
8 #include "SimG4Core/Notification/interface/BeginOfEvent.h"
9 #include "SimG4Core/Notification/interface/EndOfEvent.h"
10 #include "SimG4Core/Watcher/interface/SimProducer.h"
11 #include "SimG4Core/Notification/interface/Observer.h"
12 
13 // to retreive hits
14 #include "DataFormats/DetId/interface/DetId.h"
15 #include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
16 #include "DataFormats/ForwardDetId/interface/ForwardSubdetector.h"
17 #include "DataFormats/Math/interface/Point3D.h"
18 #include "SimDataFormats/ValidationFormats/interface/PHGCalValidInfo.h"
19 #include "SimDataFormats/CaloTest/interface/HGCalTestNumbering.h"
20 #include "DataFormats/HcalDetId/interface/HcalTestNumbering.h"
21 #include "SimG4CMS/Calo/interface/HGCNumberingScheme.h"
22 #include "SimG4CMS/Calo/interface/HGCalNumberingScheme.h"
23 
24 #include "FWCore/Framework/interface/Event.h"
25 #include "FWCore/Framework/interface/EventSetup.h"
26 #include "FWCore/Framework/interface/ESHandle.h"
27 #include "FWCore/ParameterSet/interface/ParameterSet.h"
28 #include "FWCore/MessageLogger/interface/MessageLogger.h"
29 
30 #include "Geometry/Records/interface/IdealGeometryRecord.h"
31 #include "Geometry/Records/interface/HcalSimNumberingRecord.h"
32 #include "Geometry/HcalCommonData/interface/HcalNumberingFromDDD.h"
33 #include "Geometry/HcalCommonData/interface/HcalDDDSimConstants.h"
34 #include "Geometry/HGCalCommonData/interface/HGCalGeometryMode.h"
35 #include "Geometry/HGCalCommonData/interface/HGCalDDDConstants.h"
36 
37 #include "SimG4Core/Watcher/interface/SimWatcherFactory.h"
38 #include "FWCore/PluginManager/interface/ModuleDef.h"
39 #include "FWCore/Framework/interface/MakerMacros.h"
40 
41 #include "G4Step.hh"
42 #include "G4Track.hh"
43 #include "G4NavigationHistory.hh"
44 #include "G4TouchableHistory.hh"
45 
46 #include "CLHEP/Units/GlobalSystemOfUnits.h"
47 #include "CLHEP/Units/GlobalPhysicalConstants.h"
48 
49 #include <cmath>
50 #include <iostream>
51 #include <iomanip>
52 #include <memory>
53 #include <vector>
54 #include <string>
55 
56 class SimG4HGCalValidation : public SimProducer,
57  public Observer<const BeginOfJob*>,
58  public Observer<const BeginOfEvent*>,
59  public Observer<const G4Step*> {
60 public:
61  SimG4HGCalValidation(const edm::ParameterSet& p);
62  ~SimG4HGCalValidation() override;
63 
64  void produce(edm::Event&, const edm::EventSetup&) override;
65 
66 private:
67  SimG4HGCalValidation(const SimG4HGCalValidation&); // stop default
68  const SimG4HGCalValidation& operator=(const SimG4HGCalValidation&);
69 
70  void init();
71 
72  // observer classes
73  void update(const BeginOfJob* job) override;
74  void update(const BeginOfEvent* evt) override;
75  void update(const G4Step* step) override;
76 
77  // analysis related class
78  void layerAnalysis(PHGCalValidInfo&);
79  void clear();
80 
81 private:
82  //Keep reference to instantiate HcalNumberingFromDDD later
83  HcalNumberingFromDDD* numberingFromDDD_;
84 
85  //HGCal numbering scheme
86  std::vector<HGCNumberingScheme*> hgcNumbering_;
87  std::vector<HGCalNumberingScheme*> hgcalNumbering_;
88 
89  // to read from parameter set
90  std::vector<std::string> names_;
91  std::vector<int> types_, detTypes_, subdet_;
92  std::string labelLayer_;
93 
94  // parameters from geometry
95  int levelT1_, levelT2_;
96 
97  // some private members for ananlysis
98  unsigned int count_;
99  int verbosity_;
100  double edepEE_, edepHEF_, edepHEB_;
101  std::vector<double> hgcEEedep_, hgcHEFedep_, hgcHEBedep_;
102  std::vector<unsigned int> dets_, hgchitDets_, hgchitIndex_;
103  std::vector<double> hgchitX_, hgchitY_, hgchitZ_;
104 };
105 
106 SimG4HGCalValidation::SimG4HGCalValidation(const edm::ParameterSet& p)
107  : numberingFromDDD_(nullptr), levelT1_(999), levelT2_(999), count_(0) {
108  edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("SimG4HGCalValidation");
109  names_ = m_Anal.getParameter<std::vector<std::string> >("Names");
110  types_ = m_Anal.getParameter<std::vector<int> >("Types");
111  detTypes_ = m_Anal.getParameter<std::vector<int> >("DetTypes");
112  labelLayer_ = m_Anal.getParameter<std::string>("LabelLayerInfo");
113  verbosity_ = m_Anal.getUntrackedParameter<int>("Verbosity", 0);
114 
115  produces<PHGCalValidInfo>(labelLayer_);
116 
117  if (verbosity_ > 0) {
118  edm::LogVerbatim("ValidHGCal") << "HGCalTestAnalysis:: Initialised as "
119  << "observer of begin events and of G4step "
120  << "with Parameter values: \n\tLabel : " << labelLayer_ << " and with "
121  << names_.size() << " detectors";
122  for (unsigned int k = 0; k < names_.size(); ++k)
123  edm::LogVerbatim("ValidHGCal") << " [" << k << "] " << names_[k] << " Type " << types_[k] << " DetType "
124  << detTypes_[k];
125  }
126 }
127 
128 SimG4HGCalValidation::~SimG4HGCalValidation() {
129  delete numberingFromDDD_;
130  for (auto number : hgcNumbering_)
131  delete number;
132  for (auto number : hgcalNumbering_)
133  delete number;
134 }
135 
136 void SimG4HGCalValidation::produce(edm::Event& e, const edm::EventSetup&) {
137  std::unique_ptr<PHGCalValidInfo> productLayer(new PHGCalValidInfo);
138  layerAnalysis(*productLayer);
139  e.put(std::move(productLayer), labelLayer_);
140 }
141 
142 void SimG4HGCalValidation::update(const BeginOfJob* job) {
143  const edm::EventSetup* es = (*job)();
144  for (unsigned int type = 0; type < types_.size(); ++type) {
145  int layers(0);
146  int detType = detTypes_[type];
147  G4String nameX = "HGCal";
148  if (types_[type] <= 1) {
149  if (types_[type] == 0) {
150  subdet_.emplace_back(ForwardSubdetector::ForwardEmpty);
151  if (detType == 0)
152  dets_.emplace_back((int)(DetId::HGCalEE));
153  else if (detType == 1)
154  dets_.emplace_back((int)(DetId::HGCalHSi));
155  else
156  dets_.emplace_back((int)(DetId::HGCalHSc));
157  } else {
158  dets_.push_back((unsigned int)(DetId::Forward));
159  if (detType == 0)
160  subdet_.emplace_back((int)(ForwardSubdetector::HGCEE));
161  else if (detType == 1)
162  subdet_.emplace_back((int)(ForwardSubdetector::HGCHEF));
163  else
164  subdet_.emplace_back((int)(ForwardSubdetector::HGCHEB));
165  }
166  if (detType == 0)
167  nameX = "HGCalEESensitive";
168  else if (detType == 1)
169  nameX = "HGCalHESiliconSensitive";
170  else
171  nameX = "HGCalHEScintillatorSensitive";
172  edm::ESHandle<HGCalDDDConstants> hdc;
173  es->get<IdealGeometryRecord>().get(nameX, hdc);
174  if (hdc.isValid()) {
175  levelT1_ = hdc->levelTop(0);
176  levelT2_ = hdc->levelTop(1);
177  if (hdc->tileTrapezoid()) {
178  types_[type] = -1;
179  hgcalNumbering_.emplace_back(new HGCalNumberingScheme(*hdc, (DetId::Detector)(dets_[type]), nameX));
180  } else if (hdc->waferHexagon6()) {
181  types_[type] = 1;
182  hgcNumbering_.push_back(new HGCNumberingScheme(*hdc, nameX));
183  } else {
184  types_[type] = 0;
185  hgcalNumbering_.emplace_back(new HGCalNumberingScheme(*hdc, (DetId::Detector)(dets_[type]), nameX));
186  }
187  layers = hdc->layers(false);
188  } else {
189  edm::LogError("ValidHGCal") << "Cannot find HGCalDDDConstants for " << nameX;
190  throw cms::Exception("Unknown", "ValidHGCal") << "Cannot find HGCalDDDConstants for " << nameX << "\n";
191  }
192  } else {
193  nameX = "HcalEndcap";
194  dets_.push_back((unsigned int)(DetId::Hcal));
195  subdet_.push_back((int)(HcalSubdetector::HcalEndcap));
196  edm::ESHandle<HcalDDDSimConstants> hdc;
197  es->get<HcalSimNumberingRecord>().get(hdc);
198  if (hdc.isValid()) {
199  numberingFromDDD_ = new HcalNumberingFromDDD(hdc.product());
200  layers = 18;
201  } else {
202  edm::LogError("ValidHGCal") << "Cannot find HcalDDDSimConstant";
203  throw cms::Exception("Unknown", "ValidHGCal") << "Cannot find HcalDDDSimConstant\n";
204  }
205  }
206  if (detType == 0) {
207  for (int i = 0; i < layers; ++i)
208  hgcEEedep_.push_back(0);
209  } else if (detType == 1) {
210  for (int i = 0; i < layers; ++i)
211  hgcHEFedep_.push_back(0);
212  } else {
213  for (int i = 0; i < layers; ++i)
214  hgcHEBedep_.push_back(0);
215  }
216  if (verbosity_ > 0)
217  edm::LogVerbatim("ValidHGCal") << "[" << type << "]: " << nameX << " det " << dets_[type] << " subdet "
218  << subdet_[type] << " with " << layers << " layers";
219  }
220 }
221 
222 //=================================================================== per EVENT
223 void SimG4HGCalValidation::update(const BeginOfEvent* evt) {
224  int iev = (*evt)()->GetEventID();
225  if (verbosity_ > 0)
226  edm::LogVerbatim("ValidHGCal") << "SimG4HGCalValidation: =====> Begin "
227  << "event = " << iev;
228 
229  ++count_;
230  edepEE_ = edepHEF_ = edepHEB_ = 0.;
231 
232  //HGCal variables
233  for (unsigned int i = 0; i < hgcEEedep_.size(); i++)
234  hgcEEedep_[i] = 0.;
235  for (unsigned int i = 0; i < hgcHEFedep_.size(); i++)
236  hgcHEFedep_[i] = 0.;
237  for (unsigned int i = 0; i < hgcHEBedep_.size(); i++)
238  hgcHEBedep_[i] = 0.;
239 
240  //Cache reset
241  clear();
242 }
243 
244 //=================================================================== each STEP
245 void SimG4HGCalValidation::update(const G4Step* aStep) {
246  if (aStep != nullptr) {
247  G4VPhysicalVolume* curPV = aStep->GetPreStepPoint()->GetPhysicalVolume();
248  const G4String& name = curPV->GetName();
249  G4VSensitiveDetector* curSD = aStep->GetPreStepPoint()->GetSensitiveDetector();
250  if (verbosity_ > 1)
251  edm::LogVerbatim("ValidHGCal") << "ValidHGCal::Step in " << name << " at "
252  << aStep->GetPreStepPoint()->GetPosition();
253 
254  // Only for Sensitive detector
255  if (curSD != nullptr) {
256  int type(-1);
257  for (unsigned int k = 0; k < names_.size(); ++k) {
258  if (name.find(names_[k].c_str()) != std::string::npos) {
259  type = k;
260  break;
261  }
262  }
263  int detType = (type >= 0) ? detTypes_[type] : -1;
264  if (verbosity_ > 0)
265  edm::LogVerbatim("ValidHGCal") << "ValidHGCal::In SD " << name << " type " << type << ":" << detType;
266 
267  // Right type of SD
268  if (type >= 0) {
269  //Get the 32-bit index of the hit cell
270  const G4TouchableHistory* touchable =
271  static_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
272  unsigned int index(0);
273  int layer(0);
274  G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
275  if (types_[type] <= 1) {
276  // HGCal
277  G4ThreeVector localpos = touchable->GetHistory()->GetTopTransform().TransformPoint(hitPoint);
278  float globalZ = touchable->GetTranslation(0).z();
279  int iz(globalZ > 0 ? 1 : -1);
280  int module(-1), cell(-1);
281  if (types_[type] == 1) {
282  if (touchable->GetHistoryDepth() == levelT1_) {
283  layer = touchable->GetReplicaNumber(0);
284  } else {
285  layer = touchable->GetReplicaNumber(2);
286  module = touchable->GetReplicaNumber(1);
287  cell = touchable->GetReplicaNumber(0);
288  }
289  index =
290  hgcNumbering_[type]->getUnitID((ForwardSubdetector)(subdet_[type]), layer, module, cell, iz, localpos);
291  } else {
292  if ((touchable->GetHistoryDepth() == levelT1_) || (touchable->GetHistoryDepth() == levelT2_)) {
293  layer = touchable->GetReplicaNumber(0);
294  } else {
295  layer = touchable->GetReplicaNumber(3);
296  module = touchable->GetReplicaNumber(2);
297  cell = touchable->GetReplicaNumber(1);
298  }
299  double weight(0);
300  index = hgcalNumbering_[type]->getUnitID(layer, module, cell, iz, hitPoint, weight);
301  }
302  if (verbosity_ > 1)
303  edm::LogVerbatim("ValidHGCal")
304  << "HGCal: " << name << " Layer " << layer << " Module " << module << " Cell " << cell;
305  } else {
306  // Hcal
307  int depth = (touchable->GetReplicaNumber(0)) % 10 + 1;
308  int lay = (touchable->GetReplicaNumber(0) / 10) % 100 + 1;
309  int det = (touchable->GetReplicaNumber(1)) / 1000;
310  HcalNumberingFromDDD::HcalID tmp =
311  numberingFromDDD_->unitID(det, math::XYZVectorD(hitPoint.x(), hitPoint.y(), hitPoint.z()), depth, lay);
312  index = HcalTestNumbering::packHcalIndex(tmp.subdet, tmp.zside, tmp.depth, tmp.etaR, tmp.phis, tmp.lay);
313  layer = tmp.lay;
314  if (verbosity_ > 1)
315  edm::LogVerbatim("ValidHGCal")
316  << "HCAL: " << det << ":" << depth << ":" << lay << " o/p " << tmp.subdet << ":" << tmp.zside << ":"
317  << tmp.depth << ":" << tmp.etaR << ":" << tmp.phis << ":" << tmp.lay << " point " << hitPoint << " "
318  << hitPoint.rho() << ":" << hitPoint.eta() << ":" << hitPoint.phi();
319  }
320 
321  double edeposit = aStep->GetTotalEnergyDeposit();
322  if (verbosity_ > 0)
323  edm::LogVerbatim("ValidHGCal") << "Layer " << layer << " Index " << std::hex << index << std::dec << " Edep "
324  << edeposit << " hit at " << hitPoint;
325  if (detType == 0) {
326  edepEE_ += edeposit;
327  if (layer < (int)(hgcEEedep_.size()))
328  hgcEEedep_[layer] += edeposit;
329  } else if (detType == 1) {
330  edepHEF_ += edeposit;
331  if (layer < (int)(hgcHEFedep_.size()))
332  hgcHEFedep_[layer] += edeposit;
333  } else {
334  edepHEB_ += edeposit;
335  if (layer < (int)(hgcHEBedep_.size()))
336  hgcHEBedep_[layer] += edeposit;
337  }
338  G4String nextVolume("XXX");
339  if (aStep->GetTrack()->GetNextVolume() != nullptr)
340  nextVolume = aStep->GetTrack()->GetNextVolume()->GetName();
341 
342  if (nextVolume.c_str() != name.c_str()) { //save hit when it exits cell
343  if (std::find(hgchitIndex_.begin(), hgchitIndex_.end(), index) == hgchitIndex_.end()) {
344  hgchitDets_.push_back(dets_[type]);
345  hgchitIndex_.push_back(index);
346  hgchitX_.push_back(hitPoint.x());
347  hgchitY_.push_back(hitPoint.y());
348  hgchitZ_.push_back(hitPoint.z());
349  }
350  }
351  } // it is right type of SD
352  } // it is in a SD
353  } //if aStep!=NULL
354 } //end update aStep
355 
356 //================================================================ End of EVENT
357 
358 void SimG4HGCalValidation::layerAnalysis(PHGCalValidInfo& product) {
359  if (verbosity_ > 0)
360  edm::LogVerbatim("ValidHGCal") << "\n ===>>> SimG4HGCalValidation: Energy "
361  << "deposit\n at EE : " << std::setw(6) << edepEE_ / CLHEP::MeV
362  << "\n at HEF: " << std::setw(6) << edepHEF_ / CLHEP::MeV
363  << "\n at HEB: " << std::setw(6) << edepHEB_ / CLHEP::MeV;
364 
365  //Fill HGC Variables
366  product.fillhgcHits(hgchitDets_, hgchitIndex_, hgchitX_, hgchitY_, hgchitZ_);
367  product.fillhgcLayers(edepEE_, edepHEF_, edepHEB_, hgcEEedep_, hgcHEFedep_, hgcHEBedep_);
368 }
369 
370 //---------------------------------------------------
371 void SimG4HGCalValidation::clear() {
372  hgchitDets_.erase(hgchitDets_.begin(), hgchitDets_.end());
373  hgchitIndex_.erase(hgchitIndex_.begin(), hgchitIndex_.end());
374  hgchitX_.erase(hgchitX_.begin(), hgchitX_.end());
375  hgchitY_.erase(hgchitY_.begin(), hgchitY_.end());
376  hgchitZ_.erase(hgchitZ_.begin(), hgchitZ_.end());
377 }
378 
379 DEFINE_SIMWATCHER(SimG4HGCalValidation);
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
SimG4HGCalValidation::levelT2_
int levelT2_
Definition: SimG4HGCalValidation.cc:95
SimG4HGCalValidation::hgcalNumbering_
std::vector< HGCalNumberingScheme * > hgcalNumbering_
Definition: SimG4HGCalValidation.cc:87
mps_fire.i
i
Definition: mps_fire.py:428
PHGCalValidInfo
Definition: PHGCalValidInfo.h:12
Observer
Definition: Observer.h:23
MessageLogger.h
SimG4HGCalValidation::edepEE_
double edepEE_
Definition: SimG4HGCalValidation.cc:100
HcalNumberingFromDDD::HcalID
Definition: HcalNumberingFromDDD.h:21
ForwardEmpty
Definition: ForwardSubdetector.h:5
ESHandle.h
HcalSimNumberingRecord.h
BeginOfJob.h
step
step
Definition: StallMonitor.cc:94
HGCNumberingScheme.h
ForwardSubdetector
ForwardSubdetector
Definition: ForwardSubdetector.h:4
SimProducer.h
mps_merge.weight
weight
Definition: mps_merge.py:88
HGCalGeometryMode.h
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
HcalNumberingFromDDD
Definition: HcalNumberingFromDDD.h:16
DetId::Hcal
Definition: DetId.h:28
MeV
const double MeV
EndOfEvent.h
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
SimG4HGCalValidation::dets_
std::vector< unsigned int > dets_
Definition: SimG4HGCalValidation.cc:102
HcalTestNumbering.h
SimG4HGCalValidation::hgchitIndex_
std::vector< unsigned int > hgchitIndex_
Definition: SimG4HGCalValidation.cc:102
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
createJobs.tmp
tmp
align.sh
Definition: createJobs.py:716
HcalNumberingFromDDD.h
SimG4HGCalValidation::names_
std::vector< std::string > names_
Definition: SimG4HGCalValidation.cc:90
SimG4HGCalValidation::init
void init()
ForwardSubdetector.h
Observer.h
SimG4HGCalValidation
Definition: SimG4HGCalValidation.cc:56
SimG4HGCalValidation::hgcNumbering_
std::vector< HGCNumberingScheme * > hgcNumbering_
Definition: SimG4HGCalValidation.cc:86
SimG4HGCalValidation::subdet_
std::vector< int > subdet_
Definition: SimG4HGCalValidation.cc:91
DEFINE_SIMWATCHER
#define DEFINE_SIMWATCHER(type)
Definition: SimWatcherFactory.h:13
SimG4HGCalValidation::hgchitY_
std::vector< double > hgchitY_
Definition: SimG4HGCalValidation.cc:103
SimG4HGCalValidation::hgchitX_
std::vector< double > hgchitX_
Definition: SimG4HGCalValidation.cc:103
DetId::HGCalHSi
Definition: DetId.h:33
DetId::HGCalEE
Definition: DetId.h:32
MakerMacros.h
HGCalNumberingScheme.h
HGCalDDDConstants::waferHexagon6
bool waferHexagon6() const
Definition: HGCalDDDConstants.h:126
BeginOfRun.h
edm::EventSetup::get
T get() const
Definition: EventSetup.h:80
SimG4HGCalValidation::labelLayer_
std::string labelLayer_
Definition: SimG4HGCalValidation.cc:92
contentValuesFiles.number
number
Definition: contentValuesFiles.py:53
SimG4HGCalValidation::produce
void produce(edm::Event &, const edm::EventSetup &) override
Definition: SimG4HGCalValidation.cc:136
SimG4HGCalValidation::types_
std::vector< int > types_
Definition: SimG4HGCalValidation.cc:91
SimG4HGCalValidation::hgcHEFedep_
std::vector< double > hgcHEFedep_
Definition: SimG4HGCalValidation.cc:101
edm::ESHandle
Definition: DTSurvey.h:22
SimG4HGCalValidation::operator=
const SimG4HGCalValidation & operator=(const SimG4HGCalValidation &)
SimG4HGCalValidation::hgcEEedep_
std::vector< double > hgcEEedep_
Definition: SimG4HGCalValidation.cc:101
BeginOfJob
Definition: BeginOfJob.h:8
dqmdumpme.k
k
Definition: dqmdumpme.py:60
SimG4HGCalValidation::edepHEB_
double edepHEB_
Definition: SimG4HGCalValidation.cc:100
SimG4HGCalValidation::hgchitZ_
std::vector< double > hgchitZ_
Definition: SimG4HGCalValidation.cc:103
LEDCalibrationChannels.depth
depth
Definition: LEDCalibrationChannels.py:65
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HGCalDDDConstants::tileTrapezoid
bool tileTrapezoid() const
Definition: HGCalDDDConstants.h:110
HGCEE
Definition: ForwardSubdetector.h:8
math::XYZVectorD
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > XYZVectorD
spatial vector with cartesian internal representation
Definition: Vector3D.h:8
SimG4HGCalValidation::~SimG4HGCalValidation
~SimG4HGCalValidation() override
Definition: SimG4HGCalValidation.cc:128
edm::ParameterSet
Definition: ParameterSet.h:47
HcalTestNumbering::packHcalIndex
static uint32_t packHcalIndex(int det, int z, int depth, int eta, int phi, int lay)
Definition: HcalTestNumbering.cc:7
SimG4HGCalValidation::hgchitDets_
std::vector< unsigned int > hgchitDets_
Definition: SimG4HGCalValidation.cc:102
Event.h
SimProducer
Definition: SimProducer.h:64
type
type
Definition: SiPixelVCal_PayloadInspector.cc:37
SimWatcherFactory.h
gainCalibHelper::gainCalibPI::type
type
Definition: SiPixelGainCalibHelper.h:39
BeginOfEvent.h
HGCalDDDConstants::layers
unsigned int layers(bool reco) const
Definition: HGCalDDDConstants.cc:561
ModuleDef.h
BeginOfEvent
Definition: BeginOfEvent.h:6
IdealGeometryRecord.h
DetId::Detector
Detector
Definition: DetId.h:24
edm::EventSetup
Definition: EventSetup.h:57
HcalSubdetector.h
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
get
#define get
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:44
HGCNumberingScheme
Definition: HGCNumberingScheme.h:13
PHGCalValidInfo::fillhgcLayers
void fillhgcLayers(const double edepEE, const double edepHEF, const double edepHEB, const std::vector< double > &eedep, const std::vector< double > &hefdep, const std::vector< double > &hebdep)
Definition: PHGCalValidInfo.cc:21
SimG4HGCalValidation::edepHEF_
double edepHEF_
Definition: SimG4HGCalValidation.cc:100
eostools.move
def move(src, dest)
Definition: eostools.py:511
SimG4HGCalValidation::SimG4HGCalValidation
SimG4HGCalValidation(const edm::ParameterSet &p)
Definition: SimG4HGCalValidation.cc:106
HcalEndcap
Definition: HcalAssistant.h:34
DetId.h
SimG4HGCalValidation::detTypes_
std::vector< int > detTypes_
Definition: SimG4HGCalValidation.cc:91
SimG4HGCalValidation::count_
unsigned int count_
Definition: SimG4HGCalValidation.cc:98
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
DetId::HGCalHSc
Definition: DetId.h:34
SimG4HGCalValidation::clear
void clear()
Definition: SimG4HGCalValidation.cc:371
Exception
Definition: hltDiff.cc:246
HcalNumberingFromDDD::unitID
HcalID unitID(int det, const math::XYZVectorD &pos, int depth, int lay=-1) const
Definition: HcalNumberingFromDDD.cc:30
Point3D.h
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
EventSetup.h
PHGCalValidInfo::fillhgcHits
void fillhgcHits(const std::vector< unsigned int > &hitdets, const std::vector< unsigned int > &hitindex, const std::vector< double > &hitvtxX, const std::vector< double > &hitvtxY, const std::vector< double > &hitvtxZ)
Definition: PHGCalValidInfo.cc:7
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
AlignmentPI::index
index
Definition: AlignmentPayloadInspectorHelper.h:46
HGCalDDDConstants::levelTop
int levelTop(int ind=0) const
Definition: HGCalDDDConstants.h:87
HGCHEF
Definition: ForwardSubdetector.h:9
HGCalDDDConstants.h
ParameterSet.h
SimG4HGCalValidation::hgcHEBedep_
std::vector< double > hgcHEBedep_
Definition: SimG4HGCalValidation.cc:101
SimG4HGCalValidation::verbosity_
int verbosity_
Definition: SimG4HGCalValidation.cc:99
DetId::Forward
Definition: DetId.h:30
SimG4HGCalValidation::numberingFromDDD_
HcalNumberingFromDDD * numberingFromDDD_
Definition: SimG4HGCalValidation.cc:83
HGCalNumberingScheme
Definition: HGCalNumberingScheme.h:16
edm::Event
Definition: Event.h:73
SimG4HGCalValidation::update
void update(const BeginOfJob *job) override
This routine will be called when the appropriate signal arrives.
Definition: SimG4HGCalValidation.cc:142
PHGCalValidInfo.h
edm::Log
Definition: MessageLogger.h:70
HcalSimNumberingRecord
Definition: HcalSimNumberingRecord.h:25
TauDecayModes.dec
dec
Definition: TauDecayModes.py:143
hgcalTopologyTester_cfi.layers
layers
Definition: hgcalTopologyTester_cfi.py:8
SimG4HGCalValidation::levelT1_
int levelT1_
Definition: SimG4HGCalValidation.cc:95
weight
Definition: weight.py:1
IdealGeometryRecord
Definition: IdealGeometryRecord.h:25
SimG4HGCalValidation::layerAnalysis
void layerAnalysis(PHGCalValidInfo &)
Definition: SimG4HGCalValidation.cc:358
HGCHEB
Definition: ForwardSubdetector.h:10
HcalDDDSimConstants.h
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
HGCalTestNumbering.h
Generated for CMSSW Reference Manual by doxygen 1.8.16