message: Improved parameterization of Run2TuneMCNoNcell
diff: @@ -876,41 +876,39 @@ Double_t AliPHOSTenderSupply::CorrectNonlinearity(Double_t en){
  if(fNonlinearityVersion=="Run2TuneMCNoNcell"){ //Improved Run2 tune for MC in the case of loose cluster cuts (no Ncell>2 cut)
  if(en<=0.) return 0.;
  - const Double_t xMin=0.57; //low part of the param (optimized from pi0 peak)
  + const Double_t xMin=0.850; //low part of the param (optimized from pi0 peak)
  const Double_t xMax=5.17; //Upper part of the param (optimized from pi0 peak)
  //middle part param
  const Double_t a= 1.02165 ;
  const Double_t b=-2.548e-01 ;
  - const Double_t c= 0.6413 ;
  - const Double_t d=-0.4760 ;
  - const Double_t e= 0.1165 ;
  + const Double_t c= 0.6483 ;
  + const Double_t d=-0.4980 ;
  + const Double_t e= 0.1245 ;
  Double_t ecorr=0.;
  if(en<xMin){
  - const Double_t beta = a+0.5*b/sqrt(xMin)-0.5*d/(xMin*sqrt(xMin))-e/(xMin*xMin);
  - const Double_t alpha = a*xMin+b*sqrt(xMin)+c+d/sqrt(xMin)+e/xMin-beta*xMin ;
  - ecorr= 1.007102*(alpha+beta*en) ;
  + const Double_t gamma = 0.150 ;
  + const Double_t beta = 0.5*(0.5*b*sqrt(xMin)+c+1.5*d/sqrt(xMin)+2.*e/xMin)*TMath::Power((xMin*xMin+gamma*gamma),2)/
  + (xMin*xMin*xMin);
  + const Double_t alpha = (a*xMin+b*sqrt(xMin)+c+d/sqrt(xMin)+e/xMin-beta*xMin/(xMin*xMin+gamma*gamma))/xMin ;
  + ecorr= 1.0165898*(alpha*en+beta*en/(en*en+gamma*gamma)) ;
  }
  else{
  if(en<xMax){
  - ecorr= 1.007102*(a*en+b*sqrt(en)+c+d/sqrt(en)+e/en) ;
  + ecorr= 1.0165898*(a*en+b*sqrt(en)+c+d/sqrt(en)+e/en) ;
  }
  else{
  const Double_t beta= b+2.*c/sqrt(xMax)+3.*d/xMax+4.*e/xMax/sqrt(xMax) ;
  const Double_t alpha = a+b/sqrt(xMax)+c/xMax+d/xMax/sqrt(xMax)+e/(xMax*xMax)-beta/sqrt(xMax) ;
  - ecorr= 1.007102*(alpha*en+beta*sqrt(en)) ;
  + ecorr= 1.0165898*(alpha*en+beta*sqrt(en)) ;
  }
  }
  if(ecorr<0){
  return 0.;
  }
  - Double_t aM= 0.0226 ;
  - Double_t bM= 0.415 ;
  - Double_t cM= 0.09 ;
  - return ecorr*(1.+aM*(1.-TMath::TanH((en-bM)/cM))) ;
  -
  + return ecorr ;
  }
  return en ;
  