#include <DTResolutionAnalysisTest.h>

Inheritance diagram for DTResolutionAnalysisTest:

Public Member Functions
void	beginRun (const edm::Run &r, const edm::EventSetup &c) override
	BeginRun. More...

void	bookHistos (DQMStore::IBooker &)

void	bookHistos (DQMStore::IBooker &, int wh)

void	bookHistos (DQMStore::IBooker &, int wh, int sect)

	DTResolutionAnalysisTest (const edm::ParameterSet &ps)
	Constructor. More...

std::string	getMEName (const DTSuperLayerId &slID)
	Get the ME name. More...

	~DTResolutionAnalysisTest () override
	Destructor. More...

Public Member Functions inherited from DQMEDHarvester
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

void	beginJob () override

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &) final

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

virtual void	dqmAnalyze (DQMStore::IBooker &, DQMStore::IGetter &, edm::Event const &, edm::EventSetup const &)

	DQMEDHarvester ()

	DQMEDHarvester (edm::ParameterSet const &iConfig)

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

virtual void	dqmEndRun (DQMStore::IBooker &, DQMStore::IGetter &, edm::Run const &, edm::EventSetup const &)

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

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &es) final

void	endProcessBlockProduce (edm::ProcessBlock &) final

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &es) final

	~DQMEDHarvester () override=default

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

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

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

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

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

Private Member Functions
double	meanInRange (double mean) const

void	resetMEs ()

double	sigmaInRange (double sigma) const

int	slFromBin (int bin) const

int	stationFromBin (int bin) const

Private Attributes
bool	doCalibAnalysis

MonitorElement *	globalResSummary

double	maxGoodMeanValue

double	maxGoodSigmaValue

std::map< int, MonitorElement * >	meanDistr

std::map< std::pair< int, int >, MonitorElement * >	MeanHistos

double	minBadMeanValue

double	minBadSigmaValue

edm::ESHandle< DTGeometry >	muonGeom

int	nevents

unsigned int	nLumiSegs

int	percentual

int	prescaleFactor

int	run

std::map< int, MonitorElement * >	sigmaDistr

std::map< std::pair< int, int >, MonitorElement * >	SigmaHistos

std::string	topHistoFolder

std::map< int, MonitorElement * >	wheelMeanHistos

std::map< int, std::map< int, std::map< std::string, MonitorElement * > > >	wheelRingHistos

std::map< int, MonitorElement * >	wheelSigmaHistos

Additional Inherited Members
Public Types inherited from DQMEDHarvester
typedef dqm::harvesting::DQMStore	DQMStore

typedef dqm::harvesting::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

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)

Protected Attributes inherited from DQMEDHarvester
DQMStore *	dqmstore_

edm::GetterOfProducts< DQMToken >	jobmegetter_

edm::EDPutTokenT< DQMToken >	jobToken_

edm::GetterOfProducts< DQMToken >	lumimegetter_

edm::EDPutTokenT< DQMToken >	lumiToken_

edm::GetterOfProducts< DQMToken >	runmegetter_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

DQM Test Client

Author
G. Mila - INFN Torino

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

Definition at line 30 of file DTResolutionAnalysisTest.h.

Constructor & Destructor Documentation

◆ DTResolutionAnalysisTest()

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

Constructor.

Definition at line 32 of file DTResolutionAnalysisTest.cc.

                                                                          {
   LogTrace("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") << "[DTResolutionAnalysisTest]: Constructor";
  
   prescaleFactor = ps.getUntrackedParameter<int>("diagnosticPrescale", 1);
   // permitted test range
   maxGoodMeanValue = ps.getUntrackedParameter<double>("maxGoodMeanValue", 0.02);
   minBadMeanValue = ps.getUntrackedParameter<double>("minBadMeanValue", 0.04);
   maxGoodSigmaValue = ps.getUntrackedParameter<double>("maxGoodSigmaValue", 0.08);
   minBadSigmaValue = ps.getUntrackedParameter<double>("minBadSigmaValue", 0.16);
   // top folder for the histograms in DQMStore
   topHistoFolder = ps.getUntrackedParameter<string>("topHistoFolder", "DT/02-Segments");
  
   doCalibAnalysis = ps.getUntrackedParameter<bool>("doCalibAnalysis", false);
 }

References dtResolutionTestFinalCalib_cfi::doCalibAnalysis, edm::ParameterSet::getUntrackedParameter(), LogTrace, dtResolutionAnalysisTest_cfi::maxGoodMeanValue, dtResolutionAnalysisTest_cfi::maxGoodSigmaValue, dtResolutionAnalysisTest_cfi::minBadMeanValue, dtResolutionAnalysisTest_cfi::minBadSigmaValue, Skims_PDWG_cff::prescaleFactor, and ALCARECODTCalibSynchCosmicsDQMClient_cff::topHistoFolder.

