analysis::tools::Analysis Class Reference

#include <Analysis.h>

Public Member Functions
void	addBtagEfficiencies (const std::string &)
template<class Object >
std::shared_ptr< Collection< Object > >	addCollection (const std::string &unique_name)
template<class Object >
std::shared_ptr< Collection< Object > >	addCollection (const Collection< Object > &collection)
template<class Object >
std::shared_ptr< Collection< Object > >	addCollection (const std::vector< Object > &objects, const std::string &unique_name)
template<class Object >
std::shared_ptr< PhysicsObjectTree< Object > >	addTree (const std::string &unique_name, const std::string &path)
	Analysis (const std::string &inputFilelist, const std::string &evtinfo="MssmHbb/Events/EventInfo")
std::shared_ptr< BTagCalibrationReader >	btagCalibration (const std::string &tagger, const std::string &filename, const std::string &wp, const std::string &sysType="central", const std::vector< std::string > &otherSysTypes={"up","down"})
std::shared_ptr< BTagCalibrationReader >	btagCalibration ()
void	btagEfficienciesAlgo (const std::string &)
void	btagEfficienciesFlavour (const std::string &)
float	btagEfficiency (const analysis::tools::Jet &, const int &rank=0)
template<class Object >
std::shared_ptr< Collection< Object > >	collection (const std::string &unique_name)
double	crossSection ()
double	crossSection (const std::string &title)
int	crossSections (const std::string &path)
template<class Object >
void	defaultCollection (const std::string &unique_name)
template<class Object >
std::string	defaultCollection ()
void	event (const int &event, const bool &addCollections=true)
int	event ()
FilterResults	eventFilter (const std::string &path)
std::string	fileFullName ()
std::string	fileName ()
FilterResults	generatorFilter (const std::string &path)
double	genScale ()
double	genWeight ()
float	instantLumi ()
bool	isMC ()
std::shared_ptr< JetResolutionInfo >	jetResolutionInfo (const std::string &, const std::string &)
void	listCrossSections ()
void	listGeneratorFilter ()
double	luminosity ()
double	luminosity (const std::string &title)
float	lumiPileup ()
int	lumiSection ()
template<class Object1 , class Object2 >
void	match (const std::string &collection, const std::string &match_collection, const float &deltaR=0.5)
template<class Object1 , class Object2 >
void	match (const std::string &collection, const std::vector< std::string > &match_collections, const float &deltaR=0.5)
std::shared_ptr< MuonIdWeight >	muonIDWeights (const std::string &)
int	nPileup ()
float	nTruePileup ()
int	numberEvents ()
PDF	pdf ()
std::shared_ptr< PileupWeight >	pileupWeights (const std::string &)
int	processJsonFile (const std::string &fileName="goodJson.txt")
double	rho ()
	fixedGridRhoAll More...
int	run ()
float	scaleLuminosity (const float &lumi)
int	seed (const std::string &)
	seed for random number generator read from a txt file given as a parameter More...
bool	selectJson ()
int	size ()
void	tag (const std::string &)
std::string	tag ()
template<class Object >
std::shared_ptr< PhysicsObjectTree< Object > >	tree (const std::string &unique_name)
int	triggerPrescale (const std::string &trig)
std::map< std::string, int >	triggerPrescale (const std::vector< std::string > &trigs)
bool	triggerResult (const std::string &trig)
bool	triggerResults (const std::string &path)
	~Analysis ()

Protected Member Functions
void	treeInit_ (const std::string &unique_name, const std::string &path)

Protected Attributes
std::shared_ptr< BTagCalibration >	btagcalib_
std::shared_ptr< BTagCalibrationReader >	btagcalibread_
std::string	btageff_algo_
std::string	btageff_flavour_
std::map< std::string, std::any >	c_any_
std::string	defaultGenParticle_
int	event_
FilterResults	evtfilter_
TFile *	fileBtagEff_
TFileCollection *	fileCollection_
TCollection *	fileList_
FilterResults	genfilter_
double	genScale_
double	genWeight_
std::map< int, std::vector< std::string > >	goodLumi_
std::map< std::string, TH2F * >	h2_btageff_
std::string	inputFilelist_
float	inst_lumi_
bool	is_mc_
std::shared_ptr< JetResolutionInfo >	jerinfo_
std::map< int, std::vector< int > >	json_
int	lumi_
float	lumi_pu_
std::shared_ptr< MuonIdWeight >	muonIDweights_
float	mylumi_
int	n_pu_
float	n_true_pu_
int	nevents_
PDF	pdf_
std::shared_ptr< PileupWeight >	puweights_
double	rho_
int	run_
std::map< std::string, std::any >	t_any_
TChain *	t_event_
TChain *	t_evtfilter_
TChain *	t_genfilter_
TChain *	t_triggerResults_
std::map< std::string, std::string >	t_type_
TChain *	t_xsection_
std::string	tag_
std::map< std::string, TChain * >	tree_
std::map< std::string, bool >	triggerResults_
std::map< std::string, int >	triggerResultsPS_
std::map< std::string, double >	xsections_

Detailed Description

Definition at line 60 of file Analysis.h.

Constructor & Destructor Documentation

Analysis::Analysis	(	const std::string &	inputFilelist,
		const std::string &	evtinfo = "MssmHbb/Events/EventInfo"
	)

Definition at line 25 of file Analysis.cc.

