trigger names changed (Markus)
[u/mrichter/AliRoot.git] / PWG1 / TRD / runQAesds.C
1 Int_t ExtractRunNumber(TString filename);
2 void Plot(Int_t run);
3 void runQAesds(const char *esdlist)
4 {
5   ifstream input(esdlist);
6   string str;
7   getline(input, str);
8   TString tstr(str.c_str());
9   if(tstr.Contains("alien://")){
10     printf("Data from alien\n");
11     gSystem->Load("libRAliEn");
12     TGrid::Connect("alien://");
13   }
14
15   Int_t runnr = ExtractRunNumber(tstr);
16   printf("runnr %d\n", runnr);
17
18   gSystem->Load("libANALYSIS");
19   gSystem->Load("libANALYSISalice");
20   gSystem->Load("libTENDER");
21   gSystem->Load("libPWG1");
22   AliLog::SetGlobalLogLevel(AliLog::kError);
23   gROOT->LoadMacro(Form("%s/PWG0/CreateESDChain.C", gSystem->ExpandPathName("$ALICE_ROOT")));
24   TChain *chain = CreateESDChain(esdlist, -1);
25   chain->SetBranchStatus("*FMD*",0);
26   chain->SetBranchStatus("*Calo*",0);
27   chain->SetBranchStatus("Tracks", 1);
28   chain->SetBranchStatus("ESDfriend*",1);
29   chain->Lookup();
30   chain->GetListOfFiles()->Print();
31   printf("\n ----> CHAIN HAS %d ENTRIES <----\n\n", (Int_t)chain->GetEntries());
32     
33   AliAnalysisManager *TRDqa = new AliAnalysisManager("TRD Cosmics QA");
34   TRDqa->SetInputEventHandler(new AliESDInputHandler);
35   AliTRDinfoGen *info = new AliTRDinfoGen();
36   AliTRDcheckDET *detector = new AliTRDcheckDET();
37   detector->SetDebugLevel(10);
38   detector->SetMCdata(kFALSE);
39   detector->UseClustersOutsideChamber();
40   TRDqa->AddTask(info);
41   TRDqa->AddTask(detector);
42   
43   AliAnalysisDataContainer *co[3];
44   TRDqa->ConnectInput( info, 0, TRDqa->GetCommonInputContainer());
45   co[0] = TRDqa->CreateContainer("trackInfo", TObjArray::Class(), AliAnalysisManager::kExchangeContainer);
46   co[1] = TRDqa->CreateContainer("eventInfo", AliTRDeventInfo::Class(), AliAnalysisManager::kExchangeContainer);
47   co[2] = TRDqa->CreateContainer("v0Info", TObjArray::Class(), AliAnalysisManager::kExchangeContainer);
48   TRDqa->ConnectOutput(info, 0, co[0]);
49   TRDqa->ConnectOutput(info, 1, co[1]);
50   TRDqa->ConnectOutput(info, 2, co[2]);
51   TRDqa->ConnectInput(detector, 0, co[0]);
52   TRDqa->ConnectInput(detector, 1, co[1]);
53   TRDqa->ConnectOutput(detector, 0, TRDqa->CreateContainer(detector->GetName(), TObjArray::Class(), AliAnalysisManager::kOutputContainer, Form("TRD.Perf%d.root", runnr)));
54   
55   if(TRDqa->InitAnalysis()){
56           TRDqa->Print();
57           TRDqa->StartAnalysis("local", chain);
58   }
59   gSystem->Exec(Form("mv TRD.DebugPerformance.root TRD.DbgPerf%d.root", runnr));
60   Plot(runnr);
61 }
62
63 void Plot(Int_t run){
64   TFile *f=new TFile(Form("TRD.Perf%d.root",run));
65   TH1F *hNtrSec = (TH1F*)checkDET->FindObject("hNtrksSector");
66   TH1F *hNtrks = (TH1F*)checkDET->FindObject("hNtrks"); // tracks/event
67   TH1F *hNtlsSTA = (TH1F*)checkDET->FindObject("hNtlsSTA");
68   TH1F *hNtlsBAR = (TH1F*)checkDET->FindObject("hNtlsBAR");
69   TH1F *hNclTrk = (TH1F*)checkDET->FindObject("hNclTls");
70   TProfile *hPHp = (TProfile*) checkDET->FindObject("<PH>")->FindObject("hPHt");
71   f->Close();
72
73   TH2F *hPH = new TH2F("hPH", "PH", 30, 0, 30, 100, 0, 1000);
74   TString filename = Form("TRD.DbgPerf%d.root",run);
75   TChain *debchain = new TChain("PHt");
76   if(filename.EndsWith(".root")) {  
77     debchain->AddFile(filename.Data());
78     debchain->Draw("Charge:Timebin>>hPH", "","colz");
79   } else{
80     cout << "No such file... " << filename <<endl;
81   }
82
83   const Float_t Ls1=0.04;
84   TCanvas *tc = new TCanvas ("Monitor", "Monitor",50,50,600,900);
85   tc->Divide(2,3);
86
87   tc->cd(1); 
88   pad = tc->cd(1); pad->SetLogy();
89   hNtrks->GetXaxis()->SetTitle("Tracks per event");
90   hNtrks->GetXaxis()->SetTitleOffset(.9);
91   hNtrks->GetXaxis()->SetTitleSize(0.05);
92   hNtrks->GetYaxis()->SetTitle("");
93   //hNtrks->SetTitle(Form("%d",run));
94   hNtrks->SetLineWidth(2);
95   hNtrks->SetLabelSize(Ls1,"xy");
96   hNtrks->GetXaxis()->SetRangeUser(0,20);    
97   hNtrks->Draw();
98   int Ntracks=hNtrks->Integral();
99   TLatex *text=new TLatex(8,10,Form("%d tracks",Ntracks));
100   text->SetTextSize(0.045); text->SetTextFont(42);
101   text->Draw("same");
102
103   tc->cd(2);
104   hNtrSec->GetXaxis()->SetTitle("Sector");
105   hNtrSec->GetXaxis()->SetTitleOffset(.9);
106   hNtrSec->GetXaxis()->SetTitleSize(0.05);
107   hNtrSec->GetYaxis()->SetTitle("Tracks");
108   //hNtrSec->GetYaxis()->SetTitleSize(0.05);
109   hNtrSec->GetYaxis()->SetTitleOffset(1.2);
110   hNtrSec->SetLineWidth(2);
111   hNtrSec->SetLabelSize(Ls1,"xy");
112   hNtrSec->Draw();
113
114   tc->cd(3);
115   hNtlsBAR->GetXaxis()->SetTitle("Tracklets per track");
116   hNtlsBAR->GetXaxis()->SetTitleOffset(.9);
117   hNtlsBAR->GetXaxis()->SetTitleSize(0.05);
118   hNtlsBAR->GetYaxis()->SetTitle("");
119   //hNtlsBAR->SetStats(1111);
120   hNtlsBAR->SetLabelSize(Ls1,"xy");
121   hNtlsBAR->SetLineWidth(2);
122   //hNtlsBAR->GetYaxis()->SetRangeUser(0,500);    
123   hNtlsSTA->GetXaxis()->SetTitle("Tracklets per track");
124   hNtlsSTA->GetYaxis()->SetTitle("");
125   hNtlsBAR->SetLineColor(2); hNtlsBAR->Draw("");
126   hNtlsSTA->SetLineColor(4); hNtlsSTA->Draw("same");
127   leg = new TLegend(0.14,0.35,0.65,0.5);
128   leg->AddEntry(hNtlsBAR,"Barrel","l");
129   leg->AddEntry(hNtlsSTA,"Standalone","l");
130   leg->SetBorderSize(0);
131   leg->SetFillStyle(0);
132   leg->SetMargin(0.15); //separation symbol-text
133   leg->SetEntrySeparation(0.1); //this seems not to do much...
134   leg->Draw();
135
136   tc->cd(4);
137   hNclTrk->GetXaxis()->SetTitle("Clusters per tracklet");
138   hNclTrk->GetXaxis()->SetTitleOffset(.9);
139   hNclTrk->GetXaxis()->SetTitleSize(0.05);
140   hNclTrk->GetYaxis()->SetTitle("");
141   hNclTrk->SetLabelSize(Ls1,"xy");
142   hNclTrk->SetLineWidth(2);
143   hNclTrk->Draw();
144
145   tc->cd(5);
146   hPHp->GetXaxis()->SetTitle("Time bin (100 ns)");
147   hPHp->GetXaxis()->SetTitleOffset(.9);
148   hPHp->GetXaxis()->SetTitleSize(0.05);
149   hPHp->GetYaxis()->SetTitle("<PH> (ADC ch.)");
150   hPHp->GetYaxis()->SetTitleOffset(1.3);
151   hPHp->SetLabelSize(Ls1,"xy");
152   hPHp->SetLineWidth(2);
153   hPHp->Draw();
154
155   tc->cd(6);
156   hPH->GetXaxis()->SetTitle("Time bin (100 ns)");
157   hPH->GetXaxis()->SetTitleOffset(.9);
158   hPH->GetXaxis()->SetTitleSize(0.05);
159   hPH->GetYaxis()->SetTitle("Pulse Height (ADC ch.)");
160   hPH->GetYaxis()->SetTitleOffset(1.3);
161   hPH->SetLabelSize(Ls1,"xy");
162   hPH->SetLineWidth(2);
163   hPH->Draw("colz");
164
165   //tc->SaveAs(Form("Mon_run%dp%d.eps",run,pass));
166   tc->SaveAs(Form("Mon_run%d.gif",run));
167 }
168
169 Int_t ExtractRunNumber(TString filename){
170   TObjArray *toks = filename.Tokenize("/");
171   TObjString *file = dynamic_cast<TObjString *>(toks->At(toks->GetEntries() - 2));
172   TString &filestr = file->String();
173   TString runstr = filestr(2,9);
174   return runstr.Atoi();
175 }