DTLocalTriggerTest.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author C. Battilana S. Marcellini - INFN Bologna
 */


// This class header
#include "DQM/DTMonitorClient/src/DTLocalTriggerTest.h"

// Framework headers
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "DQMServices/Core/interface/DQMStore.h"

// Geometry
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "Geometry/DTGeometry/interface/DTGeometry.h"

// Root
#include "TF1.h"
#include "TProfile.h"


//C++ headers
#include <iostream>
#include <sstream>

using namespace edm;
using namespace std;


DTLocalTriggerTest::DTLocalTriggerTest(const edm::ParameterSet& ps){

  setConfig(ps,"DTLocalTrigger");
  baseFolderDCC = "DT/03-LocalTrigger-DCC/";
  baseFolderDDU = "DT/04-LocalTrigger-DDU/";
  nMinEvts  = ps.getUntrackedParameter<int>("nEventsCert", 5000);

}


DTLocalTriggerTest::~DTLocalTriggerTest(){

}

void DTLocalTriggerTest::beginJob(){
  
  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","Summaries");
    }
    else {
      bookWheelHistos(wh,"CorrFractionSummary","Summaries");
      bookWheelHistos(wh,"2ndFractionSummary","Summaries");
    }
      }
      if (hwSource=="COM") {
    bookCmsHistos("MatchingSummary","Summaries");
      }
      else {
    bookCmsHistos("CorrFractionSummary");
    bookCmsHistos("2ndFractionSummary");
      }
      if (hwSource=="DCC") {
    bookCmsHistos("TrigGlbSummary","",true);
      }
       
    }   
  }

}


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

}


void DTLocalTriggerTest::runClientDiagnostic() {

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

  fillGlobalSummary();

}

void DTLocalTriggerTest::fillGlobalSummary() {

  float glbPerc[5] = { 1., 0.9, 0.6, 0.3, 0.01 };
  trigSource = "";
  hwSource = "DCC";  

  int nSecReadout = 0;

  for (int wh=-2; wh<=2; ++wh) {
    for (int sect=1; sect<=12; ++sect) {

      float maxErr = 8.;
      int corr   = cmsME.find(fullName("CorrFractionSummary"))->second->getBinContent(sect,wh+3);
      int second = cmsME.find(fullName("2ndFractionSummary"))->second->getBinContent(sect,wh+3);
      int lut=0;
      MonitorElement * lutsME = dbe->get(topFolder(hwSource=="DCC") + "Summaries/TrigLutSummary");
      if (lutsME) {
    lut = lutsME->getBinContent(sect,wh+3);
    maxErr+=4;
      } else {
    LogTrace(category()) << "[" << testName 
     << "Test]: DCC Lut test Summary histo not found." << endl;
      }
      (corr <5 || second<5) && nSecReadout++;
      int errcode = ((corr<5 ? corr : 4) + (second<5 ? second : 4) + (lut<5 ? lut : 4) );
      errcode = min(int((errcode/maxErr + 0.01)*5),5);
      cmsME.find("TrigGlbSummary")->second->setBinContent(sect,wh+3,glbPerc[errcode]);
    
    }
  }

  if (!nSecReadout) 
    cmsME.find("TrigGlbSummary")->second->Reset(); // white histo id DCC is not RO
  
  string nEvtsName = "DT/EventInfo/Counters/nProcessedEventsTrigger";
  MonitorElement * meProcEvts = dbe->get(nEvtsName);

  if (meProcEvts) {
    int nProcEvts = meProcEvts->getFloatValue();
    cmsME.find("TrigGlbSummary")->second->setEntries(nProcEvts < nMinEvts ? 10. : nProcEvts);
  } else {
    cmsME.find("TrigGlbSummary")->second->setEntries(nMinEvts + 1);
    LogVerbatim (category()) << "[" << testName 
     << "Test]: ME: " <<  nEvtsName << " not found!" << endl;
  }

}