{
enum {
  kKT=1<<0,
  kAKT=1<<1,
  kFullJet=1<<2,
  kChargedJet=1<<3,
  kNeutralJet=1<<4,
  kR020Jet=1<<5,
  kR030Jet=1<<6,
  kR040Jet=1<<7,
  kRX1Jet=1<<8,  // user defined radii, use SetRadius(Double_t)
  kRX2Jet=1<<9,
  kRX3Jet=1<<10
};

Float_t kJetLeadingTrackBias = 0;
Int_t kNcent = 1;

Int_t kNjetResp = 1;

if (kJetLeadingTrackBias > 1)
  kNjetResp = 2;

if (kJetLeadingTrackBias > 5)
  kNjetResp = 3;

kNjetResp *= kNcent;

if (1) {
  UInt_t kPythiaR020Charged = 0; 
  UInt_t kPythiaR030Charged = 0;
  UInt_t kPythiaR040Charged = 0;

  UInt_t kPythiaR020Full = 0;
  UInt_t kPythiaR030Full = 0;
}
else {
  UInt_t kPythiaR020Charged = kAKT|kR020Jet|kChargedJet;
  UInt_t kPythiaR030Charged = kAKT|kR030Jet|kChargedJet;
  UInt_t kPythiaR040Charged = kAKT|kR040Jet|kChargedJet;

  UInt_t kPythiaR020Full = kAKT|kR020Jet|kFullJet;
  UInt_t kPythiaR030Full = kAKT|kR030Jet|kFullJet;
}

TString kPYTHIAPath = "alien:///alice/sim/2012/LHC12a15e_fix/%d/%d/AOD129/%04d/AliAOD.root";
Int_t nAODFiles = 140;
Int_t kNpTHardBins = 11;
Double_t kPtHardBinsScaling[11] = {0, 0, 3.27317e-05, 2.57606e-06, 2.5248e-07, 2.92483e-08, 4.15631e-09, 6.6079e-10, 1.49042e-10, 3.5571e-11, 1.29474e-11};
kPtHardBinsScaling[0] = 0;
kPtHardBinsScaling[1] = 0;
kPtHardBinsScaling[2] = 0;
kPtHardBinsScaling[3] = 0;

Double_t kMinPythiaJetPt = 10;

TString kMCTracksName = "MCSelectedParticles";
TString kClusMCOnlyName = "CaloClustersMCOnly";
TString kCorrClusMCOnlyName = "CaloClustersMCOnlyCorr";
TString kEmcalClusMCOnlyName = "EmcalClustersMCOnly";

TString kTracksName = "PicoTracks";
TString kInTracksName = "tracks";
TString kClusName = "caloClusters";
TString kEMCALCellsName="emcalCells";
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 kTrackEmbSpectrumName = "PicoTracksEmbeddedSpectrum";

TString kTpcKtJetsName(Form("Jet_KTChargedR020_%s_pT0150",kTracksName.Data()));
TString kPythiaEmcalKtJetsName(Form("JetEmcal_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()));

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()));

// Toy model Single Particle embedding
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()));

//Toy model Pythia embedding
TString kMCChargedJets02Name(Form("Jet_AKTChargedR020_%s_pT0000",kMCTracksName.Data()));
TString kMCChargedJets03Name(Form("Jet_AKTChargedR030_%s_pT0000",kMCTracksName.Data()));
TString kMCChargedJets04Name(Form("Jet_AKTChargedR040_%s_pT0000",kMCTracksName.Data()));
TString kMCFullJets02Name(Form("Jet_AKTFullR020_%s_pT0000",kMCTracksName.Data()));
TString kMCFullJets03Name(Form("Jet_AKTFullR030_%s_pT0000",kMCTracksName.Data()));

TString kTpcJetsMConlyR02Name(Form("JetMConly_AKTChargedR020_%s_pT0150",kTracksName.Data()));
TString kTpcJetsMConlyR03Name(Form("JetMConly_AKTChargedR030_%s_pT0150",kTracksName.Data()));
TString kTpcJetsMConlyR04Name(Form("JetMConly_AKTChargedR040_%s_pT0150",kTracksName.Data()));
TString kEmcalJetsMConly02Name(Form("JetMConly_AKTFullR020_%s_pT0150_%s_ET0300",kTracksName.Data(),kCorrClusMCOnlyName.Data()));
TString kEmcalJetsMConly03Name(Form("JetMConly_AKTFullR030_%s_pT0150_%s_ET0300",kTracksName.Data(),kCorrClusMCOnlyName.Data()));

TString kTpcKtPythiaEmbJetsName(Form("JetPythiaEmb_KTChargedR020_%s_pT0150",kTracksName.Data()));

UInt_t kMatching = 1; //1=geometrical, 2=MClabel, 3=same collections
Double_t kMaxDistance02a = 1;
Double_t kMaxDistance03a = 1;
Double_t kMaxDistance04a = 1;
Double_t kMaxDistance02b = 1;
Double_t kMaxDistance03b = 1;
Double_t kMaxDistance04b = 1;
if (kMatching == 1) {
  kMaxDistance02a = 0.125;
  kMaxDistance03a = 0.18;
  kMaxDistance04a = 0.25;
  kMaxDistance02b = 0.125;
  kMaxDistance03b = 0.18;
  kMaxDistance04b = 0.25;
}
else if (kMatching == 3) {
  kMaxDistance02a = 1;
  kMaxDistance03a = 1;
  kMaxDistance04a = 1;
  kMaxDistance02b = 0.5;
  kMaxDistance03b = 0.5;
  kMaxDistance04b = 0.5;
}

Double_t kMaxGeoDistance02 = 0.125;
Double_t kMaxGeoDistance03 = 0.18;
Double_t kMaxGeoDistance04 = 0.25;

TString kTpcRhoName = "TpcRho";
TString kTpcRhoNameExLJ = "TpcRho_ExLJ";
TString kTpcRhoSmallName = "TpcRho_Small";
TString kTpcRho1GeVNameExLJ = "TpcRho_1GeV";
TString kTpcRhoSmall1GeVNameExLJ = "TpcRho_Small_1GeV";

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

//Toy model Pythia embedding
TString kTpcPythiaEmbRhoNameExLJ = "TpcPythiaEmbRho_ExLJ";
TString kTpcPythiaEmbRhoSmallName = "TpcPythiaEmbRho_Small";
TString kEmcalPythiaEmbRhoMeth2Name = "TpcPythiaEmbRho_Small_Scaled";

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 kEmcalMinEta = -0.7;
Double_t kEmcalMaxEta = 0.7;
Double_t kEmcalMinPhi = 80*TMath::DegToRad();
Double_t kEmcalMaxPhi = 180*TMath::DegToRad();

Double_t kPi = TMath::Pi();
Double_t kGhostArea = 0.005;
Double_t kKtJetRadius = 0.2;
Double_t kJetAreaCut = 0.557;
Double_t kTrackPtCut = 0.15;
Double_t kClusPtCut = 0.30;
Float_t kPartLevPtCut = 0;
Double_t kJetPtCut = 1;
Double_t kJetBiasTrack = 5;
Double_t kJetBiasClus = 1000;
UInt_t kEmcalAna = 1;
UInt_t kTpcAna = 0;
UInt_t kUserAna = 2;
Double_t kHadCorr = 2;
Bool_t kDoReclusterize = kTRUE;
UInt_t kClusterizerType = AliEMCALRecParam::kClusterizerv2;

Double_t kPropDist = 440;

UInt_t kPhysSel = 1<<31; //AliEmcalPhysicsSelection::kEmcalOk
//UInt_t kPhysSel = AliVEvent::kAnyINT + AliVEvent::kSemiCentral + AliVEvent::kCentral;
//UInt_t kPhysSel = AliVEvent::kAny;
//UInt_t kPhysSel = AliVEvent::kMB;

UInt_t prePhysSel = AliVEvent::kAnyINT + AliVEvent::kSemiCentral + AliVEvent::kCentral; 
//UInt_t prePhysSel = AliVEvent::kEMCEGA;
//UInt_t prePhysSel = AliVEvent::kEMCEJE;

UInt_t kMBPhysSel = AliVEvent::kMB;

UInt_t kNoPhysSel = AliVEvent::kAny;

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

if (kHadCorr > 1.9) {
  if (kClusterizerType==AliEMCALRecParam::kClusterizerv2) {
    // 300 MeV, had corr = 2
    sfunc->SetParameter(2,1.8078);
    sfunc->SetParameter(1,-0.0115442);
    sfunc->SetParameter(0,0.000109722);
  }
  else {
    // 3x3, 300 MeV, had corr = 2
    sfunc->SetParameter(2, 2.02415);
    sfunc->SetParameter(1, -8.13464e-03);
    sfunc->SetParameter(0,2.28091e-05);
  }
}
else if (kHadCorr > 1.6) {
  // 300 MeV, had corr = 1.7
  sfunc->SetParameter(2, 1.79148);
  sfunc->SetParameter(1, -1.13795e-02);
  sfunc->SetParameter(0, 1.14528e-04);
}
else if (kHadCorr > 1.2) {
  // 300 MeV, had corr = 1.3
  sfunc->SetParameter(2,1.90655);
  sfunc->SetParameter(1,-0.0159274);
  sfunc->SetParameter(0,0.0001706891);
}
else {
  // 300 MeV, had corr = 0
  sfunc->SetParameter(2,2.01265);
  sfunc->SetParameter(1,-0.0128864);
  sfunc->SetParameter(0,0.000100894);
}

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

}