NtupleReader/NtupleReader/MySystematics.h

#ifndef MySystematicsH
#define MySystematicsH

//::::::::::::::::::::::
//:: CLASS MySystematics ::
//::::::::::::::::::::::


//::::::::::::::::::
//:: HEADER FILES ::
//::::::::::::::::::

// standard C++ //
#include <iostream>
#include <vector>

// ROOT //
#include <TRandom.h>
#include "TLorentzVector.h"

// HiggsAnalysis //
#include "MyParticle.h"
#include "MyMissingEt.h"
#include "MyTools.h"
#include "NtupleReader/MyDatasetInfo.h"



class MySystematics {

 public:
   
   MySystematics( unsigned int seed ){
        init(seed);
    };

   ~MySystematics(void){
        destruct();
    }
;

 public:

  Double_t egamma_resolution( const Double_t et );

  Double_t muon_resolution( const Double_t pt );

  Double_t jet_resolution( const Double_t e, const Double_t eta ); 

  Double_t taujet_resolution( const Double_t e );

  Double_t energy_scale ( const Double_t pT, const Double_t f );

  Bool_t   is_detected( const Double_t f );

  Bool_t   is_bjet( const Double_t f_minus );

  Double_t metCorrection_x(const Double_t aodPx, const Double_t correctedPx);
  Double_t metCorrection_y(const Double_t aodPy, Double_t correctedPy);
  Double_t metCorrection_sum(const Double_t aodSum, Double_t correctedsum);
  Double_t metScale ( const Double_t aodM, const Double_t correctedM, Double_t factor=1.0);

  void apply_syst( std::vector<MyParticle*> * vp_syst, MyMissingEt * met_syst, int f_reso, int f_eScale_minus, int f_eScale_plus, int f_effi, int f_btag);

 private:

  // The MyDatasetInfo object 
  MyDatasetInfo *m_dataset_info;

  // The MyTools object 
  MyTools m_tools;


  TRandom * m_random;

  void init(unsigned int seed);
         // initialization method (variables as defined
         // in the constructor of this class)
  void destruct(void);
         // destruction method

 
};

//::::::::::::::::::::::::::::::::::::::
//:: IMPLEMENTATION OF INLINE METHODS ::
//::::::::::::::::::::::::::::::::::::::

#include "MySystematics.ixx"

#endif

---