#include <anautil.h>

Public Member Functions
void	addCut (TString n, std::function< bool()> pass, std::function< float()> weight)

void	addCutSyst (TString syst, std::vector< TString > pattern, std::vector< TString > vetopattern=std::vector< TString >())

void	addCutToLastActiveCut (TString n, std::function< bool()> pass, std::function< float()> weight)

void	addCutToSkipCutflowList (TString n)

void	addToCutTreeMap (TString n)

void	addWgtSyst (TString syst, std::function< float()>)

void	book2DHistogram (TString, std::pair< std::pair< TString, TString >, std::tuple< unsigned, float, float, unsigned, float, float, std::function< float()>, std::function< float()>>>, TString="")

void	book2DVecHistogram (TString, std::pair< std::pair< TString, TString >, std::tuple< unsigned, float, float, unsigned, float, float, std::function< std::vector< float >()>, std::function< std::vector< float >()>, std::function< std::vector< float >()>>>, TString="")

void	book2DVecHistogram (TString, std::pair< std::pair< TString, TString >, std::tuple< std::vector< float >, unsigned, float, float, std::function< std::vector< float >()>, std::function< std::vector< float >()>, std::function< std::vector< float >()>>>, TString="")

void	bookCutflowHistograms ()

void	bookCutflowHistograms_v1 ()

void	bookCutflowHistograms_v2 ()

void	bookCutflows ()

void	bookCutflowsForRegions (std::vector< TString > regions)

void	bookCutflowTree ()

void	bookEventLists ()

void	bookHistogram (TString, std::pair< TString, std::tuple< unsigned, float, float, std::function< float()>>>, TString="")

void	bookHistogram (TString, std::pair< TString, std::tuple< std::vector< float >, std::function< float()>>>, TString="")

void	bookHistograms (Histograms &histograms)

void	bookHistograms (Histograms &histograms, std::vector< TString > cutlist)

void	bookHistogramsForCut (Histograms &histograms, TString)

void	bookHistogramsForCutAndAbove (Histograms &histograms, TString)

void	bookHistogramsForCutAndBelow (Histograms &histograms, TString)

void	bookHistogramsForEndCuts (Histograms &histograms)

void	bookVecHistogram (TString, std::pair< TString, std::tuple< unsigned, float, float, std::function< std::vector< float >()>, std::function< std::vector< float >()>>>, TString="")

void	bookVecHistogram (TString, std::pair< TString, std::tuple< std::vector< float >, std::function< std::vector< float >()>, std::function< std::vector< float >()>>>, TString="")

void	copyAndEditCuts (TString, std::map< TString, TString >)

void	createWgtSystBranches ()

	Cutflow ()

	Cutflow (TFile *o)

void	fill ()

void	fillCutflow (std::vector< TString > &cutlist, TH1D h, TH1D hraw, float wgtsyst=1)

void	fillCutflow_v2 (std::vector< CutTree > &cutlist, TH1D h, TH1D *hraw, float wgtsyst=1)

void	fillCutflows (TString syst="", bool iswgtsyst=true)

void	fillCutflows_v1 (TString syst="", bool iswgtsyst=true)

void	fillCutflows_v2 (TString syst="", bool iswgtsyst=true)

void	fillHistograms (TString syst="", bool iswgtsyst=true)

void	filterCuts (std::vector< TString > ns)

CutTree &	getCut (TString n)

void	printCuts ()

void	printSetFunctionError (TString msg)

void	removeCut (TString n)

void	saveCutflows ()

void	saveHistograms ()

void	saveOutput ()

void	saveTTreeX ()

void	setCut (TString cutname, std::function< bool()> pass, std::function< float()> weight)

void	setCutLists (std::vector< TString > regions)

void	setCutSyst (TString cutname, TString syst, std::function< bool()> pass, std::function< float()> weight)

void	setEventID (int, int, unsigned long long)

void	setHistsAxesExtendable ()

void	setLastActiveCut (TString n)

void	setSaveTTreeX (bool=true)

void	setSkipSystematicHistograms (bool=true)

void	setTFile (TFile *o)

void	setVariable (TString varname, float)

void	setWgtSyst (TString syst, float weight)

	~Cutflow ()

Public Attributes
std::map< std::tuple< std::string, std::string, std::string, std::string >, TH2F * >	booked_2dhistograms

std::vector< std::tuple< std::string, std::string, std::string, std::string > >	booked_2dhistograms_nominal_keys

std::map< std::tuple< std::string, std::string, std::string >, TH1D * >	booked_histograms

std::vector< std::tuple< std::string, std::string, std::string > >	booked_histograms_nominal_keys

bool	cutflow_booked

std::map< TString, TH1D * >	cutflow_histograms

std::vector< std::tuple< TH1D , std::vector< int >, std::vector< float * >, std::function< float()> > >	cutflow_histograms_v2

std::vector< TString >	cutflow_nofill_cut_list

std::map< TString, std::vector< TString > >	cutlists

std::vector< TString >	cutsysts

CutTree	cuttree

std::map< TString, std::vector< CutTree * > >	cuttreelists

std::map< std::string, CutTree * >	cuttreemap

bool	dosavettreex

bool	doskipsysthist

bool	iseventlistbooked

CutTree *	last_active_cut

TFile *	ofile

std::map< TString, TH1D * >	rawcutflow_histograms

std::vector< std::tuple< TH1D , std::vector< int > > >	rawcutflow_histograms_v2

int	seterrorcount

std::vector< TString >	systs

std::map< TString, std::function< float()> >	systs_funcs

TTree *	t

TTreeX *	tx

Detailed Description

Definition at line 113 of file anautil.h.

Constructor & Destructor Documentation

◆ Cutflow() [1/2]

RooUtil::Cutflow::Cutflow ( )

Definition at line 7 of file anautil.cc.

References cuttree, and cuttreemap.

     : cuttree("Root"),
       last_active_cut(0),
       ofile(0),
       t(0),
       tx(0),
       iseventlistbooked(false),
       seterrorcount(0),
       doskipsysthist(0),
       dosavettreex(0),
       cutflow_booked(false) {
   cuttreemap["Root"] = &cuttree;
 }

◆ Cutflow() [2/2]

RooUtil::Cutflow::Cutflow ( TFile * o )

Definition at line 22 of file anautil.cc.

References cuttree, and cuttreemap.

     : cuttree("Root"),
       last_active_cut(0),
       ofile(o),
       t(0),
       tx(0),
       iseventlistbooked(false),
       seterrorcount(0),
       doskipsysthist(0),
       dosavettreex(0),
       cutflow_booked(false) {
   cuttreemap["Root"] = &cuttree;
 }

◆ ~Cutflow()

RooUtil::Cutflow::~Cutflow ( )

Definition at line 37 of file anautil.cc.

References isotrackApplyRegressor::k, and findQualityFiles::v.

                          {
   if (booked_histograms.size() > 0) {
     for (auto& [k, v] : booked_histograms) {
       delete v;
     }
   }
   if (booked_2dhistograms.size() > 0) {
     for (auto& [k, v] : booked_2dhistograms) {
       delete v;
     }
   }
 }

Member Function Documentation

◆ addCut()

void RooUtil::Cutflow::addCut	(	TString	n,
		std::function< bool()>	pass,
		std::function< float()>	weight
	)

Definition at line 73 of file anautil.cc.

