{
//boost::v1_43_0 cgal::v3.6 fastjet::v2.4.2
//boost::v1_53_0 cgal::v4.4 fastjet::v3.0.6_1.012
/*************************************/
/* Some definitions                  */
/*************************************/

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

//UInt_t kMBPhysSel = AliVEvent::kMB;
//UInt_t kNoPhysSel = AliVEvent::kAny;

Double_t kPi = TMath::Pi();

UInt_t kTpcAna = 0;
UInt_t kEmcalAna = 1;
UInt_t kUserAna = 2;

/*************************************/

/*************************************/
/* Trigger & Centralities            */
/*************************************/

Int_t nCentBins = 4;
Double_t minCent = 0;
Double_t maxCent = 100;

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;

/*************************************/

/*************************************/
/* Analysis Cuts & Settings          */
/*************************************/

Bool_t kDoTender = kTRUE;
Double_t kClusPtCut = 0.30;
Double_t kTrackPtCut = 0.15;
Double_t kMinPtHadCorr = 0.15;

Double_t kPropDist = 440;
Double_t kHadCorr = 2.0;
Bool_t kDoEmbedding = kFALSE;
Bool_t kDoReclusterize = kTRUE;
if (kDoEmbedding) kDoReclusterize = kTRUE;
UInt_t kClusterizerType = AliEMCALRecParam::kClusterizerv2;

Double_t kEMCtimeMin = -50e-6;
Double_t kEMCtimeMax = 100e-6;
Double_t kEMCtimeCut = 75e-6;

Double_t kGhostArea = 0.005;
Double_t kKtJetRadius = 0.2;
Double_t kJetAreaCut = 0.557;

Double_t kJetPtCut = 1;
Double_t kJetBiasTrack = 5;
Double_t kJetBiasClus = 1000;

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

Int_t kLeadHadType = 0; // 0 = charged, 1 = neutral, 2 = both
Int_t kHistoType = 1; // 1 = THnSparse, 0 = TH1/TH2/TH3

/*************************************/

/*************************************/
/* Scale Factors                     */
/*************************************/

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

if (kHadCorr > 1.9) { // had corr = 2
  if (kClusterizerType==AliEMCALRecParam::kClusterizerv2) { // v2 clusterizer
    // 300 MeV, had corr = 2
    // used for train n. 374
    sfunc->SetParameter(2, 1.76458);
    sfunc->SetParameter(1, -0.0111656);
    sfunc->SetParameter(0, 0.000107296);
    // fit from full train n. 374
    //sfunc->SetParameter(2, 1.78067);
    //sfunc->SetParameter(1, -0.0116032);
    //sfunc->SetParameter(0, 0.000114951);
  } // end v2 clusterizer
  else { // 3x3 clusterizer
    // used for train n. 408 (fit test train n.408)
    //sfunc->SetParameter(2, 1.85043);
    //sfunc->SetParameter(1, -0.0116461);
    //sfunc->SetParameter(0, 9.65308e-05);
    // fit from full train n.408
    sfunc->SetParameter(2, 1.8528);
    sfunc->SetParameter(1, -0.0121015);
    sfunc->SetParameter(0, 0.000109308);
  } // end 3x3 clusterizer
} // end had corr = 2.0
else if (kHadCorr > 1.6) { // had corr = 1.7
  // fit from test run of train n. 387
  sfunc->SetParameter(2, 1.8171);
  sfunc->SetParameter(1, -0.0150141);
  sfunc->SetParameter(0, 0.000193964);
} // end had corr = 1.7
else if (kHadCorr > 1.2) { // had corr = 1.3
  // fit from test run of train n. 389
  sfunc->SetParameter(2, 1.90496);
  sfunc->SetParameter(1, -0.016978);
  sfunc->SetParameter(0, 0.000215038);
} // end had corr = 1.3
else {   // had corr = 0
  if (kClusterizerType==AliEMCALRecParam::kClusterizerv2) { // v2 clusterizer
    // 300 MeV, had corr = 2
    // fit from full train n. 401
    sfunc->SetParameter(2, 2.0418);
    sfunc->SetParameter(1, -0.0157904);
    sfunc->SetParameter(0, 0.000148058);
  } // end v2 clusterizer
  else { // 3x3 clusterizer
    // fit from full train n. 407
    sfunc->SetParameter(2, 2.13631);
    sfunc->SetParameter(1, -0.0165245);
    sfunc->SetParameter(0, 0.000142767);
  } // end 3x3 clusterizer
} // end had corr = 0

if (prePhysSel == AliVEvent::kEMCEGA) { // GA trigger
  if (kHadCorr > 1.9) { // had corr = 2
    // old
    //sfunc->SetParameter(2,1.92213);
    //sfunc->SetParameter(1,-0.0175260);
    //sfunc->SetParameter(0,0.000152385);

    sfunc->SetParameter(2, 1.81208);
    sfunc->SetParameter(1, -0.0105506);
    sfunc->SetParameter(0, 0.000145519);
  } // end had corr = 2.0
} //end ga trigger

/*************************************/

/*************************************/
/* Containers' names                 */
/*************************************/

TString kTracksName = "PicoTracks";
TString kInTracksName = "HybridTracks";
TString kClusName = "EmcCaloClusters";
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 kTpcRhoName = "TpcRho";
TString kTpcRhoNameExLJ = "TpcRho_ExLJ";
TString kTpcRhoSmallName = "TpcRho_Small";
TString kEmcalRhoMeth2Name = "TpcRho_Small_Scaled";

TString kTpcKtJetsName(Form("Jet_KTChargedR020_%s_pT0150_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 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)));

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

/*************************************/

/*************************************/
/* Trigger Tasks                     */
/*************************************/

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

/*************************************/

/*************************************/
/* PYTHIA Embedding                  */
/*************************************/

Double_t kJetLeadingTrackBias = 5;
Int_t kNcent = 3;

Int_t kNjetResp = 1;
const char* kUsedMCParticles = "MCParticlesSelected";
if (kJetLeadingTrackBias > 1)
  kNjetResp = 2;

if (kJetLeadingTrackBias > 5)
  kNjetResp = 3;

kNjetResp *= kNcent;
Bool_t kMakeTrigger = kFALSE;
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 kMinPythiaJetPt = 0;

// Containers' names

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

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()));
TString kPythiaEmcalKtJetsName(Form("JetEmcal_KTChargedR020_%s_pT0150_pt_scheme",kTracksName.Data()));

TString kTpcPythiaEmbRhoNameExLJ = "TpcPythiaEmbRho_ExLJ";
TString kTpcPythiaEmbRhoSmallName = "TpcPythiaEmbRho_Small";
TString kEmcalPythiaEmbRhoMeth2Name = "TpcPythiaEmbRho_Small_Scaled";

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;

Double_t kPartLevPtCut = 0;

/*************************************/

/*************************************/
/* Redmer                            */
/*************************************/

AliAnalysisTaskSE* kRhoVnMod[10];

/*************************************/

/*************************************/
/* Darius                            */
/*************************************/

std::cout << "_____________________________________________________________________" << std::endl << "Printing all global variables" << std::endl << "_____________________________________________________________________" << std::endl;
gSystem->Exec("printenv");
std::cout << "_____________________________________________________________________" << std::endl;

/*************************************/
std::cout << "Print enviroment variables..." << std::endl;
std::cout << "ALIEN_JDL_LPMPRODUCTIONTYPE: " << gSystem->Getenv("ALIEN_JDL_LPMPRODUCTIONTYPE") << std::endl;
}