#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 (edm::ParameterSet const &iConfig)

	DQMEDHarvester ()

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

	EDProducer (const EDProducer &)=delete

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

const EDProducer &	operator= (const EDProducer &)=delete

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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

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<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
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)

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

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

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

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

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

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

ProducesCollector	producesCollector ()

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

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

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

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

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

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

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

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)

void	resetItemsToGetFrom (BranchType iType)

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

const DTGeometry *	muonGeom

edm::ESGetToken< DTGeometry, MuonGeometryRecord >	muonGeomToken_

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
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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 29 of file DTResolutionAnalysisTest.h.

Constructor & Destructor Documentation

◆ DTResolutionAnalysisTest()

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

Constructor.

Definition at line 31 of file DTResolutionAnalysisTest.cc.

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

     : muonGeomToken_(esConsumes<edm::Transition::BeginRun>()) {
   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);
 }

◆ ~DTResolutionAnalysisTest()

DTResolutionAnalysisTest::~DTResolutionAnalysisTest ( )

override

Destructor.

Definition at line 47 of file DTResolutionAnalysisTest.cc.

References LogTrace, and nevents.

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

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.

References visDQMUpload::context, LogTrace, muonGeom, muonGeomToken_, and nevents.

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

◆ bookHistos() [1/3]

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

Definition at line 61 of file DTResolutionAnalysisTest.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), globalResSummary, fftjetimagerecorder_cfi::histoLabel, HltBtagPostValidation_cff::histoName, meanDistr, dqm::implementation::NavigatorBase::setCurrentFolder(), sigmaDistr, relativeConstraints::station, topHistoFolder, makeMuonMisalignmentScenario::wheel, wheelMeanHistos, wheelRingHistos, and wheelSigmaHistos.

Referenced by dqmEndJob().

                                                                   {
   // 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);
     }
   }
 }

◆ bookHistos() [2/3]

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

Definition at line 299 of file DTResolutionAnalysisTest.cc.

References dqm::implementation::IBooker::book2D(), HltBtagPostValidation_cff::histoName, dqm::implementation::NavigatorBase::setCurrentFolder(), topHistoFolder, makeMuonMisalignmentScenario::wheel, wheelMeanHistos, and wheelSigmaHistos.

                                                                           {
   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);
 }

◆ bookHistos() [3/3]

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

Definition at line 340 of file DTResolutionAnalysisTest.cc.

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

                                                                                     {
   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);
   }
 }

◆ dqmEndJob()

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

overrideprotectedvirtual

Implements DQMEDHarvester.

Definition at line 141 of file DTResolutionAnalysisTest.cc.

References funct::abs(), bookHistos(), DTGeometry::chambers(), dqm::implementation::IGetter::dirExists(), mps_splice::entry, dqm::impl::MonitorElement::Fill(), HcalObjRepresent::Fill(), dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getMean(), getMEName(), dqm::impl::MonitorElement::getRMS(), dqm::legacy::MonitorElement::getTH1F(), globalResSummary, LogTrace, SiStripPI::mean, meanDistr, MeanHistos, meanInRange(), muonGeom, resetMEs(), DTChamberId::sector(), sigmaDistr, SigmaHistos, sigmaInRange(), DTChamberId::station(), DTSuperLayerId::superLayer(), DTSuperLayerId::superlayer(), topHistoFolder, mps_merge::weight, DTChamberId::wheel(), wheelMeanHistos, wheelRingHistos, and wheelSigmaHistos.

                                                                                            {
   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
 }

◆ getMEName()

string DTResolutionAnalysisTest::getMEName ( const DTSuperLayerId & slID )

Get the ME name.

Definition at line 408 of file DTResolutionAnalysisTest.cc.

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

Referenced by dqmEndJob().

                                                                      {
   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;
 }

◆ meanInRange()

double DTResolutionAnalysisTest::meanInRange ( double mean ) const

private

Definition at line 445 of file DTResolutionAnalysisTest.cc.

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

Referenced by dqmEndJob().

                                                               {
   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;
 }

◆ resetMEs()

void DTResolutionAnalysisTest::resetMEs ( )

private

Definition at line 469 of file DTResolutionAnalysisTest.cc.

References globalResSummary, timingPdfMaker::histo, meanDistr, dqm::impl::MonitorElement::Reset(), sigmaDistr, wheelMeanHistos, and wheelSigmaHistos.

Referenced by dqmEndJob().

                                         {
   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();
   }
 }

◆ sigmaInRange()

double DTResolutionAnalysisTest::sigmaInRange ( double sigma ) const

private

Definition at line 457 of file DTResolutionAnalysisTest.cc.

