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

Inheritance diagram for EtlSimHitsValidation:

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

	~EtlSimHitsValidation () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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< edm::PSimHitContainer >	etlSimHitsToken_

const std::string	folder_

const float	hitMinEnergy_

MonitorElement *	meHitEnergy_ [2]

MonitorElement *	meHitEta_ [2]

MonitorElement *	meHitEvsEta_ [2]

MonitorElement *	meHitEvsPhi_ [2]

MonitorElement *	meHitPhi_ [2]

MonitorElement *	meHitTime_ [2]

MonitorElement *	meHitTvsE_ [2]

MonitorElement *	meHitTvsEta_ [2]

MonitorElement *	meHitTvsPhi_ [2]

MonitorElement *	meHitX_ [2]

MonitorElement *	meHitXlocal_ [2]

MonitorElement *	meHitY_ [2]

MonitorElement *	meHitYlocal_ [2]

MonitorElement *	meHitZ_ [2]

MonitorElement *	meHitZlocal_ [2]

MonitorElement *	meNhits_ [2]

MonitorElement *	meNtrkPerCell_ [2]

MonitorElement *	meOccupancy_ [2]

Detailed Description

Description: ETL SIM hits validation

Implementation: [Notes on implementation]

Definition at line 41 of file EtlSimHitsValidation.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 89 of file EtlSimHitsValidation.cc.

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

     : folder_(iConfig.getParameter<std::string>("folder")),
       hitMinEnergy_(iConfig.getParameter<double>("hitMinimumEnergy")) {
   etlSimHitsToken_ = consumes<edm::PSimHitContainer>(iConfig.getParameter<edm::InputTag>("inputTag"));
 }

EtlSimHitsValidation::~EtlSimHitsValidation ( )

override

Definition at line 95 of file EtlSimHitsValidation.cc.

95 {}

Member Function Documentation

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

overrideprivate

Definition at line 98 of file EtlSimHitsValidation.cc.

References geant_units::operators::convertMmToCm(), geant_units::operators::convertUnitsTo(), TauDecayModes::dec, etlSimHitsToken_, Exception, MonitorElement::Fill(), ETLDetId::geographicalId(), relativeConstraints::geom, edm::EventSetup::get(), edm::Event::getHandle(), hitMinEnergy_, MTDGeometry::idToDet(), edm::makeValid(), meHitEnergy_, meHitEta_, meHitEvsEta_, meHitEvsPhi_, meHitPhi_, meHitTime_, meHitTvsE_, meHitTvsEta_, meHitTvsPhi_, meHitX_, meHitXlocal_, meHitY_, meHitYlocal_, meHitZ_, meHitZlocal_, meNhits_, meNtrkPerCell_, meOccupancy_, edm::ESHandle< T >::product(), DetId::rawId(), rpcPointValidation_cfi::simHit, findQualityFiles::size, ntuplemaker::time, GeomDet::toGlobal(), x, y, and z.

                                                                                       {
   using namespace edm;
   using namespace geant_units::operators;
 
   edm::ESHandle<MTDGeometry> geometryHandle;
   iSetup.get<MTDDigiGeometryRecord>().get(geometryHandle);
   const MTDGeometry* geom = geometryHandle.product();
 
   auto etlSimHitsHandle = makeValid(iEvent.getHandle(etlSimHitsToken_));
 
   std::unordered_map<uint32_t, MTDHit> m_etlHits[2];
   std::unordered_map<uint32_t, std::set<int> > m_etlTrkPerCell[2];
 
   // --- Loop over the BLT SIM hits
   for (auto const& simHit : *etlSimHitsHandle) {
     // --- Use only hits compatible with the in-time bunch-crossing
     if (simHit.tof() < 0 || simHit.tof() > 25.)
       continue;
 
     ETLDetId id = simHit.detUnitId();
 
     int idet = (id.zside() + 1) / 2;
 
     m_etlTrkPerCell[idet][id.rawId()].insert(simHit.trackId());
 
     auto simHitIt = m_etlHits[idet].emplace(id.rawId(), MTDHit()).first;
 
     // --- Accumulate the energy (in MeV) of SIM hits in the same detector cell
     (simHitIt->second).energy += convertUnitsTo(0.001_MeV, simHit.energyLoss());
 
     // --- Get the time of the first SIM hit in the cell
     if ((simHitIt->second).time == 0 || simHit.tof() < (simHitIt->second).time) {
       (simHitIt->second).time = simHit.tof();
 
       auto hit_pos = simHit.entryPoint();
       (simHitIt->second).x = hit_pos.x();
       (simHitIt->second).y = hit_pos.y();
       (simHitIt->second).z = hit_pos.z();
     }
 
   }  // simHit loop
 
   // ==============================================================================
   //  Histogram filling
   // ==============================================================================
 
   for (int idet = 0; idet < 2; ++idet) {
     meNhits_[idet]->Fill(m_etlHits[idet].size());
 
     for (auto const& hit : m_etlTrkPerCell[idet])
       meNtrkPerCell_[idet]->Fill((hit.second).size());
 
     for (auto const& hit : m_etlHits[idet]) {
       if ((hit.second).energy < hitMinEnergy_)
         continue;
 
       // --- Get the SIM hit global position
       ETLDetId detId(hit.first);
 
       DetId geoId = detId.geographicalId();
       const MTDGeomDet* thedet = geom->idToDet(geoId);
       if (thedet == nullptr)
         throw cms::Exception("EtlSimHitsValidation") << "GeographicalID: " << std::hex << geoId.rawId() << " ("
                                                      << detId.rawId() << ") is invalid!" << std::dec << std::endl;
 
       Local3DPoint local_point(
           convertMmToCm((hit.second).x), convertMmToCm((hit.second).y), convertMmToCm((hit.second).z));
       const auto& global_point = thedet->toGlobal(local_point);
 
       // --- Fill the histograms
 
       meHitEnergy_[idet]->Fill((hit.second).energy);
       meHitTime_[idet]->Fill((hit.second).time);
 
       meHitXlocal_[idet]->Fill((hit.second).x);
       meHitYlocal_[idet]->Fill((hit.second).y);
       meHitZlocal_[idet]->Fill((hit.second).z);
 
       meOccupancy_[idet]->Fill(global_point.x(), global_point.y());
 
       meHitX_[idet]->Fill(global_point.x());
       meHitY_[idet]->Fill(global_point.y());
       meHitZ_[idet]->Fill(global_point.z());
       meHitPhi_[idet]->Fill(global_point.phi());
       meHitEta_[idet]->Fill(global_point.eta());
 
       meHitTvsE_[idet]->Fill((hit.second).energy, (hit.second).time);
       meHitEvsPhi_[idet]->Fill(global_point.phi(), (hit.second).energy);
       meHitEvsEta_[idet]->Fill(global_point.eta(), (hit.second).energy);
       meHitTvsPhi_[idet]->Fill(global_point.phi(), (hit.second).time);
       meHitTvsEta_[idet]->Fill(global_point.eta(), (hit.second).time);
 
     }  // hit loop
 
   }  // idet loop
 }

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

