#include <DTLocalTriggerSynchTest.h>

Inheritance diagram for DTLocalTriggerSynchTest:

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

virtual	~DTLocalTriggerSynchTest ()
	Destructor. More...

Public Member Functions inherited from DTLocalTriggerBaseTest
	DTLocalTriggerBaseTest ()
	Constructor. More...

virtual	~DTLocalTriggerBaseTest ()
	Destructor. More...

Public Member Functions inherited from DQMEDHarvester
virtual void	analyze (edm::Event const &, edm::EventSetup const &) final

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

	DQMEDHarvester (void)

virtual void	endJob () final

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

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

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

	EDConsumerBase ()

	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 &run, const edm::EventSetup &c)
	begin Run More...

void	bookChambHistos (DQMStore::IBooker &, DTChamberId chambId, std::string htype, std::string subfolder="")
	Book the new MEs (for each chamber) More...

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

void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)
	Perform client diagnostic in online. More...

float	getFloatFromME (DQMStore::IGetter &, DTChamberId chId, std::string meType)
	Get float MEs. More...

void	makeRatioME (TH1F numerator, TH1F denominator, MonitorElement *result)
	Compute efficiency plots. More...

void	runClientDiagnostic (DQMStore::IBooker &, DQMStore::IGetter &)
	DQM Client Diagnostic. More...

Protected Member Functions inherited from DTLocalTriggerBaseTest
void	beginRun (edm::Run const &run, edm::EventSetup const &context)
	BeginRun. More...

void	bookCmsHistos (DQMStore::IBooker &, std::string hTag, std::string folder="", bool isGlb=false)
	Book the new MEs (CMS summary) More...

void	bookSectorHistos (DQMStore::IBooker &, int wheel, int sector, std::string hTag, std::string folder="")
	Book the new MEs (for each sector) More...

void	bookWheelHistos (DQMStore::IBooker &, int wheel, std::string hTag, std::string folder="")
	Book the new MEs (for each wheel) More...

std::string	category ()
	Get message logger name. More...

void	endRun (edm::Run const &run, edm::EventSetup const &context)
	Perform client diagnostic in offline. More...

std::string	fullName (std::string htype)
	Create fullname from histo partial name. More...

template<class T >
T *	getHisto (MonitorElement *me)
	Convert ME to Histogram fo type T. More...

std::string	getMEName (std::string histoTag, std::string subfolder, const DTChamberId &chambid)
	Get the ME name (by chamber) More...

std::string	getMEName (std::string histoTag, std::string subfolder, int wh)
	Get the ME name (by wheel) More...

std::pair< float, float >	phiRange (const DTChamberId &id)
	Calculate phi range for histograms. More...

void	setConfig (const edm::ParameterSet &ps, std::string name)
	Set configuration variables. More...

std::string &	topFolder (bool isTM)
	Get top folder name. More...

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
bool	bookingdone

double	bxTime

std::map< uint32_t, std::map< std::string, MonitorElement * > >	chambME

std::string	denHistoTag

int	minEntries

int	nBXHigh

int	nBXLow

std::string	numHistoTag

bool	rangeInBX

std::string	ratioHistoTag

DTTPGParameters	wPhaseMap

bool	writeDB

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from DTLocalTriggerBaseTest
std::string	baseFolderDDU

std::string	baseFolderTM

std::map< std::string, MonitorElement * >	cmsME

std::string	hwSource

std::vector< std::string >	hwSources

edm::ESHandle< DTGeometry >	muonGeom

int	nevents

unsigned int	nLumiSegs

edm::ParameterSet	parameters

int	prescaleFactor

int	run

bool	runOnline

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

std::string	sourceFolder

std::string	testName

std::string	trigSource

std::vector< std::string >	trigSources

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

Detailed Description

DQM Test Client

Author
C. Battilana - CIEMAT

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

Definition at line 22 of file DTLocalTriggerSynchTest.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 45 of file DTLocalTriggerSynchTest.cc.

                                                                           {
 
   setConfig(ps,"DTLocalTriggerSynch");
   baseFolderTM = "DT/90-LocalTriggerSynch/";
   baseFolderDDU = "DT/90-LocalTriggerSynch/";
 
   bookingdone = 0;
 
 }

DTLocalTriggerSynchTest::~DTLocalTriggerSynchTest ( )

virtual

Destructor.

Definition at line 56 of file DTLocalTriggerSynchTest.cc.

                                                  {
 
 }

Member Function Documentation

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

protected

begin Run

Definition at line 61 of file DTLocalTriggerSynchTest.cc.

