analysis::tools::JetAnalyser Class Reference

#include <JetAnalyser.h>

Inheritance diagram for analysis::tools::JetAnalyser:

Public Member Functions
virtual void	actionApplyBjetRegression ()
virtual void	actionApplyBjetRegression (const int &)
virtual float	actionApplyBtagSF (const int &, const bool &global_weight=true)
virtual void	actionApplyJER ()
virtual bool	analysisWithJets ()
float	btag (const Jet &jet, const std::string &algo)
	Returns the btag value of the jet for a given btag algorithm. More...
virtual ScaleFactors	btagSF (const int &, const std::string &)
virtual void	fillJetHistograms (const std::string &label="x")
	Fill the pre-defined histograms created by the jetHistograms() method. More...
virtual void	fillJetHistograms (const int &r=-1, const std::string &label="x", const float &sf=1., const bool &workflow=false)
	For experts only (fills pre-defined histograms for a specific jet, applying weights) More...
virtual float	getBtagSF (const int &)
	JetAnalyser ()
	constructors More...
	JetAnalyser (int argc, char *argv[])
virtual bool	jetCorrections ()
	multiple actions: apply JER and b-tag regression corrections More...
virtual void	jetHistograms (const int &n, const std::string &label="x")
	Creates pre-defined histograms in directory 'label' for analysis with 'n' jets. More...
virtual void	jetHistograms (const std::string &label="x")
	Creates pre-defined histograms in directory 'label' for analysis with config->nJetsMin() More...
std::vector< std::shared_ptr< Jet > >	jets ()
	vector of pointers of all jets from the "Jets" collection More...
void	jets (const std::string &col)
	Creates the "Jets" collection. More...
virtual void	jetSwap (const int &, const int &)
virtual bool	onlineBJetMatching (const int &)
virtual bool	onlineJetMatching (const int &)
std::vector< std::shared_ptr< Jet > >	selectedJets ()
	vector of pointers of the selectedJets More...
virtual bool	selectionBJet (const int &)
virtual bool	selectionBJetProbB (const int &)
virtual bool	selectionBJetProbBB (const int &)
virtual bool	selectionBJetProbC (const int &)
virtual bool	selectionBJetProbG (const int &)
virtual bool	selectionBJetProbLepB (const int &)
virtual bool	selectionBJetProbLight (const int &)
virtual bool	selectionDiJetMass (const int &, const int &)
virtual bool	selectionJet (const int &r)
	Given the ranking 'r' of a jet, it returns whether the jet passes the pt_min and |eta_max|, optionally pt_max, where the values of the thresholds pt_min and |eta_max|, pt_max are passed by configuration file. More...
virtual bool	selectionJet (const int &r, const float &pt_min, const float &eta_max, const float &pt_max=-1.)
	Given the ranking 'r' of a jet, it returns whether the jet passes the pt_min and |eta_max|, optionally pt_max, where the values of the thresholds pt_min and |eta_max|, pt_max (opt) are passed as arguments. More...
virtual bool	selectionJetDeta (const int &r1, const int &r2)
	Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_eta selection; the threshold delta is passed by the configuration file. More...
virtual bool	selectionJetDeta (const int &r1, const int &r2, const float &delta)
	Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_eta selection; the threshold delta is passed as an argument (if delta <= 0, |delta| = delta_min; else |delta| = delta_max ) More...
virtual bool	selectionJetDphi (const int &r1, const int &r2)
	Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_phi selection; the threshold delta is passed by the configuration file. More...
virtual bool	selectionJetDphi (const int &r1, const int &r2, const float &delta)
	Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_phi selection; the threshold delta is passed as an argument (if delta <= 0, |delta| = delta_min; else |delta| = delta_max ) More...
virtual bool	selectionJetDr (const int &r1, const int &r2)
	Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_R selection; the threshold delta is passed by the configuration file. More...
virtual bool	selectionJetDr (const int &r1, const int &r2, const float &delta)
	Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_R selection; the threshold delta is passed as an argument (if delta <= 0, |delta| = delta_min; else |delta| = delta_max ) More...
virtual bool	selectionJetId ()
virtual bool	selectionJetPileupId ()
virtual bool	selectionJetPtImbalance (const int &r1, const int &r2)
	Given the rankings r1 and r2 of two jets, it returns whether the jets passes the pt imbalance selection; the threshold delta is passed by the configuration file. More...
virtual bool	selectionJetPtImbalance (const int &r1, const int &r2, const float &delta)
	Given the rankings r1 and r2 of two jets, it returns whether the jets passes the pt imbalance selection; the threshold delta is passed as an argument (if delta <= 0, |delta| = delta_min; else |delta| = delta_max ) More...
virtual bool	selectionJetQGlikelihood (const int &, const float &)
virtual bool	selectionJetQGlikelihood (const int &)
virtual bool	selectionNJets ()
virtual bool	selectionNonBJet (const int &)
	~JetAnalyser ()
	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	applyjer_
std::map< std::string, std::shared_ptr< BTagCalibrationReader > >	bsf_reader_
std::vector< std::string >	flavours_
std::shared_ptr< JetResolutionInfo >	jerinfo_
std::vector< std::shared_ptr< Jet > >	jets_
bool	jetsanalysis_
int	n_hjets_
std::vector< std::shared_ptr< Jet > >	selectedJets_
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 JetAnalyser.h.

Constructor & Destructor Documentation

JetAnalyser::JetAnalyser

(

)

constructors

Definition at line 21 of file JetAnalyser.cc.

{
}

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

Definition at line 26 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_, applyjer_, bsf_reader_, analysis::tools::BaseAnalyser::config_, flavours_, analysis::tools::BaseAnalyser::genJetsAnalysis(), jerinfo_, and jetsanalysis_.

                                                : BaseAnalyser(argc,argv)
{
   // Physics objects
   // Jets
   jetsanalysis_  = ( analysis_->addTree<Jet> ("Jets",config_->jetsCollection()) != nullptr );
   
   applyjer_ = false;
   
   if ( config_->btagScaleFactors() != "" )
   {
      bsf_reader_["loose"]  = analysis_->btagCalibration(config_->btagAlgorithm(), config_->btagScaleFactors(), "loose"); 
      bsf_reader_["medium"] = analysis_->btagCalibration(config_->btagAlgorithm(), config_->btagScaleFactors(), "medium");
      bsf_reader_["tight"]  = analysis_->btagCalibration(config_->btagAlgorithm(), config_->btagScaleFactors(), "tight");
   }
   
   if ( config_->jerPtRes() != "" && config_->jerSF() != "" && this->genJetsAnalysis() ) // FIXME: check if files exist
   {
      jerinfo_ = analysis_->jetResolutionInfo(config_->jerPtRes(),config_->jerSF());
      applyjer_ = ( jerinfo_ != nullptr && jetsanalysis_ );
   }
   
   if ( config_->isMC() )
   {
      flavours_ = {"udsg","c","b"};
      if ( config_->useJetsExtendedFlavour() )
      {
         flavours_.push_back("cc");
         flavours_.push_back("bb");
      }
   }
//   histograms("jet",config_->nJetsMin());

}

JetAnalyser::~JetAnalyser

(

)

desctructor

Definition at line 60 of file JetAnalyser.cc.

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

Member Function Documentation

void JetAnalyser::actionApplyBjetRegression ( )

virtual

Definition at line 1044 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), and selectedJets_.

Referenced by jetCorrections(), and main().

{
   if ( ! config_->bRegression() ) return;
   
   for ( auto & j : selectedJets_ )
      j -> applyBjetRegression();

   cutflow("b jet energy regression");
}

void JetAnalyser::actionApplyBjetRegression ( const int &  r )

virtual

Definition at line 1054 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), and selectedJets_.

{
   if ( ! config_->bRegression() ) return;
   
   int j = r-1;
   selectedJets_[j]-> applyBjetRegression();
   
   cutflow(Form("Jet %d: b jet energy regression",r));
}

float JetAnalyser::actionApplyBtagSF ( const int &  r, const bool &  global_weight = true  )

virtual

Definition at line 1010 of file JetAnalyser.cc.

References btagSF(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), ReHLT_L1Trepack_HLTUser_DataNtuple::isMC, and analysis::tools::BaseAnalyser::weight_.

Referenced by analysis::mssmhbb::MssmHbbAnalyser::btagSelection(), and main().

{
   float sf = 1.;
   if ( ! config_-> isMC() || config_->btagScaleFactors() == "" ) return sf;  // will not apply btag SF
   if ( ! config_->signalRegion() && r == config_->revBtagJet() ) return sf;
   
   int j = r-1;
   std::string label = Form("Jet %d: btag SF applied (%s %s WP)",r,config_->btagAlgorithm().c_str(),config_->jetsBtagWP()[j].c_str());

   if ( config_->jetsBtagWP()[j] == "xxx" )
         label = Form("Jet %d: btag SF = 1 applied (%s %s WP)",r,config_->btagAlgorithm().c_str(),config_->jetsBtagWP()[j].c_str());
   
   if ( global_weight || config_->jetsBtagWP()[j] != "xxx" )  sf = this->btagSF(r,config_->jetsBtagWP()[j]).nominal;
   
   weight_ *= sf;
   
   cutflow(label);
   
   return sf;
   
}

void JetAnalyser::actionApplyJER ( )

virtual

Definition at line 992 of file JetAnalyser.cc.

References applyjer_, analysis::tools::BaseAnalyser::basename(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), analysis::tools::BaseAnalyser::isMC_, jerinfo_, jetsanalysis_, and selectedJets_.

Referenced by jetCorrections(), and main().

{
   if ( ! jetsanalysis_ ||  ! isMC_ ) return;
   
   std::string label = "WARNING: NO JER smearing (*** missing JER Info and/or GenJet collection ***)";
   
   if ( applyjer_ )
   {
      std::string bnpt = basename(config_->jerPtRes());
      std::string bnsf = basename(config_->jerSF());
      label = Form("JER smearing (%s,%s)",bnpt.c_str(),bnsf.c_str());
      for ( auto & j : selectedJets_ )
         j -> applyJER(*jerinfo_,0.2);
   }
   
   cutflow(label);
}

bool JetAnalyser::analysisWithJets ( )

virtual

Definition at line 72 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_, analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), analysis::tools::BaseAnalyser::genJetsAnalysis(), analysis::tools::BaseAnalyser::genParticlesAnalysis(), jets(), jets_, jetsanalysis_, selectedJets_, and analysis::tools::BaseAnalyser::triggerEmulation().

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