overrideprivate

Definition at line 196 of file EtlSimHitsValidation.cc.

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

                                                                        {
   ibook.setCurrentFolder(folder_);
 
   // --- histograms booking
 
   meNhits_[1] = ibook.book1D("EtlNhitsZpos", "Number of ETL cells with SIM hits (+Z);N_{ETL cells}", 100, 0., 5000.);
   meNhits_[0] = ibook.book1D("EtlNhitsZneg", "Number of ETL cells with SIM hits (-Z);N_{ETL cells}", 100, 0., 5000.);
   meNtrkPerCell_[1] = ibook.book1D("EtlNtrkPerCellZpos", "Number of tracks per ETL sensor (+Z);N_{trk}", 10, 0., 10.);
   meNtrkPerCell_[0] = ibook.book1D("EtlNtrkPerCellZneg", "Number of tracks per ETL sensor (-Z);N_{trk}", 10, 0., 10.);
 
   meHitEnergy_[1] = ibook.book1D("EtlHitEnergyZpos", "ETL SIM hits energy (+Z);E_{SIM} [MeV]", 100, 0., 3.);
   meHitEnergy_[0] = ibook.book1D("EtlHitEnergyZneg", "ETL SIM hits energy (-Z);E_{SIM} [MeV]", 100, 0., 3.);
   meHitTime_[1] = ibook.book1D("EtlHitTimeZpos", "ETL SIM hits ToA (+Z);ToA_{SIM} [ns]", 100, 0., 25.);
   meHitTime_[0] = ibook.book1D("EtlHitTimeZneg", "ETL SIM hits ToA (-Z);ToA_{SIM} [ns]", 100, 0., 25.);
 
   meHitXlocal_[1] = ibook.book1D("EtlHitXlocalZpos", "ETL SIM local X (+Z);X_{SIM}^{LOC} [mm]", 100, -25., 25.);
   meHitXlocal_[0] = ibook.book1D("EtlHitXlocalZneg", "ETL SIM local X (-Z);X_{SIM}^{LOC} [mm]", 100, -25., 25.);
   meHitYlocal_[1] = ibook.book1D("EtlHitYlocalZpos", "ETL SIM local Y (+Z);Y_{SIM}^{LOC} [mm]", 100, -48., 48.);
   meHitYlocal_[0] = ibook.book1D("EtlHitYlocalZneg", "ETL SIM local Y (-Z);Y_{SIM}^{LOC} [mm]", 100, -48., 48.);
   meHitZlocal_[1] = ibook.book1D("EtlHitZlocalZpos", "ETL SIM local Z (+Z);Z_{SIM}^{LOC} [mm]", 80, -0.16, 0.16);
   meHitZlocal_[0] = ibook.book1D("EtlHitZlocalZneg", "ETL SIM local Z (-Z);Z_{SIM}^{LOC} [mm]", 80, -0.16, 0.16);
 
   meOccupancy_[1] = ibook.book2D(
       "EtlOccupancyZpos", "ETL SIM hits occupancy (+Z);X_{SIM} [cm];Y_{SIM} [cm]", 135, -135., 135., 135, -135., 135.);
   meOccupancy_[0] = ibook.book2D(
       "EtlOccupancyZneg", "ETL SIM hits occupancy (-Z);X_{SIM} [cm];Y_{SIM} [cm]", 135, -135., 135., 135, -135., 135.);
 
   meHitX_[1] = ibook.book1D("EtlHitXZpos", "ETL SIM hits X (+Z);X_{SIM} [cm]", 100, -130., 130.);
   meHitX_[0] = ibook.book1D("EtlHitXZneg", "ETL SIM hits X (-Z);X_{SIM} [cm]", 100, -130., 130.);
   meHitY_[1] = ibook.book1D("EtlHitYZpos", "ETL SIM hits Y (+Z);Y_{SIM} [cm]", 100, -130., 130.);
   meHitY_[0] = ibook.book1D("EtlHitYZneg", "ETL SIM hits Y (-Z);Y_{SIM} [cm]", 100, -130., 130.);
   meHitZ_[1] = ibook.book1D("EtlHitZZpos", "ETL SIM hits Z (+Z);Z_{SIM} [cm]", 100, 303.4, 304.2);
   meHitZ_[0] = ibook.book1D("EtlHitZZneg", "ETL SIM hits Z (-Z);Z_{SIM} [cm]", 100, -304.2, -303.4);
 
   meHitPhi_[1] = ibook.book1D("EtlHitPhiZpos", "ETL SIM hits #phi (+Z);#phi_{SIM} [rad]", 100, -3.15, 3.15);
   meHitPhi_[0] = ibook.book1D("EtlHitPhiZneg", "ETL SIM hits #phi (-Z);#phi_{SIM} [rad]", 100, -3.15, 3.15);
   meHitEta_[1] = ibook.book1D("EtlHitEtaZpos", "ETL SIM hits #eta (+Z);#eta_{SIM}", 100, 1.56, 3.2);
   meHitEta_[0] = ibook.book1D("EtlHitEtaZneg", "ETL SIM hits #eta (-Z);#eta_{SIM}", 100, -3.2, -1.56);
 
   meHitTvsE_[1] = ibook.bookProfile(
       "EtlHitTvsEZpos", "ETL SIM time vs energy (+Z);E_{SIM} [MeV];T_{SIM} [ns]", 50, 0., 2., 0., 100.);
   meHitTvsE_[0] = ibook.bookProfile(
       "EtlHitTvsEZneg", "ETL SIM time vs energy (-Z);E_{SIM} [MeV];T_{SIM} [ns]", 50, 0., 2., 0., 100.);
   meHitEvsPhi_[1] = ibook.bookProfile(
       "EtlHitEvsPhiZpos", "ETL SIM energy vs #phi (+Z);#phi_{SIM} [rad];E_{SIM} [MeV]", 50, -3.15, 3.15, 0., 100.);
   meHitEvsPhi_[0] = ibook.bookProfile(
       "EtlHitEvsPhiZneg", "ETL SIM energy vs #phi (-Z);#phi_{SIM} [rad];E_{SIM} [MeV]", 50, -3.15, 3.15, 0., 100.);
   meHitEvsEta_[1] = ibook.bookProfile(
       "EtlHitEvsEtaZpos", "ETL SIM energy vs #eta (+Z);#eta_{SIM};E_{SIM} [MeV]", 50, 1.56, 3.2, 0., 100.);
   meHitEvsEta_[0] = ibook.bookProfile(
       "EtlHitEvsEtaZneg", "ETL SIM energy vs #eta (-Z);#eta_{SIM};E_{SIM} [MeV]", 50, -3.2, -1.56, 0., 100.);
   meHitTvsPhi_[1] = ibook.bookProfile(
       "EtlHitTvsPhiZpos", "ETL SIM time vs #phi (+Z);#phi_{SIM} [rad];T_{SIM} [ns]", 50, -3.15, 3.15, 0., 100.);
   meHitTvsPhi_[0] = ibook.bookProfile(
       "EtlHitTvsPhiZneg", "ETL SIM time vs #phi (-Z);#phi_{SIM} [rad];T_{SIM} [ns]", 50, -3.15, 3.15, 0., 100.);
   meHitTvsEta_[1] = ibook.bookProfile(
       "EtlHitTvsEtaZpos", "ETL SIM time vs #eta (+Z);#eta_{SIM};T_{SIM} [ns]", 50, 1.56, 3.2, 0., 100.);
   meHitTvsEta_[0] = ibook.bookProfile(
       "EtlHitTvsEtaZneg", "ETL SIM time vs #eta (-Z);#eta_{SIM};T_{SIM} [ns]", 50, -3.2, -1.56, 0., 100.);
 }

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

static

Definition at line 260 of file EtlSimHitsValidation.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), DEFINE_FWK_MODULE, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                       {
   edm::ParameterSetDescription desc;
 
   desc.add<std::string>("folder", "MTD/ETL/SimHits");
   desc.add<edm::InputTag>("inputTag", edm::InputTag("g4SimHits", "FastTimerHitsEndcap"));
   desc.add<double>("hitMinimumEnergy", 0.1);  // [MeV]
 
   descriptions.add("etlSimHits", desc);
 }