standalone Makefile, preliminary (towards PROOF par archive)

[u/mrichter/AliRoot.git] / PWG0 / dNdEta / testAnalysis2.C

//
// Script to test the dN/dEta analysis using the dNdEtaAnalysis and
// dNdEtaCorrection classes. Note that there is a cut on the events,
// so the measurement will be biassed.
//
// implementation with TSelector
//

#include "../CreateESDChain.C"

testAnalysis2(Char_t* dataDir, Int_t nRuns=20, Bool_t aMC = kFALSE)
{
  gSystem->Load("libPWG0base");

  TChain* chain = CreateESDChain(dataDir, nRuns);

  // ########################################################
  // selection of esd tracks
  AliESDtrackCuts* esdTrackCuts = new AliESDtrackCuts();
  esdTrackCuts->DefineHistograms(1);

  esdTrackCuts->SetMinNClustersTPC(50);
  esdTrackCuts->SetMaxChi2PerClusterTPC(3.5);
  esdTrackCuts->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
  esdTrackCuts->SetRequireTPCRefit(kTRUE);

  esdTrackCuts->SetMinNsigmaToVertex(3);
  esdTrackCuts->SetAcceptKingDaughters(kFALSE);

  chain->GetUserInfo()->Add(esdTrackCuts);

  dNdEtaCorrection* dNdEtaCorrection = new dNdEtaCorrection();
  dNdEtaCorrection->LoadHistograms("correction_map.root","dndeta_correction");
  dNdEtaCorrection->RemoveEdges(2,0,2);

  chain->GetUserInfo()->Add(dNdEtaCorrection);

  TString selectorName = ((aMC == kFALSE) ? "AlidNdEtaAnalysisSelector" : "AlidNdEtaAnalysisMCSelector");

  AliLog::SetClassDebugLevel(selectorName, AliLog::kInfo);

  chain->Process(selectorName + ".cxx+");
}