CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
EtlSimHitsValidation Class Reference

#include <Validation/MtdValidation/plugins/EtlSimHitsValidation.cc>

Inheritance diagram for EtlSimHitsValidation:
DQMEDAnalyzer edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >

Public Member Functions

 EtlSimHitsValidation (const edm::ParameterSet &)
 
 ~EtlSimHitsValidation () override
 
- Public Member Functions inherited from DQMEDAnalyzer
void accumulate (edm::Event const &event, edm::EventSetup const &setup) final
 
void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void beginRun (edm::Run const &run, edm::EventSetup const &setup) final
 
void beginStream (edm::StreamID id) final
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer ()
 
void endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
 
void endRun (edm::Run const &run, edm::EventSetup const &setup) final
 
virtual bool getCanSaveByLumi ()
 
- Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProduceroperator= (const EDProducer &)=delete
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static void globalEndJob (DQMEDAnalyzerGlobalCache const *)
 
static void globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
 
static void globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
 
static std::unique_ptr< DQMEDAnalyzerGlobalCacheinitializeGlobalCache (edm::ParameterSet const &)
 

Private Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 

Private Attributes

edm::EDGetTokenT< CrossingFrame< PSimHit > > etlSimHitsToken_
 
const std::string folder_
 
const float hitMinEnergy2Dis_
 
MonitorElementmeHitEnergy_ [4]
 
MonitorElementmeHitEta_ [4]
 
MonitorElementmeHitEvsEta_ [4]
 
MonitorElementmeHitEvsPhi_ [4]
 
MonitorElementmeHitPhi_ [4]
 
MonitorElementmeHitTime_ [4]
 
MonitorElementmeHitTvsE_ [4]
 
MonitorElementmeHitTvsEta_ [4]
 
MonitorElementmeHitTvsPhi_ [4]
 
MonitorElementmeHitX_ [4]
 
MonitorElementmeHitXlocal_ [4]
 
MonitorElementmeHitY_ [4]
 
MonitorElementmeHitYlocal_ [4]
 
MonitorElementmeHitZ_ [4]
 
MonitorElementmeHitZlocal_ [4]
 
MonitorElementmeNhits_ [4]
 
MonitorElementmeNtrkPerCell_ [4]
 
MonitorElementmeOccupancy_ [4]
 
edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecordmtdgeoToken_
 
edm::ESGetToken< MTDTopology, MTDTopologyRcdmtdtopoToken_
 

Additional Inherited Members

- Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore DQMStore
 
typedef dqm::reco::MonitorElement MonitorElement
 
- Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using CacheTypes = CacheContexts< T... >
 
using GlobalCache = typename CacheTypes::GlobalCache
 
using HasAbility = AbilityChecker< T... >
 
using InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
 
using LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
 
using LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
 
using LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
 
using RunCache = typename CacheTypes::RunCache
 
using RunContext = RunContextT< RunCache, GlobalCache >
 
using RunSummaryCache = typename CacheTypes::RunSummaryCache
 
- Protected Member Functions inherited from DQMEDAnalyzer
uint64_t meId () const
 
- Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMTokenlumiToken_
 
edm::EDPutTokenT< DQMTokenrunToken_
 
unsigned int streamId_
 

Detailed Description

Description: ETL SIM hits validation

Implementation: [Notes on implementation]

Definition at line 39 of file EtlSimHitsValidation.cc.

Constructor & Destructor Documentation

◆ EtlSimHitsValidation()

EtlSimHitsValidation::EtlSimHitsValidation ( const edm::ParameterSet iConfig)
explicit

Definition at line 89 of file EtlSimHitsValidation.cc.

References etlSimHitsToken_, edm::ParameterSet::getParameter(), mtdgeoToken_, and mtdtopoToken_.

90  : folder_(iConfig.getParameter<std::string>("folder")),
91  hitMinEnergy2Dis_(iConfig.getParameter<double>("hitMinimumEnergy2Dis")) {
92  etlSimHitsToken_ = consumes<CrossingFrame<PSimHit> >(iConfig.getParameter<edm::InputTag>("inputTag"));
93  mtdgeoToken_ = esConsumes<MTDGeometry, MTDDigiGeometryRecord>();
94  mtdtopoToken_ = esConsumes<MTDTopology, MTDTopologyRcd>();
95 }
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecord > mtdgeoToken_
edm::EDGetTokenT< CrossingFrame< PSimHit > > etlSimHitsToken_
const std::string folder_
edm::ESGetToken< MTDTopology, MTDTopologyRcd > mtdtopoToken_

◆ ~EtlSimHitsValidation()

EtlSimHitsValidation::~EtlSimHitsValidation ( )
override

Definition at line 97 of file EtlSimHitsValidation.cc.

97 {}

Member Function Documentation

◆ analyze()

