#include <GEMCoPadDigiValidation.h>

Inheritance diagram for GEMCoPadDigiValidation:

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

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

	GEMCoPadDigiValidation (const edm::ParameterSet &)

	~GEMCoPadDigiValidation ()

Public Member Functions inherited from GEMBaseValidation
MonitorElement *	BookHistXY (DQMStore::IBooker &, const char name, const char label, unsigned int region_num, unsigned int station_num, unsigned int layer_num=99)

MonitorElement *	BookHistZR (DQMStore::IBooker &, const char name, const char label, unsigned int region_num, unsigned int station_num, unsigned int layer_num=99)

	GEMBaseValidation (const edm::ParameterSet &ps)

MonitorElement *	getDCEta (DQMStore::IBooker &, const GEMStation *, TString, TString)

MonitorElement *	getSimpleZR (DQMStore::IBooker &, TString, TString)

std::string	getStationLabel (int i)

std::string	getSuffixName (int region, int station, int layer)

std::string	getSuffixName (int region, int station)

std::string	getSuffixName (int region)

std::string	getSuffixTitle (int region, int station, int layer)

std::string	getSuffixTitle (int region, int station)

std::string	getSuffixTitle (int region)

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

unsigned int	nPart ()

unsigned int	nRegion ()

unsigned int	nStation ()

unsigned int	nStationForLabel ()

void	setNStationForLabel (unsigned int number)

virtual	~GEMBaseValidation ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
bool	detailPlot_

edm::EDGetToken	InputTagToken_

int	maxBXGEM_

int	minBXGEM_

int	nBinXY_

std::unordered_map< UInt_t, MonitorElement * >	theCoPad_dcEta

std::unordered_map< UInt_t, MonitorElement * >	theCoPad_simple_zr

MonitorElement *	theCSCCoPad [2][3]

MonitorElement *	theCSCCoPad_bx [2][3]

MonitorElement *	theCSCCoPad_phipad [2][3]

MonitorElement *	theCSCCoPad_xy [2][3]

std::unordered_map< UInt_t, MonitorElement * >	theCSCCoPad_xy_ch

MonitorElement *	theCSCCoPad_zr [2][3]

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes inherited from GEMBaseValidation
int	nBinXY_

std::vector< double >	nBinZR_

std::vector< double >	RangeZR_

Detailed Description

Definition at line 8 of file GEMCoPadDigiValidation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 4 of file GEMCoPadDigiValidation.cc.

References detailPlot_, edm::ParameterSet::getParameter(), InputTagToken_, maxBXGEM_, and minBXGEM_.

                                                                         : GEMBaseValidation(cfg)
 {
   InputTagToken_ = consumes<GEMCoPadDigiCollection>(cfg.getParameter<edm::InputTag>("CopadLabel"));
   detailPlot_ = cfg.getParameter<bool>("detailPlot");
   minBXGEM_ = cfg.getParameter<int>("minBXGEM"); 
   maxBXGEM_ = cfg.getParameter<int>("maxBXGEM"); 
 }

GEMCoPadDigiValidation::~GEMCoPadDigiValidation ( )

Definition at line 89 of file GEMCoPadDigiValidation.cc.

                                                 {
 
 }

Member Function Documentation

void GEMCoPadDigiValidation::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Implements GEMBaseValidation.

Definition at line 95 of file GEMCoPadDigiValidation.cc.