◆ ~DTResolutionAnalysisTest()

DTResolutionAnalysisTest::~DTResolutionAnalysisTest ( )

override

Destructor.

Definition at line 47 of file DTResolutionAnalysisTest.cc.

                                                     {
   LogTrace("DTDQM|DTMonitorClient|DTResolutionAnalysisTest")
       << "DTResolutionAnalysisTest: analyzed " << nevents << " events";
 }

References LogTrace, and cmsHarvester::nevents.

Member Function Documentation

◆ beginRun()

void DTResolutionAnalysisTest::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

override

BeginRun.

Definition at line 52 of file DTResolutionAnalysisTest.cc.

                                                                                  {
   LogTrace("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") << "[DTResolutionAnalysisTest]: BeginRun";
  
   nevents = 0;
  
   // Get the geometry
   context.get<MuonGeometryRecord>().get(muonGeom);
 }

References edm::EventSetup::get(), edm::get(), LogTrace, and cmsHarvester::nevents.

◆ bookHistos() [1/3]

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

Definition at line 61 of file DTResolutionAnalysisTest.cc.

                                                                   {
   // global residual summary
  
   ibooker.setCurrentFolder(topHistoFolder);
   globalResSummary =
       ibooker.book2D("ResidualsGlbSummary", "# of SLs with good mean and good sigma of residuals", 12, 1, 13, 5, -2, 3);
  
   // book summaries for mean and sigma
   ibooker.setCurrentFolder(topHistoFolder + "/00-MeanRes");
   meanDistr[-2] = ibooker.book1D("MeanDistr", "Mean value of the residuals all (cm)", 100, -0.1, 0.1);
   meanDistr[-1] = ibooker.book1D("MeanDistr_Phi", "Mean value of the residuals #phi SL (cm)", 100, -0.1, 0.1);
   meanDistr[0] =
       ibooker.book1D("MeanDistr_ThetaWh0", "Mean values of the residuals #theta SL Wh 0 (cm)", 100, -0.1, 0.1);
   meanDistr[1] =
       ibooker.book1D("MeanDistr_ThetaWh1", "Mean value of the residuals #theta SL Wh +/-1 (cm)", 100, -0.1, 0.1);
   meanDistr[2] =
       ibooker.book1D("MeanDistr_ThetaWh2", "Mean value of the residuals #theta SL Wh +/-2 (cm)", 100, -0.1, 0.1);
  
   for (int wh = -2; wh <= 2; wh++) {
     stringstream wheel;
     wheel << wh;
     for (int st = 1; st <= 4; st++) {
       stringstream station;
       station << st;
       for (int sl = 1; sl <= 3; sl++) {
         if (sl == 2 && st == 4)
           continue;
         ibooker.setCurrentFolder(topHistoFolder + "/00-MeanRes/Wheel" + wheel.str());
         string histoLabel = "MeanDistr";
         string histoName = histoLabel + "_W" + wheel.str() + "_MB" + station.str() + "_Phi";
         string histoType = histoLabel + "_Phi";
         if (sl == 2) {
           histoName = histoLabel + "_W" + wheel.str() + "_MB" + station.str() + "_Theta";
           histoType = histoLabel + "_Theta";
         }
         (wheelRingHistos[wh][st])[histoType] = ibooker.book1D(histoName, histoLabel, 100, -0.1, 0.1);
  
         ibooker.setCurrentFolder(topHistoFolder + "/01-SigmaRes/Wheel" + wheel.str());
         histoLabel = "SigmaDistr";
         histoName = histoLabel + "_W" + wheel.str() + "_MB" + station.str() + "_Phi";
         histoType = histoLabel + "_Phi";
         if (sl == 2) {
           histoName = histoLabel + "_W" + wheel.str() + "_MB" + station.str() + "_Theta";
           histoType = histoLabel + "_Theta";
         }
         (wheelRingHistos[wh][st])[histoType] = ibooker.book1D(histoName, histoLabel, 50, 0.0, 0.2);
       }
     }
   }
  
   ibooker.setCurrentFolder(topHistoFolder + "/00-MeanRes");
   string histoTitle = "# of SLs with good mean of residuals";
   wheelMeanHistos[3] = ibooker.book2D("MeanResGlbSummary", histoTitle.c_str(), 12, 1, 13, 5, -2, 3);
   wheelMeanHistos[3]->setAxisTitle("Sector", 1);
   wheelMeanHistos[3]->setAxisTitle("Wheel", 2);
  
   ibooker.setCurrentFolder(topHistoFolder + "/01-SigmaRes");
   sigmaDistr[-2] = ibooker.book1D("SigmaDistr", "Sigma value of the residuals all (cm)", 50, 0.0, 0.2);
   sigmaDistr[-1] = ibooker.book1D("SigmaDistr_Phi", "Sigma value of the residuals #phi SL (cm)", 50, 0.0, 0.2);
   sigmaDistr[0] =
       ibooker.book1D("SigmaDistr_ThetaWh0", "Sigma value of the residuals #theta SL Wh 0 (cm)", 50, 0.0, 0.2);
   sigmaDistr[1] =
       ibooker.book1D("SigmaDistr_ThetaWh1", "Sigma value of the residuals #theta SL Wh +/-1 (cm)", 50, 0.0, 0.2);
   sigmaDistr[2] =
       ibooker.book1D("SigmaDistr_ThetaWh2", "Sigma value of the residuals #theta SL Wh +/-2 (cm)", 50, 0.0, 0.2);
  
   histoTitle = "# of SLs with good sigma of residuals";
   wheelSigmaHistos[3] = ibooker.book2D("SigmaResGlbSummary", histoTitle.c_str(), 12, 1, 13, 5, -2, 3);
   wheelSigmaHistos[3]->setAxisTitle("Sector", 1);
   wheelSigmaHistos[3]->setAxisTitle("Wheel", 2);
  
   // loop over all the CMS wheels, sectors & book the summary histos
   for (int wheel = -2; wheel <= 2; wheel++) {
     bookHistos(ibooker, wheel);
     for (int sector = 1; sector <= 12; sector++) {
       bookHistos(ibooker, wheel, sector);
     }
   }
 }

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), bookHistos(), fftjetimagerecorder_cfi::histoLabel, HltBtagPostValidation_cff::histoName, dqm::impl::MonitorElement::setAxisTitle(), dqm::implementation::NavigatorBase::setCurrentFolder(), relativeConstraints::station, ALCARECODTCalibSynchCosmicsDQMClient_cff::topHistoFolder, and makeMuonMisalignmentScenario::wheel.