References DMR_cfg::cut, and create_idmaps::n.

                                                                                            {
   cuttree.addCut(n);
   addToCutTreeMap(n);
   setLastActiveCut(n);
   setCut(n, cut, weight);
 }

◆ addCutSyst()

void RooUtil::Cutflow::addCutSyst	(	TString	syst,
		std::vector< TString >	pattern,
		std::vector< TString >	vetopattern = `std::vector<TString>()`
	)

Definition at line 412 of file anautil.cc.

References topSingleLeptonDQM_PU_cfi::pattern.

                                                                                                         {
   cutsysts.push_back(syst);
   cuttree.addSyst(syst, pattern, vetopattern);
 }

◆ addCutToLastActiveCut()

void RooUtil::Cutflow::addCutToLastActiveCut	(	TString	n,
		std::function< bool()>	pass,
		std::function< float()>	weight
	)

Definition at line 81 of file anautil.cc.

References DMR_cfg::cut, and create_idmaps::n.

                                                                                                           {
   last_active_cut->addCut(n);
   addToCutTreeMap(n);
   setLastActiveCut(n);
   setCut(n, cut, weight);
 }

◆ addCutToSkipCutflowList()

void RooUtil::Cutflow::addCutToSkipCutflowList ( TString n )

Definition at line 159 of file anautil.cc.

References spr::find(), and create_idmaps::n.

                                                       {
   if (std::find(cutflow_nofill_cut_list.begin(), cutflow_nofill_cut_list.end(), n) == cutflow_nofill_cut_list.end())
     cutflow_nofill_cut_list.push_back(n);
 }

◆ addToCutTreeMap()

void RooUtil::Cutflow::addToCutTreeMap ( TString n )

Definition at line 51 of file anautil.cc.

References RooUtil::error(), and create_idmaps::n.

                                               {
   if (cuttreemap.find(n.Data()) == cuttreemap.end())
     cuttreemap[n.Data()] = cuttree.getCutPointer(n);
   else
     error(TString::Format("Cut %s already exists! no duplicate cut names allowed!", n.Data()));
 }

◆ addWgtSyst()

void RooUtil::Cutflow::addWgtSyst	(	TString	syst,
		std::function< float()>	weight
	)

Definition at line 419 of file anautil.cc.

References mps_merge::weight.

                                                                          {
   systs.push_back(syst);
   systs_funcs[syst] = weight;
 }

◆ book2DHistogram()

void RooUtil::Cutflow::book2DHistogram	(	TString	cut,
		std::pair< std::pair< TString, TString >, std::tuple< unsigned, float, float, unsigned, float, float, std::function< float()>, std::function< float()>>>	key,
		TString	syst = `""`
	)

Definition at line 786 of file anautil.cc.

References DMR_cfg::cut, submitPVResolutionJobs::key, WZElectronSkims53X_cff::max, and SiStripPI::min.

                   {
   TString varname = key.first.first;
   TString varnamey = key.first.second;
   unsigned int nbin = std::get<0>(key.second);
   float min = std::get<1>(key.second);
   float max = std::get<2>(key.second);
   unsigned int nbiny = std::get<3>(key.second);
   float miny = std::get<4>(key.second);
   float maxy = std::get<5>(key.second);
   std::function<float()> varxdef = std::get<6>(key.second);
   std::function<float()> varydef = std::get<7>(key.second);
   TString histname = cut + syst + "__" + varname + "_v_" + varnamey;
   if (booked_2dhistograms.find(std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())) ==
       booked_2dhistograms.end()) {
     booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())] =
         new TH2F(histname, "", nbin, min, max, nbiny, miny, maxy);
     booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())]->SetDirectory(0);
     booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())]->Sumw2();
     if (syst.IsNull()) {
       booked_2dhistograms_nominal_keys.push_back(
           std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data()));
     }
     cuttreemap[cut.Data()]->addHist2D(
         booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())],
         varxdef,
         varydef,
         syst);
   }
 }

◆ book2DVecHistogram() [1/2]

void RooUtil::Cutflow::book2DVecHistogram	(	TString	cut,
		std::pair< std::pair< TString, TString >, std::tuple< unsigned, float, float, unsigned, float, float, std::function< std::vector< float >()>, std::function< std::vector< float >()>, std::function< std::vector< float >()>>>	key,
		TString	syst = `""`
	)

Definition at line 822 of file anautil.cc.

References DMR_cfg::cut, submitPVResolutionJobs::key, WZElectronSkims53X_cff::max, and SiStripPI::min.

                                                         {
   TString varname = key.first.first;
   TString varnamey = key.first.second;
   unsigned int nbin = std::get<0>(key.second);
   float min = std::get<1>(key.second);
   float max = std::get<2>(key.second);
   unsigned int nbiny = std::get<3>(key.second);
   float miny = std::get<4>(key.second);
   float maxy = std::get<5>(key.second);
   std::function<std::vector<float>()> varxdef = std::get<6>(key.second);
   std::function<std::vector<float>()> varydef = std::get<7>(key.second);
   std::function<std::vector<float>()> elemwgt = std::get<8>(key.second);
   TString histname = cut + syst + "__" + varname + "_v_" + varnamey;
   if (booked_2dhistograms.find(std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())) ==
       booked_2dhistograms.end()) {
     booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())] =
         new TH2F(histname, "", nbin, min, max, nbiny, miny, maxy);
     booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())]->SetDirectory(0);
     booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())]->Sumw2();
     if (syst.IsNull()) {
       booked_2dhistograms_nominal_keys.push_back(
           std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data()));
     }
     cuttreemap[cut.Data()]->addHist2DVec(
         booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())],
         varxdef,
         varydef,
         elemwgt,
         syst);
   }
 }

◆ book2DVecHistogram() [2/2]

void RooUtil::Cutflow::book2DVecHistogram	(	TString	cut,
		std::pair< std::pair< TString, TString >, std::tuple< std::vector< float >, unsigned, float, float, std::function< std::vector< float >()>, std::function< std::vector< float >()>, std::function< std::vector< float >()>>>	key,
		TString	syst = `""`
	)

Definition at line 865 of file anautil.cc.

References DMR_cfg::cut, mps_fire::i, and submitPVResolutionJobs::key.

                                                         {
   TString varname = key.first.first;
   TString varnamey = key.first.second;
   std::vector<float> xboundaries = std::get<0>(key.second);
   unsigned int nbiny = std::get<1>(key.second);
   float miny = std::get<2>(key.second);
   float maxy = std::get<3>(key.second);
   std::function<std::vector<float>()> varxdef = std::get<4>(key.second);
   std::function<std::vector<float>()> varydef = std::get<5>(key.second);
   std::function<std::vector<float>()> elemwgt = std::get<6>(key.second);
   TString histname = cut + syst + "__" + varname + "_v_" + varnamey;
   if (booked_2dhistograms.find(std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())) ==
       booked_2dhistograms.end()) {
     Double_t xbounds[xboundaries.size()];
     for (unsigned int i = 0; i < xboundaries.size(); ++i)
       xbounds[i] = xboundaries[i];
     booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())] =
         new TH2F(histname, "", xboundaries.size() - 1, xbounds, nbiny, miny, maxy);
     booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())]->SetDirectory(0);
     booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())]->Sumw2();
     if (syst.IsNull()) {
       booked_2dhistograms_nominal_keys.push_back(
           std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data()));
     }
     cuttreemap[cut.Data()]->addHist2DVec(
         booked_2dhistograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data(), varnamey.Data())],
         varxdef,
         varydef,
         elemwgt,
         syst);
   }
 }