void EtlSimHitsValidation::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
)
overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 100 of file EtlSimHitsValidation.cc.

References angle_units::operators::convertMmToCm(), geant_units::operators::convertUnitsTo(), TauDecayModes::dec, hcalRecHitTable_cff::detId, hcalRecHitTable_cff::energy, etlSimHitsToken_, Exception, dqm::impl::MonitorElement::Fill(), relativeConstraints::geom, edm::EventSetup::getTransientHandle(), hitMinEnergy2Dis_, l1ctLayer2EG_cff::id, iEvent, LogDebug, edm::makeValid(), meHitEnergy_, meHitEta_, meHitEvsEta_, meHitEvsPhi_, meHitPhi_, meHitTime_, meHitTvsE_, meHitTvsEta_, meHitTvsPhi_, meHitX_, meHitXlocal_, meHitY_, meHitYlocal_, meHitZ_, meHitZlocal_, meNhits_, meNtrkPerCell_, meOccupancy_, mtdgeoToken_, DetId::rawId(), nano_mu_digi_cff::rawId, rpcPointValidation_cfi::simHit, findQualityFiles::size, hcalRecHitTable_cff::time, GeomDet::toGlobal(), mps_merge::weight, x, y, z, and ecaldqm::zside().

100  {
101  using namespace edm;
102  using namespace geant_units::operators;
103 
104  auto geometryHandle = iSetup.getTransientHandle(mtdgeoToken_);
105  const MTDGeometry* geom = geometryHandle.product();
106 
107  auto etlSimHitsHandle = makeValid(iEvent.getHandle(etlSimHitsToken_));
108  MixCollection<PSimHit> etlSimHits(etlSimHitsHandle.product());
109 
110  std::unordered_map<uint32_t, MTDHit> m_etlHits[4];
111  std::unordered_map<uint32_t, std::set<int> > m_etlTrkPerCell[4];
112 
113  // --- Loop over the ETL SIM hits
114 
115  int idet = 999;
116 
117  size_t index(0);
118 
119  for (auto const& simHit : etlSimHits) {
120  index++;
121  LogDebug("EtlSimHitsValidation") << "SimHit # " << index << " detId " << simHit.detUnitId() << " ene "
122  << simHit.energyLoss() << " tof " << simHit.tof() << " tId " << simHit.trackId();
123 
124  // --- Use only hits compatible with the in-time bunch-crossing
125  if (simHit.tof() < 0 || simHit.tof() > 25.)
126  continue;
127 
128  ETLDetId id = simHit.detUnitId();
129  if ((id.zside() == -1) && (id.nDisc() == 1)) {
130  idet = 0;
131  } else if ((id.zside() == -1) && (id.nDisc() == 2)) {
132  idet = 1;
133  } else if ((id.zside() == 1) && (id.nDisc() == 1)) {
134  idet = 2;
135  } else if ((id.zside() == 1) && (id.nDisc() == 2)) {
136  idet = 3;
137  } else {
138  edm::LogWarning("EtlSimHitsValidation") << "Unknown ETL DetId configuration: " << id;
139  continue;
140  }
141 
142  m_etlTrkPerCell[idet][id.rawId()].insert(simHit.trackId());
143 
144  auto simHitIt = m_etlHits[idet].emplace(id.rawId(), MTDHit()).first;
145 
146  // --- Accumulate the energy (in MeV) of SIM hits in the same detector cell
147  (simHitIt->second).energy += convertUnitsTo(0.001_MeV, simHit.energyLoss());
148 
149  // --- Get the time of the first SIM hit in the cell
150  if ((simHitIt->second).time == 0 || simHit.tof() < (simHitIt->second).time) {
151  (simHitIt->second).time = simHit.tof();
152 
153  auto hit_pos = simHit.localPosition();
154  (simHitIt->second).x = hit_pos.x();
155  (simHitIt->second).y = hit_pos.y();
156  (simHitIt->second).z = hit_pos.z();
157  }
158  LogDebug("EtlSimHitsValidation") << "Registered in idet " << idet;
159 
160  } // simHit loop
161 
162  // ==============================================================================
163  // Histogram filling
164  // ==============================================================================
165 
166  for (int idet = 0; idet < 4; ++idet) { //two disks per side
167  meNhits_[idet]->Fill(m_etlHits[idet].size());
168  LogDebug("EtlSimHitsValidation") << "idet " << idet << " #hits " << m_etlHits[idet].size();
169 
170  for (auto const& hit : m_etlTrkPerCell[idet]) {
171  meNtrkPerCell_[idet]->Fill((hit.second).size());
172  }
173 
174  for (auto const& hit : m_etlHits[idet]) {
175  double weight = 1.0;
176  if ((hit.second).energy < hitMinEnergy2Dis_)
177  continue;
178  // --- Get the SIM hit global position
179  ETLDetId detId(hit.first);
180  DetId geoId = detId.geographicalId();
181  const MTDGeomDet* thedet = geom->idToDet(geoId);
182  if (thedet == nullptr)
183  throw cms::Exception("EtlSimHitsValidation") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
184  << detId.rawId() << ") is invalid!" << std::dec << std::endl;
185 
186  Local3DPoint local_point(
187  convertMmToCm((hit.second).x), convertMmToCm((hit.second).y), convertMmToCm((hit.second).z));
188  const auto& global_point = thedet->toGlobal(local_point);
189 
190  if (detId.discSide() == 1) {
191  weight = -weight;
192  }
193 
194  // --- Fill the histograms
195 
196  meHitEnergy_[idet]->Fill((hit.second).energy);
197  meHitTime_[idet]->Fill((hit.second).time);
198  meHitXlocal_[idet]->Fill((hit.second).x);
199  meHitYlocal_[idet]->Fill((hit.second).y);
200  meHitZlocal_[idet]->Fill((hit.second).z);
201  meOccupancy_[idet]->Fill(global_point.x(), global_point.y(), weight);
202  meHitX_[idet]->Fill(global_point.x());
203  meHitY_[idet]->Fill(global_point.y());
204  meHitZ_[idet]->Fill(global_point.z());
205  meHitPhi_[idet]->Fill(global_point.phi());
206  meHitEta_[idet]->Fill(global_point.eta());
207  meHitTvsE_[idet]->Fill((hit.second).energy, (hit.second).time);
208  meHitEvsPhi_[idet]->Fill(global_point.phi(), (hit.second).energy);
209  meHitEvsEta_[idet]->Fill(global_point.eta(), (hit.second).energy);
210  meHitTvsPhi_[idet]->Fill(global_point.phi(), (hit.second).time);
211  meHitTvsEta_[idet]->Fill(global_point.eta(), (hit.second).time);
212 
213  } // hit loop
214 
215  } // idet loop
216 }
size
Write out results.
edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecord > mtdgeoToken_
MonitorElement * meNhits_[4]
MonitorElement * meHitYlocal_[4]
edm::EDGetTokenT< CrossingFrame< PSimHit > > etlSimHitsToken_
Definition: MTDHit.h:4
Definition: weight.py:1
MonitorElement * meOccupancy_[4]
int zside(DetId const &)
MonitorElement * meHitY_[4]
MonitorElement * meHitTvsEta_[4]
constexpr NumType convertUnitsTo(double desiredUnits, NumType val)
Definition: GeantUnits.h:73
void Fill(long long x)
MonitorElement * meHitXlocal_[4]
int iEvent
Definition: GenABIO.cc:224
MonitorElement * meHitTvsPhi_[4]
MonitorElement * meHitX_[4]
MonitorElement * meHitEvsEta_[4]
MonitorElement * meHitEta_[4]
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:49
MonitorElement * meHitTime_[4]
Definition: DetId.h:17
MonitorElement * meHitPhi_[4]
MonitorElement * meHitTvsE_[4]
constexpr NumType convertMmToCm(NumType millimeters)
Definition: angle_units.h:44
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
Detector identifier class for the Endcap Timing Layer.
Definition: ETLDetId.h:16
HLT enums.
MonitorElement * meHitZlocal_[4]
ESTransientHandle< T > getTransientHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:141
Log< level::Warning, false > LogWarning
MonitorElement * meNtrkPerCell_[4]
MonitorElement * meHitEnergy_[4]
auto makeValid(const U &iOtherHandleType) noexcept(false)
Definition: ValidHandle.h:52
MonitorElement * meHitEvsPhi_[4]
#define LogDebug(id)
MonitorElement * meHitZ_[4]

