Updates in SDD macros to plot QA results
[u/mrichter/AliRoot.git] / ITS / macrosSDD / TrendQAtrainSDD.C
CommitLineData
341adf04 1#if !defined(__CINT__) || defined(__MAKECINT__)
2#include <TH1F.h>
3#include <TH2F.h>
4#include <TF1.h>
5#include <TPad.h>
6#include <TGraphErrors.h>
7#include <TROOT.h>
8#include <TFile.h>
9#include <TTree.h>
10#include <TGrid.h>
11#include <TGridResult.h>
12#include <TMath.h>
13#include <TSystem.h>
14#include <TNtuple.h>
15#include <TCanvas.h>
16#include <TStyle.h>
17#include <TLatex.h>
18#include <TLegend.h>
19#include <TLegendEntry.h>
20#include "AliITSgeomTGeo.h"
21#endif
22
23Double_t LangausFun(Double_t *x, Double_t *par);
24void MakePlots(TString ntupleFileName);
25
26
27void TrendQAtrainSDD(TString period,
28 TString recoPass,
193da5be 29 TString qaTrain="QA",
341adf04 30 Int_t firstRun=0,
31 Int_t lastRun=999999999,
32 TString fileName="QAresults.root"){
33
34 gStyle->SetOptStat(0);
35
36
37 TGrid::Connect("alien:");
38 Int_t year=0;
39 if(period.Contains("LHC09")) year=2009;
40 else if(period.Contains("LHC10")) year=2010;
4d8feb56 41 else if(period.Contains("LHC11")) year=2011;
341adf04 42
193da5be 43 TString outFilNam=Form("TrendingSDD_%s_%s_%s.root",period.Data(),recoPass.Data(),qaTrain.Data());
341adf04 44
45
193da5be 46 const Int_t nVariables=31;
47 TNtuple* ntsdd=new TNtuple("ntsdd","SDD trending","nrun:fracTrackWithClu1:errfracTrackWithClu1:fracTrackWithClu2:errfracTrackWithClu2:fracTrackWithClu3:errfracTrackWithClu3:fracTrackWithClu4:errfracTrackWithClu4:fracTrackWithClu5:errfracTrackWithClu5:fracTrackWithClu6:errfracTrackWithClu6:meanTrPts3:errmeanTrPts3:meanTrPts4:errmeanTrPts4:minDrTime:errminDrTime:meanDrTime:errmeanDrTime:fracExtra:errfracExtra:meandEdxTB0:errmeandEdxTB0:meandEdxTB5:errmeandEdxTB5:nMod95:nMod80:nMod60:nModEmpty");
341adf04 48 Float_t xnt[nVariables];
49
50 TBits* readRun=new TBits(999999);
51 readRun->ResetAllBits();
52 if(!gSystem->Exec(Form("ls -l %s > /dev/null 2>&1",outFilNam.Data()))){
53 TFile* oldfil=new TFile(outFilNam.Data());
54 TNtuple* ntmp=(TNtuple*)oldfil->Get("ntsdd");
55 Bool_t isOK=kFALSE;
56 if(ntmp){
57 if(ntmp->GetNvar()==ntsdd->GetNvar()){
58 isOK=kTRUE;
59 TObjArray* arr1=(TObjArray*)ntsdd->GetListOfBranches();
60 TObjArray* arr2=(TObjArray*)ntmp->GetListOfBranches();
61 for(Int_t iV=0; iV<ntmp->GetNvar(); iV++){
62 TString vnam1=arr1->At(iV)->GetName();
63 TString vnam2=arr2->At(iV)->GetName();
64 if(vnam1!=vnam2) isOK=kFALSE;
65 ntmp->SetBranchAddress(vnam2.Data(),&xnt[iV]);
66 }
67 if(isOK){
68 for(Int_t nE=0; nE<ntmp->GetEntries(); nE++){
69 ntmp->GetEvent(nE);
70 Int_t theRun=(Int_t)(xnt[0]+0.0001);
71 readRun->SetBitNumber(theRun);
72 ntsdd->Fill(xnt);
73 }
74 }
75 }
76 }
77 if(!isOK){
78 printf("Ntuple in local file not OK -> will be recreated\n");
79 }
80 oldfil->Close();
81 delete oldfil;
82 }
83
84 if(!gGrid||!gGrid->IsConnected()) {
85 printf("gGrid not found! exit macro\n");
86 return;
87 }
88
193da5be 89 TString path=Form("/alice/data/%d/%s/",year,period.Data());
341adf04 90 TGridResult *gr = gGrid->Query(path,fileName);
91 Int_t nFiles = gr->GetEntries();
92 printf("================>%d files found\n", nFiles);
93 if (nFiles < 1) return;
94
95
96 if (nFiles > 1){
97 for (Int_t iFil = 0; iFil <nFiles ; iFil++) {
98 TString fileNameLong=Form("%s",gr->GetKey(iFil,"turl"));
99 if(!fileNameLong.Contains(recoPass.Data())) continue;
193da5be 100 if(!fileNameLong.Contains(qaTrain.Data())) continue;
341adf04 101 TString runNumber=fileNameLong;
102 runNumber.ReplaceAll(Form("alien:///alice/data/%d/%s/",year,period.Data()),"");
103 runNumber.Remove(9,runNumber.Sizeof());
104
105 Int_t iRun=atoi(runNumber.Data());
106 if(readRun->TestBitNumber(iRun)){
107 printf("Run %d aleady in local ntuple -> skipping it\n",iRun);
108 continue;
109 }
110 if(iRun<firstRun) continue;
111 if(iRun>lastRun) continue;
112
193da5be 113 TString isMerged=fileNameLong;
114 isMerged.Remove(isMerged.Sizeof()-16);
115 isMerged.Remove(0,isMerged.Sizeof()-5);
116 if(!isMerged.Contains("QA")) continue;
117 printf("Open File %s Run %d\n",fileNameLong.Data(),iRun);
118
119
120
341adf04 121 TFile* f=TFile::Open(fileNameLong.Data());
122
123 TDirectoryFile* df=(TDirectoryFile*)f->Get("SDD_Performance");
124 if(!df){
125 printf("Run %d SDD_Performance MISSING -> Exit\n",iRun);
126 continue;
127 }
128 TList* l=(TList*)df->Get("coutputRP");
129 if(!df){
130 printf("Run %d coutputRP TList MISSING -> Exit\n",iRun);
131 continue;
e1a5dc0a 132 }
133 TH1F* hcllay=(TH1F*)l->FindObject("hCluInLay");
134 Float_t fracT[6]={0.,0.,0.,0.,0.,0.};
135 Float_t efracT[6]={0.,0.,0.,0.,0.,0.};
136 if(hcllay->GetBinContent(1)>0){
137 for(Int_t iLay=0; iLay<6; iLay++){
138 fracT[iLay]=hcllay->GetBinContent(iLay+2)/hcllay->GetBinContent(1);
139 efracT[iLay]=TMath::Sqrt(fracT[iLay]*(1-fracT[iLay])/hcllay->GetBinContent(1));
140 }
341adf04 141 }
142 TH1F* hmodT=(TH1F*)l->FindObject("hTPMod");
143 TH1F* hgamod=(TH1F*)l->FindObject("hGAMod");
144 Int_t bestMod=0;
145 for(Int_t iMod=0; iMod<260;iMod++){
146 Int_t gda=(Int_t)hgamod->GetBinContent(iMod+1);
147 if(gda>bestMod) bestMod=gda;
148 }
149 Int_t nChunks=1;
150 if(bestMod>512){
151 nChunks=(Int_t)(bestMod/512.+0.5);
152 }
153 hgamod->Scale(1./nChunks);
154
155 TH1F* hev=(TH1F*)l->FindObject("hNEvents");
156 Int_t nTotEvents=hev->GetBinContent(2);
157 Int_t nTrigEvents=hev->GetBinContent(3);
158 Int_t nEvents=nTotEvents;
159 printf("Run %d Number of Events = %d Triggered=%d\n",iRun,nTotEvents,nTrigEvents);
160 if(nTrigEvents>0){
161 nEvents=nTrigEvents;
162 }
193da5be 163 if(nTotEvents==0) continue;
341adf04 164 Int_t nModGood3=0;
165 Int_t nModGood4=0;
166 Int_t nModBadAn=0;
167 Float_t sumtp3=0;
168 Float_t sumtp4=0;
169 Float_t sumEtp3=0;
170 Float_t sumEtp4=0;
171 for(Int_t iMod=0; iMod<260; iMod++){
172 Float_t tps=hmodT->GetBinContent(iMod+1);
173 Float_t ga=hgamod->GetBinContent(iMod+1);
174 if(ga<500) nModBadAn++;
175 Float_t tpsN=0.;
176 Float_t etpsN=0.;
177 if(ga>0){
178 tpsN=tps/ga/(Float_t)nEvents;
179 etpsN=TMath::Sqrt(tps)/ga/(Float_t)nEvents;
180 if(iMod<84){
181 sumtp3+=tpsN;
182 sumEtp3+=(etpsN*etpsN);
183 nModGood3++;
184 }
185 else{
186 sumtp4+=tpsN;
187 sumEtp4+=(etpsN*etpsN);
188 nModGood4++;
189 }
190 }
191 }
192
193da5be 193 TH1F* hapmod=(TH1F*)l->FindObject("hAllPmod");
194 TH1F* hgpmod=(TH1F*)l->FindObject("hGoodPmod");
195 // TH1F* hmpmod=(TH1F*)l->FindObject("hMissPmod");
196 TH1F* hbrmod=(TH1F*)l->FindObject("hBadRegmod");
197 TH1F* hskmod=(TH1F*)l->FindObject("hSkippedmod");
198 TH1F* hoamod=(TH1F*)l->FindObject("hOutAccmod");
199 TH1F* hnrmod=(TH1F*)l->FindObject("hNoRefitmod");
200 Int_t nBelow95=0;
201 Int_t nBelow80=0;
202 Int_t nBelow60=0;
203 Int_t nZeroP=0;
204 for(Int_t imod=0; imod<260;imod++){
205 Float_t numer=hgpmod->GetBinContent(imod+1)+hbrmod->GetBinContent(imod+1)+hoamod->GetBinContent(imod+1)+hnrmod->GetBinContent(imod+1)+hskmod->GetBinContent(imod+1);
206 Float_t denom=hapmod->GetBinContent(imod+1);
207 if(denom>0){
208 Float_t eff=numer/denom;
209 if(eff<0.95) nBelow95++;
210 if(eff<0.80) nBelow80++;
211 if(eff<0.60) nBelow60++;
212 }
213 if(hmodT->GetBinContent(imod+1)<1.){
214 nZeroP++;
215 }
216 }
217
341adf04 218 TH1F* htimT=(TH1F*)l->FindObject("hDrTimTPAll");
219 TH1F* htimTe=(TH1F*)l->FindObject("hDrTimTPExtra");
220
221 Double_t fracExtra=0.;
222 Double_t errFracExtra=0.;
223 if(htimT->GetEntries()>0){
224 fracExtra=htimTe->GetEntries()/htimT->GetEntries();
225 errFracExtra=TMath::Sqrt(htimTe->GetEntries())/htimT->GetEntries();
226 }
227 Double_t averPoints=0.;
228 Double_t cntBins=0.;
229 for(Int_t iBin=1; iBin<=htimT->GetNbinsX(); iBin++){
230 Float_t tim=htimT->GetBinCenter(iBin);
231 if(tim>2000. && tim<4000.){
232 averPoints+=htimT->GetBinContent(iBin);
233 cntBins+=1;
234 }
235 }
236 Double_t minTime=-999.;
237 Double_t errMinTime=0.;
238 if(cntBins>0){
239 averPoints/=cntBins;
240 for(Int_t iBin=1; iBin<=htimT->GetNbinsX(); iBin++){
241 if(htimT->GetBinContent(iBin)>0.5*averPoints){
242 minTime=htimT->GetBinCenter(iBin);
243 errMinTime=0.5*htimT->GetBinWidth(iBin);
244 break;
245 }
246 }
247 }
248 TH1F* hSigTim0=(TH1F*)l->FindObject("hSigTimeInt0");
249 TH1F* hSigTim5=(TH1F*)l->FindObject("hSigTimeInt5");
250
e1a5dc0a 251 Int_t index=0;
252 xnt[index++]=(Float_t)iRun;
253 xnt[index++]=fracT[0];
254 xnt[index++]=efracT[0];
255 xnt[index++]=fracT[1];
256 xnt[index++]=efracT[1];
257 xnt[index++]=fracT[2];
258 xnt[index++]=efracT[2];
259 xnt[index++]=fracT[3];
260 xnt[index++]=efracT[3];
261 xnt[index++]=fracT[4];
262 xnt[index++]=efracT[4];
263 xnt[index++]=fracT[5];
264 xnt[index++]=efracT[5];
265 xnt[index++]=sumtp3/nModGood3;
266 xnt[index++]=TMath::Sqrt(sumEtp3)/nModGood3;
267 xnt[index++]=sumtp4/nModGood4;
268 xnt[index++]=TMath::Sqrt(sumEtp4)/nModGood4;
269 xnt[index++]=minTime;
270 xnt[index++]=errMinTime;
271 xnt[index++]=htimT->GetMean();
272 xnt[index++]=htimT->GetMeanError();
273 xnt[index++]=fracExtra;
274 xnt[index++]=errFracExtra;
275 xnt[index++]=hSigTim0->GetMean();
276 xnt[index++]=hSigTim0->GetMeanError();
277 xnt[index++]=hSigTim5->GetMean();
278 xnt[index++]=hSigTim5->GetMeanError();
193da5be 279 xnt[index++]=(Float_t)nBelow95;
280 xnt[index++]=(Float_t)nBelow80;
281 xnt[index++]=(Float_t)nBelow60;
282 xnt[index++]=(Float_t)nZeroP;
341adf04 283 ntsdd->Fill(xnt);
284 }
285 }
286
287 TFile* outfil=new TFile(outFilNam.Data(),"recreate");
288 outfil->cd();
289 ntsdd->Write();
290 outfil->Close();
291
292 MakePlots(outFilNam);
293
294}
295
296void MakePlots(TString ntupleFileName){
297 TFile* fil=new TFile(ntupleFileName.Data(),"read");
298 if(!fil){
299 printf("File with ntuple does not exist\n");
300 return;
301 }
302 TNtuple* ntsdd=(TNtuple*)fil->Get("ntsdd");
303
304 Float_t nrun;
305 Float_t meanTrPts3,errmeanTrPts3,meanTrPts4,errmeanTrPts4;
306 Float_t minDrTime,errminDrTime,meanDrTime,errmeanDrTime;
e1a5dc0a 307 Float_t fracTrackWithClu1,fracTrackWithClu2,errfracTrackWithClu1,errfracTrackWithClu2;
308 Float_t fracTrackWithClu3,fracTrackWithClu4,errfracTrackWithClu3,errfracTrackWithClu4;
309 Float_t fracTrackWithClu5,fracTrackWithClu6,errfracTrackWithClu5,errfracTrackWithClu6;
341adf04 310 Float_t fracExtra,errfracExtra;
311 Float_t meandEdxTB0,errmeandEdxTB0,meandEdxTB5,errmeandEdxTB5;
193da5be 312 Float_t nMod95,nMod80,nMod60,nModEmpty;
313
341adf04 314 ntsdd->SetBranchAddress("nrun",&nrun);
e1a5dc0a 315 ntsdd->SetBranchAddress("fracTrackWithClu1",&fracTrackWithClu1);
316 ntsdd->SetBranchAddress("errfracTrackWithClu1",&errfracTrackWithClu1);
317 ntsdd->SetBranchAddress("fracTrackWithClu2",&fracTrackWithClu2);
318 ntsdd->SetBranchAddress("errfracTrackWithClu2",&errfracTrackWithClu2);
319 ntsdd->SetBranchAddress("fracTrackWithClu3",&fracTrackWithClu3);
320 ntsdd->SetBranchAddress("errfracTrackWithClu3",&errfracTrackWithClu3);
321 ntsdd->SetBranchAddress("fracTrackWithClu4",&fracTrackWithClu4);
322 ntsdd->SetBranchAddress("errfracTrackWithClu4",&errfracTrackWithClu4);
323 ntsdd->SetBranchAddress("fracTrackWithClu5",&fracTrackWithClu5);
324 ntsdd->SetBranchAddress("errfracTrackWithClu5",&errfracTrackWithClu5);
325 ntsdd->SetBranchAddress("fracTrackWithClu6",&fracTrackWithClu6);
326 ntsdd->SetBranchAddress("errfracTrackWithClu6",&errfracTrackWithClu6);
193da5be 327 ntsdd->SetBranchAddress("nMod95",&nMod95);
328 ntsdd->SetBranchAddress("nMod80",&nMod80);
329 ntsdd->SetBranchAddress("nMod60",&nMod60);
330 ntsdd->SetBranchAddress("nModEmpty",&nModEmpty);
e1a5dc0a 331
341adf04 332 ntsdd->SetBranchAddress("meanTrPts3",&meanTrPts3);
333 ntsdd->SetBranchAddress("errmeanTrPts3",&errmeanTrPts3);
334 ntsdd->SetBranchAddress("meanTrPts4",&meanTrPts4);
335 ntsdd->SetBranchAddress("errmeanTrPts4",&errmeanTrPts4);
336 ntsdd->SetBranchAddress("minDrTime",&minDrTime);
337 ntsdd->SetBranchAddress("errminDrTime",&errminDrTime);
338 ntsdd->SetBranchAddress("meanDrTime",&meanDrTime);
339 ntsdd->SetBranchAddress("errmeanDrTime",&errmeanDrTime);
340 ntsdd->SetBranchAddress("fracExtra",&fracExtra);
341 ntsdd->SetBranchAddress("errfracExtra",&errfracExtra);
342 ntsdd->SetBranchAddress("meandEdxTB0",&meandEdxTB0);
343 ntsdd->SetBranchAddress("errmeandEdxTB0",&errmeandEdxTB0);
344 ntsdd->SetBranchAddress("meandEdxTB5",&meandEdxTB5);
345 ntsdd->SetBranchAddress("errmeandEdxTB5",&errmeandEdxTB5);
346
347 TH1F* histotrp3=new TH1F("histotrp3","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
348 TH1F* histotrp4=new TH1F("histotrp4","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
349 TH1F* histominTime=new TH1F("histominTime","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
350 TH1F* histomeanTime=new TH1F("histomeanTime","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
351 TH1F* histofracExtra=new TH1F("histofracExtra","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
352 TH1F* histodEdxTB0=new TH1F("histodEdxTB0","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
353 TH1F* histodEdxTB5=new TH1F("histodEdxTB5","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
e1a5dc0a 354 TH1F* histoTrackClu1=new TH1F("histoTrackClu1","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
355 TH1F* histoTrackClu2=new TH1F("histoTrackClu2","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
356 TH1F* histoTrackClu3=new TH1F("histoTrackClu3","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
357 TH1F* histoTrackClu4=new TH1F("histoTrackClu4","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
358 TH1F* histoTrackClu5=new TH1F("histoTrackClu5","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
359 TH1F* histoTrackClu6=new TH1F("histoTrackClu6","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
360
193da5be 361 TH1F* histoNmodEffBelow95=new TH1F("histoNmodEffBelow95","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
362 TH1F* histoNmodEffBelow80=new TH1F("histoNmodEffBelow80","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
363 TH1F* histoNmodEffBelow60=new TH1F("histoNmodEffBelow60","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
364 TH1F* histoNmodEmpty=new TH1F("histoNmodEmpty","",(Int_t)ntsdd->GetEntries(),0.,ntsdd->GetEntries());
341adf04 365
366 for(Int_t i=0; i<ntsdd->GetEntries();i++){
367 ntsdd->GetEvent(i);
e1a5dc0a 368 histoTrackClu1->SetBinContent(i+1,fracTrackWithClu1);
369 histoTrackClu1->SetBinError(i+1,errfracTrackWithClu1);
370 histoTrackClu1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
371 histoTrackClu2->SetBinContent(i+1,fracTrackWithClu2);
372 histoTrackClu2->SetBinError(i+1,errfracTrackWithClu2);
373 histoTrackClu2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
374 histoTrackClu3->SetBinContent(i+1,fracTrackWithClu3);
375 histoTrackClu3->SetBinError(i+1,errfracTrackWithClu3);
376 histoTrackClu3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
377 histoTrackClu4->SetBinContent(i+1,fracTrackWithClu4);
378 histoTrackClu4->SetBinError(i+1,errfracTrackWithClu4);
379 histoTrackClu4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
380 histoTrackClu5->SetBinContent(i+1,fracTrackWithClu5);
381 histoTrackClu5->SetBinError(i+1,errfracTrackWithClu5);
382 histoTrackClu5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
383 histoTrackClu6->SetBinContent(i+1,fracTrackWithClu6);
384 histoTrackClu6->SetBinError(i+1,errfracTrackWithClu6);
385 histoTrackClu6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
341adf04 386 histominTime->SetBinContent(i+1,minDrTime);
387 histominTime->SetBinError(i+1,errminDrTime);
388 histominTime->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
389 histomeanTime->SetBinContent(i+1,meanDrTime);
390 histomeanTime->SetBinError(i+1,errmeanDrTime);
391 histomeanTime->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
392 histotrp3->SetBinContent(i+1,meanTrPts3);
393 histotrp3->SetBinError(i+1,errmeanTrPts3);
394 histotrp3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
395 histotrp4->SetBinContent(i+1,meanTrPts4);
396 histotrp4->SetBinError(i+1,errmeanTrPts3);
397 histotrp4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
398 histofracExtra->SetBinContent(i+1,fracExtra);
399 histofracExtra->SetBinError(i+1,errfracExtra);
400 histofracExtra->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
401 histodEdxTB0->SetBinContent(i+1,meandEdxTB0);
402 histodEdxTB0->SetBinError(i+1,errmeandEdxTB0);
403 histodEdxTB0->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
404 histodEdxTB5->SetBinContent(i+1,meandEdxTB5);
405 histodEdxTB5->SetBinError(i+1,errmeandEdxTB5);
406 histodEdxTB5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
193da5be 407
408 histoNmodEffBelow95->SetBinContent(i+1,nMod95);
409 histoNmodEffBelow95->SetBinError(i+1,0.0001);
410 histoNmodEffBelow95->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
411 histoNmodEffBelow80->SetBinContent(i+1,nMod80);
412 histoNmodEffBelow80->SetBinError(i+1,0.0001);
413 histoNmodEffBelow80->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
414 histoNmodEffBelow60->SetBinContent(i+1,nMod60);
415 histoNmodEffBelow60->SetBinError(i+1,0.0001);
416 histoNmodEffBelow60->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
417 histoNmodEmpty->SetBinContent(i+1,nModEmpty);
418 histoNmodEmpty->SetBinError(i+1,0.0001);
419 histoNmodEmpty->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
341adf04 420 }
421
422 gStyle->SetOptStat(0);
423
424
425 TCanvas* c1=new TCanvas("c1","Occupancy");
426 histotrp3->SetLineColor(1);
427 histotrp3->SetMarkerStyle(20);
428 histotrp3->Draw();
429 histotrp3->SetMinimum(0.);
430 histotrp4->SetLineColor(2);
431 histotrp4->SetMarkerColor(2);
432 histotrp4->SetMarkerStyle(22);
433 histotrp4->Draw("same");
434 histotrp3->GetYaxis()->SetTitle("Track Point Occupancy");
435 TLegend* leg=new TLegend(0.7,0.15,0.88,0.35);
436 TLegendEntry* ent=leg->AddEntry(histotrp3,"Layer3","PL");
437 ent=leg->AddEntry(histotrp4,"Layer4","PL");
438 ent->SetTextColor(histotrp4->GetMarkerColor());
439 leg->SetFillStyle(0);
440 leg->Draw();
441 c1->Update();
442
443 TCanvas* c2=new TCanvas("c2","DriftTime",600,900);
444 c2->Divide(1,2);
445 c2->cd(1);
446 histominTime->Draw();
447 histominTime->SetMinimum(250);
448 histominTime->SetMaximum(1000);
449 histominTime->GetYaxis()->SetTitle("Minimum Drift Time (ns)");
450 c2->cd(2);
451 histomeanTime->Draw();
452 histomeanTime->GetYaxis()->SetTitle("Average Drift Time (ns)");
453 c2->Update();
454
455 TCanvas* c3=new TCanvas("c3","ExtraClusters");
456 histofracExtra->Draw();
457 histofracExtra->GetYaxis()->SetTitle("Fraction of Extra Clusters");
458 c3->Update();
459
460
461 TCanvas* c4=new TCanvas("c4","Charge");
462 histodEdxTB0->SetLineColor(1);
463 histodEdxTB0->SetMarkerStyle(20);
464 histodEdxTB0->Draw();
465 histodEdxTB0->SetMinimum(0.);
466 histodEdxTB5->SetLineColor(4);
467 histodEdxTB5->SetMarkerColor(4);
468 histodEdxTB5->SetMarkerStyle(23);
469 histodEdxTB5->Draw("same");
e1a5dc0a 470 histodEdxTB0->GetYaxis()->SetTitle("<dE/dx> (keV/300 #mum)");
341adf04 471 TLegend* leg2=new TLegend(0.6,0.15,0.88,0.35);
472 ent=leg2->AddEntry(histodEdxTB0,"Small drift time","PL");
473 ent=leg2->AddEntry(histodEdxTB5,"Large drift time","PL");
474 ent->SetTextColor(histodEdxTB5->GetMarkerColor());
475 leg2->SetFillStyle(0);
476 leg2->Draw();
477 c4->Update();
478
e1a5dc0a 479 TCanvas* c5=new TCanvas("c5","TrackWithSDD");
480 histoTrackClu3->Draw();
481 histoTrackClu3->SetLineColor(1);
482 histoTrackClu3->SetMarkerStyle(20);
483 histoTrackClu3->Draw();
484 histoTrackClu3->SetMinimum(0.);
485 histoTrackClu4->SetLineColor(2);
486 histoTrackClu4->SetMarkerColor(2);
487 histoTrackClu4->SetMarkerStyle(22);
488 histoTrackClu4->Draw("same");
489 histoTrackClu1->SetLineColor(kGray+1);
490 histoTrackClu1->SetMarkerColor(kGray+1);
491 histoTrackClu1->SetMarkerStyle(24);
492 histoTrackClu1->Draw("same");
493 histoTrackClu2->SetLineColor(kGray+2);
494 histoTrackClu2->SetMarkerColor(kGray+2);
495 histoTrackClu2->SetMarkerStyle(26);
496 histoTrackClu2->Draw("same");
497 histoTrackClu5->SetLineColor(4);
498 histoTrackClu5->SetMarkerColor(4);
499 histoTrackClu5->SetMarkerStyle(29);
500 histoTrackClu5->Draw("same");
501 histoTrackClu6->SetLineColor(kBlue+1);
502 histoTrackClu6->SetMarkerColor(kBlue+1);
503 histoTrackClu6->SetMarkerStyle(30);
504 histoTrackClu6->Draw("same");
505 histoTrackClu3->GetYaxis()->SetTitle("Fraction of Tracks with Point in SDD Layers");
506 TLegend* leg3=new TLegend(0.7,0.15,0.88,0.35);
507 ent=leg3->AddEntry(histoTrackClu1,"Layer1","PL");
508 ent->SetTextColor(histoTrackClu1->GetMarkerColor());
509 ent=leg3->AddEntry(histoTrackClu2,"Layer2","PL");
510 ent->SetTextColor(histoTrackClu2->GetMarkerColor());
511 ent=leg3->AddEntry(histoTrackClu3,"Layer3","PL");
512 ent->SetTextColor(histoTrackClu3->GetMarkerColor());
513 ent=leg3->AddEntry(histoTrackClu4,"Layer4","PL");
514 ent->SetTextColor(histoTrackClu4->GetMarkerColor());
515 ent=leg3->AddEntry(histoTrackClu5,"Layer5","PL");
516 ent->SetTextColor(histoTrackClu5->GetMarkerColor());
517 ent=leg3->AddEntry(histoTrackClu6,"Layer6","PL");
518 ent->SetTextColor(histoTrackClu6->GetMarkerColor());
519
520 leg3->SetFillStyle(0);
521 leg3->Draw();
522 c5->Update();
341adf04 523
193da5be 524 TCanvas* c6=new TCanvas("c6","Modules with low eff",800,1000);
525 c6->Divide(1,4);
526 c6->cd(1);
527 histoNmodEffBelow95->Draw("E");
528 c6->cd(2);
529 histoNmodEffBelow80->Draw("E");
530 c6->cd(3);
531 histoNmodEffBelow60->Draw("E");
532 c6->cd(4);
533 histoNmodEmpty->Draw("E");
534
341adf04 535}
536
537Double_t LangausFun(Double_t *x, Double_t *par) {
538
539 //Fit parameters:
540 //par[0]=Width (scale) parameter of Landau density
541 //par[1]=Most Probable (MP, location) parameter of Landau density
542 //par[2]=Total area (integral -inf to inf, normalization constant)
543 //par[3]=Width (sigma) of convoluted Gaussian function
544 //
545 //In the Landau distribution (represented by the CERNLIB approximation),
546 //the maximum is located at x=-0.22278298 with the location parameter=0.
547 //This shift is corrected within this function, so that the actual
548 //maximum is identical to the MP parameter.
549
550 // Numeric constants
551 Double_t invsq2pi = 0.3989422804014; // (2 pi)^(-1/2)
552 Double_t mpshift = -0.22278298; // Landau maximum location
553
554 // Control constants
555 Double_t np = 100.0; // number of convolution steps
556 Double_t sc = 5.0; // convolution extends to +-sc Gaussian sigmas
557
558 // Variables
559 Double_t xx;
560 Double_t mpc;
561 Double_t fland;
562 Double_t sum = 0.0;
563 Double_t xlow,xupp;
564 Double_t step;
565 Double_t i;
566
567
568 // MP shift correction
569 mpc = par[1] - mpshift * par[0];
570
571 // Range of convolution integral
572 xlow = x[0] - sc * par[3];
573 xupp = x[0] + sc * par[3];
574
575 step = (xupp-xlow) / np;
576
577 // Convolution integral of Landau and Gaussian by sum
578 for(i=1.0; i<=np/2; i++) {
579 xx = xlow + (i-.5) * step;
580 fland = TMath::Landau(xx,mpc,par[0]) / par[0];
581 sum += fland * TMath::Gaus(x[0],xx,par[3]);
582
583 xx = xupp - (i-.5) * step;
584 fland = TMath::Landau(xx,mpc,par[0]) / par[0];
585 sum += fland * TMath::Gaus(x[0],xx,par[3]);
586 }
587
588 return (par[2] * step * sum * invsq2pi / par[3]);
589}