◆ bookHistos() [2/3]

void DTResolutionAnalysisTest::bookHistos	(	DQMStore::IBooker &	ibooker,
		int	wh
	)

Definition at line 299 of file DTResolutionAnalysisTest.cc.

                                                                           {
   stringstream wheel;
   wheel << wh;
  
   ibooker.setCurrentFolder(topHistoFolder + "/00-MeanRes");
   string histoName = "MeanSummaryRes_W" + wheel.str();
   string histoTitle = "# of SLs with wrong mean of residuals (Wheel " + wheel.str() + ")";
  
   wheelMeanHistos[wh] = ibooker.book2D(histoName.c_str(), histoTitle.c_str(), 12, 1, 13, 11, 1, 12);
   wheelMeanHistos[wh]->setAxisTitle("Sector", 1);
   wheelMeanHistos[wh]->setBinLabel(1, "MB1_SL1", 2);
   wheelMeanHistos[wh]->setBinLabel(2, "MB1_SL2", 2);
   wheelMeanHistos[wh]->setBinLabel(3, "MB1_SL3", 2);
   wheelMeanHistos[wh]->setBinLabel(4, "MB2_SL1", 2);
   wheelMeanHistos[wh]->setBinLabel(5, "MB2_SL2", 2);
   wheelMeanHistos[wh]->setBinLabel(6, "MB2_SL3", 2);
   wheelMeanHistos[wh]->setBinLabel(7, "MB3_SL1", 2);
   wheelMeanHistos[wh]->setBinLabel(8, "MB3_SL2", 2);
   wheelMeanHistos[wh]->setBinLabel(9, "MB3_SL3", 2);
   wheelMeanHistos[wh]->setBinLabel(10, "MB4_SL1", 2);
   wheelMeanHistos[wh]->setBinLabel(11, "MB4_SL3", 2);
  
   ibooker.setCurrentFolder(topHistoFolder + "/01-SigmaRes");
   histoName = "SigmaSummaryRes_W" + wheel.str();
   histoTitle = "# of SLs with wrong sigma of residuals (Wheel " + wheel.str() + ")";
  
   wheelSigmaHistos[wh] = ibooker.book2D(histoName.c_str(), histoTitle.c_str(), 12, 1, 13, 11, 1, 12);
   wheelSigmaHistos[wh]->setAxisTitle("Sector", 1);
   wheelSigmaHistos[wh]->setBinLabel(1, "MB1_SL1", 2);
   wheelSigmaHistos[wh]->setBinLabel(2, "MB1_SL2", 2);
   wheelSigmaHistos[wh]->setBinLabel(3, "MB1_SL3", 2);
   wheelSigmaHistos[wh]->setBinLabel(4, "MB2_SL1", 2);
   wheelSigmaHistos[wh]->setBinLabel(5, "MB2_SL2", 2);
   wheelSigmaHistos[wh]->setBinLabel(6, "MB2_SL3", 2);
   wheelSigmaHistos[wh]->setBinLabel(7, "MB3_SL1", 2);
   wheelSigmaHistos[wh]->setBinLabel(8, "MB3_SL2", 2);
   wheelSigmaHistos[wh]->setBinLabel(9, "MB3_SL3", 2);
   wheelSigmaHistos[wh]->setBinLabel(10, "MB4_SL1", 2);
   wheelSigmaHistos[wh]->setBinLabel(11, "MB4_SL3", 2);
 }