◆ bookHistograms()

void EtlSimHitsValidation::bookHistograms ( DQMStore::IBooker ibook,
edm::Run const &  run,
edm::EventSetup const &  iSetup 
)
overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 219 of file EtlSimHitsValidation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), folder_, meHitEnergy_, meHitEta_, meHitEvsEta_, meHitEvsPhi_, meHitPhi_, meHitTime_, meHitTvsE_, meHitTvsEta_, meHitTvsPhi_, meHitX_, meHitXlocal_, meHitY_, meHitYlocal_, meHitZ_, meHitZlocal_, meNhits_, meNtrkPerCell_, meOccupancy_, and dqm::implementation::NavigatorBase::setCurrentFolder().

221  {
222  ibook.setCurrentFolder(folder_);
223 
224  // --- histograms booking
225 
226  meNhits_[0] = ibook.book1D("EtlNhitsZnegD1",
227  "Number of ETL cells with SIM hits (-Z, Single(topo1D)/First(topo2D) disk);N_{ETL cells}",
228  100,
229  0.,
230  5000.);
231  meNhits_[1] = ibook.book1D(
232  "EtlNhitsZnegD2", "Number of ETL cells with SIM hits (-Z, Second disk);N_{ETL cells}", 100, 0., 5000.);
233  meNhits_[2] = ibook.book1D("EtlNhitsZposD1",
234  "Number of ETL cells with SIM hits (+Z, Single(topo1D)/First(topo2D) disk);N_{ETL cells}",
235  100,
236  0.,
237  5000.);
238  meNhits_[3] = ibook.book1D(
239  "EtlNhitsZposD2", "Number of ETL cells with SIM hits (+Z, Second Disk);N_{ETL cells}", 100, 0., 5000.);
240  meNtrkPerCell_[0] = ibook.book1D("EtlNtrkPerCellZnegD1",
241  "Number of tracks per ETL sensor (-Z, Single(topo1D)/First(topo2D) disk);N_{trk}",
242  10,
243  0.,
244  10.);
245  meNtrkPerCell_[1] =
246  ibook.book1D("EtlNtrkPerCellZnegD2", "Number of tracks per ETL sensor (-Z, Second disk);N_{trk}", 10, 0., 10.);
247  meNtrkPerCell_[2] = ibook.book1D("EtlNtrkPerCellZposD1",
248  "Number of tracks per ETL sensor (+Z, Single(topo1D)/First(topo2D) disk);N_{trk}",
249  10,
250  0.,
251  10.);
252  meNtrkPerCell_[3] =
253  ibook.book1D("EtlNtrkPerCellZposD2", "Number of tracks per ETL sensor (+Z, Second disk);N_{trk}", 10, 0., 10.);
254  meHitEnergy_[0] = ibook.book1D(
255  "EtlHitEnergyZnegD1", "ETL SIM hits energy (-Z, Single(topo1D)/First(topo2D) disk);E_{SIM} [MeV]", 100, 0., 3.);
256  meHitEnergy_[1] =
257  ibook.book1D("EtlHitEnergyZnegD2", "ETL SIM hits energy (-Z, Second disk);E_{SIM} [MeV]", 100, 0., 3.);
258  meHitEnergy_[2] = ibook.book1D(
259  "EtlHitEnergyZposD1", "ETL SIM hits energy (+Z, Single(topo1D)/First(topo2D) disk);E_{SIM} [MeV]", 100, 0., 3.);
260  meHitEnergy_[3] =
261  ibook.book1D("EtlHitEnergyZposD2", "ETL SIM hits energy (+Z, Second disk);E_{SIM} [MeV]", 100, 0., 3.);
262  meHitTime_[0] = ibook.book1D(
263  "EtlHitTimeZnegD1", "ETL SIM hits ToA (-Z, Single(topo1D)/First(topo2D) disk);ToA_{SIM} [ns]", 100, 0., 25.);
264  meHitTime_[1] = ibook.book1D("EtlHitTimeZnegD2", "ETL SIM hits ToA (-Z, Second disk);ToA_{SIM} [ns]", 100, 0., 25.);
265  meHitTime_[2] = ibook.book1D(
266  "EtlHitTimeZposD1", "ETL SIM hits ToA (+Z, Single(topo1D)/First(topo2D) disk);ToA_{SIM} [ns]", 100, 0., 25.);
267  meHitTime_[3] = ibook.book1D("EtlHitTimeZposD2", "ETL SIM hits ToA (+Z, Second disk);ToA_{SIM} [ns]", 100, 0., 25.);
268 
269  meHitXlocal_[0] = ibook.book1D("EtlHitXlocalZnegD1",
270  "ETL SIM local X (-Z, Single(topo1D)/First(topo2D) disk);X_{SIM}^{LOC} [mm]",
271  100,
272  -25.,
273  25.);
274  meHitXlocal_[1] =
275  ibook.book1D("EtlHitXlocalZnegD2", "ETL SIM local X (-Z, Second disk);X_{SIM}^{LOC} [mm]", 100, -25., 25.);
276  meHitXlocal_[2] = ibook.book1D("EtlHitXlocalZposD1",
277  "ETL SIM local X (+Z, Single(topo1D)/First(topo2D) disk);X_{SIM}^{LOC} [mm]",
278  100,
279  -25.,
280  25.);
281  meHitXlocal_[3] =
282  ibook.book1D("EtlHitXlocalZposD2", "ETL SIM local X (+Z, Second disk);X_{SIM}^{LOC} [mm]", 100, -25., 25.);
283 
284  meHitYlocal_[0] = ibook.book1D("EtlHitYlocalZnegD1",
285  "ETL SIM local Y (-Z, Single(topo1D)/First(topo2D) disk);Y_{SIM}^{LOC} [mm]",
286  100,
287  -48.,
288  48.);
289  meHitYlocal_[1] =
290  ibook.book1D("EtlHitYlocalZnegD2", "ETL SIM local Y (-Z, Second Disk);Y_{SIM}^{LOC} [mm]", 100, -48., 48.);
291  meHitYlocal_[2] = ibook.book1D("EtlHitYlocalZposD1",
292  "ETL SIM local Y (+Z, Single(topo1D)/First(topo2D) disk);Y_{SIM}^{LOC} [mm]",
293  100,
294  -48.,
295  48.);
296  meHitYlocal_[3] =
297  ibook.book1D("EtlHitYlocalZposD2", "ETL SIM local Y (+Z, Second disk);Y_{SIM}^{LOC} [mm]", 100, -48., 48.);
298  meHitZlocal_[0] = ibook.book1D("EtlHitZlocalZnegD1",
299  "ETL SIM local Z (-Z, Single(topo1D)/First(topo2D) disk);Z_{SIM}^{LOC} [mm]",
300  80,
301  -0.16,
302  0.16);
303  meHitZlocal_[1] =
304  ibook.book1D("EtlHitZlocalZnegD2", "ETL SIM local Z (-Z, Second disk);Z_{SIM}^{LOC} [mm]", 80, -0.16, 0.16);
305  meHitZlocal_[2] = ibook.book1D("EtlHitZlocalZposD1",
306  "ETL SIM local Z (+Z, Single(topo1D)/First(topo2D) disk);Z_{SIM}^{LOC} [mm]",
307  80,
308  -0.16,
309  0.16);
310  meHitZlocal_[3] =
311  ibook.book1D("EtlHitZlocalZposD2", "ETL SIM local Z (+Z, Second disk);Z_{SIM}^{LOC} [mm]", 80, -0.16, 0.16);
312 
313  meOccupancy_[0] =
314  ibook.book2D("EtlOccupancyZnegD1",
315  "ETL SIM hits occupancy (-Z, Single(topo1D)/First(topo2D) disk);X_{SIM} [cm];Y_{SIM} [cm]",
316  135,
317  -135.,
318  135.,
319  135,
320  -135.,
321  135.);
322  meOccupancy_[1] = ibook.book2D("EtlOccupancyZnegD2",
323  "ETL SIM hits occupancy (-Z, Second disk);X_{SIM} [cm];Y_{SIM} [cm]",
324  135,
325  -135.,
326  135.,
327  135,
328  -135.,
329  135.);
330  meOccupancy_[2] =
331  ibook.book2D("EtlOccupancyZposD1",
332  "ETL SIM hits occupancy (+Z, Single(topo1D)/First(topo2D) disk);X_{SIM} [cm];Y_{SIM} [cm]",
333  135,
334  -135.,
335  135.,
336  135,
337  -135.,
338  135.);
339  meOccupancy_[3] = ibook.book2D("EtlOccupancyZposD2",
340  "ETL SIM hits occupancy (+Z, Second disk);X_{SIM} [cm];Y_{SIM} [cm]",
341  135,
342  -135.,
343  135.,
344  135,
345  -135.,
346  135.);
347 
348  meHitX_[0] = ibook.book1D(
349  "EtlHitXZnegD1", "ETL SIM hits X (+Z, Single(topo1D)/First(topo2D) disk);X_{SIM} [cm]", 100, -130., 130.);
350  meHitX_[1] = ibook.book1D("EtlHitXZnegD2", "ETL SIM hits X (-Z, Second disk);X_{SIM} [cm]", 100, -130., 130.);
351  meHitX_[2] = ibook.book1D(
352  "EtlHitXZposD1", "ETL SIM hits X (+Z, Single(topo1D)/First(topo2D) disk);X_{SIM} [cm]", 100, -130., 130.);
353  meHitX_[3] = ibook.book1D("EtlHitXZposD2", "ETL SIM hits X (+Z, Second disk);X_{SIM} [cm]", 100, -130., 130.);
354  meHitY_[0] = ibook.book1D(
355  "EtlHitYZnegD1", "ETL SIM hits Y (-Z, Single(topo1D)/First(topo2D) disk);Y_{SIM} [cm]", 100, -130., 130.);
356  meHitY_[1] = ibook.book1D("EtlHitYZnegD2", "ETL SIM hits Y (-Z, Second disk);Y_{SIM} [cm]", 100, -130., 130.);
357  meHitY_[2] = ibook.book1D(
358  "EtlHitYZposD1", "ETL SIM hits Y (+Z, Single(topo1D)/First(topo2D) disk);Y_{SIM} [cm]", 100, -130., 130.);
359  meHitY_[3] = ibook.book1D("EtlHitYZposD2", "ETL SIM hits Y (+Z, Second disk);Y_{SIM} [cm]", 100, -130., 130.);
360  meHitZ_[0] = ibook.book1D(
361  "EtlHitZZnegD1", "ETL SIM hits Z (-Z, Single(topo1D)/First(topo2D) disk);Z_{SIM} [cm]", 100, -302., -298.);
362  meHitZ_[1] = ibook.book1D("EtlHitZZnegD2", "ETL SIM hits Z (-Z, Second disk);Z_{SIM} [cm]", 100, -304., -300.);
363  meHitZ_[2] = ibook.book1D(
364  "EtlHitZZposD1", "ETL SIM hits Z (+Z, Single(topo1D)/First(topo2D) disk);Z_{SIM} [cm]", 100, 298., 302.);
365  meHitZ_[3] = ibook.book1D("EtlHitZZposD2", "ETL SIM hits Z (+Z, Second disk);Z_{SIM} [cm]", 100, 300., 304.);
366 
367  meHitPhi_[0] = ibook.book1D(
368  "EtlHitPhiZnegD1", "ETL SIM hits #phi (-Z, Single(topo1D)/First(topo2D) disk);#phi_{SIM} [rad]", 100, -3.15, 3.15);
369  meHitPhi_[1] =
370  ibook.book1D("EtlHitPhiZnegD2", "ETL SIM hits #phi (-Z, Second disk);#phi_{SIM} [rad]", 100, -3.15, 3.15);
371  meHitPhi_[2] = ibook.book1D(
372  "EtlHitPhiZposD1", "ETL SIM hits #phi (+Z, Single(topo1D)/First(topo2D) disk);#phi_{SIM} [rad]", 100, -3.15, 3.15);
373  meHitPhi_[3] =
374  ibook.book1D("EtlHitPhiZposD2", "ETL SIM hits #phi (+Z, Second disk);#phi_{SIM} [rad]", 100, -3.15, 3.15);
375  meHitEta_[0] = ibook.book1D(
376  "EtlHitEtaZnegD1", "ETL SIM hits #eta (-Z, Single(topo1D)/First(topo2D) disk);#eta_{SIM}", 100, -3.2, -1.56);
377  meHitEta_[1] = ibook.book1D("EtlHitEtaZnegD2", "ETL SIM hits #eta (-Z, Second disk);#eta_{SIM}", 100, -3.2, -1.56);
378  meHitEta_[2] = ibook.book1D(
379  "EtlHitEtaZposD1", "ETL SIM hits #eta (+Z, Single(topo1D)/First(topo2D) disk);#eta_{SIM}", 100, 1.56, 3.2);
380  meHitEta_[3] = ibook.book1D("EtlHitEtaZposD2", "ETL SIM hits #eta (+Z, Second disk);#eta_{SIM}", 100, 1.56, 3.2);
381 
382  meHitTvsE_[0] =
383  ibook.bookProfile("EtlHitTvsEZnegD1",
384  "ETL SIM time vs energy (-Z, Single(topo1D)/First(topo2D) disk);E_{SIM} [MeV];T_{SIM} [ns]",
385  50,
386  0.,
387  2.,
388  0.,
389  100.);
390  meHitTvsE_[1] = ibook.bookProfile(
391  "EtlHitTvsEZnegD2", "ETL SIM time vs energy (-Z, Second disk);E_{SIM} [MeV];T_{SIM} [ns]", 50, 0., 2., 0., 100.);
392  meHitTvsE_[2] =
393  ibook.bookProfile("EtlHitTvsEZposD1",
394  "ETL SIM time vs energy (+Z, Single(topo1D)/First(topo2D) disk);E_{SIM} [MeV];T_{SIM} [ns]",
395  50,
396  0.,
397  2.,
398  0.,
399  100.);
400  meHitTvsE_[3] = ibook.bookProfile(
401  "EtlHitTvsEZposD2", "ETL SIM time vs energy (+Z, Second disk);E_{SIM} [MeV];T_{SIM} [ns]", 50, 0., 2., 0., 100.);
402  meHitEvsPhi_[0] =
403  ibook.bookProfile("EtlHitEvsPhiZnegD1",
404  "ETL SIM energy vs #phi (-Z, Single(topo1D)/First(topo2D) disk);#phi_{SIM} [rad];E_{SIM} [MeV]",
405  50,
406  -3.15,
407  3.15,
408  0.,
409  100.);
410  meHitEvsPhi_[1] = ibook.bookProfile("EtlHitEvsPhiZnegD2",
411  "ETL SIM energy vs #phi (-Z, Second disk);#phi_{SIM} [rad];E_{SIM} [MeV]",
412  50,
413  -3.15,
414  3.15,
415  0.,
416  100.);
417  meHitEvsPhi_[2] =
418  ibook.bookProfile("EtlHitEvsPhiZposD1",
419  "ETL SIM energy vs #phi (+Z, Single(topo1D)/First(topo2D) disk);#phi_{SIM} [rad];E_{SIM} [MeV]",
420  50,
421  -3.15,
422  3.15,
423  0.,
424  100.);
425  meHitEvsPhi_[3] = ibook.bookProfile("EtlHitEvsPhiZposD2",
426  "ETL SIM energy vs #phi (+Z, Second disk);#phi_{SIM} [rad];E_{SIM} [MeV]",
427  50,
428  -3.15,
429  3.15,
430  0.,
431  100.);
432  meHitEvsEta_[0] =
433  ibook.bookProfile("EtlHitEvsEtaZnegD1",
434  "ETL SIM energy vs #eta (-Z, Single(topo1D)/First(topo2D) disk);#eta_{SIM};E_{SIM} [MeV]",
435  50,
436  -3.2,
437  -1.56,
438  0.,
439  100.);
440  meHitEvsEta_[1] = ibook.bookProfile("EtlHitEvsEtaZnegD2",
441  "ETL SIM energy vs #eta (-Z, Second disk);#eta_{SIM};E_{SIM} [MeV]",
442  50,
443  -3.2,
444  -1.56,
445  0.,
446  100.);
447  meHitEvsEta_[2] =
448  ibook.bookProfile("EtlHitEvsEtaZposD1",
449  "ETL SIM energy vs #eta (+Z, Single(topo1D)/First(topo2D) disk);#eta_{SIM};E_{SIM} [MeV]",
450  50,
451  1.56,
452  3.2,
453  0.,
454  100.);
455  meHitEvsEta_[3] = ibook.bookProfile("EtlHitEvsEtaZposD2",
456  "ETL SIM energy vs #eta (+Z, Second disk);#eta_{SIM};E_{SIM} [MeV]",
457  50,
458  1.56,
459  3.2,
460  0.,
461  100.);
462  meHitTvsPhi_[0] =
463  ibook.bookProfile("EtlHitTvsPhiZnegD1",
464  "ETL SIM time vs #phi (-Z, Single(topo1D)/First(topo2D) disk);#phi_{SIM} [rad];T_{SIM} [ns]",
465  50,
466  -3.15,
467  3.15,
468  0.,
469  100.);
470  meHitTvsPhi_[1] = ibook.bookProfile("EtlHitTvsPhiZnegD2",
471  "ETL SIM time vs #phi (-Z, Second disk);#phi_{SIM} [rad];T_{SIM} [ns]",
472  50,
473  -3.15,
474  3.15,
475  0.,
476  100.);
477  meHitTvsPhi_[2] =
478  ibook.bookProfile("EtlHitTvsPhiZposD1",
479  "ETL SIM time vs #phi (+Z, Single(topo1D)/First(topo2D) disk);#phi_{SIM} [rad];T_{SIM} [ns]",
480  50,
481  -3.15,
482  3.15,
483  0.,
484  100.);
485  meHitTvsPhi_[3] = ibook.bookProfile("EtlHitTvsPhiZposD2",
486  "ETL SIM time vs #phi (+Z, Second disk);#phi_{SIM} [rad];T_{SIM} [ns]",
487  50,
488  -3.15,
489  3.15,
490  0.,
491  100.);
492  meHitTvsEta_[0] =
493  ibook.bookProfile("EtlHitTvsEtaZnegD1",
494  "ETL SIM time vs #eta (-Z, Single(topo1D)/First(topo2D) disk);#eta_{SIM};T_{SIM} [ns]",
495  50,
496  -3.2,
497  -1.56,
498  0.,
499  100.);
500  meHitTvsEta_[1] = ibook.bookProfile(
501  "EtlHitTvsEtaZnegD2", "ETL SIM time vs #eta (-Z, Second disk);#eta_{SIM};T_{SIM} [ns]", 50, -3.2, -1.56, 0., 100.);
502  meHitTvsEta_[2] =
503  ibook.bookProfile("EtlHitTvsEtaZposD1",
504  "ETL SIM time vs #eta (+Z, Single(topo1D)/First(topo2D) disk);#eta_{SIM};T_{SIM} [ns]",
505  50,
506  1.56,
507  3.2,
508  0.,
509  100.);
510  meHitTvsEta_[3] = ibook.bookProfile(
511  "EtlHitTvsEtaZposD2", "ETL SIM time vs #eta (+Z, Second disk);#eta_{SIM};T_{SIM} [ns]", 50, 1.56, 3.2, 0., 100.);
512 }
MonitorElement * meNhits_[4]
MonitorElement * meHitYlocal_[4]
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:36
MonitorElement * meOccupancy_[4]
MonitorElement * meHitY_[4]
MonitorElement * meHitTvsEta_[4]
MonitorElement * meHitXlocal_[4]
MonitorElement * bookProfile(TString const &name, TString const &title, int nchX, double lowX, double highX, int, double lowY, double highY, char const *option="s", FUNC onbooking=NOOP())
Definition: DQMStore.h:408
MonitorElement * meHitTvsPhi_[4]
MonitorElement * meHitX_[4]
MonitorElement * meHitEvsEta_[4]
MonitorElement * meHitEta_[4]
MonitorElement * meHitTime_[4]
MonitorElement * meHitPhi_[4]
MonitorElement * meHitTvsE_[4]
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:221
MonitorElement * meHitZlocal_[4]
MonitorElement * meNtrkPerCell_[4]
MonitorElement * meHitEnergy_[4]
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
const std::string folder_
MonitorElement * meHitEvsPhi_[4]
MonitorElement * meHitZ_[4]

