DTTriggerEfficiencyTest Class Reference

#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
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

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	makeEfficiencyME2D (TH2F *numerator, TH2F *denominator, MonitorElement *result)
	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="")
	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::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Attributes
std::map< uint32_t, std::map < std::string, MonitorElement * > >	chambME
bool	detailedPlots
DTTrigGeomUtils *	trigGeomUtils

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
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

  Date:

  2009/07/29 11:10:52

  Revision:

  1.1

  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 39 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 49 of file DTTriggerEfficiencyTest.cc.

                                                 {
  
}

Member Function Documentation

void DTTriggerEfficiencyTest::beginJob ( void  )

protectedvirtual

BeginJob.

Reimplemented from DTLocalTriggerBaseTest.

Definition at line 54 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 61 of file DTTriggerEfficiencyTest.cc.

References DTLocalTriggerBaseTest::beginRun(), 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
    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 185 of file DTTriggerEfficiencyTest.cc.

References category(), newFWLiteAna::fullName, LogTrace, edm::max(), edm::min(), RecoTauCommonJetSelections_cfi::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::makeEfficiencyME2D ( TH2F *  numerator, TH2F *  denominator, MonitorElement *  result  )

protected

Compute 2D efficiency plots.

Definition at line 158 of file DTTriggerEfficiencyTest.cc.

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

                                                                                                          {
  
  TH2F* efficiency = result->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 spr::find(), newFWLiteAna::fullName, 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);
      // Loop over the TriggerUnits
      for (int wh=-2; wh<=2; ++wh){
    
    TH2F * TrigEffDenum   = getHisto<TH2F>(dbe->get(getMEName("TrigEffDenum","",wh)));
    TH2F * TrigEffNum     = getHisto<TH2F>(dbe->get(getMEName("TrigEffNum","",wh)));
    TH2F * TrigEffCorrNum = getHisto<TH2F>(dbe->get(getMEName("TrigEffCorrNum","",wh)));
    
    if (TrigEffDenum && TrigEffNum && TrigEffCorrNum && TrigEffDenum->GetEntries()>1) {
      
      if( whME[wh].find(fullName("TrigEffPhi")) == whME[wh].end() ){
        bookWheelHistos(wh,"TrigEffPhi");  
        bookWheelHistos(wh,"TrigEffCorrPhi");  
      }
      std::map<std::string,MonitorElement*> *innerME = &(whME[wh]);
      makeEfficiencyME2D(TrigEffNum,TrigEffDenum,innerME->find(fullName("TrigEffPhi"))->second);
      makeEfficiencyME2D(TrigEffCorrNum,TrigEffDenum,innerME->find(fullName("TrigEffCorrPhi"))->second);
      
    }

    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]);
        makeEfficiencyME2D(TrackPosvsAngleAnyQual,TrackPosvsAngle,innerME->find(fullName("TrigEffPosvsAnglePhi"))->second);
        makeEfficiencyME2D(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 51 of file DTTriggerEfficiencyTest.h.

bool DTTriggerEfficiencyTest::detailedPlots

private

Definition at line 53 of file DTTriggerEfficiencyTest.h.

DTTrigGeomUtils* DTTriggerEfficiencyTest::trigGeomUtils

private

Definition at line 52 of file DTTriggerEfficiencyTest.h.