#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
edm::EDGetTokenT< GEMCoPadDigiCollection >	copad_token_

int	gem_bx_max_

int	gem_bx_min_

MEMap2Ids	me_detail_bx_

MEMap2Ids	me_detail_occ_pad_

MEMap2Ids	me_detail_occ_phi_pad_

MEMap2Ids	me_detail_occ_xy_

MEMap2Ids	me_occ_det_

MEMap1Ids	me_occ_zr_

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

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

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

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

Int_t	getDetOccBinX (Int_t chamber_id, Int_t layer_id)

const GEMGeometry *	initGeometry (const edm::EventSetup &)

Bool_t	isMuonSimHit (const PSimHit &)

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_

const std::string	kLogCategory_

const Int_t	kMuonPDGId_ = 13

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.

     : 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);
  
   gem_bx_min_ = copad_pset.getParameter<int>("minBX");
   gem_bx_max_ = copad_pset.getParameter<int>("maxBX");
 }

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

◆ ~GEMCoPadDigiValidation()

GEMCoPadDigiValidation::~GEMCoPadDigiValidation ( )

override

Definition at line 92 of file GEMCoPadDigiValidation.cc.

92 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements GEMBaseValidation.

Definition at line 94 of file GEMCoPadDigiValidation.cc.

                                                                                       {
   const GEMGeometry* gem = initGeometry(setup);
  
   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 (auto range_iter = copad_collection->begin(); range_iter != copad_collection->end(); range_iter++) {
     GEMDetId gemid = (*range_iter).first;
     const GEMCoPadDigiCollection::Range& range = (*range_iter).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();
  
     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.
       me_occ_zr_[region_id]->Fill(std::fabs(g_z1), g_r1);
       me_occ_zr_[region_id]->Fill(std::fabs(g_z2), g_r2);
  
       Int_t bin_x = getDetOccBinX(chamber_id, layer_id);
       me_occ_det_[key2]->Fill(bin_x, roll_id);
       me_occ_det_[key2]->Fill(bin_x + 1, roll_id);
  
       // Fill detail plots.
       if (detail_plot_) {
         me_detail_occ_xy_[key2]->Fill(g_x, g_y);
         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
 }

References GEMEtaPartition::centreOfPad(), GEMSuperChamber::chamber(), relativeConstraints::chamber, GEMDetId::chamber(), copad_token_, GEMBaseValidation::detail_plot_, gem_bx_max_, GEMBaseValidation::getDetOccBinX(), GEMBaseValidation::initGeometry(), edm::HandleBase::isValid(), GEMBaseValidation::kLogCategory_, GEMDetId::layer(), me_detail_bx_, me_detail_occ_pad_, me_detail_occ_phi_pad_, me_detail_occ_xy_, me_occ_det_, me_occ_zr_, 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().

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 14 of file GEMCoPadDigiValidation.cc.

                                                                         {
   const GEMGeometry* gem = initGeometry(setup);
  
   // NOTE Occupancy
   booker.setCurrentFolder("MuonGEMDigisV/GEMDigisTask/CoPad/Occupancy");
  
   for (const auto& region : gem->regions()) {
     Int_t region_id = region->region();
  
     me_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};
  
       me_occ_det_[key2] = bookDetectorOccupancy(booker, key2, station, "copad", "CoPad");
  
       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");
       }
     }  // station loop
   }    // region loop
  
   // NOTE Bunch Crossing
   if (detail_plot_) {
     booker.setCurrentFolder("MuonGEMDigisV/GEMDigisTask/CoPad/BunchCrossing");
  
     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 28 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and GEMCoPadDigiValidation().

◆ gem_bx_min_

int GEMCoPadDigiValidation::gem_bx_min_

private

Definition at line 28 of file GEMCoPadDigiValidation.h.

Referenced by 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_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_occ_det_

MEMap2Ids GEMCoPadDigiValidation::me_occ_det_

private

Definition at line 17 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and bookHistograms().

◆ me_occ_zr_

MEMap1Ids GEMCoPadDigiValidation::me_occ_zr_

private

Definition at line 16 of file GEMCoPadDigiValidation.h.

Referenced by analyze(), and bookHistograms().

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_

◆ me_detail_bx_

◆ me_detail_occ_pad_

◆ me_detail_occ_phi_pad_

◆ me_detail_occ_xy_

◆ me_occ_det_

◆ me_occ_zr_