References btageff_algo_, btageff_flavour_, event_, fileCollection_, fileList_, genScale_, genWeight_, analysis::tools::PDF::id, inputFilelist_, inst_lumi_, is_mc_, lumi_, lumi_pu_, mylumi_, n_pu_, n_true_pu_, nevents_, pdf_, rho_, run_, t_event_, and analysis::tools::PDF::x.

{
   inputFilelist_  = inputFilelist;
   fileCollection_ = new TFileCollection("fileCollection","",inputFilelist.c_str());
   fileList_ = (TCollection*) fileCollection_->GetList();

   // event info (must be in the tree always)
   t_event_ = new TChain(evtinfo.c_str());
   t_event_ -> AddFileInfoList(fileList_);
   
   std::vector<std::string> branches;
   TObjArray * treeBranches = t_event_->GetListOfBranches();
   for ( int i = 0 ; i < treeBranches->GetEntries() ; ++i )
      branches.push_back(treeBranches->At(i)->GetName());
   
   t_event_ -> SetBranchAddress("event", &event_);
   t_event_ -> SetBranchAddress("run", &run_);
   t_event_ -> SetBranchAddress("lumisection", &lumi_);
   
   // For backward compatibility
   std::vector<std::string>::iterator it;
   it = std::find(branches.begin(),branches.end(),"nPileup");      if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &n_pu_);
   it = std::find(branches.begin(),branches.end(),"nTruePileup");  if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &n_true_pu_);
   
   it = std::find(branches.begin(),branches.end(),"lumiPileup");   if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &lumi_pu_);
   it = std::find(branches.begin(),branches.end(),"instantLumi");  if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &inst_lumi_);
   
   it = std::find(branches.begin(),branches.end(),"genWeight");    if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &genWeight_);
   it = std::find(branches.begin(),branches.end(),"genScale");     if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &genScale_);
   it = std::find(branches.begin(),branches.end(),"pdfid1");       if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &pdf_.id.first);
   it = std::find(branches.begin(),branches.end(),"pdfid2");       if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &pdf_.id.second);
   it = std::find(branches.begin(),branches.end(),"pdfx1");        if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &pdf_.x.first);
   it = std::find(branches.begin(),branches.end(),"pdfx2");        if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &pdf_.x.second);
   
   it = std::find(branches.begin(),branches.end(),"rho");          if ( it != branches.end() ) t_event_  -> SetBranchAddress( (*it).c_str(), &rho_);
   
//   t_event_ -> SetBranchAddress("nPileup", &n_pu_);
//   t_event_ -> SetBranchAddress("nTruePileup", &n_true_pu_);

   nevents_ = t_event_ -> GetEntries();

   t_event_ -> GetEntry(1); // Check whether it's mc simulation
   if (run_ == 1) {         // A bit stupid, but it's only solution that I found for the moment
     is_mc_ = true;
   } else {
     is_mc_ = false;
   }
   
   btageff_algo_    = "";
   btageff_flavour_ = "";
   
   mylumi_= -1.;

   
   //if(is_mc_) crossSection();

}

Analysis::~Analysis

(

)

Definition at line 83 of file Analysis.cc.

