#include <TriggerAnalyser.h>

Inheritance diagram for analysis::tools::TriggerAnalyser:

Public Member Functions
virtual bool	analysisWithTrigger ()

bool	l1tJetsAnalysis () const

bool	l1tMuonsAnalysis () const

virtual bool	selectionHLT ()

virtual bool	selectionL1 ()

virtual bool	selectionTrigger ()

virtual bool	selectionTriggerEmulated (const bool &, const bool &, const std::string &, const int &, const float &, const float &)

	TriggerAnalyser ()
	constructors More...

	TriggerAnalyser (int argc, char *argv[])

std::vector< std::shared_ptr< TriggerObject > >	triggerObjectsCaloJets ()

std::vector< std::shared_ptr< TriggerObject > >	triggerObjectsL1Jets ()

std::vector< std::shared_ptr< TriggerObject > >	triggerObjectsL1Muons ()

std::vector< std::shared_ptr< TriggerObject > >	triggerObjectsL3Muons ()

std::vector< std::shared_ptr< TriggerObject > >	triggerObjectsPFJets ()

	~TriggerAnalyser ()
	desctructor More...

Public Member Functions inherited from analysis::tools::BaseAnalyser
void	actionApplyPileupWeight (const int &var=0)
	apply pileup weight given a systematic variation More...

void	analyserTree ()
	root tree More...

std::shared_ptr< Analysis >	analysis ()
	returns pointer to Analysis object More...

	BaseAnalyser ()
	constructor More...

	BaseAnalyser (int argc, char *argv[])
	constructor More...

std::string	basename (const std::string &)
	returns the basename of a path More...

std::map< std::string, std::shared_ptr< TGraphAsymmErrors > >	btagEfficiencies () const

std::shared_ptr< Config >	config ()
	returns pointer to Config object More...

float	crossSection () const
	returns cross section More...

int	cutflow ()
	get cutflow index More...

void	cutflow (const int &)
	set cutflow index More...

void	cutflow (const std::string &label, const bool &ok=true)
	create and update cutflow entry in the cutflow histogram More...

virtual bool	event (const int &)
	event entry to be readout and processed More...

void	fillAnalyserTree ()
	fill root tree More...

void	fillPileupHistogram ()
	fills pileup histogram More...

void	generatorWeight ()
	generator weight More...

bool	genJetsAnalysis () const
	returns whether analysis of gen jets can be done More...

bool	genParticlesAnalysis () const
	returns whether analysis of gen particles can be done More...

std::shared_ptr< TH1F >	histogram (const std::string &)
	returns a given 1D histogram More...

std::map< std::string, std::shared_ptr< TH1F > >	histograms ()
	returns all 1D histograms More...

virtual void	histograms (const std::string &, const int &)
	create n histograms of a given type More...

int	nEvents ()
	number of events to be processed More...

std::shared_ptr< TFile >	output ()
	returns pointer to output root file More...

void	pileupHistogram ()
	creates pileup histogram More...

float	pileupWeight (const float &truepu, const int &var) const
	returns pileup weight given the true pileup and uncertainty variation in values of sigma More...

std::shared_ptr< PileupWeight >	pileupWeights () const
	returns pointer to pileup weights (MC-only) More...

void	scale (const float &)
	sets a scale More...

int	seed ()
	returns a seed for random number generator More...

int	seed (const std::string &)
	reads a file containing a seed and returns the seed or -1 if fails More...

void	seed (const int &)
	sets a seed for random number generators More...

bool	triggerEmulated (const std::string &)
	return status of an emulated l1 trigger More...

virtual bool	triggerEmulation (const std::string &, const int &, const float &, const float &, const std::string &)
	emulate l1 trigger More...

std::map< std::string, bool >	triggersEmulated ()
	return status of all emulated l1 triggers More...

float	trueInteractions () const
	returns true number of interactions More...

float	weight ()
	returns event weight More...

void	weight (const float &)
	sets event weight More...

virtual	~BaseAnalyser ()
	desctructor More...

Protected Attributes
bool	l1tjetsanalysis_

bool	l1tmuonsanalysis_

bool	triggeranalysis_

Protected Attributes inherited from analysis::tools::BaseAnalyser
std::shared_ptr< TTree >	analyser_tree_
	output tree More...

std::shared_ptr< Analysis >	analysis_
	Analysis objects. More...