References dqm::implementation::IBooker::book2D(), HltBtagPostValidation_cff::histoName, dqm::impl::MonitorElement::setAxisTitle(), dqm::implementation::NavigatorBase::setCurrentFolder(), ALCARECODTCalibSynchCosmicsDQMClient_cff::topHistoFolder, and makeMuonMisalignmentScenario::wheel.

◆ bookHistos() [3/3]

void DTResolutionAnalysisTest::bookHistos	(	DQMStore::IBooker &	ibooker,
		int	wh,
		int	sect
	)

Definition at line 340 of file DTResolutionAnalysisTest.cc.

                                                                                     {
   stringstream wheel;
   wheel << wh;
   stringstream sector;
   sector << sect;
  
   string MeanHistoName = "MeanTest_W" + wheel.str() + "_Sec" + sector.str();
   string SigmaHistoName = "SigmaTest_W" + wheel.str() + "_Sec" + sector.str();
  
   string folder = topHistoFolder + "/Wheel" + wheel.str() + "/Sector" + sector.str();
   ibooker.setCurrentFolder(folder);
  
   if (sect != 4 && sect != 10) {
     MeanHistos[make_pair(wh, sect)] =
  
         ibooker.book1D(MeanHistoName.c_str(), "Mean (from gaussian fit) of the residuals distribution", 11, 1, 12);
   } else {
     MeanHistos[make_pair(wh, sect)] =
         ibooker.book1D(MeanHistoName.c_str(), "Mean (from gaussian fit) of the residuals distribution", 13, 1, 14);
   }
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(1, "MB1_SL1", 1);
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(2, "MB1_SL2", 1);
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(3, "MB1_SL3", 1);
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(4, "MB2_SL1", 1);
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(5, "MB2_SL2", 1);
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(6, "MB2_SL3", 1);
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(7, "MB3_SL1", 1);
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(8, "MB3_SL2", 1);
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(9, "MB3_SL3", 1);
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(10, "MB4_SL1", 1);
   (MeanHistos[make_pair(wh, sect)])->setBinLabel(11, "MB4_SL3", 1);
   if (sect == 4) {
     (MeanHistos[make_pair(wh, sect)])->setBinLabel(12, "MB4S13_SL1", 1);
     (MeanHistos[make_pair(wh, sect)])->setBinLabel(13, "MB4S13_SL3", 1);
   }
   if (sect == 10) {
     (MeanHistos[make_pair(wh, sect)])->setBinLabel(12, "MB4S14_SL1", 1);
     (MeanHistos[make_pair(wh, sect)])->setBinLabel(13, "MB4S14_SL3", 1);
   }
  
   if (sect != 4 && sect != 10) {
     SigmaHistos[make_pair(wh, sect)] =
         ibooker.book1D(SigmaHistoName.c_str(), "Sigma (from gaussian fit) of the residuals distribution", 11, 1, 12);
   } else {
     SigmaHistos[make_pair(wh, sect)] =
         ibooker.book1D(SigmaHistoName.c_str(), "Sigma (from gaussian fit) of the residuals distribution", 13, 1, 14);
   }
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(1, "MB1_SL1", 1);
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(2, "MB1_SL2", 1);
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(3, "MB1_SL3", 1);
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(4, "MB2_SL1", 1);
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(5, "MB2_SL2", 1);
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(6, "MB2_SL3", 1);
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(7, "MB3_SL1", 1);
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(8, "MB3_SL2", 1);
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(9, "MB3_SL3", 1);
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(10, "MB4_SL1", 1);
   (SigmaHistos[make_pair(wh, sect)])->setBinLabel(11, "MB4_SL3", 1);
   if (sect == 4) {
     (SigmaHistos[make_pair(wh, sect)])->setBinLabel(12, "MB4S13_SL1", 1);
     (SigmaHistos[make_pair(wh, sect)])->setBinLabel(13, "MB4S13_SL3", 1);
   }
   if (sect == 10) {
     (SigmaHistos[make_pair(wh, sect)])->setBinLabel(12, "MB4S14_SL1", 1);
     (SigmaHistos[make_pair(wh, sect)])->setBinLabel(13, "MB4S14_SL3", 1);
   }
 }