References maxGoodSigmaValue, minBadSigmaValue, and relativeConstraints::value.

Referenced by dqmEndJob().

                                                                 {
   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;
 }

◆ slFromBin()

int DTResolutionAnalysisTest::slFromBin ( int bin ) const

private

Definition at line 437 of file DTResolutionAnalysisTest.cc.

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

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

◆ stationFromBin()

int DTResolutionAnalysisTest::stationFromBin ( int bin ) const

private

Definition at line 435 of file DTResolutionAnalysisTest.cc.

References newFWLiteAna::bin.

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

newFWLiteAna.bin

bin

Definition: newFWLiteAna.py:161

Member Data Documentation

◆ doCalibAnalysis

bool DTResolutionAnalysisTest::doCalibAnalysis

private

Definition at line 65 of file DTResolutionAnalysisTest.h.

Referenced by DTResolutionAnalysisTest(), and getMEName().

◆ globalResSummary

MonitorElement* DTResolutionAnalysisTest::globalResSummary

private

Definition at line 91 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), dqmEndJob(), and resetMEs().

◆ maxGoodMeanValue

double DTResolutionAnalysisTest::maxGoodMeanValue

private

Definition at line 60 of file DTResolutionAnalysisTest.h.

Referenced by DTResolutionAnalysisTest(), and meanInRange().

◆ maxGoodSigmaValue

double DTResolutionAnalysisTest::maxGoodSigmaValue

private

Definition at line 62 of file DTResolutionAnalysisTest.h.

Referenced by DTResolutionAnalysisTest(), and sigmaInRange().

◆ meanDistr

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

private

Definition at line 77 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), dqmEndJob(), and resetMEs().

◆ MeanHistos

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

private

Definition at line 71 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), and dqmEndJob().

◆ minBadMeanValue

double DTResolutionAnalysisTest::minBadMeanValue

private

Definition at line 61 of file DTResolutionAnalysisTest.h.

Referenced by DTResolutionAnalysisTest(), and meanInRange().

◆ minBadSigmaValue

double DTResolutionAnalysisTest::minBadSigmaValue

private

Definition at line 63 of file DTResolutionAnalysisTest.h.

Referenced by DTResolutionAnalysisTest(), and sigmaInRange().

◆ muonGeom

const DTGeometry* DTResolutionAnalysisTest::muonGeom

private

Definition at line 68 of file DTResolutionAnalysisTest.h.

Referenced by beginRun(), and dqmEndJob().

◆ muonGeomToken_

edm::ESGetToken<DTGeometry, MuonGeometryRecord> DTResolutionAnalysisTest::muonGeomToken_

private

Definition at line 67 of file DTResolutionAnalysisTest.h.

Referenced by beginRun().

◆ nevents

int DTResolutionAnalysisTest::nevents

private

Definition at line 53 of file DTResolutionAnalysisTest.h.

Referenced by beginRun(), and ~DTResolutionAnalysisTest().

◆ nLumiSegs

unsigned int DTResolutionAnalysisTest::nLumiSegs

private

Definition at line 54 of file DTResolutionAnalysisTest.h.

◆ percentual

int DTResolutionAnalysisTest::percentual

private

Definition at line 57 of file DTResolutionAnalysisTest.h.

◆ prescaleFactor

int DTResolutionAnalysisTest::prescaleFactor

private

Definition at line 55 of file DTResolutionAnalysisTest.h.

Referenced by DTResolutionAnalysisTest().

◆ run

int DTResolutionAnalysisTest::run

private

Definition at line 56 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.

Referenced by bookHistos(), dqmEndJob(), and resetMEs().

◆ SigmaHistos

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

private

Definition at line 72 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), and dqmEndJob().

◆ topHistoFolder

std::string DTResolutionAnalysisTest::topHistoFolder

private

Definition at line 94 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), dqmEndJob(), DTResolutionAnalysisTest(), and getMEName().

◆ wheelMeanHistos

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

private

Definition at line 74 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), dqmEndJob(), and resetMEs().

◆ wheelRingHistos

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

private

Definition at line 81 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), and dqmEndJob().

◆ wheelSigmaHistos

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

private

Definition at line 75 of file DTResolutionAnalysisTest.h.

Referenced by bookHistos(), dqmEndJob(), and resetMEs().

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

◆ muonGeomToken_

◆ nevents

◆ nLumiSegs

◆ percentual

◆ prescaleFactor

◆ run

◆ sigmaDistr

◆ SigmaHistos

◆ topHistoFolder

◆ wheelMeanHistos

◆ wheelRingHistos

◆ wheelSigmaHistos