◆ bookCutflowHistograms()

void RooUtil::Cutflow::bookCutflowHistograms ( )

Definition at line 181 of file anautil.cc.

References RooUtil::error().

                                            {
   if (cutflow_booked) {
     error("bookCutflowHistograms():: cutflows already booked! yet you called again to book it! Check your user code!");
   }
   cutflow_booked = true;
   bookCutflowHistograms_v1();
 }

◆ bookCutflowHistograms_v1()

void RooUtil::Cutflow::bookCutflowHistograms_v1 ( )

Definition at line 259 of file anautil.cc.

References RooUtil::CutflowUtil::createCutflowHistograms(), mps_check::msg, and RooUtil::print().

                                               {
   ofile->cd();
 
   std::tie(cutflow_histograms, rawcutflow_histograms) = RooUtil::CutflowUtil::createCutflowHistograms(cutlists);
 
   for (auto& syst : systs) {
     std::map<CUTFLOWMAPSTRING, THist*> cutflow_histograms_tmp;
     std::map<CUTFLOWMAPSTRING, THist*> rawcutflow_histograms_tmp;
     std::tie(cutflow_histograms_tmp, rawcutflow_histograms_tmp) =
         RooUtil::CutflowUtil::createCutflowHistograms(cutlists, syst);
     cutflow_histograms.insert(cutflow_histograms_tmp.begin(), cutflow_histograms_tmp.end());
     rawcutflow_histograms.insert(rawcutflow_histograms_tmp.begin(), rawcutflow_histograms_tmp.end());
   }
 
   for (auto& cutsyst : cutsysts) {
     std::map<CUTFLOWMAPSTRING, THist*> cutflow_histograms_tmp;
     std::map<CUTFLOWMAPSTRING, THist*> rawcutflow_histograms_tmp;
     std::tie(cutflow_histograms_tmp, rawcutflow_histograms_tmp) =
         RooUtil::CutflowUtil::createCutflowHistograms(cutlists, cutsyst);
     cutflow_histograms.insert(cutflow_histograms_tmp.begin(), cutflow_histograms_tmp.end());
     rawcutflow_histograms.insert(rawcutflow_histograms_tmp.begin(), rawcutflow_histograms_tmp.end());
   }
 
   for (auto& cutflow_histogram : cutflow_histograms) {
     TString msg = "Booked cutflow histogram for cut = " + cutflow_histogram.first;
     print(msg);
   }
   for (auto& rawcutflow_histogram : rawcutflow_histograms) {
     TString msg = "Booked rawcutflow histogram for cut = " + rawcutflow_histogram.first;
     print(msg);
   }
 }

◆ bookCutflowHistograms_v2()

void RooUtil::Cutflow::bookCutflowHistograms_v2 ( )

Definition at line 190 of file anautil.cc.

References hgcalPlots::ct, h, mps_fire::i, THist, UNITY(), and hltDeepSecondaryVertexTagInfosPFPuppi_cfi::weights.

                                               {
   ofile->cd();
 
   // Nominal
   for (auto& cuttreelist : cuttreelists) {
     THist* h = new THist(cuttreelist.first + "_cutflow", "", cuttreelist.second.size(), 0, cuttreelist.second.size());
     THist* h_raw =
         new THist(cuttreelist.first + "_rawcutflow", "", cuttreelist.second.size(), 0, cuttreelist.second.size());
     h->Sumw2();
     h_raw->Sumw2();
     h->SetDirectory(0);
     h_raw->SetDirectory(0);
     std::vector<int*> passes;
     std::vector<float*> weights;
     for (auto& ct : cuttreelist.second) {
       passes.push_back(&(ct->pass));
       weights.push_back(&(ct->weight));
     }
     cutflow_histograms_v2.push_back(std::make_tuple(h, passes, weights, UNITY));
     rawcutflow_histograms_v2.push_back(std::make_tuple(h_raw, passes));
   }
 
   // Wgt varying systematic cutflow
   for (auto& syst : systs) {
     for (auto& cuttreelist : cuttreelists) {
       THist* h =
           new THist(cuttreelist.first + syst + "_cutflow", "", cuttreelist.second.size(), 0, cuttreelist.second.size());
       THist* h_raw = new THist(
           cuttreelist.first + syst + "_rawcutflow", "", cuttreelist.second.size(), 0, cuttreelist.second.size());
       h->Sumw2();
       h_raw->Sumw2();
       h->SetDirectory(0);
       h_raw->SetDirectory(0);
       std::vector<int*> passes;
       std::vector<float*> weights;
       for (auto& ct : cuttreelist.second) {
         passes.push_back(&(ct->pass));
         weights.push_back(&(ct->weight));
       }
       cutflow_histograms_v2.push_back(std::make_tuple(h, passes, weights, systs_funcs[syst]));
       rawcutflow_histograms_v2.push_back(std::make_tuple(h_raw, passes));
     }
   }
 
   // Nominal
   for (unsigned int i = 0; i < cutsysts.size(); ++i) {
     TString cutsyst = cutsysts[i];
     for (auto& cuttreelist : cuttreelists) {
       THist* h = new THist(
           cuttreelist.first + cutsyst + "_cutflow", "", cuttreelist.second.size(), 0, cuttreelist.second.size());
       THist* h_raw = new THist(
           cuttreelist.first + cutsyst + "_rawcutflow", "", cuttreelist.second.size(), 0, cuttreelist.second.size());
       h->Sumw2();
       h_raw->Sumw2();
       h->SetDirectory(0);
       h_raw->SetDirectory(0);
       std::vector<int*> passes;
       std::vector<float*> weights;
       for (auto& ct : cuttreelist.second) {
         passes.push_back(&(ct->systpasses[i]));
         weights.push_back(&(ct->systweights[i]));
       }
       cutflow_histograms_v2.push_back(std::make_tuple(h, passes, weights, UNITY));
       rawcutflow_histograms_v2.push_back(std::make_tuple(h_raw, passes));
     }
   }
 }

◆ bookCutflows()

void RooUtil::Cutflow::bookCutflows ( )

Definition at line 300 of file anautil.cc.

                                   {
   std::vector<TString> regions = cuttree.getEndCuts();
   setCutLists(regions);
   bookCutflowTree();
   bookCutflowHistograms();
 }

◆ bookCutflowsForRegions()

void RooUtil::Cutflow::bookCutflowsForRegions ( std::vector< TString > regions )

Definition at line 293 of file anautil.cc.

                                                                       {
   setCutLists(regions);
   bookCutflowTree();
   bookCutflowHistograms();
 }

◆ bookCutflowTree()

void RooUtil::Cutflow::bookCutflowTree ( )

Definition at line 165 of file anautil.cc.

References RooUtil::CutflowUtil::createCutflowBranches(), and submitPVValidationJobs::t.

                                      {
   if (!t) {
     ofile->cd();
     t = new TTree("cut_tree", "cut_tree");
     t->SetDirectory(0);
   }
   if (!tx) {
     ofile->cd();
     tx = new TTreeX(t);
   }
   RooUtil::CutflowUtil::createCutflowBranches(cutlists, *tx);
   createWgtSystBranches();
   //    t->Print();
 }

