NtupleReader/NtupleReader/MySystematics.ixx

//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
// Jamuary 2008      Ketevi A. Assamagan
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

//::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
//:: IMPLEMENTATION OF INLINE METHODS DEFINED IN THE CLASS MySystematics ::
//::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::

//*****************************************************************************

  //:::::::::::::::::
  //:: METHOD egamma_resolution ::
  //:::::::::::::::::
inline Double_t MySystematics::egamma_resolution( const Double_t et ) {
  Double_t sigma = 0.0073 * et; 
  return m_random->Gaus( et,sigma );
}

//*****************************************************************************

  //:::::::::::::::::
  //:: METHOD muon_resolution ::
  //:::::::::::::::::
inline Double_t MySystematics::muon_resolution( const Double_t pt ) {

  if (pt==0.0) return 0.0;

  Double_t inv_pt=1.0/pt ;      
  Double_t comp1 = 0.011*inv_pt ;
  Double_t comp2 = 0.00017 ;
  Double_t sigma = sqrt( comp1*comp1 + comp2*comp2 ) ;
  Double_t new_inv_pt= m_random->Gaus( inv_pt,sigma ) ;
  Double_t new_pt=1.0/new_inv_pt ;
  return new_pt;


/* (till 27.1.08) smearing in pt with some sigma specified/recommended */
//  Double_t comp1 = 0.011*pt;
//  Double_t comp2 = 0.00017 * (pt*pt);
//  Double_t sigma = sqrt( comp1*comp1 + comp2*comp2 );
//  return m_random->Gaus( pt,sigma );

}

//*****************************************************************************

  //:::::::::::::::::
  //:: METHOD jet_resolution ::
  //:::::::::::::::::
inline Double_t MySystematics::jet_resolution( const Double_t e,const Double_t eta ) {
  Double_t sigma = 0.45 * sqrt(e);
  if (TMath::Abs(eta)>3.2) sigma = 0.63 * sqrt(e);
  return m_random->Gaus( e,sigma );
}

//*****************************************************************************

  //:::::::::::::::::
  //:: METHOD taujet_resolution ::
  //:::::::::::::::::
inline Double_t MySystematics::taujet_resolution( const Double_t e ) {
  Double_t sigma = 0.45 * sqrt(e);
  return m_random->Gaus( e,sigma );
}

//*****************************************************************************

  //:::::::::::::::::
  //:: METHOD energy_scale ::
  //:::::::::::::::::
inline Double_t MySystematics::energy_scale ( const Double_t pT, const Double_t f ) {
  return f*pT;
} 

//*****************************************************************************

  //:::::::::::::::::
  //:: METHOD is_detected ::
  //:::::::::::::::::
inline Bool_t MySystematics::is_detected( const Double_t f ) {
  Double_t ran = m_random->Uniform();
  return ( ran>f );
}

//*****************************************************************************

  //:::::::::::::::::
  //:: METHOD is_bjet ::
  //:::::::::::::::::
inline Bool_t MySystematics::is_bjet( const Double_t f_minus ) {
  Double_t ran = m_random->Uniform();
  return ( ran>f_minus );
}

//*****************************************************************************

  //:::::::::::::::::
  //:: METHOD metCorrection_x ::
  //:::::::::::::::::
inline Double_t MySystematics::metCorrection_x(const Double_t aodPx, const Double_t correctedPx) {
  return ( +aodPx-correctedPx );
}

  //:::::::::::::::::
  //:: METHOD metCorrection_y ::
  //:::::::::::::::::
inline Double_t MySystematics::metCorrection_y(const Double_t aodPy, Double_t correctedPy) {
  return ( +aodPy-correctedPy );
} 

  //:::::::::::::::::
  //:: METHOD metCorrection_sum ::
  //:::::::::::::::::
inline Double_t MySystematics::metCorrection_sum(const Double_t aodSum, Double_t correctedSum) {
  return ( +aodSum-correctedSum );
} 

  //:::::::::::::::::
  //:: METHOD metScale ::
  //S.M.: I do not use it at the moment!
  //:::::::::::::::::
inline Double_t MySystematics::metScale ( const Double_t aodM, const Double_t correctedM, Double_t factor)
{
  return (correctedM-factor*aodM);
}

//*****************************************************************************

---