std::map< std::string, std::shared_ptr< TGraphAsymmErrors > >	btageff_

std::shared_ptr< Config >	config_
	Config objects. More...

int	cutflow_
	Cutflow index. More...

bool	genjetsanalysis_

bool	genpartsanalysis_

std::map< std::string, std::shared_ptr< TH1F > >	h1_
	1D histograms mapping More...

std::map< std::string, std::shared_ptr< TH2F > >	h2_
	2D histograms mapping More...

std::shared_ptr< TFile >	hout_
	output root file More...

bool	isData_
	flag for DATA sample More...

bool	isMC_
	flag for MC sample More...

bool	istruepu_
	true pileup More...

std::string	puw_label_
	pileup weight label More...

std::shared_ptr< PileupWeight >	puweights_
	pileup weight More...

float	scale_
	overall scaling More...

int	seed_
	seed value More...

std::map< std::string, bool >	trg_emul_
	emulated triggers More...

float	weight_
	event weight More...

float	xsection_
	cross section More...

Detailed Description

Definition at line 35 of file TriggerAnalyser.h.

Constructor & Destructor Documentation

TriggerAnalyser::TriggerAnalyser ( )

constructors

Definition at line 20 of file TriggerAnalyser.cc.

21 {

22 }

TriggerAnalyser::TriggerAnalyser	(	int	argc,
		char *	argv[]
	)

Definition at line 25 of file TriggerAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_, analysis::tools::BaseAnalyser::config_, l1tjetsanalysis_, l1tmuonsanalysis_, and triggeranalysis_.

                                                         : BaseAnalyser(argc,argv)
 {
 
    triggeranalysis_ = false;
    
    if ( config_->triggerResults() != "" )
       triggeranalysis_  = analysis_->triggerResults(config_->triggerResults());
    
    if ( config_->triggerObjectsDir() != "" )
    {
       // online jets
       if ( config_->triggerObjectsL1Jets() != "l1tJets")
          analysis_->addTree<TriggerObject> (config_->triggerObjectsL1Jets()  ,Form("%s/%s", config_->triggerObjectsDir().c_str(),config_->triggerObjectsL1Jets().c_str()));
       analysis_->addTree<TriggerObject> (config_->triggerObjectsCaloJets(),Form("%s/%s", config_->triggerObjectsDir().c_str(),config_->triggerObjectsCaloJets().c_str()));
       analysis_->addTree<TriggerObject> (config_->triggerObjectsPFJets()  ,Form("%s/%s", config_->triggerObjectsDir().c_str(),config_->triggerObjectsPFJets().c_str()));
       // online b jets
       analysis_->addTree<TriggerObject> (config_->triggerObjectsBJets(),Form("%s/%s", config_->triggerObjectsDir().c_str(),config_->triggerObjectsBJets().c_str()));
       // online muons
       analysis_->addTree<TriggerObject> (config_->triggerObjectsL1Muons(),Form("%s/%s",config_->triggerObjectsDir().c_str(),config_->triggerObjectsL1Muons().c_str()));
       analysis_->addTree<TriggerObject> (config_->triggerObjectsL3Muons(),Form("%s/%s",config_->triggerObjectsDir().c_str(),config_->triggerObjectsL3Muons().c_str()));
    }
    if ( config_ -> l1tJetsCollection() != "")
    {
       l1tjetsanalysis_ = ( analysis_ -> addTree<L1TJet> ("l1tJets",config_ -> l1tJetsCollection()) != nullptr );
    }
    if ( config_ -> l1tMuonsCollection() != "")
    {
       l1tmuonsanalysis_ = ( analysis_ -> addTree<L1TMuon> ("l1tMuons",config_ -> l1tMuonsCollection()) != nullptr );
    }
 
 }

TriggerAnalyser::~TriggerAnalyser ( )

desctructor

Definition at line 57 of file TriggerAnalyser.cc.

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

Member Function Documentation

bool TriggerAnalyser::analysisWithTrigger ( )

virtual

Definition at line 184 of file TriggerAnalyser.cc.

References triggeranalysis_.

 {
    return triggeranalysis_;
 }

bool TriggerAnalyser::l1tJetsAnalysis ( ) const

Definition at line 235 of file TriggerAnalyser.cc.

References l1tjetsanalysis_.

