#include <DTDCSByLumiSummary.h>

Inheritance diagram for DTDCSByLumiSummary:

Public Member Functions
	DTDCSByLumiSummary (const edm::ParameterSet &pset)
	Constructor.
virtual	~DTDCSByLumiSummary ()
	Destructor.
Private Member Functions
virtual void	analyze (const edm::Event &event, const edm::EventSetup &setup)
virtual void	beginJob ()
virtual void	beginLuminosityBlock (const edm::LuminosityBlock &lumi, const edm::EventSetup &setup)
virtual void	endJob ()
virtual void	endLuminosityBlock (const edm::LuminosityBlock &lumi, const edm::EventSetup &setup)
virtual void	endRun (const edm::Run &run, const edm::EventSetup &setup)
Private Attributes
std::map< int, std::vector < float > >	dcsFracPerLumi
MonitorElement *	globalHVSummary
std::vector < DTTimeEvolutionHisto * >	hDCSFracTrend
DQMStore *	theDQMStore
MonitorElement *	totalDCSFraction
std::vector< MonitorElement * >	totalDCSFractionWh

Detailed Description

No description available.

Date:: 2011/03/02 14:00:06

Revision:: 1.1

Author:: C. Battilana - CIEMAT; P. Bellan - INFN PD; A. Branca = INFN PD

Definition at line 26 of file DTDCSByLumiSummary.h.

Constructor & Destructor Documentation

DTDCSByLumiSummary::DTDCSByLumiSummary ( const edm::ParameterSet & pset )

Constructor.

Definition at line 30 of file DTDCSByLumiSummary.cc.

{}

DTDCSByLumiSummary::~DTDCSByLumiSummary ( ) [virtual]

Destructor.

Definition at line 33 of file DTDCSByLumiSummary.cc.

{}

Member Function Documentation

void DTDCSByLumiSummary::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)		`[private, virtual]`

Implements edm::EDAnalyzer.

Definition at line 201 of file DTDCSByLumiSummary.cc.

                                                                           { 

}

void DTDCSByLumiSummary::beginJob ( void ) [private, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 36 of file DTDCSByLumiSummary.cc.

References cmsCodeRules::cppFunctionSkipper::operator, and MonitorElement::setLumiFlag().

                                 {

  theDQMStore = Service<DQMStore>().operator->();

  // book the ME
  theDQMStore->setCurrentFolder("DT/EventInfo/DCSContents");

  totalDCSFraction = theDQMStore->bookFloat("DTDCSSummary");  
  totalDCSFraction->setLumiFlag(); // set LumiFlag to DCS content value (save it by lumi)

  globalHVSummary = theDQMStore->book2D("HVGlbSummary","HV Status Summary",1,1,13,5,-2,3);
  globalHVSummary->setAxisTitle("Sectors",1);
  globalHVSummary->setAxisTitle("Wheel",2);

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

    stringstream wheel_str; wheel_str << wh;

    MonitorElement* FractionWh = theDQMStore->bookFloat("DT_Wheel"+wheel_str.str());  
    FractionWh->setLumiFlag(); // set LumiFlag to DCS content value (save it by lumi)

    totalDCSFractionWh.push_back(FractionWh);
  }

  globalHVSummary->Reset();

}