◆ bookEventLists()

void RooUtil::Cutflow::bookEventLists ( )

Definition at line 459 of file anautil.cc.

References RooUtil::error().

                                     {
   if (!tx)
     error("bookEventLists():: No TTreeX has been set. Forgot to call bookCutflows()?");
   if (!tx->hasBranch<int>("run"))
     tx->createBranch<int>("run");
   if (!tx->hasBranch<int>("lumi"))
     tx->createBranch<int>("lumi");
   if (!tx->hasBranch<unsigned long long>("evt"))
     tx->createBranch<unsigned long long>("evt");
   iseventlistbooked = true;
 }

◆ bookHistogram() [1/2]

void RooUtil::Cutflow::bookHistogram	(	TString	cut,
		std::pair< TString, std::tuple< unsigned, float, float, std::function< float()>>>	key,
		TString	syst = `""`
	)

Definition at line 683 of file anautil.cc.

References DMR_cfg::cut, submitPVResolutionJobs::key, WZElectronSkims53X_cff::max, SiStripPI::min, and THist.

                                                    {
   TString varname = key.first;
   unsigned int nbin = std::get<0>(key.second);
   float min = std::get<1>(key.second);
   float max = std::get<2>(key.second);
   std::function<float()> vardef = std::get<3>(key.second);
   TString histname = cut + syst + "__" + varname;
   if (booked_histograms.find(std::make_tuple(cut.Data(), syst.Data(), varname.Data())) == booked_histograms.end()) {
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())] =
         new THist(histname, "", nbin, min, max);
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())]->SetDirectory(0);
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())]->Sumw2();
     if (syst.IsNull()) {
       booked_histograms_nominal_keys.push_back(std::make_tuple(cut.Data(), syst.Data(), varname.Data()));
     }
     cuttreemap[cut.Data()]->addHist1D(
         booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())], vardef, syst);
   }
 }

◆ bookHistogram() [2/2]

void RooUtil::Cutflow::bookHistogram	(	TString	cut,
		std::pair< TString, std::tuple< std::vector< float >, std::function< float()>>>	key,
		TString	syst = `""`
	)

Definition at line 734 of file anautil.cc.

References DMR_cfg::cut, mps_fire::i, submitPVResolutionJobs::key, and THist.

                                                    {
   TString varname = key.first;
   std::vector<float> boundaries = std::get<0>(key.second);
   std::function<float()> vardef = std::get<1>(key.second);
   TString histname = cut + syst + "__" + varname;
   if (booked_histograms.find(std::make_tuple(cut.Data(), syst.Data(), varname.Data())) == booked_histograms.end()) {
     Float_t bounds[boundaries.size()];
     for (unsigned int i = 0; i < boundaries.size(); ++i)
       bounds[i] = boundaries[i];
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())] =
         new THist(histname, "", boundaries.size() - 1, bounds);
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())]->SetDirectory(0);
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())]->Sumw2();
     if (syst.IsNull()) {
       booked_histograms_nominal_keys.push_back(std::make_tuple(cut.Data(), syst.Data(), varname.Data()));
     }
     cuttreemap[cut.Data()]->addHist1D(
         booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())], vardef, syst);
   }
 }

◆ bookHistograms() [1/2]

void RooUtil::Cutflow::bookHistograms ( Histograms & histograms )

Definition at line 666 of file anautil.cc.

References bookHistograms(), SimL1EmulatorRepack_CalouGT_cff::cutlist, and nano_mu_digi_cff::region.

Referenced by main().

                                                           {
   std::vector<TString> regions = cuttree.getEndCuts();
   for (auto& region : regions) {
     std::vector<TString> cutlist = cuttree.getCutList(region);
     bookHistograms(histograms, cutlist);
   }
 }

◆ bookHistograms() [2/2]

void RooUtil::Cutflow::bookHistograms	(	Histograms &	histograms,
		std::vector< TString >	cutlist
	)

Definition at line 675 of file anautil.cc.

References DMR_cfg::cut, and SimL1EmulatorRepack_CalouGT_cff::cutlist.

                                                                                       {
   for (auto& cut : cutlist) {
     bookHistogramsForCut(histograms, cut);
   }
 }

◆ bookHistogramsForCut()

void RooUtil::Cutflow::bookHistogramsForCut	(	Histograms &	histograms,
		TString	cut
	)

Definition at line 998 of file anautil.cc.

References DMR_cfg::cut, and submitPVResolutionJobs::key.

                                                                              {
   for (auto& key : histograms.th1fs)
     bookHistogram(cut, key);
   for (auto& key : histograms.th1fs_varbin)
     bookHistogram(cut, key);
   for (auto& key : histograms.th1vecfs)
     bookVecHistogram(cut, key);
   for (auto& key : histograms.th1vecfs_varbin)
     bookVecHistogram(cut, key);
   for (auto& key : histograms.th2fs)
     book2DHistogram(cut, key);
   for (auto& key : histograms.th2vecfs)
     book2DVecHistogram(cut, key);
   for (auto& key : histograms.th2vecfs_xvarbin)
     book2DVecHistogram(cut, key);
   if (not doskipsysthist) {
     for (auto& syst : systs) {
       for (auto& key : histograms.th1fs)
         bookHistogram(cut, key, syst);
       for (auto& key : histograms.th1fs_varbin)
         bookHistogram(cut, key, syst);
       for (auto& key : histograms.th1vecfs)
         bookVecHistogram(cut, key, syst);
       for (auto& key : histograms.th1vecfs_varbin)
         bookVecHistogram(cut, key, syst);
       for (auto& key : histograms.th2fs)
         book2DHistogram(cut, key, syst);
       for (auto& key : histograms.th2vecfs)
         book2DVecHistogram(cut, key, syst);
       for (auto& key : histograms.th2vecfs_xvarbin)
         book2DVecHistogram(cut, key, syst);
     }
     for (auto& cutsyst : cutsysts) {
       for (auto& key : histograms.th1fs)
         bookHistogram(cut, key, cutsyst);
       for (auto& key : histograms.th1fs_varbin)
         bookHistogram(cut, key, cutsyst);
       for (auto& key : histograms.th1vecfs)
         bookVecHistogram(cut, key, cutsyst);
       for (auto& key : histograms.th1vecfs_varbin)
         bookVecHistogram(cut, key, cutsyst);
       for (auto& key : histograms.th2fs)
         book2DHistogram(cut, key, cutsyst);
       for (auto& key : histograms.th2vecfs)
         book2DVecHistogram(cut, key, cutsyst);
       for (auto& key : histograms.th2vecfs_xvarbin)
         book2DVecHistogram(cut, key, cutsyst);
     }
   }
 }

◆ bookHistogramsForCutAndAbove()

void RooUtil::Cutflow::bookHistogramsForCutAndAbove	(	Histograms &	histograms,
		TString	cut
	)

Definition at line 1058 of file anautil.cc.

References RooUtil::error().

                                                                                      {
   error("bookHistogramsForCutAndAbove not yet implemented");
 }

◆ bookHistogramsForCutAndBelow()

void RooUtil::Cutflow::bookHistogramsForCutAndBelow	(	Histograms &	histograms,
		TString	cut
	)

Definition at line 1050 of file anautil.cc.

References DummyCfis::c, DMR_cfg::cut, and SimL1EmulatorRepack_CalouGT_cff::cutlist.

                                                                                      {
   std::vector<TString> cutlist = cuttree.getCutListBelow(cut);
   for (auto& c : cutlist) {
     bookHistogramsForCut(histograms, c);
   }
 }

◆ bookHistogramsForEndCuts()

void RooUtil::Cutflow::bookHistogramsForEndCuts ( Histograms & histograms )

Definition at line 1063 of file anautil.cc.

References nano_mu_digi_cff::region.

                                                                     {
   std::vector<TString> regions = cuttree.getEndCuts();
   for (auto& region : regions) {
     bookHistogramsForCut(histograms, region);
   }
 }

◆ bookVecHistogram() [1/2]

void RooUtil::Cutflow::bookVecHistogram	(	TString	cut,
		std::pair< TString, std::tuple< unsigned, float, float, std::function< std::vector< float >()>, std::function< std::vector< float >()>>>	key,
		TString	syst = `""`
	)

Definition at line 706 of file anautil.cc.

References DMR_cfg::cut, submitPVResolutionJobs::key, WZElectronSkims53X_cff::max, SiStripPI::min, and THist.

                   {
   TString varname = key.first;
   unsigned int nbin = std::get<0>(key.second);
   float min = std::get<1>(key.second);
   float max = std::get<2>(key.second);
   std::function<std::vector<float>()> vardef = std::get<3>(key.second);
   std::function<std::vector<float>()> wgtdef = std::get<4>(key.second);
   TString histname = cut + syst + "__" + varname;
   if (booked_histograms.find(std::make_tuple(cut.Data(), syst.Data(), varname.Data())) == booked_histograms.end()) {
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())] =
         new THist(histname, "", nbin, min, max);
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())]->SetDirectory(0);
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())]->Sumw2();
     if (syst.IsNull()) {
       booked_histograms_nominal_keys.push_back(std::make_tuple(cut.Data(), syst.Data(), varname.Data()));
     }
     cuttreemap[cut.Data()]->addHist1DVec(
         booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())], vardef, wgtdef, syst);
   }
 }

◆ bookVecHistogram() [2/2]

void RooUtil::Cutflow::bookVecHistogram	(	TString	cut,
		std::pair< TString, std::tuple< std::vector< float >, std::function< std::vector< float >()>, std::function< std::vector< float >()>>>	key,
		TString	syst = `""`
	)

Definition at line 758 of file anautil.cc.

References DMR_cfg::cut, mps_fire::i, submitPVResolutionJobs::key, and THist.

                   {
   TString varname = key.first;
   std::vector<float> boundaries = std::get<0>(key.second);
   std::function<std::vector<float>()> vardef = std::get<1>(key.second);
   std::function<std::vector<float>()> wgtdef = std::get<2>(key.second);
   TString histname = cut + syst + "__" + varname;
   if (booked_histograms.find(std::make_tuple(cut.Data(), syst.Data(), varname.Data())) == booked_histograms.end()) {
     Float_t bounds[boundaries.size()];
     for (unsigned int i = 0; i < boundaries.size(); ++i)
       bounds[i] = boundaries[i];
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())] =
         new THist(histname, "", boundaries.size() - 1, bounds);
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())]->SetDirectory(0);
     booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())]->Sumw2();
     if (syst.IsNull()) {
       booked_histograms_nominal_keys.push_back(std::make_tuple(cut.Data(), syst.Data(), varname.Data()));
     }
     cuttreemap[cut.Data()]->addHist1DVec(
         booked_histograms[std::make_tuple(cut.Data(), syst.Data(), varname.Data())], vardef, wgtdef, syst);
   }
 }

◆ copyAndEditCuts()

void RooUtil::Cutflow::copyAndEditCuts	(	TString	,
		std::map< TString, TString >
	)

◆ createWgtSystBranches()

void RooUtil::Cutflow::createWgtSystBranches ( )

Definition at line 429 of file anautil.cc.

                                            {
   for (auto& syst : systs) {
     if (!tx->hasBranch<float>(syst))
       tx->createBranch<float>(syst);
   }
 }

◆ fill()

void RooUtil::Cutflow::fill ( void )

Definition at line 472 of file anautil.cc.

References RPCNoise_cfi::fillHistograms, and mps_check::msg.

Referenced by main().

                           {
 #ifdef USE_TTREEX
   if (!tx) {
     TString msg = "No TTreeX object set, fill()";
     printSetFunctionError(msg);
     return;
   }
   tx->setBranch<bool>("Root", 1);          // Root is internally set
   tx->setBranch<float>("Root_weight", 1);  // Root is internally set
 #else
   cuttreemap["Root"]->pass_this_cut = 1;
   cuttreemap["Root"]->weight_this_cut = 1;
 #endif
 
   // Evaluate nominal selection cutflows (the non cut varying selections)
   cuttree.evaluate(*tx, "", iseventlistbooked);
 
   // Nominal cutflow
   fillCutflows();
 
   // Wgt systematic variations
   for (auto& syst : systs)
     fillCutflows(syst);
 
   // Fill nominal histograms
   fillHistograms();
 
   if (not doskipsysthist) {
     // Wgt systematic variations
     for (auto& syst : systs)
       fillHistograms(syst);
   }
 
   for (auto& cutsyst : cutsysts) {
     cuttree.evaluate(*tx, cutsyst, iseventlistbooked);
     fillCutflows(cutsyst, false);
     if (not doskipsysthist)
       fillHistograms(cutsyst, false);
   }
 
   if (tx) {
     if (dosavettreex)
       tx->fill();
 
     tx->clear();
   }
 }

◆ fillCutflow()

void RooUtil::Cutflow::fillCutflow	(	std::vector< TString > &	cutlist,
		TH1D *	h,
		TH1D *	hraw,
		float	wgtsyst = `1`
	)

Definition at line 565 of file anautil.cc.

References SimL1EmulatorRepack_CalouGT_cff::cutlist, h, and mps_fire::i.

                                                                                                   {
   for (unsigned int i = 0; i < cutlist.size(); ++i) {
     int& pass = cuttreemap[cutlist[i].Data()]->pass;
     if (pass) {
       float& weight = cuttreemap[cutlist[i].Data()]->weight;
       h->Fill(i, weight * wgtsyst);
       hraw->Fill(i, 1);
     } else {
       return;
     }
   }
 }

◆ fillCutflow_v2()

void RooUtil::Cutflow::fillCutflow_v2	(	std::vector< CutTree *> &	cutlist,
		TH1D *	h,
		TH1D *	hraw,
		float	wgtsyst = `1`
	)

Definition at line 579 of file anautil.cc.

References hgcalPlots::ct, h, and mps_fire::i.

                                                      {
   for (unsigned int i = 0; i < cuttreelist.size(); ++i) {
     RooUtil::CutTree* ct = cuttreelist[i];
     // if (std::find(cutflow_nofill_cut_list.begin(), cutflow_nofill_cut_list.end(), ct->name) != cutflow_nofill_cut_list.end())
     //     continue;
     int& pass = ct->pass;
     if (pass) {
       float& weight = ct->weight;
       h->Fill(i, weight * wgtsyst);
       hraw->Fill(i, 1);
     } else {
       return;
     }
   }
 }

