|
|
|
#include <FourVectorHLTClient.h>
Public Member Functions | |
| FourVectorHLTClient (const edm::ParameterSet &ps) | |
| Constructor. | |
| virtual | ~FourVectorHLTClient () |
| Destructor. | |
Protected Member Functions | |
| void | analyze (const edm::Event &e, const edm::EventSetup &c) |
| Fake Analyze. | |
| void | beginJob () |
| BeginJob. | |
| void | beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context) |
| void | beginRun (const edm::Run &r, const edm::EventSetup &c) |
| BeginRun. | |
| void | endJob () |
| Endjob. | |
| void | endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c) |
| DQM Client Diagnostic. | |
| void | endRun (const edm::Run &r, const edm::EventSetup &c) |
| EndRun. | |
| void | normalizeHLTMatrix () |
Private Member Functions | |
| void | calculateRatio (TH1F *effHist, TH1F *denHist) |
| TH1F * | get1DHisto (std::string meName, DQMStore *dbi) |
| TProfile * | get1DProfile (std::string meName, DQMStore *dbi) |
| TH2F * | get2DHisto (std::string meName, DQMStore *dbi) |
| TProfile2D * | get2DProfile (std::string meName, DQMStore *dbi) |
| void | initialize () |
Private Attributes | |
| TString | clientDir_ |
| int | counterEvt_ |
| counter | |
| int | counterLS_ |
| TString | customEffDir_ |
| std::vector< std::pair < std::string, std::string > > | custompathnamepairs_ |
| DQMStore * | dbe_ |
| HLTConfigProvider | hltConfig_ |
| std::vector< TString > | hltMEName |
| std::vector< MonitorElement * > | hltMEs |
| std::vector< TString > | hltPathName |
| MonitorElement * | klmgrvTest_ |
| int | nChannels |
| prescale on number of events | |
| edm::ParameterSet | parameters_ |
| std::string | pathsSummaryFilterCountsFolder_ |
| std::string | pathsSummaryFilterEfficiencyFolder_ |
| std::string | pathsSummaryFolder_ |
| std::string | pathsSummaryHLTCorrelationsFolder_ |
| int | prescaleEvt_ |
| units of lumi sections | |
| int | prescaleLS_ |
| counter | |
| std::string | processname_ |
| Float_t | reportSummary |
| MonitorElement * | reportSummary_ |
| MonitorElement * | reportSummaryContent_ [20] |
| MonitorElement * | reportSummaryMap_ |
| TString | sourceDir_ |
| Float_t | summaryContent [20] |
| Float_t | summarySum |
| MonitorElement * | testHLTEff_ |
| std::vector< MonitorElement * > | v_ME_HLTPass_Normalized_Any_ |
| std::vector< MonitorElement * > | v_ME_HLTPassPass_ |
| std::vector< MonitorElement * > | v_ME_HLTPassPass_Normalized_ |
Definition at line 31 of file FourVectorHLTClient.h.
| FourVectorHLTClient::FourVectorHLTClient | ( | const edm::ParameterSet & | ps | ) |
Constructor.
Definition at line 43 of file FourVectorHLTClient.cc.
{
parameters_=ps;
initialize();
}
| FourVectorHLTClient::~FourVectorHLTClient | ( | ) | [virtual] |
Destructor.
Definition at line 51 of file FourVectorHLTClient.cc.
References LogDebug.
{
LogDebug("FourVectorHLTClient")<< "FourVectorHLTClient: ending...." ;
}
| void FourVectorHLTClient::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
Fake Analyze.
Implements edm::EDAnalyzer.
Definition at line 159 of file FourVectorHLTClient.cc.
References LogDebug.
{
counterEvt_++;
if (prescaleEvt_<1) return;
if (prescaleEvt_>0 && counterEvt_%prescaleEvt_ != 0) return;
LogDebug("FourVectorHLTClient")<<"analyze..." << endl;
}
| void FourVectorHLTClient::beginJob | ( | void | ) | [protected, virtual] |
BeginJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 116 of file FourVectorHLTClient.cc.
References dbe_, LogDebug, and cmsCodeRules::cppFunctionSkipper::operator.
{
LogDebug("FourVectorHLTClient")<<"[FourVectorHLTClient]: beginJob" << endl;
// get backendinterface
dbe_ = Service<DQMStore>().operator->();
}
| void FourVectorHLTClient::beginLuminosityBlock | ( | const edm::LuminosityBlock & | lumiSeg, |
| const edm::EventSetup & | context | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 149 of file FourVectorHLTClient.cc.
{
// optionally reset histograms here
}
| void FourVectorHLTClient::beginRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
BeginRun.
Reimplemented from edm::EDAnalyzer.
Definition at line 125 of file FourVectorHLTClient.cc.
References LogDebug.
{
LogDebug("FourVectorHLTClient")<<"[FourVectorHLTClient]: beginRun" << endl;
// HLT config does not change within runs!
bool changed=false;
if (!hltConfig_.init(r, context, processname_, changed)) {
processname_ = "FU";
if (!hltConfig_.init(r, context, processname_, changed)){
LogDebug("FourVectorHLTOffline") << "HLTConfigProvider failed to initialize.";
}
// check if trigger name in (new) config
// cout << "Available TriggerNames are: " << endl;
// hltConfig_.dump("Triggers");
}
}
| void FourVectorHLTClient::calculateRatio | ( | TH1F * | effHist, |
| TH1F * | denHist | ||
| ) | [private] |
Definition at line 643 of file FourVectorHLTClient.cc.
References i, gen::k, MultiGaussianStateTransform::N, and mathSSE::sqrt().
{
//cout << " NUMERATOR histogram name = " << effHist->GetName() << endl;
int nBins= effHist->GetNbinsX();
for (int i=0;i<=nBins;i++) {
float k = effHist->GetBinContent(i); // number of pass
float N = denHist->GetBinContent(i); // number of total
float ratio;
float err;
//ratio = N ? k / N : 0;
if(N>0) ratio = k / N;
else ratio = 0;
if(N > 0) {
if(ratio <= 1) {
err = sqrt(ratio*(1-ratio)/N);
}
//else if(ratio == 1) {
//err = 1./sqrt(2*N);
//}
else {
err = 0;
}
effHist->SetBinContent(i,ratio);
effHist->SetBinError(i,err);
}
}
}
| void FourVectorHLTClient::endJob | ( | void | ) | [protected, virtual] |
Endjob.
Reimplemented from edm::EDAnalyzer.
Definition at line 581 of file FourVectorHLTClient.cc.
{
}
| void FourVectorHLTClient::endLuminosityBlock | ( | const edm::LuminosityBlock & | lumiSeg, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
DQM Client Diagnostic.
Reimplemented from edm::EDAnalyzer.
Definition at line 153 of file FourVectorHLTClient.cc.
{
}
| void FourVectorHLTClient::endRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
EndRun.
add dataset name and thier triggers to the list
Reimplemented from edm::EDAnalyzer.
Definition at line 170 of file FourVectorHLTClient.cc.
References newFWLiteAna::bin, DQMStore::book1D(), DQMStore::book2D(), dbe_, DQMStore::get(), DQMStore::getContents(), edm::getName(), DQMStore::getSubdirs(), MonitorElement::getTH1F(), MonitorElement::getTH2F(), i, gen::k, prof2calltree::l, LogDebug, m, AlCaRecoCosmics_cfg::name, DQMStore::setCurrentFolder(), and findQualityFiles::size.
{
LogDebug("FourVectorHLTClient")<<"FourVectorHLTClient:: endLuminosityBlock " << endl;
// Work with PathsSummary ME's
dbe_->setCurrentFolder(pathsSummaryFolder_.c_str());
std::vector<MonitorElement*> summaryME;
summaryME = dbe_->getContents(pathsSummaryFilterCountsFolder_.c_str());
vector<float> hltPassCount;
for(unsigned int i=0; i< summaryME.size(); i++) {
TString nameME = (summaryME[i])->getName();
LogDebug("FourVectorHLTClient") << "ME = " << nameME << endl;
TString fullPathToME = pathsSummaryFilterCountsFolder_ + nameME;
LogDebug("FourVectorHLTClient") << "fullPathME = " << fullPathToME << endl;
if(nameME.Contains("Filters_")) {
TH1F* MEHist = summaryME[i]->getTH1F() ;
LogDebug("FourVectorHLTClient") << "this is TH1 first bin label = " << MEHist->GetXaxis()->GetBinLabel(1) << endl;
LogDebug("FourVectorHLTClient") << "is this is TH1 = " << fullPathToME << endl;
TString nameMEeff = nameME+"_Eff";
TString titleMEeff = nameME+" Efficiency wrt input to L1";
TH1F* effHist = (TH1F*) MEHist->Clone(nameMEeff.Data());
effHist->SetTitle(titleMEeff.Data());
float firstBinContent = MEHist->GetBinContent(1);
if(firstBinContent !=0 ) {
effHist->Scale(1./firstBinContent);
}
else {
unsigned int nBins = effHist->GetNbinsX();
for(unsigned int bin = 0; bin < nBins+1; bin++) effHist->SetBinContent(bin, 0);
}
dbe_->setCurrentFolder(pathsSummaryFilterEfficiencyFolder_.c_str());
dbe_->book1D(nameMEeff.Data(), effHist);
} // end if "Filters_" ME
} // end for MEs
// Normalize PassPass and PassFail Matrices
vector<string> name;
name.push_back("All");
name.push_back("Muon");
name.push_back("Egamma");
name.push_back("Tau");
name.push_back("JetMET");
name.push_back("Rest");
name.push_back("Special");
vector<string> datasetNames = hltConfig_.datasetNames() ;
for (unsigned int i=0;i<datasetNames.size();i++) {
name.push_back(datasetNames[i]);
}
string fullPathToME;
for (unsigned int i=0;i<name.size();i++) {
fullPathToME = pathsSummaryFolder_ +"HLT_"+name[i]+"_PassPass";
v_ME_HLTPassPass_.push_back( dbe_->get(fullPathToME));
fullPathToME = pathsSummaryHLTCorrelationsFolder_+"HLT_"+name[i]+"_PassPass_Normalized";
v_ME_HLTPassPass_Normalized_.push_back( dbe_->get(fullPathToME));
fullPathToME = pathsSummaryHLTCorrelationsFolder_+"HLT_"+name[i]+"_Pass_Normalized_Any";
v_ME_HLTPass_Normalized_Any_.push_back( dbe_->get(fullPathToME));
}
normalizeHLTMatrix();
// HLT paths
TObjArray *hltPathNameColl = new TObjArray(100);
hltPathNameColl->SetOwner();
std::vector<std::string> fullPathHLTFolders;
// get subdir of the sourceDir_ to get HLT path names
dbe_->setCurrentFolder(sourceDir_.Data());
fullPathHLTFolders = dbe_->getSubdirs();
LogDebug("FourVectorHLTClient")<<"endRun: sourceDir_(" << sourceDir_.Data() << ") has " << fullPathHLTFolders.size() << " folders. " << endl;
for(unsigned int i=0;i<fullPathHLTFolders.size();i++) {
LogDebug("FourVectorHLTClient")<<"endRun: sourceDir_(" << sourceDir_.Data() << ") folder["<< i << "] = " << fullPathHLTFolders[i] << endl;
TString hltPath = fullPathHLTFolders[i].substr(fullPathHLTFolders[i].rfind('/')+1, fullPathHLTFolders[i].size());
// Efficiencies
TString currEffFolder = clientDir_ + "/" + hltPath + "/" + customEffDir_ + "/";
LogDebug("FourVectorHLTClient")<< "Custom Efficiencies dir path = " << currEffFolder << endl;
dbe_->setCurrentFolder(currEffFolder.Data());
hltMEs = dbe_->getContents(fullPathHLTFolders[i]);
LogDebug("FourVectorHLTClient")<< "Number of MEs for this HLT path = " << hltMEs.size() << endl;
// custom efficiencies
for (std::vector<std::pair<std::string, std::string> >::iterator custompathnamepair = custompathnamepairs_.begin(); custompathnamepair != custompathnamepairs_.end(); ++custompathnamepair)
{
//TString numPathName=TString(custompathnamepair->first);
//if(!hltPath.Contains(numPathName,TString::kExact)) continue;
//TString numPathName=hltPath;
//if(!hltPath.Contains(numPathName,TString::kExact)) continue;
//LogDebug("FourVectorHLTClient")<< "hltPath = " << hltPath << " matchString = " << custompathnamepair->first << endl;
//if(!hltPath.Contains(TString(custompathnamepair->first),TString::kExact)) continue;
if(!hltPath.Contains(TString(custompathnamepair->first))) continue;
TString numPathName=hltPath;
TString denPathName=TString(custompathnamepair->second);
vector<TString> vStage;
vStage.push_back(TString("L1"));
vStage.push_back(TString("On"));
vStage.push_back(TString("Off"));
vStage.push_back(TString("Mc"));
vStage.push_back(TString("L1UM"));
vStage.push_back(TString("OnUM"));
vStage.push_back(TString("OffUM"));
vStage.push_back(TString("McUM"));
vector<TString> vObj;
vObj.push_back(TString("l1Et"));
vObj.push_back(TString("onEt"));
vObj.push_back(TString("offEt"));
vObj.push_back(TString("mcEt"));
// TH1Fs
for (unsigned int k=0; k<vObj.size();k++) {
for (unsigned int l=0; l<vStage.size();l++) {
for (unsigned int m=0; m<vStage.size();m++) {
TString oldHistPathNum;
TString oldHistPathNumBck;
TString oldHistPathDen;
// not the differeence b/w Num and NumBck, as in OnOff vs OffOn
oldHistPathNum = sourceDir_+"/"+hltPath+"/"+numPathName+"_wrt_"+denPathName+"_"+vObj[k]+vStage[l]+vStage[m];
oldHistPathNumBck = sourceDir_+"/"+hltPath+"/"+numPathName+"_wrt_"+denPathName+"_"+vObj[k]+vStage[m]+vStage[l];
// In the deominator hist name, we don't have any "UM" substrings, so remove them
TString tempDenString = vStage[l];
tempDenString.ReplaceAll("UM","") ;
oldHistPathDen = sourceDir_+"/"+hltPath+"/"+numPathName+"_wrt_"+denPathName+"_"+vObj[k]+tempDenString;
MonitorElement *numME = dbe_->get(oldHistPathNum.Data());
MonitorElement *numMEBck = dbe_->get(oldHistPathNumBck.Data());
MonitorElement *denME = dbe_->get(oldHistPathDen.Data());
LogDebug("FourVectorHLTClient")<< " oldHistPathNum = " << oldHistPathNum << endl;
LogDebug("FourVectorHLTClient")<< " oldHistPathNumBck = " << oldHistPathNumBck << endl;
LogDebug("FourVectorHLTClient")<< " oldHistPathDen = " << oldHistPathDen << endl;
//check if HLTOffline histogram exist
if ( numME && denME ) {
LogDebug("FourVectorHLTClient")<< "DID find NUM and DEN histograms to derive eff " << vStage[m]<<"To"<<vStage[l]<<" using:" <<endl
<< " NUM = " << oldHistPathNum << endl
<< " DEN = " << oldHistPathDen << endl;
}
else {
LogDebug("FourVectorHLTClient")<< "Cannot find NUM and DEN histograms to derive eff " << vStage[m]<<"To"<<vStage[l]<<" using:" <<endl
<< " NUM = " << oldHistPathNum << endl
<< " DEN = " << oldHistPathDen << endl;
if ( numMEBck && denME) {
LogDebug("FourVectorHLTClient")<< "DID find NUM and DEN histograms to derive eff " << vStage[m]<<"To"<<vStage[l]<<" using:" <<endl
<< " NUM = " << oldHistPathNumBck << endl
<< " DEN = " << oldHistPathDen << endl;
numME = numMEBck;
}
else {
LogDebug("FourVectorHLTClient")<< "Cannot find NUM and DEN histograms to derive eff " << vStage[m]<<"To"<<vStage[l]<<" using:" <<endl
<< " NUM = " << oldHistPathNumBck << endl
<< " DEN = " << oldHistPathDen << endl;
continue;
}
}
// EtaPhi histos are TH2s, the rest are TH1s
TH1F * numHist = numME->getTH1F();
TH1F * denHist = denME->getTH1F();
// build names and title for efficiency histogram
TString newHistName = hltPath +"_wrt_" + denPathName +"_"+vObj[k]+"_Eff_"+vStage[m]+"To"+vStage[l];
// if there is "UM", remove it first,then append it to the end of the name
if(newHistName.Contains("UM")) {
newHistName.ReplaceAll("UM","");
newHistName.Append("_UM");
}
TString newHistTitle = numPathName+" given " + denPathName +" "+vObj[k]+" Eff "+vStage[m]+"To"+vStage[l];
// if there is "UM", remove it first,then append it to the end of the name
if(newHistTitle.Contains("UM")) {
newHistTitle.ReplaceAll("UM","");
newHistTitle.Append("_UM");
}
TString newHistPath = currEffFolder+newHistName;
LogDebug("FourVectorHLTClient")<< "Will make efficiency histogram " << newHistPath << endl;
TH1F* effHist = (TH1F*) numHist->Clone(newHistName.Data());
effHist->SetTitle(newHistTitle.Data());
//effHist->Sumw2();
//denHist->Sumw2();
//effHist->Divide(numHist,denHist);
//effHist->Divide(numHist,denHist,1.,1.,"B");
calculateRatio(effHist, denHist);
/*
effHist->Divide(numHist,denHist,1.,1.);
for (int i=0;i<=effHist->GetNbinsX();i++) {
float err = 0;
if(numHist->GetBinContent(i)>0) err = 1./sqrt(numHist->GetBinContent(i)); // 1/sqrt(number of total)
effHist->SetBinError(i,err);
}
*/
LogDebug("FourVectorHLTClient")<< "Numerator hist " << numHist->GetName() << endl;
LogDebug("FourVectorHLTClient")<< "Denominator hist " << denHist->GetName() << endl;
LogDebug("FourVectorHLTClient")<< "Booking efficiency histogram path " << newHistPath << endl;
LogDebug("FourVectorHLTClient")<< "Booking efficiency histogram name " << newHistName << endl;
// book eff histos
dbe_->book1D(newHistName.Data(), effHist);
} // end for Stage m
} // end for Stage l
} //end for obj k
vObj.clear();
vObj.push_back(TString("l1Etal1Phi"));
vObj.push_back(TString("onEtaonPhi"));
vObj.push_back(TString("offEtaoffPhi"));
vObj.push_back(TString("mcEtamcPhi"));
for (unsigned int k=0; k<vObj.size();k++) {
for (unsigned int l=0; l<vStage.size();l++) {
for (unsigned int m=0; m<vStage.size();m++) {
TString oldHistPathNum;
TString oldHistPathNumBck;
TString oldHistPathDen;
// not the differeence b/w Num and NumBck, as in OnOff vs OffOn
oldHistPathNum = sourceDir_+"/"+hltPath+"/"+numPathName+"_wrt_"+denPathName+"_"+vObj[k]+vStage[l]+vStage[m];
oldHistPathNumBck = sourceDir_+"/"+hltPath+"/"+numPathName+"_wrt_"+denPathName+"_"+vObj[k]+vStage[m]+vStage[l];
// In the deominator hist name, we don't have any "UM" substrings, so remove them
TString tempDenString = vStage[l];
tempDenString.ReplaceAll("UM","") ;
oldHistPathDen = sourceDir_+"/"+hltPath+"/"+numPathName+"_wrt_"+denPathName+"_"+vObj[k]+tempDenString;
MonitorElement *numME = dbe_->get(oldHistPathNum.Data());
MonitorElement *numMEBck = dbe_->get(oldHistPathNumBck.Data());
MonitorElement *denME = dbe_->get(oldHistPathDen.Data());
LogDebug("FourVectorHLTClient")<< " oldHistPathNum = " << oldHistPathNum << endl;
LogDebug("FourVectorHLTClient")<< " oldHistPathNumBck = " << oldHistPathNumBck << endl;
LogDebug("FourVectorHLTClient")<< " oldHistPathDen = " << oldHistPathDen << endl;
//check if HLTOffline histogram exist
if ( numME && denME ) {
LogDebug("FourVectorHLTClient")<< "DID find NUM and DEN histograms to derive eff " << vStage[m]<<"To"<<vStage[l]<<" using:" <<endl
<< " NUM = " << oldHistPathNum << endl
<< " DEN = " << oldHistPathDen << endl;
}
else {
LogDebug("FourVectorHLTClient")<< "Cannot find NUM and DEN histograms to derive eff " << vStage[m]<<"To"<<vStage[l]<<" using:" <<endl
<< " NUM = " << oldHistPathNum << endl
<< " DEN = " << oldHistPathDen << endl;
if ( numMEBck && denME) {
LogDebug("FourVectorHLTClient")<< "DID find NUM and DEN histograms to derive eff " << vStage[m]<<"To"<<vStage[l]<<" using:" <<endl
<< " NUM = " << oldHistPathNumBck << endl
<< " DEN = " << oldHistPathDen << endl;
numME = numMEBck;
}
else {
LogDebug("FourVectorHLTClient")<< "Cannot find NUM and DEN histograms to derive eff " << vStage[m]<<"To"<<vStage[l]<<" using:" <<endl
<< " NUM = " << oldHistPathNumBck << endl
<< " DEN = " << oldHistPathDen << endl;
continue;
}
}
TH2F* numHist = numME->getTH2F();
TH2F* denHist = denME->getTH2F();
// build names and title for efficiency histogram
TString newHistName = hltPath +"_wrt_" + denPathName +"_"+vObj[k]+"_Eff_"+vStage[m]+"To"+vStage[l];
// if there is "UM", remove it first,then append it to the end of the name
if(newHistName.Contains("UM")) {
newHistName.ReplaceAll("UM","");
newHistName.Append("_UM");
}
TString newHistTitle = numPathName+" given " + denPathName +" "+vObj[k]+" Eff "+vStage[m]+"To"+vStage[l];
// if there is "UM", remove it first,then append it to the end of the name
if(newHistTitle.Contains("UM")) {
newHistTitle.ReplaceAll("UM","");
newHistTitle.Append("_UM");
}
TString newHistPath = currEffFolder+newHistName;
LogDebug("FourVectorHLTClient")<< "Will make efficiency histogram " << newHistPath << endl;
TH2F* effHist = (TH2F*) numHist->Clone(newHistName.Data());
effHist->SetTitle(newHistTitle.Data());
effHist->Sumw2();
//denHist->Sumw2();
//effHist->Divide(numHist,denHist,1.,1.,"B");
effHist->Divide(numHist,denHist,1.,1.);
//reportSummaryMap_ = dbe_->book2D(numHist->Divide(denHist));
LogDebug("FourVectorHLTClient")<< "Booking efficiency histogram path " << newHistPath << endl;
LogDebug("FourVectorHLTClient")<< "Booking efficiency histogram name " << newHistName << endl;
// book eff histos
dbe_->book2D(newHistName.Data(), effHist);
} // end for Stage m
} // end for Stage l
} //end for obj k
} // end for custompathpair
} // end loop over folders for i
hltPathNameColl->Delete();
}
| TH1F* FourVectorHLTClient::get1DHisto | ( | std::string | meName, |
| DQMStore * | dbi | ||
| ) | [private] |
| TProfile* FourVectorHLTClient::get1DProfile | ( | std::string | meName, |
| DQMStore * | dbi | ||
| ) | [private] |
| TH2F* FourVectorHLTClient::get2DHisto | ( | std::string | meName, |
| DQMStore * | dbi | ||
| ) | [private] |
| TProfile2D* FourVectorHLTClient::get2DProfile | ( | std::string | meName, |
| DQMStore * | dbi | ||
| ) | [private] |
| void FourVectorHLTClient::initialize | ( | ) | [private] |
Definition at line 58 of file FourVectorHLTClient.cc.
References dbe_, LogDebug, and cmsCodeRules::cppFunctionSkipper::operator.
{
counterLS_=0;
counterEvt_=0;
// get back-end interface
dbe_ = Service<DQMStore>().operator->();
processname_ = parameters_.getUntrackedParameter<std::string>("processname","HLT");
// base folder for the contents of this job
sourceDir_ = TString(parameters_.getUntrackedParameter<string>("hltSourceDir",""));
// remove trainling "/" from dirname
while(sourceDir_.Last('/') == sourceDir_.Length()-1) {
sourceDir_.Remove(sourceDir_.Length()-1);
}
LogDebug("FourVectorHLTClient")<< "Source dir = " << sourceDir_ << endl;
clientDir_ = TString(parameters_.getUntrackedParameter<string>("hltClientDir",""));
// remove trainling "/" from dirname
while(clientDir_.Last('/') == clientDir_.Length()-1) {
clientDir_.Remove(clientDir_.Length()-1);
}
LogDebug("FourVectorHLTClient")<< "Client dir = " << clientDir_ << endl;
prescaleLS_ = parameters_.getUntrackedParameter<int>("prescaleLS", -1);
LogDebug("FourVectorHLTClient")<< "DQM lumi section prescale = " << prescaleLS_ << " lumi section(s)"<< endl;
prescaleEvt_ = parameters_.getUntrackedParameter<int>("prescaleEvt", -1);
LogDebug("FourVectorHLTClient")<< "DQM event prescale = " << prescaleEvt_ << " events(s)"<< endl;
customEffDir_ = TString(parameters_.getUntrackedParameter<string>("customEffDir","custom-efficiencies"));
LogDebug("FourVectorHLTClient")<< "Custom Efficiencies dir = " << customEffDir_ << endl;
std::vector<edm::ParameterSet> effpaths = parameters_.getParameter<std::vector<edm::ParameterSet> >("effpaths");
std::pair<std::string, std::string> custompathnamepair;
for(std::vector<edm::ParameterSet>::iterator pathconf = effpaths.begin() ; pathconf != effpaths.end();
pathconf++) {
custompathnamepair.first =pathconf->getParameter<std::string>("pathname");
custompathnamepair.second = pathconf->getParameter<std::string>("denompathname");
custompathnamepairs_.push_back(custompathnamepair);
// customdenompathnames_.push_back(pathconf->getParameter<std::string>("denompathname"));
// custompathnames_.push_back(pathconf->getParameter<std::string>("pathname"));
}
pathsSummaryFolder_ = parameters_.getUntrackedParameter ("pathsSummaryFolder",std::string("HLT/FourVector/PathsSummary/"));
pathsSummaryHLTCorrelationsFolder_ = parameters_.getUntrackedParameter ("hltCorrelationsFolder",std::string("HLT/FourVector/PathsSummary/HLT Correlations/"));
pathsSummaryFilterCountsFolder_ = parameters_.getUntrackedParameter ("filterCountsFolder",std::string("HLT/FourVector/PathsSummary/Filters Counts/"));
pathsSummaryFilterEfficiencyFolder_ = parameters_.getUntrackedParameter ("filterEfficiencyFolder",std::string("HLT/FourVector/PathsSummary/Filters Efficiencies/"));
}
| void FourVectorHLTClient::normalizeHLTMatrix | ( | ) | [protected] |
Definition at line 681 of file FourVectorHLTClient.cc.
References MonitorElement::getTH1F(), MonitorElement::getTH2F(), and i.
{
for (unsigned int i =0;i<v_ME_HLTPassPass_.size();i++) {
MonitorElement* ME_HLTPassPass_ = v_ME_HLTPassPass_[i];
MonitorElement* ME_HLTPassPass_Normalized_ = v_ME_HLTPassPass_Normalized_[i];
MonitorElement* ME_HLTPass_Normalized_Any_ = v_ME_HLTPass_Normalized_Any_[i];
if(!ME_HLTPassPass_ || !ME_HLTPassPass_Normalized_ || !ME_HLTPass_Normalized_Any_) return;
float passCount = 0;
unsigned int nBinsX = ME_HLTPassPass_->getTH2F()->GetNbinsX();
unsigned int nBinsY = ME_HLTPassPass_->getTH2F()->GetNbinsY();
for(unsigned int binX = 0; binX < nBinsX+1; binX++) {
passCount = ME_HLTPassPass_->getTH2F()->GetBinContent(binX,binX);
for(unsigned int binY = 0; binY < nBinsY+1; binY++) {
if(passCount != 0) {
// normalize each bin to number of passCount
float normalizedBinContentPassPass = (ME_HLTPassPass_->getTH2F()->GetBinContent(binX,binY))/passCount;
//float normalizedBinContentPassFail = (ME_HLTPassFail_->getTH2F()->GetBinContent(binX,binY))/passCount;
ME_HLTPassPass_Normalized_->getTH2F()->SetBinContent(binX,binY,normalizedBinContentPassPass);
//ME_HLTPassFail_Normalized_->getTH2F()->SetBinContent(binX,binY,normalizedBinContentPassFail);
if(binX == nBinsX) {
ME_HLTPass_Normalized_Any_->getTH1F()->SetBinContent(binY,normalizedBinContentPassPass);
}
}
else {
ME_HLTPassPass_Normalized_->getTH2F()->SetBinContent(binX,binY,0);
//ME_HLTPassFail_Normalized_->getTH2F()->SetBinContent(binX,binY,0);
} // end if else
} // end for binY
} // end for binX
} // end for i
}
TString FourVectorHLTClient::clientDir_ [private] |
Definition at line 80 of file FourVectorHLTClient.h.
int FourVectorHLTClient::counterEvt_ [private] |
counter
Definition at line 100 of file FourVectorHLTClient.h.
int FourVectorHLTClient::counterLS_ [private] |
Definition at line 99 of file FourVectorHLTClient.h.
TString FourVectorHLTClient::customEffDir_ [private] |
Definition at line 81 of file FourVectorHLTClient.h.
std::vector<std::pair<std::string, std::string> > FourVectorHLTClient::custompathnamepairs_ [private] |
Definition at line 88 of file FourVectorHLTClient.h.
DQMStore* FourVectorHLTClient::dbe_ [private] |
Definition at line 78 of file FourVectorHLTClient.h.
Definition at line 108 of file FourVectorHLTClient.h.
std::vector<TString> FourVectorHLTClient::hltMEName [private] |
Definition at line 84 of file FourVectorHLTClient.h.
std::vector<MonitorElement*> FourVectorHLTClient::hltMEs [private] |
Definition at line 87 of file FourVectorHLTClient.h.
std::vector<TString> FourVectorHLTClient::hltPathName [private] |
Definition at line 85 of file FourVectorHLTClient.h.
MonitorElement* FourVectorHLTClient::klmgrvTest_ [private] |
Definition at line 113 of file FourVectorHLTClient.h.
int FourVectorHLTClient::nChannels [private] |
prescale on number of events
Definition at line 103 of file FourVectorHLTClient.h.
Definition at line 75 of file FourVectorHLTClient.h.
std::string FourVectorHLTClient::pathsSummaryFilterCountsFolder_ [private] |
Definition at line 96 of file FourVectorHLTClient.h.
std::string FourVectorHLTClient::pathsSummaryFilterEfficiencyFolder_ [private] |
Definition at line 97 of file FourVectorHLTClient.h.
std::string FourVectorHLTClient::pathsSummaryFolder_ [private] |
Definition at line 94 of file FourVectorHLTClient.h.
std::string FourVectorHLTClient::pathsSummaryHLTCorrelationsFolder_ [private] |
Definition at line 95 of file FourVectorHLTClient.h.
int FourVectorHLTClient::prescaleEvt_ [private] |
units of lumi sections
Definition at line 102 of file FourVectorHLTClient.h.
int FourVectorHLTClient::prescaleLS_ [private] |
counter
Definition at line 101 of file FourVectorHLTClient.h.
std::string FourVectorHLTClient::processname_ [private] |
Definition at line 82 of file FourVectorHLTClient.h.
Float_t FourVectorHLTClient::reportSummary [private] |
Definition at line 104 of file FourVectorHLTClient.h.
Definition at line 109 of file FourVectorHLTClient.h.
MonitorElement* FourVectorHLTClient::reportSummaryContent_[20] [private] |
Definition at line 110 of file FourVectorHLTClient.h.
Definition at line 111 of file FourVectorHLTClient.h.
TString FourVectorHLTClient::sourceDir_ [private] |
Definition at line 79 of file FourVectorHLTClient.h.
Float_t FourVectorHLTClient::summaryContent[20] [private] |
Definition at line 106 of file FourVectorHLTClient.h.
Float_t FourVectorHLTClient::summarySum [private] |
Definition at line 105 of file FourVectorHLTClient.h.
MonitorElement* FourVectorHLTClient::testHLTEff_ [private] |
Definition at line 112 of file FourVectorHLTClient.h.
std::vector<MonitorElement*> FourVectorHLTClient::v_ME_HLTPass_Normalized_Any_ [private] |
Definition at line 92 of file FourVectorHLTClient.h.
std::vector<MonitorElement*> FourVectorHLTClient::v_ME_HLTPassPass_ [private] |
Definition at line 90 of file FourVectorHLTClient.h.
std::vector<MonitorElement*> FourVectorHLTClient::v_ME_HLTPassPass_Normalized_ [private] |
Definition at line 91 of file FourVectorHLTClient.h.
1.7.3