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