#include <DTChamberEfficiencyClient.h>

Inheritance diagram for DTChamberEfficiencyClient:

Public Member Functions
	DTChamberEfficiencyClient (const edm::ParameterSet &ps)
	Constructor. More...

virtual	~DTChamberEfficiencyClient ()
	Destructor. More...

Public Member Functions inherited from DQMEDHarvester
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override final

virtual void	beginRun (edm::Run const &, edm::EventSetup const &) override

	DQMEDHarvester (void)

virtual void	endJob () override final

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override final

virtual void	endRun (edm::Run const &, edm::EventSetup const &) override

virtual void	produce (edm::Event &, edm::EventSetup const &) override final

virtual	~DQMEDHarvester ()=default

Public Member Functions inherited from edm::one::EDProducer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources, edm::EndLuminosityBlockProducer >
	EDProducer ()=default

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override

void	bookHistos (DQMStore::IBooker &)
	book the report summary More...

void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override

void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &) override

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)

Private Attributes
std::map< int, MonitorElement * >	EffDistrPerWh

MonitorElement *	globalEffDistr

MonitorElement *	globalEffSummary

edm::ESHandle< DTGeometry >	muonGeom

unsigned int	nLumiSegs

int	prescaleFactor

MonitorElement *	summaryHistos [5][2]

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

DQM Test Client

Author
M. Pelliccioni - INFN Torino

threadsafe version (//-) oct/nov 2014 - WATWanAbdullah ncpp-um-my

Definition at line 42 of file DTChamberEfficiencyClient.h.

Constructor & Destructor Documentation

DTChamberEfficiencyClient::DTChamberEfficiencyClient ( const edm::ParameterSet & ps )

Constructor.

Definition at line 30 of file DTChamberEfficiencyClient.cc.

References edm::ParameterSet::getUntrackedParameter().

 {
 
   LogVerbatim ("DTDQM|DTMonitorClient|DTChamberEfficiencyClient")
     << "DTChamberEfficiencyClient: Constructor called";
 
   prescaleFactor = pSet.getUntrackedParameter<int>("diagnosticPrescale", 1);
 
 }

DTChamberEfficiencyClient::~DTChamberEfficiencyClient ( )

virtual

Destructor.

Definition at line 40 of file DTChamberEfficiencyClient.cc.

 {
    LogVerbatim ("DTDQM|DTMonitorClient|DTChamberEfficiencyClient")
      << "DTChamberEfficiencyClient: Destructor called";
 }

Member Function Documentation

void DTChamberEfficiencyClient::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

overrideprotected

Definition at line 46 of file DTChamberEfficiencyClient.cc.

References edm::EventSetup::get().

                                                                                       {
 
   // Get the DT Geometry
   setup.get<MuonGeometryRecord>().get(muonGeom);
 
 }

void DTChamberEfficiencyClient::bookHistos ( DQMStore::IBooker & ibooker )

protected

book the report summary

Definition at line 192 of file DTChamberEfficiencyClient.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), MonitorElement::setAxisTitle(), DQMStore::IBooker::setCurrentFolder(), and makeMuonMisalignmentScenario::wheel.

 {
 
   ibooker.setCurrentFolder("DT/05-ChamberEff");
 
   globalEffSummary = ibooker.book2D("EfficiencyGlbSummary","Efficiency Summary",12,1,13,5,-2,3);
   globalEffSummary->setAxisTitle("sector",1);
   globalEffSummary->setAxisTitle("wheel",2);
 
   globalEffDistr = ibooker.book1D("TotalEfficiency","Total efficiency",51,0.,1.02);
   globalEffDistr -> setAxisTitle("Eff",1);
 
   for(int wh=-2; wh<=2; wh++){
     stringstream wheel; wheel << wh;
     string histoNameAll =  "EfficiencyMap_All_W" + wheel.str();
     string histoTitleAll =  "Efficiency map for all segments for wheel " + wheel.str();
 
     string histoNameQual =  "EfficiencyMap_Qual_W" + wheel.str();
     string histoTitleQual =  "Efficiency map for quality segments for wheel " + wheel.str();
 
     string histoNameEff =  "Efficiency_W" + wheel.str();
     string histoTitleEff =  "Segment efficiency, wheel " + wheel.str();
 
     ibooker.setCurrentFolder("DT/05-ChamberEff");
 
     summaryHistos[wh+2][0] = ibooker.book2D(histoNameAll.c_str(),histoTitleAll.c_str(),14,1.,15.,4,1.,5.);
     summaryHistos[wh+2][0]->setAxisTitle("Sector",1);
     summaryHistos[wh+2][0]->setBinLabel(1,"MB1",2);
     summaryHistos[wh+2][0]->setBinLabel(2,"MB2",2);
     summaryHistos[wh+2][0]->setBinLabel(3,"MB3",2);
     summaryHistos[wh+2][0]->setBinLabel(4,"MB4",2);
 
     EffDistrPerWh[wh+2] = ibooker.book1D(histoNameEff.c_str(),histoTitleEff.c_str(),51,0.,1.02);
     EffDistrPerWh[wh+2] -> setAxisTitle("Eff",1);
 
     ibooker.setCurrentFolder("DT/05-ChamberEff/HighQual");
 
     summaryHistos[wh+2][1] = ibooker.book2D(histoNameQual.c_str(),histoTitleQual.c_str(),14,1.,15.,4,1.,5.);
     summaryHistos[wh+2][1]->setAxisTitle("Sector",1);
     summaryHistos[wh+2][1]->setBinLabel(1,"MB1",2);
     summaryHistos[wh+2][1]->setBinLabel(2,"MB2",2);
     summaryHistos[wh+2][1]->setBinLabel(3,"MB3",2);
     summaryHistos[wh+2][1]->setBinLabel(4,"MB4",2);
    
   }
 
   return;
 }

