#include <GEMCoPadDigiValidation.h>

Inheritance diagram for GEMCoPadDigiValidation:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

	GEMCoPadDigiValidation (const edm::ParameterSet &)

	~GEMCoPadDigiValidation () override

Public Member Functions inherited from GEMBaseValidation
	GEMBaseValidation (const edm::ParameterSet &, std::string)

	~GEMBaseValidation () override=0

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

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

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 &)

Protected Member Functions inherited from GEMBaseValidation
template<typename T >
dqm::impl::MonitorElement *	bookDetectorOccupancy (DQMStore::IBooker &booker, const T &key, const GEMStation station, const char name_prfix, const char *title_prefix)

template<typename T >
dqm::impl::MonitorElement *	bookHist1D (DQMStore::IBooker &booker, const T &key, const char name, const char title, Int_t nbinsx, Double_t xlow, Double_t xup, const char x_title="", const char y_title="Entries")

template<typename T >
dqm::impl::MonitorElement *	bookHist2D (DQMStore::IBooker &booker, const T &key, const char name, const char title, Int_t nbinsx, Double_t xlow, Double_t xup, Int_t nbinsy, Double_t ylow, Double_t yup, const char x_title="", const char y_title="")

template<typename T >
dqm::impl::MonitorElement *	bookPIDHist (DQMStore::IBooker &booker, const T &key, const char name, const char title)

template<typename T >
dqm::impl::MonitorElement *	bookPIDHist (DQMStore::IBooker &booker, const T &key, Int_t ieta, const char name, const char title)

template<typename T >
dqm::impl::MonitorElement *	bookPolarOccupancy (DQMStore::IBooker &booker, const T &key, const char name_prefix, const char title_prefix)

template<typename T >
dqm::impl::MonitorElement *	bookXYOccupancy (DQMStore::IBooker &booker, const T &key, const char name_prefix, const char title_prefix)

dqm::impl::MonitorElement *	bookZROccupancy (DQMStore::IBooker &booker, Int_t region_id, const char name_prfix, const char title_prefix)

template<typename T >
dqm::impl::MonitorElement *	bookZROccupancy (DQMStore::IBooker &booker, const T &key, const char name_prfix, const char title_prefix)

Int_t	getDetOccBinX (Int_t num_layers, Int_t chamber_id, Int_t layer_id)

Int_t	getPidIdx (Int_t pid)

Bool_t	isMuonSimHit (const PSimHit &)

Float_t	toDegree (Float_t radian)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from GEMBaseValidation
Bool_t	detail_plot_

std::vector< Double_t >	eta_range_

edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomToken_

edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomTokenBeginRun_

const std::string	kLogCategory_

const Int_t	kMuonPDGId_ = 13

std::vector< Int_t >	pid_list_

Int_t	xy_occ_num_bins_

std::vector< Int_t >	zr_occ_num_bins_

std::vector< Double_t >	zr_occ_range_

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 7 of file GEMCoPadDigiValidation.h.

Constructor & Destructor Documentation

◆ GEMCoPadDigiValidation()

GEMCoPadDigiValidation::GEMCoPadDigiValidation ( const edm::ParameterSet & pset )

explicit

Definition at line 3 of file GEMCoPadDigiValidation.cc.

References copad_token_, gem_bx_max_, gem_bx_min_, geomToken_, geomTokenBeginRun_, and muonDTDigis_cfi::pset.

     : GEMBaseValidation(pset, "GEMCoPadDigiValidation") {
   const auto& copad_pset = pset.getParameterSet("gemCoPadDigi");
 
   const auto& copad_tag = copad_pset.getParameter<edm::InputTag>("inputTag");
   copad_token_ = consumes<GEMCoPadDigiCollection>(copad_tag);
   geomToken_ = esConsumes<GEMGeometry, MuonGeometryRecord>();
   geomTokenBeginRun_ = esConsumes<GEMGeometry, MuonGeometryRecord, edm::Transition::BeginRun>();
 
   gem_bx_min_ = copad_pset.getParameter<int>("minBX");
   gem_bx_max_ = copad_pset.getParameter<int>("maxBX");
 }

◆ ~GEMCoPadDigiValidation()

GEMCoPadDigiValidation::~GEMCoPadDigiValidation ( )

override

Definition at line 98 of file GEMCoPadDigiValidation.cc.

98 {}

Member Function Documentation

◆ analyze()

void GEMCoPadDigiValidation::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overridevirtual

Implements GEMBaseValidation.

Definition at line 100 of file GEMCoPadDigiValidation.cc.

References GEMEtaPartition::centreOfPad(), GEMSuperChamber::chamber(), relativeConstraints::chamber, GEMDetId::chamber(), copad_token_, GEMBaseValidation::detail_plot_, mixOne_premix_on_sim_cfi::gem, gem_bx_max_, geomToken_, GEMBaseValidation::getDetOccBinX(), edm::HandleBase::isValid(), GEMBaseValidation::kLogCategory_, GEMDetId::layer(), me_detail_bx_, me_detail_occ_det_, me_detail_occ_pad_, me_detail_occ_phi_pad_, me_detail_occ_xy_, me_detail_occ_zr_, GEMDetId::nlayers(), compareTotals::pad1, compareTotals::pad2, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), FastTimerService_cff::range, GEMDetId::region(), GEMDetId::ring(), alignCSCRings::ring_id, singleTopDQM_cfi::setup, GEMDetId::station(), GeomDet::surface(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                       {
   const auto& gemH = setup.getHandle(geomToken_);
   if (!gemH.isValid()) {
     edm::LogError(kLogCategory_) << "Failed to initialize GEM geometry.";
     return;
   }
   const GEMGeometry* gem = gemH.product();
 
   edm::Handle<GEMCoPadDigiCollection> copad_collection;
   event.getByToken(copad_token_, copad_collection);
   if (not copad_collection.isValid()) {
     edm::LogError(kLogCategory_) << "Cannot get pads by token." << std::endl;
     return;
   }
 
   // GEMCoPadDigiCollection::DigiRangeIterator
   for (const auto& copad_pair : *copad_collection) {
     GEMDetId gemid = copad_pair.first;
     const auto& range = copad_pair.second;
 
     Int_t region_id = gemid.region();
     Int_t station_id = gemid.station();
     Int_t ring_id = gemid.ring();
     Int_t layer_id = gemid.layer();
     Int_t chamber_id = gemid.chamber();
     Int_t num_layers = gemid.nlayers();
 
     ME2IdsKey key2{region_id, station_id};
 
     for (auto digi = range.first; digi != range.second; ++digi) {
       // GEM copads are stored per super chamber!
       // layer_id = 0, roll_id = 0
       GEMDetId super_chamber_id{region_id, ring_id, station_id, 0, chamber_id, 0};
       Int_t roll_id = (*digi).roll();
 
       const GeomDet* geom_det = gem->idToDet(super_chamber_id);
       if (geom_det == nullptr) {
         edm::LogError(kLogCategory_) << super_chamber_id << " : This detId cannot be "
                                      << "loaded from GEMGeometry // Original" << gemid << " station : " << station_id
                                      << std::endl
                                      << "Getting DetId failed. Discard this gem copad hit." << std::endl;
         continue;
       }
 
       const BoundPlane& surface = geom_det->surface();
       const GEMSuperChamber* super_chamber = gem->superChamber(super_chamber_id);
       if (super_chamber == nullptr) {
         edm::LogError(kLogCategory_) << "Super chamber is null for super chamber ID = " << super_chamber_id;
         continue;
       }
       Int_t pad1 = digi->pad(1);
       Int_t pad2 = digi->pad(2);
       Int_t bx1 = digi->bx(1);
       Int_t bx2 = digi->bx(2);
 
       // Filtered using BX
       if (bx1 < gem_bx_min_ or bx1 > gem_bx_max_)
         continue;
       if (bx2 < gem_bx_min_ or bx2 > gem_bx_max_)
         continue;
 
       const Int_t padArray[] = {pad1, pad2};
       LocalPoint lpArray[2];
       int ptCounter = 0;
       for (; ptCounter < 2; ++ptCounter) {
         const GEMChamber* const chamber = super_chamber->chamber(ptCounter + 1);  // Fetch chambers 1 and 2
         if (chamber == nullptr) {
           edm::LogError(kLogCategory_) << "Chamber " << (ptCounter + 1)
                                        << " is null for super chamber = " << super_chamber_id;
           break;
         }
         const GEMEtaPartition* const etaPartition = chamber->etaPartition(roll_id);
         if (etaPartition == nullptr) {
           edm::LogError(kLogCategory_) << "Eta partition " << roll_id << " is null for chamber, super chamber = ("
                                        << (ptCounter + 1) << ", " << super_chamber_id << ")";
           break;
         }
         lpArray[ptCounter] = etaPartition->centreOfPad(padArray[ptCounter]);
       }                     // end for
       if (ptCounter < 2) {  // Broke out of "for" loop
         edm::LogError(kLogCategory_) << "Skipping a digi due to bad chamber " << (ptCounter + 1);
         continue;
       }
       const GlobalPoint& gp1 = surface.toGlobal(lpArray[0]);
       const GlobalPoint& gp2 = surface.toGlobal(lpArray[1]);
 
       Float_t g_r1 = gp1.perp();
       Float_t g_r2 = gp2.perp();
 
       Float_t g_z1 = gp1.z();
       Float_t g_z2 = gp2.z();
 
       Float_t g_phi = gp1.phi();
       Float_t g_x = gp1.x();
       Float_t g_y = gp1.y();
 
       // Fill normal plots.
 
       Int_t bin_x = getDetOccBinX(num_layers, chamber_id, layer_id);
 
       // Fill detail plots.
       if (detail_plot_) {
         me_detail_occ_zr_[region_id]->Fill(std::fabs(g_z1), g_r1);
         me_detail_occ_zr_[region_id]->Fill(std::fabs(g_z2), g_r2);
         me_detail_occ_xy_[key2]->Fill(g_x, g_y);
         me_detail_occ_det_[key2]->Fill(bin_x, roll_id);
         me_detail_occ_det_[key2]->Fill(bin_x + 1, roll_id);
         me_detail_occ_phi_pad_[key2]->Fill(g_phi, pad1);
         me_detail_occ_pad_[key2]->Fill(pad1);
         me_detail_bx_[key2]->Fill(bx1);
         me_detail_bx_[key2]->Fill(bx2);
       }  // detail_plot_
     }    // loop over digis
   }      // loop over range iters
 }

◆ bookHistograms()

void GEMCoPadDigiValidation::bookHistograms	(	DQMStore::IBooker &	booker,
		edm::Run const &	run,
		edm::EventSetup const &	setup
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 16 of file GEMCoPadDigiValidation.cc.

References GEMBaseValidation::bookDetectorOccupancy(), GEMBaseValidation::bookHist1D(), GEMBaseValidation::bookHist2D(), GEMBaseValidation::bookXYOccupancy(), GEMBaseValidation::bookZROccupancy(), GEMBaseValidation::detail_plot_, mixOne_premix_on_sim_cfi::gem, geomTokenBeginRun_, GEMBaseValidation::kLogCategory_, M_PI, me_detail_bx_, me_detail_occ_det_, me_detail_occ_pad_, me_detail_occ_phi_pad_, me_detail_occ_xy_, me_detail_occ_zr_, HLT_2022v15_cff::region, dqm::implementation::NavigatorBase::setCurrentFolder(), singleTopDQM_cfi::setup, and relativeConstraints::station.

                                                                         {
   const auto& gemH = setup.getHandle(geomTokenBeginRun_);
   if (!gemH.isValid()) {
     edm::LogError(kLogCategory_) << "Failed to initialize GEM geometry.";
     return;
   }
   const GEMGeometry* gem = gemH.product();
 
   // NOTE Occupancy
   booker.setCurrentFolder("GEM/CoPad");
 
   for (const auto& region : gem->regions()) {
     Int_t region_id = region->region();
 
     if (detail_plot_)
       me_detail_occ_zr_[region_id] = bookZROccupancy(booker, region_id, "copad", "CoPad");
 
     for (const auto& station : region->stations()) {
       Int_t station_id = station->station();
       const auto& superChamberVec = station->superChambers();
       if (superChamberVec.empty() || superChamberVec[0] == nullptr) {
         edm::LogError(kLogCategory_) << "Super chambers missing or null for region = " << region_id
                                      << " and station = " << station_id;
         continue;
       }
       const auto& chamberVec = superChamberVec[0]->chambers();
       if (chamberVec.empty() || chamberVec[0] == nullptr) {
         edm::LogError(kLogCategory_) << "Chambers missing or null for region, station, super chamber = (" << region_id
                                      << ", " << station_id << ", " << superChamberVec[0]->id() << ")";
         continue;
       }
       const auto& etaPartitionVec = chamberVec[0]->etaPartitions();
       if (etaPartitionVec.empty() || etaPartitionVec[0] == nullptr) {
         edm::LogError(kLogCategory_) << "Eta partition missing or null for region, station, super chamber, chamber = ("
                                      << region_id << ", " << station_id << ", " << superChamberVec[0]->id() << ", "
                                      << chamberVec[0]->id() << ")";
         continue;
       }
       Int_t num_pads = etaPartitionVec[0]->npads();
       ME2IdsKey key2{region_id, station_id};
 
       if (detail_plot_) {
         me_detail_occ_xy_[key2] = bookXYOccupancy(booker, key2, "copad", "CoPad");
 
         me_detail_occ_phi_pad_[key2] = bookHist2D(booker,
                                                   key2,
                                                   "copad_occ_phi_pad",
                                                   "CoPad Occupancy",
                                                   280,
                                                   -M_PI,
                                                   M_PI,
                                                   num_pads / 2,
                                                   0,
                                                   num_pads,
                                                   "#phi [rad]",
                                                   "Pad number");
 
         me_detail_occ_pad_[key2] = bookHist1D(
             booker, key2, "copad_occ_pad", "CoPad Ocupancy per pad number", num_pads, 0.5, num_pads + 0.5, "Pad number");
 
         me_detail_occ_det_[key2] = bookDetectorOccupancy(booker, key2, station, "copad", "CoPad");
       }
     }  // station loop
   }    // region loop
 
   // NOTE Bunch Crossing
   if (detail_plot_) {
     for (const auto& region : gem->regions()) {
       Int_t region_id = region->region();
       for (const auto& station : region->stations()) {
         Int_t station_id = station->station();
         ME2IdsKey key2{region_id, station_id};
 
         me_detail_bx_[key2] =
             bookHist1D(booker, key2, "copad_bx", "CoPad Bunch Crossing", 5, -2.5, 2.5, "Bunch crossing");
       }  // station loop
     }    // region loop
   }      // detail plot
 }

Member Data Documentation

◆ copad_token_

edm::EDGetTokenT<GEMCoPadDigiCollection> GEMCoPadDigiValidation::copad_token_

private

Definition at line 25 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and GEMCoPadDigiValidation().

◆ gem_bx_max_

int GEMCoPadDigiValidation::gem_bx_max_

private

Definition at line 30 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and GEMCoPadDigiValidation().

◆ gem_bx_min_

int GEMCoPadDigiValidation::gem_bx_min_

private

Definition at line 30 of file GEMCoPadDigiValidation.h.

Referenced by GEMCoPadDigiValidation().

◆ geomToken_

edm::ESGetToken<GEMGeometry, MuonGeometryRecord> GEMCoPadDigiValidation::geomToken_

private

Definition at line 26 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and GEMCoPadDigiValidation().

◆ geomTokenBeginRun_

edm::ESGetToken<GEMGeometry, MuonGeometryRecord> GEMCoPadDigiValidation::geomTokenBeginRun_

private

Definition at line 27 of file GEMCoPadDigiValidation.h.

Referenced by bookHistograms(), and GEMCoPadDigiValidation().

◆ me_detail_bx_

MEMap2Ids GEMCoPadDigiValidation::me_detail_bx_

private

Definition at line 22 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and bookHistograms().

◆ me_detail_occ_det_

MEMap2Ids GEMCoPadDigiValidation::me_detail_occ_det_

private

Definition at line 17 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and bookHistograms().

◆ me_detail_occ_pad_

MEMap2Ids GEMCoPadDigiValidation::me_detail_occ_pad_

private

Definition at line 20 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and bookHistograms().

◆ me_detail_occ_phi_pad_

MEMap2Ids GEMCoPadDigiValidation::me_detail_occ_phi_pad_

private

Definition at line 19 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and bookHistograms().

◆ me_detail_occ_xy_

MEMap2Ids GEMCoPadDigiValidation::me_detail_occ_xy_

private

Definition at line 18 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and bookHistograms().

◆ me_detail_occ_zr_

MEMap1Ids GEMCoPadDigiValidation::me_detail_occ_zr_

private

Definition at line 16 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and bookHistograms().

Private Attributes
edm::EDGetTokenT< GEMCoPadDigiCollection >	copad_token_

int	gem_bx_max_

int	gem_bx_min_

edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomToken_

edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomTokenBeginRun_

MEMap2Ids	me_detail_bx_

MEMap2Ids	me_detail_occ_det_

MEMap2Ids	me_detail_occ_pad_

MEMap2Ids	me_detail_occ_phi_pad_

MEMap2Ids	me_detail_occ_xy_

MEMap1Ids	me_detail_occ_zr_

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ GEMCoPadDigiValidation()

◆ ~GEMCoPadDigiValidation()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

Member Data Documentation

◆ copad_token_

◆ gem_bx_max_

◆ gem_bx_min_

◆ geomToken_

◆ geomTokenBeginRun_

◆ me_detail_bx_

◆ me_detail_occ_det_

◆ me_detail_occ_pad_

◆ me_detail_occ_phi_pad_

◆ me_detail_occ_xy_

◆ me_detail_occ_zr_