//
//  ConfigFemtoAnalysis.C
//

AliFemtoManager*
ConfigFemtoAnalysis()
{
  const double PionMass = 0.13956995;

  //
  // Setup Event Reader
  //
  AliFemtoEventReaderAOD *rdr = new AliFemtoEventReaderAODChain();

  rdr->SetFilterBit(7);
  rdr->SetEPVZERO(kTRUE);
  rdr->SetUseMultiplicity(AliFemtoEventReaderAOD::kCentrality);
  rdr->SetCentralityFlattening(kFALSE);
  rdr->SetReadV0(0);


  AliFemtoManager *manager = new AliFemtoManager();

  manager->SetEventReader(rdr);


  // Analysis

  AliFemtoSimpleAnalysis *analysis = new AliFemtoSimpleAnalysis();
  analysis->SetNumEventsToMix(5);
  analysis->SetVerboseMode(kFALSE);

  // Event Cut

  AliFemtoEventCut *ev_cut = new AliFemtoBasicEventCut();
  ev_cut->AddCutMonitor(
    new AliFemtoCutMonitorEventMult("EvPass"),
    new AliFemtoCutMonitorEventMult("EvFail")
  );
  analysis->SetEventCut(ev_cut);

  // Particle Cut 1

  AliFemtoBasicTrackCut* particle_cut = new AliFemtoBasicTrackCut();
  particle_cut->SetPt(0.5, 2);
  particle_cut->SetRapidity(-0.8, 0.8);
  particle_cut->SetCharge(1);
  particle_cut->AddCutMonitor(
    new AliFemtoCutMonitorParticleYPt("particle1_pass", PionMass),
    new AliFemtoCutMonitorParticleYPt("particle1_fail", PionMass)
  );
  analysis->SetFirstParticleCut(particle_cut);

  // Particle Cut 2

  particle_cut = new AliFemtoBasicTrackCut();
  particle_cut->SetPt(2, 4);
  particle_cut->SetRapidity(-0.8, 0.8);
  particle_cut->SetCharge(1);
  particle_cut->AddCutMonitor(
    new AliFemtoCutMonitorParticleYPt("particle2_pass", PionMass),
    new AliFemtoCutMonitorParticleYPt("particle2_fail", PionMass)
  );
  analysis->SetSecondParticleCut(particle_cut);

  // Pair Cut

  AliFemtoShareQualityPairCut *pair_cut = new AliFemtoShareQualityPairCut();
  pair_cut->SetShareFractionMax(0.2);
  analysis->SetPairCut(pair_cut);

  // Correlation Functions

  analysis->AddCorrFctn(new AliFemtoQinvCorrFctn(
    "Qinv",
    144, 0.0, 2.0
  ));


  analysis->AddCorrFctn(new AliFemtoAvgSepCorrFctn(
    "AvgSep",
    144, 0, 15
  ));

  manager->AddAnalysis(analysis);

  return manager;
}