Referenced by analysis::tools::Analyser::event().

 {
    return l1tjetsanalysis_;
 }

bool TriggerAnalyser::l1tMuonsAnalysis ( ) const

Definition at line 240 of file TriggerAnalyser.cc.

References l1tmuonsanalysis_.

Referenced by analysis::tools::Analyser::event().

 {
    return l1tmuonsanalysis_;
 }

bool TriggerAnalyser::selectionHLT ( )

virtual

Definition at line 133 of file TriggerAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_, analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, and analysis::tools::BaseAnalyser::weight_.

Referenced by main(), and selectionTrigger().

 {
    if ( config_->hltPath_ == "" ) return true;
    
    ++cutflow_;
    if ( ! analysis_->triggerResult(config_->hltPath_) ) return false;
    
    if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
       h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,(config_->hltPath_).c_str());
    
    h1_["cutflow"] -> Fill(cutflow_,weight_);
 
    return true;
 }

bool TriggerAnalyser::selectionL1 ( )

virtual

Definition at line 148 of file TriggerAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_, analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, and analysis::tools::BaseAnalyser::weight_.

Referenced by main(), and selectionTrigger().

 {
    if ( config_->l1Seed_ == "" ) return true;
    
    ++cutflow_;
    if ( ! analysis_->triggerResult(config_->l1Seed_)  ) return false;
    
    if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
       h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,(config_->l1Seed_).c_str());
    
    h1_["cutflow"] -> Fill(cutflow_,weight_);
 
    return true;
 }

bool TriggerAnalyser::selectionTrigger ( )

virtual

Emulated triggers

Definition at line 70 of file TriggerAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, selectionHLT(), selectionL1(), and selectionTriggerEmulated().

Referenced by main().

 {
    bool hlt = selectionHLT();
    bool l1  = selectionL1();
    
    // L1 muon trigger
    bool l1muon = true;
    if ( config_->triggerEmulateL1Muons() != "" &&  config_->triggerEmulateL1MuonsNMin() > 0 && config_->triggerObjectsL1Muons() != "" )
    {
       int nmin = config_->triggerEmulateL1MuonsNMin();
       float ptmin = config_->triggerEmulateL1MuonsPtMin();
       float etamax = config_->triggerEmulateL1MuonsEtaMax();
       l1muon = selectionTriggerEmulated(l1,hlt,config_->triggerEmulateL1Muons(),nmin,ptmin,etamax);
    }
    // L3 muon trigger
    bool l3muon = true;
    if ( config_->triggerEmulateL3Muons() != "" &&  config_->triggerEmulateL3MuonsNMin() > 0  && config_->triggerObjectsL3Muons() != "" )
    {
       int nmin = config_->triggerEmulateL3MuonsNMin();
       float ptmin = config_->triggerEmulateL3MuonsPtMin();
       float etamax = config_->triggerEmulateL3MuonsEtaMax();
       l3muon = selectionTriggerEmulated(l1,hlt,config_->triggerEmulateL3Muons(),nmin,ptmin,etamax);
    }
       
    // L1 jet trigger
    bool l1jet = true;
    if ( config_->triggerEmulateL1Jets() != "" &&  config_->triggerEmulateL1JetsNMin() > 0 && config_->triggerObjectsL1Jets() != "" )
    {
       int nmin = config_->triggerEmulateL1JetsNMin();
       float ptmin = config_->triggerEmulateL1JetsPtMin();
       float etamax = config_->triggerEmulateL1JetsEtaMax();
       l1jet = selectionTriggerEmulated(l1,hlt,config_->triggerEmulateL1Jets(),nmin,ptmin,etamax);
    }
    
    // Calo jet trigger
    bool calojet = true;
    if ( config_->triggerEmulateCaloJets() != "" &&  config_->triggerEmulateCaloJetsNMin() > 0 && config_->triggerObjectsCaloJets() != "" )
    {
       int nmin = config_->triggerEmulateCaloJetsNMin();
       float ptmin = config_->triggerEmulateCaloJetsPtMin();
       float etamax = config_->triggerEmulateCaloJetsEtaMax();
       calojet = selectionTriggerEmulated(l1,hlt,config_->triggerEmulateCaloJets(),nmin,ptmin,etamax);
    }
    
    
    // PF jet trigger
    bool pfjet = true;
    if ( config_->triggerEmulatePFJets() != "" &&  config_->triggerEmulatePFJetsNMin() > 0 && config_->triggerObjectsPFJets() != "" )
    {
       int nmin = config_->triggerEmulatePFJetsNMin();
       float ptmin = config_->triggerEmulatePFJetsPtMin();
       float etamax = config_->triggerEmulatePFJetsEtaMax();
       pfjet = selectionTriggerEmulated(l1,hlt,config_->triggerEmulatePFJets(),nmin,ptmin,etamax);
    }
    
    
    bool emul = l1muon && l3muon && l1jet && calojet && pfjet;
    
    return (hlt && l1 && emul);
    
 }

