AtlfastB.h

Go to the documentation of this file.

// ================================================
// AtlfastB class description
// ================================================
//
// Namespace Atlfast::
//
// class: AtlfastB
//
// It is based upon atlfastb.f by E.Richter-Was
//
//  Author    : Mireia Dosil
//
//
// Description: 
//
// The original jet flavour labels from JetMaker are rearranged here to 
// simulate the overall efficiencies measured in full simulation.
// Correction factors are also computed to match different tagging hypotheses.
// 
//
// ................................................................
//

#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

#include "CLHEP/Random/RandomEngine.h"
#include "CLHEP/Random/RandFlat.h"

#include "AtlfastEvent/Interpolator.h"
#include "AtlfastEvent/CollectionDefs.h"

class MsgStream;
class ISvcLocator;

namespace Atlfast {


  class TesIO;
  class Jet;

  using std::string;

  class AtlfastB : public Algorithm
{

 private:
  
  int m_atlfBNSet;  
  
  double m_epsib;   
  double m_epsic;   
  double m_epsij;   
  double m_epsitau; 

  double m_epsitau1P; 
  double m_epsitau3P; 

  string m_Tau1P3Pcorrfile;  

  double m_corr1P3P[20];        
  double m_corr_prob1P[6][20];  
  double m_corr_prob3P[6][20];  

  //input files for jet correction factors

  string m_corrjfile;  
  string m_corrcfile;  

  bool   m_useTDRBParam;         
  bool   m_interpolateBTagging;  

  double m_epsilonBjet;          

  double m_correctionFactor;

  Interpolator* m_beff_interpolator;     
  Interpolator* m_brejpu_interpolator;   
  Interpolator* m_brejnpu_interpolator;  
  Interpolator* m_brejtau_interpolator;  
  Interpolator* m_brejc_interpolator;    

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

  int m_indtauveto;  
  


  //switches for the different methods
  bool m_AtlfBJetSwitch;     
  bool m_AtlfCalSwitch;      
  bool m_AtlfTauSwitch;      
  bool m_AtlfTauVetoSwitch;  
  bool m_AtlfTrigMuoSwitch;  
  bool m_AtlfTau1P3PSwitch;  
  
  bool m_corrtaumom;         

  
  TesIO* m_tesIO;     

  HepRandomEngine* m_pRandomEngine;  
  RandFlat* m_pRandFlatGenerator;    
  
  const JetCollection* m_Jets;         
  JetCollection::const_iterator m_it;  

  Interpolator* m_taueff_interpolator;  
  Interpolator* m_taurej_interpolator;  

  
  std::string m_inputLocation;    

  std::string m_mcLocation; 



 public:

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


  StatusCode atlfBje_TDR();           
  StatusCode atlfBje();               
  StatusCode atlfCal();               
  StatusCode atlfTau(double epsitau); 
  StatusCode atlfTau1P3P();           
  StatusCode atlfTauVeto(int ind);    
  StatusCode atlfTrigMuo();           

  StatusCode tautag(double pt, double eta, double efftau, double &rjet, int &iflag);


  StatusCode tauveto(double pt, int ind, double &efftau, double &effjet);  

  void       makeInterpolator(Interpolator* &intptr, string intname, bool contbounds);
  
  
  double fitcoreb(double pt); 
  double fitcoreu(double pt); 


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

  
}  //end of namespace bracket

#endif

---