|
|
|
#include <HLTMonPhotonClient.h>
Public Member Functions | |
| HLTMonPhotonClient (const edm::ParameterSet &) | |
| ~HLTMonPhotonClient () | |
Private Member Functions | |
| virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
| virtual void | beginJob () |
| virtual void | endJob () |
Private Attributes | |
| MonitorElement * | cumFilterEff |
| DQMStore * | dbe |
| std::string | dirname_ |
| MonitorElement * | eventCounter |
| ofstream | logFile_ |
| bool | monitorDaemon_ |
| std::string | outputFile_ |
| MonitorElement * | pixelhistosEt [4] |
| MonitorElement * | pixelhistosEta [4] |
| MonitorElement * | pixelhistosEtaOut [4] |
| MonitorElement * | pixelhistosEtOut [4] |
| MonitorElement * | pixelhistosPhi [4] |
| MonitorElement * | pixelhistosPhiOut [4] |
| MonitorElement * | relFilterEff |
| std::string | sourcedirname_ |
| edm::InputTag | sourcetag_ |
| std::vector< edm::InputTag > | theHLTCollectionLabels |
Definition at line 30 of file HLTMonPhotonClient.h.
| HLTMonPhotonClient::HLTMonPhotonClient | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 14 of file HLTMonPhotonClient.cc.
References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), LogDebug, NULL, and cmsCodeRules::cppFunctionSkipper::operator.
{
LogDebug("HLTMonPhotonClient") << "constructor...." ;
logFile_.open("HLTMonPhotonClient.log");
dbe = NULL;
if (iConfig.getUntrackedParameter < bool > ("DQMStore", false)) {
dbe = Service < DQMStore > ().operator->();
dbe->setVerbose(0);
}
outputFile_ =
iConfig.getUntrackedParameter <std::string>("outputFile", "");
if (outputFile_.size() != 0) {
LogInfo("HLTMonPhotonClient") << "Photon Trigger Monitoring histograms will be saved to "
<< outputFile_ ;
}
else {
outputFile_ = "PhotonDQM.root";
}
bool disable =
iConfig.getUntrackedParameter < bool > ("disableROOToutput", false);
if (disable) {
outputFile_ = "";
}
sourcetag_=iConfig.getParameter<edm::InputTag>("SourceTag");
theHLTCollectionLabels = iConfig.getParameter<std::vector<edm::InputTag> >("theHLTCollectionLabels");
dirname_="HLT/HLTMonPhoton/"+iConfig.getParameter<std::string>("@module_label");
sourcedirname_="HLT/HLTMonPhoton/"+sourcetag_.label();
if (dbe != NULL) {
dbe->setCurrentFolder(dirname_);
}
}
| HLTMonPhotonClient::~HLTMonPhotonClient | ( | ) |
Definition at line 57 of file HLTMonPhotonClient.cc.
{
// do anything here that needs to be done at desctruction time
// (e.g. close files, deallocate resources etc.)
}
| void HLTMonPhotonClient::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDAnalyzer.
Definition at line 72 of file HLTMonPhotonClient.cc.
{
TH1F* refhisto = eventCounter->getTH1F();
if(refhisto->GetBinContent(1) != 0){
// relFilterEff->setBinContent(6,refhisto->GetBinContent(5)/refhisto->GetBinContent(2));
// cumFilterEff->setBinContent(5,refhisto->GetBinContent(5)/refhisto->GetBinContent(1));
cumFilterEff->setBinContent(4,refhisto->GetBinContent(4)/refhisto->GetBinContent(1));
cumFilterEff->setBinContent(3,refhisto->GetBinContent(3)/refhisto->GetBinContent(1));
cumFilterEff->setBinContent(2,refhisto->GetBinContent(2)/refhisto->GetBinContent(1));
}else{
relFilterEff->setBinContent(6,0);
cumFilterEff->setBinContent(5,0);
cumFilterEff->setBinContent(4,0);
cumFilterEff->setBinContent(3,0);
cumFilterEff->setBinContent(2,0);
}
for(int i = 0; i<4; i++){
if(refhisto->GetBinContent(4-i) != 0){
relFilterEff->setBinContent(5-i,refhisto->GetBinContent(5-i)/refhisto->GetBinContent(4-i));
}else{
relFilterEff->setBinContent(5-i,0);
}
}
if(refhisto->GetBinContent(1) != 0){
relFilterEff->setBinContent(1,refhisto->GetBinContent(1)/refhisto->GetBinContent(1));
cumFilterEff->setBinContent(2,refhisto->GetBinContent(2)/refhisto->GetBinContent(1));
cumFilterEff->setBinContent(1,refhisto->GetBinContent(1)/refhisto->GetBinContent(1));
}else{
relFilterEff->setBinContent(1,0);
cumFilterEff->setBinContent(2,0);
cumFilterEff->setBinContent(1,0);
}
TH1F* num;
TH1F* denom;
for(int i=0; i<3; i++){
denom = pixelhistosEt[0]->getTH1F();
num = pixelhistosEt[i+1]->getTH1F();
for(int j=1; j <= pixelhistosEtOut[i]->getNbinsX();j++ ){
if(denom->GetBinContent(j)!=0){
pixelhistosEtOut[i]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
}
else{
pixelhistosEtOut[i]->setBinContent(j,0.);
}
}
denom = pixelhistosEta[0]->getTH1F();
num = pixelhistosEta[i+1]->getTH1F();
for(int j=1; j <= pixelhistosEtaOut[i]->getNbinsX();j++ ){
if(denom->GetBinContent(j)!=0)
pixelhistosEtaOut[i]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
else
pixelhistosEtaOut[i]->setBinContent(j,0.);
}
denom = pixelhistosPhi[0]->getTH1F();
num = pixelhistosPhi[i+1]->getTH1F();
for(int j=1; j <= pixelhistosPhiOut[i]->getNbinsX();j++ ){
if(denom->GetBinContent(j)!=0)
pixelhistosPhiOut[i]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
else
pixelhistosPhiOut[i]->setBinContent(j,0.);
}
}
denom = pixelhistosEt[0]->getTH1F();
num = pixelhistosEt[3]->getTH1F();
for(int j=1; j <= pixelhistosEtOut[3]->getNbinsX();j++ ){
if(denom->GetBinContent(j)!=0){
pixelhistosEtOut[3]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
}
else{
pixelhistosEtOut[3]->setBinContent(j,0.);
}
}
denom = pixelhistosEta[0]->getTH1F();
num = pixelhistosEta[3]->getTH1F();
for(int j=1; j <= pixelhistosEtaOut[3]->getNbinsX();j++ ){
if(denom->GetBinContent(j)!=0){
pixelhistosEtaOut[3]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
}
else{
pixelhistosEtaOut[3]->setBinContent(j,0.);
}
}
denom = pixelhistosPhi[0]->getTH1F();
num = pixelhistosPhi[3]->getTH1F();
for(int j=1; j <= pixelhistosPhiOut[3]->getNbinsX();j++ ){
if(denom->GetBinContent(j)!=0){
pixelhistosPhiOut[3]->setBinContent(j,num->GetBinContent(j)/denom->GetBinContent(j));
}
else{
pixelhistosPhiOut[3]->setBinContent(j,0.);
}
}
}
| void HLTMonPhotonClient::beginJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 175 of file HLTMonPhotonClient.cc.
References DQMStore::book1D(), DQMStore::get(), MonitorElement::getTH1F(), i, cmsCodeRules::cppFunctionSkipper::operator, DQMStore::rmdir(), MonitorElement::setBinLabel(), and DQMStore::setCurrentFolder().
{
DQMStore *dbe = 0;
dbe = Service < DQMStore > ().operator->();
if (dbe) {
dbe->setCurrentFolder(dirname_);
dbe->rmdir(dirname_);
}
if (dbe) {
dbe->setCurrentFolder(dirname_);
std::string tmpname = "";
for(int i = 0; i<4; i++){
LogInfo("HLTMonPhotonClient") << "loop iteration: "<<i ;
tmpname = sourcedirname_ + "/" + theHLTCollectionLabels[i].label() + " Eta Dist";
pixelhistosEta[i]=dbe->get(tmpname);
tmpname = sourcedirname_ + "/" + theHLTCollectionLabels[i].label() + " Et Dist";
pixelhistosEt[i]=dbe->get(tmpname);
tmpname = sourcedirname_ + "/" + theHLTCollectionLabels[i].label() + " Phi Dist";
pixelhistosPhi[i]=dbe->get(tmpname);
}
TH1F* refhist;
int nBins;
double xMin;
double xMax;
refhist = pixelhistosEt[0]->getTH1F();
nBins = refhist->GetNbinsX();
xMin = refhist->GetXaxis()->GetXmin();
xMax = refhist->GetXaxis()->GetXmax();
pixelhistosEtOut[0] =dbe->book1D("EcalIsol Eff vs Et","EcalIsol Eff vs Et",nBins,xMin,xMax);
pixelhistosEtOut[1] =dbe->book1D("HcalIsol Eff vs Et","HcalIsol Eff vs Et",nBins,xMin,xMax);
pixelhistosEtOut[2] =dbe->book1D("Track Isol Eff vs Et","Track Isol Eff vs Et",nBins,xMin,xMax);
pixelhistosEtOut[3] =dbe->book1D("Total Eff vs Et","Total Eff vs Et",nBins,xMin,xMax);
refhist = pixelhistosEta[0]->getTH1F();
nBins = refhist->GetNbinsX();
xMin = refhist->GetXaxis()->GetXmin();
xMax = refhist->GetXaxis()->GetXmax();
pixelhistosEtaOut[0] =dbe->book1D("EcalIsol Eff vs Eta","EcalIsol Eff vs Eta",nBins,xMin,xMax);
pixelhistosEtaOut[1] =dbe->book1D("HcalIsol Eff vs Eta","HcalIsol Eff vs Eta",nBins,xMin,xMax);
pixelhistosEtaOut[2] =dbe->book1D("Track Isol Eff vs Eta","Track Isol Eff vs Eta",nBins,xMin,xMax);
pixelhistosEtaOut[3] =dbe->book1D("Total Eff vs Eta","Total Eff vs Eta",nBins,xMin,xMax);
refhist = pixelhistosPhi[0]->getTH1F();
nBins = refhist->GetNbinsX();
xMin = refhist->GetXaxis()->GetXmin();
xMax = refhist->GetXaxis()->GetXmax();
pixelhistosPhiOut[0] =dbe->book1D("EcalIsol Eff vs Phi","EcalIsol Eff vs Phi",nBins,xMin,xMax);
pixelhistosPhiOut[1] =dbe->book1D("HcalIsol Eff vs Phi","HcalIsol Eff vs Phi",nBins,xMin,xMax);
pixelhistosPhiOut[2] =dbe->book1D("Track Isol Eff vs Phi","Track Isol Eff vs Phi",nBins,xMin,xMax);
pixelhistosPhiOut[3] =dbe->book1D("Total Eff vs Phi","Total Eff vs Phi",nBins,xMin,xMax);
tmpname = sourcedirname_ + "/Evts Passing Filters";
eventCounter = dbe->get(tmpname);
relFilterEff = dbe->book1D("Relative Filter Effs","Relative Filter Effs",5,0,5);
relFilterEff->setBinLabel(1,"HLTIsoEtFilterEff");
relFilterEff->setBinLabel(2,"HLTIsoEcalIsolEff");
relFilterEff->setBinLabel(3,"HLTIsoHcalIsolEff");
relFilterEff->setBinLabel(4,"HLTIsoTrackIsolEff");
// relFilterEff->setBinLabel(5,"FIXME");
// relFilterEff->setBinLabel(6,"FIXME");
cumFilterEff = dbe->book1D("Cumulative Filter Effs","Cumulative Filter Effs",5,0,5);
cumFilterEff->setBinLabel(1,"HLTIsoEtFilterEff");
cumFilterEff->setBinLabel(2,"HLTIsoEcalIsolEff");
cumFilterEff->setBinLabel(3,"HLTIsoHcalIsolEff");
cumFilterEff->setBinLabel(4,"HLTIsoTrackIsolEff");
// cumFilterEff->setBinLabel(5,"TrackIsolEff");
} // end "if(dbe)"
}
| void HLTMonPhotonClient::endJob | ( | void | ) | [private, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 256 of file HLTMonPhotonClient.cc.
{
// std::cout << "HLTMonPhotonClient: end job...." << std::endl;
if (outputFile_.size() != 0 && dbe)
dbe->save(outputFile_);
return;
}
MonitorElement* HLTMonPhotonClient::cumFilterEff [private] |
Definition at line 55 of file HLTMonPhotonClient.h.
DQMStore* HLTMonPhotonClient::dbe [private] |
Definition at line 43 of file HLTMonPhotonClient.h.
std::string HLTMonPhotonClient::dirname_ [private] |
Definition at line 57 of file HLTMonPhotonClient.h.
MonitorElement* HLTMonPhotonClient::eventCounter [private] |
Definition at line 53 of file HLTMonPhotonClient.h.
ofstream HLTMonPhotonClient::logFile_ [private] |
Definition at line 60 of file HLTMonPhotonClient.h.
bool HLTMonPhotonClient::monitorDaemon_ [private] |
Definition at line 59 of file HLTMonPhotonClient.h.
std::string HLTMonPhotonClient::outputFile_ [private] |
Definition at line 61 of file HLTMonPhotonClient.h.
MonitorElement* HLTMonPhotonClient::pixelhistosEt[4] [private] |
Definition at line 47 of file HLTMonPhotonClient.h.
MonitorElement* HLTMonPhotonClient::pixelhistosEta[4] [private] |
Definition at line 48 of file HLTMonPhotonClient.h.
MonitorElement* HLTMonPhotonClient::pixelhistosEtaOut[4] [private] |
Definition at line 51 of file HLTMonPhotonClient.h.
MonitorElement* HLTMonPhotonClient::pixelhistosEtOut[4] [private] |
Definition at line 50 of file HLTMonPhotonClient.h.
MonitorElement* HLTMonPhotonClient::pixelhistosPhi[4] [private] |
Definition at line 49 of file HLTMonPhotonClient.h.
MonitorElement* HLTMonPhotonClient::pixelhistosPhiOut[4] [private] |
Definition at line 52 of file HLTMonPhotonClient.h.
MonitorElement* HLTMonPhotonClient::relFilterEff [private] |
Definition at line 54 of file HLTMonPhotonClient.h.
std::string HLTMonPhotonClient::sourcedirname_ [private] |
Definition at line 58 of file HLTMonPhotonClient.h.
edm::InputTag HLTMonPhotonClient::sourcetag_ [private] |
Definition at line 45 of file HLTMonPhotonClient.h.
std::vector<edm::InputTag> HLTMonPhotonClient::theHLTCollectionLabels [private] |
Definition at line 62 of file HLTMonPhotonClient.h.
1.7.3