DTLocalTriggerTest Class Reference

#include <DTLocalTriggerTest.h>

Inheritance diagram for DTLocalTriggerTest:

Public Member Functions
	DTLocalTriggerTest (const edm::ParameterSet &ps)
	Constructor. More...
virtual	~DTLocalTriggerTest ()
	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	runClientDiagnostic ()
	Run client analysis. 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

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:

  2008/10/07 14:26:43

  Revision:

  1.13

  Author

  C. Battilana S. Marcellini - INFN Bologna

Definition at line 20 of file DTLocalTriggerTest.h.

Constructor & Destructor Documentation

DTLocalTriggerTest::DTLocalTriggerTest

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 36 of file DTLocalTriggerTest.cc.

                                                               {

  setConfig(ps,"DTLocalTrigger");
  baseFolderDCC = "DT/03-LocalTrigger-DCC/";
  baseFolderDDU = "DT/04-LocalTrigger-DDU/";

}

DTLocalTriggerTest::~DTLocalTriggerTest ( )

virtual

Destructor.

Definition at line 45 of file DTLocalTriggerTest.cc.

                                       {

}

Member Function Documentation

void DTLocalTriggerTest::beginJob ( void  )

protectedvirtual

BeginJob.

Reimplemented from DTLocalTriggerBaseTest.

Definition at line 49 of file DTLocalTriggerTest.cc.

References DTLocalTriggerBaseTest::beginJob(), and Parameters::parameters.

                                 {
  
  DTLocalTriggerBaseTest::beginJob();


  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 (hwSource=="COM") {
        bookWheelHistos(wh,"MatchingPhi");
      } 
      else { 
        for (int sect=1; sect<=12; ++sect){
          bookSectorHistos(wh,sect,"BXDistribPhi");
          bookSectorHistos(wh,sect,"QualDistribPhi");
        }
        bookWheelHistos(wh,"CorrectBXPhi");
        bookWheelHistos(wh,"ResidualBXPhi");
        bookWheelHistos(wh,"CorrFractionPhi");
        bookWheelHistos(wh,"2ndFractionPhi");
        bookWheelHistos(wh,"TriggerInclusivePhi");
        bookWheelHistos(wh,"CorrectBXTheta");
        if (hwSource=="DDU") {
          bookWheelHistos(wh,"HFractionTheta");
        }
      }
    }
      }
    }
  }
  // Summary test histo booking (only static)
  for (iTr = trigSources.begin(); iTr != trEnd; ++iTr){
    trigSource = (*iTr);
    for (iHw = hwSources.begin(); iHw != hwSources.end(); ++iHw){
      hwSource = (*iHw);
      // Loop over the TriggerUnits
      for (int wh=-2; wh<=2; ++wh){
    if (hwSource=="COM") {
      bookWheelHistos(wh,"MatchingSummary","Matching");
    }
    else {
      bookWheelHistos(wh,"CorrFractionSummary");
      bookWheelHistos(wh,"2ndFractionSummary");
    }
      }
      if (hwSource=="COM") {
    bookCmsHistos("MatchingSummary","Matching");
      }
      else {
    bookCmsHistos("CorrFractionSummary");
    bookCmsHistos("2ndFractionSummary");
      } 
    }   
  }

}

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

protectedvirtual

BeginRun.

Reimplemented from DTLocalTriggerBaseTest.

Definition at line 118 of file DTLocalTriggerTest.cc.

References DTLocalTriggerBaseTest::beginRun().

                                                                        {
  
  DTLocalTriggerBaseTest::beginRun(r,c);

}

void DTLocalTriggerTest::runClientDiagnostic ( )

protectedvirtual

Run client analysis.

Implements DTLocalTriggerBaseTest.

Definition at line 125 of file DTLocalTriggerTest.cc.