References DTLocalTriggerBaseTest::beginRun(), hgcalDigitizer_cfi::bxTime, and metProducer_cfi::parameters.

                                                                                   {  
 
   numHistoTag   = parameters.getParameter<string>("numHistoTag");
   denHistoTag   = parameters.getParameter<string>("denHistoTag");
   ratioHistoTag = parameters.getParameter<string>("ratioHistoTag");
   bxTime        = parameters.getParameter<double>("bxTimeInterval");
   rangeInBX     = parameters.getParameter<bool>("rangeWithinBX");
   nBXLow        = parameters.getParameter<int>("nBXLow");
   nBXHigh       = parameters.getParameter<int>("nBXHigh");
   minEntries    = parameters.getParameter<int>("minEntries");
 
   DTLocalTriggerBaseTest::beginRun(run,c);
 }

void DTLocalTriggerSynchTest::bookChambHistos	(	DQMStore::IBooker &	,
		DTChamberId	chambId,
		std::string	htype,
		std::string	subfolder = `""`
	)

protected

Book the new MEs (for each chamber)

Definition at line 261 of file DTLocalTriggerSynchTest.cc.

References DQMStore::IBooker::book1D(), hgcalDigitizer_cfi::bxTime, python.rootplot.argparse::category, printsummarytable::folder, edm::max(), edm::min(), pileupCalc::nbins, DetId::rawId(), DTChamberId::sector(), DQMStore::IBooker::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

                                                                                                       {
   
   stringstream wheel; wheel << chambId.wheel();
   stringstream station; station << chambId.station();   
   stringstream sector; sector << chambId.sector();
 
   string fullType  = fullName(htype);
   bool isTM = hwSource=="TM" ;
   string HistoName = fullType + "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
 
   string folder = topFolder(isTM) + "Wheel" + wheel.str() + "/Sector" + sector.str() + "/Station" + station.str();
   if ( subfolder!="") { folder += "7" + subfolder; }
 
   ibooker.setCurrentFolder(folder);
 
   LogPrint(category()) << "[" << testName << "Test]: booking " << folder << "/" <<HistoName;
 
   
   uint32_t indexChId = chambId.rawId();
   float min = rangeInBX ?      0 : nBXLow*bxTime;
   float max = rangeInBX ? bxTime : nBXHigh*bxTime;
   int nbins = static_cast<int>(ceil( rangeInBX ? bxTime : (nBXHigh-nBXLow)*bxTime));
 
   chambME[indexChId][fullType] = ibooker.book1D(HistoName.c_str(),"All/HH ratio vs Muon Arrival Time",nbins,min,max);
 
 }

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

protectedvirtual

Reimplemented from DTLocalTriggerBaseTest.

Definition at line 165 of file DTLocalTriggerSynchTest.cc.

References DTTPGParameters::begin(), hgcalDigitizer_cfi::bxTime, python.rootplot.argparse::category, DTLocalTriggerBaseTest::dqmEndJob(), DTTPGParameters::end(), DQMStore::IGetter::get(), DTTimeUnits::ns, metProducer_cfi::parameters, DTTPGParameters::set(), and DTCalibDBUtils::writeToDB().

                                                                                            {
 
   DTLocalTriggerBaseTest::dqmEndJob(ibooker,igetter);
 
   if ( parameters.getParameter<bool>("writeDB")) {
     LogVerbatim(category()) << "[" << testName 
                 << "Test]: writeDB flag set to true. Producing peak position database." << endl;
 
     DTTPGParameters* delayMap = new DTTPGParameters();
     hwSource =  parameters.getParameter<bool>("dbFromTM") ? "TM" : "DDU";
     std::vector<const DTChamber*>::const_iterator chambIt  = muonGeom->chambers().begin();
     std::vector<const DTChamber*>::const_iterator chambEnd = muonGeom->chambers().end();
       for (; chambIt!=chambEnd; ++chambIt) { 
 
     DTChamberId chId = (*chambIt)->id();
     float fineDelay = 0;
     int coarseDelay = static_cast<int>((getFloatFromME(igetter,chId,"tTrig_SL1") + 
                                                                         getFloatFromME(igetter,chId,"tTrig_SL3"))*0.5/bxTime);
 
     bool fineDiff   = parameters.getParameter<bool>("fineParamDiff");
     bool coarseDiff = parameters.getParameter<bool>("coarseParamDiff");
 
     TH1F *ratioH     = getHisto<TH1F>(igetter.get(getMEName(ratioHistoTag,"", chId)));    
     if (ratioH->GetEntries()>minEntries) {        
       TF1 *fitF=ratioH->GetFunction("mypol8");
       if (fitF) { fineDelay=fitF->GetMaximumX(0,bxTime); }
     } else {
       LogInfo(category()) << "[" << testName 
                   << "Test]: Ratio histogram for chamber " << chId
                   << " is empty. Worst Phase value set to 0." << endl;
     }
 
     if (fineDiff || coarseDiff) {
       float wFine;
       int wCoarse;
       wPhaseMap.get(chId,wCoarse,wFine,DTTimeUnits::ns);
       if (fineDiff)   { fineDelay = wFine - fineDelay; }
       if (coarseDiff) { coarseDelay = wCoarse - coarseDelay; }
     } 
     delayMap->set(chId,coarseDelay,fineDelay,DTTimeUnits::ns);
       }
 
       std::vector< std::pair<DTTPGParametersId,DTTPGParametersData> >::const_iterator dbIt  = delayMap->begin();
       std::vector< std::pair<DTTPGParametersId,DTTPGParametersData> >::const_iterator dbEnd = delayMap->end();
       for (; dbIt!=dbEnd; ++dbIt) {
     LogVerbatim(category()) << "[" << testName << "Test]: DB entry for Wh " << (*dbIt).first.wheelId 
                 << " Sec " << (*dbIt).first.sectorId 
                 << " St " << (*dbIt).first.stationId 
                 << " has coarse " << (*dbIt).second.nClock
                 << " and phase " << (*dbIt).second.tPhase << std::endl;
       }
       
       string delayRecord = "DTTPGParametersRcd";
       DTCalibDBUtils::writeToDB(delayRecord,delayMap);
       
   }
 
 }

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

protectedvirtual

Perform client diagnostic in online.

Reimplemented from DTLocalTriggerBaseTest.

Definition at line 76 of file DTLocalTriggerSynchTest.cc.

References python.rootplot.argparse::category, edm::EventSetup::get(), and metProducer_cfi::parameters.

                                                                                       {
 
   if (bookingdone) return;
 
   vector<string>::const_iterator iTr   = trigSources.begin();
   vector<string>::const_iterator trEnd = trigSources.end();
   vector<string>::const_iterator iHw   = hwSources.begin();
   vector<string>::const_iterator hwEnd = hwSources.end();
 
   //Booking
   if(parameters.getUntrackedParameter<bool>("staticBooking", true)){
     for (; iTr != trEnd; ++iTr){
       trigSource = (*iTr);
       for (; iHw != hwEnd; ++iHw){
     hwSource = (*iHw);
     std::vector<const DTChamber*>::const_iterator chambIt  = muonGeom->chambers().begin();
     std::vector<const DTChamber*>::const_iterator chambEnd = muonGeom->chambers().end();
     for (; chambIt!=chambEnd; ++chambIt) { 
       DTChamberId chId = ((*chambIt)->id());
       bookChambHistos(ibooker,chId,ratioHistoTag);
     }
       }
     }
   }
 
   LogVerbatim(category()) << "[" << testName << "Test]: book Histograms" << endl;
 
   if (parameters.getParameter<bool>("fineParamDiff")) {
     ESHandle<DTTPGParameters> wPhaseHandle;
     c.get<DTTPGParametersRcd>().get(wPhaseHandle);
     wPhaseMap = (*wPhaseHandle);
   }
 
   bookingdone = 1; 
 
 }

float DTLocalTriggerSynchTest::getFloatFromME	(	DQMStore::IGetter &	igetter,
		DTChamberId	chId,
		std::string	meType
	)

protected

Get float MEs.

Definition at line 233 of file DTLocalTriggerSynchTest.cc.

References python.rootplot.argparse::category, DQMStore::IGetter::get(), MonitorElement::getFloatValue(), DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

                                                                                      {
    
    stringstream wheel; wheel << chId.wheel();
    stringstream station; station << chId.station();
    stringstream sector; sector << chId.sector();
 
    string folderName = topFolder(hwSource=="TM") + "Wheel" +  wheel.str() +
      "/Sector" + sector.str() + "/Station" + station.str() + "/" ; 
 
    string histoname = sourceFolder + folderName 
      + meType
      + "_W" + wheel.str()  
      + "_Sec" + sector.str()
      + "_St" + station.str();
 
    MonitorElement* me = igetter.get(histoname);
    if (me) { 
      return me->getFloatValue(); 
    }
    else { 
      LogProblem(category()) << "[" << testName << "Test]: " << histoname << " is not a valid ME. 0 returned" << std::endl;
    }
    
    return 0;
 
  }

void DTLocalTriggerSynchTest::makeRatioME	(	TH1F *	numerator,
		TH1F *	denominator,
		MonitorElement *	result
	)

protected

Compute efficiency plots.

