PhotonDataCertification Class Reference

#include <PhotonDataCertification.h>

Inheritance diagram for PhotonDataCertification:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	beginRun (const edm::Run &, const edm::EventSetup &)
virtual void	endJob ()
virtual void	endRun (const edm::Run &, const edm::EventSetup &)
	PhotonDataCertification (const edm::ParameterSet &pset)
	~PhotonDataCertification ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Attributes
DQMStore *	dbe_
edm::ParameterSet	parameters_
bool	verbose_

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 Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 49 of file PhotonDataCertification.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 19 of file PhotonDataCertification.cc.

References gather_cfg::cout, dbe_, cmsCodeRules.cppFunctionSkipper::operator, MultipleCompare::pset, and DQMStore::setVerbose().

{

  dbe_ = 0;
  dbe_ = edm::Service<DQMStore>().operator->();
  dbe_->setVerbose(0);
  parameters_ = pset;
  verbose_ = parameters_.getParameter<bool>("verbose");

  if(verbose_) cout << ">>> Constructor (PhotonDataCertification) <<<" << endl;

}

PhotonDataCertification::~PhotonDataCertification

(

)

Definition at line 34 of file PhotonDataCertification.cc.

{
}

Member Function Documentation

void PhotonDataCertification::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::EDAnalyzer.

Definition at line 45 of file PhotonDataCertification.cc.

References gather_cfg::cout.

{
   using namespace edm;

   if(verbose_) std::cout << ">>> Analyze (PhotonDataCertification) <<<" << std::endl;

}

void PhotonDataCertification::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 55 of file PhotonDataCertification.cc.

{
}

void PhotonDataCertification::beginRun ( const edm::Run &  run, const edm::EventSetup &  c  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 67 of file PhotonDataCertification.cc.

References gather_cfg::cout, and edm::RunBase::id().

{

  if(verbose_) std::cout << ">>> BeginRun (PhotonDataCertification) <<<" << std::endl;
  if(verbose_) std::cout << ">>> "<< run.id() << std::endl;

}

void PhotonDataCertification::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 61 of file PhotonDataCertification.cc.

{
}

void PhotonDataCertification::endRun ( const edm::Run &  run, const edm::EventSetup &  c  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 77 of file PhotonDataCertification.cc.

References DQMStore::book2D(), DQMStore::bookFloat(), gather_cfg::cout, dbe_, MonitorElement::Fill(), DQMStore::get(), MonitorElement::getFloatValue(), MonitorElement::getName(), MonitorElement::getQReports(), MonitorElement::getTH2F(), i, and DQMStore::setCurrentFolder().

{
  if(verbose_) std::cout << ">>> EndRun (PhotonDataCertification) <<<" << std::endl;

  std::vector<string> histoNameVector;


  //booking histograms according to naming conventions

  if(!dbe_) return;

  dbe_->setCurrentFolder("Egamma/EventInfo/");

  MonitorElement*  reportSummary = dbe_->bookFloat("reportSummary");
  MonitorElement*  CertificationSummary = dbe_->bookFloat("CertificationSummary");

  //for now these will hold values from eta/phi tests for spikes/holes
  MonitorElement*  reportSummaryMap = dbe_->book2D("reportSummaryMap","reportSummaryMap",2,0,2,2,0,2);
  MonitorElement*  CertificationSummaryMap = dbe_->book2D("CertificationSummaryMap","CertificationSummaryMap",2,0,2,2,0,2);

  TH2 * reportSummaryMapTH2 = reportSummaryMap->getTH2F();

  reportSummaryMapTH2->GetXaxis()->SetBinLabel(1,"Eta");
  reportSummaryMapTH2->GetXaxis()->SetBinLabel(2,"Phi");
  reportSummaryMapTH2->GetYaxis()->SetBinLabel(1,"SpikeTest");
  reportSummaryMapTH2->GetYaxis()->SetBinLabel(2,"HoleTest");

  TH2 * CertificationSummaryMapTH2 = CertificationSummaryMap->getTH2F();

  CertificationSummaryMapTH2->GetXaxis()->SetBinLabel(1,"Eta");
  CertificationSummaryMapTH2->GetXaxis()->SetBinLabel(2,"Phi");
  CertificationSummaryMapTH2->GetYaxis()->SetBinLabel(1,"SpikeTest");
  CertificationSummaryMapTH2->GetYaxis()->SetBinLabel(2,"HoleTest");



  //histograms to be tested and results put into Egamma/EventInfo/reportSummaryContents
  histoNameVector.push_back("nPhoAllEcal");
  histoNameVector.push_back("r9AllEcal");
  histoNameVector.push_back("phoEta");
  histoNameVector.push_back("phoPhi");


//   histoNameVector.push_back("phoEtAllEcal");
//   histoNameVector.push_back("EfficiencyVsEtHLT");
//   histoNameVector.push_back("EfficiencyVsEtaHLT");
//   histoNameVector.push_back("EfficiencyVsEtLoose");
//   histoNameVector.push_back("EfficiencyVsEtaLoose");


  // to do:  what do we want in certification contents?
  //  dbe_->setCurrentFolder("Egamma/EventInfo/CertificationContents/");


  //looping over histograms to be tested
  if(verbose_) std::cout << "\n>>> looping over histograms to be tested <<<\n\n";
  for(std::vector<string>::iterator it=histoNameVector.begin();it!=histoNameVector.end();++it){

    string HistoName = (*it);
    if(verbose_) std::cout << ">>> " << HistoName;

    MonitorElement * TestHist=0;

    if(HistoName.find("Efficiency")!=std::string::npos){
      TestHist = dbe_->get("Egamma/PhotonAnalyzer/Efficiencies/"+HistoName);
    }
    else{
      TestHist = dbe_->get("Egamma/PhotonAnalyzer/GoodCandidatePhotons/Et above 0 GeV/"+HistoName);
    }
    bool validMe = TestHist!=0;
    if(verbose_)  std::cout << " is valid? " << validMe << "\n";
    if(!validMe) continue;

    if(verbose_)  std::cout << ">>> TestHist Name: " << TestHist->getName() << "\n\n";



    //get QReports associated to each ME
    std::vector<QReport *> myQReports = TestHist->getQReports();
    if(verbose_)  cout << TestHist->getName() <<": myQReports.size() = " << myQReports.size() << "\n\n";
    for(uint i=0;i!=myQReports.size();++i) {

      std::string qtname = myQReports[i]->getQRName() ; // get QT name
      float qtresult = myQReports[i]->getQTresult(); // get QT result value
      int qtstatus = myQReports[i]->getStatus() ; // get QT status value

      if(verbose_) std::cout << "\tTest " << i << ":  qtname: " << qtname   << "\n";
      if(verbose_) std::cout << "\tTest " << i << ":  qtresult: " << qtresult  << std::endl;
      if(verbose_) std::cout << "\tTest " << i << ":  qtstatus: " << qtstatus    << "\n\n";


      //book and fill float for each test done
      dbe_->setCurrentFolder("Egamma/EventInfo/reportSummaryContents/");
      MonitorElement * Float  = dbe_->bookFloat(HistoName+"_"+qtname);
      Float->Fill(qtresult);

    }
  }



  //filling summaries based on geometric spikes/holes

  float etaSpikeTestResult=0;
  float phiSpikeTestResult=0;
  float etaHoleTestResult=0;
  float phiHoleTestResult=0;

  if(dbe_->get("Egamma/EventInfo/reportSummaryContents/phoEta_SpikeTest")!=0)
    etaSpikeTestResult = dbe_->get("Egamma/EventInfo/reportSummaryContents/phoEta_SpikeTest")->getFloatValue();
  else  etaSpikeTestResult = -1;

  if(dbe_->get("Egamma/EventInfo/reportSummaryContents/phoPhi_SpikeTest")!=0)
    phiSpikeTestResult = dbe_->get("Egamma/EventInfo/reportSummaryContents/phoPhi_SpikeTest")->getFloatValue();
  else  phiSpikeTestResult = -1;

  if(dbe_->get("Egamma/EventInfo/reportSummaryContents/phoEta_HoleTest")!=0)
    etaHoleTestResult = dbe_->get("Egamma/EventInfo/reportSummaryContents/phoEta_HoleTest")->getFloatValue();
  else  etaHoleTestResult = -1;

  if(dbe_->get("Egamma/EventInfo/reportSummaryContents/phoPhi_HoleTest")!=0)
    phiHoleTestResult = dbe_->get("Egamma/EventInfo/reportSummaryContents/phoPhi_HoleTest")->getFloatValue();
  else  phiHoleTestResult = -1;



  if(verbose_) std::cout << ">>>  Results of tests to be put into Summary Maps  <<<\n\n";
  if(verbose_) std::cout << "\tetaSpikeTestResult= " << etaSpikeTestResult << "\n";
  if(verbose_) std::cout << "\tphiSpikeTestResult= " << phiSpikeTestResult << "\n";
  if(verbose_) std::cout << "\tetaHoleTestResult= " << etaHoleTestResult << "\n";
  if(verbose_) std::cout << "\tphiHoleTestResult= " << phiHoleTestResult << "\n\n";


  //fill reportSummary & CertificationSummary with average of hole/spike tests
  float reportSummaryFloat = (etaSpikeTestResult+etaHoleTestResult+phiSpikeTestResult+phiHoleTestResult)/4.;
  if(reportSummaryFloat<0) reportSummaryFloat = -1;
  reportSummary->Fill(reportSummaryFloat);
  CertificationSummary->Fill(reportSummaryFloat);

  reportSummaryMap->Fill(0,0,etaSpikeTestResult);
  reportSummaryMap->Fill(0,1,etaHoleTestResult);
  reportSummaryMap->Fill(1,0,phiSpikeTestResult);
  reportSummaryMap->Fill(1,1,phiHoleTestResult);

  CertificationSummaryMap->Fill(0,0,etaSpikeTestResult);
  CertificationSummaryMap->Fill(0,1,etaHoleTestResult);
  CertificationSummaryMap->Fill(1,0,phiSpikeTestResult);
  CertificationSummaryMap->Fill(1,1,phiHoleTestResult);




}

Member Data Documentation

DQMStore* PhotonDataCertification::dbe_

private

Definition at line 66 of file PhotonDataCertification.h.

edm::ParameterSet PhotonDataCertification::parameters_

private

Definition at line 67 of file PhotonDataCertification.h.

Referenced by Mixins._TypedParameterizable::clone(), Types.PSet::clone(), and Mixins._TypedParameterizable::copy().

bool PhotonDataCertification::verbose_

private

Definition at line 69 of file PhotonDataCertification.h.