---

AtlfastB.h

Go to the documentation of this file.

// ================================================
// AtlfastB class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
//
// This version....
//
//
// Namespace Atlfast::
//
// class: AtlfastB
//
//It is based upon atlfastb.f by E.Richter-Was
//
//  Author    : Mireia Dosil
//
//
// Description: 
//
//  Algorithm for jet tagging and jet tagging efficiencies computation
//  Translated into c++ athena from the fortran code atlfastb.f
// 
//
// ................................................................
//

#ifndef ATLFAST_B
#define ATLFAST__B

// STL
#ifndef STD_VECTOR_H
#include <vector>
#define STD_VECTOR_H
#endif

#ifndef STD_STRING_H
#include <string>
#define STD_STRING_H
#endif

// Gaudi 
#ifndef GAUDIKERNEL_ALGORITHM_H
#include "GaudiKernel/Algorithm.h"
#define GAUDIKERNEL_ALGORITHM_H
#endif



//*************************************
//Default parameters private to this algorithm
//*************************************

#define DEFAULT_AtlfBjetSwitch     1
#define DEFAULT_AtlfCalSwitch      1
#define DEFAULT_AtlfTauSwitch      1
#define DEFAULT_AtlfTauVetoSwitch  0
#define DEFAULT_AtlfTrigMuoSwitch  0
#define DEFAULT_AtlfBNSet          1  //Choose the set of efficiencies and pt dependent factors
#define DEFAULT_TauEff             0.5 //Tau jet tagging efficiency
#define DEFAULT_TauVetoOption      1  //For the atlfTauVeto(int ind) method 
#define DEFAULT_JetCorrFile        "ptcorrectionjet.dat"
#define DEFAULT_CJetCorrFile        "ptcorrectioncjet.dat"
#define DEFAULT_RandomSeed        1234567

//*****************************************
// Default paths in the Transient event store to get/put entities
//
//******************************************

#define DEFAULT_inputLocation           "/Event/AtlfastJets" 
#define DEFAULT_outputLocation          "/Event/AtlfastCalJets" 

//***************************************************************
// AtlfastB class declaration
//
// A "Gaudi algorithm" is something which gets scheduled and controlled
// by the framework. In its simplest definition it is something which
//  - gets clled for each event 
//  - can get anything it wants out of the Transient Event Store"
//  - can put anything it makes into the store.

//
//****************************************************************
class MsgStream;
class ISvcLocator;

//using HepMC::GenEvent;

class RandFlat;
class HepRandomEngine;

namespace Atlfast {
  class TesIO;
  class Jet;

  class AtlfastB : public Algorithm
{

 private:
  
  
  //--------------------------------
  // Parameters of this algorithm
  // (descriptions above)
  //--------------------------------

  
  // Parameters for the different atlfastB methods
  // parametrization prepared by E. Ros, updated November'99, low luminosity
  // apply randomized b-tagging with efficiencies specified
  // by one of given sets, pT dependent correction factors included
  // m_atlfBNSet=1  --> m_epsib=0.5
  // m_atlfBNSet=2  --> m_epsib=0.6
  // m_atlfBNSet=3  --> m_epsib=0.7
  //                                          
  // canonical for low luminosity performance,no pT dependence
  // m_atlfBNSet=5 --> m_epsib=0.60 epsi_c=0.10, R=100
  //                                          
  // Inner Detector TDR, April'97, low luminosity, pT dependent correction
  // factors included, parametrization prepared by E. Ros
  // m_atlfBNSet=11 --> m_epsib=0.33
  // m_atlfBNSet=12 --> m_epsib=0.43
  // m_atlfBNSet=13 --> m_epsib=0.53
  // m_atlfBNSet=14 --> m_epsib=0.624
  //                                          


  int m_atlfBNSet;
  
  //Jet tagging efficiencies
  float m_epsib;
  float m_epsic;
  float m_epsij;
  float m_epsitau;

  //input files for jet correction factors

  string m_corrjfile;
  string m_corrcfile;

  //jet correction factors
  float m_corrj[5][15];
  float m_corrc[5][15];

  //
  int m_indtauveto;  
  
  // for an analysis where you want to reject taus use atlfVeto      
  // apply veto for tau-jet 
  // (m_indtauveto=1): tau-background fixed to 5%, atlfVeto gives the corresponding jet-efficiency
  // (m_indtauveto=2): jet-efficiency fixed at 90%, atlfVeto gives the corresponding tau-background


//switch for the different methods
  int m_iAtlfBje;
  int m_iAtlfCal;
  int m_iAtlfTau;
  int m_iAtlfTauVeto;
  int m_iAtlfTrigMuo;
  
  //other variables
  
  bool m_notInit;
  TesIO* m_tesIO;
  int m_randomSeed;
  HepRandomEngine* m_pRandomEngine;
  RandFlat* m_pRandFlatGenerator;
  
  std::vector<Jet*>  m_Jets;
  std::vector<Jet*>::iterator m_it;



  
  //--------------------------------------------------
  // Paths in the Transient Event store to get/put things
  // These are currently set in member variables so that
  // they can be overwritten by the job options service
  //--------------------------------------------------

  std::string m_inputLocation ;
  std::string m_outputLocation ;

 public:

  //-------------------------
  // Constructors/Destructors
  //
  // Gaudi requires that the constructor takes certain arguments
  // (and passes them directly to the constructor of the base class)
  //-------------------------

  AtlfastB( const std::string& name, ISvcLocator* pSvcLocator ) ; 
   virtual ~AtlfastB();




  //-------------------------------
  // helper
  //-------------------------------

  

  StatusCode atlfBje(int nset);
  StatusCode atlfCal();
  StatusCode atlfTau(float epsitau);
  StatusCode tautag(double pt, double eta, double efftau, double &rjet);
  StatusCode atlfTauVeto(int ind);
  StatusCode tauveto(double pt, int ind, double &efftau, double &effjet);  
  StatusCode atlfTrigMuo();
  
  
  double fitcoreb(double pt);
  double fitcoreu(double pt);

  //------------------------------------------------------
  // Methods used by Gaudi to run the algorithm
  //------------------------------------------------------

   StatusCode initialize() ;
   StatusCode execute() ;
   StatusCode finalize() ;

};

  
}  //end of namespace bracket

#endif

---