References dqm::implementation::IBooker::book1D(), printsummarytable::folder, dqm::implementation::NavigatorBase::setCurrentFolder(), ALCARECODTCalibSynchCosmicsDQMClient_cff::topHistoFolder, and makeMuonMisalignmentScenario::wheel.

◆ dqmEndJob()

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

overrideprotectedvirtual

Implements DQMEDHarvester.

Definition at line 141 of file DTResolutionAnalysisTest.cc.

                                                                                            {
   if (!igetter.dirExists(topHistoFolder)) {
     LogTrace("DTDQM|DTMonitorClient|DTResolutionAnalysisTest")
         << "[DTResolutionAnalysisTest]: Base folder " << topHistoFolder << " does not exist. Skipping client operation."
         << endl;
     return;
   }
  
   bookHistos(ibooker);  // histos booked only if top histo folder exist
   // as Standard/AlcaReco Harvest is performed in the same step
  
   LogTrace("DTDQM|DTMonitorClient|DTResolutionAnalysisTest")
       << "[DTResolutionAnalysisTest]: End of Run transition, performing the DQM client operation" << endl;
  
   // reset the ME with fixed scale
   resetMEs();
  
   for (vector<const DTChamber*>::const_iterator ch_it = muonGeom->chambers().begin();
        ch_it != muonGeom->chambers().end();
        ++ch_it) {  // loop over the chambers
  
     DTChamberId chID = (*ch_it)->id();
  
     // Fill the test histos
     for (vector<const DTSuperLayer*>::const_iterator sl_it = (*ch_it)->superLayers().begin();
          sl_it != (*ch_it)->superLayers().end();
          ++sl_it) {  // loop over SLs
  
       DTSuperLayerId slID = (*sl_it)->id();
       MonitorElement* res_histo = igetter.get(getMEName(slID));
  
       if (res_histo) {  // Gaussian Fit
         float statMean = res_histo->getMean(1);
         float statSigma = res_histo->getRMS(1);
         double mean = -1;
         double sigma = -1;
         TH1F* histo_root = res_histo->getTH1F();
  
         // fill the summaries
         int entry = (chID.station() - 1) * 3;
         int binSect = slID.sector();
         if (slID.sector() == 13)
           binSect = 4;
         else if (slID.sector() == 14)
           binSect = 10;
         int binSL = entry + slID.superLayer();
         if (chID.station() == 4 && slID.superLayer() == 3)
           binSL--;
         if ((slID.sector() == 13 || slID.sector() == 14) && slID.superLayer() == 1)
           binSL = 12;
         if ((slID.sector() == 13 || slID.sector() == 14) && slID.superLayer() == 3)
           binSL = 13;
  
         if (histo_root->GetEntries() > 20) {
           TF1* gfit = new TF1("Gaussian", "gaus", (statMean - (2 * statSigma)), (statMean + (2 * statSigma)));
           try {
             histo_root->Fit(gfit, "Q0 SERIAL", "", -0.1, 0.1);
           } catch (cms::Exception& iException) {
             LogWarning("DTDQM|DTMonitorModule|DTResolutionAnalysisTask")
                 << "[DTResolutionAnalysisTask]: Exception when fitting SL : " << slID;
             // FIXME: the SL is set as OK in the summary
             double weight = 1 / 11.;
             if ((binSect == 4 || binSect == 10) && slID.station() == 4)
               weight = 1 / 22.;
             globalResSummary->Fill(binSect, slID.wheel(), weight);
             continue;
           }
  
           if (gfit) {
             // get the mean and the sigma of the distribution
             mean = gfit->GetParameter(1);
             sigma = gfit->GetParameter(2);
  
             // fill the distributions
             meanDistr[-2]->Fill(mean);
             sigmaDistr[-2]->Fill(sigma);
             if (slID.superlayer() == 2) {
               meanDistr[abs(slID.wheel())]->Fill(mean);
               sigmaDistr[abs(slID.wheel())]->Fill(sigma);
             } else {
               meanDistr[-1]->Fill(mean);
               sigmaDistr[-1]->Fill(sigma);
             }
  
             string histoType = "MeanDistr_Phi";
             if (slID.superlayer() == 2)
               histoType = "MeanDistr_Theta";
             (wheelRingHistos[slID.wheel()][slID.station()])[histoType]->Fill(mean);
  
             histoType = "SigmaDistr_Phi";
             if (slID.superlayer() == 2)
               histoType = "SigmaDistr_Theta";
             (wheelRingHistos[slID.wheel()][slID.station()])[histoType]->Fill(sigma);
  
             // sector summaries
             MeanHistos[make_pair(slID.wheel(), binSect)]->setBinContent(binSL, mean);
             SigmaHistos[make_pair(slID.wheel(), binSect)]->setBinContent(binSL, sigma);
  
             if ((slID.sector() == 13 || slID.sector() == 14) && binSL == 12)
               binSL = 10;
             if ((slID.sector() == 13 || slID.sector() == 14) && binSL == 13)
               binSL = 11;
  
             if ((slID.sector() == 13 || slID.sector() == 14)) {
               double MeanVal = wheelMeanHistos[slID.wheel()]->getBinContent(binSect, binSL);
               double MeanBinVal = (MeanVal > 0. && MeanVal < meanInRange(mean)) ? MeanVal : meanInRange(mean);
               wheelMeanHistos[slID.wheel()]->setBinContent(binSect, binSL, MeanBinVal);
  
               double SigmaVal = wheelSigmaHistos[slID.wheel()]->getBinContent(binSect, binSL);
               double SigmaBinVal = (SigmaVal > 0. && SigmaVal < sigmaInRange(sigma)) ? SigmaVal : sigmaInRange(sigma);
               wheelSigmaHistos[slID.wheel()]->setBinContent(binSect, binSL, SigmaBinVal);
  
             } else {
               wheelMeanHistos[slID.wheel()]->setBinContent(binSect, binSL, meanInRange(mean));
               wheelSigmaHistos[slID.wheel()]->setBinContent(binSect, binSL, sigmaInRange(sigma));
             }
  
             // set the weight
             double weight = 1 / 11.;
             if ((binSect == 4 || binSect == 10) && slID.station() == 4)
               weight = 1 / 22.;
  
             // test the values of mean and sigma
             if ((meanInRange(mean) > 0.85) && (sigmaInRange(sigma) > 0.85)) {  // sigma and mean ok
               globalResSummary->Fill(binSect, slID.wheel(), weight);
               wheelMeanHistos[3]->Fill(binSect, slID.wheel(), weight);
               wheelSigmaHistos[3]->Fill(binSect, slID.wheel(), weight);
             } else {
               if ((meanInRange(mean) < 0.85) && (sigmaInRange(sigma) > 0.85)) {  // only sigma ok
                 wheelSigmaHistos[3]->Fill(binSect, slID.wheel(), weight);
               }
               if ((meanInRange(mean) > 0.85) && (sigmaInRange(sigma) < 0.85)) {  // only mean ok
                 wheelMeanHistos[3]->Fill(binSect, slID.wheel(), weight);
               }
             }
           }
           delete gfit;
         } else {
           LogVerbatim("DTDQM|DTMonitorModule|DTResolutionAnalysisTask")
               << "[DTResolutionAnalysisTask] Fit of " << slID << " not performed because # entries < 20 ";
           // FIXME: the SL is set as OK in the summary
           double weight = 1 / 11.;
           if ((binSect == 4 || binSect == 10) && slID.station() == 4)
             weight = 1 / 22.;
           globalResSummary->Fill(binSect, slID.wheel(), weight);
           wheelMeanHistos[3]->Fill(binSect, slID.wheel(), weight);
           wheelSigmaHistos[3]->Fill(binSect, slID.wheel(), weight);
           wheelMeanHistos[slID.wheel()]->setBinContent(binSect, binSL, 1.);
           wheelSigmaHistos[slID.wheel()]->setBinContent(binSect, binSL, 1.);
         }
       } else {
         LogWarning("DTDQM|DTMonitorModule|DTResolutionAnalysisTask")
             << "[DTResolutionAnalysisTask] Histo: " << getMEName(slID) << " not found" << endl;
       }
     }  // loop on SLs
   }    // Loop on Stations
 }