◆ fillCutflows()

void RooUtil::Cutflow::fillCutflows	(	TString	syst = `""`,
		bool	iswgtsyst = `true`
	)

Definition at line 522 of file anautil.cc.

522 { fillCutflows_v2(syst, iswgtsyst); }

RooUtil::Cutflow::fillCutflows_v2

void fillCutflows_v2(TString syst="", bool iswgtsyst=true)

Definition: anautil.cc:538

◆ fillCutflows_v1()

void RooUtil::Cutflow::fillCutflows_v1	(	TString	syst = `""`,
		bool	iswgtsyst = `true`
	)

Definition at line 525 of file anautil.cc.

References SimL1EmulatorRepack_CalouGT_cff::cutlist, and RooUtil::CutflowUtil::fillCutflow().

                                                                  {
   for (auto& pair : cutlists) {
     const TString& region_name = pair.first;
     std::vector<TString>& cutlist = pair.second;
     float wgtsyst = (!syst.IsNull() and iswgtsyst) ? systs_funcs[syst]() : 1;
     fillCutflow(cutlist,
                 cutflow_histograms[(region_name + syst).Data()],
                 rawcutflow_histograms[(region_name + syst).Data()],
                 wgtsyst);
   }
 }

◆ fillCutflows_v2()

void RooUtil::Cutflow::fillCutflows_v2	(	TString	syst = `""`,
		bool	iswgtsyst = `true`
	)

Definition at line 538 of file anautil.cc.

                                                                  {
   for (auto& pair : cuttreelists) {
     const TString& region_name = pair.first;
     std::vector<RooUtil::CutTree*>& cuttreelist = pair.second;
     float wgtsyst = (!syst.IsNull() and iswgtsyst) ? systs_funcs[syst]() : 1;
     fillCutflow_v2(cuttreelist,
                    cutflow_histograms[(region_name + syst).Data()],
                    rawcutflow_histograms[(region_name + syst).Data()],
                    wgtsyst);
   }
 }

◆ fillHistograms()

void RooUtil::Cutflow::fillHistograms	(	TString	syst = `""`,
		bool	iswgtsyst = `true`
	)

Definition at line 599 of file anautil.cc.

                                                                 {
 #ifdef USE_CUTLAMBDA
   float wgtsyst = (!syst.IsNull() and iswgtsyst) ? systs_funcs[syst]() : 1;
   cuttree.fillHistograms(syst, wgtsyst);
 #else
   float wgtsyst = (!syst.IsNull() and iswgtsyst) ? tx->getBranch<float>(syst) : 1.;
   cuttree.fillHistograms(*tx, syst, wgtsyst);
 #endif
 
   //    for (auto& key1d : booked_histograms_nominal_keys)
   //    {
   //        TString cutname = std::get<0>(key1d);
   //        std::get<1>(key1d) = syst;
   //        TString varname = std::get<2>(key1d);
   //        bool& passed = cuttreemap[cutname.Data()]->pass;
   //        if (!passed)
   //            continue;
   //        float& weight = cuttreemap[cutname.Data()]->weight;
   //        float wgtsyst = (!syst.IsNull() and iswgtsyst)? tx->getBranch<float>(syst) : 1.;
   //        THist* h = booked_histograms[key1d];
   //        h->Fill(tx->getBranch<float>(varname), weight * wgtsyst);
   //    }
   //
   //    for (auto& key2d : booked_2dhistograms_nominal_keys)
   //    {
   //        TString cutname = std::get<0>(key2d);
   //        std::get<1>(key2d) = syst;
   //        TString varname = std::get<2>(key2d);
   //        TString varnamey = std::get<3>(key2d);
   //        bool& passed = cuttreemap[cutname.Data()]->pass;
   //        if (!passed)
   //            continue;
   //        float& weight = cuttreemap[cutname.Data()]->weight;
   //        float wgtsyst = (!syst.IsNull() and iswgtsyst)? tx->getBranch<float>(syst) : 1.;
   //        TH2F* h = booked_2dhistograms[key2d];
   //        h->Fill(tx->getBranch<float>(varname), tx->getBranch<float>(varnamey), weight * wgtsyst);
   //    }
 
   //    for (auto& pair : booked_histograms)
   //    {
   //        TString cutname = std::get<0>(pair.first);
   //        TString sysname = std::get<1>(pair.first);
   //        TString varname = std::get<2>(pair.first);
   //        THist* h = pair.second;
   //        bool& passed = cuttreemap[cutname.DATA()]->pass;
   //        float& weight = cuttreemap[cutname.DATA()]->weight;
   //        float wgtsyst = sysname.IsNull() ? 1. : tx->getBranch<float>(sysname);
   //        if (passed)
   //            h->Fill(tx->getBranch<float>(varname), weight * wgtsyst);
   //    }
   //
   //    for (auto& pair : booked_2dhistograms)
   //    {
   //        TString cutname = std::get<0>(pair.first);
   //        TString sysname = std::get<1>(pair.first);
   //        TString varname = std::get<2>(pair.first);
   //        TString varnamey = std::get<3>(pair.first);
   //        TH2F* h = pair.second;
   //        bool& passed = cuttreemap[cutname.DATA()]->pass;
   //        float& weight = cuttreemap[cutname.DATA()]->weight;
   //        float wgtsyst = sysname.IsNull() ? 1. : tx->getBranch<float>(sysname);
   //        if (passed)
   //            h->Fill(tx->getBranch<float>(varname), tx->getBranch<float>(varnamey), weight * wgtsyst);
   //    }
 }

◆ filterCuts()

void RooUtil::Cutflow::filterCuts ( std::vector< TString > ns )

Definition at line 114 of file anautil.cc.

References submitPVValidationJobs::child, DMR_cfg::cut, SimL1EmulatorRepack_CalouGT_cff::cutlist, spr::find(), mps_fire::i, and create_idmaps::n.

                                                      {
   //
   std::vector<TString> to_not_remove;
   for (auto& n : ns) {
     std::vector<TString> cutlist = cuttree.getCutList(n);
     for (auto& cut : cutlist) {
       to_not_remove.push_back(cut);
     }
   }
 
   for (auto& n : ns) {
     std::vector<TString> cutlist = cuttree.getCutList(n);
     for (unsigned int i = 0; i < cutlist.size() - 1; ++i) {
       RooUtil::CutTree* cut = cuttree.getCutPointer(cutlist[i]);
       std::vector<RooUtil::CutTree*> toremove;
       for (auto& child : cut->children) {
         if (not child->name.EqualTo(cutlist[i + 1])) {
           if (std::find(to_not_remove.begin(), to_not_remove.end(), child->name) == to_not_remove.end())
             toremove.push_back(child);
         }
       }
       for (auto& child : toremove) {
         cut->children.erase(std::find(cut->children.begin(), cut->children.end(), child));
         cuttreemap.erase(cuttreemap.find(child->name.Data()));
       }
     }
   }
 }

◆ getCut()

RooUtil::CutTree & RooUtil::Cutflow::getCut ( TString n )

Definition at line 65 of file anautil.cc.

References DummyCfis::c, and create_idmaps::n.

                                                 {
   RooUtil::CutTree& c = cuttree.getCut(n);
   setLastActiveCut(n);
   return c;
 }

◆ printCuts()

void RooUtil::Cutflow::printCuts ( )

