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.
8 // implementation with TSelector
11 #include "../CreateESDChain.C"
12 #include "../PWG0Helper.C"
14 void testAnalysis2(Char_t* data, Int_t nRuns=20, Int_t offset=0, Bool_t aMC = kFALSE, Bool_t aDebug = kFALSE, Bool_t aProof = kFALSE, const char* correctionMapFile = "correction_map.root", const char* correctionMapFolder = "dndeta_correction", const char* option = "", const char* proofServer = "jgrosseo@lxb6046")
18 connectProof(proofServer);
19 gProof->AddInput(new TParameter<long>("PROOF_MaxSlavesPerNode", (long)2));
20 gProof->AddInput(new TNamed("PROOF_Packetizer", "TAdaptivePacketizer"));
23 TString libraries("libEG;libGeom;libESD;libPWG0base");
24 TString packages("adaptivepacketizer;PWG0base");
26 if (!prepareQuery(libraries, packages, kTRUE))
29 //TODO somehow prevent loading several times
30 gROOT->ProcessLine(".L CreateCuts.C");
31 gROOT->ProcessLine(".L drawPlots.C");
33 TChain* chain = CreateESDChain(data, nRuns, offset);
39 // selection of esd tracks
40 AliESDtrackCuts* esdTrackCuts = CreateTrackCuts();
43 printf("ERROR: esdTrackCuts could not be created\n");
47 inputList.Add(esdTrackCuts);
50 TString selectorName = ((aMC == kFALSE) ? "AlidNdEtaAnalysisESDSelector" : "AlidNdEtaAnalysisMCSelector");
51 AliLog::SetClassDebugLevel(selectorName, AliLog::kInfo);
53 selectorName += ".cxx+";
58 Int_t result = executeQuery(chain, &inputList, selectorName, option);
64 dNdEtaAnalysis* fdNdEtaAnalysis = new dNdEtaAnalysis("dndeta", "dndeta");
66 TFile* file = TFile::Open("analysis_mc.root");
70 cout << "Error. File not found" << endl;
73 fdNdEtaAnalysis->LoadHistograms();
74 fdNdEtaAnalysis->DrawHistograms(kTRUE);
77 FinishAnalysisAll("analysis_esd_raw.root", "analysis_esd.root", correctionMapFile, correctionMapFolder);
81 void FinishAnalysisAll(const char* dataInput = "analysis_esd_raw.root", const char* dataOutput = "analysis_esd.root", const char* correctionMapFile = "correction_map.root", const char* correctionMapFolder = "dndeta_correction")
83 gSystem->Load("libPWG0base");
85 AlidNdEtaCorrection* dNdEtaCorrection = new AlidNdEtaCorrection(correctionMapFolder, correctionMapFolder);
86 TFile::Open(correctionMapFile);
87 dNdEtaCorrection->LoadHistograms();
89 TFile* file = TFile::Open(dataInput);
93 cout << "Error. File not found" << endl;
97 dNdEtaAnalysis* fdNdEtaAnalysis = new dNdEtaAnalysis("dndeta", "dndeta");
98 fdNdEtaAnalysis->LoadHistograms("dndeta");
99 fdNdEtaAnalysis->Finish(dNdEtaCorrection, 0.3, AlidNdEtaCorrection::kINEL);
100 fdNdEtaAnalysis->DrawHistograms(kTRUE);
101 TFile* file2 = TFile::Open(dataOutput, "RECREATE");
102 fdNdEtaAnalysis->SaveHistograms();
105 fdNdEtaAnalysis = new dNdEtaAnalysis("dndetaTr", "dndetaTr");
106 fdNdEtaAnalysis->LoadHistograms("dndeta");
107 fdNdEtaAnalysis->Finish(dNdEtaCorrection, 0.3, AlidNdEtaCorrection::kVertexReco);
108 fdNdEtaAnalysis->DrawHistograms(kTRUE);
110 fdNdEtaAnalysis->SaveHistograms();
113 fdNdEtaAnalysis = new dNdEtaAnalysis("dndetaTrVtx", "dndetaTrVtx");
114 fdNdEtaAnalysis->LoadHistograms("dndeta");
115 fdNdEtaAnalysis->Finish(dNdEtaCorrection, 0.3, AlidNdEtaCorrection::kTrack2Particle);
116 fdNdEtaAnalysis->DrawHistograms(kTRUE);
118 fdNdEtaAnalysis->SaveHistograms();
121 void FinishAnalysis(const char* analysisFile = "analysis_esd.root", const char* analysisDir = "dndeta", const char* correctionMapFile = "correction_map.root", const char* correctionMapFolder = "dndeta_correction", Bool_t useUncorrected = kFALSE, Bool_t simple = kFALSE)
123 gSystem->Load("libPWG0base");
125 TFile* file = TFile::Open(analysisFile);
127 dNdEtaAnalysis* fdNdEtaAnalysis = new dNdEtaAnalysis(analysisDir, analysisDir);
128 fdNdEtaAnalysis->LoadHistograms();
130 if (correctionMapFile)
132 AlidNdEtaCorrection* dNdEtaCorrection = new AlidNdEtaCorrection(correctionMapFolder, correctionMapFolder);
133 TFile::Open(correctionMapFile);
134 dNdEtaCorrection->LoadHistograms();
136 fdNdEtaAnalysis->Finish(dNdEtaCorrection, 0.3);
139 fdNdEtaAnalysis->DrawHistograms(simple);
141 TH1* hist = fdNdEtaAnalysis->GetdNdEtaHistogram(1);
142 Int_t binLeft = hist->GetXaxis()->FindBin(-0.5);
143 Int_t binRight = hist->GetXaxis()->FindBin(0.5);
144 Float_t value1 = hist->Integral(binLeft, binRight);
146 hist = fdNdEtaAnalysis->GetdNdEtaHistogram(2);
147 Float_t value2 = hist->Integral(binLeft, binRight);
149 printf("Ratio is %f, values are %f %f\n", value1 / value2, value1, value2);