void DTChamberEfficiencyClient::dqmEndJob	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

overrideprotectedvirtual

Implements DQMEDHarvester.

Definition at line 61 of file DTChamberEfficiencyClient.cc.

References bookHistos(), HcalObjRepresent::Fill(), DQMStore::IGetter::get(), MonitorElement::getBinContent(), MonitorElement::getBinError(), MonitorElement::getTH2F(), gen::k, LogTrace, anotherprimaryvertexanalyzer_cfi::nBinX, anotherprimaryvertexanalyzer_cfi::nBinY, funct::pow(), mathSSE::sqrt(), relativeConstraints::station, and makeMuonMisalignmentScenario::wheel.

 {
   LogVerbatim ("DTDQM|DTMonitorClient|DTChamberEfficiencyClient")
     << "DTChamberEfficiencyClient: endRun";
 
   bookHistos(ibooker);  
 
   // reset the global summary
   globalEffSummary->Reset();
 
   //Loop over the wheels
   for(int wheel=-2;wheel<=2;wheel++){
     stringstream wheel_str; wheel_str << wheel;
 
     // Get the ME produced by EfficiencyTask Source
     // All means no selection on segments, Qual means segments with at least 12 hits
     
     MonitorElement* MECountAll = igetter.get("DT/05-ChamberEff/Task/hCountSectVsChamb_All_W" + wheel_str.str());    
     MonitorElement* MECountQual = igetter.get("DT/05-ChamberEff/Task/hCountSectVsChamb_Qual_W" + wheel_str.str());  
     MonitorElement* MEExtrap = igetter.get("DT/05-ChamberEff/Task/hExtrapSectVsChamb_W" + wheel_str.str()); 
 
     //get the TH2F
     if(!MECountAll || !(MECountAll->getTH2F())) {
       edm::LogWarning("DTChamberEfficiencyClient") << "ME not available" << std::endl;
       return;
     }
       
     TH2F* hCountAll = MECountAll->getTH2F();
     TH2F* hCountQual = MECountQual->getTH2F();
     TH2F* hExtrap = MEExtrap->getTH2F();
 
 
 
     const int nBinX = summaryHistos[wheel+2][0]->getNbinsX();
     const int nBinY = summaryHistos[wheel+2][0]->getNbinsY();
 
     for(int j=1;j<=nBinX;j++){
       for(int k=1;k<=nBinY;k++){
     summaryHistos[wheel+2][0]->setBinContent(j,k,0.);
     summaryHistos[wheel+2][1]->setBinContent(j,k,0.);
 
     const float numerAll = hCountAll->GetBinContent(j,k);
     const float numerQual = hCountQual->GetBinContent(j,k);
     const float denom = hExtrap->GetBinContent(j,k);
 
     if(denom != 0.){
       const float effAll= numerAll/denom;
       const float eff_error_All = sqrt((effAll+effAll*effAll)/denom);
 
       const float effQual= numerQual/denom;
       const float eff_error_Qual = sqrt((effQual+effQual*effQual)/denom);
 
       summaryHistos[wheel+2][0]->setBinContent(j,k,effAll);
       summaryHistos[wheel+2][0]->setBinError(j,k,eff_error_All);
 
       summaryHistos[wheel+2][1]->setBinContent(j,k,effQual);
       summaryHistos[wheel+2][1]->setBinError(j,k,eff_error_Qual);
 
         // Fill 1D eff distributions
         globalEffDistr -> Fill(effAll);
         EffDistrPerWh[wheel+2] -> Fill(effAll);
 
     }
       }
     }
   }
 
   // fill the global eff. summary
   // problems at a granularity smaller than the chamber are ignored
   for(int wheel=-2; wheel<=2; wheel++) { // loop over wheels
     // retrieve the chamber efficiency summary
     MonitorElement * segmentWheelSummary = summaryHistos[wheel+2][0];
     if(segmentWheelSummary != 0) {
 
       for(int sector=1; sector<=12; sector++) { // loop over sectors
         float nFailingChambers = 0.;
 
     double meaneff = 0.;
     double errorsum = 0.;
 
     for(int station = 1; station != 5; ++station) { // loop over stations
       
       const double tmpefficiency = segmentWheelSummary->getBinContent(sector, station);
       const double tmpvariance = pow(segmentWheelSummary->getBinError(sector, station),2);
       
       if(tmpefficiency < 0.2 || tmpvariance == 0){
         nFailingChambers++;
         continue;
       }
 
       meaneff += tmpefficiency/tmpvariance;
       errorsum += 1./tmpvariance;
 
       LogTrace("DTDQM|DTMonitorClient|DTChamberEfficiencyClient")
         << "Wheel: " << wheel << " Stat: " << station
         << " Sect: " << sector << " status: " << meaneff/errorsum << endl;
     }
 
     if(sector == 4 || sector == 10) {
       int whichSector = (sector == 4) ? 13 : 14;
 
       const double tmpefficiency = segmentWheelSummary->getBinContent(whichSector, 4);
       const double tmpvariance = pow(segmentWheelSummary->getBinError(whichSector, 4),2);
 
       if(tmpefficiency > 0.2 && tmpvariance != 0) {
         meaneff += tmpefficiency/tmpvariance;
         errorsum += 1./tmpvariance;
       }
       else nFailingChambers++;
 
     }
 
     double eff_result = 0;
         if(errorsum != 0) eff_result = meaneff/errorsum;
 
     if(nFailingChambers != 0) {
       if(sector != 4 && sector != 10) eff_result = eff_result*(4.-nFailingChambers)/4.;
       else eff_result = eff_result*(5.-nFailingChambers)/5.;
     }
 
     if(eff_result > 0.7) globalEffSummary->Fill(sector,wheel,1.);
     else if(eff_result < 0.7 && eff_result > 0.5) globalEffSummary->Fill(sector,wheel,0.6);
     else if(eff_result < 0.5 && eff_result > 0.3) globalEffSummary->Fill(sector,wheel,0.4);
     else if(eff_result < 0.3 && eff_result > 0.) globalEffSummary->Fill(sector,wheel,0.15);
 
       }
     }
   }
   return;
 }