{
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

void Analysis::addBtagEfficiencies

(

const std::string &

filename

)

Definition at line 416 of file Analysis.cc.

References btageff_algo_, btageff_flavour_, fileBtagEff_, and h2_btageff_.

Referenced by main().

{
   fileBtagEff_ = new TFile(filename.c_str(),"OLD");
   std::string ftitle = fileBtagEff_->GetTitle();
   std::vector<std::string> field;
   if ( ftitle != "" )
   {
      boost::split(field, ftitle, boost::is_any_of(":"));
      btageff_flavour_ = field[0];
      btageff_algo_    = field[1];
   }
   
   TList * mylist = fileBtagEff_->GetListOfKeys();
   for ( int i = 0 ; i < mylist->GetSize() ; ++i  )
   {
      std::string className = ((TKey*) mylist -> At(i)) -> GetClassName();
      std::string objName   = ((TKey*) mylist -> At(i)) -> GetName();
      if ( className == "TH2F" )
         h2_btageff_[objName] = (TH2F*) fileBtagEff_->Get(objName.c_str());
//      if ( className == "TH2D" )
//         h2_btageff_[objName] = (TH2D*) fileBtagEff_->Get(objName.c_str());
   }
   
}

template<class Object >

std::shared_ptr< Collection< Object > > Analysis::addCollection

(

const std::string &

unique_name

)

Definition at line 297 of file Analysis.h.

      {
         // Still need to see how to deal with collections not originating from the ntuple,
         // e.g. a selected jets collection from the ntuple jets collection.
         
         // If tree does not exist, return NULL
         std::map<std::string, std::any >::iterator it = t_any_.find(unique_name);
         if ( it == t_any_.end() )
            return nullptr;
         
         auto tree = std::any_cast< std::shared_ptr< PhysicsObjectTree<Object> > > (t_any_[unique_name]);
         c_any_[unique_name] = std::shared_ptr< Collection<Object> > ( new Collection<Object>(tree -> collection()));
         std::shared_ptr< Collection<Object> > ret = std::any_cast< std::shared_ptr< Collection<Object> > > (c_any_[unique_name]);
         
         return ret;
      }

template<class Object >

std::shared_ptr< Collection< Object > > Analysis::addCollection

(

const Collection< Object > &

collection

)

Definition at line 315 of file Analysis.h.

References analysis::tools::Collection< Object >::name().

      {
         std::string unique_name = collection.name();
         t_any_[unique_name] = nullptr;
         c_any_[unique_name] = std::shared_ptr< Collection<Object> > ( new Collection<Object>(collection) );
         std::shared_ptr< Collection<Object> > ret = std::any_cast< std::shared_ptr< Collection<Object> > > (c_any_[unique_name]);
         
         return ret;
      }

template<class Object >

std::shared_ptr< Collection< Object > > Analysis::addCollection	(	const std::vector< Object > &	objects,
		const std::string &	unique_name
	)

Definition at line 326 of file Analysis.h.

      {
         Collection<Object> collection(objects,unique_name);
         t_any_[unique_name] = nullptr;
         c_any_[unique_name] = std::shared_ptr< Collection<Object> > ( new Collection<Object>(collection) );
         std::shared_ptr< Collection<Object> > ret = std::any_cast< std::shared_ptr< Collection<Object> > > (c_any_[unique_name]);
         return ret;
      }

template<class Object >

std::shared_ptr< PhysicsObjectTree< Object > > Analysis::addTree	(	const std::string &	unique_name,
		const std::string &	path
	)

Definition at line 273 of file Analysis.h.

References nano_94X_2017_rereco31Mar18-v1::name, and naf_mult_submit::type.

Referenced by main().

      {
         if ( path == "" || unique_name == "" ) return nullptr;
         this->treeInit_(unique_name,path);
         t_any_[unique_name] = std::shared_ptr< PhysicsObjectTree<Object> > ( new PhysicsObjectTree<Object>(tree_[unique_name], unique_name) );
         std::string type = boost::core::demangle(typeid(Object).name());
         std::vector<std::string> tmp;
         boost::split( tmp, type, boost::is_any_of("::"));
         t_type_[unique_name] = tmp.back();
         return std::any_cast< std::shared_ptr< PhysicsObjectTree<Object> > > (t_any_[unique_name]);
      }

std::shared_ptr< BTagCalibrationReader > Analysis::btagCalibration	(	const std::string &	tagger,
		const std::string &	filename,
		const std::string &	wp,
		const std::string &	sysType = "central",
		const std::vector< std::string > &	otherSysTypes = {"up", "down"}
	)

Definition at line 511 of file Analysis.cc.

References btagcalib_, btagcalibread_, BTagEntry::FLAV_B, BTagEntry::FLAV_C, BTagEntry::FLAV_UDSG, BTagEntry::OP_LOOSE, BTagEntry::OP_MEDIUM, BTagEntry::OP_RESHAPING, and BTagEntry::OP_TIGHT.

Referenced by main().

{
   std::string wps = wp;
   std::transform(wps.begin(), wps.end(), wps.begin(), ::tolower);
   
   BTagEntry::OperatingPoint op = BTagEntry::OP_MEDIUM;
   if ( wps == "loose" )    op = BTagEntry::OP_LOOSE;
   if ( wps == "medium" )   op = BTagEntry::OP_MEDIUM;
   if ( wps == "tight" )    op = BTagEntry::OP_TIGHT;
   if ( wps == "reshape" )  op = BTagEntry::OP_RESHAPING;
   
   btagcalib_     = std::shared_ptr<BTagCalibration>      ( new BTagCalibration(tagger,filename));
   btagcalibread_ = std::shared_ptr<BTagCalibrationReader>( new BTagCalibrationReader(op,sysType,otherSysTypes) );
   
   
   btagcalibread_ -> load(*btagcalib_,             // calibration instance
                     BTagEntry::FLAV_B,           // btag flavour - B
                     "comb");                     // measurement type   
   
   btagcalibread_ -> load(*btagcalib_,             // calibration instance
                     BTagEntry::FLAV_C,           // btag flavour - C
                     "comb");                     // measurement type   
   
   btagcalibread_ -> load(*btagcalib_,             // calibration instance
                     BTagEntry::FLAV_UDSG,        // btag flavour - UDSG
                     "incl");                     // measurement type   
   
   
   return btagcalibread_;
   
}

std::shared_ptr< BTagCalibrationReader > Analysis::btagCalibration

(

)

Definition at line 548 of file Analysis.cc.

References btagcalibread_.

{
   return btagcalibread_;
}

void Analysis::btagEfficienciesAlgo ( const std::string &  algo )

inline

Definition at line 394 of file Analysis.h.

Referenced by main().

{ btageff_algo_    = algo; }

void Analysis::btagEfficienciesFlavour ( const std::string &  flavour )

inline

Definition at line 395 of file Analysis.h.

Referenced by main().

{ btageff_flavour_ = flavour; }

float Analysis::btagEfficiency	(	const analysis::tools::Jet &	jet,
		const int &	rank = 0
	)

Definition at line 441 of file Analysis.cc.

References btageff_flavour_, analysis::tools::Candidate::eta(), analysis::tools::Jet::extendedFlavour(), analysis::tools::Jet::flavour(), h2_btageff_, and analysis::tools::Candidate::pt().

Referenced by main().

{
   float eff = 0.;
   std::string flav;
   std::string srank;
   if ( rank < 1 ) srank = "";
   else            srank = Form("%i",rank);
   if ( btageff_flavour_ == "Extended" || btageff_flavour_ == "extended" )
   {
      flav = jet.extendedFlavour();
      if ( flav == "udsg" ) flav = "l";
   }
   else
   {
      int iflav = jet.flavour(btageff_flavour_);
      if ( abs(iflav) == 5 )                flav = "b";
      if ( abs(iflav) == 4 )                flav = "c";
      if ( abs(iflav) < 4 || iflav == 21 )  flav = "l";
   }
   float pt = jet.pt();
   float eta = fabs(jet.eta());
   
   std::string hname = Form("h_%sjet%s_eff_pt_eta",flav.c_str(),srank.c_str());

   int bin = h2_btageff_[hname] -> FindBin(pt,eta);
   
   eff = h2_btageff_[hname] -> GetBinContent(bin);
         
   return eff;
}

template<class Object >

std::shared_ptr< Collection< Object > > Analysis::collection

(

const std::string &

unique_name

)

Definition at line 336 of file Analysis.h.

Referenced by main().

      {
         std::shared_ptr< Collection<Object> > ret = std::any_cast< std::shared_ptr< Collection<Object> > > (c_any_[unique_name]);
         return ret;
      }

double Analysis::crossSection

(

)

Definition at line 241 of file Analysis.cc.

Referenced by luminosity().

{
   return this -> crossSection("crossSection");
}

double Analysis::crossSection

(

const std::string &

title

)

Definition at line 245 of file Analysis.cc.

References t_xsection_, and xsections_.

{
   if ( t_xsection_ == NULL ) return -1.;
   return xsections_[xs];
}

int Analysis::crossSections

(

const std::string &

path

)

Definition at line 219 of file Analysis.cc.

References fileList_, t_xsection_, and xsections_.

{
   if ( path == "" ) return -2;
   t_xsection_  = new TChain(path.c_str());
   int ok = t_xsection_ -> AddFileInfoList(fileList_);
   if ( ok == 0 )
   {
      std::cout << "tree does not exist" << std::endl;
      return -1;
   }
   TObjArray * xsecBranches = t_xsection_->GetListOfBranches();
   for ( int i = 0 ; i < xsecBranches->GetEntries() ; ++i )
   {
      std::string branch = xsecBranches->At(i)->GetName();
      if ( branch == "run" ) continue;
      xsections_[branch] = 0;
      t_xsection_ -> SetBranchAddress(branch.c_str(), &xsections_[branch]);
   }
   t_xsection_ -> GetEntry(0);
   return 0;
}

template<class Object >

void Analysis::defaultCollection

(

const std::string &

unique_name

)

Definition at line 358 of file Analysis.h.

      { 
         if ( std::is_same<Object,GenParticle>::value ) defaultGenParticle_ = unique_name; 
      }

template<class Object >

std::string Analysis::defaultCollection

(

)

Definition at line 363 of file Analysis.h.

      { 
         std::string ret;
         if ( std::is_same<Object,GenParticle>::value ) ret = defaultGenParticle_ ;
         return ret; 
      }

void Analysis::event	(	const int &	event,
		const bool &	addCollections = true
	)

Definition at line 94 of file Analysis.cc.

References genScale_, genWeight_, analysis::tools::PDF::id, n_pu_, n_true_pu_, nano_94X_2017_rereco31Mar18-v1::name, pdf_, t_any_, t_event_, t_type_, tree(), naf_mult_submit::type, and analysis::tools::PDF::x.

Referenced by main().

{
   // Initialisation for backward compatibility
   n_pu_ = -1;
   n_true_pu_ = -1;
   
   genWeight_ = -1.;
   genScale_  = -1.;
   pdf_.id.first  = 0;
   pdf_.id.second = 0;
   pdf_.x.first  = -1.;
   pdf_.x.second = -1.;
   
   t_event_ -> GetEntry(event);
   if ( !addCollections) return;
   
   for ( auto & tree : t_any_ )
   {
      std::string name = tree.first;
      std::string type = t_type_[name];
      if ( type == "Jet" )            this->addCollection<Jet>(name);
      if ( type == "Muon" )           this->addCollection<Muon>(name);
      if ( type == "GenJet" )         this->addCollection<GenJet>(name);
      if ( type == "MET" )            this->addCollection<MET>(name);
      if ( type == "Vertex" )         this->addCollection<Vertex>(name);
      if ( type == "TriggerObject" )  this->addCollection<TriggerObject>(name);
      if ( type == "GenParticle" )    this->addCollection<GenParticle>(name);
      if ( type == "Candidate" )      this->addCollection<Candidate>(name);
      if ( type == "JetTag" )         this->addCollection<JetTag>(name);
      if ( type == "L1TMuon" )        this->addCollection<L1TMuon>(name);
      if ( type == "L1TJet" )         this->addCollection<L1TJet>(name);
      if ( type == "RecoMuon" )       this->addCollection<RecoMuon>(name);
      if ( type == "RecoTrack" )      this->addCollection<RecoTrack>(name);
   }
   
}

int Analysis::event ( )

inline

Definition at line 378 of file Analysis.h.

References event_.

{ return event_;     }

FilterResults Analysis::eventFilter

(

const std::string &

path

)

Definition at line 348 of file Analysis.cc.

References analysis::tools::FilterResults::efficiency, evtfilter_, fileList_, analysis::tools::FilterResults::filtered, t_evtfilter_, and analysis::tools::FilterResults::total.

Referenced by main().

{
   t_evtfilter_  = new TChain(path.c_str());
   t_evtfilter_ -> AddFileInfoList(fileList_);

   unsigned int ntotal;
   unsigned int nfiltered;
   unsigned int sumtotal = 0;
   unsigned int sumfiltered = 0;

   t_evtfilter_ -> SetBranchAddress("nEventsTotal", &ntotal);
   t_evtfilter_ -> SetBranchAddress("nEventsFiltered", &nfiltered);

   for ( int i = 0; i < t_evtfilter_->GetEntries(); ++i )
   {
      t_evtfilter_ -> GetEntry(i);
      sumtotal += ntotal;
      sumfiltered += nfiltered;
   }


   evtfilter_.total = sumtotal;
   evtfilter_.filtered = sumfiltered;
   evtfilter_.efficiency = float(sumfiltered)/sumtotal;

   return evtfilter_;
}

std::string Analysis::fileFullName ( )

inline

Definition at line 397 of file Analysis.h.

Referenced by fileName().

{ return std::string(t_event_ -> GetFile() -> GetName()) ;    }

std::string Analysis::fileName

(

)

Definition at line 553 of file Analysis.cc.

References fileFullName().

Referenced by main().

{
   std::string filename;
   filename = boost::filesystem::basename(this->fileFullName())+boost::filesystem::extension(this->fileFullName());
   return filename ;
}

FilterResults Analysis::generatorFilter

(

const std::string &

path

)

Definition at line 296 of file Analysis.cc.

References analysis::tools::FilterResults::efficiency, fileList_, analysis::tools::FilterResults::filtered, genfilter_, t_genfilter_, and analysis::tools::FilterResults::total.

{
   t_genfilter_  = new TChain(path.c_str());
   t_genfilter_ -> AddFileInfoList(fileList_);

   unsigned int ntotal;
   unsigned int nfiltered;
   unsigned int sumtotal = 0;
   unsigned int sumfiltered = 0;

   t_genfilter_ -> SetBranchAddress("nEventsTotal", &ntotal);
   t_genfilter_ -> SetBranchAddress("nEventsFiltered", &nfiltered);

   for ( int i = 0; i < t_genfilter_->GetEntries(); ++i )
   {
      t_genfilter_ -> GetEntry(i);
      sumtotal += ntotal;
      sumfiltered += nfiltered;
   }


   genfilter_.total = sumtotal;
   genfilter_.filtered = sumfiltered;
   genfilter_.efficiency = float(sumfiltered)/sumtotal;

   return genfilter_;
}

double Analysis::genScale ( )

inline

Definition at line 390 of file Analysis.h.

{ return genScale_;  }

double Analysis::genWeight ( )

inline

Definition at line 389 of file Analysis.h.

{ return genWeight_; }

float Analysis::instantLumi ( )

inline

Definition at line 387 of file Analysis.h.

{ return inst_lumi_;}

bool Analysis::isMC ( )

inline

Definition at line 381 of file Analysis.h.

{ return is_mc_ ;    }

std::shared_ptr< JetResolutionInfo > Analysis::jetResolutionInfo	(	const std::string &	f_jer,
		const std::string &	f_jersf
	)

Definition at line 489 of file Analysis.cc.

References jerinfo_.

Referenced by main().

{
   JetResolution res = JetResolution(f_jer);
   JetResolutionScaleFactor sf = JetResolutionScaleFactor(f_jersf);
   jerinfo_ = std::make_shared<JetResolutionInfo>(JetResolutionInfo{res,sf});
   return jerinfo_;
}

void Analysis::listCrossSections

(

)

Definition at line 273 of file Analysis.cc.

References t_xsection_, and xsections_.

{
   std::cout << "=======================================================" << std::endl;
   std::cout << "  CROSS SECTIONS" << std::endl;
   std::cout << "=======================================================" << std::endl;
   if ( t_xsection_ == NULL )
   {
      std::cout << "No cross section tree has been declared." << std::endl;
      std::cout << "=======================================================" << std::endl;
      std::cout << std::endl;
      std::cout << std::endl;
      return;
   }

   for ( auto& xs : xsections_ )
   {
      std::cout << xs.first << " = " << xs.second << " pb " << std::endl;
   }
   std::cout << "=======================================================" << std::endl;
   std::cout << std::endl;
   std::cout << std::endl;
}

void Analysis::listGeneratorFilter

(

)

Definition at line 324 of file Analysis.cc.

References analysis::tools::FilterResults::efficiency, analysis::tools::FilterResults::filtered, genfilter_, t_genfilter_, and analysis::tools::FilterResults::total.

{
   std::cout << "=======================================================" << std::endl;
   std::cout << "  GENERATOR FILTER" << std::endl;
   std::cout << "=======================================================" << std::endl;
   if ( t_genfilter_ == NULL )
   {
      std::cout << "No generator tree has been declared." << std::endl;
      std::cout << "=======================================================" << std::endl;
      std::cout << std::endl;
      std::cout << std::endl;
      return;
   }
   std::cout << "Total generated events      = " << genfilter_.total << std::endl;
   std::cout << "Filtered generated events   = " << genfilter_.filtered << std::endl;
   std::cout << "Generator Filter Efficiency = " << genfilter_.efficiency << std::endl;

   std::cout << "=======================================================" << std::endl;
   std::cout << std::endl;
   std::cout << std::endl;


}

double Analysis::luminosity

(

)

Definition at line 251 of file Analysis.cc.

References crossSection(), and nevents_.

Referenced by scaleLuminosity().

{
    return (nevents_ / this -> crossSection() );
}

double Analysis::luminosity

(

const std::string &

title

)

Definition at line 256 of file Analysis.cc.

References crossSection(), nevents_, and t_xsection_.

{
    if ( t_xsection_ == NULL ) return -1.;
    return (nevents_ / this -> crossSection(xs));
}

float Analysis::lumiPileup ( )

inline

Definition at line 386 of file Analysis.h.

{ return lumi_pu_;   }

int Analysis::lumiSection ( )

inline

Definition at line 380 of file Analysis.h.

References lumi_.

Referenced by main().

{ return lumi_ ;     }

template<class Object1 , class Object2 >

void Analysis::match	(	const std::string &	collection,
		const std::string &	match_collection,
		const float &	deltaR = 0.5
	)

Definition at line 343 of file Analysis.h.

Referenced by main().

      {
         if ( match_collection == "" ) return;
         auto o1 = std::any_cast< std::shared_ptr< Collection<Object1> > > (c_any_[collection]);
         auto o2 = std::any_cast< std::shared_ptr< Collection<Object2> > > (c_any_[match_collection]);
         o1->matchTo(o2->vectorCandidates(),o2->name(), deltaR);
      }

template<class Object1 , class Object2 >

void Analysis::match	(	const std::string &	collection,
		const std::vector< std::string > &	match_collections,
		const float &	deltaR = 0.5
	)

Definition at line 352 of file Analysis.h.

      {
         for ( auto & mc : match_collections )
            this->match<Object1,Object2>(collection, mc, deltaR);
      }

std::shared_ptr< MuonIdWeight > Analysis::muonIDWeights

(

const std::string &

f_muID

)

Definition at line 504 of file Analysis.cc.

References muonIDweights_.

{
   muonIDweights_ = std::make_shared<MuonIdWeight>(MuonIdWeight(f_muID));
   return muonIDweights_;
}

int Analysis::nPileup ( )

inline

Definition at line 383 of file Analysis.h.

Referenced by main().

{ return n_pu_;      }

float Analysis::nTruePileup ( )

inline

Definition at line 384 of file Analysis.h.

{ return n_true_pu_; }

int Analysis::numberEvents ( )

inline

Definition at line 376 of file Analysis.h.

{ return nevents_;   }

PDF Analysis::pdf ( )

inline

Definition at line 391 of file Analysis.h.

{ return pdf_;       }

std::shared_ptr< PileupWeight > Analysis::pileupWeights

(

const std::string &

f_pu

)

Definition at line 497 of file Analysis.cc.

References puweights_.

{
   puweights_ = std::make_shared<PileupWeight>(PileupWeight(f_pu));
   return puweights_;
}

int Analysis::processJsonFile

(

const std::string &

fileName = "goodJson.txt"

)

Definition at line 378 of file Analysis.cc.

References json_, and run().

Referenced by main().

{
   using boost::property_tree::ptree;
   ptree pt;
   read_json(fileName , pt);
   
   for (auto & element: pt)
   {
      int run = std::stoi(element.first);
      std::vector<int> lumiranges;
      for ( auto & property_array : element.second )
      {
         for ( auto & property : property_array.second )
         {
            lumiranges.push_back(property.second.get_value<int>());
         }
      }
      if ( lumiranges.size()%2 != 0 ) return -1;
//      std::sort(lumiranges.begin(), lumiranges.end());  // not really needed
      
      json_[run] = lumiranges;
   }
   return 0;
}

double Analysis::rho ( )

inline

fixedGridRhoAll

Definition at line 392 of file Analysis.h.

{ return rho_; }

int Analysis::run ( )

inline

Definition at line 379 of file Analysis.h.

References run_.

Referenced by main(), and processJsonFile().

{ return run_  ;     }

float Analysis::scaleLuminosity

(

const float &

lumi

)

Definition at line 262 of file Analysis.cc.

References luminosity(), and mylumi_.

{
   float lumiScale = 1.;
   mylumi_ = lumi;
   if ( mylumi_ < 0 ) return lumiScale;
   
   lumiScale  = mylumi_/this->luminosity();
   return lumiScale;
   
}

int Analysis::seed

(

const std::string &

name

)

seed for random number generator read from a txt file given as a parameter

Definition at line 561 of file Analysis.cc.

References f.

{
   int seed = 1;
   std::ifstream f(name.c_str(),std::ios_base::in);
   if ( ! f.good() )   return -1;
   
   f >> seed;
   f.close();
   if ( seed < 1 )     return -1;
   
   return seed;
}

bool Analysis::selectJson

(

)

Definition at line 403 of file Analysis.cc.

References json_, lumi_, and run_.

Referenced by main().

{
    bool isGood = false;
   
   std::vector<int> lumiranges = json_[run_];
   for ( size_t i = 0 ; i< lumiranges.size() ; i += 2 )
   {
      if ( lumi_ >= lumiranges.at(i) && lumi_ <= lumiranges.at(i+1) ) isGood = true;
   }
    return isGood;
}

int Analysis::size ( )

inline

Definition at line 377 of file Analysis.h.

Referenced by main().

{ return nevents_;   }

void Analysis::tag ( const std::string &  t )

inline

Definition at line 371 of file Analysis.h.

                                                          { std::cout << "Tag " << t << " has been defined." << std::endl ;
                                                                                 tag_ = t;     }

std::string Analysis::tag ( )

inline

Definition at line 373 of file Analysis.h.

{ return tag_ ; }

template<class Object >

std::shared_ptr< PhysicsObjectTree< Object > > Analysis::tree

(

const std::string &

unique_name

)

Definition at line 286 of file Analysis.h.

Referenced by event().

      {
         // If tree does not exist, return NULL
         std::map<std::string, std::any >::iterator it = t_any_.find(unique_name);
         if ( it == t_any_.end() )
            return nullptr;
         return std::any_cast< std::shared_ptr< PhysicsObjectTree<Object> > > (t_any_[unique_name]);
      }

void Analysis::treeInit_ ( const std::string &  unique_name, const std::string &  path  )

protected

Definition at line 136 of file Analysis.cc.

References fileList_, t_event_, and tree_.

{
   std::string treeTitle = ((TTree*) t_event_->GetFile()->Get(path.c_str())) -> GetTitle();
   tree_[unique_name] = new TChain(path.c_str(),treeTitle.c_str());
   tree_[unique_name] -> AddFileInfoList(fileList_);
   t_event_ -> AddFriend(tree_[unique_name]);

   treeTitle.erase(std::remove(treeTitle.begin(),treeTitle.end(),' '),treeTitle.end());
   std::string classname = treeTitle.substr(0,treeTitle.find_first_of("|"));
   std::string inputTag  = treeTitle.substr(treeTitle.find_first_of("|")+1);

}

int Analysis::triggerPrescale

(

const std::string &

trig

)

Definition at line 197 of file Analysis.cc.

References triggerResultsPS_.

Referenced by main(), and triggerPrescale().

{
//   if ( t_triggerResults_ == NULL ) return -1.;
   return triggerResultsPS_["ps_"+trig];
}

std::map< std::string, int > Analysis::triggerPrescale

(

const std::vector< std::string > &

trigs

)

Definition at line 203 of file Analysis.cc.

References triggerPrescale().

{
   std::map<std::string,int> ps;
   for ( auto & trig : trigs )
   {
      ps[trig] = triggerPrescale(trig);
   }
   return ps;
}

bool Analysis::triggerResult

(

const std::string &

trig

)

Definition at line 190 of file Analysis.cc.

References t_triggerResults_, and triggerResults_.

Referenced by main().

{
   if ( trig == "" ) return true;
   if ( t_triggerResults_ == NULL ) return false;
   return triggerResults_[trig];
}

bool Analysis::triggerResults

(

const std::string &

path

)

Definition at line 162 of file Analysis.cc.

References fileList_, t_event_, t_triggerResults_, triggerResults_, and triggerResultsPS_.

Referenced by main().

{
   t_triggerResults_  = new TChain(path.c_str());
   int ok = t_triggerResults_ -> AddFileInfoList(fileList_);
   t_event_ -> AddFriend(t_triggerResults_);
   if ( ok == 0 )
   {
      std::cout << "tree does not exist" << std::endl;
      return false;
   }
   TObjArray * triggerBranches = t_triggerResults_ -> GetListOfBranches();
   for ( int i = 0 ; i < triggerBranches->GetEntries() ; ++i )
   {
      std::string branch = triggerBranches->At(i)->GetName();
      if ( TString(branch).BeginsWith("ps") )
      {
         triggerResultsPS_[branch] = 1;
         t_triggerResults_ -> SetBranchAddress(branch.c_str(), &triggerResultsPS_[branch]);
      }
      else
      {
         triggerResults_[branch] = 0;
         t_triggerResults_ -> SetBranchAddress(branch.c_str(), &triggerResults_[branch]);
      }
   }
   return true;
}

Member Data Documentation

std::shared_ptr<BTagCalibration> analysis::tools::Analysis::btagcalib_

protected

Definition at line 191 of file Analysis.h.

Referenced by btagCalibration().

std::shared_ptr<BTagCalibrationReader> analysis::tools::Analysis::btagcalibread_

protected

Definition at line 192 of file Analysis.h.

Referenced by btagCalibration().

std::string analysis::tools::Analysis::btageff_algo_

protected

Definition at line 189 of file Analysis.h.

Referenced by addBtagEfficiencies(), and Analysis().

std::string analysis::tools::Analysis::btageff_flavour_

protected

Definition at line 188 of file Analysis.h.

Referenced by addBtagEfficiencies(), Analysis(), and btagEfficiency().

std::map<std::string, std::any > analysis::tools::Analysis::c_any_

protected

Definition at line 249 of file Analysis.h.

std::string analysis::tools::Analysis::defaultGenParticle_

protected

Definition at line 211 of file Analysis.h.

int analysis::tools::Analysis::event_

protected

Definition at line 213 of file Analysis.h.

Referenced by Analysis().

FilterResults analysis::tools::Analysis::evtfilter_

protected

Definition at line 208 of file Analysis.h.

Referenced by eventFilter().

TFile* analysis::tools::Analysis::fileBtagEff_

protected

Definition at line 186 of file Analysis.h.

Referenced by addBtagEfficiencies().

TFileCollection* analysis::tools::Analysis::fileCollection_

protected

Definition at line 178 of file Analysis.h.

Referenced by Analysis().

TCollection* analysis::tools::Analysis::fileList_

protected

Definition at line 179 of file Analysis.h.

Referenced by Analysis(), crossSections(), eventFilter(), generatorFilter(), treeInit_(), and triggerResults().

FilterResults analysis::tools::Analysis::genfilter_

protected

Definition at line 207 of file Analysis.h.

Referenced by generatorFilter(), and listGeneratorFilter().

double analysis::tools::Analysis::genScale_

protected

Definition at line 225 of file Analysis.h.

Referenced by Analysis(), and event().

double analysis::tools::Analysis::genWeight_

protected

Definition at line 224 of file Analysis.h.

Referenced by Analysis(), and event().

std::map<int,std::vector<std::string> > analysis::tools::Analysis::goodLumi_

protected

Definition at line 206 of file Analysis.h.

std::map<std::string,TH2F *> analysis::tools::Analysis::h2_btageff_

protected

Definition at line 187 of file Analysis.h.

Referenced by addBtagEfficiencies(), and btagEfficiency().

std::string analysis::tools::Analysis::inputFilelist_

protected

Definition at line 180 of file Analysis.h.

Referenced by Analysis().

float analysis::tools::Analysis::inst_lumi_

protected

Definition at line 222 of file Analysis.h.

Referenced by Analysis().

bool analysis::tools::Analysis::is_mc_

protected

Definition at line 216 of file Analysis.h.

Referenced by Analysis().

std::shared_ptr<JetResolutionInfo> analysis::tools::Analysis::jerinfo_

protected

Definition at line 194 of file Analysis.h.

Referenced by jetResolutionInfo().

std::map<int,std::vector<int> > analysis::tools::Analysis::json_

protected

Definition at line 255 of file Analysis.h.

Referenced by processJsonFile(), and selectJson().

int analysis::tools::Analysis::lumi_

protected

Definition at line 215 of file Analysis.h.

Referenced by Analysis(), and selectJson().

float analysis::tools::Analysis::lumi_pu_

protected

Definition at line 221 of file Analysis.h.

Referenced by Analysis().

std::shared_ptr<MuonIdWeight> analysis::tools::Analysis::muonIDweights_

protected

Definition at line 200 of file Analysis.h.

Referenced by muonIDWeights().

float analysis::tools::Analysis::mylumi_

protected

Definition at line 252 of file Analysis.h.

Referenced by Analysis(), and scaleLuminosity().

int analysis::tools::Analysis::n_pu_

protected

Definition at line 218 of file Analysis.h.

Referenced by Analysis(), and event().

float analysis::tools::Analysis::n_true_pu_

protected

Definition at line 219 of file Analysis.h.

Referenced by Analysis(), and event().

int analysis::tools::Analysis::nevents_

protected

Definition at line 230 of file Analysis.h.

Referenced by Analysis(), and luminosity().

PDF analysis::tools::Analysis::pdf_

protected

Definition at line 226 of file Analysis.h.

Referenced by Analysis(), and event().

std::shared_ptr<PileupWeight> analysis::tools::Analysis::puweights_

protected

Definition at line 197 of file Analysis.h.

Referenced by pileupWeights().

double analysis::tools::Analysis::rho_

protected

Definition at line 228 of file Analysis.h.

Referenced by Analysis().

int analysis::tools::Analysis::run_

protected

Definition at line 214 of file Analysis.h.

Referenced by Analysis(), and selectJson().

std::map<std::string, std::any > analysis::tools::Analysis::t_any_

protected

Definition at line 245 of file Analysis.h.

Referenced by event().

TChain* analysis::tools::Analysis::t_event_

protected

Definition at line 237 of file Analysis.h.

Referenced by Analysis(), event(), treeInit_(), and triggerResults().

TChain* analysis::tools::Analysis::t_evtfilter_

protected

Definition at line 236 of file Analysis.h.

Referenced by eventFilter().

TChain* analysis::tools::Analysis::t_genfilter_

protected

Definition at line 235 of file Analysis.h.

Referenced by generatorFilter(), and listGeneratorFilter().

TChain* analysis::tools::Analysis::t_triggerResults_

protected

Definition at line 238 of file Analysis.h.

Referenced by triggerResult(), and triggerResults().

std::map<std::string, std::string > analysis::tools::Analysis::t_type_

protected

Definition at line 246 of file Analysis.h.

Referenced by event().

TChain* analysis::tools::Analysis::t_xsection_

protected

Definition at line 234 of file Analysis.h.

Referenced by crossSection(), crossSections(), listCrossSections(), and luminosity().

std::string analysis::tools::Analysis::tag_

protected

Definition at line 183 of file Analysis.h.

std::map<std::string, TChain*> analysis::tools::Analysis::tree_

protected

Definition at line 242 of file Analysis.h.

Referenced by treeInit_().

std::map<std::string, bool> analysis::tools::Analysis::triggerResults_

protected

Definition at line 204 of file Analysis.h.

Referenced by triggerResult(), and triggerResults().

std::map<std::string, int> analysis::tools::Analysis::triggerResultsPS_

protected

Definition at line 205 of file Analysis.h.

Referenced by triggerPrescale(), and triggerResults().

std::map<std::string, double> analysis::tools::Analysis::xsections_

protected

Definition at line 203 of file Analysis.h.

Referenced by crossSection(), crossSections(), and listCrossSections().

---

The documentation for this class was generated from the following files:

• Analysis.h
• Analysis.cc