#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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 *	meHitEnergy_

MonitorElement *	meHitEta_

MonitorElement *	meHitEvsEta_

MonitorElement *	meHitEvsPhi_

MonitorElement *	meHitEvsZ_

MonitorElement *	meHitLogEnergy_

MonitorElement *	meHitPhi_

MonitorElement *	meHitTime_

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_

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 >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	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 50 of file BtlSimHitsValidation.cc.

Constructor & Destructor Documentation

◆ BtlSimHitsValidation()

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

explicit

Definition at line 102 of file BtlSimHitsValidation.cc.

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

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

◆ ~BtlSimHitsValidation()

BtlSimHitsValidation::~BtlSimHitsValidation ( )

override

Definition at line 108 of file BtlSimHitsValidation.cc.

108 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 111 of file BtlSimHitsValidation.cc.

                                                                                       {
   using namespace edm;
   using namespace geant_units::operators;
  
   edm::ESHandle<MTDGeometry> geometryHandle;
   iSetup.get<MTDDigiGeometryRecord>().get(geometryHandle);
   const MTDGeometry* geom = geometryHandle.product();
  
   edm::ESHandle<MTDTopology> topologyHandle;
   iSetup.get<MTDTopologyRcd>().get(topologyHandle);
   const MTDTopology* topology = topologyHandle.product();
  
   auto btlSimHitsHandle = makeValid(iEvent.getHandle(btlSimHitsToken_));
   MixCollection<PSimHit> btlSimHits(btlSimHitsHandle.product());
  
   std::unordered_map<uint32_t, MTDHit> m_btlHits;
   std::unordered_map<uint32_t, std::set<int> > m_btlTrkPerCell;
  
   // --- Loop over the BLT SIM hits
   for (auto const& simHit : btlSimHits) {
     // --- Use only hits compatible with the in-time bunch-crossing
     if (simHit.tof() < 0 || simHit.tof() > 25.)
       continue;
  
     DetId id = simHit.detUnitId();
  
     m_btlTrkPerCell[id.rawId()].insert(simHit.trackId());
  
     auto simHitIt = m_btlHits.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
   // ==============================================================================
  
   if (!m_btlHits.empty())
     meNhits_->Fill(log10(m_btlHits.size()));
  
   for (auto const& hit : m_btlTrkPerCell)
     meNtrkPerCell_->Fill((hit.second).size());
  
   for (auto const& hit : m_btlHits) {
     meHitLogEnergy_->Fill(log10((hit.second).energy));
  
     if ((hit.second).energy < hitMinEnergy_)
       continue;
  
     // --- Get the SIM hit global position
     BTLDetId detId(hit.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((hit.second).x), convertMmToCm((hit.second).y), convertMmToCm((hit.second).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((hit.second).energy);
     meHitTime_->Fill((hit.second).time);
  
     meHitXlocal_->Fill((hit.second).x);
     meHitYlocal_->Fill((hit.second).y);
     meHitZlocal_->Fill((hit.second).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((hit.second).energy, (hit.second).time);
     meHitEvsPhi_->Fill(global_point.phi(), (hit.second).energy);
     meHitEvsEta_->Fill(global_point.eta(), (hit.second).energy);
     meHitEvsZ_->Fill(global_point.z(), (hit.second).energy);
     meHitTvsPhi_->Fill(global_point.phi(), (hit.second).time);
     meHitTvsEta_->Fill(global_point.eta(), (hit.second).time);
     meHitTvsZ_->Fill(global_point.z(), (hit.second).time);
  
   }  // hit loop
  
   // --- This is to count the number of processed events, needed in the harvesting step
   meNevents_->Fill(0.5);
 }

References btlSimHitsToken_, BTLDetId::column(), geant_units::operators::convertMmToCm(), geant_units::operators::convertUnitsTo(), MTDTopologyMode::crysLayoutFromTopoMode(), TauDecayModes::dec, HCALHighEnergyHPDFilter_cfi::energy, Exception, dqm::impl::MonitorElement::Fill(), BTLDetId::geographicalId(), relativeConstraints::geom, edm::EventSetup::get(), get, hitMinEnergy_, iEvent, edm::makeValid(), meHitEnergy_, meHitEta_, meHitEvsEta_, meHitEvsPhi_, meHitEvsZ_, meHitLogEnergy_, meHitPhi_, meHitTime_, meHitTvsE_, meHitTvsEta_, meHitTvsPhi_, meHitTvsZ_, meHitX_, meHitXlocal_, meHitY_, meHitYlocal_, meHitZ_, meHitZlocal_, meNevents_, meNhits_, meNtrkPerCell_, meOccupancy_, edm::ESHandle< T >::product(), DetId::rawId(), BTLDetId::row(), rpcPointValidation_cfi::simHit, ProxyMTDTopology::specificTopology(), ntuplemaker::time, GeomDet::toGlobal(), GeomDet::topology(), ecaldqm::topology(), x, y, and z.

◆ bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 219 of file BtlSimHitsValidation.cc.

                                                                        {
   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.);
 }

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

◆ fillDescriptions()

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

static

Definition at line 265 of file BtlSimHitsValidation.cc.

                                                                                       {
   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("btlSimHits", desc);
 }