{
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 = "TpcKtJets";
TString kTpcKtJets1GeVName = "TpcKtJets1GeV";
TString kTpcJetsR02Name = "TpcJetsR02";
TString kTpcJetsR03Name = "TpcJetsR03";
TString kTpcJetsR04Name = "TpcJetsR04";
TString kTpcEmbJetsR02Name = "TpcEmbJetsR02";
TString kTpcEmbJetsR03Name = "TpcEmbJetsR03";
TString kTpcEmbJetsR04Name = "TpcEmbJetsR04";

// Toy model
TString kTpcKtEmbJetsName = "TpcKtEmbJetsR02";
TString kTpcKtRandJetsName = "TpcKtRandJets";
TString kTpcEmbSpectrumJets02Name = "TpcEmbSpectrumJetsR02";
TString kTpcEmbSpectrumJets03Name = "TpcEmbSpectrumJetsR03";
TString kTpcEmbSpectrumJets04Name = "TpcEmbSpectrumJetsR04";
TString kTpcRandJets02Name = "TpcRandJetsR02";
TString kTpcRandJets03Name = "TpcRandJetsR03";
TString kTpcRandJets04Name = "TpcRandJetsR04";
TString kTpcRandJets02_10GeVName = "TpcRandJetsR02_10GeV";
TString kTpcRandJets03_10GeVName = "TpcRandJetsR03_10GeV";
TString kTpcRandJets04_10GeVName = "TpcRandJetsR04_10GeV";
TString kTpcRandEmbSpectrumJets02Name = "TpcRandEmbSpectrumJetsR02";
TString kTpcRandEmbSpectrumJets03Name = "TpcRandEmbSpectrumJetsR03";
TString kTpcRandEmbSpectrumJets04Name = "TpcRandEmbSpectrumJetsR04";
TString kTpcRandEmbJets02Name = "TpcRandEmbJetsR02";
TString kTpcRandEmbJets03Name = "TpcRandEmbJetsR03";
TString kTpcRandEmbJets04Name = "TpcRandEmbJetsR04";

TString kEmcalKtJetsName = "EmcalKtJets";
TString kEmcalKtJets1GeVName = "EmcalKtJets1GeV";
TString kEmcalJets025Name = "EmcalJetsR025";
TString kEmcalEmbJets025Name = "EmcalEmbJetsR025";
TString kEmcalJets015Name = "EmcalJetsR015";
TString kEmcalEmbJets015Name = "EmcalEmbJetsR015";
TString kEmcalJets02Name = "EmcalJetsR02";
TString kEmcalEmbJets02Name = "EmcalEmbJetsR02";
TString kEmcalJets03Name = "EmcalJetsR03";
TString kEmcalEmbJets03Name = "EmcalEmbJetsR03";
TString kEmcalJetsClusOnlyName = "EmcalJetsClusOnly";
TString kEmcalEmbJetsClusOnlyName = "EmcalEmbJetsClusOnly";
TString kEmcalJets025_1GeVName = "EmcalJetsR025_1GeV";
TString kEmcalEmbJets025_1GeVName = "EmcalEmbJetsR025_1GeV";
TString kEmcalJets015_1GeVName = "EmcalJetsR015_1GeV";
TString kEmcalEmbJets015_1GeVName = "EmcalEmbJetsR015_1GeV";
TString kEmcalJets02_1GeVName = "EmcalJetsR02_1GeV";
TString kEmcalEmbJets02_1GeVName = "EmcalEmbJetsR02_1GeV";
TString kEmcalJets03_1GeVName = "EmcalJetsR03_1GeV";
TString kEmcalEmbJets03_1GeVName = "EmcalEmbJetsR03_1GeV";

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 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

/*
// 150 MeV
TF1* sfunc = new TF1("sfunc","[0]*x*x+[1]*x+[2]",-1,100);
  sfunc->SetParameter(2,1.97944);
  sfunc->SetParameter(1,-0.0125036);
  sfunc->SetParameter(0,0.000189736);
*/

/*
// 300 MeV
TF1* sfunc = new TF1("sfunc","[0]*x*x+[1]*x+[2]",-1,100);
  sfunc->SetParameter(2,1.71502);
  sfunc->SetParameter(1,-8.47854e-03);
  sfunc->SetParameter(0,6.28200e-05);
*/
/*
// 300 MeV, had corr = 1.7
TF1* sfunc = new TF1("sfunc","[0]*x*x+[1]*x+[2]",-1,100);
  sfunc->SetParameter(2,1.73158);
  sfunc->SetParameter(1,-5.26873e-03);
  sfunc->SetParameter(0,-3.73229e-05);
*/

// 300 MeV, had corr = 1.3
TF1* sfunc = new TF1("sfunc","[0]*x*x+[1]*x+[2]",-1,100);
  sfunc->SetParameter(2,1.7954);
  sfunc->SetParameter(1, -6.29823e-03);
  sfunc->SetParameter(0,-3.42213e-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);
  sfunc1GeV->SetParameter(2,1.52519);
  sfunc1GeV->SetParameter(1,-1.04701e-02);
  sfunc1GeV->SetParameter(0, 8.99280e-05);
}