◆ fillDescriptions()

void EtlSimHitsValidation::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 515 of file EtlSimHitsValidation.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

515  {
517 
518  desc.add<std::string>("folder", "MTD/ETL/SimHits");
519  desc.add<edm::InputTag>("inputTag", edm::InputTag("mix", "g4SimHitsFastTimerHitsEndcap"));
520  desc.add<double>("hitMinimumEnergy2Dis", 0.001); // [MeV]
521 
522  descriptions.add("etlSimHitsValid", desc);
523 }
void add(std::string const &label, ParameterSetDescription const &psetDescription)

Member Data Documentation

◆ etlSimHitsToken_

edm::EDGetTokenT<CrossingFrame<PSimHit> > EtlSimHitsValidation::etlSimHitsToken_
private

Definition at line 56 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and EtlSimHitsValidation().

◆ folder_

const std::string EtlSimHitsValidation::folder_
private

Definition at line 53 of file EtlSimHitsValidation.cc.

Referenced by bookHistograms().

◆ hitMinEnergy2Dis_

const float EtlSimHitsValidation::hitMinEnergy2Dis_
private

Definition at line 54 of file EtlSimHitsValidation.cc.

Referenced by analyze().

◆ meHitEnergy_

MonitorElement* EtlSimHitsValidation::meHitEnergy_[4]
private

