TagProbeFitTreeAnalyzer Class Reference

Inheritance diagram for TagProbeFitTreeAnalyzer:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override
void	calculateEfficiency (string name, const edm::ParameterSet &pset)
void	endRun (const edm::Run &run, const edm::EventSetup &setup) override
	TagProbeFitTreeAnalyzer (const edm::ParameterSet &pset)
	~TagProbeFitTreeAnalyzer () override
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
TagProbeFitter	fitter
unsigned int	split_mode

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 14 of file TagProbeFitTreeAnalyzer.cc.

Constructor & Destructor Documentation

TagProbeFitTreeAnalyzer::TagProbeFitTreeAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 26 of file TagProbeFitTreeAnalyzer.cc.

References TagProbeFitter::addCategory(), TagProbeFitter::addExpression(), TagProbeFitter::addPdf(), TagProbeFitter::addThresholdCategory(), TagProbeFitter::addVariable(), writedatasetfile::args, calculateEfficiency(), eostools::cat(), myMessageLogger_cff::categories, egammaForCoreTracking_cff::cuts, electronProducer_cfi::efficiencies, cppFunctionSkipper::exception, Exception, edm::ParameterSet::existsAs(), jets_cff::expr, cmsPerfPublish::fail(), fitter, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNamesForType(), edm::ParameterSet::getUntrackedParameter(), dataset::name, TagProbeFitter::setBinnedFit(), TagProbeFitter::setBinsForMassPlots(), TagProbeFitter::setQuiet(), TagProbeFitter::setSaveDistributionsPlot(), TagProbeFitter::setSplitMode(), TagProbeFitter::setWeightVar(), split_mode, AlCaHLTBitMon_QueryRunRegistry::string, JetChargeProducer_cfi::var, and objects.autophobj::variables.

                                                                           :
  fitter( pset.getParameter<vector<string> >("InputFileNames"),
          pset.getParameter<string>("InputDirectoryName"),
          pset.getParameter<string>("InputTreeName"),
          pset.getParameter<string>("OutputFileName"),
          pset.existsAs<unsigned int>("NumCPU")?pset.getParameter<unsigned int>("NumCPU"):1,
          pset.existsAs<bool>("SaveWorkspace")?pset.getParameter<bool>("SaveWorkspace"):false,
      pset.existsAs<bool>("floatShapeParameters")?pset.getParameter<bool>("floatShapeParameters"):true,
      pset.existsAs<vector<string> >("fixVars")?pset.getParameter<vector<string> >("fixVars"):vector<string>()
      ),
  split_mode( pset.existsAs<unsigned int>("SplitMode")?pset.getParameter<unsigned int>("SplitMode"):0 )
{
  fitter.setQuiet(pset.getUntrackedParameter("Quiet",false));
  fitter.setSplitMode( split_mode );

  if (pset.existsAs<bool>("binnedFit")) {
    bool binned = pset.getParameter<bool>("binnedFit");
    fitter.setBinnedFit(binned, binned ? pset.getParameter<uint32_t>("binsForFit") : 0);
  } else if (pset.existsAs<uint32_t>("binsForMassPlots")) {
    fitter.setBinsForMassPlots(pset.getParameter<uint32_t>("binsForMassPlots"));
  }

  if (pset.existsAs<bool>("saveDistributionsPlot")) {
    fitter.setSaveDistributionsPlot(pset.getParameter<bool>("saveDistributionsPlot"));
  }
  if (pset.existsAs<std::string>("WeightVariable")) {
    fitter.setWeightVar(pset.getParameter<std::string>("WeightVariable"));
  }
  const ParameterSet variables = pset.getParameter<ParameterSet>("Variables");
  vector<string> variableNames = variables.getParameterNamesForType<vector<string> >();
  for (vector<string>::const_iterator name = variableNames.begin(); name != variableNames.end(); name++) {
    vector<string> var = variables.getParameter<vector<string> >(*name);
    double lo, hi;
    if(var.size()>=4 && !(istringstream(var[1])>>lo).fail() && !(istringstream(var[2])>>hi).fail()){
      fitter.addVariable(*name, var[0], lo, hi, var[3]);
    }else{
      LogError("TagProbeFitTreeAnalyzer")<<"Could not create variable: "<<*name<<
      ". Example: pt = cms.vstring(\"Probe pT\", \"1.0\", \"100.0\", \"GeV/c\") ";
    }
  }

  const ParameterSet categories = pset.getParameter<ParameterSet>("Categories");
  vector<string> categoryNames = categories.getParameterNamesForType<vector<string> >();
  for (vector<string>::const_iterator name = categoryNames.begin(); name != categoryNames.end(); name++) {
    vector<string> cat = categories.getParameter<vector<string> >(*name);
    if(cat.size()==2){
      fitter.addCategory(*name, cat[0], cat[1]);
    }else{
      LogError("TagProbeFitTreeAnalyzer")<<"Could not create category: "<<*name<<
      ". Example: mcTrue = cms.vstring(\"MC True\", \"dummy[true=1,false=0]\") ";
    }
  }

  if (pset.existsAs<ParameterSet>("Expressions")) {
    const ParameterSet exprs = pset.getParameter<ParameterSet>("Expressions");
    vector<string> exprNames = exprs.getParameterNamesForType<vector<string> >();
    for (vector<string>::const_iterator name = exprNames.begin(); name != exprNames.end(); name++) {
        vector<string> expr = exprs.getParameter<vector<string> >(*name);
        if(expr.size()>=2){
            vector<string> args(expr.begin()+2,expr.end());
            fitter.addExpression(*name, expr[0], expr[1], args);
        }else{
            LogError("TagProbeFitTreeAnalyzer")<<"Could not create expr: "<<*name<<
                ". Example: qop = cms.vstring(\"qOverP\", \"charge/p\", \"charge\", \"p\") ";
        }
    }
  }


  if (pset.existsAs<ParameterSet>("Cuts")) {
    const ParameterSet cuts = pset.getParameter<ParameterSet>("Cuts");
    vector<string> cutNames = cuts.getParameterNamesForType<vector<string> >();
    for (vector<string>::const_iterator name = cutNames.begin(); name != cutNames.end(); name++) {
        vector<string> cat = cuts.getParameter<vector<string> >(*name);
        if(cat.size()==3){
            fitter.addThresholdCategory(*name, cat[0], cat[1], atof(cat[2].c_str()));
        }else{
            LogError("TagProbeFitTreeAnalyzer")<<"Could not create cut: "<<*name<<
                ". Example: matched = cms.vstring(\"Matched\", \"deltaR\", \"0.5\") ";
        }
    }
  }

  if(pset.existsAs<ParameterSet>("PDFs")){
    const ParameterSet pdfs = pset.getParameter<ParameterSet>("PDFs");
    vector<string> pdfNames = pdfs.getParameterNamesForType<vector<string> >();
    for (vector<string>::const_iterator name = pdfNames.begin(); name != pdfNames.end(); name++) {
      vector<string> pdf = pdfs.getParameter<vector<string> >(*name);
      fitter.addPdf(*name, pdf);
    }
  }

  const ParameterSet efficiencies = pset.getParameter<ParameterSet>("Efficiencies");
  vector<string> efficiencyNames = efficiencies.getParameterNamesForType<ParameterSet>();
  for (vector<string>::const_iterator name = efficiencyNames.begin(); name != efficiencyNames.end(); name++) {
    try {
        calculateEfficiency(*name, efficiencies.getParameter<ParameterSet>(*name));
    } catch (std::exception &ex) {
        throw cms::Exception("Error", ex.what());
    }
  }
}