References funct::abs(), bookHistos(), dqm::implementation::IGetter::dirExists(), mps_splice::entry, HcalObjRepresent::Fill(), dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getMean(), dqm::impl::MonitorElement::getRMS(), dqm::legacy::MonitorElement::getTH1F(), LogTrace, SiStripPI::mean, DTChamberId::sector(), DTChamberId::station(), DTSuperLayerId::superLayer(), DTSuperLayerId::superlayer(), ALCARECODTCalibSynchCosmicsDQMClient_cff::topHistoFolder, mps_merge::weight, and DTChamberId::wheel().

◆ getMEName()

string DTResolutionAnalysisTest::getMEName ( const DTSuperLayerId & slID )

Get the ME name.

Definition at line 408 of file DTResolutionAnalysisTest.cc.

                                                                      {
   stringstream wheel;
   wheel << slID.wheel();
   stringstream station;
   station << slID.station();
   stringstream sector;
   sector << slID.sector();
   stringstream superLayer;
   superLayer << slID.superlayer();
  
   string folderName =
       topHistoFolder + "/Wheel" + wheel.str() + "/Sector" + sector.str() + "/Station" + station.str() + "/";
  
   if (doCalibAnalysis)
     folderName =
         "DT/DTCalibValidation/Wheel" + wheel.str() + "/Station" + station.str() + "/Sector" + sector.str() + "/";
  
   string histoname = folderName + "hResDist" + "_W" + wheel.str() + "_St" + station.str() + "_Sec" + sector.str() +
                      "_SL" + superLayer.str();
  
   if (doCalibAnalysis)
     histoname = folderName + "hResDist_STEP3" + "_W" + wheel.str() + "_St" + station.str() + "_Sec" + sector.str() +
                 "_SL" + superLayer.str();
  
   return histoname;
 }