Definition at line 66 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitEta_

MonitorElement* EtlSimHitsValidation::meHitEta_[4]
private

Definition at line 79 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitEvsEta_

MonitorElement* EtlSimHitsValidation::meHitEvsEta_[4]
private

Definition at line 83 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitEvsPhi_

MonitorElement* EtlSimHitsValidation::meHitEvsPhi_[4]
private

Definition at line 82 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitPhi_

MonitorElement* EtlSimHitsValidation::meHitPhi_[4]
private

Definition at line 78 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitTime_

MonitorElement* EtlSimHitsValidation::meHitTime_[4]
private

Definition at line 67 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitTvsE_

MonitorElement* EtlSimHitsValidation::meHitTvsE_[4]
private

Definition at line 81 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitTvsEta_

MonitorElement* EtlSimHitsValidation::meHitTvsEta_[4]
private

Definition at line 85 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitTvsPhi_

MonitorElement* EtlSimHitsValidation::meHitTvsPhi_[4]
private

Definition at line 84 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitX_

MonitorElement* EtlSimHitsValidation::meHitX_[4]
private

Definition at line 75 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitXlocal_

MonitorElement* EtlSimHitsValidation::meHitXlocal_[4]
private

Definition at line 69 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitY_

MonitorElement* EtlSimHitsValidation::meHitY_[4]
private