void DTChamberEfficiencyClient::dqmEndLuminosityBlock	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter,
		edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	setup
	)

overrideprotectedvirtual

Reimplemented from DQMEDHarvester.

Definition at line 54 of file DTChamberEfficiencyClient.cc.

 {
   LogVerbatim ("DTDQM|DTMonitorClient|DTChamberEfficiencyClient")
     << "DTChamberEfficiencyClient: endluminosityBlock";
 }  

Member Data Documentation

std::map<int, MonitorElement*> DTChamberEfficiencyClient::EffDistrPerWh

private

Definition at line 74 of file DTChamberEfficiencyClient.h.

MonitorElement* DTChamberEfficiencyClient::globalEffDistr

private

Definition at line 73 of file DTChamberEfficiencyClient.h.

MonitorElement* DTChamberEfficiencyClient::globalEffSummary

private

Definition at line 71 of file DTChamberEfficiencyClient.h.

edm::ESHandle<DTGeometry> DTChamberEfficiencyClient::muonGeom

private

Definition at line 67 of file DTChamberEfficiencyClient.h.

unsigned int DTChamberEfficiencyClient::nLumiSegs

private

Definition at line 64 of file DTChamberEfficiencyClient.h.

int DTChamberEfficiencyClient::prescaleFactor

private

Definition at line 65 of file DTChamberEfficiencyClient.h.

MonitorElement* DTChamberEfficiencyClient::summaryHistos[5][2]

private

Definition at line 70 of file DTChamberEfficiencyClient.h.

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation