#include <DTTriggerEfficiencyTest.h>

Inheritance diagram for DTTriggerEfficiencyTest:

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

virtual	~DTTriggerEfficiencyTest ()
	Destructor. More...

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

virtual	~DTLocalTriggerBaseTest ()
	Destructor. More...

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

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

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
void	beginJob ()
	BeginJob. More...

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

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

void	bookHistos (std::string hTag, std::string folder)
	Book the new MEs (global) More...

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

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

void	makeEfficiencyME (TH2F numerator, TH2F denominator, MonitorElement result2DWh, MonitorElement result1DWh, MonitorElement *result1D)
	Compute 1D/2D efficiency plots. More...

void	makeEfficiencyME (TH2F numerator, TH2F denominator, MonitorElement *result2DWh)
	Compute 2D efficiency plots. More...

void	runClientDiagnostic ()
	DQM Client Diagnostic. More...

Protected Member Functions inherited from DTLocalTriggerBaseTest
void	analyze (const edm::Event &e, const edm::EventSetup &c)
	Analyze. More...

void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
	Perform begin lumiblock operations. More...

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

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

void	bookWheelHistos (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	endJob ()
	EndJob. More...

void	endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
	Perform client diagnostic in online. 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 isDCC)
	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
std::map< uint32_t, std::map < std::string, MonitorElement * > >	chambME

bool	detailedPlots

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

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

DTTrigGeomUtils *	trigGeomUtils

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Attributes inherited from DTLocalTriggerBaseTest
std::string	baseFolderDCC

std::string	baseFolderDDU

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

DQMStore *	dbe

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

Definition at line 20 of file DTTriggerEfficiencyTest.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 37 of file DTTriggerEfficiencyTest.cc.

References edm::ParameterSet::getUntrackedParameter().

                                                                          {
 
   setConfig(ps,"DTTriggerEfficiency");
   baseFolderDCC = "DT/03-LocalTrigger-DCC/";
   baseFolderDDU = "DT/04-LocalTrigger-DDU/";
   detailedPlots = ps.getUntrackedParameter<bool>("detailedAnalysis",true);
 
 }

DTTriggerEfficiencyTest::~DTTriggerEfficiencyTest ( )

virtual

Destructor.

Definition at line 47 of file DTTriggerEfficiencyTest.cc.

                                                  {
 
 }

Member Function Documentation

void DTTriggerEfficiencyTest::beginJob ( void )

protectedvirtual

BeginJob.

Reimplemented from DTLocalTriggerBaseTest.

Definition at line 52 of file DTTriggerEfficiencyTest.cc.

References DTLocalTriggerBaseTest::beginJob().

                                       {
 
   DTLocalTriggerBaseTest::beginJob();
 
 }

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

protectedvirtual

BeginRun.

Reimplemented from DTLocalTriggerBaseTest.

Definition at line 59 of file DTTriggerEfficiencyTest.cc.

References DTLocalTriggerBaseTest::beginRun(), bookHistos(), and Parameters::parameters.

                                                                             {
 
   DTLocalTriggerBaseTest::beginRun(r,c);
   trigGeomUtils = new DTTrigGeomUtils(muonGeom);
 
   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);
         // Loop over the TriggerUnits
         bookHistos("TrigEffPhi","");
         bookHistos("TrigEffCorrPhi","");
         for (int wh=-2; wh<=2; ++wh){
           if (detailedPlots) {
             for (int sect=1; sect<=12; ++sect){
               for (int stat=1; stat<=4; ++stat){
                 DTChamberId chId(wh,stat,sect);
                 bookChambHistos(chId,"TrigEffPosvsAnglePhi","Segment");
                 bookChambHistos(chId,"TrigEffPosvsAngleCorrPhi","Segment");
               }
             }
           }
           bookWheelHistos(wh,"TrigEffPhi","");  
           bookWheelHistos(wh,"TrigEffCorrPhi","");  
         }
       }
     }
   }
 
 }

void DTTriggerEfficiencyTest::bookChambHistos	(	DTChamberId	chambId,
		std::string	htype,
		std::string	folder = `""`
	)

protected

Book the new MEs (for each chamber)

Definition at line 313 of file DTTriggerEfficiencyTest.cc.

References python.rootplot.argparse::category, newFWLiteAna::fullName, LogTrace, edm::max(), edm::min(), pileupCalc::nbins, DetId::rawId(), DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, and DTChamberId::wheel().

                                                                                               {
 
   stringstream wheel; wheel << chambId.wheel();
   stringstream station; station << chambId.station();   
   stringstream sector; sector << chambId.sector();
 
   string fullType  = fullName(htype);
   bool isDCC = hwSource=="DCC" ;
   string HistoName = fullType + "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
 
   dbe->setCurrentFolder(topFolder(isDCC) + 
       "Wheel" + wheel.str() +
       "/Sector" + sector.str() +
       "/Station" + station.str() + 
       "/" + folder + "/");
 
   LogTrace(category()) << "[" << testName << "Test]: booking " + topFolder(isDCC) + "Wheel" << wheel.str() 
     <<"/Sector" << sector.str() << "/Station" << station.str() << "/" + folder + "/" << HistoName;
 
 
   uint32_t indexChId = chambId.rawId();
   float min, max;
   int nbins;
   trigGeomUtils->phiRange(chambId,min,max,nbins,20);
   if (htype.find("TrigEffPosvsAnglePhi") == 0 ){
     chambME[indexChId][fullType] = dbe->book2D(HistoName.c_str(),"Trigger efficiency (any qual.) position vs angle (Phi)",12,-30.,30.,nbins,min,max);
     return;
   }
   if (htype.find("TrigEffPosvsAngleCorrPhi") == 0 ){
     chambME[indexChId][fullType] = dbe->book2D(HistoName.c_str(),"Trigger efficiency (correlated) pos vs angle (Phi)",12,-30.,30.,nbins,min,max);
     return;
   }
 
 }

void DTTriggerEfficiencyTest::bookHistos	(	std::string	hTag,
		std::string	folder
	)

protected

Book the new MEs (global)

Definition at line 239 of file DTTriggerEfficiencyTest.cc.

References newFWLiteAna::fullName.

                                                                   {
 
   string basedir;  
   bool isDCC = hwSource=="DCC" ;  
   basedir = topFolder(isDCC);   //Book summary histo outside Task directory 
 
   if (folder != "") {
     basedir += folder +"/" ;
   }
   dbe->setCurrentFolder(basedir);
 
   string fullTag = fullName(hTag);
   string hname = fullTag + "_All";
 
   globalEffDistr[fullTag] = dbe->book1D(hname.c_str(),hname.c_str(),51,0.,1.02);
   globalEffDistr[fullTag] ->setAxisTitle("Trig Eff",1);
 
 }

void DTTriggerEfficiencyTest::bookWheelHistos	(	int	wheel,
		std::string	hTag,
		std::string	folder
	)

protected

Book the new MEs (for each wheel)

Definition at line 258 of file DTTriggerEfficiencyTest.cc.

References python.rootplot.argparse::category, newFWLiteAna::fullName, LogTrace, MonitorElement::setAxisTitle(), and MonitorElement::setBinLabel().

                                                                                  {
 
   stringstream wh; wh << wheel;
   string basedir;  
   bool isDCC = hwSource=="DCC" ;  
   if (hTag.find("Summary") != string::npos) {
     basedir = topFolder(isDCC);   //Book summary histo outside wheel directories
   } else {
     basedir = topFolder(isDCC) + "Wheel" + wh.str() + "/" ;
 
   }
   if (folder != "") {
     basedir += folder +"/" ;
   }
   dbe->setCurrentFolder(basedir);
 
   string fullTag = fullName(hTag);
   string hname    = fullTag+ "_W" + wh.str();
 
   string hnameAll = fullTag+ "_All_W" + wh.str();
 
   LogTrace(category()) << "[" << testName << "Test]: booking "<< basedir << hname;
 
   (EffDistrPerWh[wheel])[fullTag] = dbe->book1D(hnameAll.c_str(),hnameAll.c_str(),51,0.,1.02);
 
   if (hTag.find("Phi")!= string::npos ||
       hTag.find("Summary") != string::npos ){    
     MonitorElement* me = dbe->book2D(hname.c_str(),hname.c_str(),12,1,13,4,1,5);
 
     //     setLabelPh(me);
     me->setBinLabel(1,"MB1",2);
     me->setBinLabel(2,"MB2",2);
     me->setBinLabel(3,"MB3",2);
     me->setBinLabel(4,"MB4",2);
     me->setAxisTitle("Sector",1);
 
     whME[wheel][fullTag] = me;
     return;
   }
 
   if (hTag.find("Theta") != string::npos){
     MonitorElement* me =dbe->book2D(hname.c_str(),hname.c_str(),12,1,13,3,1,4);
 
     //     setLabelTh(me);
     me->setBinLabel(1,"MB1",2);
     me->setBinLabel(2,"MB2",2);
     me->setBinLabel(3,"MB3",2);
     me->setAxisTitle("Sector",1);
 
     whME[wheel][fullTag] = me;
     return;
   }
 
 }