bool TriggerAnalyser::selectionTriggerEmulated	(	const bool &	l1,
		const bool &	hlt,
		const std::string &	name,
		const int &	nmin,
		const float &	ptmin,
		const float &	etamax
	)

virtual

Definition at line 163 of file TriggerAnalyser.cc.

References analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, analysis::tools::BaseAnalyser::triggerEmulated(), and analysis::tools::BaseAnalyser::weight_.

Referenced by selectionTrigger().

 {
    
    ++cutflow_;
    
    if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
       h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Emulated: %s (n >= %d, pT >= %4.1f GeV, |eta| <= %4.1f)",name.c_str(),nmin,ptmin,etamax));
    
    if ( ! ( l1 && hlt ) ) return false;
    if ( ! triggerEmulated(name) ) return false;
    
    h1_["cutflow"] -> Fill(cutflow_,weight_);
 
    return true;
 }

std::vector< std::shared_ptr< TriggerObject > > TriggerAnalyser::triggerObjectsCaloJets ( )

Definition at line 199 of file TriggerAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_, and analysis::tools::BaseAnalyser::config_.

 {
    auto collection = analysis_->collection<TriggerObject>(config_->triggerObjectsCaloJets());
    std::vector< std::shared_ptr<TriggerObject> > objects;
    for ( int j = 0 ; j < collection->size() ; ++j )
       objects.push_back(std::make_shared<TriggerObject>(collection->at(j)));
    return objects;
 }

std::vector< std::shared_ptr< TriggerObject > > TriggerAnalyser::triggerObjectsL1Jets ( )

Definition at line 190 of file TriggerAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_, and analysis::tools::BaseAnalyser::config_.

Referenced by main().

 {
    auto collection = analysis_->collection<TriggerObject>(config_->triggerObjectsL1Jets());
    std::vector< std::shared_ptr<TriggerObject> > objects;
    for ( int j = 0 ; j < collection->size() ; ++j )
       objects.push_back(std::make_shared<TriggerObject>(collection->at(j)));
    return objects;
 }

std::vector< std::shared_ptr< TriggerObject > > TriggerAnalyser::triggerObjectsL1Muons ( )

Definition at line 217 of file TriggerAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_, and analysis::tools::BaseAnalyser::config_.

 {
    auto collection = analysis_->collection<TriggerObject>(config_->triggerObjectsL1Muons());
    std::vector< std::shared_ptr<TriggerObject> > objects;
    for ( int j = 0 ; j < collection->size() ; ++j )
       objects.push_back(std::make_shared<TriggerObject>(collection->at(j)));
    return objects;
 }

std::vector< std::shared_ptr< TriggerObject > > TriggerAnalyser::triggerObjectsL3Muons ( )

Definition at line 226 of file TriggerAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_, and analysis::tools::BaseAnalyser::config_.

 {
    auto collection = analysis_->collection<TriggerObject>(config_->triggerObjectsL3Muons());
    std::vector< std::shared_ptr<TriggerObject> > objects;
    for ( int j = 0 ; j < collection->size() ; ++j )
       objects.push_back(std::make_shared<TriggerObject>(collection->at(j)));
    return objects;
 }

std::vector< std::shared_ptr< TriggerObject > > TriggerAnalyser::triggerObjectsPFJets ( )

Definition at line 208 of file TriggerAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_, and analysis::tools::BaseAnalyser::config_.

 {
    auto collection = analysis_->collection<TriggerObject>(config_->triggerObjectsPFJets());
    std::vector< std::shared_ptr<TriggerObject> > objects;
    for ( int j = 0 ; j < collection->size() ; ++j )
       objects.push_back(std::make_shared<TriggerObject>(collection->at(j)));
    return objects;
 }