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

Inheritance diagram for BtlSimHitsValidation:

Public Member Functions
	BtlSimHitsValidation (const edm::ParameterSet &)

	~BtlSimHitsValidation () 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 EDProducer &	operator= (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< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (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 > >	btlSimHitsToken_

const std::string	folder_

const float	hitMinEnergy_

MonitorElement *	meHitDeltaE12vsE1_

MonitorElement *	meHitDeltaT2_

MonitorElement *	meHitDeltaT3_

MonitorElement *	meHitDeltaT4_

MonitorElement *	meHitE1overEcellBulk1_

MonitorElement *	meHitE1overEcellBulk2_

MonitorElement *	meHitE1overEcellBulk3_

MonitorElement *	meHitE1overEcellBulk4_

MonitorElement *	meHitE1overEcellTail1_

MonitorElement *	meHitE1overEcellTail2_

MonitorElement *	meHitE1overEcellTail3_

MonitorElement *	meHitE1overEcellTail4_

MonitorElement *	meHitEnergy_

MonitorElement *	meHitEta_

MonitorElement *	meHitEvsEta_

MonitorElement *	meHitEvsPhi_

MonitorElement *	meHitEvsZ_

MonitorElement *	meHitLogEnergy_

MonitorElement *	meHitPhi_

MonitorElement *	meHitTime_

MonitorElement *	meHitTrkID1_

MonitorElement *	meHitTrkID2_

MonitorElement *	meHitTrkID3_

MonitorElement *	meHitTrkID4_

MonitorElement *	meHitTvsE_

MonitorElement *	meHitTvsEta_

MonitorElement *	meHitTvsPhi_

MonitorElement *	meHitTvsZ_

MonitorElement *	meHitX_

MonitorElement *	meHitXlocal_

MonitorElement *	meHitY_

MonitorElement *	meHitYlocal_

MonitorElement *	meHitZ_

MonitorElement *	meHitZlocal_

MonitorElement *	meNevents_

MonitorElement *	meNhits_

MonitorElement *	meNtrkPerCell_

MonitorElement *	meOccupancy_

edm::ESGetToken< MTDGeometry, MTDDigiGeometryRecord >	mtdgeoToken_

edm::ESGetToken< MTDTopology, MTDTopologyRcd >	mtdtopoToken_

Static Private Attributes
static constexpr float	BXTime_ = 25.

static constexpr float	cellEneCut_ = 10.

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< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Description: BTL SIM hits validation

Implementation: [Notes on implementation]

Definition at line 44 of file BtlSimHitsValidation.cc.

Constructor & Destructor Documentation

◆ BtlSimHitsValidation()

BtlSimHitsValidation::BtlSimHitsValidation ( const edm::ParameterSet & iConfig )

explicit

Definition at line 123 of file BtlSimHitsValidation.cc.

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

     : folder_(iConfig.getParameter<std::string>("folder")),
       hitMinEnergy_(iConfig.getParameter<double>("hitMinimumEnergy")) {
   btlSimHitsToken_ = consumes<CrossingFrame<PSimHit> >(iConfig.getParameter<edm::InputTag>("inputTag"));
   mtdgeoToken_ = esConsumes<MTDGeometry, MTDDigiGeometryRecord>();
   mtdtopoToken_ = esConsumes<MTDTopology, MTDTopologyRcd>();
 }

◆ ~BtlSimHitsValidation()

BtlSimHitsValidation::~BtlSimHitsValidation ( )

override

Definition at line 131 of file BtlSimHitsValidation.cc.

131 {}

Member Function Documentation

◆ analyze()

void BtlSimHitsValidation::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 134 of file BtlSimHitsValidation.cc.

References btlSimHitsToken_, BXTime_, cellEneCut_, angle_units::operators::convertGeVToMeV(), angle_units::operators::convertMmToCm(), MTDTopologyMode::crysLayoutFromTopoMode(), TauDecayModes::dec, hcalRecHitTable_cff::detId, Exception, dqm::impl::MonitorElement::Fill(), relativeConstraints::geom, edm::EventSetup::getTransientHandle(), hitMinEnergy_, iEvent, createfilelist::int, edm::makeValid(), meHitDeltaE12vsE1_, meHitDeltaT2_, meHitDeltaT3_, meHitDeltaT4_, meHitE1overEcellBulk1_, meHitE1overEcellBulk2_, meHitE1overEcellBulk3_, meHitE1overEcellBulk4_, meHitE1overEcellTail1_, meHitE1overEcellTail2_, meHitE1overEcellTail3_, meHitE1overEcellTail4_, meHitEnergy_, meHitEta_, meHitEvsEta_, meHitEvsPhi_, meHitEvsZ_, meHitLogEnergy_, meHitPhi_, meHitTime_, meHitTrkID1_, meHitTrkID2_, meHitTrkID3_, meHitTrkID4_, meHitTvsE_, meHitTvsEta_, meHitTvsPhi_, meHitTvsZ_, meHitX_, meHitXlocal_, meHitY_, meHitYlocal_, meHitZ_, meHitZlocal_, meNevents_, meNhits_, meNtrkPerCell_, meOccupancy_, mtdgeoToken_, mtdtopoToken_, RectangularMTDTopology::nrows(), RectangularMTDTopology::pixelToModuleLocalPoint(), DetId::rawId(), nano_mu_digi_cff::rawId, rpcPointValidation_cfi::simHit, ProxyMTDTopology::specificTopology(), edm::swap(), GeomDet::toGlobal(), GeomDet::topology(), HLT_2024v14_cff::topology, hit::x, hit::y, and hit::z.

                                                                                       {
   using namespace edm;
   using namespace geant_units::operators;
 
   auto geometryHandle = iSetup.getTransientHandle(mtdgeoToken_);
   const MTDGeometry* geom = geometryHandle.product();
 
   auto topologyHandle = iSetup.getTransientHandle(mtdtopoToken_);
   const MTDTopology* topology = topologyHandle.product();
 
   auto btlSimHitsHandle = makeValid(iEvent.getHandle(btlSimHitsToken_));
   MixCollection<PSimHit> btlSimHits(btlSimHitsHandle.product());
 
   std::unordered_map<uint32_t, std::unordered_map<uint64_t, MTDHit> > m_btlHitsPerCell;
 
   // --- Loop over the BLT SIM hits and accumulate the hits with the same track ID in each cell
   for (auto const& simHit : btlSimHits) {
     // --- Use only SIM hits compatible with the in-time bunch-crossing
     if (simHit.tof() < 0. || simHit.tof() > BXTime_) {
       continue;
     }
 
     DetId id = simHit.detUnitId();
 
     // --- Build a global track ID by combining the SIM hit's eventId and trackId
     uint64_t globalTrkID = ((uint64_t)simHit.eventId().rawId() << 32) | simHit.trackId();
 
     // --- Sum the energies of SIM hits with the same track ID in the same cell
     m_btlHitsPerCell[id.rawId()][globalTrkID].energy += convertGeVToMeV(simHit.energyLoss());
 
     // --- Assign the time and position of the first SIM hit in time to the accumulated hit
     if (m_btlHitsPerCell[id.rawId()][globalTrkID].time == 0. ||
         simHit.tof() < m_btlHitsPerCell[id.rawId()][globalTrkID].time) {
       m_btlHitsPerCell[id.rawId()][globalTrkID].time = simHit.tof();
 
       auto hit_pos = simHit.localPosition();
       m_btlHitsPerCell[id.rawId()][globalTrkID].x = hit_pos.x();
       m_btlHitsPerCell[id.rawId()][globalTrkID].y = hit_pos.y();
       m_btlHitsPerCell[id.rawId()][globalTrkID].z = hit_pos.z();
     }
 
   }  // simHit loop
 
   // ==============================================================================
   //  Histogram filling
   // ==============================================================================
 
   if (!m_btlHitsPerCell.empty()) {
     meNhits_->Fill(log10(m_btlHitsPerCell.size()));
   }
 
   // --- Loop over the BTL cells
   for (auto const& cell : m_btlHitsPerCell) {
     // --- Get the map of the hits in the cell
     const auto& m_hits = cell.second;
 
     // --- Skip cells with no hits
     if (m_hits.empty()) {
       continue;
     }
 
     // --- Loop over the hits in the cell, sum the hit energies and store the hit IDs in a vector
     std::vector<uint64_t> v_hitID;
     float ene_tot_cell = 0.;
 
     for (auto const& hit : m_hits) {
       ene_tot_cell += hit.second.energy;
       v_hitID.push_back(hit.first);
     }
 
     meHitLogEnergy_->Fill(log10(ene_tot_cell));
 
     // --- Skip cells with a total anergy less than hitMinEnergy_
     if (ene_tot_cell < hitMinEnergy_) {
       continue;
     }
 
     // --- Order the hits in time
     bool swapped;
     for (unsigned int ihit = 0; ihit < v_hitID.size() - 1; ++ihit) {
       swapped = false;
       for (unsigned int jhit = 0; jhit < v_hitID.size() - ihit - 1; ++jhit) {
         if (m_hits.at(v_hitID[jhit]).time > m_hits.at(v_hitID[jhit + 1]).time) {
           std::swap(v_hitID[jhit], v_hitID[jhit + 1]);
           swapped = true;
         }
       }
       if (swapped == false) {
         break;
       }
     }
 
     // --- Get the longer time interval between the hits in the cell
     float deltaT_max = 0.;
     for (unsigned int ihit = 0; ihit < v_hitID.size() - 1; ++ihit) {
       float deltaT = m_hits.at(v_hitID[ihit + 1]).time - m_hits.at(v_hitID[ihit]).time;
 
       if (deltaT > deltaT_max) {
         deltaT_max = deltaT;
       }
     }
 
     // --- Get the hit global position
     BTLDetId detId(cell.first);
     DetId geoId = detId.geographicalId(MTDTopologyMode::crysLayoutFromTopoMode(topology->getMTDTopologyMode()));
     const MTDGeomDet* thedet = geom->idToDet(geoId);
     if (thedet == nullptr) {
       throw cms::Exception("BtlSimHitsValidation") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                                    << detId.rawId() << ") is invalid!" << std::dec << std::endl;
     }
     const ProxyMTDTopology& topoproxy = static_cast<const ProxyMTDTopology&>(thedet->topology());
     const RectangularMTDTopology& topo = static_cast<const RectangularMTDTopology&>(topoproxy.specificTopology());
 
     Local3DPoint local_point(convertMmToCm(m_hits.at(v_hitID[0]).x),
                              convertMmToCm(m_hits.at(v_hitID[0]).y),
                              convertMmToCm(m_hits.at(v_hitID[0]).z));
 
     local_point = topo.pixelToModuleLocalPoint(local_point, detId.row(topo.nrows()), detId.column(topo.nrows()));
     const auto& global_point = thedet->toGlobal(local_point);
 
     // --- Fill the histograms
 
     meHitEnergy_->Fill(ene_tot_cell);
     meHitTime_->Fill(m_hits.at(v_hitID[0]).time);
 
     meHitXlocal_->Fill(m_hits.at(v_hitID[0]).x);
     meHitYlocal_->Fill(m_hits.at(v_hitID[0]).y);
     meHitZlocal_->Fill(m_hits.at(v_hitID[0]).z);
 
     meOccupancy_->Fill(global_point.z(), global_point.phi());
 
     meHitX_->Fill(global_point.x());
     meHitY_->Fill(global_point.y());
     meHitZ_->Fill(global_point.z());
     meHitPhi_->Fill(global_point.phi());
     meHitEta_->Fill(global_point.eta());
 
     meHitTvsE_->Fill(ene_tot_cell, m_hits.at(v_hitID[0]).time);
     meHitEvsPhi_->Fill(global_point.phi(), ene_tot_cell);
     meHitEvsEta_->Fill(global_point.eta(), ene_tot_cell);
     meHitEvsZ_->Fill(global_point.z(), ene_tot_cell);
     meHitTvsPhi_->Fill(global_point.phi(), m_hits.at(v_hitID[0]).time);
     meHitTvsEta_->Fill(global_point.eta(), m_hits.at(v_hitID[0]).time);
     meHitTvsZ_->Fill(global_point.z(), m_hits.at(v_hitID[0]).time);
 
     meNtrkPerCell_->Fill(v_hitID.size());
 
     // --- First hit in the cell
     int trackID = (int)(v_hitID[0] & 0xFFFFFFFF) / 100000000;
     meHitTrkID1_->Fill(trackID);
 
     // cells in the bulk of energy distribution
     if (ene_tot_cell < cellEneCut_) {
       meHitE1overEcellBulk1_->Fill(m_hits.at(v_hitID[0]).energy / ene_tot_cell);
     }
     // cells in the tail of energy distribution
     else {
       meHitE1overEcellTail1_->Fill(m_hits.at(v_hitID[0]).energy / ene_tot_cell);
     }
 
     // --- Second hit in the cell
     if (v_hitID.size() == 2) {
       trackID = (int)(v_hitID[1] & 0xFFFFFFFF) / 100000000;
       meHitTrkID2_->Fill(trackID);
 
       meHitDeltaT2_->Fill(deltaT_max);
 
       // cells in the bulk of energy distribution
       if (ene_tot_cell < cellEneCut_) {
         meHitE1overEcellBulk2_->Fill(m_hits.at(v_hitID[1]).energy / ene_tot_cell);
       }
       // cells in the tail of energy distribution
       else {
         meHitE1overEcellTail2_->Fill(m_hits.at(v_hitID[1]).energy / ene_tot_cell);
       }
 
       meHitDeltaE12vsE1_->Fill(m_hits.at(v_hitID[0]).energy,
                                m_hits.at(v_hitID[0]).energy - m_hits.at(v_hitID[1]).energy);
 
     }
     // --- Third hit in the cell
     else if (v_hitID.size() == 3) {
       trackID = (int)(v_hitID[2] & 0xFFFFFFFF) / 100000000;
       meHitTrkID3_->Fill(trackID);
 
       meHitDeltaT3_->Fill(deltaT_max);
 
       // cells in the bulk of energy distribution
       if (ene_tot_cell < cellEneCut_) {
         meHitE1overEcellBulk3_->Fill(m_hits.at(v_hitID[2]).energy / ene_tot_cell);
       }
       // cells in the tail of energy distribution
       else {
         meHitE1overEcellTail3_->Fill(m_hits.at(v_hitID[2]).energy / ene_tot_cell);
       }
 
     }
     // --- Fourth hit in the cell and next ones
     else if (v_hitID.size() >= 4) {
       for (unsigned int ihit = 3; ihit < v_hitID.size(); ++ihit) {
         trackID = (int)(v_hitID[ihit] & 0xFFFFFFFF) / 100000000;
         meHitTrkID4_->Fill(trackID);
 
         // cells in the bulk of energy distribution
         if (ene_tot_cell < cellEneCut_) {
           meHitE1overEcellBulk4_->Fill(m_hits.at(v_hitID[ihit]).energy / ene_tot_cell);
         }
         // cells in the tail of energy distribution
         else {
           meHitE1overEcellTail4_->Fill(m_hits.at(v_hitID[ihit]).energy / ene_tot_cell);
         }
       }
 
       meHitDeltaT4_->Fill(deltaT_max);
     }
 
   }  // cell loop
 
   // --- This is to count the number of processed events, needed in the harvesting step
   meNevents_->Fill(0.5);
 }

◆ bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 357 of file BtlSimHitsValidation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book1I(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), folder_, meHitDeltaE12vsE1_, meHitDeltaT2_, meHitDeltaT3_, meHitDeltaT4_, meHitE1overEcellBulk1_, meHitE1overEcellBulk2_, meHitE1overEcellBulk3_, meHitE1overEcellBulk4_, meHitE1overEcellTail1_, meHitE1overEcellTail2_, meHitE1overEcellTail3_, meHitE1overEcellTail4_, meHitEnergy_, meHitEta_, meHitEvsEta_, meHitEvsPhi_, meHitEvsZ_, meHitLogEnergy_, meHitPhi_, meHitTime_, meHitTrkID1_, meHitTrkID2_, meHitTrkID3_, meHitTrkID4_, meHitTvsE_, meHitTvsEta_, meHitTvsPhi_, meHitTvsZ_, meHitX_, meHitXlocal_, meHitY_, meHitYlocal_, meHitZ_, meHitZlocal_, meNevents_, meNhits_, meNtrkPerCell_, meOccupancy_, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                        {
   ibook.setCurrentFolder(folder_);
 
   // --- histograms booking
 
   meNevents_ = ibook.book1D("BtlNevents", "Number of events", 1, 0., 1.);
 
   meNhits_ = ibook.book1D("BtlNhits", "Number of BTL cells with SIM hits;log_{10}(N_{BTL cells})", 100, 0., 5.25);
   meNtrkPerCell_ = ibook.book1D("BtlNtrkPerCell", "Number of tracks per BTL cell;N_{trk}", 10, 0., 10.);
 
   meHitEnergy_ = ibook.book1D("BtlHitEnergy", "BTL SIM hits energy;E_{SIM} [MeV]", 100, 0., 20.);
   meHitLogEnergy_ = ibook.book1D("BtlHitLogEnergy", "BTL SIM hits energy;log_{10}(E_{SIM} [MeV])", 200, -6., 3.);
   meHitTime_ = ibook.book1D("BtlHitTime", "BTL SIM hits ToA;ToA_{SIM} [ns]", 100, 0., 25.);
 
   meHitXlocal_ = ibook.book1D("BtlHitXlocal", "BTL SIM local X;X_{SIM}^{LOC} [mm]", 100, -30., 30.);
   meHitYlocal_ = ibook.book1D("BtlHitYlocal", "BTL SIM local Y;Y_{SIM}^{LOC} [mm]", 100, -1.65, 1.65);
   meHitZlocal_ = ibook.book1D("BtlHitZlocal", "BTL SIM local Z;Z_{SIM}^{LOC} [mm]", 100, -2., 2.);
 
   meOccupancy_ = ibook.book2D(
       "BtlOccupancy", "BTL SIM hits occupancy;z_{SIM} [cm];#phi_{SIM} [rad]", 130, -260., 260., 200, -3.15, 3.15);
 
   meHitX_ = ibook.book1D("BtlHitX", "BTL SIM hits X;X_{SIM} [cm]", 100, -120., 120.);
   meHitY_ = ibook.book1D("BtlHitY", "BTL SIM hits Y;Y_{SIM} [cm]", 100, -120., 120.);
   meHitZ_ = ibook.book1D("BtlHitZ", "BTL SIM hits Z;Z_{SIM} [cm]", 100, -260., 260.);
   meHitPhi_ = ibook.book1D("BtlHitPhi", "BTL SIM hits #phi;#phi_{SIM} [rad]", 200, -3.15, 3.15);
   meHitEta_ = ibook.book1D("BtlHitEta", "BTL SIM hits #eta;#eta_{SIM}", 100, -1.55, 1.55);
 
   meHitTvsE_ =
       ibook.bookProfile("BtlHitTvsE", "BTL SIM time vs energy;E_{SIM} [MeV];T_{SIM} [ns]", 50, 0., 20., 0., 100.);
   meHitEvsPhi_ = ibook.bookProfile(
       "BtlHitEvsPhi", "BTL SIM energy vs #phi;#phi_{SIM} [rad];E_{SIM} [MeV]", 50, -3.15, 3.15, 0., 100.);
   meHitEvsEta_ =
       ibook.bookProfile("BtlHitEvsEta", "BTL SIM energy vs #eta;#eta_{SIM};E_{SIM} [MeV]", 50, -1.55, 1.55, 0., 100.);
   meHitEvsZ_ =
       ibook.bookProfile("BtlHitEvsZ", "BTL SIM energy vs Z;Z_{SIM} [cm];E_{SIM} [MeV]", 50, -260., 260., 0., 100.);
   meHitTvsPhi_ = ibook.bookProfile(
       "BtlHitTvsPhi", "BTL SIM time vs #phi;#phi_{SIM} [rad];T_{SIM} [ns]", 50, -3.15, 3.15, 0., 100.);
   meHitTvsEta_ =
       ibook.bookProfile("BtlHitTvsEta", "BTL SIM time vs #eta;#eta_{SIM};T_{SIM} [ns]", 50, -1.55, 1.55, 0., 100.);
   meHitTvsZ_ =
       ibook.bookProfile("BtlHitTvsZ", "BTL SIM time vs Z;Z_{SIM} [cm];T_{SIM} [ns]", 50, -260., 260., 0., 100.);
 
   meHitTrkID1_ = ibook.book1I("BtlHitTrkID1", "Category of the 1^{st} hit in time;Hit category", 8, 0, 8);
   meHitTrkID2_ = ibook.book1I("BtlHitTrkID2", "Category of the 2^{nd} hit in time;Hit category", 8, 0, 8);
   meHitTrkID3_ = ibook.book1I("BtlHitTrkID3", "Category of the 3^{rd} hit in time;Hit category", 8, 0, 8);
   meHitTrkID4_ = ibook.book1I("BtlHitTrkID4", "Category of the #geq4^{th} hit in time;Hit category", 8, 0, 8);
 
   meHitDeltaT2_ = ibook.book1D(
       "BtlHitDeltaT2", "Time interval between hits in the same cell (2 hits);#DeltaT_{2} [ns]", 100., 0., 25.);
   meHitDeltaT3_ = ibook.book1D(
       "BtlHitDeltaT3", "Max time interval between hits in the same cell (3 hits);#DeltaT_{3} [ns]", 100., 0., 25.);
   meHitDeltaT4_ = ibook.book1D("BtlHitDeltaT4",
                                "Max time interval between hits in the same cell (#geq4 hits);#DeltaT_{#geq4} [ns]",
                                100.,
                                0.,
                                25.);
 
   meHitDeltaE12vsE1_ = ibook.bookProfile(
       "BtlHitDeltaE12vsE", "E_{1}-E_{2} vs E_{1};E_{1} [MeV];E_{1}-E_{2} [MeV]", 100, 0., 100., 0., 100.);
 
   meHitE1overEcellBulk1_ = ibook.book1D("BtlHitE1overEtotBulk1",
                                         "Energy fraction of the 1^{st} hit in time (E_{cell}<10 MeV);E_{1} / E_{cell}",
                                         100,
                                         0.,
                                         1.);
   meHitE1overEcellBulk2_ = ibook.book1D("BtlHitE1overEtotBulk2",
                                         "Energy fraction of the 2^{nd} hit in time (E_{cell}<10 MeV);E_{2} / E_{cell}",
                                         100,
                                         0.,
                                         1.);
   meHitE1overEcellBulk3_ = ibook.book1D("BtlHitE1overEtotBulk3",
                                         "Energy fraction of the 3^{rd} hit in time (E_{cell}<10 MeV);E_{3} / E_{cell}",
                                         100,
                                         0.,
                                         1.);
   meHitE1overEcellBulk4_ =
       ibook.book1D("BtlHitE1overEtotBulk4",
                    "Energy fraction of the #geq4^{th} hit in time (E_{cell}<10 MeV);E_{#geq4} / E_{cell}",
                    100,
                    0.,
                    1.);
   meHitE1overEcellTail1_ = ibook.book1D("BtlHitE1overEtotTail1",
                                         "Energy fraction of the 1^{st} hit in time (E_{cell}>10 MeV);E_{1} / E_{cell}",
                                         100,
                                         0.,
                                         1.);
   meHitE1overEcellTail2_ = ibook.book1D("BtlHitE1overEtotTail2",
                                         "Energy fraction of the 2^{nd} hit in time (E_{cell}>10 MeV);E_{2} / E_{cell}",
                                         100,
                                         0.,
                                         1.);
   meHitE1overEcellTail3_ = ibook.book1D("BtlHitE1overEtotTail3",
                                         "Energy fraction of the 3^{rd} hit in time (E_{cell}>10 MeV);E_{3} / E_{cell}",
                                         100,
                                         0.,
                                         1.);
   meHitE1overEcellTail4_ =
       ibook.book1D("BtlHitE1overEtotTail4",
                    "Energy fraction of the #geq4^{th} hit in time (E_{cell}>10 MeV);E_{#geq4} / E_{cell}",
                    100,
                    0.,
                    1.);
 }

◆ fillDescriptions()

void BtlSimHitsValidation::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 464 of file BtlSimHitsValidation.cc.

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

                                                                                       {
   edm::ParameterSetDescription desc;
 
   desc.add<std::string>("folder", "MTD/BTL/SimHits");
   desc.add<edm::InputTag>("inputTag", edm::InputTag("mix", "g4SimHitsFastTimerHitsBarrel"));
   desc.add<double>("hitMinimumEnergy", 1.);  // [MeV]
 
   descriptions.add("btlSimHitsValid", desc);
 }