References GEMEtaPartition::centreOfPad(), detailPlot_, GEMGeometry::etaPartition(), MonitorElement::Fill(), edm::EventSetup::get(), GEMBaseValidation::getSuffixName(), GEMGeometry::idToDet(), InputTagToken_, edm::HandleBase::isValid(), LogDebug, maxBXGEM_, minBXGEM_, GEMEtaPartition::npads(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), GEMDetId::region(), reco::return(), relativeConstraints::ring, GEMDetId::roll(), AlCaHLTBitMon_QueryRunRegistry::string, GeomDet::surface(), theCoPad_dcEta, theCoPad_simple_zr, theCSCCoPad, theCSCCoPad_bx, theCSCCoPad_phipad, theCSCCoPad_xy, theCSCCoPad_xy_ch, theCSCCoPad_zr, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

 {
   const GEMGeometry* GEMGeometry_;
   try {
     edm::ESHandle<GEMGeometry> hGeom;
     iSetup.get<MuonGeometryRecord>().get(hGeom);
     GEMGeometry_ = &*hGeom;
   }
   catch( edm::eventsetup::NoProxyException<GEMGeometry>& e) {
     edm::LogError("GEMCoPadDigiValidation") << "+++ Error : GEM geometry is unavailable on event loop. +++\n";
     return;
   }
   edm::Handle<GEMCoPadDigiCollection> gem_digis;
   e.getByToken(InputTagToken_, gem_digis);
   if (!gem_digis.isValid()) {
     edm::LogError("GEMCoPadDigiValidation") << "Cannot get pads by token.";
     return ;
   }
 
 
   // Alert! copad's station2 means GEMGeometry's station3. 
   for (GEMCoPadDigiCollection::DigiRangeIterator cItr=gem_digis->begin(); cItr!=gem_digis->end(); cItr++) {
     GEMDetId id = (*cItr).first;
     int re = id.region();
     int st = id.station();
     int la = id.layer();
     Short_t chamber = (Short_t) id.chamber();
     GEMCoPadDigiCollection::const_iterator digiItr;
     //loop over digis of given roll
     for (digiItr = (*cItr ).second.first; digiItr != (*cItr ).second.second; ++digiItr)
     {
       if ( st ==2 ) st = 3; // due to gap GEMGeometry and Copad's keep information.
 
       GEMDetId roId = GEMDetId(re, id.ring(), st, la, chamber, digiItr->roll());
       Short_t nroll = roId.roll();  
       LogDebug("GEMCoPadDigiValidation")<<"roId : "<<roId;
       const GeomDet* gdet = GEMGeometry_->idToDet(roId);
       if ( gdet == nullptr) {
         edm::LogError("GEMCoPadDigiValidation")<<roId<<" : This part can not load from GEMGeometry // Original"<<id<<" station : "<<st;
         edm::LogError("GEMCoPadDigiValidation")<<"Getting DetId failed. Discard this gem copad hit.Maybe it comes from unmatched geometry between GEN and DIGI.";
         continue; 
       }
       const BoundPlane & surface = gdet->surface();
       const GEMEtaPartition * roll = GEMGeometry_->etaPartition(roId);
       LogDebug("GEMCoPadDigiValidation")<<" roll's n pad : "<<roll->npads();
 
       Short_t pad1 = (Short_t) digiItr->pad(1);
       Short_t pad2 = (Short_t) digiItr->pad(2);
       Short_t bx1  = (Short_t) digiItr->bx(1);
       Short_t bx2  = (Short_t) digiItr->bx(2);
       LogDebug("GEMCoPadDigiValidation")<<" copad #1 pad : "<<pad1<<"  bx : "<<bx1;
       LogDebug("GEMCoPadDigiValidation")<<" copad #2 pad : "<<pad2<<"  bx : "<<bx2;
 
       // Filtered using BX
       if ( bx1 < (Short_t)minBXGEM_ || bx1 > (Short_t)maxBXGEM_) continue;
       if ( bx2 < (Short_t)minBXGEM_ || bx2 > (Short_t)maxBXGEM_) continue;
 
       LocalPoint lp1 = roll->centreOfPad(pad1);
       LocalPoint lp2 = roll->centreOfPad(pad2);
 
       GlobalPoint gp1 = surface.toGlobal(lp1);
       GlobalPoint gp2 = surface.toGlobal(lp2);
       Float_t g_r1 = (Float_t) gp1.perp();
       Float_t g_r2 = (Float_t) gp2.perp();
       Float_t g_z1 = (Float_t) gp1.z();
       Float_t g_z2 = (Float_t) gp2.z();
 
       Float_t g_phi = (Float_t) gp1.phi();
       Float_t g_x = (Float_t) gp1.x();
       Float_t g_y = (Float_t) gp1.y();
 
       int region_num=0;
       if ( re == -1 ) region_num = 0 ; 
       else if (re == 1 ) region_num = 1; 
       else {
         edm::LogError("GEMCoPadDIGIValidation")<<"region : "<<re<<std::endl;
       }
       int binX = (chamber-1)*2+(la-1);
       int binY = nroll;
       int station_num = st-1;
 
       // Fill normal plots.
       TString histname_suffix = getSuffixName( re); 
       TString simple_zr_histname = TString::Format("copad_simple_zr%s",histname_suffix.Data());
       theCoPad_simple_zr[simple_zr_histname.Hash()]->Fill( fabs(g_z1), g_r1);
       theCoPad_simple_zr[simple_zr_histname.Hash()]->Fill( fabs(g_z2), g_r2);
 
       histname_suffix = getSuffixName( re, st ) ;
       TString dcEta_histname = TString::Format("copad_dcEta%s",histname_suffix.Data());
       theCoPad_dcEta[dcEta_histname.Hash()]->Fill( binX, binY); 
       theCoPad_dcEta[dcEta_histname.Hash()]->Fill( binX+1, binY); 
 
       // Fill detail plots.
       if ( detailPlot_) {
         theCSCCoPad_xy[region_num][station_num]->Fill(g_x,g_y);     
         theCSCCoPad_phipad[region_num][station_num]->Fill(g_phi,pad1);
         theCSCCoPad[region_num][station_num]->Fill(pad1);
         theCSCCoPad_bx[region_num][station_num]->Fill(bx1);
         theCSCCoPad_bx[region_num][station_num]->Fill(bx2);
         theCSCCoPad_zr[region_num][station_num]->Fill(g_z1,g_r1);
         theCSCCoPad_zr[region_num][station_num]->Fill(g_z2,g_r2);
         std::string name_prefix = getSuffixName( re, st) ;
         TString hname;
         if ( chamber %2 == 0 ) { hname = TString::Format("copad_dg_xy%s_even",name_prefix.c_str()); }
         else { hname = TString::Format("copad_dg_xy%s_odd",name_prefix.c_str()); }
         theCSCCoPad_xy_ch[hname.Hash()]->Fill(g_x,g_y);
       }
     }
   }
 }

