AnalyserBJetsExample.cc File Reference

#include "Analysis/Tools/interface/Analyser.h"

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 10 of file AnalyserBJetsExample.cc.

References analysis::tools::JetAnalyser::actionApplyBjetRegression(), analysis::tools::JetAnalyser::actionApplyJER(), submitCrab3Data::config, analysis::tools::BaseAnalyser::config(), analysis::tools::Analyser::event(), analysis::tools::JetAnalyser::fillJetHistograms(), analysis::tools::JetAnalyser::jetHistograms(), analysis::tools::BaseAnalyser::nEvents(), analysis::tools::JetAnalyser::selectionBJet(), analysis::tools::JetAnalyser::selectionJet(), analysis::tools::JetAnalyser::selectionJetDphi(), analysis::tools::JetAnalyser::selectionJetId(), analysis::tools::JetAnalyser::selectionJetPileupId(), and analysis::tools::JetAnalyser::selectionNJets().

{
   Analyser analyser(argc,argv);
   auto config = analyser.config();
   
// HISTOGRAMS definitions  
   // create some predefined jet histograms
   // if not defined, the number of jets is nJetMin from the configurations
   analyser.jetHistograms("selection01");
   analyser.jetHistograms("selection02");
   analyser.jetHistograms("selection03");
   
   for ( int i = 0 ; i < analyser.nEvents() ; ++i )
   {
      if ( ! analyser.event(i)                  )   continue;
      
   // JETS
      analyser.actionApplyJER();                               // correction : jet energy resolution smearing
      analyser.actionApplyBjetRegression();                    // correction : jet energy regression (for b jets)
      if ( ! analyser.selectionJetId()          )   continue;  // selection  : jet identification 
      if ( ! analyser.selectionJetPileupId()    )   continue;  // selection  : jet Pileup identification 
      if ( ! analyser.selectionNJets()          )   continue;  // selection  : number of jets 
      
      // the lines below will only be executed if rank <= nJetsMin
      if ( ! analyser.selectionJet(1)           )   continue;  // selection  : jet1 pt and eta 
      if ( ! analyser.selectionJet(2)           )   continue;  // selection  : jet2 pt and eta 
      if ( ! analyser.selectionJetDphi(1,2)     )   continue;  // selection  : delta_phi_jets (1,2) [or  MIN(neg): analyser.selectionJetDphi(1,2,-2.0) / MAX(pos): analyser.selectionJetDphi(1,2,+2.0)]
   // HISTOGRAMS
      analyser.fillJetHistograms("selection01");               // histograms : jets fill
   // BTAG
      if ( ! analyser.selectionBJet(1)          )   continue;
      analyser.fillJetHistograms("selection02");               // histograms : jets fill
      if ( ! analyser.selectionBJet(2)          )   continue;
      analyser.fillJetHistograms("selection03");               // histograms : jets fill
   // HISTOGRAMS
      
   }  //end event loop
   

} // end main