References end, python.tagInventory::entries, spr::find(), newFWLiteAna::fullName, and Parameters::parameters.

                                             {

  // 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 stat=1; stat<=4; ++stat){
    for (int wh=-2; wh<=2; ++wh){
      for (int sect=1; sect<=12; ++sect){
        DTChamberId chId(wh,stat,sect);
        int sector_id = (wh+3)+(sect-1)*5;
        // uint32_t indexCh = chId.rawId();
        
        if (hwSource=="COM") {
          // Perform DCC-DDU matching test and generates summaries (Phi view)
          TH2F * DDUvsDCC = getHisto<TH2F>(dbe->get(getMEName("QualDDUvsQualDCC","LocalTriggerPhi", chId)));
          if (DDUvsDCC) {
        
        int matchSummary   = 1;
        
        if (DDUvsDCC->GetEntries()>1) {
          
          double entries     = DDUvsDCC->GetEntries();
          double corrEntries = 0;
          for (int ibin=2; ibin<=8; ++ibin) {
            corrEntries += DDUvsDCC->GetBinContent(ibin,ibin);
          }
          double corrRatio   = corrEntries/entries;
          
          if (corrRatio < parameters.getUntrackedParameter<double>("matchingFracError",.65)){
            matchSummary = 2;
          }
          else if (corrRatio < parameters.getUntrackedParameter<double>("matchingFracWarning",.85)){
            matchSummary = 3;
          }
          else {
            matchSummary = 0;
          }
          
          if( whME[wh].find(fullName("MatchingPhi")) == whME[wh].end() ){
            bookWheelHistos(wh,"MatchingPhi");
          }
          
          whME[wh].find(fullName("MatchingPhi"))->second->setBinContent(sect,stat,corrRatio);
          
        }
        
        whME[wh].find(fullName("MatchingSummary"))->second->setBinContent(sect,stat,matchSummary);

          }
        }
        else {
          // Perform DCC/DDU common plot analysis (Phi ones)
          TH2F * BXvsQual      = getHisto<TH2F>(dbe->get(getMEName("BXvsQual","LocalTriggerPhi", chId)));
          TH1F * BestQual      = getHisto<TH1F>(dbe->get(getMEName("BestQual","LocalTriggerPhi", chId)));
          TH2F * Flag1stvsQual = getHisto<TH2F>(dbe->get(getMEName("Flag1stvsQual","LocalTriggerPhi", chId))); 
          if (BXvsQual && Flag1stvsQual && BestQual) {

        int corrSummary   = 1;
        int secondSummary = 1;
        
        if (BestQual->GetEntries()>1) {
          
          TH1D* BXHH    = BXvsQual->ProjectionY("",6,7,"");
          TH1D* Flag1st = Flag1stvsQual->ProjectionY();
          int BXOK_bin  = BXHH->GetEntries()>=1 ? BXHH->GetMaximumBin() : 51;
          double BXMean = BXHH->GetEntries()>=1 ? BXHH->GetMean() : 51;
          double BX_OK  = BXvsQual->GetYaxis()->GetBinCenter(BXOK_bin);
          double trigsFlag2nd = Flag1st->GetBinContent(2);
          double trigs = Flag1st->GetEntries();
          double besttrigs = BestQual->GetEntries();
          double besttrigsCorr = BestQual->Integral(5,7,"");
          delete BXHH;
          delete Flag1st;
          
          double corrFrac   = besttrigsCorr/besttrigs;
          double secondFrac = trigsFlag2nd/trigs;
          if (corrFrac < parameters.getUntrackedParameter<double>("corrFracError",.5)){
            corrSummary = 2;
          }
          else if (corrFrac < parameters.getUntrackedParameter<double>("corrFracWarning",.6)){
            corrSummary = 3;
          }
          else {
            corrSummary = 0;
          }
          if (secondFrac > parameters.getUntrackedParameter<double>("secondFracError",.2)){
            secondSummary = 2;
          }
          else if (secondFrac > parameters.getUntrackedParameter<double>("secondFracWarning",.1)){
            secondSummary = 3;
          }
          else {
            secondSummary = 0;
          }
          
          if( secME[sector_id].find(fullName("BXDistribPhi")) == secME[sector_id].end() ){
            bookSectorHistos(wh,sect,"QualDistribPhi");
            bookSectorHistos(wh,sect,"BXDistribPhi");
          }

          TH1D* BXDistr   = BXvsQual->ProjectionY();
          TH1D* QualDistr = BXvsQual->ProjectionX();
          std::map<std::string,MonitorElement*> *innerME = &(secME[sector_id]);
          
          int nbinsBX        = BXDistr->GetNbinsX();
          int firstBinCenter = static_cast<int>(BXDistr->GetBinCenter(1));
          int lastBinCenter  = static_cast<int>(BXDistr->GetBinCenter(nbinsBX));
          int iMin = firstBinCenter>-4 ? firstBinCenter : -4;
          int iMax = lastBinCenter<20  ? lastBinCenter  : 20;
          for (int ibin=iMin+5;ibin<=iMax+5; ++ibin) {
            innerME->find(fullName("BXDistribPhi"))->second->setBinContent(ibin,stat,BXDistr->GetBinContent(ibin-5-firstBinCenter+1));
          }
          for (int ibin=1;ibin<=7;++ibin) {
            innerME->find(fullName("QualDistribPhi"))->second->setBinContent(ibin,stat,QualDistr->GetBinContent(ibin));
          }

          delete BXDistr;
          delete QualDistr;

          if( whME[wh].find(fullName("CorrectBXPhi")) == whME[wh].end() ){
            bookWheelHistos(wh,"ResidualBXPhi");
            bookWheelHistos(wh,"CorrectBXPhi");
            bookWheelHistos(wh,"CorrFractionPhi");
            bookWheelHistos(wh,"2ndFractionPhi");
            bookWheelHistos(wh,"TriggerInclusivePhi");
          }
          
          innerME = &(whME[wh]);
          innerME->find(fullName("CorrectBXPhi"))->second->setBinContent(sect,stat,BX_OK+0.00001);
          innerME->find(fullName("ResidualBXPhi"))->second->setBinContent(sect,stat,round(25.*(BXMean-BX_OK))+0.00001);
          innerME->find(fullName("CorrFractionPhi"))->second->setBinContent(sect,stat,corrFrac);
          innerME->find(fullName("TriggerInclusivePhi"))->second->setBinContent(sect,stat,besttrigs);
          innerME->find(fullName("2ndFractionPhi"))->second->setBinContent(sect,stat,secondFrac);
          
        }

        whME[wh].find(fullName("CorrFractionSummary"))->second->setBinContent(sect,stat,corrSummary);
        whME[wh].find(fullName("2ndFractionSummary"))->second->setBinContent(sect,stat,secondSummary);

          }

          if (hwSource=="DDU") {
        // Perform DDU plot analysis (Theta ones)       
        TH2F * ThetaBXvsQual = getHisto<TH2F>(dbe->get(getMEName("ThetaBXvsQual","LocalTriggerTheta", chId)));
        TH1F * ThetaBestQual = getHisto<TH1F>(dbe->get(getMEName("ThetaBestQual","LocalTriggerTheta", chId)));
    
        // no theta triggers in stat 4!
        if (ThetaBXvsQual && ThetaBestQual && stat<4 && ThetaBestQual->GetEntries()>1) {
          TH1D* BXH       = ThetaBXvsQual->ProjectionY("",4,4,"");
          int    BXOK_bin = BXH->GetEffectiveEntries()>=1 ? BXH->GetMaximumBin(): 10;
          double BX_OK    = ThetaBXvsQual->GetYaxis()->GetBinCenter(BXOK_bin);
          double trigs    = ThetaBestQual->GetEntries(); 
          double trigsH   = ThetaBestQual->GetBinContent(4);
          delete BXH; 
        
          // if( secME[sector_id].find(fullName("HFractionTheta")) == secME[sector_id].end() ){
          //        // bookSectorHistos(wh,sect,"CorrectBXTheta");
          //        bookSectorHistos(wh,sect,"HFractionTheta");
          //          }
          //std::map<std::string,MonitorElement*> *innerME = &(secME.find(sector_id)->second);
          // innerME->find(fullName("CorrectBXTheta"))->second->setBinContent(stat,BX_OK);
          //innerME->find(fullName("HFractionTheta"))->second->setBinContent(stat,trigsH/trigs);
        
          if( whME[wh].find(fullName("HFractionTheta")) == whME[wh].end() ){
            bookWheelHistos(wh,"CorrectBXTheta");
            bookWheelHistos(wh,"HFractionTheta");
          }
          std::map<std::string,MonitorElement*> *innerME = &(whME.find(wh)->second);
          innerME->find(fullName("CorrectBXTheta"))->second->setBinContent(sect,stat,BX_OK+0.00001);
          innerME->find(fullName("HFractionTheta"))->second->setBinContent(sect,stat,trigsH/trigs);
        
        }
          }
          else if (hwSource=="DCC") {
        // Perform DCC plot analysis (Theta ones)       
        TH2F * ThetaPosvsBX = getHisto<TH2F>(dbe->get(getMEName("PositionvsBX","LocalTriggerTheta", chId)));
          
        // no theta triggers in stat 4!
        if (ThetaPosvsBX && stat<4 && ThetaPosvsBX->GetEntries()>1) {
          TH1D* BX        = ThetaPosvsBX->ProjectionX();
          int    BXOK_bin = BX->GetEffectiveEntries()>=1 ? BX->GetMaximumBin(): 10;
          double BX_OK    = ThetaPosvsBX->GetXaxis()->GetBinCenter(BXOK_bin);
          delete BX; 
        
          if( whME[wh].find(fullName("CorrectBXTheta")) == whME[wh].end() ){
            bookWheelHistos(wh,"CorrectBXTheta");
          }
          std::map<std::string,MonitorElement*> *innerME = &(whME.find(wh)->second);
          innerME->find(fullName("CorrectBXTheta"))->second->setBinContent(sect,stat,BX_OK+0.00001);
        
        }
          }
        }

      }
    }
      }
    }
  } 

  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);  
      for (int wh=-2; wh<=2; ++wh){
    std::map<std::string,MonitorElement*> *innerME = &(whME[wh]);
    if(hwSource=="COM") {
      TH2F* matchWhSummary   = getHisto<TH2F>(innerME->find(fullName("MatchingSummary"))->second);
      for (int sect=1; sect<=12; ++sect){
        int matchErr      = 0;
        int matchNoData   = 0;
        for (int stat=1; stat<=4; ++stat){
          switch (static_cast<int>(matchWhSummary->GetBinContent(sect,stat))) {
          case 1:
        matchNoData++;
          case 2:
        matchErr++;
          }
        }
        if (matchNoData == 4)   matchErr   = 5;
        cmsME.find(fullName("MatchingSummary"))->second->setBinContent(sect,wh+3,matchErr);
      }
    }
    else {
      TH2F* corrWhSummary   = getHisto<TH2F>(innerME->find(fullName("CorrFractionSummary"))->second);
      TH2F* secondWhSummary = getHisto<TH2F>(innerME->find(fullName("2ndFractionSummary"))->second);
      for (int sect=1; sect<=12; ++sect){
        int corrErr      = 0;
        int secondErr    = 0;
        int corrNoData   = 0;
        int secondNoData = 0;
        for (int stat=1; stat<=4; ++stat){
          switch (static_cast<int>(corrWhSummary->GetBinContent(sect,stat))) {
          case 1:
        corrNoData++;
          case 2:
        corrErr++;
          }
          switch (static_cast<int>(secondWhSummary->GetBinContent(sect,stat))) {
          case 1:
        secondNoData++;
          case 2:
        secondErr++;
          }
        }
        if (corrNoData == 4)   corrErr   = 5;
        if (secondNoData == 4) secondErr = 5;
        cmsME.find(fullName("CorrFractionSummary"))->second->setBinContent(sect,wh+3,corrErr);
        cmsME.find(fullName("2ndFractionSummary"))->second->setBinContent(sect,wh+3,secondErr);
      }
    }
      }
    }
  }

}