string DTTriggerEfficiencyTest::getMEName	(	std::string	histoTag,
		std::string	folder,
		int	wh
	)

protected

Get the ME name (by wheel)

Definition at line 226 of file DTTriggerEfficiencyTest.cc.

References newFWLiteAna::fullName.

                                                                                 {
 
   stringstream wheel; wheel << wh;
 
   string folderName =  topFolder(hwSource=="DCC") + folder + "/";
 
   string histoname = sourceFolder + folderName 
     + fullName(histoTag) + "_W" + wheel.str();
 
   return histoname;
 
 }

void DTTriggerEfficiencyTest::makeEfficiencyME	(	TH2F *	numerator,
		TH2F *	denominator,
		MonitorElement *	result2DWh,
		MonitorElement *	result1DWh,
		MonitorElement *	result1D
	)

protected

Compute 1D/2D efficiency plots.

Definition at line 171 of file DTTriggerEfficiencyTest.cc.

References postValidation_cfi::efficiency, relativeConstraints::error, MonitorElement::Fill(), MonitorElement::getTH2F(), and mathSSE::sqrt().

                                                                                                                                                                   {
 
   TH2F* efficiency = result2DWh->getTH2F();
   efficiency->Divide(numerator,denominator,1,1,"");
 
   int nbinsx = efficiency->GetNbinsX();
   int nbinsy = efficiency->GetNbinsY();
   for (int binx=1; binx<=nbinsx; ++binx){
     for (int biny=1; biny<=nbinsy; ++biny){
       float error = 0;
       float bineff = efficiency->GetBinContent(binx,biny);
 
       result1DWh->Fill(bineff);
       result1D->Fill(bineff);
 
       if (denominator->GetBinContent(binx,biny)){
         error = sqrt(bineff*(1-bineff)/denominator->GetBinContent(binx,biny));
       }
       else {
         error = 1;
         efficiency->SetBinContent(binx,biny,0.);
       }
 
       efficiency->SetBinError(binx,biny,error);
     }
   }
 
 }

void DTTriggerEfficiencyTest::makeEfficiencyME	(	TH2F *	numerator,
		TH2F *	denominator,
		MonitorElement *	result2DWh
	)

protected

Compute 2D efficiency plots.

Definition at line 200 of file DTTriggerEfficiencyTest.cc.

References postValidation_cfi::efficiency, relativeConstraints::error, MonitorElement::getTH2F(), and mathSSE::sqrt().

                                                                                                             {
 
   TH2F* efficiency = result2DWh->getTH2F();
   efficiency->Divide(numerator,denominator,1,1,"");
 
   int nbinsx = efficiency->GetNbinsX();
   int nbinsy = efficiency->GetNbinsY();
   for (int binx=1; binx<=nbinsx; ++binx){
     for (int biny=1; biny<=nbinsy; ++biny){
       float error = 0;
       float bineff = efficiency->GetBinContent(binx,biny);
 
       if (denominator->GetBinContent(binx,biny)){
         error = sqrt(bineff*(1-bineff)/denominator->GetBinContent(binx,biny));
       }
       else {
         error = 1;
         efficiency->SetBinContent(binx,biny,0.);
       }
 
       efficiency->SetBinError(binx,biny,error);
     }
   }
 
 }    

void DTTriggerEfficiencyTest::runClientDiagnostic ( )

protectedvirtual

DQM Client Diagnostic.

Implements DTLocalTriggerBaseTest.

Definition at line 99 of file DTTriggerEfficiencyTest.cc.