Definition at line 62 of file anautil.cc.

62 { cuttree.printCuts(); }

RooUtil::CutTree::printCuts

void printCuts(int indent=0, std::vector< int > multichild=std::vector< int >())

Definition: cutflowutil.h:119

RooUtil::Cutflow::cuttree

CutTree cuttree

Definition: anautil.h:115

◆ printSetFunctionError()

void RooUtil::Cutflow::printSetFunctionError ( TString msg )

Definition at line 1077 of file anautil.cc.

References mps_check::msg, and RooUtil::print().

                                                       {
   if (seterrorcount < 100) {
     print(msg);
     seterrorcount++;
   } else if (seterrorcount == 100) {
     print(msg);
     print("Suppressing Cutflow::set\"Func\"() errors ... ");
     seterrorcount++;
   } else {
     return;
   }
 }

◆ removeCut()

void RooUtil::Cutflow::removeCut ( TString n )

Definition at line 107 of file anautil.cc.

References DummyCfis::c, spr::find(), and create_idmaps::n.

                                         {
   RooUtil::CutTree* c = cuttree.getCutPointer(n);
   c->parent->children.erase(std::find(c->parent->children.begin(), c->parent->children.end(), c));
   cuttreemap.erase(cuttreemap.find(n.Data()));
 }

◆ saveCutflows()

void RooUtil::Cutflow::saveCutflows ( )

Definition at line 318 of file anautil.cc.

References RooUtil::CutflowUtil::saveCutflowHistograms().

                                   {
   // Save cutflow histograms
   ofile->cd();
   RooUtil::CutflowUtil::saveCutflowHistograms(cutflow_histograms, rawcutflow_histograms);
 }

◆ saveHistograms()

void RooUtil::Cutflow::saveHistograms ( )

Definition at line 325 of file anautil.cc.

                                     {
   ofile->cd();
   for (auto& pair : booked_histograms)
     pair.second->Write();
   for (auto& pair : booked_2dhistograms)
     pair.second->Write();
 }

◆ saveOutput()

void RooUtil::Cutflow::saveOutput ( )

Definition at line 308 of file anautil.cc.

References corrVsCorr::filename, mps_check::msg, RooUtil::print(), and ntuplePlotting::saveHistograms().

Referenced by main(), and run_lst().

                                 {
   saveCutflows();
   saveHistograms();
   saveTTreeX();
   TString filename = ofile->GetName();
   TString msg = "Wrote output to " + filename;
   print(msg);
 }

◆ saveTTreeX()

void RooUtil::Cutflow::saveTTreeX ( )

Definition at line 334 of file anautil.cc.

                                 {
   if (dosavettreex) {
     ofile->cd();
     tx->save(ofile);
   }
 }

◆ setCut()

void RooUtil::Cutflow::setCut	(	TString	cutname,
		std::function< bool()>	pass,
		std::function< float()>	weight
	)

Definition at line 343 of file anautil.cc.

References mps_merge::weight.

                                                                                                   {
   cuttreemap[cutname.Data()]->pass_this_cut_func = pass;
   cuttreemap[cutname.Data()]->weight_this_cut_func = weight;
 }

◆ setCutLists()

void RooUtil::Cutflow::setCutLists ( std::vector< TString > regions )

Definition at line 144 of file anautil.cc.

References nano_mu_digi_cff::region.

                                                            {
   for (auto& region : regions) {
     cutlists[region] = cuttree.getCutList(region);
     //        std::cout << region << std::endl;
     //        for (auto& cutname : cutlists[region])
     //            std::cout << cutname << std::endl;
   }
   for (auto& region : regions) {
     for (auto& cutname : cutlists[region]) {
       cuttreelists[region].push_back(cuttreemap[cutname.Data()]);
     }
   }
 }

◆ setCutSyst()

void RooUtil::Cutflow::setCutSyst	(	TString	cutname,
		TString	syst,
		std::function< bool()>	pass,
		std::function< float()>	weight
	)

Definition at line 349 of file anautil.cc.

References RooUtil::error(), and mps_merge::weight.

                                                                {
   if (cuttreemap[cutname.Data()]->systs.find(syst) == cuttreemap[cutname.Data()]->systs.end()) {
     error(
         TString::Format("setCutSyst():: Did not find syst=%s from the cut=%s! Did you actually book this syst for the "
                         "cut properly using addCutSyst() ?",
                         syst.Data(),
                         cutname.Data()));
   }
   cuttreemap[cutname.Data()]->systs[syst]->pass_this_cut_func = pass;
   cuttreemap[cutname.Data()]->systs[syst]->weight_this_cut_func = weight;
 }

◆ setEventID()

void RooUtil::Cutflow::setEventID	(	int	run,
		int	lumi,
		unsigned long long	evt
	)

Definition at line 447 of file anautil.cc.

References BXlumiParameters_cfi::lumi, mps_check::msg, and writedatasetfile::run.

                                                                          {
   if (!tx) {
     TString msg = "No TTreeX object set, setEventID()";
     printSetFunctionError(msg);
     return;
   }
   tx->setBranch<int>("run", run, false, true);
   tx->setBranch<int>("lumi", lumi, false, true);
   tx->setBranch<unsigned long long>("evt", evt, false, true);
 }

◆ setHistsAxesExtendable()

void RooUtil::Cutflow::setHistsAxesExtendable ( )

Definition at line 1091 of file anautil.cc.

                                             {
   for (auto& pair : booked_histograms)
     pair.second->SetCanExtend(TH1::kAllAxes);
 }

◆ setLastActiveCut()

void RooUtil::Cutflow::setLastActiveCut ( TString n )

Definition at line 59 of file anautil.cc.

References create_idmaps::n.

59 { last_active_cut = cuttree.getCutPointer(n); }

create_idmaps.n

n

Definition: create_idmaps.py:21

RooUtil::Cutflow::last_active_cut

CutTree * last_active_cut

Definition: anautil.h:116

RooUtil::Cutflow::cuttree

CutTree cuttree

Definition: anautil.h:115

RooUtil::CutTree::getCutPointer

CutTree * getCutPointer(TString n)

Definition: cutflowutil.h:248

◆ setSaveTTreeX()

void RooUtil::Cutflow::setSaveTTreeX ( bool v = true )

Definition at line 1074 of file anautil.cc.

References findQualityFiles::v.

1074 { dosavettreex = v; }

findQualityFiles.v

v

Definition: findQualityFiles.py:179

RooUtil::Cutflow::dosavettreex

bool dosavettreex

Definition: anautil.h:143

◆ setSkipSystematicHistograms()

void RooUtil::Cutflow::setSkipSystematicHistograms ( bool v = true )

Definition at line 1071 of file anautil.cc.

References findQualityFiles::v.

1071 { doskipsysthist = v; }

findQualityFiles.v

v

Definition: findQualityFiles.py:179

RooUtil::Cutflow::doskipsysthist

bool doskipsysthist

Definition: anautil.h:142

◆ setTFile()

void RooUtil::Cutflow::setTFile ( TFile * o )

inline

Definition at line 148 of file anautil.h.

References EcalTangentSkim_cfg::o, and ofile.

Referenced by initializeInputsAndOutputs(), and main().

148 { ofile = o; }

EcalTangentSkim_cfg.o

o

Definition: EcalTangentSkim_cfg.py:42

RooUtil::Cutflow::ofile

