|
|
|
#include <L1TRate.h>
Public Member Functions | |
| L1TRate (const edm::ParameterSet &ps) | |
| virtual | ~L1TRate () |
Protected Member Functions | |
| void | analyze (const edm::Event &e, const edm::EventSetup &c) |
| void | beginJob () |
| virtual void | beginLuminosityBlock (edm::LuminosityBlock const &lumiBlock, edm::EventSetup const &c) |
| void | beginRun (const edm::Run &run, const edm::EventSetup &iSetup) |
| void | endJob () |
| virtual void | endLuminosityBlock (edm::LuminosityBlock const &lumiBlock, edm::EventSetup const &c) |
| void | endRun (const edm::Run &run, const edm::EventSetup &iSetup) |
Private Member Functions | |
| bool | getXSexFitsOMDS (const edm::ParameterSet &ps) |
| bool | getXSexFitsPython (const edm::ParameterSet &ps) |
Private Attributes | |
| DQMStore * | dbe |
| std::map< TString, int > | m_algoBit |
| MonitorElement * | m_ErrorMonitor |
| std::map< std::string, bool > | m_inputCategories |
| edm::InputTag | m_l1GtDataDaqInputTag |
| const std::vector< std::vector < int > > * | m_listsPrescaleFactors |
| std::map< int, double > | m_lsLuminosity |
| std::map< int, int > | m_lsPrescaleIndex |
| std::map< int, std::map < TString, double > > | m_lsRates |
| int | m_lsShiftGTRates |
| int | m_maxNbins |
| std::string | m_outputFile |
| edm::ParameterSet | m_parameters |
| int | m_refPrescaleSet |
| edm::InputTag | m_scalersSource |
| std::map< std::string, std::string > | m_selectedTriggers |
| std::map< TString, TF1 * > | m_templateFunctions |
| bool | m_verbose |
| std::map< TString, MonitorElement * > | m_xSecObservedToExpected |
| std::map< TString, MonitorElement * > | m_xSecVsInstLumi |
| L1TRate::L1TRate | ( | const edm::ParameterSet & | ps | ) |
Definition at line 41 of file L1TRate.cc.
References gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), NULL, and cppFunctionSkipper::operator.
{
m_maxNbins = 2500; // Maximum LS for each run (for binning purposes)
m_parameters = ps;
// Mapping parameter input variables
m_scalersSource = m_parameters.getParameter <InputTag>("inputTagScalersResults");
m_l1GtDataDaqInputTag = m_parameters.getParameter <InputTag>("inputTagL1GtDataDaq");
m_verbose = m_parameters.getUntrackedParameter<bool> ("verbose",false);
m_refPrescaleSet = m_parameters.getParameter <int> ("refPrescaleSet");
m_lsShiftGTRates = m_parameters.getUntrackedParameter<int> ("lsShiftGTRates",0);
// Getting which categories to monitor
ParameterSet Categories = ps.getParameter<ParameterSet>("categories");
m_inputCategories["Mu"] = Categories.getUntrackedParameter<bool>("Mu");
m_inputCategories["EG"] = Categories.getUntrackedParameter<bool>("EG");
m_inputCategories["IsoEG"] = Categories.getUntrackedParameter<bool>("IsoEG");
m_inputCategories["Jet"] = Categories.getUntrackedParameter<bool>("Jet");
m_inputCategories["CenJet"] = Categories.getUntrackedParameter<bool>("CenJet");
m_inputCategories["ForJet"] = Categories.getUntrackedParameter<bool>("ForJet");
m_inputCategories["TauJet"] = Categories.getUntrackedParameter<bool>("TauJet");
m_inputCategories["ETM"] = Categories.getUntrackedParameter<bool>("ETM");
m_inputCategories["ETT"] = Categories.getUntrackedParameter<bool>("ETT");
m_inputCategories["HTT"] = Categories.getUntrackedParameter<bool>("HTT");
m_inputCategories["HTM"] = Categories.getUntrackedParameter<bool>("HTM");
// Inicializing Variables
dbe = NULL;
if (ps.getUntrackedParameter < bool > ("dqmStore", false)) {
dbe = Service < DQMStore > ().operator->();
dbe->setVerbose(0);
}
// What to do if we want our output to be saved to a external file
m_outputFile = ps.getUntrackedParameter < string > ("outputFile", "");
if (m_outputFile.size() != 0) {
cout << "L1T Monitoring histograms will be saved to " << m_outputFile.c_str() << endl;
}
bool disable = ps.getUntrackedParameter < bool > ("disableROOToutput", false);
if (disable) {m_outputFile = "";}
if (dbe != NULL) {dbe->setCurrentFolder("L1T/L1TRate");}
}
| L1TRate::~L1TRate | ( | ) | [virtual] |
Definition at line 90 of file L1TRate.cc.
{}
| void L1TRate::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
Implements edm::EDAnalyzer.
Definition at line 334 of file L1TRate.cc.
References gather_cfg::cout, edm::Event::getByLabel(), Level1TriggerRates::gtAlgoCountsRate(), i, edm::EventBase::id(), edm::HandleBase::isValid(), edm::EventID::luminosityBlock(), and edm::EventID::run().
{
edm::Handle<L1GlobalTriggerReadoutRecord> gtReadoutRecordData;
edm::Handle<Level1TriggerScalersCollection> triggerScalers;
edm::Handle<LumiScalersCollection> colLScal;
iEvent.getByLabel(m_l1GtDataDaqInputTag, gtReadoutRecordData);
iEvent.getByLabel(m_scalersSource , colLScal);
iEvent.getByLabel(m_scalersSource , triggerScalers);
// Integers
int EventRun = iEvent.id().run();
unsigned int eventLS = iEvent.id().luminosityBlock();
// Getting the trigger trigger rates from GT and buffering it
if(triggerScalers.isValid()){
Level1TriggerScalersCollection::const_iterator itL1TScalers = triggerScalers->begin();
Level1TriggerRates trigRates(*itL1TScalers,EventRun);
int gtLS = (*itL1TScalers).lumiSegmentNr()+m_lsShiftGTRates;
// If we haven't got the data from this LS yet get it
if(m_lsRates.find(gtLS)==m_lsRates.end()){
if (m_verbose) {cout << "[L1TRate:] Buffering GT Rates for LS=" << gtLS << endl;}
map<TString,double> bufferRate;
// Buffer the rate informations for all selected bits
for(map<string,string>::const_iterator i=m_selectedTriggers.begin() ; i!=m_selectedTriggers.end() ; i++){
string tTrigger = (*i).second;
// If trigger name is defined we store the rate
if(tTrigger != "Undefined"){
unsigned int trigBit = m_algoBit[tTrigger];
double trigRate = trigRates.gtAlgoCountsRate()[trigBit];
bufferRate[tTrigger] = trigRate;
}
}
m_lsRates[gtLS] = bufferRate;
}
}
// Getting from the SCAL the luminosity information and buffering it
if(colLScal.isValid() && colLScal->size()){
LumiScalersCollection::const_iterator itLScal = colLScal->begin();
unsigned int scalLS = itLScal->sectionNumber();
// If we haven't got the data from this SCAL LS yet get it
if(m_lsLuminosity.find(scalLS)==m_lsLuminosity.end()){
if (m_verbose) {cout << "[L1TRate:] Buffering SCAL-HF Lumi for LS=" << scalLS << endl;}
double instLumi = itLScal->instantLumi(); // Getting Instant Lumi from HF (via SCAL)
double deadTimeNormHF = itLScal->deadTimeNormalization(); // Getting Dead Time Normalization from HF (via SCAL)
// If HF Dead Time Corrections is requested we apply it
// NOTE: By default this is assumed false since for now WbM fits do NOT assume this correction
if(m_parameters.getUntrackedParameter<bool>("useHFDeadTimeNormalization",false)){
// Protecting for deadtime = 0
if(deadTimeNormHF==0){instLumi = 0;}
else {instLumi = instLumi/deadTimeNormHF;}
}
// Buffering the luminosity information
m_lsLuminosity[scalLS]=instLumi;
}
}
// Getting the prescale index used when this event was triggered
if(gtReadoutRecordData.isValid()){
// If we haven't got the data from this LS yet get it
if(m_lsPrescaleIndex.find(eventLS)==m_lsPrescaleIndex.end()){
if (m_verbose) {cout << "[L1TRate:] Buffering Prescale Index for LS=" << eventLS << endl;}
// Getting Final Decision Logic (FDL) Data from GT
const vector<L1GtFdlWord>& gtFdlVectorData = gtReadoutRecordData->gtFdlVector();
// Getting the index for the fdl data for this event
int indexFDL=0;
for(unsigned int i=0; i<gtFdlVectorData.size(); i++){
if(gtFdlVectorData[i].bxInEvent()==0){indexFDL=i; break;}
}
int CurrentPrescalesIndex = gtFdlVectorData[indexFDL].gtPrescaleFactorIndexAlgo();
m_lsPrescaleIndex[eventLS] = CurrentPrescalesIndex;
}
}
}
| void L1TRate::beginJob | ( | void | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 93 of file L1TRate.cc.
References gather_cfg::cout, cppFunctionSkipper::operator, DQMStore::rmdir(), and DQMStore::setCurrentFolder().
| void L1TRate::beginLuminosityBlock | ( | edm::LuminosityBlock const & | lumiBlock, |
| edm::EventSetup const & | c | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 244 of file L1TRate.cc.
References gather_cfg::cout, edm::LuminosityBlockBase::id(), and edm::LuminosityBlockID::luminosityBlock().
| void L1TRate::beginRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | iSetup | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 122 of file L1TRate.cc.
References gather_cfg::cout, edm::EventSetup::get(), L1TMenuHelper::getEtaRangeByAlias(), L1TMenuHelper::getLUSOTrigger(), L1TMenuHelper::getPrescaleByAlias(), L1TMenuHelper::getQualityAlias(), L1GtTriggerMenu::gtAlgorithmMap(), L1GtPrescaleFactors::gtPrescaleFactors(), i, relval_steps::menu, and edm::ESHandle< T >::product().
{
if (m_verbose) {cout << "[L1TRate:] Called beginRun." << endl;}
ESHandle<L1GtTriggerMenu> menuRcd;
ESHandle<L1GtPrescaleFactors> l1GtPfAlgo;
iSetup.get<L1GtTriggerMenuRcd>() .get(menuRcd);
iSetup.get<L1GtPrescaleFactorsAlgoTrigRcd>().get(l1GtPfAlgo);
const L1GtTriggerMenu* menu = menuRcd .product();
const L1GtPrescaleFactors* m_l1GtPfAlgo = l1GtPfAlgo.product();
// Initializing DQM Monitor Elements
dbe->setCurrentFolder("L1T/L1TRate");
m_ErrorMonitor = dbe->book1D("ErrorMonitor", "ErrorMonitor",5,0,5);
m_ErrorMonitor->setBinLabel(1,"WARNING_DB_CONN_FAILED"); // Errors from L1TOMDSHelper
m_ErrorMonitor->setBinLabel(2,"WARNING_DB_QUERY_FAILED"); // Errors from L1TOMDSHelper
m_ErrorMonitor->setBinLabel(3,"WARNING_DB_INCORRECT_NBUNCHES"); // Errors from L1TOMDSHelper
m_ErrorMonitor->setBinLabel(4,"WARNING_PY_MISSING_FIT");
m_ErrorMonitor->setBinLabel(5,"UNKNOWN");
// Retriving the list of prescale sets
m_listsPrescaleFactors = &(m_l1GtPfAlgo->gtPrescaleFactors());
// Getting Lowest Prescale Single Object Triggers from the menu
L1TMenuHelper myMenuHelper = L1TMenuHelper(iSetup);
m_selectedTriggers = myMenuHelper.getLUSOTrigger(m_inputCategories,m_refPrescaleSet);
//-> Getting template fits for the algLo cross sections
int srcAlgoXSecFit = m_parameters.getParameter<int>("srcAlgoXSecFit");
if (srcAlgoXSecFit == 0){getXSexFitsOMDS (m_parameters);}
else if(srcAlgoXSecFit == 1){getXSexFitsPython(m_parameters);}
for (CItAlgo algo = menu->gtAlgorithmMap().begin(); algo!=menu->gtAlgorithmMap().end(); ++algo){
m_algoBit[(algo->second).algoAlias()] = (algo->second).algoBitNumber();
}
double minInstantLuminosity = m_parameters.getParameter<double>("minInstantLuminosity");
double maxInstantLuminosity = m_parameters.getParameter<double>("maxInstantLuminosity");
// Initializing DQM Monitor Elements
for(map<string,string>::const_iterator i=m_selectedTriggers.begin() ; i!=m_selectedTriggers.end() ; i++){
TString tCategory = (*i).first;
TString tTrigger = (*i).second;
TString tErrorMessage = "";
TF1* tTestFunction;
if(tTrigger != "Undefined" && m_templateFunctions.find(tTrigger) != m_templateFunctions.end()){
tTestFunction = m_templateFunctions[tTrigger];
}
else if(tTrigger == "Undefined"){
TString tFunc = "-1";
tTestFunction = new TF1("FitParametrization_"+tTrigger,tFunc,0,double(m_maxNbins)-0.5);
}
else if(m_templateFunctions.find(tTrigger) == m_templateFunctions.end()){
TString tFunc = "-1";
tTestFunction = new TF1("FitParametrization_"+tTrigger,tFunc,0,double(m_maxNbins)-0.5);
tErrorMessage = " (Undefined Test Function)";
}
else{
TString tFunc = "-1";
tTestFunction = new TF1("FitParametrization_"+tTrigger,tFunc,0,double(m_maxNbins)-0.5);
}
if(tTrigger != "Undefined"){
if(myMenuHelper.getPrescaleByAlias(tCategory,tTrigger) != 1){
tErrorMessage += " WARNING: Default Prescale = ";
tErrorMessage += myMenuHelper.getPrescaleByAlias(tCategory,tTrigger);
}
if (tCategory == "Mu" && myMenuHelper.getEtaRangeByAlias(tCategory,tTrigger) != 4294967295){
tErrorMessage += " WARNING: Eta Range = ";
tErrorMessage += myMenuHelper.getEtaRangeByAlias(tCategory,tTrigger);
}
else if(tCategory == "EG" && myMenuHelper.getEtaRangeByAlias(tCategory,tTrigger) != 32639){
tErrorMessage += " WARNING: Eta Range = ";
tErrorMessage += myMenuHelper.getEtaRangeByAlias(tCategory,tTrigger);
}
else if(tCategory == "IsoEG" && myMenuHelper.getEtaRangeByAlias(tCategory,tTrigger) != 32639){
tErrorMessage += " WARNING: Eta Range = ";
tErrorMessage += myMenuHelper.getEtaRangeByAlias(tCategory,tTrigger);
}
if(tCategory == "Mu" && myMenuHelper.getQualityAlias(tCategory,tTrigger) != 240){
tErrorMessage += " WARNING: Quality = ";
tErrorMessage += myMenuHelper.getQualityAlias(tCategory,tTrigger);
}
}
dbe->setCurrentFolder("L1T/L1TRate/TriggerCrossSections");
m_xSecVsInstLumi[tTrigger] = dbe->bookProfile(tCategory,
"Cross Sec. vs Inst. Lumi Algo: "+tTrigger+tErrorMessage,
m_maxNbins,
minInstantLuminosity,
maxInstantLuminosity,0,500);
m_xSecVsInstLumi[tTrigger] ->setAxisTitle("Instantaneous Luminosity [10^{30}cm^{-2}s^{-1}]" ,1);
m_xSecVsInstLumi[tTrigger] ->setAxisTitle("Algorithm #sigma [#mu b]" ,2);
m_xSecVsInstLumi[tTrigger] ->getTProfile()->GetListOfFunctions()->Add(tTestFunction);
m_xSecVsInstLumi[tTrigger] ->getTProfile()->SetMarkerStyle(23);
dbe->setCurrentFolder("L1T/L1TRate/Certification");
m_xSecObservedToExpected[tTrigger] = dbe->book1D(tCategory, "Algo: "+tTrigger+tErrorMessage,m_maxNbins,-0.5,double(m_maxNbins)-0.5);
m_xSecObservedToExpected[tTrigger] ->setAxisTitle("Lumi Section" ,1);
m_xSecObservedToExpected[tTrigger] ->setAxisTitle("#sigma_{obs} / #sigma_{exp}" ,2);
}
}
| void L1TRate::endJob | ( | void | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 109 of file L1TRate.cc.
References gather_cfg::cout.
{
if (m_verbose) {cout << "[L1TRate:] Called endJob." << endl;}
if (m_outputFile.size() != 0 && dbe)
dbe->save(m_outputFile);
return;
}
| void L1TRate::endLuminosityBlock | ( | edm::LuminosityBlock const & | lumiBlock, |
| edm::EventSetup const & | c | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 251 of file L1TRate.cc.
References gather_cfg::cout, getTProfile(), i, edm::LuminosityBlockBase::id(), python::rootplot::utilities::ls(), fjr2json::lumi, and edm::LuminosityBlockID::luminosityBlock().
{
int eventLS = lumiBlock.id().luminosityBlock();
if (m_verbose) {cout << "[L1TRate:] Called endLuminosityBlock at LS=" << eventLS << endl;}
// We can certify LS -1 since we should have available:
// gt rates: (current LS)-1
// prescale: current LS
// lumi : current LS
//eventLS--;
// Checking if all necessary quantities are defined for our calculations
bool isDefRate,isDefLumi,isDefPrescaleIndex;
map<TString,double>* rates=0;
double lumi=0;
int prescalesIndex=0;
// Reseting MonitorElements so we can refill them
for(map<string,string>::const_iterator i=m_selectedTriggers.begin() ; i!=m_selectedTriggers.end() ; i++){
string tTrigger = (*i).second;
m_xSecObservedToExpected[tTrigger]->getTH1()->Reset("ICE");
m_xSecVsInstLumi [tTrigger]->getTH1()->Reset("ICE");
}
for(map<int,map<TString,double> >::iterator i=m_lsRates.begin() ; i!=m_lsRates.end() ; i++){
unsigned int ls = (*i).first;
rates = &(*i).second;
isDefRate=true;
if(m_lsLuminosity.find(ls)==m_lsLuminosity.end()){isDefLumi=false;}
else{
isDefLumi=true;
lumi=m_lsLuminosity[ls];
}
if(m_lsPrescaleIndex.find(ls)==m_lsPrescaleIndex.end()){isDefPrescaleIndex=false;}
else{
isDefPrescaleIndex=true;
prescalesIndex=m_lsPrescaleIndex[ls];
}
if(isDefRate && isDefLumi && isDefPrescaleIndex){
const vector<int>& currentPrescaleFactors = (*m_listsPrescaleFactors).at(prescalesIndex);
for(map<string,string>::const_iterator i=m_selectedTriggers.begin() ; i!=m_selectedTriggers.end() ; i++){
string tTrigger = (*i).second;
TF1* tTestFunction = (TF1*) m_xSecVsInstLumi[tTrigger]->getTProfile()->GetListOfFunctions()->First();
// If trigger name is defined we get the rate fit parameters
if(tTrigger != "Undefined"){
unsigned int trigBit = m_algoBit[tTrigger];
double trigPrescale = currentPrescaleFactors[trigBit];
double trigRate = (*rates)[tTrigger];
if(lumi!=0 && trigPrescale!=0 && trigRate!=0){
double AlgoXSec = (trigPrescale*trigRate)/lumi;
double TemplateFunctionValue = tTestFunction->Eval(lumi);
// Checking against Template function
int ibin = m_xSecObservedToExpected[tTrigger]->getTH1()->FindBin(ls);
m_xSecObservedToExpected[tTrigger]->setBinContent(ibin,AlgoXSec/TemplateFunctionValue);
m_xSecVsInstLumi [tTrigger]->Fill(lumi,AlgoXSec);
if(m_verbose){cout<<"[L1TRate:] ls="<<ls<<" Algo="<<tTrigger<<" XSec="<<AlgoXSec<<" Test="<<AlgoXSec/TemplateFunctionValue<<endl;}
}
else{
int ibin = m_xSecObservedToExpected[tTrigger]->getTH1()->FindBin(ls);
m_xSecObservedToExpected[tTrigger]->setBinContent(ibin,0.000001);
if(m_verbose){cout << "[L1TRate:] Algo="<< tTrigger<< " XSec=Failed" << endl;}
}
}
}
}
}
}
| void L1TRate::endRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | iSetup | ||
| ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 239 of file L1TRate.cc.
References gather_cfg::cout.
| bool L1TRate::getXSexFitsOMDS | ( | const edm::ParameterSet & | ps | ) | [private] |
Definition at line 438 of file L1TRate.cc.
References a, WbMTriggerXSecFit::bitName, L1TOMDSHelper::connect(), WbMTriggerXSecFit::fitFunction, edm::ParameterSet::getParameter(), L1TOMDSHelper::NO_ERROR, WbMTriggerXSecFit::p0, WbMTriggerXSecFit::p1, WbMTriggerXSecFit::p2, and WbMTriggerXSecFit::pm1.
{
bool noError = true;
string oracleDB = ps.getParameter<string>("oracleDB");
string pathCondDB = ps.getParameter<string>("pathCondDB");
L1TOMDSHelper myOMDSHelper;
int conError;
myOMDSHelper.connect(oracleDB,pathCondDB,conError);
map<string,WbMTriggerXSecFit> wbmFits;
if(conError == L1TOMDSHelper::NO_ERROR){
int errorRetrive;
wbmFits = myOMDSHelper.getWbMAlgoXsecFits(errorRetrive);
// Filling errors if they exist
if(errorRetrive != L1TOMDSHelper::NO_ERROR){
noError = false;
string eName = myOMDSHelper.enumToStringError(errorRetrive);
m_ErrorMonitor->Fill(eName);
}
}else{
noError = false;
string eName = myOMDSHelper.enumToStringError(conError);
m_ErrorMonitor->Fill(eName);
}
double minInstantLuminosity = m_parameters.getParameter<double>("minInstantLuminosity");
double maxInstantLuminosity = m_parameters.getParameter<double>("maxInstantLuminosity");
// Getting rate fit parameters for all input triggers
for(map<string,string>::const_iterator a=m_selectedTriggers.begin() ; a!=m_selectedTriggers.end() ; a++){
string tTrigger = (*a).second;
// If trigger name is defined we get the rate fit parameters
if(tTrigger != "Undefined"){
if(wbmFits.find(tTrigger) != wbmFits.end()){
WbMTriggerXSecFit tWbMParameters = wbmFits[tTrigger];
vector<double> tParameters;
tParameters.push_back(tWbMParameters.pm1);
tParameters.push_back(tWbMParameters.p0);
tParameters.push_back(tWbMParameters.p1);
tParameters.push_back(tWbMParameters.p2);
// Retriving and populating the m_templateFunctions array
m_templateFunctions[tTrigger] = new TF1("FitParametrization_"+tWbMParameters.bitName,
tWbMParameters.fitFunction,
minInstantLuminosity,maxInstantLuminosity);
m_templateFunctions[tTrigger] ->SetParameters(&tParameters[0]);
m_templateFunctions[tTrigger] ->SetLineWidth(1);
m_templateFunctions[tTrigger] ->SetLineColor(kRed);
}else{noError = false;}
}
}
return noError;
}
| bool L1TRate::getXSexFitsPython | ( | const edm::ParameterSet & | ps | ) | [private] |
Definition at line 513 of file L1TRate.cc.
References a, b, and edm::ParameterSet::getParameter().
{
// error meaning
bool noError = true;
// Getting fit parameters
std::vector<edm::ParameterSet> m_fitParameters = ps.getParameter< vector<ParameterSet> >("fitParameters");
double minInstantLuminosity = m_parameters.getParameter<double>("minInstantLuminosity");
double maxInstantLuminosity = m_parameters.getParameter<double>("maxInstantLuminosity");
// Getting rate fit parameters for all input triggers
for(map<string,string>::const_iterator a=m_selectedTriggers.begin() ; a!=m_selectedTriggers.end() ; a++){
string tTrigger = (*a).second;
// If trigger name is defined we get the rate fit parameters
if(tTrigger != "Undefined"){
bool foundFit = false;
for(unsigned int b=0 ; b<m_fitParameters.size() ; b++){
if(tTrigger == m_fitParameters[b].getParameter<string>("AlgoName")){
TString tAlgoName = m_fitParameters[b].getParameter< string > ("AlgoName");
TString tTemplateFunction = m_fitParameters[b].getParameter< string > ("TemplateFunction");
vector<double> tParameters = m_fitParameters[b].getParameter< vector<double> >("Parameters");
// Retriving and populating the m_templateFunctions array
m_templateFunctions[tTrigger] = new TF1("FitParametrization_"+tAlgoName,tTemplateFunction,
minInstantLuminosity,maxInstantLuminosity);
m_templateFunctions[tTrigger] ->SetParameters(&tParameters[0]);
m_templateFunctions[tTrigger] ->SetLineWidth(1);
m_templateFunctions[tTrigger] ->SetLineColor(kRed);
foundFit = true;
break;
}
if(!foundFit){
noError = false;
string eName = "WARNING_PY_MISSING_FIT";
m_ErrorMonitor->Fill(eName);
}
}
}
}
return noError;
}
DQMStore* L1TRate::dbe [private] |
std::map<TString,int> L1TRate::m_algoBit [private] |
MonitorElement* L1TRate::m_ErrorMonitor [private] |
std::map<std::string,bool> L1TRate::m_inputCategories [private] |
edm::InputTag L1TRate::m_l1GtDataDaqInputTag [private] |
const std::vector< std::vector<int> >* L1TRate::m_listsPrescaleFactors [private] |
std::map<int,double> L1TRate::m_lsLuminosity [private] |
std::map<int,int> L1TRate::m_lsPrescaleIndex [private] |
std::map<int,std::map<TString,double> > L1TRate::m_lsRates [private] |
int L1TRate::m_lsShiftGTRates [private] |
int L1TRate::m_maxNbins [private] |
std::string L1TRate::m_outputFile [private] |
edm::ParameterSet L1TRate::m_parameters [private] |
int L1TRate::m_refPrescaleSet [private] |
edm::InputTag L1TRate::m_scalersSource [private] |
std::map<std::string,std::string> L1TRate::m_selectedTriggers [private] |
std::map<TString,TF1*> L1TRate::m_templateFunctions [private] |
bool L1TRate::m_verbose [private] |
std::map<TString,MonitorElement*> L1TRate::m_xSecObservedToExpected [private] |
std::map<TString,MonitorElement*> L1TRate::m_xSecVsInstLumi [private] |
1.7.3