{
   jets_.clear();
   selectedJets_.clear();
   
   // trigger emulation
   // L1 jets
   std::string triggerObjectsL1Jets;
   if ( config_->triggerObjectsL1Jets() != "" )
   {
      triggerObjectsL1Jets = config_->triggerObjectsL1Jets();
      if ( config_->triggerEmulateL1Jets() != "" &&  config_->triggerEmulateL1JetsNMin() > 0 )
      {
         int nmin = config_->triggerEmulateL1JetsNMin();
         float ptmin = config_->triggerEmulateL1JetsPtMin();
         float etamax = config_->triggerEmulateL1JetsEtaMax();
         std::string newL1Jets = config_->triggerEmulateL1Jets();
         triggerEmulation(triggerObjectsL1Jets,nmin,ptmin,etamax,newL1Jets);
         triggerObjectsL1Jets = newL1Jets;
      }
   }
   
   // Calo jets
   std::string triggerObjectsCaloJets;
   if ( config_->triggerObjectsCaloJets() != "" )
   {
      triggerObjectsCaloJets = config_->triggerObjectsCaloJets();
      if ( config_->triggerEmulateCaloJets() != "" &&  config_->triggerEmulateCaloJetsNMin() > 0 )
      {
         int nmin = config_->triggerEmulateCaloJetsNMin();
         float ptmin = config_->triggerEmulateCaloJetsPtMin();
         float etamax = config_->triggerEmulateCaloJetsEtaMax();
         std::string newCaloJets = config_->triggerEmulateCaloJets();
         triggerEmulation(triggerObjectsCaloJets,nmin,ptmin,etamax,newCaloJets);
         triggerObjectsCaloJets = newCaloJets;
      }
   }

   // PF jets
   std::string triggerObjectsPFJets;
   if ( config_->triggerObjectsPFJets() != "" )
   {
      triggerObjectsPFJets = config_->triggerObjectsPFJets();
      if ( config_->triggerEmulatePFJets() != "" &&  config_->triggerEmulatePFJetsNMin() > 0 )
      {
         int nmin = config_->triggerEmulatePFJetsNMin();
         float ptmin = config_->triggerEmulatePFJetsPtMin();
         float etamax = config_->triggerEmulatePFJetsEtaMax();
         std::string newPFJets = config_->triggerEmulatePFJets();
         triggerEmulation(triggerObjectsPFJets,nmin,ptmin,etamax,newPFJets);
         triggerObjectsPFJets = newPFJets;
      }
   }

   
   if ( ! jetsanalysis_ ) return false;
   
   cutflow(Form("Using Jet collection: %s",(config_->jetsCollection()).c_str()));
   
   
   if ( config_->triggerObjectsL1Jets() != "" )
      analysis_->match<Jet,TriggerObject>("Jets",triggerObjectsL1Jets,config_-> triggerMatchL1JetsDrMax());
   if ( config_->triggerObjectsCaloJets() != "" )
      analysis_->match<Jet,TriggerObject>("Jets",triggerObjectsCaloJets,config_-> triggerMatchCaloJetsDrMax());
   if ( config_->triggerObjectsPFJets() != "" )
      analysis_->match<Jet,TriggerObject>("Jets",triggerObjectsPFJets,config_-> triggerMatchPFJetsDrMax());
   analysis_->match<Jet,TriggerObject>("Jets",config_->triggerObjectsBJets(),config_-> triggerMatchCaloBJetsDrMax());

   // std::shared_ptr< Collection<Jet> >
   auto jets = analysis_->collection<Jet>("Jets");
   
   if ( this->genParticlesAnalysis() && config_ -> useJetsExtendedFlavour() )
   {
      auto particles = analysis_->collection<GenParticle>("GenParticles");
      jets->associatePartons(particles,0.4,1.,config_->pythia8());
   }
   
   if ( this->genJetsAnalysis() )
   {
      auto genjets = analysis_->collection<GenJet>("GenJets");
      jets->addGenJets(genjets);
   }
   
   for ( int j = 0 ; j < jets->size() ; ++j )  jets_.push_back(std::make_shared<Jet>(jets->at(j)));
   
   selectedJets_ = jets_;
   
   return true;
}

float JetAnalyser::btag	(	const Jet &	jet,
		const std::string &	algo
	)

Returns the btag value of the jet for a given btag algorithm.

Definition at line 366 of file JetAnalyser.cc.

References analysis::tools::Jet::btag(), and analysis::tools::BaseAnalyser::config_.

Referenced by fillJetHistograms(), selectionBJet(), selectionBJetProbB(), selectionBJetProbBB(), selectionBJetProbC(), selectionBJetProbG(), selectionBJetProbLepB(), selectionBJetProbLight(), and selectionNonBJet().

{
   float btag;
   if ( config_->btagAlgorithm() == "csvivf" || config_->btagAlgorithm() == "csv" )   {      btag = jet.btag("btag_csvivf");   }
   else if ( config_->btagAlgorithm() == "deepcsv" )                            {      btag = jet.btag("btag_deepb") + jet.btag("btag_deepbb");   }
   else if ( config_->btagAlgorithm() == "deepbvsall" )                         {      btag = jet.btag("btag_deepbvsall");   }
   else if ( config_->btagAlgorithm() == "deepflavour" || config_->btagAlgorithm() == "deepjet" ) {      btag = jet.btag("btag_dfb") + jet.btag("btag_dfbb") + jet.btag("btag_dflepb");   }
   else                                                                   {      btag = -9999;   }
   
   return btag;
}

ScaleFactors JetAnalyser::btagSF ( const int &  r, const std::string &  wp  )

virtual

Definition at line 980 of file JetAnalyser.cc.

References bsf_reader_, analysis::tools::ScaleFactors::down, analysis::tools::ScaleFactors::nominal, selectedJets_, and analysis::tools::ScaleFactors::up.

Referenced by actionApplyBtagSF(), and getBtagSF().

{
   ScaleFactors sf;
   int j = r-1;
   sf.nominal = selectedJets_[j]->btagSF(bsf_reader_[wp]);
   sf.up      = selectedJets_[j]->btagSFup(bsf_reader_[wp]);
   sf.down    = selectedJets_[j]->btagSFdown(bsf_reader_[wp]);
   
   return sf;
}

void JetAnalyser::fillJetHistograms ( const std::string &  label = "x" )

virtual

Fill the pre-defined histograms created by the jetHistograms() method.

Definition at line 795 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), analysis::tools::Composite< A, B >::deltaEta(), analysis::tools::Composite< A, B >::deltaR(), analysis::tools::Candidate::eta(), analysis::tools::BaseAnalyser::h1_, analysis::tools::Candidate::m(), n_hjets_, analysis::tools::BaseAnalyser::output(), analysis::tools::Candidate::phi(), analysis::tools::Candidate::pt(), selectedJets_, and analysis::tools::BaseAnalyser::weight_.

Referenced by analysis::mssmhbb::MssmHbbAnalyser::endSelection(), and main().