void GEMCoPadDigiValidation::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	Run,
		edm::EventSetup const &	iSetup
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 11 of file GEMCoPadDigiValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), GEMBaseValidation::BookHistXY(), GEMBaseValidation::BookHistZR(), detailPlot_, GEMBaseValidation::getDCEta(), GEMBaseValidation::getSimpleZR(), GEMBaseValidation::getSuffixName(), GEMBaseValidation::getSuffixTitle(), GEMBaseValidation::initGeometry(), GEMBaseValidation::nStation(), PI, Pi, HLT_25ns10e33_v2_cff::region, GEMGeometry::regions(), relativeConstraints::station, AlCaHLTBitMon_QueryRunRegistry::string, theCoPad_dcEta, theCoPad_simple_zr, theCSCCoPad, theCSCCoPad_bx, theCSCCoPad_phipad, theCSCCoPad_xy, theCSCCoPad_xy_ch, and theCSCCoPad_zr.

                                                                                                                         {
   const GEMGeometry* GEMGeometry_ = initGeometry(iSetup);
   //if ( nStationForLabel() ==3 ) setNStationForLabel(2) ;
 
   const double PI = TMath::Pi();
 
   int npadsGE11 = GEMGeometry_->regions()[0]->stations()[0]->superChambers()[0]->chambers()[0]->etaPartitions()[0]->npads();
   int npadsGE21 = 0;
   int nPads = 0;
 
   if ( nStation() > 1 ) {
     //const std::vector<const GEMStation*>&  sts = GEMGeometry_->regions()[0]->stations();
     std::vector<const GEMStation*>::const_iterator it = GEMGeometry_->regions()[0]->stations().end();
     --it;
     npadsGE21  = (*it)->superChambers()[0]->chambers()[0]->etaPartitions()[0]->npads();
   }
 
   for( auto& region : GEMGeometry_->regions()  ){
     int re = region->region();
     TString title_suffix = getSuffixTitle(re); 
     TString histname_suffix = getSuffixName( re) ;
     TString simpleZR_title    = TString::Format("Copad ZR Occupancy%s; |Z|(cm) ; R(cm)",title_suffix.Data());
     TString simpleZR_histname = TString::Format("copad_simple_zr%s",histname_suffix.Data());
 
     auto* simpleZR = getSimpleZR(ibooker, simpleZR_title, simpleZR_histname);
     if ( simpleZR != nullptr) {
       theCoPad_simple_zr[simpleZR_histname.Hash() ] = simpleZR;
     }
     for( auto& station : region->stations()) {
       int st = station->station();
       //if ( nStation()==3 && st ==2 ) continue;  // skip st2 short
       //st = (station->station()==1) ? 1 : 2;   // 1 = station 1, 3 = station2l. Should be 2.
       TString title_suffix2 = getSuffixTitle( re , st) ;
       TString histname_suffix2 = getSuffixName( re, st) ;
 
       TString dcEta_title    = TString::Format("Copad's occupancy for detector component %s; # of sub-detector ;#eta-partition",title_suffix2.Data());
       TString dcEta_histname = TString::Format("copad_dcEta%s",histname_suffix2.Data());
 
       auto* dcEta = getDCEta(ibooker, station, dcEta_title, dcEta_histname);
       if ( dcEta != nullptr) {
         theCoPad_dcEta[ dcEta_histname.Hash() ] = dcEta;
       }
     }
 
     if ( detailPlot_) {
       for( auto&  region : GEMGeometry_->regions() ) {
         int re = region->region();
         int region_num = (re+1)/2;
         std::string name_prefix  = getSuffixName( re);
         std::string label_prefix = getSuffixTitle(re);
         for( auto& station : region->stations() ) {
           int st = station->station();
           int station_num = st-1;
           //if ( nStation()==3 && st ==2 ) continue;  // skip st2 short
           //st = (station->station()==1) ? 1 : 2;   // 1 = station 1, 3 = station2l. Should be 2.
             
           if ( st == 1 ) nPads = npadsGE11;
           else nPads = npadsGE21;
           name_prefix  = getSuffixName( re, st) ;
           label_prefix = getSuffixTitle( re, st) ;
           theCSCCoPad_phipad[region_num][station_num] = ibooker.book2D( ("copad_dg_phipad"+name_prefix).c_str(), ("Digi occupancy: "+label_prefix+"; phi [rad]; Pad number").c_str(), 280,-PI,PI, nPads/2,0,nPads );
           theCSCCoPad[region_num][station_num] = ibooker.book1D( ("copad_dg"+name_prefix).c_str(), ("Digi occupancy per pad number: "+label_prefix+";Pad number; entries").c_str(), nPads,0.5,nPads+0.5);
           theCSCCoPad_bx[region_num][station_num] = ibooker.book1D( ("copad_dg_bx"+name_prefix).c_str(), ("Bunch crossing: "+label_prefix+"; bunch crossing ; entries").c_str(), 11,-5.5,5.5);
           theCSCCoPad_zr[region_num][station_num] = BookHistZR( ibooker, "copad_dg","CoPad Digi",region_num,station_num);
           theCSCCoPad_xy[region_num][station_num] = BookHistXY( ibooker, "copad_dg","CoPad Digi",region_num,station_num);
           TString xy_name = TString::Format("copad_dg_xy%s_odd",name_prefix.c_str());
           TString xy_title = TString::Format("Digi XY occupancy %s at odd chambers",label_prefix.c_str());
           theCSCCoPad_xy_ch[ xy_name.Hash()] = ibooker.book2D(xy_name, xy_title, 360, -360,360, 360, -360, 360);
           xy_name = TString::Format("copad_dg_xy%s_even",name_prefix.c_str());
           xy_title = TString::Format("Digi XY occupancy %s at even chambers",label_prefix.c_str());
           theCSCCoPad_xy_ch[ xy_name.Hash()] = ibooker.book2D(xy_name, xy_title, 360, -360,360, 360, -360, 360);
         }
       }
     }
   }
 }