message: Option to use dN/dy for upper and lower FONLL predictions
diff: enum EDDecay{kD0Kpi,kDplusKpipi,kDstarD0pi,kDsKKpi,kLcpKpi,kLcK0Sp,kLcpiL,kDplusKKpi};
  enum EFidY{kFixedY,kPtDepY};
  enum EPtShape{kFlatPt,kPtFONLL8TeV,kPtFONLL8TeVfeeddown,kPtFONLL7TeV,kPtPythia7TeV,kPtFONLL5TeV,kPtFONLL13TeVprompt,kPtFONLL13TeVfeeddown,kPtPythia13TeVprompt,kPtPythia13TeVfeeddown};
  -enum EYShape{kFlatY,kYFONLL5TeV};
  +enum EYShape{kFlatY,kYFONLL5TeV,kYFONLLMax5TeV,kYFONLLMin5TeV};
  // Configuration
  Int_t fDDecay=kD0Kpi;
  @@ -31,7 +31,7 @@ Double_t fEtaMaxDau=0.9;
  Int_t fOptionYFiducial=kPtDepY;
  Double_t fYMaxFidAccCut=0.8;
  Int_t fPtShape=kPtFONLL5TeV;
  -Int_t fYShape=kFlatY;
  +Int_t fYShape=kYFONLL5TeV;
  TString fDecayTableFileName="$ALICE_PHYSICS/PWGHF/vertexingHF/macros/decaytable_acc.dat";
  Int_t fDebugLevel=0;
  Int_t totTrials=10000000;
  @@ -252,6 +252,14 @@ void ComputeAcceptance(){
  funcY=new TF1("fsigyfonll5","[0]+((x>0.5)*[1]/sqrt(x))+((x<0.5)*([1]*sqrt(2)))",0.,50);
  funcY->SetParameters(7.22252e-01,5.06914);
  outFileName.Append("FONLLy.root");
  + }else if(fYShape==kYFONLLMax5TeV){
  + funcY=new TF1("fsigyfonll5","TMath::Min(8.,[0]+[1]/sqrt(x))",0.,50);
  + funcY->SetParameters(7.05870e-01,5.11006);
  + outFileName.Append("FONLLyMax.root");
  + }else if(fYShape==kYFONLLMin5TeV){
  + funcY=new TF1("fsigyfonll5","TMath::Min(8.,[0]+[1]/sqrt(x)+[2]/(x*x*x))",0.,50);
  + funcY->SetParameters(6.99901e-01,5.27825,1.54156e+01);
  + outFileName.Append("FONLLyMin.root");
  }else{
  outFileName.Append("flaty.root");
  }
  