---

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
//***************************************************************
// 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;

  using std::string;
  
  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
    bool m_AtlfBJetSwitch;
    bool m_AtlfCalSwitch;
    bool m_AtlfTauSwitch;
    bool m_AtlfTauVetoSwitch;
    bool m_AtlfTrigMuoSwitch;
    
    //other variables
    
    bool m_notInit;
    TesIO* m_tesIO;
    HepRandomEngine* m_pRandomEngine;
    RandFlat* m_pRandFlatGenerator;
    
    std::vector<Jet*>  m_Jets;
    std::vector<Jet*>::const_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

---