#include "AliJetHadronCorrectionv0.h"
const Int_t maxVariant = 8; // size eta grid
-const Int_t nVec = 10; // size momentum grid
-const Int_t nPol = 4; // number coefficients of polinom
+const Int_t nVec = 11; // size momentum grid
+const Int_t nPol = 4; // number coefficients of polynom
static Double_t etaGrid[maxVariant]={ 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.67};
-static Double_t pGrid[nVec]={ 0.2, 0.5, 1.0, 2.0, 3.0, 5.0,10.0,15.0,25.0,40.0};
+static Double_t pGrid[nVec]={ 0.5, 1.0, 2.0, 3.0, 5.0, 10.0, 20.0, 30.0, 50.0, 70.0, 100.0};
// c[][] - first index for eta, second for momentum
static Double_t c[maxVariant][nPol] ={
{1.305705e-01, 3.725653e-01, -1.219962e-02, 1.806235e-04},
Double_t
AliJetHadronCorrectionv0::GetEnergy(Double_t pmom, Double_t eta, Int_t /*gid*/)
{
+// Get energy for momentum pmom and pseudorapidity eta
Int_t iEta=0; // index
Double_t etaw = TMath::Abs(eta);
if(etaw > etaGrid[maxVariant-1]) etaw = etaGrid[maxVariant-1];