TagProbeFitTreeAnalyzer::~TagProbeFitTreeAnalyzer ( )

inlineoverride

Definition at line 17 of file TagProbeFitTreeAnalyzer.cc.

{};

Member Function Documentation

void TagProbeFitTreeAnalyzer::analyze ( const edm::Event &  event, const edm::EventSetup &  eventSetup  )

inlineoverride

Definition at line 18 of file TagProbeFitTreeAnalyzer.cc.

{};

void TagProbeFitTreeAnalyzer::calculateEfficiency	(	string	name,
		const edm::ParameterSet &	pset
	)

Definition at line 129 of file TagProbeFitTreeAnalyzer.cc.

References RecoTauDiscriminantConfiguration::binning, TagProbeFitter::calculateEfficiency(), gather_cfg::cout, DEFINE_FWK_MODULE, edm::ParameterSet::existsAs(), cmsRelvalreport::exit, fitter, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNamesForType(), mps_fire::i, genParticles_cff::map, gen::n, and JetChargeProducer_cfi::var.

Referenced by TagProbeFitTreeAnalyzer().

                                                                                         {
  vector<string> effCatState = pset.getParameter<vector<string> >("EfficiencyCategoryAndState");
  if(effCatState.empty() ||  (effCatState.size() % 2 == 1)){
    cout<<"EfficiencyCategoryAndState must be a even-sized list of category names and states of that category (cat1, state1, cat2, state2, ...)."<<endl;
    exit(1);
  }

  vector<string> unbinnedVariables;
  if(pset.existsAs<vector<string> >("UnbinnedVariables")){
    unbinnedVariables = pset.getParameter<vector<string> >("UnbinnedVariables");
  }

  const ParameterSet binVars = pset.getParameter<ParameterSet>("BinnedVariables");
  map<string, vector<double> >binnedVariables;
  vector<string> variableNames = binVars.getParameterNamesForType<vector<double> >();
  for (vector<string>::const_iterator var = variableNames.begin(); var != variableNames.end(); var++) {
    vector<double> binning = binVars.getParameter<vector<double> >(*var);
    binnedVariables[*var] = binning;
  }
  map<string, vector<string> >mappedCategories;
  vector<string> categoryNames = binVars.getParameterNamesForType<vector<string> >();
  for (vector<string>::const_iterator var = categoryNames.begin(); var != categoryNames.end(); var++) {
    vector<string> map = binVars.getParameter<vector<string> >(*var);
    mappedCategories[*var] = map;
  }

  vector<string> binToPDFmap;
  if(pset.existsAs<vector<string> >("BinToPDFmap")){
    binToPDFmap = pset.getParameter<vector<string> >("BinToPDFmap");
  }
  if((!binToPDFmap.empty()) && (binToPDFmap.size()%2 == 0)){
    cout<<"BinToPDFmap must have odd size, first string is the default, followed by binRegExp - PDFname pairs!"<<endl;
    exit(2);
  }

  vector<string> effCats, effStates;
  for (size_t i = 0, n = effCatState.size()/2; i < n; ++i) {
    effCats.push_back(effCatState[2*i]);
    effStates.push_back(effCatState[2*i+1]);
  }

  fitter.calculateEfficiency(name, effCats, effStates, unbinnedVariables, binnedVariables, mappedCategories, binToPDFmap);
}

void TagProbeFitTreeAnalyzer::endRun ( const edm::Run &  run, const edm::EventSetup &  setup  )

inlineoverride

Definition at line 19 of file TagProbeFitTreeAnalyzer.cc.

References dataset::name.

{};

Member Data Documentation

TagProbeFitter TagProbeFitTreeAnalyzer::fitter

private

Definition at line 22 of file TagProbeFitTreeAnalyzer.cc.

Referenced by calculateEfficiency(), and TagProbeFitTreeAnalyzer().

unsigned int TagProbeFitTreeAnalyzer::split_mode

private

Definition at line 23 of file TagProbeFitTreeAnalyzer.cc.

Referenced by TagProbeFitTreeAnalyzer().