TFile * ofile

Definition: anautil.h:132

◆ setVariable()

void RooUtil::Cutflow::setVariable	(	TString	varname,
		float	val
	)

Definition at line 437 of file anautil.cc.

References mps_check::msg, and heppy_batch::val.

                                                            {
   if (!tx) {
     TString msg = "No TTreeX object set, setVariable() for " + varname;
     printSetFunctionError(msg);
     return;
   }
   tx->setBranch<float>(varname, val, false, true);
 }

◆ setWgtSyst()

void RooUtil::Cutflow::setWgtSyst	(	TString	syst,
		float	weight
	)

Definition at line 402 of file anautil.cc.

References mps_check::msg, and mps_merge::weight.

                                                           {
   if (!tx) {
     TString msg = "No TTreeX object set, setWgtSyst() for " + syst;
     printSetFunctionError(msg);
     return;
   }
   tx->setBranch<float>(syst, weight, false, true);
 }

Member Data Documentation

◆ booked_2dhistograms

std::map<std::tuple< std::string , std::string , std::string , std::string >, TH2F*> RooUtil::Cutflow::booked_2dhistograms

Definition at line 123 of file anautil.h.

◆ booked_2dhistograms_nominal_keys

std::vector<std::tuple< std::string , std::string , std::string , std::string > > RooUtil::Cutflow::booked_2dhistograms_nominal_keys

Definition at line 127 of file anautil.h.

◆ booked_histograms

std::map<std::tuple< std::string , std::string , std::string >, TH1D *> RooUtil::Cutflow::booked_histograms

Definition at line 121 of file anautil.h.

◆ booked_histograms_nominal_keys

std::vector<std::tuple< std::string , std::string , std::string > > RooUtil::Cutflow::booked_histograms_nominal_keys

Definition at line 125 of file anautil.h.

◆ cutflow_booked

bool RooUtil::Cutflow::cutflow_booked

Definition at line 144 of file anautil.h.

◆ cutflow_histograms

std::map< TString , TH1D *> RooUtil::Cutflow::cutflow_histograms

Definition at line 118 of file anautil.h.

◆ cutflow_histograms_v2

std::vector<std::tuple< TH1D *, std::vector<int*>, std::vector<float*>, std::function<float()> > > RooUtil::Cutflow::cutflow_histograms_v2

Definition at line 129 of file anautil.h.

◆ cutflow_nofill_cut_list

std::vector<TString> RooUtil::Cutflow::cutflow_nofill_cut_list

Definition at line 131 of file anautil.h.

◆ cutlists

std::map<TString, std::vector<TString> > RooUtil::Cutflow::cutlists

Definition at line 138 of file anautil.h.

◆ cutsysts

std::vector<TString> RooUtil::Cutflow::cutsysts

Definition at line 135 of file anautil.h.

◆ cuttree

CutTree RooUtil::Cutflow::cuttree

Definition at line 115 of file anautil.h.

Referenced by Cutflow().

◆ cuttreelists

std::map<TString, std::vector<CutTree*> > RooUtil::Cutflow::cuttreelists

Definition at line 139 of file anautil.h.

◆ cuttreemap

std::map< std::string , CutTree*> RooUtil::Cutflow::cuttreemap

Definition at line 117 of file anautil.h.

Referenced by Cutflow().

◆ dosavettreex

bool RooUtil::Cutflow::dosavettreex

Definition at line 143 of file anautil.h.

◆ doskipsysthist

bool RooUtil::Cutflow::doskipsysthist

Definition at line 142 of file anautil.h.

◆ iseventlistbooked

bool RooUtil::Cutflow::iseventlistbooked

Definition at line 140 of file anautil.h.

◆ last_active_cut

CutTree* RooUtil::Cutflow::last_active_cut

Definition at line 116 of file anautil.h.

◆ ofile

TFile* RooUtil::Cutflow::ofile

Definition at line 132 of file anautil.h.

Referenced by setTFile().

◆ rawcutflow_histograms

std::map< TString , TH1D *> RooUtil::Cutflow::rawcutflow_histograms

Definition at line 119 of file anautil.h.

◆ rawcutflow_histograms_v2

std::vector<std::tuple< TH1D *, std::vector<int*> > > RooUtil::Cutflow::rawcutflow_histograms_v2

Definition at line 130 of file anautil.h.

◆ seterrorcount

int RooUtil::Cutflow::seterrorcount

Definition at line 141 of file anautil.h.

◆ systs

std::vector<TString> RooUtil::Cutflow::systs

Definition at line 136 of file anautil.h.

◆ systs_funcs

std::map<TString, std::function<float()> > RooUtil::Cutflow::systs_funcs

Definition at line 137 of file anautil.h.

◆ t

TTree* RooUtil::Cutflow::t

Definition at line 133 of file anautil.h.

◆ tx

TTreeX* RooUtil::Cutflow::tx

Definition at line 134 of file anautil.h.

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Cutflow() [1/2]

◆ Cutflow() [2/2]

◆ ~Cutflow()

Member Function Documentation

◆ addCut()

◆ addCutSyst()

◆ addCutToLastActiveCut()

◆ addCutToSkipCutflowList()

◆ addToCutTreeMap()

◆ addWgtSyst()

◆ book2DHistogram()

◆ book2DVecHistogram() [1/2]

◆ book2DVecHistogram() [2/2]

◆ bookCutflowHistograms()

◆ bookCutflowHistograms_v1()

◆ bookCutflowHistograms_v2()

◆ bookCutflows()

◆ bookCutflowsForRegions()

◆ bookCutflowTree()

◆ bookEventLists()

◆ bookHistogram() [1/2]

◆ bookHistogram() [2/2]

◆ bookHistograms() [1/2]

◆ bookHistograms() [2/2]

◆ bookHistogramsForCut()

◆ bookHistogramsForCutAndAbove()

◆ bookHistogramsForCutAndBelow()

◆ bookHistogramsForEndCuts()

◆ bookVecHistogram() [1/2]

◆ bookVecHistogram() [2/2]

◆ copyAndEditCuts()

◆ createWgtSystBranches()

◆ fill()

◆ fillCutflow()

◆ fillCutflow_v2()

◆ fillCutflows()

◆ fillCutflows_v1()

◆ fillCutflows_v2()

◆ fillHistograms()

◆ filterCuts()

◆ getCut()

◆ printCuts()

◆ printSetFunctionError()

◆ removeCut()

◆ saveCutflows()

◆ saveHistograms()

◆ saveOutput()

◆ saveTTreeX()

◆ setCut()

◆ setCutLists()

◆ setCutSyst()

◆ setEventID()

◆ setHistsAxesExtendable()

◆ setLastActiveCut()

◆ setSaveTTreeX()

◆ setSkipSystematicHistograms()

◆ setTFile()

◆ setVariable()

◆ setWgtSyst()

Member Data Documentation

◆ booked_2dhistograms

◆ booked_2dhistograms_nominal_keys

◆ booked_histograms

◆ booked_histograms_nominal_keys

◆ cutflow_booked

◆ cutflow_histograms

◆ cutflow_histograms_v2

◆ cutflow_nofill_cut_list

◆ cutlists

◆ cutsysts

◆ cuttree

◆ cuttreelists

◆ cuttreemap

◆ dosavettreex

◆ doskipsysthist

◆ iseventlistbooked