]>
Commit | Line | Data |
---|---|---|
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) | |
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 | ||
45 | AliLog::SetClassDebugLevel(selectorName, AliLog::kInfo); | |
46 | ||
47 | TStopwatch timer; | |
48 | timer.Start(); | |
49 | ||
50 | chain->Process(selectorName + ".cxx+"); | |
51 | ||
52 | timer.Stop(); | |
53 | timer.Print(); | |
54 | } |