{
/*************************************/
/* 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
};

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 kPhysSelCentral = AliVEvent::kCentral;
//UInt_t kComPhysSel = AliVEvent::kMB + AliVEvent::kSemiCentral + AliVEvent::kCentral + AliVEvent::kEMCEGA; 
//UInt_t kComPhysSel = AliVEvent::kAny;
//UInt_t kComPhysSel = AliVEvent::kEMCEJE;
UInt_t kComPhysSel = AliVEvent::kMB + AliVEvent::kSemiCentral + AliVEvent::kCentral; 
//UInt_t kComPhysSel = AliVEvent::kMB + AliVEvent::kCentral; 
//UInt_t kComPhysSel = AliVEvent::kEMCEJE;
//UInt_t kComPhysSel = AliVEvent::kMB;

// Some trigger definitions for Joel
UInt_t kPhysSelJetHadEMCEJE = AliVEvent::kMB + AliVEvent::kSemiCentral + AliVEvent::kCentral + AliVEvent::kEMCEJE;
UInt_t kPhysSelJetHadEMCEJGA = AliVEvent::kMB + AliVEvent::kSemiCentral + AliVEvent::kCentral + AliVEvent::kEMCEJE + AliVEvent::kEMCEGA;
UInt_t kPhysSelJetHadEMCEGA = AliVEvent::kMB + AliVEvent::kSemiCentral + AliVEvent::kCentral + AliVEvent::kEMCEGA; 
UInt_t kPhysSelJetHadMB = AliVEvent::kMB + AliVEvent::kSemiCentral + AliVEvent::kCentral; 

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

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

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

Bool_t kMakeEmcalTriggers = kFALSE;
const char* kEmcalTriggers = "";
if(kMakeEmcalTriggers) kEmcalTriggers = "EmcalTriggers";

Double_t kTrackEff = 1.0; 

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

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

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

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

Double_t kJetLeadingTrackBias = 0;
Int_t kNcent = 1;

Int_t kNjetResp = 1;
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/AOD149/%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};

Double_t kPtHardBinsScaling[11] = {0.000000E+00, 5.206101E-05, 5.859497E-06, 4.444755E-07, 4.344664E-08, 5.154750E-09, 6.956634E-10, 1.149828E-10, 2.520137E-11, 6.222240E-12, 2.255832E-12}

kPtHardBinsScaling[0] = 0;
kPtHardBinsScaling[1] = 0;
kPtHardBinsScaling[2] = 0;

//Separate train for eacht pt hard bin
for(Int_t i =1; i<11; i++) kPtHardBinsScaling[i] = 0;
kPtHardBinsScaling[8] = 1;

// Uncomment the line below to have a flat distribution of pt hard bins
//for(Int_t i =1; i<11; i++) kPtHardBinsScaling[i] = 1;

Double_t kMinPythiaJetPt = 1;//20; //changed by Marta

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

TString kTracksName = "PicoTracks";
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 kTpcRhoExLJScaledName = "TpcRho_ExLJ_Scaled";

TString kTpcKtJetsName(Form("Jet_KTChargedR020_%s_pT0150_pt_scheme",kTracksName.Data())); 
TString kTpcKtJetsR04Name(Form("Jet_KTChargedR040_%s_pT0150_pt_scheme",kTracksName.Data())); 

TString kTpcJetsR01Name(Form("Jet_AKTChargedR010_%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 kEmcalJets01Name(Form("Jet_AKTFullR010_%s_pT0150_%s_ET0%d_pt_scheme",kTracksName.Data(),kCorrClusName.Data(),TMath::FloorNint(kClusPtCut*1000+0.5)));
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 kEmalJet02UncorrName = "Jet_AKTNeutralR020_EmcCaloClusters_ET0300_pt_scheme";

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

// 1GeV constituent cut
TString kTpcKtJets1GeVName(Form("Jet_KTChargedR020_%s_pT1000_pt_scheme",kTracksName.Data()));
TString kTpcRho1GeVNameExLJ = "TpcRho_1GeV_ExLJ";
TString kTpcRho1GeVSmallName = "TpcRho_1GeV_Small";
TString kEmcalRho1GeVMeth2Name = "TpcRho_1GeV_Small_Scaled";
TString kTpcJets1GeVR02Name(Form("Jet_AKTChargedR020_%s_pT1000_pt_scheme",kTracksName.Data()));
TString kTpcJets1GeVR03Name(Form("Jet_AKTChargedR030_%s_pT1000_pt_scheme",kTracksName.Data()));
TString kEmcalJets1GeVR02Name(Form("Jet_AKTFullR020_%s_pT1000_%s_ET1000_pt_scheme",kTracksName.Data(),kCorrClusName.Data()));
TString kEmcalJets1GeVR03Name(Form("Jet_AKTFullR030_%s_pT1000_%s_ET1000_pt_scheme",kTracksName.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()));

//Jets with E recombination scheme (added by Marta)
TString kTpcKtJetsESchemeName(Form("Jet_KTChargedR020_%s_pT0150_E_scheme",kTracksName.Data())); 
TString kNeutralKtJetsESchemeName(Form("Jet_KTNeutralR020_%s_ET0300_E_scheme",kCorrClusName.Data())); 

//Jets with E recombination scheme (added by Davide&Leticia)
TString kTpcKtJetsESchemeName1GeV(Form("Jet_KTChargedR020_%s_pT1000_E_scheme",kTracksName.Data())); 
TString kTpcKtJetsESchemeNameEmb(Form("Jet_KTChargedR020_%s_pT0150_E_scheme","PicoTracksEmb")); 
TString kTpcAKtJetsESchemeNameEmb(Form("Jet_AKTChargedR040_%s_pT0150_E_scheme","PicoTracksEmb")); 
TString kTpcAKtJetsMCESchemeNameEmb(Form("JetMCOnly_AKTChargedR040_%s_pT0150_E_scheme","PicoTracksEmb")); 
Double_t kPtQG=0.15;
/*************************************/

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

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

// Containers' names
TString kMCTracksName = "MCSelectedParticles";
const char* kUsedMCParticles = "MCParticlesSelected";

TString kMCChargedJetsR02Name(Form("Jet_AKTChargedR020_%s_pT0000_pt_scheme",kMCTracksName.Data()));
TString kMCChargedJetsR03Name(Form("Jet_AKTChargedR030_%s_pT0000_pt_scheme",kMCTracksName.Data()));
TString kMCChargedJetsR04Name(Form("Jet_AKTChargedR040_%s_pT0000_pt_scheme",kMCTracksName.Data()));
TString kMCFullJetsR02Name(Form("Jet_AKTFullR020_%s_pT0000_pt_scheme",kMCTracksName.Data()));
TString kMCFullJetsR03Name(Form("Jet_AKTFullR030_%s_pT0000_pt_scheme",kMCTracksName.Data()));

//Exclude embedding
TString kOrigClusExcludeEmbName = "CaloClustersExcludeEmb";
TString kClusExcludeEmbName = "EmcCaloClustersExcludeEmb";
TString kCorrClusExcludeEmbName = "CaloClustersCorrExcludeEmb";
TString kEmcalClusExcludeEmbName = "EmcalClusters_";
kEmcalClusExcludeEmbName += kClusExcludeEmbName;
TString kEmcalTracksExcludeEmbName = "EmcalTracks_";
kEmcalTracksExcludeEmbName += "ExcludeEmb";
kEmcalTracksExcludeEmbName += kTracksName;

TString kTpcJetsR02ExcludeEmbName(Form("JetExcludeEmb_AKTChargedR020_%s_pT0150_pt_scheme",kTracksName.Data()));
TString kTpcJetsR03ExcludeEmbName(Form("JetExcludeEmb_AKTChargedR030_%s_pT0150_pt_scheme",kTracksName.Data()));
TString kTpcJetsR04ExcludeEmbName(Form("JetExcludeEmb_AKTChargedR040_%s_pT0150_pt_scheme",kTracksName.Data()));
TString kEmcalJetsR02ExcludeEmbName(Form("JetExcludeEmb_AKTFullR020_%s_pT0150_%s_ET0300_pt_scheme",kTracksName.Data(),kCorrClusExcludeEmbName.Data()));
TString kEmcalJetsR03ExcludeEmbName(Form("JetExcludeEmb_AKTFullR030_%s_pT0150_%s_ET0300_pt_scheme",kTracksName.Data(),kCorrClusExcludeEmbName.Data()));

TString kTpcKtJetsExcludeEmbName(Form("JetExcludeEmb_KTChargedR020_%s_pT0150_pt_scheme",kTracksName.Data()));

TString kTpcRhoExLJExcludeEmbName = "TpcRho_ExLJ_ExcludeEmb";
TString kTpcRhoSmallExcludeEmbName = "TpcRho_Small_ExcludeEmb";
TString kEmcalRhoMeth2ExcludeEmbName = "TpcRho_Small_ExcludeEmb_Scaled";

//Embedding only
TString kOrigClusEmbOnlyName = "CaloClustersEmbOnly";
TString kClusEmbOnlyName = "EmcCaloClustersEmbOnly";
TString kCorrClusEmbOnlyName = "CaloClustersCorrEmbOnly";
TString kEmcalClusEmbOnlyName = "EmcalClusters_";
kEmcalClusEmbOnlyName += kClusEmbOnlyName;
TString kEmcalTracksEmbOnlyName = "EmcalTracks_";
kEmcalTracksEmbOnlyName += "EmbOnly";
kEmcalTracksEmbOnlyName += kTracksName;

TString kTpcJetsR02EmbOnlyName(Form("JetEmbOnly_AKTChargedR020_%s_pT0150_pt_scheme",kTracksName.Data()));
TString kTpcJetsR03EmbOnlyName(Form("JetEmbOnly_AKTChargedR030_%s_pT0150_pt_scheme",kTracksName.Data()));
TString kTpcJetsR04EmbOnlyName(Form("JetEmbOnly_AKTChargedR040_%s_pT0150_pt_scheme",kTracksName.Data()));
TString kEmcalJetsR02EmbOnlyName(Form("JetEmbOnly_AKTFullR020_%s_pT0150_%s_ET0300_pt_scheme",kTracksName.Data(),kCorrClusEmbOnlyName.Data()));
TString kEmcalJetsR03EmbOnlyName(Form("JetEmbOnly_AKTFullR030_%s_pT0150_%s_ET0300_pt_scheme",kTracksName.Data(),kCorrClusEmbOnlyName.Data()));

//PbPb only
TString kTpcJetsR02PbPbOnlyName(Form("JetPbPbOnly_AKTChargedR020_%s_pT0150_pt_scheme",kTracksName.Data()));

// Other settings
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                            */
/*************************************/

/*************************************/
/* Marta                            */
/*************************************/
TString kRhoMassName = "RhoMass";
Float_t kChJetPhiMin = -1;
Float_t kChJetPhiMax = -1;
Float_t kChJetPhiOffset = 0.;

TF1* srhomfunc = new TF1("srhomfunc","[0]*x*x+[1]*x+[2]",-1,100);
srhomfunc->SetParameter(2, 1.68354);
srhomfunc->SetParameter(1, -2.86991e-03);
srhomfunc->SetParameter(0, -1.49981e-05);

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

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

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

if (kHadCorr > 1.9) { // had corr = 2
  if (kClusterizerType==AliEMCALRecParam::kClusterizerv2) { // v2 clusterizer
    // 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);

    // fit from test train n. 457
    sfunc1GeV->SetParameter(2, 1.63498);
    sfunc1GeV->SetParameter(1, -0.0208189);
    sfunc1GeV->SetParameter(0, 0.000303435);
  } // 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 (0) { // 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

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


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

/*************************************/
const char* kPeriod = "lhc11h";
const char* kTrackCuts = "Hybrid_LHC11h";
const char* kDatatype = "AOD";
const char* pass = "pass2";
TString kOrigClusName = "caloClusters";
TString kEMCALCellsName="emcalCells";
TString kInTracksName = "AODFilterTracks";
}