References dtResolutionTestFinalCalib_cfi::doCalibAnalysis, ALCARECODTCalibSynchCosmicsDQM_cff::folderName, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTSuperLayerId::superlayer(), ALCARECODTCalibSynchCosmicsDQMClient_cff::topHistoFolder, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

◆ meanInRange()

double DTResolutionAnalysisTest::meanInRange ( double mean ) const

private

Definition at line 445 of file DTResolutionAnalysisTest.cc.

                                                               {
   double value(0.);
   if (fabs(mean) <= maxGoodMeanValue) {
     value = 1.;
   } else if (fabs(mean) > maxGoodMeanValue && fabs(mean) < minBadMeanValue) {
     value = 0.9;
   } else if (fabs(mean) >= minBadMeanValue) {
     value = 0.1;
   }
   return value;
 }

References dtResolutionAnalysisTest_cfi::maxGoodMeanValue, SiStripPI::mean, dtResolutionAnalysisTest_cfi::minBadMeanValue, and relativeConstraints::value.

◆ resetMEs()

void DTResolutionAnalysisTest::resetMEs ( )

private

Definition at line 469 of file DTResolutionAnalysisTest.cc.

                                         {
   globalResSummary->Reset();
   // Reset the summary histo
   for (map<int, MonitorElement*>::const_iterator histo = wheelMeanHistos.begin(); histo != wheelMeanHistos.end();
        histo++) {
     (*histo).second->Reset();
   }
   for (map<int, MonitorElement*>::const_iterator histo = wheelSigmaHistos.begin(); histo != wheelSigmaHistos.end();
        histo++) {
     (*histo).second->Reset();
   }
  
   for (int indx = -2; indx != 3; ++indx) {
     meanDistr[indx]->Reset();
     sigmaDistr[indx]->Reset();
   }
 }

References timingPdfMaker::histo.

◆ sigmaInRange()

double DTResolutionAnalysisTest::sigmaInRange ( double sigma ) const

private

Definition at line 457 of file DTResolutionAnalysisTest.cc.

                                                                 {
   double value(0.);
   if (sigma <= maxGoodSigmaValue) {
     value = 1.;
   } else if (sigma > maxGoodSigmaValue && sigma < minBadSigmaValue) {
     value = 0.9;
   } else if (sigma >= minBadSigmaValue) {
     value = 0.1;
   }
   return value;
 }

References dtResolutionAnalysisTest_cfi::maxGoodSigmaValue, dtResolutionAnalysisTest_cfi::minBadSigmaValue, and relativeConstraints::value.

◆ slFromBin()

int DTResolutionAnalysisTest::slFromBin ( int bin ) const

private

Definition at line 437 of file DTResolutionAnalysisTest.cc.

                                                      {
   int ret = bin % 3;
   if (ret == 0 || bin == 11)
     ret = 3;
  
   return ret;
 }

References newFWLiteAna::bin, and runTheMatrix::ret.

◆ stationFromBin()

int DTResolutionAnalysisTest::stationFromBin ( int bin ) const

private

Definition at line 435 of file DTResolutionAnalysisTest.cc.

435 { return (int)(bin / 3.1) + 1; }

References newFWLiteAna::bin.

Member Data Documentation

◆ doCalibAnalysis

bool DTResolutionAnalysisTest::doCalibAnalysis

private

Definition at line 66 of file DTResolutionAnalysisTest.h.

◆ globalResSummary

MonitorElement* DTResolutionAnalysisTest::globalResSummary

