PWG0/dNdEta/testAnalysis2.C

   1 /* $Id$ */
   2
   3 //
   4 // Script to test the dN/dEta analysis using the dNdEtaAnalysis and
   5 // dNdEtaCorrection classes. Note that there is a cut on the events,
   6 // so the measurement will be biassed.
   7 //
   8 // implementation with TSelector
   9 //
  10
  11 #include "../CreateESDChain.C"
  12
  13 testAnalysis2(Char_t* dataDir, Int_t nRuns=20, Int_t offset=0, Bool_t aMC = kFALSE, Bool_t aDebug = kFALSE)
  14 {
  15   gSystem->Load("libPWG0base");
  16
  17   TChain* chain = CreateESDChainFromDir(dataDir, nRuns, offset);
  18
  19   // ########################################################
  20   // selection of esd tracks
  21   AliESDtrackCuts* esdTrackCuts = new AliESDtrackCuts();
  22   esdTrackCuts->DefineHistograms(1);
  23
  24   esdTrackCuts->SetMinNClustersTPC(50);
  25   esdTrackCuts->SetMaxChi2PerClusterTPC(3.5);
  26   esdTrackCuts->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
  27   esdTrackCuts->SetRequireTPCRefit(kTRUE);
  28
  29   esdTrackCuts->SetMinNsigmaToVertex(3);
  30   esdTrackCuts->SetAcceptKingDaughters(kFALSE);
  31
  32   chain->GetUserInfo()->Add(esdTrackCuts);
  33
  34   if (aMC == kFALSE)
  35   {
  36     dNdEtaCorrection* dNdEtaCorrection = new dNdEtaCorrection();
  37     dNdEtaCorrection->LoadHistograms("correction_map.root","dndeta_correction");
  38     dNdEtaCorrection->RemoveEdges(2, 0, 2);
  39
  40     chain->GetUserInfo()->Add(dNdEtaCorrection);
  41   }
  42
  43   TString selectorName = ((aMC == kFALSE) ? "AlidNdEtaAnalysisESDSelector" : "AlidNdEtaAnalysisMCSelector");
  44   AliLog::SetClassDebugLevel(selectorName, AliLog::kInfo);
  45
  46   selectorName += ".cxx++";
  47   if (aDebug != kFALSE)
  48     selectorName += "g";
  49
  50   TStopwatch timer;
  51   timer.Start();
  52
  53   chain->Process(selectorName);
  54
  55   timer.Stop();
  56   timer.Print();
  57 }