void DTDCSByLumiSummary::beginLuminosityBlock	(	const edm::LuminosityBlock &	lumi,
		const edm::EventSetup &	setup
	)		`[private, virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 65 of file DTDCSByLumiSummary.cc.

                                                                                                   {

  // CB LumiFlag marked products are reset on LS boundaries
  totalDCSFraction->Reset(); 

  for(int wh=-2;wh<=2;wh++){
    totalDCSFractionWh[wh+2]->Reset();
  }

}

void DTDCSByLumiSummary::endJob ( void ) [private, virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 196 of file DTDCSByLumiSummary.cc.

                                {  

}

void DTDCSByLumiSummary::endLuminosityBlock	(	const edm::LuminosityBlock &	lumi,
		const edm::EventSetup &	setup
	)		`[private, virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 77 of file DTDCSByLumiSummary.cc.

References MonitorElement::getBinContent(), edm::LuminosityBlockBase::id(), LogTrace, and edm::LuminosityBlockID::luminosityBlock().

                                                                                                 {

  // Get the by lumi product plot from the task
  int lumiNumber = lumi.id().luminosityBlock();

  bool null_pointer_histo(0);

  std::vector<float> wh_activeFrac;

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

    stringstream wheel_str; wheel_str << wh;    

    string hActiveUnitsPath = "DT/EventInfo/DCSContents/hActiveUnits"+wheel_str.str();
    MonitorElement *hActiveUnits = theDQMStore->get(hActiveUnitsPath);

    if (hActiveUnits) {
      float activeFrac = static_cast<float>(hActiveUnits->getBinContent(2)) /  // CB 2nd bin is # of active channels
        hActiveUnits->getBinContent(1);    // first bin is overall number of channels

      if(activeFrac < 0.) activeFrac=-1;

      wh_activeFrac.push_back(activeFrac);

      // Fill by lumi Certification ME
      totalDCSFraction->Fill(activeFrac); 
      totalDCSFractionWh[wh+2]->Fill(activeFrac);

    } else {
      LogTrace("DTDQM|DTMonitorClient|DTDCSByLumiSummary")
        << "[DTDCSByLumiSummary]: got null pointer retrieving histo at :" 
        << hActiveUnitsPath << " for lumi # " << lumiNumber
        << "client operation not performed." << endl;

      null_pointer_histo=true;
    }    

  } // end loop on wheels

  if(!null_pointer_histo) dcsFracPerLumi[lumiNumber] = wh_activeFrac; // Fill map to be used to compute trend plots

}

void DTDCSByLumiSummary::endRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)		`[private, virtual]`

Reimplemented from edm::EDAnalyzer.

Definition at line 121 of file DTDCSByLumiSummary.cc.

References csvReporter::r.

                                                                             {


  // Book trend plots ME & loop on map to fill it with by lumi info
  map<int,std::vector<float> >::const_iterator fracPerLumiIt  = dcsFracPerLumi.begin();
  map<int,std::vector<float> >::const_iterator fracPerLumiEnd = dcsFracPerLumi.end();

  int fLumi = dcsFracPerLumi.begin()->first;
  int lLumi = dcsFracPerLumi.rbegin()->first;
  theDQMStore->setCurrentFolder("DT/EventInfo/DCSContents");

  int nLumis = lLumi-fLumi + 1.;

  // trend plots
  for(int wh=-2; wh<=2; wh++) {

    stringstream wheel_str; wheel_str << wh;    

    DTTimeEvolutionHisto* trend;

    trend = new DTTimeEvolutionHisto(theDQMStore, "hDCSFracTrendWh" + wheel_str.str(), "Fraction of DT-HV ON Wh" + wheel_str.str(),
        nLumis, fLumi, 1, false, 2);

    hDCSFracTrend.push_back(trend);

  }

  float goodLSperWh[5] = {0,0,0,0,0}; 
  float badLSperWh[5] = {0,0,0,0,0};

  // fill trend plots and save infos for summaryPlot
  for(;fracPerLumiIt!=fracPerLumiEnd;++fracPerLumiIt) {

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

      std::vector<float> activeFracPerWh;
      activeFracPerWh =  fracPerLumiIt->second;

      hDCSFracTrend[wh+2]->setTimeSlotValue(activeFracPerWh[wh+2],fracPerLumiIt->first);

      if( activeFracPerWh[wh+2] > 0 ) { // we do not count the lumi were the DTs are off (no real problem), 
        // even if this can happen in the middle of a run (real problem: to be fixed)
        if( activeFracPerWh[wh+2] > 0.9 ) goodLSperWh[wh+2]++;
        else { 
          badLSperWh[wh+2]++;
        }
      } else {  // there is no HV value OR all channels OFF
        if( activeFracPerWh[wh+2] < 0 ) badLSperWh[wh+2]=-1;       // if there were no HV values, activeFrac returning -1
      }

    }

  }

  // fill summaryPlot
  for(int wh=-2; wh<=2; wh++) {

    if( goodLSperWh[wh+2] != 0 || badLSperWh[wh+2] == -1 ) {

      float r = badLSperWh[wh+2]/fabs(goodLSperWh[wh+2] + badLSperWh[wh+2]);
      if( r > 0.5 ) globalHVSummary->Fill(1,wh,0);
      else globalHVSummary->Fill(1,wh,1); 
      if( r == -1 ) globalHVSummary->Fill(1,wh,-1);    

    } else {

      globalHVSummary->Fill(1,wh,0);

    }

  }

}