private

Definition at line 91 of file DTResolutionAnalysisTest.h.

◆ maxGoodMeanValue

double DTResolutionAnalysisTest::maxGoodMeanValue

private

Definition at line 61 of file DTResolutionAnalysisTest.h.

◆ maxGoodSigmaValue

double DTResolutionAnalysisTest::maxGoodSigmaValue

private

Definition at line 63 of file DTResolutionAnalysisTest.h.

◆ meanDistr

std::map<int, MonitorElement*> DTResolutionAnalysisTest::meanDistr

private

Definition at line 77 of file DTResolutionAnalysisTest.h.

◆ MeanHistos

std::map<std::pair<int, int>, MonitorElement*> DTResolutionAnalysisTest::MeanHistos

private

Definition at line 71 of file DTResolutionAnalysisTest.h.

◆ minBadMeanValue

double DTResolutionAnalysisTest::minBadMeanValue

private

Definition at line 62 of file DTResolutionAnalysisTest.h.

◆ minBadSigmaValue

double DTResolutionAnalysisTest::minBadSigmaValue

private

Definition at line 64 of file DTResolutionAnalysisTest.h.

◆ muonGeom

edm::ESHandle<DTGeometry> DTResolutionAnalysisTest::muonGeom

private

Definition at line 68 of file DTResolutionAnalysisTest.h.

◆ nevents

int DTResolutionAnalysisTest::nevents

private

Definition at line 54 of file DTResolutionAnalysisTest.h.

◆ nLumiSegs

unsigned int DTResolutionAnalysisTest::nLumiSegs

private

Definition at line 55 of file DTResolutionAnalysisTest.h.

◆ percentual

int DTResolutionAnalysisTest::percentual

private

Definition at line 58 of file DTResolutionAnalysisTest.h.

◆ prescaleFactor

int DTResolutionAnalysisTest::prescaleFactor

private

Definition at line 56 of file DTResolutionAnalysisTest.h.

◆ run

int DTResolutionAnalysisTest::run

private

Definition at line 57 of file DTResolutionAnalysisTest.h.

Referenced by generateEDF.LumiInfo::__str__(), DTWorkflow.DTWorkflow::all(), Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), personalPlayback.Playback::do_create_lumi(), personalPlayback.Playback::do_exec(), o2olib.O2OTool::execute(), generateEDF.LumiInfo::fixXingInfo(), and dqm_interfaces.DirWalkerFile::walk().

◆ sigmaDistr

std::map<int, MonitorElement*> DTResolutionAnalysisTest::sigmaDistr

private

Definition at line 78 of file DTResolutionAnalysisTest.h.

◆ SigmaHistos

std::map<std::pair<int, int>, MonitorElement*> DTResolutionAnalysisTest::SigmaHistos

private

Definition at line 72 of file DTResolutionAnalysisTest.h.

◆ topHistoFolder

std::string DTResolutionAnalysisTest::topHistoFolder

private

Definition at line 94 of file DTResolutionAnalysisTest.h.

◆ wheelMeanHistos

std::map<int, MonitorElement*> DTResolutionAnalysisTest::wheelMeanHistos

private

Definition at line 74 of file DTResolutionAnalysisTest.h.

◆ wheelRingHistos

std::map<int, std::map<int, std::map<std::string, MonitorElement*> > > DTResolutionAnalysisTest::wheelRingHistos

private

Definition at line 81 of file DTResolutionAnalysisTest.h.

◆ wheelSigmaHistos

std::map<int, MonitorElement*> DTResolutionAnalysisTest::wheelSigmaHistos

private

Definition at line 75 of file DTResolutionAnalysisTest.h.

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DTResolutionAnalysisTest()

◆ ~DTResolutionAnalysisTest()

Member Function Documentation

◆ beginRun()

◆ bookHistos() [1/3]

◆ bookHistos() [2/3]

◆ bookHistos() [3/3]

◆ dqmEndJob()

◆ getMEName()

◆ meanInRange()

◆ resetMEs()

◆ sigmaInRange()

◆ slFromBin()

◆ stationFromBin()

Member Data Documentation

◆ doCalibAnalysis

◆ globalResSummary

◆ maxGoodMeanValue

◆ maxGoodSigmaValue

◆ meanDistr

◆ MeanHistos

◆ minBadMeanValue

◆ minBadSigmaValue

◆ muonGeom

◆ nevents

◆ nLumiSegs

◆ percentual

◆ prescaleFactor

◆ run

◆ sigmaDistr

◆ SigmaHistos

◆ topHistoFolder

◆ wheelMeanHistos

◆ wheelRingHistos

◆ wheelSigmaHistos