JetRecalibrator.h

Go to the documentation of this file.

//                                                                      //
// Reacalibrate jets following ATLFastSTL                               //
//                                                                      //
namespace Atlfast{

  class JetRecalibrator{
  public:

    JetRecalibrator(){}

    double bCalibrate(double pt){

      //Convert to GeV
      pt=pt/GeV;

      //b jets
      float core, a0, a1, a2, a3, a4, a5;
      
      if(pt < 10.){ return 0;
      } else if ( pt < 55.0) {    
        a0 =   1.2715;     
        a1 =   .12241;     
        a2 =  -.10480E-01;         
        a3 =  +.33310E-03;         
        a4 =  -.47454E-05;         
        a5 =   .25436E-07;         
        core =   a0 + a1*pt + a2*pt*pt + a3*pt*pt*pt 
          + a4*pt*pt*pt*pt + a5*pt*pt*pt*pt*pt;
        core =   core * 1.006;
        core = core * pt * GeV;
        return core;
      } else {
        a0 =   1.18;       
        a1 =   -.16672E-02;        
        a2 =    .44414E-05;        
        core = a0+a1*pt + a2*pt*pt;

        //convert back to original units
        core = core * pt* GeV;
        return core;
      }
    }

    double uCalibrate(float pt) {

      //Convert to GeV
      pt=pt/GeV;


      //  calibrate u-jets 

      float core, a0, a1, a2, a3, a4, a5;
      
      if(pt < 10.){ return 0;
      } else if ( pt < 45.0) {    
        a0 =  1.5085;      
        a1 =   .31468E-01;         
        a2 =  -.36973E-02;         
        a3 =  +.11220E-03;         
        a4 =  -.13921E-05;         
        a5 =   .61538E-08;         
        core =    a0 + a1*pt + a2*pt*pt + a3*pt*pt*pt 
          + a4*pt*pt*pt*pt + a5*pt*pt*pt*pt*pt;
        core = core * pt * GeV;
        return core;
       } else {
         a0 =   1.18;      
         a1 =   -.16672E-02;       
         a2 =    .44414E-05 ;      
         core =   a0 + a1*pt + a2*pt*pt;
         core =   core/1.25;

        //convert back to original units
         core = core * pt * GeV;
         return core;
       }
    }
  };
}

---