Definition at line 76 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitYlocal_

MonitorElement* EtlSimHitsValidation::meHitYlocal_[4]
private

Definition at line 70 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitZ_

MonitorElement* EtlSimHitsValidation::meHitZ_[4]
private

Definition at line 77 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meHitZlocal_

MonitorElement* EtlSimHitsValidation::meHitZlocal_[4]
private

Definition at line 71 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meNhits_

MonitorElement* EtlSimHitsValidation::meNhits_[4]
private

Definition at line 63 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meNtrkPerCell_

MonitorElement* EtlSimHitsValidation::meNtrkPerCell_[4]
private

Definition at line 64 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ meOccupancy_

MonitorElement* EtlSimHitsValidation::meOccupancy_[4]
private

Definition at line 73 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and bookHistograms().

◆ mtdgeoToken_

edm::ESGetToken<MTDGeometry, MTDDigiGeometryRecord> EtlSimHitsValidation::mtdgeoToken_
private

Definition at line 58 of file EtlSimHitsValidation.cc.

Referenced by analyze(), and EtlSimHitsValidation().

◆ mtdtopoToken_

edm::ESGetToken<MTDTopology, MTDTopologyRcd> EtlSimHitsValidation::mtdtopoToken_
private

Definition at line 59 of file EtlSimHitsValidation.cc.

Referenced by EtlSimHitsValidation().