References bookHistos(), spr::find(), newFWLiteAna::fullName, DetId::rawId(), and edm::second().

                                                   {
 
   // 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);
       // Loop over the TriggerUnits
       if( globalEffDistr.find(fullName("TrigEffPhi")) == globalEffDistr.end() ){
         bookHistos("TrigEffPhi","");
         bookHistos("TrigEffCorrPhi","");
       }
       for (int wh=-2; wh<=2; ++wh){
 
         TH2F * TrigEffDenum   = getHisto<TH2F>(dbe->get(getMEName("TrigEffDenum","Task",wh)));
         TH2F * TrigEffNum     = getHisto<TH2F>(dbe->get(getMEName("TrigEffNum","Task",wh)));
         TH2F * TrigEffCorrNum = getHisto<TH2F>(dbe->get(getMEName("TrigEffCorrNum","Task",wh)));
 
         if (TrigEffDenum && TrigEffNum && TrigEffCorrNum && TrigEffDenum->GetEntries()>1) {
 
           if( whME[wh].find(fullName("TrigEffPhi")) == whME[wh].end() ){
             bookWheelHistos(wh,"TrigEffPhi","");  
             bookWheelHistos(wh,"TrigEffCorrPhi","");  
           }
 
           MonitorElement* Eff1DAll_TrigEffPhi = (&globalEffDistr)->find(fullName("TrigEffPhi"))->second;
           MonitorElement* Eff1DAll_TrigEffCorrPhi = (&globalEffDistr)->find(fullName("TrigEffCorrPhi"))->second;
 
           MonitorElement* Eff1DWh_TrigEffPhi = (&(EffDistrPerWh[wh]))->find(fullName("TrigEffPhi"))->second;
           MonitorElement* Eff1DWh_TrigEffCorrPhi = (&(EffDistrPerWh[wh]))->find(fullName("TrigEffCorrPhi"))->second;
 
           MonitorElement* Eff2DWh_TrigEffPhi = (&(whME[wh]))->find(fullName("TrigEffPhi"))->second;
           MonitorElement* Eff2DWh_TrigEffCorrPhi = (&(whME[wh]))->find(fullName("TrigEffCorrPhi"))->second;
 
           makeEfficiencyME(TrigEffNum,TrigEffDenum,Eff2DWh_TrigEffPhi,Eff1DWh_TrigEffPhi,Eff1DAll_TrigEffPhi);
           makeEfficiencyME(TrigEffCorrNum,TrigEffDenum,Eff2DWh_TrigEffCorrPhi,Eff1DWh_TrigEffCorrPhi,Eff1DAll_TrigEffCorrPhi);
 
         }
 
         if (detailedPlots) {
           for (int stat=1; stat<=4; ++stat){
             for (int sect=1; sect<=12; ++sect){
               DTChamberId chId(wh,stat,sect);
               uint32_t indexCh = chId.rawId();
 
               // Perform Efficiency analysis (Phi+Segments 2D)
               TH2F * TrackPosvsAngle        = getHisto<TH2F>(dbe->get(getMEName("TrackPosvsAngle","Segment", chId)));
               TH2F * TrackPosvsAngleAnyQual = getHisto<TH2F>(dbe->get(getMEName("TrackPosvsAngleAnyQual","Segment", chId)));
               TH2F * TrackPosvsAngleCorr    = getHisto<TH2F>(dbe->get(getMEName("TrackPosvsAngleCorr","Segment", chId)));
 
               if (TrackPosvsAngle && TrackPosvsAngleAnyQual && TrackPosvsAngleCorr && TrackPosvsAngle->GetEntries()>1) {
 
                 if( chambME[indexCh].find(fullName("TrigEffAnglePhi")) == chambME[indexCh].end()){
                   bookChambHistos(chId,"TrigEffPosvsAnglePhi","Segment");
                   bookChambHistos(chId,"TrigEffPosvsAngleCorrPhi","Segment");
                 }
 
                 std::map<std::string,MonitorElement*> *innerME = &(chambME[indexCh]);
                 makeEfficiencyME(TrackPosvsAngleAnyQual,TrackPosvsAngle,innerME->find(fullName("TrigEffPosvsAnglePhi"))->second);
                 makeEfficiencyME(TrackPosvsAngleCorr,TrackPosvsAngle,innerME->find(fullName("TrigEffPosvsAngleCorrPhi"))->second);
 
               }
             }
           }
         }
       }
 
     }
   } 
 
 }

Member Data Documentation

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

private

Definition at line 65 of file DTTriggerEfficiencyTest.h.

bool DTTriggerEfficiencyTest::detailedPlots

private

Definition at line 67 of file DTTriggerEfficiencyTest.h.

std::map<int,std::map<std::string,MonitorElement*> > DTTriggerEfficiencyTest::EffDistrPerWh

private

Definition at line 64 of file DTTriggerEfficiencyTest.h.

std::map<std::string, MonitorElement*> DTTriggerEfficiencyTest::globalEffDistr

private

Definition at line 63 of file DTTriggerEfficiencyTest.h.

DTTrigGeomUtils* DTTriggerEfficiencyTest::trigGeomUtils

private

Definition at line 66 of file DTTriggerEfficiencyTest.h.

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation