{
const char* kPeriod = "lhc11h";
const char* kTrackCuts = "lhc11h";
const char* kTrackCutsRosi = "lhc10hold";
const char* kDatatype = "ESD";
const char* kUsedMCParticles = "MCParticlesSelected";

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

Bool_t kMakeTrigger = kFALSE;

Int_t kHistoType = 1;

Bool_t kDoTender = kTRUE;

Double_t kClusPtCut = 0.30;
Double_t kTrackPtCut = 0.15;
Double_t kMinPtHadCorr = 0.15;

Double_t kJetLeadingTrackBias = 5;
Int_t kNcent = 3;

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/AOD141/%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;
//Double_t kPtHardBinsScaling[11] = {0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1};

Double_t kMinPythiaJetPt = 0;

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

Int_t kLeadHadType = 0; // 0 = charged, 1 = neutral, 2 = both

TString kTracksName = "PicoTracks";
TString kInTracksName = "HybridTracks";
TString kClusName = "EmcCaloClusters";
TString kOrigClusName = "CaloClusters";
TString kEMCALCellsName="EMCALCells";
TString kCorrClusName = "CaloClustersCorr";
TString kEmcalTracksName = "EmcalTracks_";
kEmcalTracksName += kTracksName;
TString kEmcalClusName = "EmcalClusters_";
kEmcalClusName += kClusName;

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_pt_scheme",kTracksName.Data())); 
TString kPythiaEmcalKtJetsName(Form("JetEmcal_KTChargedR020_%s_pT0150_pt_scheme",kTracksName.Data()));
TString kTpcKtJets1GeVName(Form("Jet_KTChargedR020_%s_pT1000_pt_scheme",kTracksName.Data()));
TString kTpcJetsR02Name(Form("Jet_AKTChargedR020_%s_pT0150_pt_scheme",kTracksName.Data()));
TString kTpcJetsR03Name(Form("Jet_AKTChargedR030_%s_pT0150_pt_scheme",kTracksName.Data()));
TString kTpcJetsR04Name(Form("Jet_AKTChargedR040_%s_pT0150_pt_scheme",kTracksName.Data()));
TString kTpcEmbJetsR02Name(Form("Jet_AKTChargedR020_%s_pT0150_pt_scheme",kTrackEmbName.Data()));
TString kTpcEmbJetsR03Name(Form("Jet_AKTChargedR030_%s_pT0150_pt_scheme",kTrackEmbName.Data()));
TString kTpcEmbJetsR04Name(Form("Jet_AKTChargedR040_%s_pT0150_pt_scheme",kTrackEmbName.Data()));

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

TString kEmcalJets015Name(Form("Jet_AKTFullR015_%s_pT0150_%s_ET0%d_pt_scheme",kTracksName.Data(),kCorrClusName.Data(),TMath::FloorNint(kClusPtCut*1000+0.5)));
TString kEmcalJets02Name(Form("Jet_AKTFullR020_%s_pT0150_%s_ET0%d_pt_scheme",kTracksName.Data(),kCorrClusName.Data(),TMath::FloorNint(kClusPtCut*1000+0.5)));
TString kEmcalJets025Name(Form("Jet_AKTFullR025_%s_pT0150_%s_ET0%d_pt_scheme",kTracksName.Data(),kCorrClusName.Data(),TMath::FloorNint(kClusPtCut*1000+0.5)));
TString kEmcalJets03Name(Form("Jet_AKTFullR030_%s_pT0150_%s_ET0%d_pt_scheme",kTracksName.Data(),kCorrClusName.Data(),TMath::FloorNint(kClusPtCut*1000+0.5)));

TString kEmcalEmbJets015Name(Form("Jet_AKTFullR015_%s_pT0150_%s_ET0%d_pt_scheme",kTrackEmcalEmbName.Data(),kCorrClusName.Data(),TMath::FloorNint(kClusPtCut*1000+0.5)));
TString kEmcalEmbJets02Name(Form("Jet_AKTFullR020_%s_pT0150_%s_ET0%d_pt_scheme",kTrackEmcalEmbName.Data(),kCorrClusName.Data(),TMath::FloorNint(kClusPtCut*1000+0.5)));
TString kEmcalEmbJets025Name(Form("Jet_AKTFullR025_%s_pT0150_%s_ET0%d_pt_scheme",kTrackEmcalEmbName.Data(),kCorrClusName.Data(),TMath::FloorNint(kClusPtCut*1000+0.5)));
TString kEmcalEmbJets03Name(Form("Jet_AKTFullR030_%s_pT0150_%s_ET0%d_pt_scheme",kTrackEmcalEmbName.Data(),kCorrClusName.Data(),TMath::FloorNint(kClusPtCut*1000+0.5)));

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

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

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

// Toy model Single Particle embedding
TString kTpcEmbSpectrumJets02Name(Form("Jet_AKTChargedR020_%s_pT0150_pt_scheme",kTrackEmbSpectrumName.Data()));
TString kTpcEmbSpectrumJets03Name(Form("Jet_AKTChargedR030_%s_pT0150_pt_scheme",kTrackEmbSpectrumName.Data()));
TString kTpcEmbSpectrumJets04Name(Form("Jet_AKTChargedR040_%s_pT0150_pt_scheme",kTrackEmbSpectrumName.Data()));

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

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

TString kTpcKtPythiaEmbJetsName(Form("JetPythiaEmb_KTChargedR020_%s_pT0150_pt_scheme",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;
Double_t kMaxGeoDistance02 = 1.2;
Double_t kMaxGeoDistance03 = 1.2;
Double_t kMaxGeoDistance04 = 1.2;
/*
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;
}
*/

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 kEMCtimeMin = -50e-6;
Double_t kEMCtimeMax = 100e-6;
Double_t kEMCtimeCut = 75e-6;

Double_t kPi = TMath::Pi();
Double_t kGhostArea = 0.005;
Double_t kKtJetRadius = 0.2;
Double_t kJetAreaCut = 0.557;
Double_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.0;
Bool_t kDoEmbedding = kFALSE;
Bool_t kDoReclusterize = kTRUE;
UInt_t kClusterizerType = AliEMCALRecParam::kClusterizerv2;

if (kDoEmbedding) kDoReclusterize = kTRUE;

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 kPhysSel = AliVEvent::kEMCEGA;

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

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 (0) { // 430 propagation (AOD115 & lts)
   if (prePhysSel == AliVEvent::kEMCEGA) {
      // 300 MeV, had corr = 2, GA trigger, v2 cluster
      sfunc->SetParameter(2,1.9007);
      sfunc->SetParameter(1,-0.0102328);
      sfunc->SetParameter(0,0.000128441);
   } 
   else 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 (kEMCtimeMin>-200e-9) { // 440 propagation (AOD142, AOD145 & ESD) & cluster time cuts
      // 300 MeV, had corr = 2
      // fit from train n. 370, run 169838, histogram ScalevsCent in range 0-50
      sfunc->SetParameter(2, 1.74475);
      sfunc->SetParameter(1, -0.0112865);
      sfunc->SetParameter(0, 0.000113246);
  }
  else { // 440 propagation (AOD142, AOD145 & ESD) & no time cuts
      // 300 MeV, had corr = 2
      sfunc->SetParameter(2,1.76458);
      sfunc->SetParameter(1,-0.0111656);
      sfunc->SetParameter(0,0.000107296);
  }

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

const char* pass = "pass2";

// Trigger

TString kTriggerClusName = "L1TriggerClusters";
TString kTriggerClusFastORName = "L1TriggerClustersFastOR";
TF1 *kTriggerThreshold = new TF1("eth","[0] + [1]*x + [2]*x*x",0,100);
//kTriggerThreshold->SetParameter(0,137.8);
//kTriggerThreshold->SetParameter(1,-1.28);
//kTriggerThreshold->SetParameter(0,79.9);
//kTriggerThreshold->SetParameter(1,-0.701);
kTriggerThreshold->SetParameter(0,5.);
kTriggerThreshold->SetParameter(1,0.0031);
kTriggerThreshold->SetParameter(2,0.);

//JetFlow

AliAnalysisTaskSE* kRhoVnMod[10];

}