{
   this->output()->cd();
   
   this->output()->cd(label.c_str());
   
   int n = n_hjets_;
   
   if ( n > config_->nJetsMin() ) n = config_->nJetsMin();
   
   h1_[Form("jet_hist_weight_%s",label.c_str())] -> Fill(0.,weight_);
   
   for ( int j = 0; j < n; ++j )
   {
      // fill each jet with addtional SF to weight = 1
      this->fillJetHistograms(j+1,label,1.);
      
      if ( config_ -> doDijet() )
      {
         for ( int k = j+1; k < n && j < n; ++k )
         {
            Composite<Jet,Jet> c_ij(*(selectedJets_[j]),*(selectedJets_[k]));
            
            h1_[Form("dptrel_jet%d%d_%s" , j+1,k+1,label.c_str())] -> Fill(fabs(selectedJets_[j]->pt()-selectedJets_[k]->pt())/selectedJets_[j]->pt(),weight_);
            h1_[Form("dpt_jet%d%d_%s" , j+1,k+1,label.c_str())]    -> Fill(fabs(selectedJets_[j]->pt()-selectedJets_[k]->pt()),weight_);
            h1_[Form("dr_jet%d%d_%s",j+1,k+1,label.c_str())]       -> Fill(c_ij.deltaR(),weight_);
            h1_[Form("deta_jet%d%d_%s",j+1,k+1,label.c_str())]     -> Fill(c_ij.deltaEta(),weight_);
            h1_[Form("dphi_jet%d%d_%s",j+1,k+1,label.c_str())]     -> Fill(fabs(selectedJets_[j]->deltaPhi(*selectedJets_[k])));
            h1_[Form("pt_jet%d%d_%s",j+1,k+1,label.c_str())]       -> Fill(c_ij.pt(),weight_);
            h1_[Form("eta_jet%d%d_%s",j+1,k+1,label.c_str())]      -> Fill(c_ij.eta(),weight_);
            h1_[Form("phi_jet%d%d_%s",j+1,k+1,label.c_str())]      -> Fill(c_ij.phi()*180./acos(-1.),weight_);
            if ( config_->isMC() || !config_->signalRegion() )
            {
               h1_[Form("m_jet%d%d_%s",j+1,k+1,label.c_str())]  -> Fill(c_ij.m(),weight_);
            }
            else  // blind
            {
               h1_[Form("m_jet%d%d_%s",j+1,k+1,label.c_str())]  -> Fill(0.,weight_);
            }
            if ( config_->isMC() && config_->histogramJetsPerFlavour() )
            {
               std::string flv1 = "udsg";
               std::string flv2 = "udsg";
               if ( ! config_ -> useJetsExtendedFlavour() )
               {
                  if ( abs(selectedJets_[j]->flavour()) == 4 ) flv1 = "c"; 
                  if ( abs(selectedJets_[j]->flavour()) == 5 ) flv1 = "b"; 
                  if ( abs(selectedJets_[k]->flavour()) == 4 ) flv2 = "c"; 
                  if ( abs(selectedJets_[k]->flavour()) == 5 ) flv2 = "b"; 
               }
               else
               {
                  flv1 = selectedJets_[j]->extendedFlavour();
                  flv2 = selectedJets_[k]->extendedFlavour();
               }
               h1_[Form("dptrel_jet%d%d_%s_%s_%s" , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> Fill(fabs(selectedJets_[j]->pt()-selectedJets_[k]->pt())/selectedJets_[j]->pt(),weight_);
               h1_[Form("dpt_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> Fill(fabs(selectedJets_[j]->pt()-selectedJets_[k]->pt()),weight_);
               h1_[Form("dr_jet%d%d_%s_%s_%s"     , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> Fill(c_ij.deltaR(),weight_);
               h1_[Form("deta_jet%d%d_%s_%s_%s"   , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> Fill(c_ij.deltaEta(),weight_);
               h1_[Form("dphi_jet%d%d_%s_%s_%s"   , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> Fill(fabs(selectedJets_[j]->deltaPhi(*selectedJets_[k])));
               h1_[Form("pt_jet%d%d_%s_%s_%s"     , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> Fill(c_ij.pt(),weight_);
               h1_[Form("eta_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> Fill(c_ij.eta(),weight_);
               h1_[Form("phi_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> Fill(c_ij.phi()*180./acos(-1.),weight_);
               h1_[Form("m_jet%d%d_%s_%s_%s"      , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> Fill(c_ij.m(),weight_);

            }
         }
      }
   }
   this->output()->cd();
   
   cutflow(Form("*** Filling jets histograms - %s",label.c_str()));
   
}

void JetAnalyser::fillJetHistograms ( const int &  r = -1, const std::string &  label = "x", const float &  sf = 1., const bool &  workflow = false  )

virtual

For experts only (fills pre-defined histograms for a specific jet, applying weights)

Definition at line 870 of file JetAnalyser.cc.

References btag(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, analysis::tools::BaseAnalyser::h2_, n_hjets_, analysis::tools::BaseAnalyser::output(), selectedJets_, and analysis::tools::BaseAnalyser::weight_.

{
   if ( r < 1 || r > n_hjets_ ) return;
   
   if ( workflow ) // BE CAREFUL with this
   {
      this->output()->cd();
      ++ cutflow_;
      if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("*** Filling jet # %d histograms - %s",r, label.c_str()));
      this->output()->cd(label.c_str());
   }
   
   int j = r-1;
   
   // 1D histograms
   h1_[Form("pt_jet%d_%s",j+1,label.c_str())]   -> Fill(selectedJets_[j]->pt(),weight_*sf);
   // barrel and endcap pt distributions
   float eta = selectedJets_[j]->eta();
   if ( config_-> histogramJetsRegionSplit() )
   {
      if ( eta < -1.4 )               h1_[Form("pt_jet%d_me_%s"  , j+1,label.c_str())] -> Fill(selectedJets_[j]->pt(),weight_*sf);
      if ( eta >  1.4 )               h1_[Form("pt_jet%d_pe_%s"  , j+1,label.c_str())] -> Fill(selectedJets_[j]->pt(),weight_*sf);
      if ( eta <  0.0 && eta > -1.0 ) h1_[Form("pt_jet%d_mb_%s"  , j+1,label.c_str())] -> Fill(selectedJets_[j]->pt(),weight_*sf);
      if ( eta >  0.0 && eta <  1.0 ) h1_[Form("pt_jet%d_pb_%s"  , j+1,label.c_str())] -> Fill(selectedJets_[j]->pt(),weight_*sf);
      if ( eta < -1.0 && eta > -1.4 ) h1_[Form("pt_jet%d_mbe_%s" , j+1,label.c_str())] -> Fill(selectedJets_[j]->pt(),weight_*sf);
      if ( eta >  1.0 && eta <  1.4 ) h1_[Form("pt_jet%d_pbe_%s" , j+1,label.c_str())] -> Fill(selectedJets_[j]->pt(),weight_*sf);
   }
      
   //
   h1_[Form("eta_jet%d_%s",j+1,label.c_str())]  -> Fill(selectedJets_[j]->eta(),weight_*sf);
   h1_[Form("phi_jet%d_%s",j+1,label.c_str())]  -> Fill(selectedJets_[j]->phi()*180./acos(-1.),weight_*sf);
   float mybtag = btag(*selectedJets_[j],config_->btagAlgorithm());
   float mybtaglog = 1.e-7;
   if ( mybtag > 0 ) mybtaglog = -log(1.-mybtag);
   h1_[Form("btag_jet%d_%s",j+1,label.c_str())]    -> Fill(mybtag,weight_*sf);
   h1_[Form("btaglog_jet%d_%s",j+1,label.c_str())] -> Fill(mybtaglog,weight_*sf);
   h1_[Form("qglikelihood_jet%d_%s", j+1,label.c_str())] -> Fill(selectedJets_[j]->qgLikelihood(),weight_*sf);
   h1_[Form("nconstituents_jet%d_%s", j+1,label.c_str())] -> Fill(selectedJets_[j]->constituents(),weight_*sf);
      
   if ( config_->btagAlgorithm() == "deepcsv")
   {
      h1_[Form("btag_light_jet%d_%s", j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_deeplight"),weight_*sf);
      h1_[Form("btag_c_jet%d_%s"    , j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_deepc"),weight_*sf); 
      h1_[Form("btag_b_jet%d_%s"    , j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_deepb"),weight_*sf); 
      h1_[Form("btag_bb_jet%d_%s"   , j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_deepbb"),weight_*sf); 
      h1_[Form("btag_cc_jet%d_%s"   , j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_deepcc"),weight_*sf); 
   }
   if ( config_->btagAlgorithm() == "deepflavour" || config_->btagAlgorithm() == "deepjet" )
   {
      h1_[Form("btag_light_jet%d_%s", j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dflight"),weight_*sf);
      h1_[Form("btag_g_jet%d_%s"    , j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dfg")    ,weight_*sf); 
      h1_[Form("btag_c_jet%d_%s"    , j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dfc")    ,weight_*sf); 
      h1_[Form("btag_b_jet%d_%s"    , j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dfb")    ,weight_*sf); 
      h1_[Form("btag_bb_jet%d_%s"   , j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dfbb")   ,weight_*sf); 
      h1_[Form("btag_lepb_jet%d_%s" , j+1,label.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dflepb") ,weight_*sf); 
   }
   // 2D histograms
   h2_[Form("pt_eta_jet%d_%s"  , j+1,label.c_str())] -> Fill(selectedJets_[j]->pt(), selectedJets_[j]->eta(), weight_*sf);
      
      
   if ( config_->isMC() && config_->histogramJetsPerFlavour() )
   {
      std::string flv = "udsg";
      if ( ! config_ -> useJetsExtendedFlavour() )
      {
         if ( abs(selectedJets_[j]->flavour()) == 4 ) flv = "c"; 
         if ( abs(selectedJets_[j]->flavour()) == 5 ) flv = "b"; 
      }
      else
      {
         flv = selectedJets_[j]->extendedFlavour();
      }
      // 1D histograms
      h1_[Form("pt_jet%d_%s_%s"  , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->pt(),weight_*sf);
      h1_[Form("eta_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->eta(),weight_*sf);
      h1_[Form("phi_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->phi()*180./acos(-1.),weight_*sf);
      h1_[Form("btag_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())]  -> Fill(btag(*selectedJets_[j],config_->btagAlgorithm()),weight_*sf);
      h1_[Form("qglikelihood_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())] -> Fill(selectedJets_[j]->qgLikelihood(),weight_*sf);
      h1_[Form("nconstituents_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())] -> Fill(selectedJets_[j]->constituents(),weight_*sf);
      if ( config_->btagAlgorithm() == "deepcsv")
      {
         h1_[Form("btag_light_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_deeplight"),weight_*sf);
         h1_[Form("btag_c_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_deepc"),weight_*sf); 
         h1_[Form("btag_b_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_deepb"),weight_*sf); 
         h1_[Form("btag_bb_jet%d_%s_%s"   , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_deepbb"),weight_*sf); 
         h1_[Form("btag_cc_jet%d_%s_%s"   , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_deepcc"),weight_*sf); 
      }
      if ( config_->btagAlgorithm() == "deepflavour" || config_->btagAlgorithm() == "deepjet")
      {
         h1_[Form("btag_light_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dflight"),weight_*sf);
         h1_[Form("btag_g_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dfg")    ,weight_*sf); 
         h1_[Form("btag_c_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dfc")    ,weight_*sf); 
         h1_[Form("btag_b_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dfb")    ,weight_*sf); 
         h1_[Form("btag_bb_jet%d_%s_%s"   , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dfbb")   ,weight_*sf); 
         h1_[Form("btag_lepb_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str())]  -> Fill(selectedJets_[j]->btag("btag_dflepb") ,weight_*sf); 
         
         
      }
      // 2D histograms
      h2_[Form("pt_eta_jet%d_%s_%s"  , j+1,label.c_str(),flv.c_str())] -> Fill(selectedJets_[j]->pt(), selectedJets_[j]->eta(), weight_*sf);
   }
      
   this->output()->cd();
   if ( workflow ) h1_["cutflow"] -> Fill(cutflow_,weight_*sf);
   
}

float JetAnalyser::getBtagSF ( const int &  r )

virtual

Definition at line 1032 of file JetAnalyser.cc.

References btagSF(), analysis::tools::BaseAnalyser::config_, and ReHLT_L1Trepack_HLTUser_DataNtuple::isMC.

{
   float sf = 1.;
   int j = r-1;
   if ( ! config_-> isMC() || config_->btagScaleFactors() == "" ) return sf;  // will not apply btag SF
   if ( ! config_->signalRegion() && r == config_->revBtagJet() ) return sf;
   
   if ( config_->jetsBtagWP()[j] != "xxx" )  sf = this->btagSF(r,config_->jetsBtagWP()[j]).nominal;
   
   return sf;
}

bool JetAnalyser::jetCorrections ( )

virtual

multiple actions: apply JER and b-tag regression corrections

Definition at line 1429 of file JetAnalyser.cc.

References actionApplyBjetRegression(), actionApplyJER(), and analysis::tools::BaseAnalyser::config().

Referenced by main().

{
// CORRECTIONS
   // b energy regression
      if ( this->config()->bRegression() )  this->actionApplyBjetRegression();
   // Jet energy resolution smearing
      this->actionApplyJER();
      
      return true;
      
}

void JetAnalyser::jetHistograms ( const int &  n, const std::string &  label = "x"  )

virtual

Creates pre-defined histograms in directory 'label' for analysis with 'n' jets.

Definition at line 172 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, flavours_, analysis::tools::BaseAnalyser::h1_, analysis::tools::BaseAnalyser::h2_, naf_mult_submit::int, n_hjets_, and analysis::tools::BaseAnalyser::output().

Referenced by jetHistograms(), main(), and analysis::mssmhbb::MssmHbbAnalyser::MssmHbbAnalyser().

{
   this->output()->cd();
   this->output()->mkdir(label.c_str());
   this->output()->cd(label.c_str());
   
   n_hjets_ = n;
   
   h1_[Form("jet_hist_weight_%s",label.c_str())] = std::make_shared<TH1F>(Form("jet_hist_weight_%s",label.c_str()) , Form("jet_hist_weight_%s",label.c_str()) ,1 , 0. , 1. );
   
   // btag variable binning
   float min1 = 0.0;
   float max1 = 0.1;
   float size1 = 0.0001;
   int nbins1 = int((max1-min1)/size1);
   float min2 = 0.1;
   float max2 = 0.9;
   float size2 = 0.001;
   int nbins2 = int((max2-min2)/size2);
   float min3 = 0.9;
   float max3 = 1.0;
   float size3 = 0.0001;
   int nbins3 = int((max3-min3)/size3);
   int nbins_btag = nbins1+nbins2+nbins3;
   
   std::vector<float> bins_btag;
   bins_btag.clear();
   int counter = 0; 
   for ( int i = 0; i<nbins1; ++i)   {  bins_btag.push_back(min1 + size1*i);  ++counter; }
   for ( int i = 0; i<nbins2; ++i)   {  bins_btag.push_back(min2 + size2*i);  ++counter; }
   for ( int i = 0; i<nbins3+1; ++i) {  bins_btag.push_back(min3 + size3*i);  ++counter; }
   
   // uniform binning for btag (comment it out if you want the variable binning above)
   float size = 0.0002;
   nbins_btag = int(1./size);
   bins_btag.clear();
   for ( int i = 0; i<nbins_btag+1; ++i)   {  bins_btag.push_back(size*i);  ++counter; }
   
   
   for ( int j = 0; j < n; ++j ) // loop over jets
   {
      // 1D histograms
      h1_[Form("pt_jet%d_%s"  , j+1,label.c_str())]  = std::make_shared<TH1F>(Form("pt_jet%d"  , j+1) , Form("pt_jet%d_%s"  , j+1,label.c_str()) ,1500 , 0   , 1500  );
      h1_[Form("eta_jet%d_%s" , j+1,label.c_str())]  = std::make_shared<TH1F>(Form("eta_jet%d" , j+1) , Form("eta_jet%d_%s" , j+1,label.c_str()) , 600 , -3, 3 );
      h1_[Form("phi_jet%d_%s" , j+1,label.c_str())]  = std::make_shared<TH1F>(Form("phi_jet%d" , j+1) , Form("phi_jet%d_%s" , j+1,label.c_str()) , 360 , -180, 180 );
      h1_[Form("btag_jet%d_%s", j+1,label.c_str())]  = std::make_shared<TH1F>(Form("btag_jet%d", j+1) , Form("btag_jet%d_%s", j+1,label.c_str()) , nbins_btag, &bins_btag[0] );
      h1_[Form("btaglog_jet%d_%s", j+1,label.c_str())]  = std::make_shared<TH1F>(Form("btaglog_jet%d", j+1) , Form("btaglog_jet%d_%s", j+1,label.c_str()) , 200 , 1.e-6, 10 );
      h1_[Form("qglikelihood_jet%d_%s", j+1,label.c_str())]  = std::make_shared<TH1F>(Form("qglikelihood_jet%d", j+1) , Form("qglikelihood_jet%d_%s", j+1,label.c_str()) , 200 , 0, 1 );
      h1_[Form("nconstituents_jet%d_%s", j+1,label.c_str())]  = std::make_shared<TH1F>(Form("nconstituents_jet%d", j+1) , Form("nconstituents_jet%d_%s", j+1,label.c_str()) , 200 , 0, 200 );
      // pt distributions separate for barrel and endcap, minus and plus sides and overlaps
      if ( config_-> histogramJetsRegionSplit() )
      {
         h1_[Form("pt_jet%d_me_%s"  , j+1,label.c_str())]  = std::make_shared<TH1F>(Form("pt_jet%d_me"  , j+1) , Form("pt_jet%d_me_%s"  , j+1,label.c_str()) ,1500 , 0   , 1500  );
         h1_[Form("pt_jet%d_pe_%s"  , j+1,label.c_str())]  = std::make_shared<TH1F>(Form("pt_jet%d_pe"  , j+1) , Form("pt_jet%d_pe_%s"  , j+1,label.c_str()) ,1500 , 0   , 1500  );
         h1_[Form("pt_jet%d_mb_%s"  , j+1,label.c_str())]  = std::make_shared<TH1F>(Form("pt_jet%d_mb"  , j+1) , Form("pt_jet%d_mb_%s"  , j+1,label.c_str()) ,1500 , 0   , 1500  );
         h1_[Form("pt_jet%d_pb_%s"  , j+1,label.c_str())]  = std::make_shared<TH1F>(Form("pt_jet%d_pb"  , j+1) , Form("pt_jet%d_pb_%s"  , j+1,label.c_str()) ,1500 , 0   , 1500  );
         h1_[Form("pt_jet%d_mbe_%s" , j+1,label.c_str())]  = std::make_shared<TH1F>(Form("pt_jet%d_mbe" , j+1) , Form("pt_jet%d_mbe_%s"  , j+1,label.c_str()) ,1500 , 0   , 1500  );
         h1_[Form("pt_jet%d_pbe_%s" , j+1,label.c_str())]  = std::make_shared<TH1F>(Form("pt_jet%d_pbe" , j+1) , Form("pt_jet%d_pbe_%s"  , j+1,label.c_str()) ,1500 , 0   , 1500  );
      }
      
      h1_[Form("pt_jet%d_%s"  , j+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d p_{T} [GeV]",j+1));
      h1_[Form("eta_jet%d_%s" , j+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d  #eta",j+1));
      h1_[Form("phi_jet%d_%s" , j+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d  #phi",j+1));
      h1_[Form("btag_jet%d_%s", j+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d btag discriminator",j+1)); 
      h1_[Form("btaglog_jet%d_%s", j+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d -ln(1-btag discriminator)",j+1)); 
      h1_[Form("qglikelihood_jet%d_%s", j+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d q-g likelihood",j+1)); 
      h1_[Form("nconstituents_jet%d_%s", j+1,label.c_str())]-> GetXaxis() -> SetTitle(Form("Jet %d n constituents",j+1)); 
      
      if ( config_->btagAlgorithm() == "deepcsv")
      {
         h1_[Form("btag_light_jet%d_%s", j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_light_jet%d", j+1) , Form("btag_light_jet%d_%s", j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
         h1_[Form("btag_c_jet%d_%s"    , j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_c_jet%d"    , j+1) , Form("btag_c_jet%d_%s"    , j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
         h1_[Form("btag_b_jet%d_%s"    , j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_b_jet%d"    , j+1) , Form("btag_b_jet%d_%s"    , j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
         h1_[Form("btag_bb_jet%d_%s"   , j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_bb_jet%d"   , j+1) , Form("btag_bb_jet%d_%s"   , j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
         h1_[Form("btag_cc_jet%d_%s"   , j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_cc_jet%d"   , j+1) , Form("btag_cc_jet%d_%s"   , j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
      }
      if ( config_->btagAlgorithm() == "deepflavour" || config_->btagAlgorithm() == "deepjet" )
      {
         h1_[Form("btag_light_jet%d_%s", j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_light_jet%d", j+1) , Form("btag_light_jet%d_%s", j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
         h1_[Form("btag_g_jet%d_%s"    , j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_g_jet%d"    , j+1) , Form("btag_g_jet%d_%s"    , j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
         h1_[Form("btag_c_jet%d_%s"    , j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_c_jet%d"    , j+1) , Form("btag_c_jet%d_%s"    , j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
         h1_[Form("btag_b_jet%d_%s"    , j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_b_jet%d"    , j+1) , Form("btag_b_jet%d_%s"    , j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
         h1_[Form("btag_bb_jet%d_%s"   , j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_bb_jet%d"   , j+1) , Form("btag_bb_jet%d_%s"   , j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
         h1_[Form("btag_lepb_jet%d_%s" , j+1,label.c_str())] = std::make_shared<TH1F>(Form("btag_lepb_jet%d" , j+1) , Form("btag_lepb_jet%d_%s" , j+1,label.c_str()) , nbins_btag, &bins_btag[0] );         
      }
      
      // 2D histograms
      h2_[Form("pt_eta_jet%d_%s"  , j+1,label.c_str())]  = std::make_shared<TH2F>(Form("pt_eta_jet%d"  , j+1) , Form("pt_eta_jet%d_%s"  , j+1,label.c_str()) ,1500 , 0   , 1500, 600, -3, 3  );
      h2_[Form("pt_eta_jet%d_%s"  , j+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d p_{T} [GeV]",j+1));
      h2_[Form("pt_eta_jet%d_%s"  , j+1,label.c_str())] -> GetYaxis() -> SetTitle(Form("Jet %d  #eta",j+1));
      
      if ( config_->isMC() && config_->histogramJetsPerFlavour() )
      {
         for ( auto & flv : flavours_ ) // flavour dependent histograms
         {
            // 1D histograms
            h1_[Form("pt_jet%d_%s_%s"  , j+1,label.c_str(),flv.c_str())]  = std::make_shared<TH1F>(Form("pt_jet%d_%s"   , j+1, flv.c_str()) , Form("pt_jet%d_%s_%s"  , j+1,label.c_str(),flv.c_str()) ,1500 , 0   , 1500  );
            h1_[Form("eta_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str())]  = std::make_shared<TH1F>(Form("eta_jet%d_%s"  , j+1, flv.c_str()) , Form("eta_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str()) , 600 , -3, 3 );
            h1_[Form("phi_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str())]  = std::make_shared<TH1F>(Form("phi_jet%d_%s"  , j+1, flv.c_str()) , Form("phi_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str()) , 360 , -180, 180 );
            h1_[Form("btag_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())]  = std::make_shared<TH1F>(Form("btag_jet%d_%s" , j+1, flv.c_str()) , Form("btag_jet%d_%s_%s", j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );
            h1_[Form("qglikelihood_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())]  = std::make_shared<TH1F>(Form("qglikelihood_jet%d_%s" , j+1, flv.c_str()) , Form("qglikelihood_jet%d_%s_%s", j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );
            h1_[Form("nconstituents_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())]  = std::make_shared<TH1F>(Form("nconstituents_jet%d_%s" , j+1, flv.c_str()) , Form("nconstituents_jet%d_%s_%s", j+1,label.c_str(),flv.c_str()) , 200 , 0, 200 );
            
            h1_[Form("pt_jet%d_%s_%s"  , j+1,label.c_str(),flv.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d (%s) p_{T} [GeV]"        , j+1, flv.c_str()));
            h1_[Form("eta_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d (%s)  #eta"              , j+1, flv.c_str()));
            h1_[Form("phi_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d (%s)  #phi"              , j+1, flv.c_str()));
            h1_[Form("btag_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d (%s) btag discriminator" , j+1, flv.c_str())); 
            h1_[Form("qglikelihood_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d (%s) q-g likelihood" , j+1, flv.c_str())); 
            h1_[Form("nconstituents_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d (%s) n constituents" , j+1, flv.c_str())); 
            
            if ( config_->btagAlgorithm() == "deepcsv")
            {
               h1_[Form("btag_light_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_light_jet%d_%s", j+1,flv.c_str()) , Form("btag_light_jet%d_%s_%s", j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
               h1_[Form("btag_c_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_c_jet%d_%s"    , j+1,flv.c_str()) , Form("btag_c_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
               h1_[Form("btag_b_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_b_jet%d_%s"    , j+1,flv.c_str()) , Form("btag_b_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
               h1_[Form("btag_bb_jet%d_%s_%s"   , j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_bb_jet%d_%s"   , j+1,flv.c_str()) , Form("btag_bb_jet%d_%s_%s"   , j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
               h1_[Form("btag_cc_jet%d_%s_%s"   , j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_cc_jet%d_%s"   , j+1,flv.c_str()) , Form("btag_cc_jet%d_%s_%s"   , j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
            }
            if ( config_->btagAlgorithm() == "deepflavour" || config_->btagAlgorithm() == "deepjet" )
            {
               h1_[Form("btag_light_jet%d_%s_%s", j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_light_jet%d_%s", j+1,flv.c_str()) , Form("btag_light_jet%d_%s_%s", j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
               h1_[Form("btag_g_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_g_jet%d_%s"    , j+1,flv.c_str()) , Form("btag_g_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
               h1_[Form("btag_c_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_c_jet%d_%s"    , j+1,flv.c_str()) , Form("btag_c_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
               h1_[Form("btag_b_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_b_jet%d_%s"    , j+1,flv.c_str()) , Form("btag_b_jet%d_%s_%s"    , j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
               h1_[Form("btag_bb_jet%d_%s_%s"   , j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_bb_jet%d_%s"   , j+1,flv.c_str()) , Form("btag_bb_jet%d_%s_%s"   , j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
               h1_[Form("btag_lepb_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str())] = std::make_shared<TH1F>(Form("btag_lepb_jet%d_%s" , j+1,flv.c_str()) , Form("btag_lepb_jet%d_%s_%s" , j+1,label.c_str(),flv.c_str()) , nbins_btag, &bins_btag[0] );         
            }
            // 2D histograms
            h2_[Form("pt_eta_jet%d_%s_%s"  , j+1,label.c_str(),flv.c_str())]  = std::make_shared<TH2F>(Form("pt_eta_jet%d_%s" , j+1,flv.c_str()) , Form("pt_eta_jet%d_%s_%s"  , j+1,label.c_str(),flv.c_str()) ,1500 , 0   , 1500, 600, -3, 3  );
            h2_[Form("pt_eta_jet%d_%s_%s"  , j+1,label.c_str(),flv.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d (%s) p_{T} [GeV]" , j+1,flv.c_str()));
            h2_[Form("pt_eta_jet%d_%s_%s"  , j+1,label.c_str(),flv.c_str())] -> GetYaxis() -> SetTitle(Form("Jet %d (%s)  #eta" ,j+1,flv.c_str()));

         }
      }
      if ( config_->doDijet()  )  // dijet histograms
      {
         for ( int k = j+1; k < n && j < n; ++k )
         {
            h1_[Form("dptrel_jet%d%d_%s" , j+1,k+1,label.c_str())]  = std::make_shared<TH1F>(Form("dptrel_jet%d%d" , j+1,k+1) , Form("dptrel_jet%d%d_%s" , j+1,k+1,label.c_str()) ,1000 , 0,1 );
            h1_[Form("dpt_jet%d%d_%s"    , j+1,k+1,label.c_str())]  = std::make_shared<TH1F>(Form("dpt_jet%d%d" , j+1,k+1)    , Form("dpt_jet%d%d_%s"    , j+1,k+1,label.c_str()) ,1000 , 0,1000 );
            h1_[Form("dr_jet%d%d_%s"     , j+1,k+1,label.c_str())]  = std::make_shared<TH1F>(Form("dr_jet%d%d"  , j+1,k+1)    , Form("dr_jet%d%d_%s"     , j+1,k+1,label.c_str()) , 100 , 0, 5 );
            h1_[Form("deta_jet%d%d_%s"   , j+1,k+1,label.c_str())]  = std::make_shared<TH1F>(Form("deta_jet%d%d", j+1,k+1)    , Form("deta_jet%d%d_%s"   , j+1,k+1,label.c_str()) , 100 , 0,10 );
            h1_[Form("dphi_jet%d%d_%s"   , j+1,k+1,label.c_str())]  = std::make_shared<TH1F>(Form("dphi_jet%d%d", j+1,k+1)    , Form("dphi_jet%d%d_%s"   , j+1,k+1,label.c_str()) , 315 , 0, 3.15 );
            h1_[Form("pt_jet%d%d_%s"     , j+1,k+1,label.c_str())]  = std::make_shared<TH1F>(Form("pt_jet%d%d"  , j+1,k+1)    , Form("pt_jet%d%d_%s"     , j+1,k+1,label.c_str()) , 300 , 0,3000 );
            h1_[Form("eta_jet%d%d_%s"    , j+1,k+1,label.c_str())]  = std::make_shared<TH1F>(Form("eta_jet%d%d" , j+1,k+1)    , Form("eta_jet%d%d_%s"    , j+1,k+1,label.c_str()) , 200 , -10,10 );
            h1_[Form("phi_jet%d%d_%s"    , j+1,k+1,label.c_str())]  = std::make_shared<TH1F>(Form("phi_jet%d%d" , j+1,k+1)    , Form("phi_jet%d%d_%s"    , j+1,k+1,label.c_str()) , 360 , -180,180 );
            h1_[Form("m_jet%d%d_%s"      , j+1,k+1,label.c_str())]  = std::make_shared<TH1F>(Form("m_jet%d%d"   , j+1,k+1)    , Form("m_jet%d%d_%s"      , j+1,k+1,label.c_str()) ,3000 , 0,3000 );
            
            h1_[Form("dptrel_jet%d%d_%s" , j+1,k+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("#DeltaP_{T}(Jet %d, Jet %d)/Jet %d p_{T}",j+1,k+1,j+1));
            h1_[Form("dpt_jet%d%d_%s"    , j+1,k+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("#Delta p_{T}(Jet %d, Jet %d) [GeV]",j+1,k+1));
            h1_[Form("dr_jet%d%d_%s"     , j+1,k+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("#DeltaR(Jet %d, Jet %d)",j+1,k+1));
            h1_[Form("deta_jet%d%d_%s"   , j+1,k+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("#Delta#eta(Jet %d, Jet %d)",j+1,k+1));
            h1_[Form("dphi_jet%d%d_%s"   , j+1,k+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("#Delta#phi(Jet %d, Jet %d)",j+1,k+1));
            h1_[Form("pt_jet%d%d_%s"     , j+1,k+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d + Jet %d p_{T} [GeV]",j+1,k+1));
            h1_[Form("eta_jet%d%d_%s"    , j+1,k+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d + Jet %d  #eta",j+1,k+1));
            h1_[Form("phi_jet%d%d_%s"    , j+1,k+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("Jet %d + Jet %d  #phi",j+1,k+1));
            h1_[Form("m_jet%d%d_%s"      , j+1,k+1,label.c_str())] -> GetXaxis() -> SetTitle(Form("M_{%d%d} [GeV]",j+1,k+1));
            
            if ( config_->isMC() && config_->histogramJetsPerFlavour() )
            {
               for ( auto & flv1 : flavours_ )
               {
                  for ( auto & flv2 : flavours_ )
                  {
                     h1_[Form("dptrel_jet%d%d_%s_%s_%s" , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  = std::make_shared<TH1F>(Form("dptrel_jet%d%d_%s_%s" , j+1,k+1,flv1.c_str(),flv2.c_str()) , Form("dptrel_jet%d%d_%s_%s_%s" , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str()) , 100 , 0,1 );
                     h1_[Form("dpt_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  = std::make_shared<TH1F>(Form("dpt_jet%d%d_%s_%s"    , j+1,k+1,flv1.c_str(),flv2.c_str()) , Form("dpt_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str()) , 100 , 0,1000 );
                     h1_[Form("dr_jet%d%d_%s_%s_%s"     , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  = std::make_shared<TH1F>(Form("dr_jet%d%d_%s_%s"     , j+1,k+1,flv1.c_str(),flv2.c_str()) , Form("dr_jet%d%d_%s_%s_%s"     , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str()) , 100 , 0, 5 );
                     h1_[Form("deta_jet%d%d_%s_%s_%s"   , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  = std::make_shared<TH1F>(Form("deta_jet%d%d_%s_%s"   , j+1,k+1,flv1.c_str(),flv2.c_str()) , Form("deta_jet%d%d_%s_%s_%s"   , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str()) , 100 , 0,10 );
                     h1_[Form("dphi_jet%d%d_%s_%s_%s"   , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  = std::make_shared<TH1F>(Form("dphi_jet%d%d_%s_%s"   , j+1,k+1,flv1.c_str(),flv2.c_str()) , Form("dphi_jet%d%d_%s_%s_%s"   , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str()) , 315 , 0, 3.15 );
                     h1_[Form("pt_jet%d%d_%s_%s_%s"     , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  = std::make_shared<TH1F>(Form("pt_jet%d%d_%s_%s"     , j+1,k+1,flv1.c_str(),flv2.c_str()) , Form("pt_jet%d%d_%s_%s_%s"     , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str()) , 300 , 0,300 );
                     h1_[Form("eta_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  = std::make_shared<TH1F>(Form("eta_jet%d%d_%s_%s"    , j+1,k+1,flv1.c_str(),flv2.c_str()) , Form("eta_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str()) , 200 , -10,10 );
                     h1_[Form("phi_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  = std::make_shared<TH1F>(Form("phi_jet%d%d_%s_%s"    , j+1,k+1,flv1.c_str(),flv2.c_str()) , Form("phi_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str()) , 360 , -180,180 );
                     h1_[Form("m_jet%d%d_%s_%s_%s"      , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  = std::make_shared<TH1F>(Form("m_jet%d%d_%s_%s"      , j+1,k+1,flv1.c_str(),flv2.c_str()) , Form("m_jet%d%d_%s_%s_%s"      , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str()) ,3000 , 0,300 );
                     
                     h1_[Form("dptrel_jet%d%d_%s_%s_%s" , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> GetXaxis() -> SetTitle(Form("#DeltaP_{T}(Jet %d (%s), Jet %d (%s))/Jet %d p_{T}",j+1,flv1.c_str(),k+1,flv2.c_str(),j+1));
                     h1_[Form("dpt_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> GetXaxis() -> SetTitle(Form("#Delta p_{T}(Jet %d (%s), Jet %d (%s)) [GeV]",j+1,flv1.c_str(),k+1,flv2.c_str()));
                     h1_[Form("dr_jet%d%d_%s_%s_%s"     , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> GetXaxis() -> SetTitle(Form("#DeltaR(Jet %d (%s), Jet %d (%s))",j+1,flv1.c_str(),k+1,flv2.c_str()));
                     h1_[Form("deta_jet%d%d_%s_%s_%s"   , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> GetXaxis() -> SetTitle(Form("#Delta#eta(Jet %d (%s), Jet %d (%s))",j+1,flv1.c_str(),k+1,flv2.c_str()));
                     h1_[Form("dphi_jet%d%d_%s_%s_%s"   , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> GetXaxis() -> SetTitle(Form("#Delta#phi(Jet %d (%s), Jet %d (%s))",j+1,flv1.c_str(),k+1,flv2.c_str()));
                     h1_[Form("pt_jet%d%d_%s_%s_%s"     , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> GetXaxis() -> SetTitle(Form("Jet %d (%s) + Jet %d (%s) p_{T} [GeV]",j+1,flv1.c_str(),k+1,flv2.c_str()));
                     h1_[Form("eta_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> GetXaxis() -> SetTitle(Form("Jet %d (%s) + Jet %d (%s)  #eta",j+1,flv1.c_str(),k+1,flv2.c_str()));
                     h1_[Form("phi_jet%d%d_%s_%s_%s"    , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> GetXaxis() -> SetTitle(Form("Jet %d (%s) + Jet %d (%s)  #phi",j+1,flv1.c_str(),k+1,flv2.c_str()));
                     h1_[Form("m_jet%d%d_%s_%s_%s"      , j+1,k+1,label.c_str(),flv1.c_str(),flv2.c_str())]  -> GetXaxis() -> SetTitle(Form("M_{%d%d} (%s)(%s) [GeV]",j+1,k+1,flv1.c_str(),flv2.c_str()));                                          
                  }
               }
            }
         }
      }
   }
   
   this->output()->cd();
}

void JetAnalyser::jetHistograms ( const std::string &  label = "x" )

virtual

Creates pre-defined histograms in directory 'label' for analysis with config->nJetsMin()

Definition at line 168 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, and jetHistograms().

{
   this -> jetHistograms(config_->nJetsMin(), label);
}

std::vector< std::shared_ptr< Jet > > JetAnalyser::jets

(

)

vector of pointers of all jets from the "Jets" collection

Definition at line 613 of file JetAnalyser.cc.

References jets_.

Referenced by analysisWithJets().

{
   return jets_;
}

void JetAnalyser::jets

(

const std::string &

col

)

Creates the "Jets" collection.

Definition at line 163 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::analysis_.

{
   analysis_->addTree<Jet> ("Jets",col);
}

void JetAnalyser::jetSwap ( const int &  r1, const int &  r2  )

virtual

Definition at line 1104 of file JetAnalyser.cc.

References selectedJets_.

Referenced by analysis::mssmhbb::MssmHbbAnalyser::muonJet().

{
   if ( r1 == r2 ) return;
   int j1 = r1-1;
   int j2 = r2-1;
//    ++ cutflow_;
//    if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
//       h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d <--> Jet %d: jets ranking was swapped ",r1,r2));
   
   auto jet1 = selectedJets_[j1];
   auto jet2 = selectedJets_[j2];
   
   selectedJets_[j1] = jet2;
   selectedJets_[j2] = jet1;
   
//    h1_["cutflow"] -> Fill(cutflow_,weight_);
   
}

bool JetAnalyser::onlineBJetMatching ( const int &  r )

virtual

Definition at line 778 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), and selectedJets_.

Referenced by analysis::mssmhbb::MssmHbbAnalyser::btagSelection(), and main().

{
   if ( config_->triggerObjectsBJets() == "" ) return true;
   
   bool isgood = true;
   std::string label = Form("Jet %d: online b jet match (deltaR < %4.3f)",r,config_-> triggerMatchCaloBJetsDrMax());
   
   int j = r-1;
   
   std::shared_ptr<Jet> jet = selectedJets_[j];
   isgood = ( jet->matched(config_->triggerObjectsBJets()) );
   
   cutflow(label,isgood);
   
   return isgood;
}

bool JetAnalyser::onlineJetMatching ( const int &  r )

virtual

Definition at line 746 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), and selectedJets_.

Referenced by analysis::mssmhbb::MssmHbbAnalyser::jetSelection(), and main().

{
   if ( config_->triggerObjectsL1Jets() == "" && config_->triggerObjectsCaloJets() == "" && config_->triggerObjectsPFJets() == "") return true;
   if ( config_->nJetsMin() < 0 ) return true;
   
   bool isgood = true;
   std::string label = Form("Jet %d: online jet match (deltaR: L1 < %4.3f, Calo < %4.3f, PF < %4.3f)",r,config_-> triggerMatchL1JetsDrMax(),config_-> triggerMatchCaloJetsDrMax(),config_-> triggerMatchPFJetsDrMax());
   
   int j = r-1;
   
   std::string triggerObjectsL1Jets = config_->triggerObjectsL1Jets();
   if ( config_->triggerEmulateL1Jets() != "" &&  config_->triggerEmulateL1JetsNMin() > 0 )
      triggerObjectsL1Jets = config_->triggerEmulateL1Jets();
   std::string triggerObjectsCaloJets = config_->triggerObjectsCaloJets();
   if ( config_->triggerEmulateCaloJets() != "" &&  config_->triggerEmulateCaloJetsNMin() > 0 )
      triggerObjectsCaloJets = config_->triggerEmulateCaloJets();
   std::string triggerObjectsPFJets = config_->triggerObjectsPFJets();
   if ( config_->triggerEmulatePFJets() != "" &&  config_->triggerEmulatePFJetsNMin() > 0 )
      triggerObjectsPFJets = config_->triggerEmulatePFJets();
   
   
   std::shared_ptr<Jet> jet = selectedJets_[j];
   isgood = (           jet->matched(triggerObjectsL1Jets   ) );
   isgood = ( isgood && jet->matched(triggerObjectsCaloJets ) );
   isgood = ( isgood && jet->matched(triggerObjectsPFJets   ) );

   cutflow(label,isgood);
   
   return isgood;
}

std::vector< std::shared_ptr< Jet > > JetAnalyser::selectedJets

(

)

vector of pointers of the selectedJets

Definition at line 618 of file JetAnalyser.cc.

References selectedJets_.

{
   return selectedJets_;
}

bool JetAnalyser::selectionBJet ( const int &  r )

virtual

Definition at line 708 of file JetAnalyser.cc.

References btag(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), naf_mult_submit::int, selectedJets_, and selectionNonBJet().

Referenced by analysis::mssmhbb::MssmHbbAnalyser::btagSelection(), and main().

{
   if ( config_->nJetsMin() < config_->nBJetsMin() || config_->nBJetsMin() < 1 || r > config_->nBJetsMin() ||  (int)(config_->jetsBtagWP()).size() < config_->nBJetsMin() ) return true;
   
   if ( ! config_->signalRegion() && r == config_->revBtagJet() ) return this->selectionNonBJet(r);

   int j = r-1;
   if ( config_->btagWP(config_->jetsBtagWP()[j]) < 0 ) return true; // there is no selection here, so will not update the cutflow
         
   bool isgood = true;
   std::string label = Form("Jet %d: %s btag > %6.4f (%s)",r,config_->btagAlgorithm().c_str(),config_->btagWP(config_->jetsBtagWP()[j]),config_->jetsBtagWP()[j].c_str());
   
   isgood = ( btag(*selectedJets_[j],config_->btagAlgorithm()) > config_->btagWP(config_->jetsBtagWP()[j]) );
   
   cutflow(label,isgood);
   
   return isgood;
}

bool JetAnalyser::selectionBJetProbB ( const int &  r )

virtual

Definition at line 1172 of file JetAnalyser.cc.

References btag(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, selectedJets_, and analysis::tools::BaseAnalyser::weight_.

{
   if ( config_->jetsBtagProbB().size() == 0 ) return true;
   int j = r-1;
   float wp = config_->jetsBtagProbB()[j];
   std::string algo = config_->btagAlgorithm();
   if ( fabs(wp) > 1 ) return true; // there is no selection here, so will not update the cutflow
   
   ++ cutflow_;
   if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
   {
      if ( wp > 0 )
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob b < %6.4f",r,algo.c_str(),fabs(wp)));
      else
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob b > %6.4f",r,algo.c_str(),fabs(wp)));
   }
         
   
   if ( r > config_->nBJetsMin() ) 
   {
      std::cout << "* warning * -  JetAnalyser::selectionBJetProbB(): given jet rank > nbjetsmin. Returning false! " << std::endl;
      return false;
   }
   
   // jet  btag
   if ( algo == "deepcsv" )
   {
      if ( wp > 0 && selectedJets_[j]->btag("btag_deepb") > fabs(wp) ) return false;
      if ( wp < 0 && selectedJets_[j]->btag("btag_deepb") < fabs(wp) ) return false;
   }
   if ( algo == "deepflavour" || algo == "deepjet" )
   {
      if ( wp > 0 && selectedJets_[j]->btag("btag_dfb") > fabs(wp) ) return false;
      if ( wp < 0 && selectedJets_[j]->btag("btag_dfb") < fabs(wp) ) return false;
   }
   
   h1_["cutflow"] -> Fill(cutflow_,weight_);
   
   return true;
}

bool JetAnalyser::selectionBJetProbBB ( const int &  r )

virtual

Definition at line 1216 of file JetAnalyser.cc.

References btag(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, selectedJets_, and analysis::tools::BaseAnalyser::weight_.

{
   if ( config_->jetsBtagProbBB().size() == 0 ) return true;
   int j = r-1;
   float wp = config_->jetsBtagProbBB()[j];
   std::string algo = config_->btagAlgorithm();
   if ( fabs(wp) > 1 ) return true; // there is no selection here, so will not update the cutflow
   
   ++ cutflow_;
   if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
   {
      if ( wp > 0 )
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob bb < %6.4f",r,algo.c_str(),fabs(wp)));
      else
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob bb > %6.4f",r,algo.c_str(),fabs(wp)));
   }
         
   
   if ( r > config_->nBJetsMin() ) 
   {
      std::cout << "* warning * -  JetAnalyser::selectionBJetProbBB(): given jet rank > nbjetsmin. Returning false! " << std::endl;
      return false;
   }
   
   // jet  btag
   if ( algo == "deepcsv" )
   {
      if ( wp > 0 && selectedJets_[j]->btag("btag_deepbb") > fabs(wp) ) return false;
      if ( wp < 0 && selectedJets_[j]->btag("btag_deepbb") < fabs(wp) ) return false;
   }
   if ( algo == "deepflavour" || algo == "deepjet" )
   {
      if ( wp > 0 && selectedJets_[j]->btag("btag_dfbb") > fabs(wp) ) return false;
      if ( wp < 0 && selectedJets_[j]->btag("btag_dfbb") < fabs(wp) ) return false;
   }
   
   h1_["cutflow"] -> Fill(cutflow_,weight_);
   
   return true;
}

bool JetAnalyser::selectionBJetProbC ( const int &  r )

virtual

Definition at line 1295 of file JetAnalyser.cc.

References btag(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, selectedJets_, and analysis::tools::BaseAnalyser::weight_.

{
   if ( config_->jetsBtagProbC().size() == 0 ) return true;
   int j = r-1;
   float wp = config_->jetsBtagProbC()[j];
   std::string algo = config_->btagAlgorithm();
   if ( fabs(wp) > 1 ) return true; // there is no selection here, so will not update the cutflow
   
   ++ cutflow_;
   if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
   {
      if ( wp > 0 )
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob c < %6.4f",r,algo.c_str(),fabs(wp)));
      else
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob c > %6.4f",r,algo.c_str(),fabs(wp)));
   }
         
   
   if ( r > config_->nBJetsMin() ) 
   {
      std::cout << "* warning * -  JetAnalyser::selectionBJetProbC(): given jet rank > nbjetsmin. Returning false! " << std::endl;
      return false;
   }
   
   // jet  btag
   if ( algo == "deepcsv" )
   {
      if ( wp > 0 && selectedJets_[j]->btag("btag_deepc") > fabs(wp) ) return false;
      if ( wp < 0 && selectedJets_[j]->btag("btag_deepc") < fabs(wp) ) return false;
   }
   if ( algo == "deepflavour" || algo == "deepjet" )
   {
      if ( wp > 0 && selectedJets_[j]->btag("btag_dfc") > fabs(wp) ) return false;
      if ( wp < 0 && selectedJets_[j]->btag("btag_dfc") < fabs(wp) ) return false;
   }
   
   h1_["cutflow"] -> Fill(cutflow_,weight_);
   
   return true;
}

bool JetAnalyser::selectionBJetProbG ( const int &  r )

virtual

Definition at line 1336 of file JetAnalyser.cc.

References btag(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, selectedJets_, and analysis::tools::BaseAnalyser::weight_.

{
   if ( config_->jetsBtagProbG().size() == 0 ) return true;
   int j = r-1;
   float wp = config_->jetsBtagProbG()[j];
   std::string algo = config_->btagAlgorithm();
   if ( fabs(wp) > 1 || algo == "deepcsv" ) return true; // there is no selection here, so will not update the cutflow
   
   ++ cutflow_;
   if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
   {
      if ( wp > 0 )
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob g < %6.4f",r,algo.c_str(),fabs(wp)));
      else
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob g > %6.4f",r,algo.c_str(),fabs(wp)));
   }
         
   
   if ( r > config_->nBJetsMin() ) 
   {
      std::cout << "* warning * -  JetAnalyser::selectionBJetProbG(): given jet rank > nbjetsmin. Returning false! " << std::endl;
      return false;
   }
   
   // jet  btag
   if ( algo == "deepflavour" || algo == "deepjet" )
   {
      if ( wp > 0 && selectedJets_[j]->btag("btag_dfg") > fabs(wp) ) return false;
      if ( wp < 0 && selectedJets_[j]->btag("btag_dfg") < fabs(wp) ) return false;
   }
   
   h1_["cutflow"] -> Fill(cutflow_,weight_);
   
   return true;
}

bool JetAnalyser::selectionBJetProbLepB ( const int &  r )

virtual

Definition at line 1257 of file JetAnalyser.cc.

References btag(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, selectedJets_, and analysis::tools::BaseAnalyser::weight_.

{
   if ( config_->jetsBtagProbLepB().size() == 0 ) return true;
   int j = r-1;
   float wp = config_->jetsBtagProbLepB()[j];
   std::string algo = config_->btagAlgorithm();
   if ( fabs(wp) > 1 || algo == "deepcsv" ) return true; // there is no selection here, so will not update the cutflow
   
   ++ cutflow_;
   if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
   {
      if ( wp > 0 )
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob lepb < %6.4f",r,algo.c_str(),fabs(wp)));
      else
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob lepb > %6.4f",r,algo.c_str(),fabs(wp)));
   }
         
   
   if ( r > config_->nBJetsMin() ) 
   {
      std::cout << "* warning * -  JetAnalyser::selectionBJetProbLepB(): given jet rank > nbjetsmin. Returning false! " << std::endl;
      return false;
   }
   
   // jet  btag
   if ( algo == "deepflavour" || algo == "deepjet" )
   {
      if ( wp > 0 && selectedJets_[j]->btag("btag_dflepb") > fabs(wp) ) return false;
      if ( wp < 0 && selectedJets_[j]->btag("btag_dflepb") < fabs(wp) ) return false;
   }
   
   h1_["cutflow"] -> Fill(cutflow_,weight_);
   
   return true;
}

bool JetAnalyser::selectionBJetProbLight ( const int &  r )

virtual

Definition at line 1372 of file JetAnalyser.cc.

References btag(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, selectedJets_, and analysis::tools::BaseAnalyser::weight_.

{
   if ( config_->jetsBtagProbLight().size() == 0 ) return true;
   int j = r-1;
   float wp = config_->jetsBtagProbLight()[j];
   std::string algo = config_->btagAlgorithm();
   if ( fabs(wp) > 1 ) return true; // there is no selection here, so will not update the cutflow
   
   ++ cutflow_;
   if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" ) 
   {
      if ( wp > 0 )
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob light < %6.4f",r,algo.c_str(),fabs(wp)));
      else
         h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,Form("Jet %d: %s btag prob light > %6.4f",r,algo.c_str(),fabs(wp)));
   }
         
   
   if ( r > config_->nBJetsMin() ) 
   {
      std::cout << "* warning * -  JetAnalyser::selectionBJetProbLight(): given jet rank > nbjetsmin. Returning false! " << std::endl;
      return false;
   }
   
   // jet  btag
   if ( algo == "deepcsv" )
   {
      if ( wp > 0 && selectedJets_[j]->btag("btag_deeplight") > fabs(wp) ) return false;
      if ( wp < 0 && selectedJets_[j]->btag("btag_deeplight") < fabs(wp) ) return false;
   }
   if ( algo == "deepflavour" || algo == "deepjet" )
   {
      if ( wp > 0 && selectedJets_[j]->btag("btag_dflight") > fabs(wp) ) return false;
      if ( wp < 0 && selectedJets_[j]->btag("btag_dflight") < fabs(wp) ) return false;
   }
   
   h1_["cutflow"] -> Fill(cutflow_,weight_);
   
   return true;
}

bool JetAnalyser::selectionDiJetMass ( const int &  r1, const int &  r2  )

virtual

Definition at line 1066 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow_, analysis::tools::BaseAnalyser::h1_, analysis::tools::Candidate::m(), selectedJets_, and analysis::tools::BaseAnalyser::weight_.

{
   float min = config_->massMin();
   float max = config_->massMax();
   if ( min < 0. && max < 0. ) return true;
   
   ++cutflow_;
   if ( std::string(h1_["cutflow"] -> GetXaxis()-> GetBinLabel(cutflow_+1)) == "" )
   {
      std::string label = Form("M(Jet %d + Jet %d)",r1,r2);
      if ( min > 0. )
      {   
         if ( max > 0. && max > min ) label = Form("%5.1f GeV < %s < %5.1f GeV",min,label.c_str(),max);
         else                         label = Form("%s > %5.1f GeV",label.c_str(),min);
      }
      else
      {
         label = Form("%s < %5.1f GeV",label.c_str(),max);
      }
      
      h1_["cutflow"] -> GetXaxis()-> SetBinLabel(cutflow_+1,label.c_str());
   }
   
   int j1 = r1-1;
   int j2 = r2-1;
   Composite<Jet,Jet> c_j1j2(*(selectedJets_[j1]),*(selectedJets_[j2]));
   float mass = c_j1j2.m();
   
   if ( min > 0. && mass < min ) return false;
   if ( max > 0. && mass > max ) return false;
   
   h1_["cutflow"] -> Fill(cutflow_,weight_);
   
   return true;
}

bool JetAnalyser::selectionJet ( const int &  r )

virtual

Given the ranking 'r' of a jet, it returns whether the jet passes the pt_min and |eta_max|, optionally pt_max, where the values of the thresholds pt_min and |eta_max|, pt_max are passed by configuration file.

Definition at line 378 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_.

Referenced by analysis::mssmhbb::MssmHbbAnalyser::jetSelection(), and main().

{
   if ( r > config_->nJetsMin() ) return true;
   int j = r-1;
   
   float pt_min = config_->jetsPtMin()[j];
   float pt_max = -1.;
   float eta_max = config_->jetsEtaMax()[j];
   
   if ( config_->jetsPtMax().size() > 0 && config_->jetsPtMax()[j] > config_->jetsPtMin()[j] ) pt_max = config_->jetsPtMax()[j];
   
   bool isgood = this -> selectionJet(r,pt_min,eta_max,pt_max);
   
   return isgood;
}

bool JetAnalyser::selectionJet ( const int &  r, const float &  pt_min, const float &  eta_max, const float &  pt_max = -1.  )

virtual

Given the ranking 'r' of a jet, it returns whether the jet passes the pt_min and |eta_max|, optionally pt_max, where the values of the thresholds pt_min and |eta_max|, pt_max (opt) are passed as arguments.

Definition at line 394 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), and selectedJets_.

{
   if ( r > config_->nJetsMin() ) return true;
   if ( (int)selectedJets_.size() < r ) return false; // is this correct?
   
   bool isgood = true;
   
   std::string label = Form("Jet %d: pt > %5.1f GeV and |eta| < %3.1f",r,pt_min, eta_max );
   if ( pt_max > pt_min )
      label = Form("Jet %d: pt > %5.1f GeV and pt < %5.1f GeV and |eta| < %3.1f",r,pt_min, pt_max, eta_max );
   
   int j = r-1;
   
   // kinematic selection
   if ( selectedJets_[j] -> pt() < pt_min           && !(pt_min < 0) )  isgood = false;
   if ( fabs(selectedJets_[j] -> eta()) > eta_max   && !(eta_max < 0) ) isgood = false;
   if ( config_->jetsPtMax().size() > 0 )
   {
      if ( selectedJets_[j] -> pt() > pt_max && !(pt_max < pt_min ) )   isgood = false;
   }
   
   cutflow(label,isgood);
   
   return isgood;
}

bool JetAnalyser::selectionJetDeta ( const int &  r1, const int &  r2  )

virtual

Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_eta selection; the threshold delta is passed by the configuration file.

Definition at line 445 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_.

Referenced by analysis::mssmhbb::MssmHbbAnalyser::jetSelection().

{
   bool ok = true;
   if (config_->jetsDetaMax() < 0 )
   {
      ok = ok && true;
   }
   else
   {
      ok = ok && selectionJetDeta(r1,r2,config_->jetsDetaMax());
   }
   
   if (config_->jetsDetaMin() < 0 )
   {
      ok = ok && true;
   }
   else
   {
      ok = ok && selectionJetDeta(r1,r2,-1*config_->jetsDetaMin());
   }
   return ok;
   
}

bool JetAnalyser::selectionJetDeta ( const int &  r1, const int &  r2, const float &  delta  )

virtual

Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_eta selection; the threshold delta is passed as an argument (if delta <= 0, |delta| = delta_min; else |delta| = delta_max )

Definition at line 421 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), and selectedJets_.

{
   if ( r1 > config_->nJetsMin() ||  r2 > config_->nJetsMin() || delta == 0 ) return true;
   
   bool isgood = true;
   
   std::string label = Form("Deta(jet %d, jet %d) < %4.2f",r1,r2,fabs(delta));
   if ( delta < 0 )
      label = Form("Deta(jet %d, jet %d) > %4.2f",r1,r2,fabs(delta));
   
   int j1 = r1-1;
   int j2 = r2-1;
   
   if ( delta > 0 )
      isgood = ( fabs(selectedJets_[j1]->eta() - selectedJets_[j2]->eta()) < fabs(delta) );
   else
      isgood = ( fabs(selectedJets_[j1]->eta() - selectedJets_[j2]->eta()) > fabs(delta) );

   cutflow(label,isgood);    
    
   return isgood;
   
}

bool JetAnalyser::selectionJetDphi ( const int &  r1, const int &  r2  )

virtual

Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_phi selection; the threshold delta is passed by the configuration file.

Definition at line 492 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_.

Referenced by main().

{
   bool ok = true;
   if (config_->jetsDphiMax() < 0 )
   {
      ok = ok && true;
   }
   else
   {
      ok = ok && selectionJetDphi(r1,r2,config_->jetsDphiMax());
   }
   
   if (config_->jetsDphiMin() < 0 )
   {
      ok = ok && true;
   }
   else
   {
      ok = ok && selectionJetDphi(r1,r2,-1*config_->jetsDphiMin());
   }
   return ok;
   
}

bool JetAnalyser::selectionJetDphi ( const int &  r1, const int &  r2, const float &  delta  )

virtual

Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_phi selection; the threshold delta is passed as an argument (if delta <= 0, |delta| = delta_min; else |delta| = delta_max )

Definition at line 469 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), and selectedJets_.

{
   if ( r1 > config_->nJetsMin() ||  r2 > config_->nJetsMin()  || delta == 0 ) return true;

   bool isgood = true;
   
   std::string label = Form("Dphi(jet %d, jet %d) < %4.2f",r1,r2,fabs(delta));
   if ( delta < 0 )
      label = Form("Dphi(jet %d, jet %d) > %4.2f",r1,r2,fabs(delta));
   
   int j1 = r1-1;
   int j2 = r2-1;
   
   if ( delta > 0 )
      isgood = ( fabs(selectedJets_[j1]->deltaPhi(*selectedJets_[j2])) < fabs(delta) );
   else
      isgood = ( fabs(selectedJets_[j1]->deltaPhi(*selectedJets_[j2])) > fabs(delta) );
        
   cutflow(label,isgood);
    
   return isgood;
   
}

bool JetAnalyser::selectionJetDr ( const int &  r1, const int &  r2  )

virtual

Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_R selection; the threshold delta is passed by the configuration file.

Definition at line 541 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_.

Referenced by analysis::mssmhbb::MssmHbbAnalyser::jetSelection(), and main().

{
   bool ok = true;
   if (config_->jetsDrMax() < 0 )
   {
      ok = ok && true;
   }
   else
   {
      ok = ok && selectionJetDr(r1,r2,config_->jetsDrMax());
   }
   
   if (config_->jetsDrMin() < 0 )
   {
      ok = ok && true;
   }
   else
   {
      ok = ok && selectionJetDr(r1,r2,-1*config_->jetsDrMin());
   }
   return ok;
}

bool JetAnalyser::selectionJetDr ( const int &  r1, const int &  r2, const float &  delta  )

virtual

Given the rankings r1 and r2 of two jets, it returns whether the jets passes the delta_R selection; the threshold delta is passed as an argument (if delta <= 0, |delta| = delta_min; else |delta| = delta_max )

Definition at line 517 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), and selectedJets_.

{
   if ( r1 > config_->nJetsMin() ||  r2 > config_->nJetsMin() || delta == 0 ) return true;

   bool isgood = true;
   
   std::string label = Form("DR(jet %d, jet %d) < %4.2f",r1,r2,fabs(delta));
   if ( delta < 0 )
      label = Form("DR(jet %d, jet %d) > %4.2f",r1,r2,fabs(delta));
      
   int j1 = r1-1;
   int j2 = r2-1;
   
   if ( delta > 0 )
      isgood = ( selectedJets_[j1]->deltaR(*selectedJets_[j2]) < fabs(delta) );
   else
      isgood = ( selectedJets_[j1]->deltaR(*selectedJets_[j2]) > fabs(delta) );

   cutflow(label,isgood);
    
   return isgood;
   
}

bool JetAnalyser::selectionJetId ( )

virtual

Definition at line 624 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), jetsanalysis_, and selectedJets_.

Referenced by main(), and analysis::tools::Analyser::preselection().

{
   if ( ! jetsanalysis_ ) return true;
   
   bool isgood = true;
   std::string label = Form("JetId: %s",config_->jetsId().c_str());
   
   auto jet = std::begin(selectedJets_);
   while ( jet != std::end(selectedJets_) )
   {
      if ( ! (*jet)->id(config_->jetsId() ) )
         jet = selectedJets_.erase(jet);
      else
         ++jet;
   }
   isgood = ( selectedJets_.size() > 0 );
   
   cutflow(label,isgood);
   
   return isgood;
}

bool JetAnalyser::selectionJetPileupId ( )

virtual

Definition at line 646 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), jetsanalysis_, and selectedJets_.

Referenced by main(), and analysis::tools::Analyser::preselection().

{
   if ( ! jetsanalysis_ ) return true;
   
   bool isgood = true;
   std::string label = Form("JetPileupId: %s",config_->jetsPuId().c_str());
   
   auto jet = std::begin(selectedJets_);
   while ( jet != std::end(selectedJets_) )
   {
      if ( ! (*jet)->pileupJetIdFullId(config_->jetsPuId()) )
         jet = selectedJets_.erase(jet);
      else
         ++jet;
   }
   
   isgood = ( selectedJets_.size() > 0 );
   
   cutflow(label,isgood);
   
   return isgood;

}

bool JetAnalyser::selectionJetPtImbalance ( const int &  r1, const int &  r2  )

virtual

Given the rankings r1 and r2 of two jets, it returns whether the jets passes the pt imbalance selection; the threshold delta is passed by the configuration file.

Definition at line 588 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_.

{
   bool ok = true;
   if (config_->jetsPtImbalanceMax() < 0 )
   {
      ok = ok && true;
   }
   else
   {
      ok = ok && selectionJetPtImbalance(r1,r2,config_->jetsPtImbalanceMax());
   }
   
   if (config_->jetsPtImbalanceMin() < 0 )
   {
      ok = ok && true;
   }
   else
   {
      ok = ok && selectionJetPtImbalance(r1,r2,-1*config_->jetsPtImbalanceMin());
   }
   return ok;
   
}

bool JetAnalyser::selectionJetPtImbalance ( const int &  r1, const int &  r2, const float &  delta  )

virtual

Given the rankings r1 and r2 of two jets, it returns whether the jets passes the pt imbalance selection; the threshold delta is passed as an argument (if delta <= 0, |delta| = delta_min; else |delta| = delta_max )

Definition at line 564 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), jetsanalysis_, and selectedJets_.

{
   if ( ! jetsanalysis_ ) return true;
   if ( r1 > config_->nJetsMin() ||  r2 > config_->nJetsMin() || delta == 0 ) return true;
   
   bool isgood = true;
   std::string label = Form("DpT(jet %d, jet %d)/jet %d pT < %4.2f",r1,r2,r1,fabs(delta));
   if ( delta < 0 )
      label = Form("DpT(jet %d, jet %d)/jet %d pT > %4.2f",r1,r2,r1,fabs(delta));
   
   int j1 = r1-1;
   int j2 = r2-1;
   
   if ( delta > 0 )
      isgood = ( fabs(selectedJets_[j1]->pt() - selectedJets_[j2]->pt())/selectedJets_[j1]->pt() < fabs(delta) );
   else
      isgood = ( fabs(selectedJets_[j1]->pt() - selectedJets_[j2]->pt())/selectedJets_[j1]->pt() > fabs(delta) );

        
   cutflow(label,isgood);
    
   return isgood;
   
}

bool JetAnalyser::selectionJetQGlikelihood ( const int &  r, const float &  cut  )

virtual

Definition at line 1125 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::cutflow(), and selectedJets_.

Referenced by selectionJetQGlikelihood().

{
   bool isgood = true;
   std::string label = Form("Jet %d Q-G likelihood < %4.2f",r,fabs(cut));
   
   int j = r-1;
   
   if ( cut < 0 )
      label = Form("Jet %d Q-G likelihood > %4.2f",r,fabs(cut));
   
   if ( cut > 0 )
      isgood = ( selectedJets_[j]->qgLikelihood() < fabs(cut) );
   else
      isgood = ( selectedJets_[j]->qgLikelihood() > fabs(cut) );
        
   cutflow(label,isgood);
    
   return isgood;
   
}

bool JetAnalyser::selectionJetQGlikelihood ( const int &  r )

virtual

Definition at line 1146 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, naf_mult_submit::int, and selectionJetQGlikelihood().

{
   int j = r-1;
   bool ok = true;
   if ( config_->jetsQGmax().size() == 0 || config_->jetsQGmax()[j] < 0 || (int)config_->jetsQGmax().size() < r )
   {
      ok = ok && true;
   }
   else
   {
      ok = ok && selectionJetQGlikelihood(r,config_->jetsQGmax()[j]);
   }
   
   if (config_->jetsQGmin().size() == 0 || config_->jetsQGmin()[j] < 0 || (int)config_->jetsQGmin().size() < r )
   {
      ok = ok && true;
   }
   else
   {
      ok = ok && selectionJetQGlikelihood(r,-1*config_->jetsQGmin()[j] );
   }
   return ok;
   
}

bool JetAnalyser::selectionNJets ( )

virtual

Definition at line 670 of file JetAnalyser.cc.

References analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), naf_mult_submit::int, and selectedJets_.

Referenced by analysis::mssmhbb::MssmHbbAnalyser::jetSelection(), and main().

{
   if ( config_->nJetsMin() < 0 ) return true;
   
   bool isgood = true;
   std::string label;
   
   if ( config_->nJetsMax() <= 0 )
   {
      label = Form("NJets >= %d",config_->nJetsMin());
      isgood = ((int)selectedJets_.size() >= config_->nJetsMin());
   }
   else if ( config_->nJets() >= 0 )
   {
      label = Form("NJets = %d",config_->nJets());
      isgood = ((int)selectedJets_.size() == config_->nJets());
   }
   else
   {
      if ( config_->nJetsMin() == 0 )
      {
         label = Form("NJets <= %d",config_->nJetsMax());
         isgood = ((int)selectedJets_.size() <= config_->nJetsMax());
      }
      else
      {
         label = Form("%d <= NJets <= %d",config_->nJetsMin(),config_->nJetsMax());
         isgood =  ((int)selectedJets_.size() >= config_->nJetsMin() && (int)selectedJets_.size() <= config_->nJetsMax());
      }
   }
   
   cutflow(label,isgood);
   
   return isgood;
   
}

bool JetAnalyser::selectionNonBJet ( const int &  r )

virtual

Definition at line 728 of file JetAnalyser.cc.

References btag(), analysis::tools::BaseAnalyser::config_, analysis::tools::BaseAnalyser::cutflow(), and selectedJets_.

Referenced by selectionBJet().

{
   if ( config_->btagWP(config_->revBtagWP()) < 0 ) return true; // there is no selection here, so will not update the cutflow
   
   bool isgood = true;
   std::string label = Form("Jet %d: %s btag < %6.4f (%s) [reverse btag]",r,config_->btagAlgorithm().c_str(),config_->btagWP(config_->revBtagWP()),config_->revBtagWP().c_str());
   
   int j = r-1;
   
   // jet  non btag
   isgood = ( btag(*selectedJets_[j],config_->btagAlgorithm()) < config_->btagWP(config_->revBtagWP()) );
   
   cutflow(label,isgood);
   
   return isgood;
}

Member Data Documentation

bool analysis::tools::JetAnalyser::applyjer_

protected

Definition at line 53 of file JetAnalyser.h.

Referenced by actionApplyJER(), and JetAnalyser().

std::map<std::string, std::shared_ptr<BTagCalibrationReader> > analysis::tools::JetAnalyser::bsf_reader_

protected

Definition at line 55 of file JetAnalyser.h.

Referenced by btagSF(), and JetAnalyser().

std::vector<std::string> analysis::tools::JetAnalyser::flavours_

protected

Definition at line 59 of file JetAnalyser.h.

Referenced by JetAnalyser(), and jetHistograms().

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

protected

Definition at line 57 of file JetAnalyser.h.

Referenced by actionApplyJER(), and JetAnalyser().

std::vector< std::shared_ptr<Jet> > analysis::tools::JetAnalyser::jets_

protected

Definition at line 46 of file JetAnalyser.h.

Referenced by analysisWithJets(), and jets().

bool analysis::tools::JetAnalyser::jetsanalysis_

protected

Definition at line 52 of file JetAnalyser.h.

Referenced by actionApplyJER(), analysisWithJets(), JetAnalyser(), analysis::mssmhbb::MssmHbbAnalyser::muonJet(), selectionJetId(), selectionJetPileupId(), and selectionJetPtImbalance().

int analysis::tools::JetAnalyser::n_hjets_

protected

Definition at line 50 of file JetAnalyser.h.

Referenced by fillJetHistograms(), and jetHistograms().

std::vector< std::shared_ptr<Jet> > analysis::tools::JetAnalyser::selectedJets_

protected

Definition at line 47 of file JetAnalyser.h.

Referenced by actionApplyBjetRegression(), actionApplyJER(), analysisWithJets(), btagSF(), fillJetHistograms(), analysis::mssmhbb::MssmHbbAnalyser::fillMssmHbbHistograms(), analysis::mssmhbb::MssmHbbAnalyser::fillMssmHbbTree(), jetSwap(), analysis::mssmhbb::MssmHbbAnalyser::muonJet(), analysis::tools::Analyser::muonJet(), onlineBJetMatching(), onlineJetMatching(), selectedJets(), selectionBJet(), selectionBJetProbB(), selectionBJetProbBB(), selectionBJetProbC(), selectionBJetProbG(), selectionBJetProbLepB(), selectionBJetProbLight(), selectionDiJetMass(), selectionJet(), selectionJetDeta(), selectionJetDphi(), selectionJetDr(), selectionJetId(), selectionJetPileupId(), selectionJetPtImbalance(), selectionJetQGlikelihood(), selectionNJets(), and selectionNonBJet().

---

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

• JetAnalyser.h
• JetAnalyser.cc