{
TString kTracksName = "PicoTracks";
TString kInTracksName = "tracks";
TString kClusName = "caloClusters";
TString kCorrClusName = "CaloClustersCorr";
TString kEmcalTracksName = "EmcalTracks";
TString kEmcalClusName = "EmcalClusters";

TString kClusRandName = "CaloClustersCorrRandomized";
TString kClusEmbName = "CaloClustersCorrEmbedded";

TString kTrackRandName = "PicoTracksRandomized";
TString kTrackEmbName = "PicoTracksEmbedded";
TString kTrackEmcalEmbName = "PicoTracksEmcalEmbedded";
TString kTrackRandEmbName = "PicoTracksRandomizedEmbedded";
TString kTrackRandEmbSpectrumName = "PicoTracksRandomizedEmbeddedSpectrum";
TString kTrackEmbSpectrumName = "PicoTracksEmbeddedSpectrum";

TString kTpcKtJetsName(Form("Jet_KTChargedR020_%s_pT0150",kTracksName.Data()));
TString kTpcKtJets1GeVName(Form("Jet_KTChargedR020_%s_pT1000",kTracksName.Data()));
TString kTpcJetsR02Name(Form("Jet_AKTChargedR020_%s_pT0150",kTracksName.Data()));
TString kTpcJetsR03Name(Form("Jet_AKTChargedR030_%s_pT0150",kTracksName.Data()));
TString kTpcJetsR04Name(Form("Jet_AKTChargedR040_%s_pT0150",kTracksName.Data()));
TString kTpcEmbJetsR02Name(Form("Jet_AKTChargedR020_%s_pT0150",kTrackEmbName.Data()));
TString kTpcEmbJetsR03Name(Form("Jet_AKTChargedR030_%s_pT0150",kTrackEmbName.Data()));
TString kTpcEmbJetsR04Name(Form("Jet_AKTChargedR040_%s_pT0150",kTrackEmbName.Data()));

// Toy model
TString kTpcKtEmbJetsName(Form("Jet_KTChargedR020_%s_pT0150",kTrackEmbName.Data()));
TString kTpcKtRandJetsName(Form("Jet_KTChargedR020_%s_pT0150",kTrackRandName.Data()));
TString kTpcEmbSpectrumJets02Name(Form("Jet_AKTChargedR020_%s_pT0150",kTrackEmbSpectrumName.Data()));
TString kTpcEmbSpectrumJets03Name(Form("Jet_AKTChargedR030_%s_pT0150",kTrackEmbSpectrumName.Data()));
TString kTpcEmbSpectrumJets04Name(Form("Jet_AKTChargedR040_%s_pT0150",kTrackEmbSpectrumName.Data()));
TString kTpcRandJets02Name(Form("Jet_AKTChargedR020_%s_pT0150",kTrackRandName.Data()));
TString kTpcRandJets03Name(Form("Jet_AKTChargedR030_%s_pT0150",kTrackRandName.Data()));
TString kTpcRandJets04Name(Form("Jet_AKTChargedR040_%s_pT0150",kTrackRandName.Data()));
TString kTpcRandJets02_10GeVName(Form("Jet_AKTChargedR020_%s_pT10000",kTrackRandName.Data()));
TString kTpcRandJets03_10GeVName(Form("Jet_AKTChargedR030_%s_pT10000",kTrackRandName.Data()));
TString kTpcRandJets04_10GeVName(Form("Jet_AKTChargedR040_%s_pT10000",kTrackRandName.Data()));
TString kTpcRandEmbSpectrumJets02Name(Form("Jet_AKTChargedR020_%s_pT0150",kTrackRandEmbSpectrumName.Data()));
TString kTpcRandEmbSpectrumJets03Name(Form("Jet_AKTChargedR030_%s_pT0150",kTrackRandEmbSpectrumName.Data()));
TString kTpcRandEmbSpectrumJets04Name(Form("Jet_AKTChargedR040_%s_pT0150",kTrackRandEmbSpectrumName.Data()));
TString kTpcRandEmbJets02Name(Form("Jet_AKTChargedR020_%s_pT0150",kTrackRandEmbName.Data()));
TString kTpcRandEmbJets03Name(Form("Jet_AKTChargedR030_%s_pT0150",kTrackRandEmbName.Data()));
TString kTpcRandEmbJets04Name(Form("Jet_AKTChargedR040_%s_pT0150",kTrackRandEmbName.Data()));

TString kEmcalKtJetsName(Form("Jet_KTFullR020_%s_pT0150_%s_ET0300",kTracksName.Data(),kCorrClusName.Data()));
TString kEmcalKtJets1GeVName(Form("Jet_KTFullR020_%s_pT1000_%s_ET1000",kTracksName.Data(),kCorrClusName.Data()));

TString kEmcalJets015Name(Form("Jet_AKTFullR015_%s_pT0150_%s_ET0300",kTracksName.Data(),kCorrClusName.Data()));
TString kEmcalJets02Name(Form("Jet_AKTFullR020_%s_pT0150_%s_ET0300",kTracksName.Data(),kCorrClusName.Data()));
TString kEmcalJets025Name(Form("Jet_AKTFullR025_%s_pT0150_%s_ET0300",kTracksName.Data(),kCorrClusName.Data()));
TString kEmcalJets03Name(Form("Jet_AKTFullR030_%s_pT0150_%s_ET0300",kTracksName.Data(),kCorrClusName.Data()));

TString kEmcalEmbJets015Name(Form("Jet_AKTFullR015_%s_pT0150_%s_ET0300",kTrackEmcalEmbName.Data(),kCorrClusName.Data()));
TString kEmcalEmbJets02Name(Form("Jet_AKTFullR020_%s_pT0150_%s_ET0300",kTrackEmcalEmbName.Data(),kCorrClusName.Data()));
TString kEmcalEmbJets025Name(Form("Jet_AKTFullR025_%s_pT0150_%s_ET0300",kTrackEmcalEmbName.Data(),kCorrClusName.Data()));
TString kEmcalEmbJets03Name(Form("Jet_AKTFullR030_%s_pT0150_%s_ET0300",kTrackEmcalEmbName.Data(),kCorrClusName.Data()));

TString kEmcalJets015_1GeVName(Form("Jet_AKTFullR015_%s_pT1000_%s_ET1000",kTracksName.Data(),kCorrClusName.Data()));
TString kEmcalJets02_1GeVName(Form("Jet_AKTFullR020_%s_pT1000_%s_ET1000",kTracksName.Data(),kCorrClusName.Data()));
TString kEmcalJets025_1GeVName(Form("Jet_AKTFullR025_%s_pT1000_%s_ET1000",kTracksName.Data(),kCorrClusName.Data()));
TString kEmcalJets03_1GeVName(Form("Jet_AKTFullR030_%s_pT1000_%s_ET1000",kTracksName.Data(),kCorrClusName.Data()));

TString kEmcalJets02_3GeVName(Form("Jet_AKTFullR020_%s_pT3000_%s_ET3000",kTracksName.Data(),kCorrClusName.Data()));

TString kEmcalEmbJets015_1GeVName(Form("Jet_AKTFullR015_%s_pT1000_%s_ET1000",kTrackEmcalEmbName.Data(),kCorrClusName.Data()));
TString kEmcalEmbJets02_1GeVName(Form("Jet_AKTFullR020_%s_pT1000_%s_ET1000",kTrackEmcalEmbName.Data(),kCorrClusName.Data()));
TString kEmcalEmbJets025_1GeVName(Form("Jet_AKTFullR025_%s_pT1000_%s_ET1000",kTrackEmcalEmbName.Data(),kCorrClusName.Data()));
TString kEmcalEmbJets03_1GeVName(Form("Jet_AKTFullR030_%s_pT1000_%s_ET1000",kTrackEmcalEmbName.Data(),kCorrClusName.Data()));

TString kEmcalJetsClusOnlyName(Form("Jet_AKTNeutralR020_%s_ET0300",kCorrClusName.Data()));
TString kEmcalEmbJetsClusOnlyName(Form("Jet_AKTNeutralR020_%s_ET0300",kClusEmbName.Data()));

TString kTpcRhoName = "TpcRho";
TString kTpcRhoNameExLJ = "TpcRho_ExLJ";
TString kTpcRhoSmallName = "TpcRho_Small";
TString kTpcRhoMeth4Name = "TpcRhoAverage";
TString kTpcRhoSmallMeth4Name = "TpcRhoAverage_Small";
TString kTpcEmbRhoName = "TpcEmbRho";
TString kTpcRandRhoName = "TpcRandRho";
TString kTpcRho1GeVNameExLJ = "TpcRho_1GeV";
TString kTpcRhoSmall1GeVNameExLJ = "TpcRho_Small_1GeV";

TString kEmcalRhoMeth1Name = "EmcalRho";
TString kEmcalRhoMeth1NameExLJ = "EmcalRho_ExLJ";
TString kEmcalRhoClusOnlyMeth1Name = "EmcalRhoClusOnly";
TString kEmcalRho1GeVMeth1Name = "RhoMeth1_1GeV";
TString kEmcalRhoMeth2Name = "TpcRho_Small_Scaled";
TString kEmcalRho1GeVMeth2Name = "TpcRho_Small_1GeV_Scaled";
TString kEmcalRhoMeth3Name = "EmcalRhoMeth3";
TString kEmcalRhoMeth4Name = "EmcalRhoAverage";

Double_t kRhoMinEta = -0.5;
Double_t kRhoMaxEta = 0.5;
Double_t kRhoMinPhi = 30*TMath::DegToRad()+0.2;
Double_t kRhoMaxPhi = 230*TMath::DegToRad()-0.2;

Double_t kPi = TMath::Pi();
Double_t kGhostArea = 0.005;
Double_t kJetRadius = 0.2;
Double_t kKtJetRadius = 0.2;
Double_t kJetAreaCut = 0.557;
Double_t kTrackPtCut = 0.15;
Double_t kClusPtCut = 0.3;
Double_t kJetPtCut = 1;
Double_t kJetBiasTrack = 5;
Double_t kJetBiasClus = 1000;
UInt_t kEmcalAna = 1;
UInt_t kTpcAna = 0;
UInt_t kUserAna = 2;
//UInt_t kPhysSel = AliVEvent::kAnyINT + AliVEvent::kSemiCentral + AliVEvent::kCentral;
//UInt_t kPhysSel = -2147483648; //AliEmcalPhysicsSelection::kEmcalOk
UInt_t kPhysSel = AliVEvent::kAny;
UInt_t prePhysSel = AliVEvent::kAnyINT + AliVEvent::kSemiCentral + AliVEvent::kCentral; 
//UInt_t prePhysSel = AliVEvent::kEMCEGA;
//UInt_t prePhysSel = AliVEvent::kEMCEJE;

TF1* sfunc = new TF1("sfunc","[0]*x*x+[1]*x+[2]",-1,100);

// 300 MeV, had corr = 2
  sfunc->SetParameter(2,1.80183);
  sfunc->SetParameter(1,-0.0134124);
  sfunc->SetParameter(0,0.000146031);

// 300 MeV, had corr = 1.7
/*
  sfunc->SetParameter(2, 1.79148);
  sfunc->SetParameter(1, -1.13795e-02);
  sfunc->SetParameter(0, 1.14528e-04);
*/

// 300 MeV, had corr = 1.3
/*
  sfunc->SetParameter(2,1.90655);
  sfunc->SetParameter(1,-0.0159274);
  sfunc->SetParameter(0,0.0001706891);
*/

// 3x3, 300 MeV, had corr = 2
/*
  sfunc->SetParameter(2, 2.02415);
  sfunc->SetParameter(1, -8.13464e-03);
  sfunc->SetParameter(0,2.28091e-05);
*/

TF1* rfunc = new TF1("rfunc","[0]*x*x*x*x+[1]*x*x*x+[2]*x*x+[3]*x+[4]",-1,100);
  rfunc->SetParameter(0,0.00000102399);
  rfunc->SetParameter(1,-0.000584397);
  rfunc->SetParameter(2,0.125030);
  rfunc->SetParameter(3,-9.646059);
  rfunc->SetParameter(4,260.343);

TF1* sfunc1GeV = new TF1("sfunc1GeV","[0]*x*x+[1]*x+[2]",-1,100);

// 1 GeV, had corr = 2
  sfunc1GeV->SetParameter(2,1.52519);
  sfunc1GeV->SetParameter(1,-1.04701e-02);
  sfunc1GeV->SetParameter(0, 8.99280e-05);
}