Definition at line 226 of file DTLocalTriggerSynchTest.cc.

References PostProcessor_cff::efficiency, and MonitorElement::getTH1F().

                                                                                                    {
   
   TH1F* efficiency = result->getTH1F();
   efficiency->Divide(numerator,denominator,1,1,"");
   
 }

void DTLocalTriggerSynchTest::runClientDiagnostic	(	DQMStore::IBooker &	ibooker,
		DQMStore::IGetter &	igetter
	)

protectedvirtual

DQM Client Diagnostic.

Implements DTLocalTriggerBaseTest.

Definition at line 115 of file DTLocalTriggerSynchTest.cc.

References python.rootplot.argparse::category, DQMStore::IGetter::get(), MonitorElement::getName(), and DetId::rawId().

                                                                                                       {
 
   // Loop over Trig & Hw sources
   for (vector<string>::const_iterator iTr = trigSources.begin(); iTr != trigSources.end(); ++iTr){
     trigSource = (*iTr);
     for (vector<string>::const_iterator iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw){
       hwSource = (*iHw);
       std::vector<const DTChamber*>::const_iterator chambIt  = muonGeom->chambers().begin();
       std::vector<const DTChamber*>::const_iterator chambEnd = muonGeom->chambers().end();
       for (; chambIt!=chambEnd; ++chambIt) { 
     DTChamberId chId = (*chambIt)->id();
     uint32_t indexCh = chId.rawId();
 
     // Perform peak finding
 
     TH1F *numH     = getHisto<TH1F>(igetter.get(getMEName(numHistoTag,"", chId)));
     TH1F *denH     = getHisto<TH1F>(igetter.get(getMEName(denHistoTag,"", chId)));
         
     if (numH && denH && numH->GetEntries()>minEntries && denH->GetEntries()>minEntries) {         
       std::map<std::string,MonitorElement*> innerME = chambME[indexCh];
       MonitorElement* ratioH = innerME.find(fullName(ratioHistoTag))->second;
       makeRatioME(numH,denH,ratioH);
       try {
             //Need our own copy to avoid threading problems
         TF1 mypol8("mypol8","pol8");
         getHisto<TH1F>(ratioH)->Fit(&mypol8,"CQO");
       } catch (cms::Exception& iException) {
         edm::LogPrint(category()) << "[" << testName 
                      << "Test]: Error fitting " 
                      << ratioH->getName() << " returned 0" << endl;
       }
     } else { 
       if (!numH || !denH) {
         LogPrint(category()) << "[" << testName 
                 << "Test]: At least one of the required Histograms was not found for chamber " 
                 << chId << ". Peaks not computed" << endl;
       } else {
         LogPrint(category()) << "[" << testName 
                 << "Test]: Number of plots entries for " 
                 << chId << " is less than minEntries=" 
                 << minEntries <<".  Peaks not computed" << endl;
       }
     }
 
       }
     }
   } 
 
 }

Member Data Documentation

bool DTLocalTriggerSynchTest::bookingdone

private

Definition at line 72 of file DTLocalTriggerSynchTest.h.

double DTLocalTriggerSynchTest::bxTime

private

Definition at line 64 of file DTLocalTriggerSynchTest.h.

std::map<uint32_t,std::map<std::string,MonitorElement*> > DTLocalTriggerSynchTest::chambME

private

Definition at line 60 of file DTLocalTriggerSynchTest.h.

std::string DTLocalTriggerSynchTest::denHistoTag

private

Definition at line 62 of file DTLocalTriggerSynchTest.h.

int DTLocalTriggerSynchTest::minEntries

private

Definition at line 68 of file DTLocalTriggerSynchTest.h.

int DTLocalTriggerSynchTest::nBXHigh

private

Definition at line 67 of file DTLocalTriggerSynchTest.h.

int DTLocalTriggerSynchTest::nBXLow

private

Definition at line 66 of file DTLocalTriggerSynchTest.h.

std::string DTLocalTriggerSynchTest::numHistoTag

private

Definition at line 61 of file DTLocalTriggerSynchTest.h.

bool DTLocalTriggerSynchTest::rangeInBX

private

Definition at line 65 of file DTLocalTriggerSynchTest.h.

std::string DTLocalTriggerSynchTest::ratioHistoTag

private

Definition at line 63 of file DTLocalTriggerSynchTest.h.

DTTPGParameters DTLocalTriggerSynchTest::wPhaseMap

private

Definition at line 70 of file DTLocalTriggerSynchTest.h.

bool DTLocalTriggerSynchTest::writeDB

private

Definition at line 69 of file DTLocalTriggerSynchTest.h.

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation