PhotonDataCertification.cc

Go to the documentation of this file.

#include <iostream>

#include "DQMOffline/EGamma/plugins/PhotonDataCertification.h"
#include "FWCore/Framework/interface/Run.h"
#include "RooGlobalFunc.h"
#include "RooRealVar.h"
#include "RooDataSet.h"
#include "RooBreitWigner.h"
#include "RooDataHist.h"
#include "RooFitResult.h"


//
// Original Author:  Louis James Antonelli
//         Created:  Thu Jan 22 13:42:28CET 2009
//

using namespace std;

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

{

  
  parameters_ = pset;
  verbose_ = parameters_.getParameter<bool>("verbose");

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

}


PhotonDataCertification::~PhotonDataCertification()
{
}


//
// member functions
//

// ------------ method called to for each event  ------------
void
PhotonDataCertification::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
{
   //if(verbose_) std::cout << ">>> Analyze (PhotonDataCertification) <<<" << std::endl;

}

// ------------ method called once each job just before starting event loop  ------------
void
PhotonDataCertification::beginJob()
{
  dbe_ = 0;
  dbe_ = edm::Service<DQMStore>().operator->();
  dbe_->setVerbose(0);
}

// ------------ method called once each job just after ending the event loop  ------------
void
PhotonDataCertification::endJob()
{
}

// ------------ method called just before starting a new run  ------------
void
PhotonDataCertification::beginRun(const edm::Run& run, const edm::EventSetup& c)
{

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

}

// ------------ method called right after a run ends ------------
void
PhotonDataCertification::endRun(const edm::Run& run, const edm::EventSetup& c)
{
  using namespace RooFit;
  if(verbose_) std::cout << ">>> EndRun (PhotonDataCertification) <<<" << std::endl;
  if(!dbe_) return;
  
  //booking histograms according to naming conventions
  dbe_->setCurrentFolder("Egamma/EventInfo/");
  reportSummary_ = dbe_->bookFloat("reportSummary");
  reportSummaryMap_ = dbe_->book2D("reportSummaryMap","reportSummaryMap",3,0,3,1,0,1);
  
  TH2F * reportSummaryMapTH2 = reportSummaryMap_->getTH2F();
  reportSummaryMapTH2->GetXaxis()->SetBinLabel(1,"EB");
  reportSummaryMapTH2->GetXaxis()->SetBinLabel(2,"EE");
  reportSummaryMapTH2->GetXaxis()->SetBinLabel(3,"Total");
  reportSummaryMapTH2->GetYaxis()->SetBinLabel(1,"InvMassTest");
  
  float EBResult = invMassZtest("Egamma/PhotonAnalyzer/InvMass/h_02_invMassIsoPhotonsEBarrel","invMassIsolatedPhotonsEB");
  float EEResult = invMassZtest("Egamma/PhotonAnalyzer/InvMass/h_03_invMassIsoPhotonsEEndcap","invMassIsolatedPhotonsEE");
  float AllResult = invMassZtest("Egamma/PhotonAnalyzer/InvMass/h_01_invMassAllIsolatedPhotons","invMassAllIsolatedPhotons");

  if(verbose_){
    std::cout << "EBResult: " << EBResult << std::endl;
    std::cout << "EEResult: " << EEResult << std::endl;
    std::cout << "AllResult: " << AllResult << std::endl;
  }
  
  reportSummaryMapTH2->SetBinContent(1, 1, EBResult);
  reportSummaryMapTH2->SetBinContent(2, 1, EEResult);
  reportSummaryMapTH2->SetBinContent(3, 1, AllResult);
  reportSummary_->Fill(AllResult);
 
}

float PhotonDataCertification::invMassZtest(string path, TString name){
  float ZMass = 91.2;
  float ZWidth = 2.5;
  MonitorElement *TestElem=0;
  TestElem = dbe_->get(path);
  if(TestElem==0) return 0;
  TH1F *TestHist = TestElem->getTH1F();
  if(TestHist==0) return 0;
  RooRealVar mass("mass","Mass_{2#gamma}", 0, 200,"GeV");
  RooRealVar mRes("M_{Z}", "Z Mass", ZMass, 70, 110);
  RooRealVar gamma("#Gamma", "#Gamma", ZWidth, 0, 10.0);
  RooBreitWigner BreitWigner("BreitWigner","Breit-Wigner",mass,mRes,gamma);
  RooDataHist test(name, name, mass, TestHist);

  BreitWigner.fitTo(test, RooFit::Range(80,100));

  if(abs(mRes.getValV() - ZMass) < ZWidth){return 1.0;}
  else if(abs(mRes.getValV() - ZMass) < gamma.getValV()){return 0.9;}
  else{return 0.0;}
 
}