/cvmfs/atlas.cern.ch/repo/sw/ASG/AnalysisBase/2.4.31/ElectronPhotonShowerShapeFudgeTool/ElectronPhotonShowerShapeFudgeTool/TElectronMCShifterTool.h

Go to the documentation of this file.

#ifndef TELECTRONMCSHIFTERTOOL_H
#define TELECTRONMCSHIFTERTOOL_H

#include <map>
#include <vector>


namespace ElePIDNames{
  enum Var {
     DeltaPoverP
    ,TRTHighTOutliersRatio
    ,d0significance
    ,deltaeta1
    ,deltaphi2
    ,deltaphiRescaled
    ,eratio
    ,f1
    ,f3
    ,fside
    ,reta
    ,rhad
    ,rphi
    ,trackd0pvunbiased
    ,weta2
    ,ws3
    ,wstot
    ,e277
    ,DeltaE
  };
}


class TElectronMCShifterTool{
public:
  TElectronMCShifterTool();

  ~TElectronMCShifterTool();


 double getShiftDelta(double eT, double eta, ElePIDNames::Var var){
   int ibin_combined = getShifterEtBin(eT)*6+getShifterEtaBin(eta); //Convert the 2d
   return Shifts[var].at(ibin_combined);
 };


 double getShiftDeltaByBin(int eT, int eta, ElePIDNames::Var var){
   int ibin_combined = getShifterEtBin(eT)*6+getShifterEtaBin(eta); //Convert the 2d
   return Shifts[var].at(ibin_combined);
 };


 double shiftVar(double eT,double eta,ElePIDNames::Var var,double val);


 double shiftVarByBin(int eTBin,int etaBin,ElePIDNames::Var var,double val);


 void shiftAll( float eT,
                float eta,
                float& rhad1  ,
                float& rhad   ,
                float& reta   ,
                float& rphi   ,
                float& weta2  ,
                float& f1     ,
                float& f3     ,
                float& fside  ,
                float& ws3    ,
                float& wstot  ,
                float& eratio ,
                float& e277   ,
                float& DeltaE ,
                float& deltaeta1,
                float& deltaphiRescaled
              );

 //static const double Shifts[17][54]; // 17 variables x 9 eta bins x 6 et bins
 std::map<ElePIDNames::Var, std::vector<float> > Shifts;
 std::map<ElePIDNames::Var, std::vector<float> > Widths;

private :



 unsigned int getShifterEtBin(double eT) const;


 unsigned int getShifterEtaBin(double eta) const;

};

#endif

---