]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGPP/ITS/AliITSQAtrendAlign.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
better log validation
[u/mrichter/AliRoot.git] / PWGPP / ITS / AliITSQAtrendAlign.C
CommitLineData
839ef7fc 1#if !defined(__CINT__) || defined(__MAKECINT__)
2#include <TH1F.h>
3#include <TH2F.h>
4#include <TF1.h>
5#include <TProfile.h>
6#include <TPad.h>
7#include <TGraphErrors.h>
8#include <TROOT.h>
9#include <TFile.h>
10#include <TTree.h>
11#include <TGrid.h>
12#include <TGridResult.h>
13#include <TMath.h>
14#include <TSystem.h>
15#include <TNtuple.h>
16#include <TCanvas.h>
17#include <TStyle.h>
18#include <TLatex.h>
19#include <TLegend.h>
20#include <TLegendEntry.h>
21#include <Riostream.h>
22#include "AliITSgeomTGeo.h"
23#include <TLine.h>
24#include "AliCDBManager.h"
25#include "AliCDBEntry.h"
26#include "AliITSCalibration.h"
27#include "AliITSOnlineCalibrationSPDhandler.h"
28#include "AliITSBadChannelsSSD.h"
29#include "AliITSBadChannelsSSDv2.h"
30#include "AliITSgeomTGeo.h"
31#include <TObjArray.h>
32#endif
33
34TString pdfFileNames="";
35void MakePlots(Int_t run1=-1,Int_t run2=99999999,TString ntupleFileName="TrendingITSAlign_norm.root");
36void ReadOldSSDBadChannels(TObjArray *array, AliITSBadChannelsSSDv2 *badChannelsSSD);
37void AliITSQAtrendAlign(TString runListFile="QA_Trend_list.txt",TString ntupleFileName="TrendingITSAlign_norm.root");
38
39////////////////////////////////////////////////////////////////
40// Please, read this comment before using this macro
41//
42// INPUT FILE: a text file (by default QA_Trend_list.txt) which contains
43// a list of the complete path+file name of the QAresults root files
44// (without the alien:// prefix).
45// One file per line. The order is irrelevant.
46//
47// USAGE:
48//
49// Function AliITSQAtrendAlign():
50// it looks for a local root file named TrendingITSAlign_norm.root,
51// where the ntuples used to build the trending plots are
52// stored. This file is used to generate incrementally the ntuple contents:
53// when you add a new entry in the QA_Trend_list.txt file and you have a
54// local copy of the TrendingITSAlign_norm.root file, only the additional run will
55// be processed. The whole list is processed only the first time you use the
56// macro.
57// For each run, the macro performs a gaussian fit of the xlocal and zlocal
58// residuals from QA train of each module of each ITS layer;
59// then evaluates the fractions of modules
60// whose xlocal and zlocal residual mean and sigma are smaller than chosen
61// threshold values; the fractions are normalized to the number of good modules
62// directly read from OCDB fpor each run.
63// Please, bear in mind that this macro is RAM-intensive: all the
64// ntuples are kept in memory. It is better to add few runs to the list at
65// each time, according to the RAM of your computer.
66// The function AliITSQAtrendAlign does not produce any plot.
67//
68// Function MakePlots(run1,run2):
69// it produces the plots. For each canvas a PDF file is created.
70// A PDF file with all the canvases merged is also produced, named by default TrendingITSAlign_norm.root
71// The first two arguments define a range for the runs to be displayed
72// These two arguments are optional: by default, all the runs
73// found in the ntuples are displayed
74////////////////////////////////////////////////////////////////
75
76/* $Id$ */
77
78//_____________________________________________________________________//
79void ReadOldSSDBadChannels(TObjArray *array,
80 AliITSBadChannelsSSDv2 *badChannelsSSD) {
81 Int_t fNMod = array->GetEntries();
82 cout<<"Converting old calibration object for bad channels..."<<endl;
83 for (Int_t iModule = 0; iModule < fNMod; iModule++) {
84 //for (Int_t iModule = 0; iModule < 1; iModule++) {
85 AliITSBadChannelsSSD *bad = (AliITSBadChannelsSSD*) (array->At(iModule));
86 TArrayI arrayPSide = bad->GetBadPChannelsList();
87 for(Int_t iPCounter = 0; iPCounter < arrayPSide.GetSize(); iPCounter++)
88 badChannelsSSD->AddBadChannelP(iModule,
89 iPCounter,
90 (Char_t)arrayPSide.At(iPCounter));
91
92 TArrayI arrayNSide = bad->GetBadNChannelsList();
93 for(Int_t iNCounter = 0; iNCounter < arrayNSide.GetSize(); iNCounter++)
94 badChannelsSSD->AddBadChannelN(iModule,
95 iNCounter,
96 (Char_t)arrayNSide.At(iNCounter));
97
98 }//loop over modules
99}
100
101//_____________________________________________________________________//
102
103
104void AliITSQAtrendAlign(TString runListFile,TString ntupleFileName){
105
106 TGrid::Connect("alien:");
107
108 const Int_t nVariables=169;
109 TNtuple* ntits=new TNtuple("ntitsalign","ITS ALIGN trending","nrun:SDDfracModmux20_3:errSDDfracModmux20_3:SDDfracModmux50_3:errSDDfracModmux50_3:SDDfracModmux100_3:errSDDfracModmux100_3:SDDfracModsigx100_3:errSDDfracModsigx100_3:SDDfracModsigx200_3:errSDDfracModsigx200_3:SDDfracModsigx300_3:errSDDfracModsigx300_3:SDDfracModmuxlr20_3:errSDDfracModmuxlr20_3:SDDfracModmuxlr50_3:errSDDfracModmuxlr50_3:SDDfracModmuxlr100_3:errSDDfracModmuxlr100_3:SDDfracModexcx100_3:errSDDfracModexcx100_3:SDDfracModexcx200_3:errSDDfracModexcx200_3:SDDfracModexcx300_3:errSDDfracModexcx300_3:SDDfracModmuz50_3:errSDDfracModmuz50_3:SDDfracModmuz100_3:errSDDfracModmuz100_3:SDDfracModmuz300_3:errSDDfracModmuz300_3:SDDfracModsigz100_3:errSDDfracModsigz100_3:SDDfracModsigz300_3:errSDDfracModsigz300_3:SDDfracModsigz500_3:errSDDfracModsigz500_3:SDDfracModmux20_4:errSDDfracModmux20_4:SDDfracModmux50_4:errSDDfracModmux50_4:SDDfracModmux100_4:errSDDfracModmux100_4:SDDfracModsigx100_4:errSDDfracModsigx100_4:SDDfracModsigx200_4:errSDDfracModsigx200_4:SDDfracModsigx300_4:errSDDfracModsigx300_4:SDDfracModmuxlr20_4:errSDDfracModmuxlr20_4:SDDfracModmuxlr50_4:errSDDfracModmuxlr50_4:SDDfracModmuxlr100_4:errSDDfracModmuxlr100_4:SDDfracModexcx100_4:errSDDfracModexcx100_4:SDDfracModexcx200_4:errSDDfracModexcx200_4:SDDfracModexcx300_4:errSDDfracModexcx300_4:SDDfracModmuz50_4:errSDDfracModmuz50_4:SDDfracModmuz100_4:errSDDfracModmuz100_4:SDDfracModmuz300_4:errSDDfracModmuz300_4:SDDfracModsigz100_4:errSDDfracModsigz100_4:SDDfracModsigz300_4:errSDDfracModsigz300_4:SDDfracModsigz500_4:errSDDfracModsigz500_4:SPDfracModmux20_1:errSPDfracModmux20_1:SPDfracModmux50_1:errSPDfracModmux50_1:SPDfracModmux100_1:errSPDfracModmux100_1:SPDfracModsigx100_1:errSPDfracModsigx100_1:SPDfracModsigx200_1:errSPDfracModsigx200_1:SPDfracModsigx300_1:errSPDfracModsigx300_1:SPDfracModmuz50_1:errSPDfracModmuz50_1:SPDfracModmuz100_1:errSPDfracModmuz100_1:SPDfracModmuz300_1:errSPDfracModmuz300_1:SPDfracModsigz100_1:errSPDfracModsigz100_1:SPDfracModsigz300_1:errSPDfracModsigz300_1:SPDfracModsigz500_1:errSPDfracModsigz500_1:SPDfracModmux20_2:errSPDfracModmux20_2:SPDfracModmux50_2:errSPDfracModmux50_2:SPDfracModmux100_2:errSPDfracModmux100_2:SPDfracModsigx100_2:errSPDfracModsigx100_2:SPDfracModsigx200_2:errSPDfracModsigx200_2:SPDfracModsigx300_2:errSPDfracModsigx300_2:SPDfracModmuz50_2:errSPDfracModmuz50_2:SPDfracModmuz100_2:errSPDfracModmuz100_2:SPDfracModmuz300_2:errSPDfracModmuz300_2:SPDfracModsigz100_2:errSPDfracModsigz100_2:SPDfracModsigz300_2:errSPDfracModsigz300_2:SPDfracModsigz500_2:errSPDfracModsigz500_2:SSDfracModmux20_5:errSSDfracModmux20_5:SSDfracModmux50_5:errSSDfracModmux50_5:SSDfracModmux100_5:errSSDfracModmux100_5:SSDfracModsigx100_5:errSSDfracModsigx100_5:SSDfracModsigx200_5:errSSDfracModsigx200_5:SSDfracModsigx300_5:errSSDfracModsigx300_5:SSDfracModmuz50_5:errSSDfracModmuz50_5:SSDfracModmuz100_5:errSSDfracModmuz100_5:SSDfracModmuz300_5:errSSDfracModmuz300_5:SSDfracModsigz100_5:errSSDfracModsigz100_5:SSDfracModsigz300_5:errSSDfracModsigz300_5:SSDfracModsigz500_5:errSSDfracModsigz500_5:SSDfracModmux20_6:errSSDfracModmux20_6:SSDfracModmux50_6:errSSDfracModmux50_6:SSDfracModmux100_6:errSSDfracModmux100_6:SSDfracModsigx100_6:errSSDfracModsigx100_6:SSDfracModsigx200_6:errSSDfracModsigx200_6:SSDfracModsigx300_6:errSSDfracModsigx300_6:SSDfracModmuz50_6:errSSDfracModmuz50_6:SSDfracModmuz100_6:errSSDfracModmuz100_6:SSDfracModmuz300_6:errSSDfracModmuz300_6:SSDfracModsigz100_6:errSSDfracModsigz100_6:SSDfracModsigz300_6:errSSDfracModsigz300_6:SSDfracModsigz500_6:errSSDfracModsigz500_6");
110 Float_t xnt[nVariables];
111
112 // inizializzazione variabili OCDB
113
114 AliCDBManager* man = AliCDBManager::Instance();
115 man->SetDefaultStorage("alien://folder=/alice/data/2011/OCDB");
116 AliCDBEntry *entrySDD;
117 AliITSCalibration *cal = NULL;
118 TObjArray *calSDD;
119//
120
121 TBits* readRun=new TBits(999999);
122 readRun->ResetAllBits();
123
124 if(!gSystem->Exec(Form("ls -l %s > /dev/null 2>&1",ntupleFileName.Data()))){ // esiste gia' un file con l'ntupla ...
125 TFile* oldfil=new TFile(ntupleFileName.Data()); // che viene letta per aggiungere altre entries ....
126 TNtuple* ntmp=(TNtuple*)oldfil->Get("ntitsalign");
127
128 Bool_t isOK=kFALSE;
129 if(ntmp){
130 // cout << "old ntu = " << ntmp->GetNvar() << ", new ntu " << ntits->GetNvar() << endl;
131 if(ntmp->GetNvar()==ntits->GetNvar()){
132 isOK=kTRUE;
133 TObjArray* arr1=(TObjArray*)ntits->GetListOfBranches();
134 TObjArray* arr2=(TObjArray*)ntmp->GetListOfBranches();
135 for(Int_t iV=0; iV<ntmp->GetNvar(); iV++){
136 TString vnam1=arr1->At(iV)->GetName();
137 TString vnam2=arr2->At(iV)->GetName();
138 if(vnam1!=vnam2) isOK=kFALSE;
139 ntmp->SetBranchAddress(vnam2.Data(),&xnt[iV]);
140 }
141 if(isOK){
142 for(Int_t nE=0; nE<ntmp->GetEntries(); nE++){
143 ntmp->GetEvent(nE);
144 Int_t theRun=(Int_t)(xnt[0]+0.0001);
145 readRun->SetBitNumber(theRun);
146 ntits->Fill(xnt);
147 }
148 }
149 }
150 }
151 if(!isOK){
152 printf("Ntuple in local file not OK -> will be recreated\n");
153 }
154 oldfil->Close();
155 delete oldfil;
156 }
157
158#define MAX_LINES 200
159#define MAX_LINE_LEN 255
160
161 char strings[MAX_LINES][MAX_LINE_LEN];
162 ifstream in(runListFile.Data());
163 int j = 0;
164 Int_t nrun=0;
165 Int_t runNumb[MAX_LINES];
166 Bool_t goout = kFALSE;
167 while ( in ) {
168 in.getline(strings[j], MAX_LINE_LEN);
169 TString aux(strings[j]);
170 Int_t lentrail=0;
171 if(aux.Contains("LHC11h/") || aux.Contains("LHC11e/")){
172 lentrail = 27;
173 }
174 else if(aux.Contains("LHC11h_2/")){
175 lentrail = 29;
176 }
177 else {
178 if(!aux.IsNull())printf("Unrecognised path name %s \n",aux.Data());
179 goout = kTRUE;
180 }
181 if(goout)break;
182 if(aux.Length()<lentrail)continue;
183 aux=aux.Remove(0,lentrail);
184 aux=aux.Remove(6,aux.Length());
185 runNumb[j]=atoi(aux.Data());
186 printf("%d ) - path %s \n",runNumb[j],strings[j]);
187 j++;
188 nrun++;
189 }
190
191 // loop on runs
192
193 printf("\n ******************* Loop on runs *********** \n");
194
195 Int_t twice=0;
196
197 for(Int_t jru=0;jru<nrun;jru++) {
198 printf("jru=%d - run number= %d \n",jru,runNumb[jru]);
199 Int_t iRun=runNumb[jru];
200 if(readRun->TestBitNumber(iRun))printf("Run %d - already processed\n",iRun);
201 if(readRun->TestBitNumber(iRun))continue;
202 //cout << "Value from file is " <<t << endl;
203
204 if(jru>0){
205 for(Int_t jru2=0;jru2<jru;jru2++) {
206 if(runNumb[jru2]==runNumb[jru]){
207 printf("Run %d - twice in the list - skipped\n",runNumb[jru]);
208 twice=1;
209 }
210 }
211 }
212 if(twice==1){twice=0;continue;}
213
214 printf("%s\n",strings[jru]);
215
216 if(!gGrid||!gGrid->IsConnected()) {
217 printf("gGrid not found! exit macro\n");
218 return;
219 }
220 TFile *f=TFile::Open(Form("alien://%s",strings[jru]));
221
222 TDirectoryFile *df=(TDirectoryFile*)f->Get("ITSAlignQA");
223 if(!df){
224 printf("Run %d ITSAlignQA MISSING -> Exit\n",iRun);
225 continue;
226 }
227 TList *l = (TList*)df->Get("clistITSAlignQA");
228 if(!df){
229 printf("Run %d clistITSAlignQA TList MISSING -> Exit\n",iRun);
230 continue;
231 }
232
233 TH1F* hev=(TH1F*)l->FindObject("hNEvents"); // histo 1d numero di eventi processati
234 Int_t nTotEvents=hev->GetBinContent(2);
235 Int_t nTrigEvents=hev->GetBinContent(3);
236 printf("Run %d Number of Events = %d Triggered=%d\n",iRun,nTotEvents,nTrigEvents);
237 if(nTotEvents==0) continue;
238
239 // working variables
240 Int_t nModmux20_1=0;
241 Int_t nModmux50_1=0;
242 Int_t nModmux100_1=0;
243 Int_t nModsigx100_1=0;
244 Int_t nModsigx200_1=0;
245 Int_t nModsigx300_1=0;
246 Int_t nModmuz50_1=0;
247 Int_t nModmuz100_1=0;
248 Int_t nModmuz300_1=0;
249 Int_t nModsigz100_1=0;
250 Int_t nModsigz300_1=0;
251 Int_t nModsigz500_1=0;
252
253 Int_t nModmux20_2=0;
254 Int_t nModmux50_2=0;
255 Int_t nModmux100_2=0;
256 Int_t nModsigx100_2=0;
257 Int_t nModsigx200_2=0;
258 Int_t nModsigx300_2=0;
259 Int_t nModmuz50_2=0;
260 Int_t nModmuz100_2=0;
261 Int_t nModmuz300_2=0;
262 Int_t nModsigz100_2=0;
263 Int_t nModsigz300_2=0;
264 Int_t nModsigz500_2=0;
265
266 Int_t nModmux20_3=0;
267 Int_t nModmux50_3=0;
268 Int_t nModmux100_3=0;
269 Int_t nModsigx100_3=0;
270 Int_t nModsigx200_3=0;
271 Int_t nModsigx300_3=0;
272 Int_t nModmuxlr20_3=0;
273 Int_t nModmuxlr50_3=0;
274 Int_t nModmuxlr100_3=0;
275 Int_t nModexcx100_3=0;
276 Int_t nModexcx200_3=0;
277 Int_t nModexcx300_3=0;
278 Int_t nModmuz50_3=0;
279 Int_t nModmuz100_3=0;
280 Int_t nModmuz300_3=0;
281 Int_t nModsigz100_3=0;
282 Int_t nModsigz300_3=0;
283 Int_t nModsigz500_3=0;
284
285 Int_t nModmux20_4=0;
286 Int_t nModmux50_4=0;
287 Int_t nModmux100_4=0;
288 Int_t nModsigx100_4=0;
289 Int_t nModsigx200_4=0;
290 Int_t nModsigx300_4=0;
291 Int_t nModmuxlr20_4=0;
292 Int_t nModmuxlr50_4=0;
293 Int_t nModmuxlr100_4=0;
294 Int_t nModexcx100_4=0;
295 Int_t nModexcx200_4=0;
296 Int_t nModexcx300_4=0;
297 Int_t nModmuz50_4=0;
298 Int_t nModmuz100_4=0;
299 Int_t nModmuz300_4=0;
300 Int_t nModsigz100_4=0;
301 Int_t nModsigz300_4=0;
302 Int_t nModsigz500_4=0;
303
304 Int_t nModmux20_5=0;
305 Int_t nModmux50_5=0;
306 Int_t nModmux100_5=0;
307 Int_t nModsigx100_5=0;
308 Int_t nModsigx200_5=0;
309 Int_t nModsigx300_5=0;
310 Int_t nModmuz50_5=0;
311 Int_t nModmuz100_5=0;
312 Int_t nModmuz300_5=0;
313 Int_t nModsigz100_5=0;
314 Int_t nModsigz300_5=0;
315 Int_t nModsigz500_5=0;
316
317 Int_t nModmux20_6=0;
318 Int_t nModmux50_6=0;
319 Int_t nModmux100_6=0;
320 Int_t nModsigx100_6=0;
321 Int_t nModsigx200_6=0;
322 Int_t nModsigx300_6=0;
323 Int_t nModmuz50_6=0;
324 Int_t nModmuz100_6=0;
325 Int_t nModmuz300_6=0;
326 Int_t nModsigz100_6=0;
327 Int_t nModsigz300_6=0;
328 Int_t nModsigz500_6=0;
329
330 // mio pezzo lettura e fit
331
332 for(Int_t nmod=0; nmod<240; nmod++){ // loop on SPD modules
333
334 TString nameResXvsX=Form("hSPDResidX%d",nmod);
335 TH2F *histoResXvsX=(TH2F*)l->FindObject(nameResXvsX);
336 TH1D *histoResX=histoResXvsX->ProjectionY();
337
338 TString nameResZvsZ=Form("hSPDResidZ%d",nmod);
339 TH2F *histoResZvsZ=(TH2F*)l->FindObject(nameResZvsZ);
340 TH1D *histoResZ=histoResZvsZ->ProjectionY();
341
342 TF1 *gg_x = new TF1("gg_x","gaus",-0.050, 0.050);
343 TF1 *gg_z = new TF1("gg_z","gaus",-0.15,0.15);
344
345 if(histoResX->GetEntries()){
346 histoResX->Fit("gg_x","RQ0");
347 if(nmod<80){ // layer 1
348 if(TMath::Abs(gg_x->GetParameter(1))*10000<5.) nModmux20_1+=1; // micron
349 if(TMath::Abs(gg_x->GetParameter(1))*10000<10.) nModmux50_1+=1;
350 if(TMath::Abs(gg_x->GetParameter(1))*10000<50.) nModmux100_1+=1;
351 if(gg_x->GetParameter(2)*10000<50.) nModsigx100_1+=1;
352 if(gg_x->GetParameter(2)*10000<70.) nModsigx200_1+=1;
353 if(gg_x->GetParameter(2)*10000<200.) nModsigx300_1+=1;
354 }
355 else { // layer 2
356 if(TMath::Abs(gg_x->GetParameter(1))*10000<10.) nModmux20_2+=1; // micron
357 if(TMath::Abs(gg_x->GetParameter(1))*10000<20.) nModmux50_2+=1;
358 if(TMath::Abs(gg_x->GetParameter(1))*10000<50.) nModmux100_2+=1;
359 if(gg_x->GetParameter(2)*10000<130.) nModsigx100_2+=1;
360 if(gg_x->GetParameter(2)*10000<160.) nModsigx200_2+=1;
361 if(gg_x->GetParameter(2)*10000<200.) nModsigx300_2+=1;
362 }
363 }
364
365
366 if(histoResZ->GetEntries()){
367 histoResZ->Fit("gg_z","RQ0");
368 if(nmod<80){ // layer 1
369 if(TMath::Abs(gg_z->GetParameter(1))*10000<10.) nModmuz50_1+=1;
370 if(TMath::Abs(gg_z->GetParameter(1))*10000<20.) nModmuz100_1+=1;
371 if(TMath::Abs(gg_z->GetParameter(1))*10000<50.) nModmuz300_1+=1;
372 if(gg_z->GetParameter(2)*10000<150.) nModsigz100_1+=1;
373 if(gg_z->GetParameter(2)*10000<160.) nModsigz300_1+=1;
374 if(gg_z->GetParameter(2)*10000<200.) nModsigz500_1+=1;
375 }
376
377 else { // layer 2
378 if(TMath::Abs(gg_z->GetParameter(1))*10000<10.) nModmuz50_2+=1;
379 if(TMath::Abs(gg_z->GetParameter(1))*10000<20.) nModmuz100_2+=1;
380 if(TMath::Abs(gg_z->GetParameter(1))*10000<50.) nModmuz300_2+=1;
381 if(gg_z->GetParameter(2)*10000<180.) nModsigz100_2+=1;
382 if(gg_z->GetParameter(2)*10000<190.) nModsigz300_2+=1;
383 if(gg_z->GetParameter(2)*10000<250.) nModsigz500_2+=1;
384 }
385 }
386
387 } // loop on SPD modules -- end
388
389 for(Int_t nmod=240; nmod<500; nmod++){ // loop on SDD modules
390
391 TString nameResXvsX=Form("hSDDResidXvsX%d",nmod);
392 TH2F *histoResXvsX=(TH2F*)l->FindObject(nameResXvsX);
393 TH1D *histoResX=histoResXvsX->ProjectionY();
394 TH1D *histoResX_l=histoResXvsX->ProjectionY("residual xloc<0",1,20);
395 TH1D *histoResX_r=histoResXvsX->ProjectionY("residual xloc>0",21,40);
396
397 TString nameResZvsZ=Form("hSDDResidZvsZ%d",nmod);
398 TH2F *histoResZvsZ=(TH2F*)l->FindObject(nameResZvsZ);
399 TH1D *histoResZ=histoResZvsZ->ProjectionY();
400
401 // TString nameResZvsX=Form("hSDDResidZvsX%d",nmod);
402 // TH2F *histoResZvsX=(TH2F*)l->FindObject(nameResZvsX);
403
404 TF1 *gg_x = new TF1("gg_x","gaus",-0.050, 0.050);
405 TF1 *gg_l = new TF1("gg_l","gaus",-0.050, 0.050);
406 TF1 *gg_r = new TF1("gg_r","gaus",-0.050, 0.050);
407 TF1 *gg_z = new TF1("gg_z","gaus",-0.15,0.15);
408
409 if(histoResX->GetEntries()){
410 histoResX->Fit("gg_x","RQ0");
411 if(nmod<324){ // layer 3
412 if(TMath::Abs(gg_x->GetParameter(1))*10000<20.) nModmux20_3+=1; // micron
413 if(TMath::Abs(gg_x->GetParameter(1))*10000<50.) nModmux50_3+=1;
414 if(TMath::Abs(gg_x->GetParameter(1))*10000<100.) nModmux100_3+=1;
415 if(gg_x->GetParameter(2)*10000<200.) nModsigx100_3+=1;
416 if(gg_x->GetParameter(2)*10000<250.) nModsigx200_3+=1;
417 if(gg_x->GetParameter(2)*10000<400.) nModsigx300_3+=1;
418 }
419 else { // layer 4
420 if(TMath::Abs(gg_x->GetParameter(1))*10000<20.) nModmux20_4+=1; // micron
421 if(TMath::Abs(gg_x->GetParameter(1))*10000<50.) nModmux50_4+=1;
422 if(TMath::Abs(gg_x->GetParameter(1))*10000<100.) nModmux100_4+=1;
423 if(gg_x->GetParameter(2)*10000<200.) nModsigx100_4+=1;
424 if(gg_x->GetParameter(2)*10000<250.) nModsigx200_4+=1;
425 if(gg_x->GetParameter(2)*10000<400.) nModsigx300_4+=1;
426 }
427 }
428
429 if(histoResX_l->GetEntries()&&histoResX_r->GetEntries()){
430 histoResX_l->Fit("gg_l","RQ0");
431 histoResX_r->Fit("gg_r","RQ0");
432 if(nmod<324){ // layer 3
433 if(TMath::Abs(gg_l->GetParameter(1)*10000-gg_r->GetParameter(1)*10000)<20.) nModmuxlr20_3+=1;
434 if(TMath::Abs(gg_l->GetParameter(1)*10000-gg_r->GetParameter(1)*10000)<50.) nModmuxlr50_3+=1;
435 if(TMath::Abs(gg_l->GetParameter(1)*10000-gg_r->GetParameter(1)*10000)<100.) nModmuxlr100_3+=1;
436 }
437 else { // layer 4
438 if(TMath::Abs(gg_l->GetParameter(1)*10000-gg_r->GetParameter(1)*10000)<20.) nModmuxlr20_4+=1;
439 if(TMath::Abs(gg_l->GetParameter(1)*10000-gg_r->GetParameter(1)*10000)<50.) nModmuxlr50_4+=1;
440 if(TMath::Abs(gg_l->GetParameter(1)*10000-gg_r->GetParameter(1)*10000)<100.) nModmuxlr100_4+=1;
441 }
442 }
443
444 if(histoResZ->GetEntries()){
445 histoResZ->Fit("gg_z","RQ0");
446 if(nmod<324){ // layer 3
447 if(TMath::Abs(gg_z->GetParameter(1))*10000<30.) nModmuz50_3+=1;
448 if(TMath::Abs(gg_z->GetParameter(1))*10000<50.) nModmuz100_3+=1;
449 if(TMath::Abs(gg_z->GetParameter(1))*10000<100.) nModmuz300_3+=1;
450 if(gg_z->GetParameter(2)*10000<150.) nModsigz100_3+=1;
451 if(gg_z->GetParameter(2)*10000<200.) nModsigz300_3+=1;
452 if(gg_z->GetParameter(2)*10000<350.) nModsigz500_3+=1;
453 }
454
455 else { // layer 4
456 if(TMath::Abs(gg_z->GetParameter(1))*10000<30.) nModmuz50_4+=1;
457 if(TMath::Abs(gg_z->GetParameter(1))*10000<50.) nModmuz100_4+=1;
458 if(TMath::Abs(gg_z->GetParameter(1))*10000<100.) nModmuz300_4+=1;
459 if(gg_z->GetParameter(2)*10000<250.) nModsigz100_4+=1;
460 if(gg_z->GetParameter(2)*10000<300.) nModsigz300_4+=1;
461 if(gg_z->GetParameter(2)*10000<450.) nModsigz500_4+=1;
462 }
463 }
464
465
466 Double_t min=-10000.0;
467 Double_t max=10000.0;
468
469 TProfile *profResXvsX = histoResXvsX->ProfileX();
470 if(profResXvsX->GetEntries()>0){
471 max=profResXvsX->GetBinContent(1);
472 min=profResXvsX->GetBinContent(1);
473 for(Int_t jj=0;jj<39;jj++){
474 if(profResXvsX->GetBinContent(jj+1)>max) max = profResXvsX->GetBinContent(jj+1);
475 if(profResXvsX->GetBinContent(jj+1)<min) min = profResXvsX->GetBinContent(jj+1);
476 }
477
478 if(nmod<324){ // layer 3
479 if(TMath::Abs(max-min)*10000.<100) nModexcx100_3+=1;
480 if(TMath::Abs(max-min)*10000.<200) nModexcx200_3+=1;
481 if(TMath::Abs(max-min)*10000.<300) nModexcx300_3+=1;
482 }
483
484 else { // layer 4
485 if(TMath::Abs(max-min)*10000.<100) nModexcx100_4+=1;
486 if(TMath::Abs(max-min)*10000.<200) nModexcx200_4+=1;
487 if(TMath::Abs(max-min)*10000.<300) nModexcx300_4+=1;
488 }
489 }
490
491 } // loop on SDD modules -- end
492
493 for(Int_t nmod=500; nmod<2198; nmod++){ // loop on SSD modules
494
495 TString nameResXvsX=Form("hSSDResidX%d",nmod);
496 TH2F *histoResXvsX=(TH2F*)l->FindObject(nameResXvsX);
497 TH1D *histoResX=histoResXvsX->ProjectionY();
498
499 TString nameResZvsZ=Form("hSSDResidZ%d",nmod);
500 TH2F *histoResZvsZ=(TH2F*)l->FindObject(nameResZvsZ);
501 TH1D *histoResZ=histoResZvsZ->ProjectionY();
502
503 TF1 *gg_x = new TF1("gg_x","gaus",-0.050, 0.050);
504 TF1 *gg_z = new TF1("gg_z","gaus",-0.15,0.15);
505
506 if(histoResX->GetEntries()){
507 histoResX->Fit("gg_x","RQ0");
508 if(nmod<1248){ // layer 5
509 if(TMath::Abs(gg_x->GetParameter(1))*10000<20.) nModmux20_5+=1; // micron
510 if(TMath::Abs(gg_x->GetParameter(1))*10000<40.) nModmux50_5+=1;
511 if(TMath::Abs(gg_x->GetParameter(1))*10000<100.) nModmux100_5+=1;
512 if(gg_x->GetParameter(2)*10000<70.) nModsigx100_5+=1;
513 if(gg_x->GetParameter(2)*10000<100.) nModsigx200_5+=1;
514 if(gg_x->GetParameter(2)*10000<300.) nModsigx300_5+=1;
515 }
516 else { // layer 6
517 if(TMath::Abs(gg_x->GetParameter(1))*10000<20.) nModmux20_6+=1; // micron
518 if(TMath::Abs(gg_x->GetParameter(1))*10000<40.) nModmux50_6+=1;
519 if(TMath::Abs(gg_x->GetParameter(1))*10000<100.) nModmux100_6+=1;
520 if(gg_x->GetParameter(2)*10000<200.) nModsigx100_6+=1;
521 if(gg_x->GetParameter(2)*10000<300.) nModsigx200_6+=1;
522 if(gg_x->GetParameter(2)*10000<400.) nModsigx300_6+=1;
523 }
524 }
525
526
527 if(histoResZ->GetEntries()){
528 histoResZ->Fit("gg_z","RQ0");
529 if(nmod<1248){ // layer 5
530 if(TMath::Abs(gg_z->GetParameter(1))*10000<50.) nModmuz50_5+=1;
531 if(TMath::Abs(gg_z->GetParameter(1))*10000<100.) nModmuz100_5+=1;
532 if(TMath::Abs(gg_z->GetParameter(1))*10000<300.) nModmuz300_5+=1;
533 if(gg_z->GetParameter(2)*10000<1000.) nModsigz100_5+=1;
534 if(gg_z->GetParameter(2)*10000<1100.) nModsigz300_5+=1;
535 if(gg_z->GetParameter(2)*10000<1500.) nModsigz500_5+=1;
536 }
537
538 else { // layer 6
539 if(TMath::Abs(gg_z->GetParameter(1))*10000<50.) nModmuz50_6+=1;
540 if(TMath::Abs(gg_z->GetParameter(1))*10000<100.) nModmuz100_6+=1;
541 if(TMath::Abs(gg_z->GetParameter(1))*10000<300.) nModmuz300_6+=1;
542 if(gg_z->GetParameter(2)*10000<1000.) nModsigz100_6+=1;
543 if(gg_z->GetParameter(2)*10000<1100.) nModsigz300_6+=1;
544 if(gg_z->GetParameter(2)*10000<1500.) nModsigz500_6+=1;
545 }
546 }
547
548 } // loop on SSD modules -- end
549
550 Int_t count1 = 0;
551 Int_t count2 = 0;
552 Int_t count3 = 0;
553 Int_t count4 = 0;
554 Int_t count5 = 0;
555 Int_t count6 = 0;
556
557// Int_t nPSideChannelsTotal = 0, nNSideChannelsTotal = 0;
558// Int_t nBadPSideChannelsTotal = 0, nBadNSideChannelsTotal = 0;
559 Int_t nBadPSideChannels = 0, nBadNSideChannels = 0;
560 Int_t layer = 0, ladder = 0, module = 0;
561// Int_t nPSideChannelsLayer5 = 0, nNSideChannelsLayer5 = 0;
562// Int_t nPSideChannelsLayer6 = 0, nNSideChannelsLayer6 = 0;
563// Int_t nPSideChannelsLayer5Total = 0, nNSideChannelsLayer5Total = 0;
564// Int_t nPSideChannelsLayer6Total = 0, nNSideChannelsLayer6Total = 0;
565// Int_t nBadPSideChannelsLayer5Total = 0, nBadNSideChannelsLayer5Total = 0;
566// Int_t nBadPSideChannelsLayer6Total = 0, nBadNSideChannelsLayer6Total = 0;
567// Int_t badPmodule5 = 0, badNmodule5 = 0, badPmodule6 = 0, badNmodule6 = 0;
568
569 // lettura moduli SDD Ok da OCDB
570 man->SetRun(iRun);
571 entrySDD = AliCDBManager::Instance()->Get("ITS/Calib/CalibSDD");
572 calSDD = (TObjArray *)entrySDD->GetObject();
573 for(Int_t i=0;i<84;i++){cal = (AliITSCalibration*)
574 calSDD->At(i); if(cal->IsBad())count3++;}
575 for(Int_t i=84;i<calSDD->GetEntries();i++){cal = (AliITSCalibration*)
576 calSDD->At(i); if(cal->IsBad())count4++;}
577// cout << endl << endl << count3 << " " << count4 << endl;
578
579 // ntuple filling
580 Int_t index=0;
581 // 0
582 xnt[index++]=(Float_t)iRun;
583 // 1 SDD L3
584 xnt[index++]=nModmux20_3/(84.-count3);
585 xnt[index++]=TMath::Sqrt(nModmux20_3)/(84.-count3);
586 xnt[index++]=nModmux50_3/(84.-count3);
587 xnt[index++]=TMath::Sqrt(nModmux50_3)/(84.-count3);
588 xnt[index++]=nModmux100_3/(84.-count3);
589 xnt[index++]=TMath::Sqrt(nModmux100_3)/(84.-count3);
590 xnt[index++]=nModsigx100_3/(84.-count3);
591 xnt[index++]=TMath::Sqrt(nModsigx100_3)/(84.-count3);
592 xnt[index++]=nModsigx200_3/(84.-count3);
593 xnt[index++]=TMath::Sqrt(nModsigx200_3)/(84.-count3);
594 xnt[index++]=nModsigx300_3/(84.-count3);
595 xnt[index++]=TMath::Sqrt(nModsigx300_3)/(84.-count3);
596 xnt[index++]=nModmuxlr20_3/(84.-count3);
597 xnt[index++]=TMath::Sqrt(nModmuxlr20_3)/(84.-count3);
598 xnt[index++]=nModmuxlr50_3/(84.-count3);
599 xnt[index++]=TMath::Sqrt(nModmuxlr50_3)/(84.-count3);
600 xnt[index++]=nModmuxlr100_3/(84.-count3);
601 xnt[index++]=TMath::Sqrt(nModmuxlr100_3)/(84.-count3);
602 xnt[index++]=nModexcx100_3/(84.-count3);
603 xnt[index++]=TMath::Sqrt(nModexcx100_3)/(84.-count3);
604 xnt[index++]=nModexcx200_3/(84.-count3);
605 xnt[index++]=TMath::Sqrt(nModexcx200_3)/(84.-count3);
606 xnt[index++]=nModexcx300_3/(84.-count3);
607 xnt[index++]=TMath::Sqrt(nModexcx300_3)/(84.-count3);
608 xnt[index++]=nModmuz50_3/(84.-count3);
609 xnt[index++]=TMath::Sqrt(nModmuz50_3)/(84.-count3);
610 xnt[index++]=nModmuz100_3/(84.-count3);
611 xnt[index++]=TMath::Sqrt(nModmuz100_3)/(84.-count3);
612 xnt[index++]=nModmuz300_3/(84.-count3);
613 xnt[index++]=TMath::Sqrt(nModmuz300_3)/(84.-count3);
614 xnt[index++]=nModsigz100_3/(84.-count3);
615 xnt[index++]=TMath::Sqrt(nModsigz100_3)/(84.-count3);
616 xnt[index++]=nModsigz300_3/(84.-count3);
617 xnt[index++]=TMath::Sqrt(nModsigz300_3)/(84.-count3);
618 xnt[index++]=nModsigz500_3/(84.-count3);
619 xnt[index++]=TMath::Sqrt(nModsigz500_3)/(84.-count3);
620 // 37 SDD L4
621 xnt[index++]=nModmux20_4/(176.-count4);
622 xnt[index++]=TMath::Sqrt(nModmux20_4)/(176.-count4);
623 xnt[index++]=nModmux50_4/(176.-count4);
624 xnt[index++]=TMath::Sqrt(nModmux50_4)/(176.-count4);
625 xnt[index++]=nModmux100_4/(176.-count4);
626 xnt[index++]=TMath::Sqrt(nModmux100_4)/(176.-count4);
627 xnt[index++]=nModsigx100_4/(176.-count4);
628 xnt[index++]=TMath::Sqrt(nModsigx100_4)/(176.-count4);
629 xnt[index++]=nModsigx200_4/(176.-count4);
630 xnt[index++]=TMath::Sqrt(nModsigx200_4)/(176.-count4);
631 xnt[index++]=nModsigx300_4/(176.-count4);
632 xnt[index++]=TMath::Sqrt(nModsigx300_4)/(176.-count4);
633 xnt[index++]=nModmuxlr20_4/(176.-count4);
634 xnt[index++]=TMath::Sqrt(nModmuxlr20_4)/(176.-count4);
635 xnt[index++]=nModmuxlr50_4/(176.-count4);
636 xnt[index++]=TMath::Sqrt(nModmuxlr50_4)/(176.-count4);
637 xnt[index++]=nModmuxlr100_4/(176.-count4);
638 xnt[index++]=TMath::Sqrt(nModmuxlr100_4)/(176.-count4);
639 xnt[index++]=nModexcx100_4/(176.-count4);
640 xnt[index++]=TMath::Sqrt(nModexcx100_4)/(176.-count4);
641 xnt[index++]=nModexcx200_4/(176.-count4);
642 xnt[index++]=TMath::Sqrt(nModexcx200_4)/(176.-count4);
643 xnt[index++]=nModexcx300_4/(176.-count4);
644 xnt[index++]=TMath::Sqrt(nModexcx300_4)/(176.-count4);
645 xnt[index++]=nModmuz50_4/(176.-count4);
646 xnt[index++]=TMath::Sqrt(nModmuz50_4)/(176.-count4);
647 xnt[index++]=nModmuz100_4/(176.-count4);
648 xnt[index++]=TMath::Sqrt(nModmuz100_4)/(176.-count4);
649 xnt[index++]=nModmuz300_4/(176.-count4);
650 xnt[index++]=TMath::Sqrt(nModmuz300_4)/(176.-count4);
651 xnt[index++]=nModsigz100_4/(176.-count4);
652 xnt[index++]=TMath::Sqrt(nModsigz100_4)/(176.-count4);
653 xnt[index++]=nModsigz300_4/(176.-count4);
654 xnt[index++]=TMath::Sqrt(nModsigz300_4)/(176.-count4);
655 xnt[index++]=nModsigz500_4/(176.-count4);
656 xnt[index++]=TMath::Sqrt(nModsigz500_4)/(176.-count4);
657 // 73 SPD L1
658 // lettura moduli SPD bad da OCDB
659 AliITSOnlineCalibrationSPDhandler *h = new AliITSOnlineCalibrationSPDhandler();
660 h->ReadDeadFromDB(iRun,"alien://folder=/alice/data/2011/OCDB");
661// Int_t nDeadModules[2]={0,0}; // n dead inner and outer
662// Double_t fractions[2]={0.,0.}; // fraction of dead modules
663 // # dead modules inner layer
664 for(Int_t i=0; i<80; i++) count1+=h->GetNrBad(i)/(5*8192);
665 // # dead modules outer layer
666 for(Int_t jj=80; jj<240; jj++) count2+=h->GetNrBad(jj)/(5*8192);
667
668 xnt[index++]=nModmux20_1/(80.-count1);
669 xnt[index++]=TMath::Sqrt(nModmux20_1)/(80.-count1);
670 xnt[index++]=nModmux50_1/(80.-count1);
671 xnt[index++]=TMath::Sqrt(nModmux50_1)/(80.-count1);
672 xnt[index++]=nModmux100_1/(80.-count1);
673 xnt[index++]=TMath::Sqrt(nModmux100_1)/(80.-count1);
674 xnt[index++]=nModsigx100_1/(80.-count1);
675 xnt[index++]=TMath::Sqrt(nModsigx100_1)/(80.-count1);
676 xnt[index++]=nModsigx200_1/(80.-count1);
677 xnt[index++]=TMath::Sqrt(nModsigx200_1)/(80.-count1);
678 xnt[index++]=nModsigx300_1/(80.-count1);
679 xnt[index++]=TMath::Sqrt(nModsigx300_1)/(80.-count1);
680 xnt[index++]=nModmuz50_1/(80.-count1);
681 xnt[index++]=TMath::Sqrt(nModmuz50_1)/(80.-count1);
682 xnt[index++]=nModmuz100_1/(80.-count1);
683 xnt[index++]=TMath::Sqrt(nModmuz100_1)/(80.-count1);
684 xnt[index++]=nModmuz300_1/(80.-count1);
685 xnt[index++]=TMath::Sqrt(nModmuz300_1)/(80.-count1);
686 xnt[index++]=nModsigz100_1/(80.-count1);
687 xnt[index++]=TMath::Sqrt(nModsigz100_1)/(80.-count1);
688 xnt[index++]=nModsigz300_1/(80.-count1);
689 xnt[index++]=TMath::Sqrt(nModsigz300_1)/(80.-count1);
690 xnt[index++]=nModsigz500_1/(80.-count1);
691 xnt[index++]=TMath::Sqrt(nModsigz500_1)/(80.-count1);
692 // 97 SPD L2
693 xnt[index++]=nModmux20_2/(160.-count2);
694 xnt[index++]=TMath::Sqrt(nModmux20_2)/(160.-count2);
695 xnt[index++]=nModmux50_2/(160.-count2);
696 xnt[index++]=TMath::Sqrt(nModmux50_2)/(160.-count2);
697 xnt[index++]=nModmux100_2/(160.-count2);
698 xnt[index++]=TMath::Sqrt(nModmux100_2)/(160.-count2);
699 xnt[index++]=nModsigx100_2/(160.-count2);
700 xnt[index++]=TMath::Sqrt(nModsigx100_2)/(160.-count2);
701 xnt[index++]=nModsigx200_2/(160.-count2);
702 xnt[index++]=TMath::Sqrt(nModsigx200_2)/(160.-count2);
703 xnt[index++]=nModsigx300_2/(160.-count2);
704 xnt[index++]=TMath::Sqrt(nModsigx300_2)/(160.-count2);
705 xnt[index++]=nModmuz50_2/(160.-count2);
706 xnt[index++]=TMath::Sqrt(nModmuz50_2)/(160.-count2);
707 xnt[index++]=nModmuz100_2/(160.-count2);
708 xnt[index++]=TMath::Sqrt(nModmuz100_2)/(160.-count2);
709 xnt[index++]=nModmuz300_2/(160.-count2);
710 xnt[index++]=TMath::Sqrt(nModmuz300_2)/(160.-count2);
711 xnt[index++]=nModsigz100_2/(160.-count2);
712 xnt[index++]=TMath::Sqrt(nModsigz100_2)/(160.-count2);
713 xnt[index++]=nModsigz300_2/(160.-count2);
714 xnt[index++]=TMath::Sqrt(nModsigz300_2)/(160.-count2);
715 xnt[index++]=nModsigz500_2/(160.-count2);
716 xnt[index++]=TMath::Sqrt(nModsigz500_2)/(160.-count2);
717 // 121 SSD L5
718 // lettura moduli SSD bad da OCDB
719 const Int_t fgkSSDMODULES = 1698;
720 static const Int_t fgkDefaultNStripsSSD = 768;
721 AliITSBadChannelsSSDv2 *badChannelsSSD = new AliITSBadChannelsSSDv2();
722 AliCDBEntry *entryBadChannelsSSD = man->Get("ITS/Calib/BadChannelsSSD");
723 TObject *empty = (TObject *)entryBadChannelsSSD->GetObject();
724 TString objectname = empty->GetName();
725 if(objectname=="TObjArray") {
726 TObjArray *badChannelsSSDOld = (TObjArray *)entryBadChannelsSSD->GetObject();
727 ReadOldSSDBadChannels(badChannelsSSDOld, badChannelsSSD);
728 }
729 else if(objectname=="AliITSBadChannelsSSDv2") {
730 cout<<"Reading the new format of the calibration file..."<<endl;
731 badChannelsSSD = (AliITSBadChannelsSSDv2 *)entryBadChannelsSSD->GetObject();
732 }
733 nBadPSideChannels = 0;
734 nBadNSideChannels = 0;
735 layer = 0;
736 ladder = 0;
737 module = 0;
738 for(Int_t i = 0; i < fgkSSDMODULES; i++) { // loop sui moduli
739 AliITSgeomTGeo::GetModuleId(i+500,layer,ladder,module);
740 nBadPSideChannels = 0, nBadNSideChannels = 0;
741
742 Int_t badChannel = 0;
743 for(Int_t ja = 0; ja < fgkDefaultNStripsSSD; ja++) { // loop sulle strips del modulo
744 badChannel = (Int_t)(badChannelsSSD->GetBadChannelP(i,ja)); // bad channels P side
745 //cout<<"Module: "<<i+500<< " Strip: "<<j<<" - "<<badChannel<<endl;
746 if(badChannel != 0) nBadPSideChannels += 1;
747
748 badChannel = (Int_t)(badChannelsSSD->GetBadChannelN(i,ja)); // bad channels N side
749 //cout<<"Module: "<<i+500<< " Strip: "<<fgkDefaultNStripsSSD+j+1<<" - "<<badChannel<<endl;
750 if(badChannel != 0) nBadNSideChannels += 1;
751 } // end loop sulle strip del modulo
752
753 // determinazione numero moduli morti sui due lati per ogni layer
754 if(nBadNSideChannels == 768 && nBadPSideChannels == 768){
755 if(layer == 5){
756 count5++;
757 }
758 else {
759 count6++;
760 }
761 }
762 } // endl loop sui moduli
763
764 xnt[index++]=nModmux20_5/(748.-count5);
765 xnt[index++]=TMath::Sqrt(nModmux20_5)/(748.-count5);
766 xnt[index++]=nModmux50_5/(748.-count5);
767 xnt[index++]=TMath::Sqrt(nModmux50_5)/(748.-count5);
768 xnt[index++]=nModmux100_5/(748.-count5);
769 xnt[index++]=TMath::Sqrt(nModmux100_5)/(748.-count5);
770 xnt[index++]=nModsigx100_5/(748.-count5);
771 xnt[index++]=TMath::Sqrt(nModsigx100_5)/(748.-count5);
772 xnt[index++]=nModsigx200_5/(748.-count5);
773 xnt[index++]=TMath::Sqrt(nModsigx200_5)/(748.-count5);
774 xnt[index++]=nModsigx300_5/(748.-count5);
775 xnt[index++]=TMath::Sqrt(nModsigx300_5)/(748.-count5);
776 xnt[index++]=nModmuz50_5/(748.-count5);
777 xnt[index++]=TMath::Sqrt(nModmuz50_5)/(748.-count5);
778 xnt[index++]=nModmuz100_5/(748.-count5);
779 xnt[index++]=TMath::Sqrt(nModmuz100_5)/(748.-count5);
780 xnt[index++]=nModmuz300_5/(748.-count5);
781 xnt[index++]=TMath::Sqrt(nModmuz300_5)/(748.-count5);
782 xnt[index++]=nModsigz100_5/(748.-count5);
783 xnt[index++]=TMath::Sqrt(nModsigz100_5)/(748.-count5);
784 xnt[index++]=nModsigz300_5/(748.-count5);
785 xnt[index++]=TMath::Sqrt(nModsigz300_5)/(748.-count5);
786 xnt[index++]=nModsigz500_5/(748.-count5);
787 xnt[index++]=TMath::Sqrt(nModsigz500_5)/(748.-count5);
788 // 145 SSD L6
789 xnt[index++]=nModmux20_6/(950.-count6);
790 xnt[index++]=TMath::Sqrt(nModmux20_6)/(950.-count6);
791 xnt[index++]=nModmux50_6/(950.-count6);
792 xnt[index++]=TMath::Sqrt(nModmux50_6)/(950.-count6);
793 xnt[index++]=nModmux100_6/(950.-count6);
794 xnt[index++]=TMath::Sqrt(nModmux100_6)/(950.-count6);
795 xnt[index++]=nModsigx100_6/(950.-count6);
796 xnt[index++]=TMath::Sqrt(nModsigx100_6)/(950.-count6);
797 xnt[index++]=nModsigx200_6/(950.-count6);
798 xnt[index++]=TMath::Sqrt(nModsigx200_6)/(950.-count6);
799 xnt[index++]=nModsigx300_6/(950.-count6);
800 xnt[index++]=TMath::Sqrt(nModsigx300_6)/(950.-count6);
801 xnt[index++]=nModmuz50_6/(950.-count6);
802 xnt[index++]=TMath::Sqrt(nModmuz50_6)/(950.-count6);
803 xnt[index++]=nModmuz100_6/(950.-count6);
804 xnt[index++]=TMath::Sqrt(nModmuz100_6)/(950.-count6);
805 xnt[index++]=nModmuz300_6/(950.-count6);
806 xnt[index++]=TMath::Sqrt(nModmuz300_6)/(950.-count6);
807 xnt[index++]=nModsigz100_6/(950.-count6);
808 xnt[index++]=TMath::Sqrt(nModsigz100_6)/(950.-count6);
809 xnt[index++]=nModsigz300_6/(950.-count6);
810 xnt[index++]=TMath::Sqrt(nModsigz300_6)/(950.-count6);
811 xnt[index++]=nModsigz500_6/(950.-count6);
812 xnt[index++]=TMath::Sqrt(nModsigz500_6)/(950.-count6);
813 // 168 - fine SSD L6
814 ntits->Fill(xnt);
815 // cout<<"\n\nirun sDD"<<iRun<<endl<<endl;
816
817 } // loop on runs
818
819 TFile* outfil=new TFile(ntupleFileName.Data(),"recreate");
820 outfil->cd();
821 ntits->Write();
822 outfil->Close();
823 delete outfil;
824
825 // MakePlots(ntupleFileName);
826}
827
828
829
830void MakePlots(Int_t run1,Int_t run2,TString ntupleFileName){
831 if(run1>=run2){
832 printf("****** ERROR: invalid run range %d - %d\n",run1,run2);
833 return;
834 }
835 TFile* fil=new TFile(ntupleFileName.Data(),"read");
836 if(!fil){
837 printf("File with ntuple does not exist\n");
838 return;
839 }
840 TNtuple* ntits=(TNtuple*)fil->Get("ntitsalign");
841
842 Float_t nrun;
843
844 Float_t SDDfracModmux20_3,SDDfracModmux50_3,SDDfracModmux100_3,SDDfracModsigx100_3,SDDfracModsigx200_3,SDDfracModsigx300_3;
845 Float_t SDDfracModmuxlr20_3,SDDfracModmuxlr50_3,SDDfracModmuxlr100_3,SDDfracModexcx100_3,SDDfracModexcx200_3,SDDfracModexcx300_3;
846 Float_t SDDfracModmuz50_3,SDDfracModmuz100_3,SDDfracModmuz300_3,SDDfracModsigz100_3,SDDfracModsigz300_3,SDDfracModsigz500_3;
847 Float_t errSDDfracModmux20_3,errSDDfracModmux50_3,errSDDfracModmux100_3,errSDDfracModsigx100_3,errSDDfracModsigx200_3,errSDDfracModsigx300_3;
848 Float_t errSDDfracModmuxlr20_3,errSDDfracModmuxlr50_3,errSDDfracModmuxlr100_3,errSDDfracModexcx100_3,errSDDfracModexcx200_3,errSDDfracModexcx300_3;
849 Float_t errSDDfracModmuz50_3,errSDDfracModmuz100_3,errSDDfracModmuz300_3,errSDDfracModsigz100_3,errSDDfracModsigz300_3,errSDDfracModsigz500_3;
850
851 Float_t SDDfracModmux20_4,SDDfracModmux50_4,SDDfracModmux100_4,SDDfracModsigx100_4,SDDfracModsigx200_4,SDDfracModsigx300_4;
852 Float_t SDDfracModmuxlr20_4,SDDfracModmuxlr50_4,SDDfracModmuxlr100_4,SDDfracModexcx100_4,SDDfracModexcx200_4,SDDfracModexcx300_4;
853 Float_t SDDfracModmuz50_4,SDDfracModmuz100_4,SDDfracModmuz300_4,SDDfracModsigz100_4,SDDfracModsigz300_4,SDDfracModsigz500_4;
854 Float_t errSDDfracModmux20_4,errSDDfracModmux50_4,errSDDfracModmux100_4,errSDDfracModsigx100_4,errSDDfracModsigx200_4,errSDDfracModsigx300_4;
855 Float_t errSDDfracModmuxlr20_4,errSDDfracModmuxlr50_4,errSDDfracModmuxlr100_4,errSDDfracModexcx100_4,errSDDfracModexcx200_4,errSDDfracModexcx300_4;
856 Float_t errSDDfracModmuz50_4,errSDDfracModmuz100_4,errSDDfracModmuz300_4,errSDDfracModsigz100_4,errSDDfracModsigz300_4,errSDDfracModsigz500_4;
857
858 Float_t SPDfracModmux20_1,SPDfracModmux50_1,SPDfracModmux100_1,SPDfracModsigx100_1,SPDfracModsigx200_1,SPDfracModsigx300_1;
859 Float_t SPDfracModmuz50_1,SPDfracModmuz100_1,SPDfracModmuz300_1,SPDfracModsigz100_1,SPDfracModsigz300_1,SPDfracModsigz500_1;
860 Float_t errSPDfracModmux20_1,errSPDfracModmux50_1,errSPDfracModmux100_1,errSPDfracModsigx100_1,errSPDfracModsigx200_1,errSPDfracModsigx300_1;
861 Float_t errSPDfracModmuz50_1,errSPDfracModmuz100_1,errSPDfracModmuz300_1,errSPDfracModsigz100_1,errSPDfracModsigz300_1,errSPDfracModsigz500_1;
862 Float_t SPDfracModmux20_2,SPDfracModmux50_2,SPDfracModmux100_2,SPDfracModsigx100_2,SPDfracModsigx200_2,SPDfracModsigx300_2;
863 Float_t SPDfracModmuz50_2,SPDfracModmuz100_2,SPDfracModmuz300_2,SPDfracModsigz100_2,SPDfracModsigz300_2,SPDfracModsigz500_2;
864 Float_t errSPDfracModmux20_2,errSPDfracModmux50_2,errSPDfracModmux100_2,errSPDfracModsigx100_2,errSPDfracModsigx200_2,errSPDfracModsigx300_2;
865 Float_t errSPDfracModmuz50_2,errSPDfracModmuz100_2,errSPDfracModmuz300_2,errSPDfracModsigz100_2,errSPDfracModsigz300_2,errSPDfracModsigz500_2;
866
867 Float_t SSDfracModmux20_5,SSDfracModmux50_5,SSDfracModmux100_5,SSDfracModsigx100_5,SSDfracModsigx200_5,SSDfracModsigx300_5;
868 Float_t SSDfracModmuz50_5,SSDfracModmuz100_5,SSDfracModmuz300_5,SSDfracModsigz100_5,SSDfracModsigz300_5,SSDfracModsigz500_5;
869 Float_t errSSDfracModmux20_5,errSSDfracModmux50_5,errSSDfracModmux100_5,errSSDfracModsigx100_5,errSSDfracModsigx200_5,errSSDfracModsigx300_5;
870 Float_t errSSDfracModmuz50_5,errSSDfracModmuz100_5,errSSDfracModmuz300_5,errSSDfracModsigz100_5,errSSDfracModsigz300_5,errSSDfracModsigz500_5;
871 Float_t SSDfracModmux20_6,SSDfracModmux50_6,SSDfracModmux100_6,SSDfracModsigx100_6,SSDfracModsigx200_6,SSDfracModsigx300_6;
872 Float_t SSDfracModmuz50_6,SSDfracModmuz100_6,SSDfracModmuz300_6,SSDfracModsigz100_6,SSDfracModsigz300_6,SSDfracModsigz500_6;
873 Float_t errSSDfracModmux20_6,errSSDfracModmux50_6,errSSDfracModmux100_6,errSSDfracModsigx100_6,errSSDfracModsigx200_6,errSSDfracModsigx300_6;
874 Float_t errSSDfracModmuz50_6,errSSDfracModmuz100_6,errSSDfracModmuz300_6,errSSDfracModsigz100_6,errSSDfracModsigz300_6,errSSDfracModsigz500_6;
875
876 ntits->SetBranchAddress("nrun",&nrun); // SDD
877 ntits->SetBranchAddress("SDDfracModmux20_3",&SDDfracModmux20_3);
878 ntits->SetBranchAddress("errSDDfracModmux20_3",&errSDDfracModmux20_3);
879 ntits->SetBranchAddress("SDDfracModmux50_3",&SDDfracModmux50_3);
880 ntits->SetBranchAddress("errSDDfracModmux50_3",&errSDDfracModmux50_3);
881 ntits->SetBranchAddress("SDDfracModmux100_3",&SDDfracModmux100_3);
882 ntits->SetBranchAddress("errSDDfracModmux100_3",&errSDDfracModmux100_3);
883 ntits->SetBranchAddress("SDDfracModsigx100_3",&SDDfracModsigx100_3);
884 ntits->SetBranchAddress("errSDDfracModsigx100_3",&errSDDfracModsigx100_3);
885 ntits->SetBranchAddress("SDDfracModsigx200_3",&SDDfracModsigx200_3);
886 ntits->SetBranchAddress("errSDDfracModsigx200_3",&errSDDfracModsigx200_3);
887 ntits->SetBranchAddress("SDDfracModsigx300_3",&SDDfracModsigx300_3);
888 ntits->SetBranchAddress("errSDDfracModsigx300_3",&errSDDfracModsigx300_3);
889 ntits->SetBranchAddress("SDDfracModmuxlr20_3",&SDDfracModmuxlr20_3);
890 ntits->SetBranchAddress("errSDDfracModmuxlr20_3",&errSDDfracModmuxlr20_3);
891 ntits->SetBranchAddress("SDDfracModmuxlr50_3",&SDDfracModmuxlr50_3);
892 ntits->SetBranchAddress("errSDDfracModmuxlr50_3",&errSDDfracModmuxlr50_3);
893 ntits->SetBranchAddress("SDDfracModmuxlr100_3",&SDDfracModmuxlr100_3);
894 ntits->SetBranchAddress("errSDDfracModmuxlr100_3",&errSDDfracModmuxlr100_3);
895 ntits->SetBranchAddress("SDDfracModexcx100_3",&SDDfracModexcx100_3);
896 ntits->SetBranchAddress("errSDDfracModexcx100_3",&errSDDfracModexcx100_3);
897 ntits->SetBranchAddress("SDDfracModexcx200_3",&SDDfracModexcx200_3);
898 ntits->SetBranchAddress("errSDDfracModexcx200_3",&errSDDfracModexcx200_3);
899 ntits->SetBranchAddress("SDDfracModexcx300_3",&SDDfracModexcx300_3);
900 ntits->SetBranchAddress("errSDDfracModexcx300_3",&errSDDfracModexcx300_3);
901 ntits->SetBranchAddress("SDDfracModmuz50_3",&SDDfracModmuz50_3);
902 ntits->SetBranchAddress("errSDDfracModmuz50_3",&errSDDfracModmuz50_3);
903 ntits->SetBranchAddress("SDDfracModmuz100_3",&SDDfracModmuz100_3);
904 ntits->SetBranchAddress("errSDDfracModmuz100_3",&errSDDfracModmuz100_3);
905 ntits->SetBranchAddress("SDDfracModmuz300_3",&SDDfracModmuz300_3);
906 ntits->SetBranchAddress("errSDDfracModmuz300_3",&errSDDfracModmuz300_3);
907 ntits->SetBranchAddress("SDDfracModsigz100_3",&SDDfracModsigz100_3);
908 ntits->SetBranchAddress("errSDDfracModsigz100_3",&errSDDfracModsigz100_3);
909 ntits->SetBranchAddress("SDDfracModsigz300_3",&SDDfracModsigz300_3);
910 ntits->SetBranchAddress("errSDDfracModsigz300_3",&errSDDfracModsigz300_3);
911 ntits->SetBranchAddress("SDDfracModsigz500_3",&SDDfracModsigz500_3);
912 ntits->SetBranchAddress("errSDDfracModsigz500_3",&errSDDfracModsigz500_3);
913 ntits->SetBranchAddress("SDDfracModmux20_4",&SDDfracModmux20_4);
914 ntits->SetBranchAddress("errSDDfracModmux20_4",&errSDDfracModmux20_4);
915 ntits->SetBranchAddress("SDDfracModmux50_4",&SDDfracModmux50_4);
916 ntits->SetBranchAddress("errSDDfracModmux50_4",&errSDDfracModmux50_4);
917 ntits->SetBranchAddress("SDDfracModmux100_4",&SDDfracModmux100_4);
918 ntits->SetBranchAddress("errSDDfracModmux100_4",&errSDDfracModmux100_4);
919 ntits->SetBranchAddress("SDDfracModsigx100_4",&SDDfracModsigx100_4);
920 ntits->SetBranchAddress("errSDDfracModsigx100_4",&errSDDfracModsigx100_4);
921 ntits->SetBranchAddress("SDDfracModsigx200_4",&SDDfracModsigx200_4);
922 ntits->SetBranchAddress("errSDDfracModsigx200_4",&errSDDfracModsigx200_4);
923 ntits->SetBranchAddress("SDDfracModsigx300_4",&SDDfracModsigx300_4);
924 ntits->SetBranchAddress("errSDDfracModsigx300_4",&errSDDfracModsigx300_4);
925 ntits->SetBranchAddress("SDDfracModmuxlr20_4",&SDDfracModmuxlr20_4);
926 ntits->SetBranchAddress("errSDDfracModmuxlr20_4",&errSDDfracModmuxlr20_4);
927 ntits->SetBranchAddress("SDDfracModmuxlr50_4",&SDDfracModmuxlr50_4);
928 ntits->SetBranchAddress("errSDDfracModmuxlr50_4",&errSDDfracModmuxlr50_4);
929 ntits->SetBranchAddress("SDDfracModmuxlr100_4",&SDDfracModmuxlr100_4);
930 ntits->SetBranchAddress("errSDDfracModmuxlr100_4",&errSDDfracModmuxlr100_4);
931 ntits->SetBranchAddress("SDDfracModexcx100_4",&SDDfracModexcx100_4);
932 ntits->SetBranchAddress("errSDDfracModexcx100_4",&errSDDfracModexcx100_4);
933 ntits->SetBranchAddress("SDDfracModexcx200_4",&SDDfracModexcx200_4);
934 ntits->SetBranchAddress("errSDDfracModexcx200_4",&errSDDfracModexcx200_4);
935 ntits->SetBranchAddress("SDDfracModexcx300_4",&SDDfracModexcx300_4);
936 ntits->SetBranchAddress("errSDDfracModexcx300_4",&errSDDfracModexcx300_4);
937 ntits->SetBranchAddress("SDDfracModmuz50_4",&SDDfracModmuz50_4);
938 ntits->SetBranchAddress("errSDDfracModmuz50_4",&errSDDfracModmuz50_4);
939 ntits->SetBranchAddress("SDDfracModmuz100_4",&SDDfracModmuz100_4);
940 ntits->SetBranchAddress("errSDDfracModmuz100_4",&errSDDfracModmuz100_4);
941 ntits->SetBranchAddress("SDDfracModmuz300_4",&SDDfracModmuz300_4);
942 ntits->SetBranchAddress("errSDDfracModmuz300_4",&errSDDfracModmuz300_4);
943 ntits->SetBranchAddress("SDDfracModsigz100_4",&SDDfracModsigz100_4);
944 ntits->SetBranchAddress("errSDDfracModsigz100_4",&errSDDfracModsigz100_4);
945 ntits->SetBranchAddress("SDDfracModsigz300_4",&SDDfracModsigz300_4);
946 ntits->SetBranchAddress("errSDDfracModsigz300_4",&errSDDfracModsigz300_4);
947 ntits->SetBranchAddress("SDDfracModsigz500_4",&SDDfracModsigz500_4);
948 ntits->SetBranchAddress("errSDDfracModsigz500_4",&errSDDfracModsigz500_4);
949
950 ntits->SetBranchAddress("SPDfracModmux20_1",&SPDfracModmux20_1); // SPD
951 ntits->SetBranchAddress("errSPDfracModmux20_1",&errSPDfracModmux20_1);
952 ntits->SetBranchAddress("SPDfracModmux50_1",&SPDfracModmux50_1);
953 ntits->SetBranchAddress("errSPDfracModmux50_1",&errSPDfracModmux50_1);
954 ntits->SetBranchAddress("SPDfracModmux100_1",&SPDfracModmux100_1);
955 ntits->SetBranchAddress("errSPDfracModmux100_1",&errSPDfracModmux100_1);
956 ntits->SetBranchAddress("SPDfracModsigx100_1",&SPDfracModsigx100_1);
957 ntits->SetBranchAddress("errSPDfracModsigx100_1",&errSPDfracModsigx100_1);
958 ntits->SetBranchAddress("SPDfracModsigx200_1",&SPDfracModsigx200_1);
959 ntits->SetBranchAddress("errSPDfracModsigx200_1",&errSPDfracModsigx200_1);
960 ntits->SetBranchAddress("SPDfracModsigx300_1",&SPDfracModsigx300_1);
961 ntits->SetBranchAddress("errSPDfracModsigx300_1",&errSPDfracModsigx300_1);
962 ntits->SetBranchAddress("SPDfracModmuz50_1",&SPDfracModmuz50_1);
963 ntits->SetBranchAddress("errSPDfracModmuz50_1",&errSPDfracModmuz50_1);
964 ntits->SetBranchAddress("SPDfracModmuz100_1",&SPDfracModmuz100_1);
965 ntits->SetBranchAddress("errSPDfracModmuz100_1",&errSPDfracModmuz100_1);
966 ntits->SetBranchAddress("SPDfracModmuz300_1",&SPDfracModmuz300_1);
967 ntits->SetBranchAddress("errSPDfracModmuz300_1",&errSPDfracModmuz300_1);
968 ntits->SetBranchAddress("SPDfracModsigz100_1",&SPDfracModsigz100_1);
969 ntits->SetBranchAddress("errSPDfracModsigz100_1",&errSPDfracModsigz100_1);
970 ntits->SetBranchAddress("SPDfracModsigz300_1",&SPDfracModsigz300_1);
971 ntits->SetBranchAddress("errSPDfracModsigz300_1",&errSPDfracModsigz300_1);
972 ntits->SetBranchAddress("SPDfracModsigz500_1",&SPDfracModsigz500_1);
973 ntits->SetBranchAddress("errSPDfracModsigz500_1",&errSPDfracModsigz500_1);
974 ntits->SetBranchAddress("SPDfracModmux20_2",&SPDfracModmux20_2);
975 ntits->SetBranchAddress("errSPDfracModmux20_2",&errSPDfracModmux20_2);
976 ntits->SetBranchAddress("SPDfracModmux50_2",&SPDfracModmux50_2);
977 ntits->SetBranchAddress("errSPDfracModmux50_2",&errSPDfracModmux50_2);
978 ntits->SetBranchAddress("SPDfracModmux100_2",&SPDfracModmux100_2);
979 ntits->SetBranchAddress("errSPDfracModmux100_2",&errSPDfracModmux100_2);
980 ntits->SetBranchAddress("SPDfracModsigx100_2",&SPDfracModsigx100_2);
981 ntits->SetBranchAddress("errSPDfracModsigx100_2",&errSPDfracModsigx100_2);
982 ntits->SetBranchAddress("SPDfracModsigx200_2",&SPDfracModsigx200_2);
983 ntits->SetBranchAddress("errSPDfracModsigx200_2",&errSPDfracModsigx200_2);
984 ntits->SetBranchAddress("SPDfracModsigx300_2",&SPDfracModsigx300_2);
985 ntits->SetBranchAddress("errSPDfracModsigx300_2",&errSPDfracModsigx300_2);
986 ntits->SetBranchAddress("SPDfracModmuz50_2",&SPDfracModmuz50_2);
987 ntits->SetBranchAddress("errSPDfracModmuz50_2",&errSPDfracModmuz50_2);
988 ntits->SetBranchAddress("SPDfracModmuz100_2",&SPDfracModmuz100_2);
989 ntits->SetBranchAddress("errSPDfracModmuz100_2",&errSPDfracModmuz100_2);
990 ntits->SetBranchAddress("SPDfracModmuz300_2",&SPDfracModmuz300_2);
991 ntits->SetBranchAddress("errSPDfracModmuz300_2",&errSPDfracModmuz300_2);
992 ntits->SetBranchAddress("SPDfracModsigz100_2",&SPDfracModsigz100_2);
993 ntits->SetBranchAddress("errSPDfracModsigz100_2",&errSPDfracModsigz100_2);
994 ntits->SetBranchAddress("SPDfracModsigz300_2",&SPDfracModsigz300_2);
995 ntits->SetBranchAddress("errSPDfracModsigz300_2",&errSPDfracModsigz300_2);
996 ntits->SetBranchAddress("SPDfracModsigz500_2",&SPDfracModsigz500_2);
997 ntits->SetBranchAddress("errSPDfracModsigz500_2",&errSPDfracModsigz500_2);
998
999 ntits->SetBranchAddress("SSDfracModmux20_5",&SSDfracModmux20_5); // SSD
1000 ntits->SetBranchAddress("errSSDfracModmux20_5",&errSSDfracModmux20_5);
1001 ntits->SetBranchAddress("SSDfracModmux50_5",&SSDfracModmux50_5);
1002 ntits->SetBranchAddress("errSSDfracModmux50_5",&errSSDfracModmux50_5);
1003 ntits->SetBranchAddress("SSDfracModmux100_5",&SSDfracModmux100_5);
1004 ntits->SetBranchAddress("errSSDfracModmux100_5",&errSSDfracModmux100_5);
1005 ntits->SetBranchAddress("SSDfracModsigx100_5",&SSDfracModsigx100_5);
1006 ntits->SetBranchAddress("errSSDfracModsigx100_5",&errSSDfracModsigx100_5);
1007 ntits->SetBranchAddress("SSDfracModsigx200_5",&SSDfracModsigx200_5);
1008 ntits->SetBranchAddress("errSSDfracModsigx200_5",&errSSDfracModsigx200_5);
1009 ntits->SetBranchAddress("SSDfracModsigx300_5",&SSDfracModsigx300_5);
1010 ntits->SetBranchAddress("errSSDfracModsigx300_5",&errSSDfracModsigx300_5);
1011 ntits->SetBranchAddress("SSDfracModmuz50_5",&SSDfracModmuz50_5);
1012 ntits->SetBranchAddress("errSSDfracModmuz50_5",&errSSDfracModmuz50_5);
1013 ntits->SetBranchAddress("SSDfracModmuz100_5",&SSDfracModmuz100_5);
1014 ntits->SetBranchAddress("errSSDfracModmuz100_5",&errSSDfracModmuz100_5);
1015 ntits->SetBranchAddress("SSDfracModmuz300_5",&SSDfracModmuz300_5);
1016 ntits->SetBranchAddress("errSSDfracModmuz300_5",&errSSDfracModmuz300_5);
1017 ntits->SetBranchAddress("SSDfracModsigz100_5",&SSDfracModsigz100_5);
1018 ntits->SetBranchAddress("errSSDfracModsigz100_5",&errSSDfracModsigz100_5);
1019 ntits->SetBranchAddress("SSDfracModsigz300_5",&SSDfracModsigz300_5);
1020 ntits->SetBranchAddress("errSSDfracModsigz300_5",&errSSDfracModsigz300_5);
1021 ntits->SetBranchAddress("SSDfracModsigz500_5",&SSDfracModsigz500_5);
1022 ntits->SetBranchAddress("errSSDfracModsigz500_5",&errSSDfracModsigz500_5);
1023 ntits->SetBranchAddress("SSDfracModmux20_6",&SSDfracModmux20_6);
1024 ntits->SetBranchAddress("errSSDfracModmux20_6",&errSSDfracModmux20_6);
1025 ntits->SetBranchAddress("SSDfracModmux50_6",&SSDfracModmux50_6);
1026 ntits->SetBranchAddress("errSSDfracModmux50_6",&errSSDfracModmux50_6);
1027 ntits->SetBranchAddress("SSDfracModmux100_6",&SSDfracModmux100_6);
1028 ntits->SetBranchAddress("errSSDfracModmux100_6",&errSSDfracModmux100_6);
1029 ntits->SetBranchAddress("SSDfracModsigx100_6",&SSDfracModsigx100_6);
1030 ntits->SetBranchAddress("errSSDfracModsigx100_6",&errSSDfracModsigx100_6);
1031 ntits->SetBranchAddress("SSDfracModsigx200_6",&SSDfracModsigx200_6);
1032 ntits->SetBranchAddress("errSSDfracModsigx200_6",&errSSDfracModsigx200_6);
1033 ntits->SetBranchAddress("SSDfracModsigx300_6",&SSDfracModsigx300_6);
1034 ntits->SetBranchAddress("errSSDfracModsigx300_6",&errSSDfracModsigx300_6);
1035 ntits->SetBranchAddress("SSDfracModmuz50_6",&SSDfracModmuz50_6);
1036 ntits->SetBranchAddress("errSSDfracModmuz50_6",&errSSDfracModmuz50_6);
1037 ntits->SetBranchAddress("SSDfracModmuz100_6",&SSDfracModmuz100_6);
1038 ntits->SetBranchAddress("errSSDfracModmuz100_6",&errSSDfracModmuz100_6);
1039 ntits->SetBranchAddress("SSDfracModmuz300_6",&SSDfracModmuz300_6);
1040 ntits->SetBranchAddress("errSSDfracModmuz300_6",&errSSDfracModmuz300_6);
1041 ntits->SetBranchAddress("SSDfracModsigz100_6",&SSDfracModsigz100_6);
1042 ntits->SetBranchAddress("errSSDfracModsigz100_6",&errSSDfracModsigz100_6);
1043 ntits->SetBranchAddress("SSDfracModsigz300_6",&SSDfracModsigz300_6);
1044 ntits->SetBranchAddress("errSSDfracModsigz300_6",&errSSDfracModsigz300_6);
1045 ntits->SetBranchAddress("SSDfracModsigz500_6",&SSDfracModsigz500_6);
1046 ntits->SetBranchAddress("errSSDfracModsigz500_6",&errSSDfracModsigz500_6);
1047
1048 // Sort entries according to run number in the chosen range
1049 // Same order is assumed for all the subsequent ntuples
1050 Int_t nr=ntits->GetEntries();
1051 Int_t *myIndex = new Int_t [nr];
1052 Int_t *noRuns = new Int_t [nr];
1053 Int_t kRunsToPlot=0;
1054 printf("Processing runs from %d up to %d\n",run1,run2);
1055 for(Int_t i=0; i<nr;i++){
1056 ntits->GetEvent(i);
1057 Int_t intrun = static_cast<Int_t>(nrun+0.01);
1058 if(intrun>=run1 && intrun<=run2){
1059 printf("Accepting run number %d in position %d\n",intrun,kRunsToPlot);
1060 noRuns[i]=intrun;
1061 kRunsToPlot++;
1062 }
1063 else {
1064 noRuns[i]=run2+10;
1065 printf("Rejecting run number %d - out of range\n",intrun);
1066 }
1067 }
1068 TMath::Sort(nr,noRuns,myIndex,kFALSE);
1069 printf("Total number of runs accepted for display %d\n",kRunsToPlot);
1070 if(kRunsToPlot==0)return;
1071 for(Int_t i=0;i<kRunsToPlot;i++)printf("Position %d ) Run: %d\n",i,noRuns[myIndex[i]]);
1072
1073 TH1F* hfmux20_3=new TH1F("hfmux20_3","",kRunsToPlot,0.,kRunsToPlot);
1074 TH1F* hfmux50_3=new TH1F("hfmux50_3","",kRunsToPlot,0.,kRunsToPlot);
1075 TH1F* hfmux100_3=new TH1F("hfmux100_3","",kRunsToPlot,0.,kRunsToPlot);
1076 TH1F* hfsigx100_3=new TH1F("hfsigx100_3","",kRunsToPlot,0.,kRunsToPlot);
1077 TH1F* hfsigx200_3=new TH1F("hfsigx200_3","",kRunsToPlot,0.,kRunsToPlot);
1078 TH1F* hfsigx300_3=new TH1F("hfsigx300_3","",kRunsToPlot,0.,kRunsToPlot);
1079 TH1F* hfmuxlr20_3=new TH1F("hfmuxlr20_3","",kRunsToPlot,0.,kRunsToPlot);
1080 TH1F* hfmuxlr50_3=new TH1F("hfmuxlr50_3","",kRunsToPlot,0.,kRunsToPlot);
1081 TH1F* hfmuxlr100_3=new TH1F("hfmuxlr100_3","",kRunsToPlot,0.,kRunsToPlot);
1082 TH1F* hfexcx100_3=new TH1F("hfexcx100_3","",kRunsToPlot,0.,kRunsToPlot);
1083 TH1F* hfexcx200_3=new TH1F("hfexcx200_3","",kRunsToPlot,0.,kRunsToPlot);
1084 TH1F* hfexcx300_3=new TH1F("hfexcx300_3","",kRunsToPlot,0.,kRunsToPlot);
1085 TH1F* hfmuz50_3=new TH1F("hfmuz50_3","",kRunsToPlot,0.,kRunsToPlot);
1086 TH1F* hfmuz100_3=new TH1F("hfmuz100_3","",kRunsToPlot,0.,kRunsToPlot);
1087 TH1F* hfmuz300_3=new TH1F("hfmuz300_3","",kRunsToPlot,0.,kRunsToPlot);
1088 TH1F* hfsigz100_3=new TH1F("hfsigz100_3","",kRunsToPlot,0.,kRunsToPlot);
1089 TH1F* hfsigz300_3=new TH1F("hfsigz300_3","",kRunsToPlot,0.,kRunsToPlot);
1090 TH1F* hfsigz500_3=new TH1F("hfsigz500_3","",kRunsToPlot,0.,kRunsToPlot);
1091
1092 TH1F* hfmux20_4=new TH1F("hfmux20_4","",kRunsToPlot,0.,kRunsToPlot);
1093 TH1F* hfmux50_4=new TH1F("hfmux50_4","",kRunsToPlot,0.,kRunsToPlot);
1094 TH1F* hfmux100_4=new TH1F("hfmux100_4","",kRunsToPlot,0.,kRunsToPlot);
1095 TH1F* hfsigx100_4=new TH1F("hfsigx100_4","",kRunsToPlot,0.,kRunsToPlot);
1096 TH1F* hfsigx200_4=new TH1F("hfsigx200_4","",kRunsToPlot,0.,kRunsToPlot);
1097 TH1F* hfsigx300_4=new TH1F("hfsigx300_4","",kRunsToPlot,0.,kRunsToPlot);
1098 TH1F* hfmuxlr20_4=new TH1F("hfmuxlr20_4","",kRunsToPlot,0.,kRunsToPlot);
1099 TH1F* hfmuxlr50_4=new TH1F("hfmuxlr50_4","",kRunsToPlot,0.,kRunsToPlot);
1100 TH1F* hfmuxlr100_4=new TH1F("hfmuxlr100_4","",kRunsToPlot,0.,kRunsToPlot);
1101 TH1F* hfexcx100_4=new TH1F("hfexcx100_4","",kRunsToPlot,0.,kRunsToPlot);
1102 TH1F* hfexcx200_4=new TH1F("hfexcx200_4","",kRunsToPlot,0.,kRunsToPlot);
1103 TH1F* hfexcx300_4=new TH1F("hfexcx300_4","",kRunsToPlot,0.,kRunsToPlot);
1104 TH1F* hfmuz50_4=new TH1F("hfmuz50_4","",kRunsToPlot,0.,kRunsToPlot);
1105 TH1F* hfmuz100_4=new TH1F("hfmuz100_4","",kRunsToPlot,0.,kRunsToPlot);
1106 TH1F* hfmuz300_4=new TH1F("hfmuz300_4","",kRunsToPlot,0.,kRunsToPlot);
1107 TH1F* hfsigz100_4=new TH1F("hfsigz100_4","",kRunsToPlot,0.,kRunsToPlot);
1108 TH1F* hfsigz300_4=new TH1F("hfsigz300_4","",kRunsToPlot,0.,kRunsToPlot);
1109 TH1F* hfsigz500_4=new TH1F("hfsigz500_4","",kRunsToPlot,0.,kRunsToPlot);
1110
1111 TH1F* hfmux20_1=new TH1F("hfmux20_1","",kRunsToPlot,0.,kRunsToPlot);
1112 TH1F* hfmux50_1=new TH1F("hfmux50_1","",kRunsToPlot,0.,kRunsToPlot);
1113 TH1F* hfmux100_1=new TH1F("hfmux100_1","",kRunsToPlot,0.,kRunsToPlot);
1114 TH1F* hfsigx100_1=new TH1F("hfsigx100_1","",kRunsToPlot,0.,kRunsToPlot);
1115 TH1F* hfsigx200_1=new TH1F("hfsigx200_1","",kRunsToPlot,0.,kRunsToPlot);
1116 TH1F* hfsigx300_1=new TH1F("hfsigx300_1","",kRunsToPlot,0.,kRunsToPlot);
1117 TH1F* hfmuz50_1=new TH1F("hfmuz50_1","",kRunsToPlot,0.,kRunsToPlot);
1118 TH1F* hfmuz100_1=new TH1F("hfmuz100_1","",kRunsToPlot,0.,kRunsToPlot);
1119 TH1F* hfmuz300_1=new TH1F("hfmuz300_1","",kRunsToPlot,0.,kRunsToPlot);
1120 TH1F* hfsigz100_1=new TH1F("hfsigz100_1","",kRunsToPlot,0.,kRunsToPlot);
1121 TH1F* hfsigz300_1=new TH1F("hfsigz300_1","",kRunsToPlot,0.,kRunsToPlot);
1122 TH1F* hfsigz500_1=new TH1F("hfsigz500_1","",kRunsToPlot,0.,kRunsToPlot);
1123
1124 TH1F* hfmux20_2=new TH1F("hfmux20_2","",kRunsToPlot,0.,kRunsToPlot);
1125 TH1F* hfmux50_2=new TH1F("hfmux50_2","",kRunsToPlot,0.,kRunsToPlot);
1126 TH1F* hfmux100_2=new TH1F("hfmux100_2","",kRunsToPlot,0.,kRunsToPlot);
1127 TH1F* hfsigx100_2=new TH1F("hfsigx100_2","",kRunsToPlot,0.,kRunsToPlot);
1128 TH1F* hfsigx200_2=new TH1F("hfsigx200_2","",kRunsToPlot,0.,kRunsToPlot);
1129 TH1F* hfsigx300_2=new TH1F("hfsigx300_2","",kRunsToPlot,0.,kRunsToPlot);
1130 TH1F* hfmuz50_2=new TH1F("hfmuz50_2","",kRunsToPlot,0.,kRunsToPlot);
1131 TH1F* hfmuz100_2=new TH1F("hfmuz100_2","",kRunsToPlot,0.,kRunsToPlot);
1132 TH1F* hfmuz300_2=new TH1F("hfmuz300_2","",kRunsToPlot,0.,kRunsToPlot);
1133 TH1F* hfsigz100_2=new TH1F("hfsigz100_2","",kRunsToPlot,0.,kRunsToPlot);
1134 TH1F* hfsigz300_2=new TH1F("hfsigz300_2","",kRunsToPlot,0.,kRunsToPlot);
1135 TH1F* hfsigz500_2=new TH1F("hfsigz500_2","",kRunsToPlot,0.,kRunsToPlot);
1136
1137 TH1F* hfmux20_5=new TH1F("hfmux20_5","",kRunsToPlot,0.,kRunsToPlot);
1138 TH1F* hfmux50_5=new TH1F("hfmux50_5","",kRunsToPlot,0.,kRunsToPlot);
1139 TH1F* hfmux100_5=new TH1F("hfmux100_5","",kRunsToPlot,0.,kRunsToPlot);
1140 TH1F* hfsigx100_5=new TH1F("hfsigx100_5","",kRunsToPlot,0.,kRunsToPlot);
1141 TH1F* hfsigx200_5=new TH1F("hfsigx200_5","",kRunsToPlot,0.,kRunsToPlot);
1142 TH1F* hfsigx300_5=new TH1F("hfsigx300_5","",kRunsToPlot,0.,kRunsToPlot);
1143 TH1F* hfmuz50_5=new TH1F("hfmuz50_5","",kRunsToPlot,0.,kRunsToPlot);
1144 TH1F* hfmuz100_5=new TH1F("hfmuz100_5","",kRunsToPlot,0.,kRunsToPlot);
1145 TH1F* hfmuz300_5=new TH1F("hfmuz300_5","",kRunsToPlot,0.,kRunsToPlot);
1146 TH1F* hfsigz100_5=new TH1F("hfsigz100_5","",kRunsToPlot,0.,kRunsToPlot);
1147 TH1F* hfsigz300_5=new TH1F("hfsigz300_5","",kRunsToPlot,0.,kRunsToPlot);
1148 TH1F* hfsigz500_5=new TH1F("hfsigz500_5","",kRunsToPlot,0.,kRunsToPlot);
1149
1150 TH1F* hfmux20_6=new TH1F("hfmux20_6","",kRunsToPlot,0.,kRunsToPlot);
1151 TH1F* hfmux50_6=new TH1F("hfmux50_6","",kRunsToPlot,0.,kRunsToPlot);
1152 TH1F* hfmux100_6=new TH1F("hfmux100_6","",kRunsToPlot,0.,kRunsToPlot);
1153 TH1F* hfsigx100_6=new TH1F("hfsigx100_6","",kRunsToPlot,0.,kRunsToPlot);
1154 TH1F* hfsigx200_6=new TH1F("hfsigx200_6","",kRunsToPlot,0.,kRunsToPlot);
1155 TH1F* hfsigx300_6=new TH1F("hfsigx300_6","",kRunsToPlot,0.,kRunsToPlot);
1156 TH1F* hfmuz50_6=new TH1F("hfmuz50_6","",kRunsToPlot,0.,kRunsToPlot);
1157 TH1F* hfmuz100_6=new TH1F("hfmuz100_6","",kRunsToPlot,0.,kRunsToPlot);
1158 TH1F* hfmuz300_6=new TH1F("hfmuz300_6","",kRunsToPlot,0.,kRunsToPlot);
1159 TH1F* hfsigz100_6=new TH1F("hfsigz100_6","",kRunsToPlot,0.,kRunsToPlot);
1160 TH1F* hfsigz300_6=new TH1F("hfsigz300_6","",kRunsToPlot,0.,kRunsToPlot);
1161 TH1F* hfsigz500_6=new TH1F("hfsigz500_6","",kRunsToPlot,0.,kRunsToPlot);
1162
1163
1164 for(Int_t i=0; i<kRunsToPlot;i++){
1165 ntits->GetEvent(myIndex[i]);
1166
1167 hfmux20_3->SetBinContent(i+1,SDDfracModmux20_3);
1168// hfmux20_3->SetBinError(i+1,errSDDfracModmux20_3);
1169 hfmux20_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1170 hfmux50_3->SetBinContent(i+1,SDDfracModmux50_3);
1171// hfmux50_3->SetBinError(i+1,errSDDfracModmux50_3);
1172 hfmux50_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1173 hfmux100_3->SetBinContent(i+1,SDDfracModmux100_3);
1174// hfmux100_3->SetBinError(i+1,errSDDfracModmux100_3);
1175 hfmux100_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1176 hfsigx100_3->SetBinContent(i+1,SDDfracModsigx100_3);
1177// hfsigx100_3->SetBinError(i+1,errSDDfracModsigx100_3);
1178 hfsigx100_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1179 hfsigx200_3->SetBinContent(i+1,SDDfracModsigx200_3);
1180// hfsigx200_3->SetBinError(i+1,errSDDfracModsigx200_3);
1181 hfsigx200_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1182 hfsigx300_3->SetBinContent(i+1,SDDfracModsigx300_3);
1183// hfsigx300_3->SetBinError(i+1,errSDDfracModsigx300_3);
1184 hfsigx300_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1185 hfmuxlr20_3->SetBinContent(i+1,SDDfracModmuxlr20_3);
1186// hfmuxlr20_3->SetBinError(i+1,errSDDfracModmuxlr20_3);
1187 hfmuxlr20_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1188 hfmuxlr50_3->SetBinContent(i+1,SDDfracModmuxlr50_3);
1189// hfmuxlr50_3->SetBinError(i+1,errSDDfracModmux50_3);
1190 hfmuxlr50_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1191 hfmuxlr100_3->SetBinContent(i+1,SDDfracModmux100_3);
1192// hfmuxlr100_3->SetBinError(i+1,errSDDfracModmux100_3);
1193 hfmuxlr100_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1194 hfexcx100_3->SetBinContent(i+1,SDDfracModexcx100_3);
1195// hfexcx100_3->SetBinError(i+1,errSDDfracModexcx100_3);
1196 hfexcx100_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1197 hfexcx200_3->SetBinContent(i+1,SDDfracModexcx200_3);
1198// hfexcx200_3->SetBinError(i+1,errSDDfracModexcx200_3);
1199 hfexcx200_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1200 hfexcx300_3->SetBinContent(i+1,SDDfracModexcx300_3);
1201// hfexcx300_3->SetBinError(i+1,errSDDfracModexcx300_3);
1202 hfexcx300_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1203 hfmuz50_3->SetBinContent(i+1,SDDfracModmuz50_3);
1204// hfmuz50_3->SetBinError(i+1,errSDDfracModmuz50_3);
1205 hfmuz50_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1206 hfmuz100_3->SetBinContent(i+1,SDDfracModmuz100_3);
1207// hfmuz100_3->SetBinError(i+1,errSDDfracModmuz100_3);
1208 hfmuz100_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1209 hfmuz300_3->SetBinContent(i+1,SDDfracModmuz300_3);
1210// hfmuz300_3->SetBinError(i+1,errSDDfracModmuz300_3);
1211 hfmuz300_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1212 hfsigz100_3->SetBinContent(i+1,SDDfracModsigz100_3);
1213// hfsigz100_3->SetBinError(i+1,errSDDfracModsigz100_3);
1214 hfsigz100_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1215 hfsigz300_3->SetBinContent(i+1,SDDfracModsigz300_3);
1216// hfsigz300_3->SetBinError(i+1,errSDDfracModsigz300_3);
1217 hfsigz300_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1218 hfsigz500_3->SetBinContent(i+1,SDDfracModsigz500_3);
1219// hfsigz500_3->SetBinError(i+1,errSDDfracModsigz500_3);
1220 hfsigz500_3->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1221
1222 hfmux20_4->SetBinContent(i+1,SDDfracModmux20_4);
1223// hfmux20_4->SetBinError(i+1,errSDDfracModmux20_4);
1224 hfmux20_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1225 hfmux50_4->SetBinContent(i+1,SDDfracModmux50_4);
1226// hfmux50_4->SetBinError(i+1,errSDDfracModmux50_4);
1227 hfmux50_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1228 hfmux100_4->SetBinContent(i+1,SDDfracModmux100_4);
1229// hfmux100_4->SetBinError(i+1,errSDDfracModmux100_4);
1230 hfmux100_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1231 hfsigx100_4->SetBinContent(i+1,SDDfracModsigx100_4);
1232// hfsigx100_4->SetBinError(i+1,errSDDfracModsigx100_4);
1233 hfsigx100_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1234 hfsigx200_4->SetBinContent(i+1,SDDfracModsigx200_4);
1235// hfsigx200_4->SetBinError(i+1,errSDDfracModsigx200_4);
1236 hfsigx200_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1237 hfsigx300_4->SetBinContent(i+1,SDDfracModsigx300_4);
1238// hfsigx300_4->SetBinError(i+1,errSDDfracModsigx300_4);
1239 hfsigx300_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1240 hfmuxlr20_4->SetBinContent(i+1,SDDfracModmuxlr20_4);
1241// hfmuxlr20_4->SetBinError(i+1,errSDDfracModmuxlr20_4);
1242 hfmuxlr20_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1243 hfmuxlr50_4->SetBinContent(i+1,SDDfracModmuxlr50_4);
1244// hfmuxlr50_4->SetBinError(i+1,errSDDfracModmux50_4);
1245 hfmuxlr50_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1246 hfmuxlr100_4->SetBinContent(i+1,SDDfracModmux100_4);
1247// hfmuxlr100_4->SetBinError(i+1,errSDDfracModmux100_4);
1248 hfmuxlr100_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1249 hfexcx100_4->SetBinContent(i+1,SDDfracModexcx100_4);
1250// hfexcx100_4->SetBinError(i+1,errSDDfracModexcx100_4);
1251 hfexcx100_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1252 hfexcx200_4->SetBinContent(i+1,SDDfracModexcx200_4);
1253// hfexcx200_4->SetBinError(i+1,errSDDfracModexcx200_4);
1254 hfexcx200_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1255 hfexcx300_4->SetBinContent(i+1,SDDfracModexcx300_4);
1256// hfexcx300_4->SetBinError(i+1,errSDDfracModexcx300_4);
1257 hfexcx300_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1258 hfmuz50_4->SetBinContent(i+1,SDDfracModmuz50_4);
1259// hfmuz50_4->SetBinError(i+1,errSDDfracModmuz50_4);
1260 hfmuz50_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1261 hfmuz100_4->SetBinContent(i+1,SDDfracModmuz100_4);
1262// hfmuz100_4->SetBinError(i+1,errSDDfracModmuz100_4);
1263 hfmuz100_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1264 hfmuz300_4->SetBinContent(i+1,SDDfracModmuz300_4);
1265// hfmuz300_4->SetBinError(i+1,errSDDfracModmuz300_4);
1266 hfmuz300_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1267 hfsigz100_4->SetBinContent(i+1,SDDfracModsigz100_4);
1268// hfsigz100_4->SetBinError(i+1,errSDDfracModsigz100_4);
1269 hfsigz100_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1270 hfsigz300_4->SetBinContent(i+1,SDDfracModsigz300_4);
1271// hfsigz300_4->SetBinError(i+1,errSDDfracModsigz300_4);
1272 hfsigz300_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1273 hfsigz500_4->SetBinContent(i+1,SDDfracModsigz500_4);
1274// hfsigz500_4->SetBinError(i+1,errSDDfracModsigz500_4);
1275 hfsigz500_4->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1276
1277 hfmux20_1->SetBinContent(i+1,SPDfracModmux20_1);
1278// hfmux20_1->SetBinError(i+1,errSPDfracModmux20_1);
1279 hfmux20_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1280 hfmux50_1->SetBinContent(i+1,SPDfracModmux50_1);
1281// hfmux50_1->SetBinError(i+1,errSPDfracModmux50_1);
1282 hfmux50_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1283 hfmux100_1->SetBinContent(i+1,SPDfracModmux100_1);
1284// hfmux100_1->SetBinError(i+1,errSPDfracModmux100_1);
1285 hfmux100_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1286 hfsigx100_1->SetBinContent(i+1,SPDfracModsigx100_1);
1287// hfsigx100_1->SetBinError(i+1,errSPDfracModsigx100_1);
1288 hfsigx100_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1289 hfsigx200_1->SetBinContent(i+1,SPDfracModsigx200_1);
1290// hfsigx200_1->SetBinError(i+1,errSPDfracModsigx200_1);
1291 hfsigx200_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1292 hfsigx300_1->SetBinContent(i+1,SPDfracModsigx300_1);
1293// hfsigx300_1->SetBinError(i+1,errSPDfracModsigx300_1);
1294 hfsigx300_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1295 hfmuz50_1->SetBinContent(i+1,SPDfracModmuz50_1);
1296// hfmuz50_1->SetBinError(i+1,errSPDfracModmuz50_1);
1297 hfmuz50_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1298 hfmuz100_1->SetBinContent(i+1,SPDfracModmuz100_1);
1299// hfmuz100_1->SetBinError(i+1,errSPDfracModmuz100_1);
1300 hfmuz100_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1301 hfmuz300_1->SetBinContent(i+1,SPDfracModmuz300_1);
1302// hfmuz300_1->SetBinError(i+1,errSPDfracModmuz300_1);
1303 hfmuz300_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1304 hfsigz100_1->SetBinContent(i+1,SPDfracModsigz100_1);
1305// hfsigz100_1->SetBinError(i+1,errSPDfracModsigz100_1);
1306 hfsigz100_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1307 hfsigz300_1->SetBinContent(i+1,SPDfracModsigz300_1);
1308// hfsigz300_1->SetBinError(i+1,errSPDfracModsigz300_1);
1309 hfsigz300_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1310 hfsigz500_1->SetBinContent(i+1,SPDfracModsigz500_1);
1311// hfsigz500_1->SetBinError(i+1,errSPDfracModsigz500_1);
1312 hfsigz500_1->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1313
1314 hfmux20_2->SetBinContent(i+1,SPDfracModmux20_2);
1315// hfmux20_2->SetBinError(i+1,errSPDfracModmux20_2);
1316 hfmux20_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1317 hfmux50_2->SetBinContent(i+1,SPDfracModmux50_2);
1318// hfmux50_2->SetBinError(i+1,errSPDfracModmux50_2);
1319 hfmux50_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1320 hfmux100_2->SetBinContent(i+1,SPDfracModmux100_2);
1321// hfmux100_2->SetBinError(i+1,errSPDfracModmux100_2);
1322 hfmux100_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1323 hfsigx100_2->SetBinContent(i+1,SPDfracModsigx100_2);
1324// hfsigx100_2->SetBinError(i+1,errSPDfracModsigx100_2);
1325 hfsigx100_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1326 hfsigx200_2->SetBinContent(i+1,SPDfracModsigx200_2);
1327// hfsigx200_2->SetBinError(i+1,errSPDfracModsigx200_2);
1328 hfsigx200_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1329 hfsigx300_2->SetBinContent(i+1,SPDfracModsigx300_2);
1330// hfsigx300_2->SetBinError(i+1,errSPDfracModsigx300_2);
1331 hfsigx300_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1332 hfmuz50_2->SetBinContent(i+1,SPDfracModmuz50_2);
1333// hfmuz50_2->SetBinError(i+1,errSPDfracModmuz50_2);
1334 hfmuz50_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1335 hfmuz100_2->SetBinContent(i+1,SPDfracModmuz100_2);
1336// hfmuz100_2->SetBinError(i+1,errSPDfracModmuz100_2);
1337 hfmuz100_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1338 hfmuz300_2->SetBinContent(i+1,SPDfracModmuz300_2);
1339// hfmuz300_2->SetBinError(i+1,errSPDfracModmuz300_2);
1340 hfmuz300_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1341 hfsigz100_2->SetBinContent(i+1,SPDfracModsigz100_2);
1342// hfsigz100_2->SetBinError(i+1,errSPDfracModsigz100_2);
1343 hfsigz100_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1344 hfsigz300_2->SetBinContent(i+1,SPDfracModsigz300_2);
1345// hfsigz300_2->SetBinError(i+1,errSPDfracModsigz300_2);
1346 hfsigz300_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1347 hfsigz500_2->SetBinContent(i+1,SPDfracModsigz500_2);
1348// hfsigz500_2->SetBinError(i+1,errSPDfracModsigz500_2);
1349 hfsigz500_2->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1350
1351 hfmux20_5->SetBinContent(i+1,SSDfracModmux20_5);
1352// hfmux20_5->SetBinError(i+1,errSSDfracModmux20_5);
1353 hfmux20_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1354 hfmux50_5->SetBinContent(i+1,SSDfracModmux50_5);
1355// hfmux50_5->SetBinError(i+1,errSSDfracModmux50_5);
1356 hfmux50_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1357 hfmux100_5->SetBinContent(i+1,SSDfracModmux100_5);
1358// hfmux100_5->SetBinError(i+1,errSSDfracModmux100_5);
1359 hfmux100_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1360 hfsigx100_5->SetBinContent(i+1,SSDfracModsigx100_5);
1361// hfsigx100_5->SetBinError(i+1,errSSDfracModsigx100_5);
1362 hfsigx100_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1363 hfsigx200_5->SetBinContent(i+1,SSDfracModsigx200_5);
1364// hfsigx200_5->SetBinError(i+1,errSSDfracModsigx200_5);
1365 hfsigx200_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1366 hfsigx300_5->SetBinContent(i+1,SSDfracModsigx300_5);
1367// hfsigx300_5->SetBinError(i+1,errSSDfracModsigx300_5);
1368 hfsigx300_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1369 hfmuz50_5->SetBinContent(i+1,SSDfracModmuz50_5);
1370// hfmuz50_5->SetBinError(i+1,errSSDfracModmuz50_5);
1371 hfmuz50_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1372 hfmuz100_5->SetBinContent(i+1,SSDfracModmuz100_5);
1373// hfmuz100_5->SetBinError(i+1,errSSDfracModmuz100_5);
1374 hfmuz100_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1375 hfmuz300_5->SetBinContent(i+1,SSDfracModmuz300_5);
1376// hfmuz300_5->SetBinError(i+1,errSSDfracModmuz300_5);
1377 hfmuz300_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1378 hfsigz100_5->SetBinContent(i+1,SSDfracModsigz100_5);
1379// hfsigz100_5->SetBinError(i+1,errSSDfracModsigz100_5);
1380 hfsigz100_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1381 hfsigz300_5->SetBinContent(i+1,SSDfracModsigz300_5);
1382// hfsigz300_5->SetBinError(i+1,errSSDfracModsigz300_5);
1383 hfsigz300_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1384 hfsigz500_5->SetBinContent(i+1,SSDfracModsigz500_5);
1385// hfsigz500_5->SetBinError(i+1,errSSDfracModsigz500_5);
1386 hfsigz500_5->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1387
1388 hfmux20_6->SetBinContent(i+1,SSDfracModmux20_6);
1389// hfmux20_6->SetBinError(i+1,errSSDfracModmux20_6);
1390 hfmux20_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1391 hfmux50_6->SetBinContent(i+1,SSDfracModmux50_6);
1392// hfmux50_6->SetBinError(i+1,errSSDfracModmux50_6);
1393 hfmux50_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1394 hfmux100_6->SetBinContent(i+1,SSDfracModmux100_6);
1395// hfmux100_6->SetBinError(i+1,errSSDfracModmux100_6);
1396 hfmux100_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1397 hfsigx100_6->SetBinContent(i+1,SSDfracModsigx100_6);
1398// hfsigx100_6->SetBinError(i+1,errSSDfracModsigx100_6);
1399 hfsigx100_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1400 hfsigx200_6->SetBinContent(i+1,SSDfracModsigx200_6);
1401// hfsigx200_6->SetBinError(i+1,errSSDfracModsigx200_6);
1402 hfsigx200_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1403 hfsigx300_6->SetBinContent(i+1,SSDfracModsigx300_6);
1404// hfsigx300_6->SetBinError(i+1,errSSDfracModsigx300_6);
1405 hfsigx300_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1406 hfmuz50_6->SetBinContent(i+1,SSDfracModmuz50_6);
1407// hfmuz50_6->SetBinError(i+1,errSSDfracModmuz50_6);
1408 hfmuz50_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1409 hfmuz100_6->SetBinContent(i+1,SSDfracModmuz100_6);
1410// hfmuz100_6->SetBinError(i+1,errSSDfracModmuz100_6);
1411 hfmuz100_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1412 hfmuz300_6->SetBinContent(i+1,SSDfracModmuz300_6);
1413// hfmuz300_6->SetBinError(i+1,errSSDfracModmuz300_6);
1414 hfmuz300_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1415 hfsigz100_6->SetBinContent(i+1,SSDfracModsigz100_6);
1416// hfsigz100_6->SetBinError(i+1,errSSDfracModsigz100_6);
1417 hfsigz100_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1418 hfsigz300_6->SetBinContent(i+1,SSDfracModsigz300_6);
1419// hfsigz300_6->SetBinError(i+1,errSSDfracModsigz300_6);
1420 hfsigz300_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1421 hfsigz500_6->SetBinContent(i+1,SSDfracModsigz500_6);
1422// hfsigz500_6->SetBinError(i+1,errSSDfracModsigz500_6);
1423 hfsigz500_6->GetXaxis()->SetBinLabel(i+1,Form("%d",(Int_t)nrun));
1424
1425
1426 }
1427
1428 gROOT->SetStyle("Plain");
1429 // gStyle->SetOptTitle(0);
1430 gStyle->SetOptStat(0);
1431 gStyle->SetFillColor(0);
1432 gStyle->SetTextFont(32);
1433
1434
1435
1436 TCanvas* c1=new TCanvas("c1","SDD Layer 3", 1200,800);
1437 c1->Divide(1,2);
1438 c1->cd(1);
1439 hfmux20_3->SetTitle("fraction(| #mu x-residual | < threshold) - layer 3");
1440 hfmux20_3->SetLineColor(2);
1441 hfmux20_3->SetMarkerStyle(31);
1442 hfmux20_3->SetMarkerSize(2);
1443 hfmux20_3->SetMarkerColor(2);
1444 hfmux20_3->SetMinimum(0.);
1445 hfmux20_3->SetMaximum(1.2);
1446 hfmux20_3->GetXaxis()->SetTitle("run number");
1447 hfmux20_3->GetYaxis()->SetTitleOffset(1.2);
1448 // hfmux20_3->GetYaxis()->SetTitle("%");
1449 hfmux20_3->Draw("P");
1450 hfmux50_3->SetLineColor(4);
1451 hfmux50_3->SetMarkerStyle(21);
1452 hfmux50_3->SetMarkerColor(4);
1453 hfmux50_3->Draw("same P");
1454 hfmux100_3->SetLineColor(6);
1455 hfmux100_3->SetMarkerStyle(22);
1456 hfmux100_3->SetMarkerColor(6);
1457 hfmux100_3->Draw("same P");
1458
1459 // quiqui
1460 TLine *l3 = new TLine(0.,0.877,kRunsToPlot,0.877);
1461 l3->SetLineStyle(2);
1462 l3->SetLineWidth(2);
1463// l3->Draw();
1464 TLegend* leg=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
1465 TLegendEntry* enta=leg->AddEntry(hfmux20_3,"|#mu res x | < 20 #mum","PL"); // P: marker, L; line
1466 enta->SetTextColor(hfmux20_3->GetMarkerColor());
1467 TLegendEntry* entb=leg->AddEntry(hfmux50_3,"|#mu res x | < 50 #mum","PL"); // P: marker, L; line
1468 entb->SetTextColor(hfmux50_3->GetMarkerColor());
1469 TLegendEntry* entc=leg->AddEntry(hfmux100_3,"|#mu res x | < 100 #mum","PL"); // P: marker, L; line
1470 entc->SetTextColor(hfmux100_3->GetMarkerColor());
1471 leg->SetFillStyle(0);
1472 leg->Draw();
1473 // c1->Update();
1474
1475 c1->cd(2);
1476 hfmuxlr20_3->SetTitle("fraction(| #muL - #muR | < threshold) - layer 3");
1477 hfmuxlr20_3->SetLineColor(2);
1478 hfmuxlr20_3->SetMarkerStyle(31);
1479 hfmuxlr20_3->SetMarkerSize(2);
1480 hfmuxlr20_3->SetMarkerColor(2);
1481 hfmuxlr20_3->SetMinimum(0.);
1482 hfmuxlr20_3->SetMaximum(1.2);
1483 hfmuxlr20_3->GetXaxis()->SetTitle("run number");
1484 hfmuxlr20_3->GetYaxis()->SetTitleOffset(1.2);
1485 // hfmuxlr20_3->GetYaxis()->SetTitle("%");
1486 hfmuxlr20_3->Draw("P");
1487 hfmuxlr50_3->SetLineColor(4);
1488 hfmuxlr50_3->SetMarkerStyle(21);
1489 hfmuxlr50_3->SetMarkerColor(4);
1490 hfmuxlr50_3->Draw("same P");
1491 hfmuxlr100_3->SetLineColor(6);
1492 hfmuxlr100_3->SetMarkerStyle(22);
1493 hfmuxlr100_3->SetMarkerColor(6);
1494 hfmuxlr100_3->Draw("same P");
1495
1496// l3->Draw();
1497
1498 TLegend* leg2=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
1499 TLegendEntry* ent2a=leg2->AddEntry(hfmuxlr20_3,"|#muL - #muR| < 20 #mum","PL"); // P: marker, L; line
1500 ent2a->SetTextColor(hfmuxlr20_3->GetMarkerColor());
1501 TLegendEntry* ent2b=leg2->AddEntry(hfmuxlr50_3,"|#muL - #muR| < 50 #mum","PL"); // P: marker, L; line
1502 ent2b->SetTextColor(hfmuxlr50_3->GetMarkerColor());
1503 TLegendEntry* ent2c=leg2->AddEntry(hfmuxlr100_3,"|#muL - #muR| < 100 #mum","PL"); // P: marker, L; line
1504 ent2c->SetTextColor(hfmuxlr100_3->GetMarkerColor());
1505 leg2->SetFillStyle(0);
1506 leg2->Draw();
1507
1508 c1->SaveAs("SDD_resx_L3");
1509// pdfFileNames+=" SDD_resx_L3";
1510 c1->Update();
1511
1512 TCanvas* c2=new TCanvas("c2","SDD Layer 4", 1200,800);
1513 c2->Divide(1,2);
1514 c2->cd(1);
1515 hfmux20_4->SetTitle("fraction(| #mu x-residual | < threshold) - layer 4");
1516 hfmux20_4->SetLineColor(2);
1517 hfmux20_4->SetMarkerStyle(31);
1518 hfmux20_4->SetMarkerSize(2);
1519 hfmux20_4->SetMarkerColor(2);
1520 hfmux20_4->SetMinimum(0.);
1521 hfmux20_4->SetMaximum(1.2);
1522 hfmux20_4->GetXaxis()->SetTitle("run number");
1523 hfmux20_4->GetYaxis()->SetTitleOffset(1.2);
1524 // hfmux20_4->GetYaxis()->SetTitle("%");
1525 hfmux20_4->Draw("P");
1526 hfmux50_4->SetLineColor(4);
1527 hfmux50_4->SetMarkerStyle(21);
1528 hfmux50_4->SetMarkerColor(4);
1529 hfmux50_4->Draw("same P");
1530 hfmux100_4->SetLineColor(6);
1531 hfmux100_4->SetMarkerStyle(22);
1532 hfmux100_4->SetMarkerColor(6);
1533 hfmux100_4->Draw("same P");
1534 TLine *l4 = new TLine(0.,0.839,kRunsToPlot,0.839);
1535 l4->SetLineStyle(2);
1536 l4->SetLineWidth(2);
1537// l4->Draw();
1538
1539 TLegend* legg=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
1540 TLegendEntry* entta=legg->AddEntry(hfmux20_4,"|#mu res x | < 20 #mum","PL"); // P: marker, L; line
1541 entta->SetTextColor(hfmux20_4->GetMarkerColor());
1542 TLegendEntry* enttb=legg->AddEntry(hfmux50_4,"|#mu res x | < 50 #mum","PL"); // P: marker, L; line
1543 enttb->SetTextColor(hfmux50_4->GetMarkerColor());
1544 TLegendEntry* enttc=legg->AddEntry(hfmux100_4,"|#mu res x | < 100 #mum","PL"); // P: marker, L; line
1545 enttc->SetTextColor(hfmux100_4->GetMarkerColor());
1546 legg->SetFillStyle(0);
1547 legg->Draw();
1548 // c2->Update();
1549
1550 c2->cd(2);
1551 hfmuxlr20_4->SetTitle("fraction(| #muL - #muR | < threshold) - layer 4");
1552 hfmuxlr20_4->SetLineColor(2);
1553 hfmuxlr20_4->SetMarkerStyle(31);
1554 hfmuxlr20_4->SetMarkerSize(2);
1555 hfmuxlr20_4->SetMarkerColor(2);
1556 hfmuxlr20_4->SetMinimum(0.);
1557 hfmuxlr20_4->SetMaximum(1.2);
1558 hfmuxlr20_4->GetXaxis()->SetTitle("run number");
1559 hfmuxlr20_4->GetYaxis()->SetTitleOffset(1.2);
1560 // hfmuxlr20_4->GetYaxis()->SetTitle("%");
1561 hfmuxlr20_4->Draw("P");
1562 hfmuxlr50_4->SetLineColor(4);
1563 hfmuxlr50_4->SetMarkerStyle(21);
1564 hfmuxlr50_4->SetMarkerColor(4);
1565 hfmuxlr50_4->Draw("same P");
1566 hfmuxlr100_4->SetLineColor(6);
1567 hfmuxlr100_4->SetMarkerStyle(22);
1568 hfmuxlr100_4->SetMarkerColor(6);
1569 hfmuxlr100_4->Draw("same P");
1570
1571// l4->Draw();
1572
1573 TLegend* legg2=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
1574 TLegendEntry* entt2a=legg2->AddEntry(hfmuxlr20_4,"|#muL - #muR| < 20 #mum","PL"); // P: marker, L; line
1575 entt2a->SetTextColor(hfmuxlr20_4->GetMarkerColor());
1576 TLegendEntry* entt2b=legg2->AddEntry(hfmuxlr50_4,"|#muL - #muR| < 50 #mum","PL"); // P: marker, L; line
1577 entt2b->SetTextColor(hfmuxlr50_4->GetMarkerColor());
1578 TLegendEntry* entt2c=legg2->AddEntry(hfmuxlr100_4,"|#muL - #muR| < 100 #mum","PL"); // P: marker, L; line
1579 entt2c->SetTextColor(hfmuxlr100_4->GetMarkerColor());
1580 legg2->SetFillStyle(0);
1581 legg2->Draw();
1582
1583 c2->SaveAs("SDD_resx_L4");
1584// pdfFileNames+=" SDD_resx_L4";
1585 c2->Update();
1586
1587 //////////qui
1588
1589 TCanvas* c3=new TCanvas("c3","SDD z residuals", 1200,800);
1590 c3->Divide(1,2);
1591 c3->cd(1);
1592 hfmuz50_3->SetTitle("fraction(| #mu z-residual | < threshold) - layer 3");
1593 hfmuz50_3->SetLineColor(2);
1594 hfmuz50_3->SetMarkerStyle(31);
1595 hfmuz50_3->SetMarkerSize(2);
1596 hfmuz50_3->SetMarkerColor(2);
1597 hfmuz50_3->SetMinimum(0.);
1598 hfmuz50_3->SetMaximum(1.2);
1599 hfmuz50_3->GetXaxis()->SetTitle("run number");
1600 hfmuz50_3->GetYaxis()->SetTitleOffset(1.2);
1601 // hfmuz50_3->GetYaxis()->SetTitle("%");
1602 hfmuz50_3->Draw("P");
1603 hfmuz100_3->SetLineColor(4);
1604 hfmuz100_3->SetMarkerStyle(21);
1605 hfmuz100_3->SetMarkerColor(4);
1606 hfmuz100_3->Draw("same P");
1607 hfmuz300_3->SetLineColor(6);
1608 hfmuz300_3->SetMarkerStyle(22);
1609 hfmuz300_3->SetMarkerColor(6);
1610 hfmuz300_3->Draw("same P");
1611// l3->Draw();
1612 TLegend* lege=new TLegend(0.5,0.15,0.88,0.35); // x1, y1, x2, y2
1613 TLegendEntry* ent_a=lege->AddEntry(hfmuz50_3,"|#mu res z| < 30 #mum","PL"); // P: marker, L; line
1614 ent_a->SetTextColor(hfmuz50_3->GetMarkerColor());
1615 TLegendEntry* ent_b=lege->AddEntry(hfmuz100_3,"|#mu res z| < 50 #mum","PL"); // P: marker, L; line
1616 ent_b->SetTextColor(hfmuz100_3->GetMarkerColor());
1617 TLegendEntry* ent_c=lege->AddEntry(hfmuz300_3,"|#mu res z| < 100 #mum","PL"); // P: marker, L; line
1618 ent_c->SetTextColor(hfmuz300_3->GetMarkerColor());
1619 lege->SetFillStyle(0);
1620 lege->Draw();
1621 // c3->Update();
1622
1623 c3->cd(2);
1624 hfmuz50_4->SetTitle("fraction(| #mu z-residual | < threshold) - layer 4");
1625 hfmuz50_4->SetLineColor(2);
1626 hfmuz50_4->SetMarkerStyle(31);
1627 hfmuz50_4->SetMarkerSize(2);
1628 hfmuz50_4->SetMarkerColor(2);
1629 hfmuz50_4->SetMinimum(0.);
1630 hfmuz50_4->SetMaximum(1.2);
1631 hfmuz50_4->GetXaxis()->SetTitle("run number");
1632 hfmuz50_4->GetYaxis()->SetTitleOffset(1.2);
1633 // hfmuz50_4->GetYaxis()->SetTitle("%");
1634 hfmuz50_4->Draw("P");
1635 hfmuz100_4->SetLineColor(4);
1636 hfmuz100_4->SetMarkerStyle(21);
1637 hfmuz100_4->SetMarkerColor(4);
1638 hfmuz100_4->Draw("same P");
1639 hfmuz300_4->SetLineColor(6);
1640 hfmuz300_4->SetMarkerStyle(22);
1641 hfmuz300_4->SetMarkerColor(6);
1642 hfmuz300_4->Draw("same P");
1643// l4->Draw();
1644 TLegend* leggi=new TLegend(0.5,0.15,0.88,0.35); // x1, y1, x2, y2
1645 TLegendEntry* en_ta=leggi->AddEntry(hfmuz50_4,"|#mu res z| < 30 #mum","PL"); // P: marker, L; line
1646 en_ta->SetTextColor(hfmuz50_4->GetMarkerColor());
1647 TLegendEntry* en_tb=leggi->AddEntry(hfmuz100_4,"|#mu res z| < 50 #mum","PL"); // P: marker, L; line
1648 en_tb->SetTextColor(hfmuz100_4->GetMarkerColor());
1649 TLegendEntry* en_tc=leggi->AddEntry(hfmuz300_4,"|#mu res z| < 100 #mum","PL"); // P: marker, L; line
1650 en_tc->SetTextColor(hfmuz300_4->GetMarkerColor());
1651 leggi->SetFillStyle(0);
1652 leggi->Draw();
1653 c3->SaveAs("SDD_resz");
1654// pdfFileNames+=" SDD_resz";
1655 c3->Update();
1656
1657 //////////////////////////////////////////////// SPD ////////////////////////////////////////////////////
1658
1659 TCanvas* c4=new TCanvas("c4","SPD x residuals", 1200,800);
1660 c4->Divide(1,2);
1661 c4->cd(1);
1662 hfmux20_1->SetTitle("fraction(| #mu x-residual | < threshold) - layer 1");
1663 hfmux20_1->SetLineColor(2);
1664 hfmux20_1->SetMarkerStyle(31);
1665 hfmux20_1->SetMarkerSize(2);
1666 hfmux20_1->SetMarkerColor(2);
1667 hfmux20_1->SetMinimum(0.);
1668 hfmux20_1->SetMaximum(1.2);
1669 hfmux20_1->GetXaxis()->SetTitle("run number");
1670 hfmux20_1->GetYaxis()->SetTitleOffset(1.2);
1671 // hfmux20_1->GetYaxis()->SetTitle("%");
1672 hfmux20_1->Draw("P");
1673 hfmux50_1->SetLineColor(4);
1674 hfmux50_1->SetMarkerStyle(21);
1675 hfmux50_1->SetMarkerColor(4);
1676 hfmux50_1->Draw("same P");
1677 hfmux100_1->SetLineColor(6);
1678 hfmux100_1->SetMarkerStyle(22);
1679 hfmux100_1->SetMarkerColor(6);
1680 hfmux100_1->Draw("same P");
1681 TLine *l1 = new TLine(0.,0.625,kRunsToPlot,0.625);
1682 l1->SetLineStyle(2);
1683 l1->SetLineWidth(2);
1684// l1->Draw();
1685
1686 TLegend* lege1=new TLegend(0.5,0.15,0.88,0.35); // x1, y1, x2, y2
1687 TLegendEntry* ent_a1=lege1->AddEntry(hfmux20_1,"|#mu res x| < 5 #mum","PL"); // P: marker, L; line
1688 ent_a1->SetTextColor(hfmux20_1->GetMarkerColor());
1689 TLegendEntry* ent_b1=lege1->AddEntry(hfmux50_1,"|#mu res x| < 10 #mum","PL"); // P: marker, L; line
1690 ent_b1->SetTextColor(hfmux50_1->GetMarkerColor());
1691 TLegendEntry* ent_c1=lege1->AddEntry(hfmux100_1,"|#mu res x| < 50 #mum","PL"); // P: marker, L; line
1692 ent_c1->SetTextColor(hfmux100_1->GetMarkerColor());
1693 lege1->SetFillStyle(0);
1694 lege1->Draw();
1695 // c4->Update();
1696
1697 c4->cd(2);
1698 hfmux20_2->SetTitle("fraction(| #mu x-residual | < threshold) - layer 2");
1699 hfmux20_2->SetLineColor(2);
1700 hfmux20_2->SetMarkerStyle(31);
1701 hfmux20_2->SetMarkerSize(2);
1702 hfmux20_2->SetMarkerColor(2);
1703 hfmux20_2->SetMinimum(0.);
1704 hfmux20_2->SetMaximum(1.2);
1705 hfmux20_2->GetXaxis()->SetTitle("run number");
1706 hfmux20_2->GetYaxis()->SetTitleOffset(1.2);
1707 // hfmux20_2->GetYaxis()->SetTitle("%");
1708 hfmux20_2->Draw("P");
1709 hfmux50_2->SetLineColor(4);
1710 hfmux50_2->SetMarkerStyle(21);
1711 hfmux50_2->SetMarkerColor(4);
1712 hfmux50_2->Draw("same P");
1713 hfmux100_2->SetLineColor(6);
1714 hfmux100_2->SetMarkerStyle(22);
1715 hfmux100_2->SetMarkerColor(6);
1716 hfmux100_2->Draw("same P");
1717 TLine *l2 = new TLine(0.,0.631,kRunsToPlot,0.631);
1718 l2->SetLineStyle(2);
1719 l2->SetLineWidth(2);
1720// l2->Draw();
1721
1722 TLegend* leggi1=new TLegend(0.5,0.15,0.88,0.35); // x1, y1, x2, y2
1723 TLegendEntry* en_ta1=leggi1->AddEntry(hfmux20_2,"|#mu res x| < 10 #mum","PL"); // P: marker, L; line
1724 en_ta1->SetTextColor(hfmux20_2->GetMarkerColor());
1725 TLegendEntry* en_tb1=leggi1->AddEntry(hfmux50_2,"|#mu res x| < 20 #mum","PL"); // P: marker, L; line
1726 en_tb1->SetTextColor(hfmux50_2->GetMarkerColor());
1727 TLegendEntry* en_tc1=leggi1->AddEntry(hfmux100_2,"|#mu res x| < 50 #mum","PL"); // P: marker, L; line
1728 en_tc1->SetTextColor(hfmux100_2->GetMarkerColor());
1729 leggi1->SetFillStyle(0);
1730 leggi1->Draw();
1731 c4->SaveAs("SPD_resx");
1732// pdfFileNames+=" SPD_resx";
1733 c4->Update();
1734
1735 TCanvas* c5=new TCanvas("c5","SPD z residuals", 1200,800);
1736 c5->Divide(1,2);
1737 c5->cd(1);
1738 hfmuz50_1->SetTitle("fraction(| #mu z-residual | < threshold) - layer 1");
1739 hfmuz50_1->SetLineColor(2);
1740 hfmuz50_1->SetMarkerStyle(31);
1741 hfmuz50_1->SetMarkerColor(2);
1742 hfmuz50_1->SetMinimum(0.);
1743 hfmuz50_1->SetMaximum(1.2);
1744 hfmuz50_1->GetXaxis()->SetTitle("run number");
1745 hfmuz50_1->GetYaxis()->SetTitleOffset(1.2);
1746 // hfmuz50_1->GetYaxis()->SetTitle("%");
1747 hfmuz50_1->Draw("P");
1748 hfmuz100_1->SetLineColor(4);
1749 hfmuz100_1->SetMarkerStyle(21);
1750 hfmuz100_1->SetMarkerColor(4);
1751 hfmuz100_1->Draw("same P");
1752 hfmuz300_1->SetLineColor(6);
1753 hfmuz300_1->SetMarkerStyle(22);
1754 hfmuz300_1->SetMarkerColor(6);
1755 hfmuz300_1->Draw("same P");
1756// l1->Draw();
1757
1758 TLegend* lege2=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
1759 TLegendEntry* ent_a2=lege2->AddEntry(hfmuz50_1,"|#mu res z| < 10 #mum","PL"); // P: marker, L; line
1760 ent_a2->SetTextColor(hfmuz50_1->GetMarkerColor());
1761 TLegendEntry* ent_b2=lege2->AddEntry(hfmuz100_1,"|#mu res z| < 20 #mum","PL"); // P: marker, L; line
1762 ent_b2->SetTextColor(hfmuz100_1->GetMarkerColor());
1763 TLegendEntry* ent_c2=lege2->AddEntry(hfmuz300_1,"|#mu res z| < 50 #mum","PL"); // P: marker, L; line
1764 ent_c2->SetTextColor(hfmuz300_1->GetMarkerColor());
1765 lege2->SetFillStyle(0);
1766 lege2->Draw();
1767 // c3->Update();
1768
1769 c5->cd(2);
1770 hfmuz50_2->SetTitle("fraction(| #mu z-residual | < threshold) - layer 2");
1771 hfmuz50_2->SetLineColor(2);
1772 hfmuz50_2->SetMarkerStyle(31);
1773 hfmuz50_2->SetMarkerSize(2);
1774 hfmuz50_2->SetMarkerColor(2);
1775 hfmuz50_2->SetMinimum(0.);
1776 hfmuz50_2->SetMaximum(1.2);
1777 hfmuz50_2->GetXaxis()->SetTitle("run number");
1778 hfmuz50_2->GetYaxis()->SetTitleOffset(1.2);
1779 // hfmuz50_2->GetYaxis()->SetTitle("%");
1780 hfmuz50_2->Draw("P");
1781 hfmuz100_2->SetLineColor(4);
1782 hfmuz100_2->SetMarkerStyle(21);
1783 hfmuz100_2->SetMarkerColor(4);
1784 hfmuz100_2->Draw("same P");
1785 hfmuz300_2->SetLineColor(6);
1786 hfmuz300_2->SetMarkerStyle(22);
1787 hfmuz300_2->SetMarkerColor(6);
1788 hfmuz300_2->Draw("same P");
1789// l2->Draw();
1790
1791 TLegend* leggi2=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
1792 TLegendEntry* en_ta2=leggi2->AddEntry(hfmuz50_2,"|#mu res z| < 10 #mum","PL"); // P: marker, L; line
1793 en_ta2->SetTextColor(hfmuz50_2->GetMarkerColor());
1794 TLegendEntry* en_tb2=leggi2->AddEntry(hfmuz100_2,"|#mu res z| < 20 #mum","PL"); // P: marker, L; line
1795 en_tb2->SetTextColor(hfmuz100_2->GetMarkerColor());
1796 TLegendEntry* en_tc2=leggi2->AddEntry(hfmuz300_2,"|#mu res z| < 50 #mum","PL"); // P: marker, L; line
1797 en_tc2->SetTextColor(hfmuz300_2->GetMarkerColor());
1798 leggi2->SetFillStyle(0);
1799 leggi2->Draw();
1800 c5->SaveAs("SPD_resz");
1801// pdfFileNames+=" SPD_resz";
1802 c5->Update();
1803
1804 //////////////////////////////////////////////// SSD ////////////////////////////////////////////////////
1805
1806 TCanvas* c6=new TCanvas("c6","SSD x residuals", 1200,800);
1807 c6->Divide(1,2);
1808 c6->cd(1);
1809 hfmux20_5->SetTitle("fraction(| #mu x-residual | < threshold) - layer 5");
1810 hfmux20_5->SetLineColor(2);
1811 hfmux20_5->SetMarkerStyle(31);
1812 hfmux20_5->SetMarkerSize(2);
1813 hfmux20_5->SetMarkerColor(2);
1814 hfmux20_5->SetMinimum(0.);
1815 hfmux20_5->SetMaximum(1.2);
1816 hfmux20_5->GetXaxis()->SetTitle("run number");
1817 hfmux20_5->GetYaxis()->SetTitleOffset(1.2);
1818 // hfmux20_5->GetYaxis()->SetTitle("%");
1819 hfmux20_5->Draw("P");
1820 hfmux50_5->SetLineColor(4);
1821 hfmux50_5->SetMarkerStyle(21);
1822 hfmux50_5->SetMarkerColor(4);
1823 hfmux50_5->Draw("same P");
1824 hfmux100_5->SetLineColor(6);
1825 hfmux100_5->SetMarkerStyle(22);
1826 hfmux100_5->SetMarkerColor(6);
1827 hfmux100_5->Draw("same P");
1828 TLine *l5 = new TLine(0.,0.915,kRunsToPlot,0.915);
1829 l5->SetLineStyle(2);
1830 l5->SetLineWidth(2);
1831// l5->Draw();
1832
1833 TLegend* lege3=new TLegend(0.5,0.15,0.88,0.35); // x1, y1, x2, y2
1834 TLegendEntry* ent_a3=lege3->AddEntry(hfmux20_5,"|#mu res x| < 20 #mum","PL"); // P: marker, L; line
1835 ent_a3->SetTextColor(hfmux20_5->GetMarkerColor());
1836 TLegendEntry* ent_b3=lege3->AddEntry(hfmux50_5,"|#mu res x| < 40 #mum","PL"); // P: marker, L; line
1837 ent_b3->SetTextColor(hfmux50_5->GetMarkerColor());
1838 TLegendEntry* ent_c3=lege3->AddEntry(hfmux100_5,"|#mu res x| < 100 #mum","PL"); // P: marker, L; line
1839 ent_c3->SetTextColor(hfmux100_5->GetMarkerColor());
1840 lege3->SetFillStyle(0);
1841 lege3->Draw();
1842 // c6->Update();
1843
1844 c6->cd(2);
1845 hfmux20_6->SetTitle("fraction(| #mu x-residual | < threshold) - layer 6");
1846 hfmux20_6->SetLineColor(2);
1847 hfmux20_6->SetMarkerStyle(31);
1848 hfmux20_6->SetMarkerSize(2);
1849 hfmux20_6->SetMarkerColor(2);
1850 hfmux20_6->SetMinimum(0.);
1851 hfmux20_6->SetMaximum(1.2);
1852 hfmux20_6->GetXaxis()->SetTitle("run number");
1853 hfmux20_6->GetYaxis()->SetTitleOffset(1.2);
1854 // hfmux20_6->GetYaxis()->SetTitle("%");
1855 hfmux20_6->Draw("P");
1856 hfmux50_6->SetLineColor(4);
1857 hfmux50_6->SetMarkerStyle(21);
1858 hfmux50_6->SetMarkerColor(4);
1859 hfmux50_6->Draw("same P");
1860 hfmux100_6->SetLineColor(6);
1861 hfmux100_6->SetMarkerStyle(22);
1862 hfmux100_6->SetMarkerColor(6);
1863 hfmux100_6->Draw("same P");
1864 TLine *l6 = new TLine(0.,0.881,kRunsToPlot,0.881);
1865 l6->SetLineStyle(2);
1866 l6->SetLineWidth(2);
1867// l6->Draw();
1868
1869 TLegend* leggi3=new TLegend(0.5,0.15,0.88,0.35); // x1, y1, x2, y2
1870 TLegendEntry* en_ta3=leggi3->AddEntry(hfmux20_6,"|#mu res x| < 20 #mum","PL"); // P: marker, L; line
1871 en_ta3->SetTextColor(hfmux20_6->GetMarkerColor());
1872 TLegendEntry* en_tb3=leggi3->AddEntry(hfmux50_6,"|#mu res x| < 40 #mum","PL"); // P: marker, L; line
1873 en_tb3->SetTextColor(hfmux50_6->GetMarkerColor());
1874 TLegendEntry* en_tc3=leggi3->AddEntry(hfmux100_6,"|#mu res x| < 100 #mum","PL"); // P: marker, L; line
1875 en_tc3->SetTextColor(hfmux100_6->GetMarkerColor());
1876 leggi3->SetFillStyle(0);
1877 leggi3->Draw();
1878 c6->SaveAs("SSD_resx");
1879// pdfFileNames+=" SSD_resx";
1880 c6->Update();
1881
1882 TCanvas* c7=new TCanvas("c7","SSD z residuals", 1200,800);
1883 c7->Divide(1,2);
1884 c7->cd(1);
1885 hfmuz50_5->SetTitle("fraction(| #mu z-residual | < threshold) - layer 5");
1886 hfmuz50_5->SetLineColor(2);
1887 hfmuz50_5->SetMarkerStyle(31);
1888 hfmuz50_5->SetMarkerSize(2);
1889 hfmuz50_5->SetMarkerColor(2);
1890 hfmuz50_5->SetMinimum(0.);
1891 hfmuz50_5->SetMaximum(1.2);
1892 hfmuz50_5->GetXaxis()->SetTitle("run number");
1893 hfmuz50_5->GetYaxis()->SetTitleOffset(1.2);
1894 // hfmuz50_5->GetYaxis()->SetTitle("%");
1895 hfmuz50_5->Draw("P");
1896 hfmuz100_5->SetLineColor(4);
1897 hfmuz100_5->SetMarkerStyle(21);
1898 hfmuz100_5->SetMarkerColor(4);
1899 hfmuz100_5->Draw("same P");
1900 hfmuz300_5->SetLineColor(6);
1901 hfmuz300_5->SetMarkerStyle(22);
1902 hfmuz300_5->SetMarkerColor(6);
1903 hfmuz300_5->Draw("same P");
1904// l5->Draw();
1905
1906 TLegend* lege4=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
1907 TLegendEntry* ent_a4=lege4->AddEntry(hfmuz50_5,"|#mu res z| < 50 #mum","PL"); // P: marker, L; line
1908 ent_a4->SetTextColor(hfmuz50_5->GetMarkerColor());
1909 TLegendEntry* ent_b4=lege4->AddEntry(hfmuz100_5,"|#mu res z| < 100 #mum","PL"); // P: marker, L; line
1910 ent_b4->SetTextColor(hfmuz100_5->GetMarkerColor());
1911 TLegendEntry* ent_c4=lege4->AddEntry(hfmuz300_5,"|#mu res z| < 300 #mum","PL"); // P: marker, L; line
1912 ent_c4->SetTextColor(hfmuz300_5->GetMarkerColor());
1913 lege4->SetFillStyle(0);
1914 lege4->Draw();
1915 // c7->Update();
1916
1917 c7->cd(2);
1918 hfmuz50_6->SetTitle("fraction(| #mu z-residual | < threshold) - layer 6");
1919 hfmuz50_6->SetLineColor(2);
1920 hfmuz50_6->SetMarkerStyle(31);
1921 hfmuz50_6->SetMarkerSize(2);
1922 hfmuz50_6->SetMarkerColor(2);
1923 hfmuz50_6->SetMinimum(0.);
1924 hfmuz50_6->SetMaximum(1.2);
1925 hfmuz50_6->GetXaxis()->SetTitle("run number");
1926 hfmuz50_6->GetYaxis()->SetTitleOffset(1.2);
1927 // hfmuz50_6->GetYaxis()->SetTitle("%");
1928 hfmuz50_6->Draw("P");
1929 hfmuz100_6->SetLineColor(4);
1930 hfmuz100_6->SetMarkerStyle(21);
1931 hfmuz100_6->SetMarkerColor(4);
1932 hfmuz100_6->Draw("same P");
1933 hfmuz300_6->SetLineColor(6);
1934 hfmuz300_6->SetMarkerStyle(22);
1935 hfmuz300_6->SetMarkerColor(6);
1936 hfmuz300_6->Draw("same P");
1937// l6->Draw();
1938 TLegend* leggi4=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
1939 TLegendEntry* en_ta4=leggi4->AddEntry(hfmuz50_6,"|#mu res z| < 50 #mum","PL"); // P: marker, L; line
1940 en_ta4->SetTextColor(hfmuz50_6->GetMarkerColor());
1941 TLegendEntry* en_tb4=leggi4->AddEntry(hfmuz100_6,"|#mu res z| < 100 #mum","PL"); // P: marker, L; line
1942 en_tb4->SetTextColor(hfmuz100_6->GetMarkerColor());
1943 TLegendEntry* en_tc4=leggi4->AddEntry(hfmuz300_6,"|#mu res z| < 300 #mum","PL"); // P: marker, L; line
1944 en_tc4->SetTextColor(hfmuz300_6->GetMarkerColor());
1945 leggi4->SetFillStyle(0);
1946 leggi4->Draw();
1947 c7->SaveAs("SSD_resz");
1948// pdfFileNames+=" SSD_resz";
1949 c7->Update();
1950
1951
1952 // -------------------------- plot sigma --------------------------------------
1953 // -------------------------- SPD x -------------------------------------
1954
1955 TCanvas* c8=new TCanvas("c8","SPD x res width", 1200,800);
1956 c8->Divide(1,2);
1957 c8->cd(1);
1958 hfsigx100_1->SetTitle("fraction(#sigma x-residual < threshold) - layer 1");
1959 hfsigx100_1->SetLineColor(2);
1960 hfsigx100_1->SetMarkerStyle(31);
1961 hfsigx100_1->SetMarkerSize(2);
1962 hfsigx100_1->SetMarkerColor(2);
1963 hfsigx100_1->SetMinimum(0.);
1964 hfsigx100_1->SetMaximum(1.2);
1965 hfsigx100_1->GetXaxis()->SetTitle("run number");
1966 hfsigx100_1->GetYaxis()->SetTitleOffset(1.2);
1967 // hfsigx100_1->GetYaxis()->SetTitle("%");
1968 hfsigx100_1->Draw("P");
1969 hfsigx200_1->SetLineColor(4);
1970 hfsigx200_1->SetMarkerStyle(21);
1971 hfsigx200_1->SetMarkerColor(4);
1972 hfsigx200_1->Draw("same P");
1973 hfsigx300_1->SetLineColor(6);
1974 hfsigx300_1->SetMarkerStyle(22);
1975 hfsigx300_1->SetMarkerColor(6);
1976 hfsigx300_1->Draw("same P");
1977// l1->Draw();
1978
1979 TLegend* lege2s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
1980 TLegendEntry* ent_a2s=lege2s->AddEntry(hfsigx100_1,"|#sigma res x| < 50 #mum","PL"); // P: marker, L; line
1981 ent_a2s->SetTextColor(hfsigx100_1->GetMarkerColor());
1982 TLegendEntry* ent_b2s=lege2s->AddEntry(hfsigx200_1,"|#sigma res x| < 70 #mum","PL"); // P: marker, L; line
1983 ent_b2s->SetTextColor(hfsigx200_1->GetMarkerColor());
1984 TLegendEntry* ent_c2s=lege2s->AddEntry(hfsigx300_1,"|#sigma res x| < 200 #mum","PL"); // P: marker, L; line
1985 ent_c2s->SetTextColor(hfsigx300_1->GetMarkerColor());
1986 lege2s->SetFillStyle(0);
1987 lege2s->Draw();
1988
1989 c8->cd(2);
1990 hfsigx100_2->SetTitle("fraction(#sigma x-residual < threshold) - layer 2");
1991 hfsigx100_2->SetLineColor(2);
1992 hfsigx100_2->SetMarkerStyle(31);
1993 hfsigx100_2->SetMarkerSize(2);
1994 hfsigx100_2->SetMarkerColor(2);
1995 hfsigx100_2->SetMinimum(0.);
1996 hfsigx100_2->SetMaximum(1.2);
1997 hfsigx100_2->GetXaxis()->SetTitle("run number");
1998 hfsigx100_2->GetYaxis()->SetTitleOffset(1.2);
1999 // hfsigx100_2->GetYaxis()->SetTitle("%");
2000 hfsigx100_2->Draw("P");
2001 hfsigx200_2->SetLineColor(4);
2002 hfsigx200_2->SetMarkerStyle(21);
2003 hfsigx200_2->SetMarkerColor(4);
2004 hfsigx200_2->Draw("same P");
2005 hfsigx300_2->SetLineColor(6);
2006 hfsigx300_2->SetMarkerStyle(22);
2007 hfsigx300_2->SetMarkerColor(6);
2008 hfsigx300_2->Draw("same P");
2009// l2->Draw();
2010
2011 TLegend* leggi2s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2012 TLegendEntry* en_ta2s=leggi2s->AddEntry(hfsigx100_2,"#sigma res x < 130 #mum","PL"); // P: marker, L; line
2013 en_ta2s->SetTextColor(hfsigx100_2->GetMarkerColor());
2014 TLegendEntry* en_tb2s=leggi2s->AddEntry(hfsigx200_2,"#sigma res x < 160 #mum","PL"); // P: marker, L; line
2015 en_tb2s->SetTextColor(hfsigx200_2->GetMarkerColor());
2016 TLegendEntry* en_tc2s=leggi2s->AddEntry(hfsigx300_2,"#sigma res x < 200 #mum","PL"); // P: marker, L; line
2017 en_tc2s->SetTextColor(hfsigx300_2->GetMarkerColor());
2018 leggi2s->SetFillStyle(0);
2019 leggi2s->Draw();
2020 c8->SaveAs("SPD_sigx");
2021// pdfFileNames+=" SPD_sigx";
2022 c8->Update();
2023
2024 // -------------------------- SPD z -------------------------------------
2025
2026 TCanvas* c9=new TCanvas("c9","SPD z res width", 1200,800);
2027 c9->Divide(1,2);
2028 c9->cd(1);
2029 hfsigz100_1->SetTitle("fraction(#sigma z-residual < threshold) - layer 1");
2030 hfsigz100_1->SetLineColor(2);
2031 hfsigz100_1->SetMarkerStyle(31);
2032 hfsigz100_1->SetMarkerSize(2);
2033 hfsigz100_1->SetMarkerColor(2);
2034 hfsigz100_1->SetMinimum(0.);
2035 hfsigz100_1->SetMaximum(1.2);
2036 hfsigz100_1->GetXaxis()->SetTitle("run number");
2037 hfsigz100_1->GetYaxis()->SetTitleOffset(1.2);
2038 // hfsigz100_1->GetYaxis()->SetTitle("%");
2039 hfsigz100_1->Draw("P");
2040 hfsigz300_1->SetLineColor(4);
2041 hfsigz300_1->SetMarkerStyle(21);
2042 hfsigz300_1->SetMarkerColor(4);
2043 hfsigz300_1->Draw("same P");
2044 hfsigz500_1->SetLineColor(6);
2045 hfsigz500_1->SetMarkerStyle(22);
2046 hfsigz500_1->SetMarkerColor(6);
2047 hfsigz500_1->Draw("same P");
2048// l1->Draw();
2049 TLegend* lege20s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2050 TLegendEntry* ent_a20s=lege20s->AddEntry(hfsigz100_1,"#sigma res z < 130 #mum","PL"); // P: marker, L; line
2051 ent_a20s->SetTextColor(hfsigz100_1->GetMarkerColor());
2052 TLegendEntry* ent_b20s=lege20s->AddEntry(hfsigz300_1,"#sigma res z < 160 #mum","PL"); // P: marker, L; line
2053 ent_b20s->SetTextColor(hfsigz300_1->GetMarkerColor());
2054 TLegendEntry* ent_c20s=lege20s->AddEntry(hfsigz500_1,"#sigma res z < 200 #mum","PL"); // P: marker, L; line
2055 ent_c20s->SetTextColor(hfsigz500_1->GetMarkerColor());
2056 lege20s->SetFillStyle(0);
2057 lege20s->Draw();
2058
2059 c9->cd(2);
2060 hfsigz100_2->SetTitle("fraction(#sigma z-residual < threshold) - layer 2");
2061 hfsigz100_2->SetLineColor(2);
2062 hfsigz100_2->SetMarkerStyle(31);
2063 hfsigz100_2->SetMarkerSize(2);
2064 hfsigz100_2->SetMarkerColor(2);
2065 hfsigz100_2->SetMinimum(0.);
2066 hfsigz100_2->SetMaximum(1.2);
2067 hfsigz100_2->GetXaxis()->SetTitle("run number");
2068 hfsigz100_2->GetYaxis()->SetTitleOffset(1.2);
2069 // hfsigz100_2->GetYaxis()->SetTitle("%");
2070 hfsigz100_2->Draw("P");
2071 hfsigz300_2->SetLineColor(4);
2072 hfsigz300_2->SetMarkerStyle(21);
2073 hfsigz300_2->SetMarkerColor(4);
2074 hfsigz300_2->Draw("same P");
2075 hfsigz500_2->SetLineColor(6);
2076 hfsigz500_2->SetMarkerStyle(22);
2077 hfsigz500_2->SetMarkerColor(6);
2078 hfsigz500_2->Draw("same P");
2079// l2->Draw();
2080 TLegend* leggi21s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2081 TLegendEntry* en_ta21s=leggi21s->AddEntry(hfsigz100_2,"#sigma res z < 180 #mum","PL"); // P: marker, L; line
2082 en_ta21s->SetTextColor(hfsigz100_2->GetMarkerColor());
2083 TLegendEntry* en_tb21s=leggi21s->AddEntry(hfsigz300_2,"#sigma res z < 190 #mum","PL"); // P: marker, L; line
2084 en_tb21s->SetTextColor(hfsigz300_2->GetMarkerColor());
2085 TLegendEntry* en_tc21s=leggi21s->AddEntry(hfsigz500_2,"#sigma res z < 250 #mum","PL"); // P: marker, L; line
2086 en_tc21s->SetTextColor(hfsigz500_2->GetMarkerColor());
2087 leggi21s->SetFillStyle(0);
2088 leggi21s->Draw();
2089 c9->SaveAs("SPD_sigz");
2090// pdfFileNames+=" SPD_sigz";
2091 c9->Update();
2092
2093 //---------------------------- SDD x ----------------------------------------
2094
2095 TCanvas* c10=new TCanvas("c10","SDD x res width", 1200,800);
2096 c10->Divide(1,2);
2097 c10->cd(1);
2098 hfsigx100_3->SetTitle("fraction(#sigma x-residual < threshold) - layer 3");
2099 hfsigx100_3->SetLineColor(2);
2100 hfsigx100_3->SetMarkerStyle(31);
2101 hfsigx100_3->SetMarkerSize(2);
2102 hfsigx100_3->SetMarkerColor(2);
2103 hfsigx100_3->SetMinimum(0.);
2104 hfsigx100_3->SetMaximum(1.2);
2105 hfsigx100_3->GetXaxis()->SetTitle("run number");
2106 hfsigx100_3->GetYaxis()->SetTitleOffset(1.2);
2107 // hfsigx100_3->GetYaxis()->SetTitle("%");
2108 hfsigx100_3->Draw("P");
2109 hfsigx200_3->SetLineColor(4);
2110 hfsigx200_3->SetMarkerStyle(21);
2111 hfsigx200_3->SetMarkerColor(4);
2112 hfsigx200_3->Draw("same P");
2113 hfsigx300_3->SetLineColor(6);
2114 hfsigx300_3->SetMarkerStyle(22);
2115 hfsigx300_3->SetMarkerColor(6);
2116 hfsigx300_3->Draw("same P");
2117// l3->Draw();
2118 TLegend* lege22s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2119 TLegendEntry* ent_a22s=lege22s->AddEntry(hfsigx100_3,"#sigma res x < 200 #mum","PL"); // P: marker, L; line
2120 ent_a22s->SetTextColor(hfsigx100_3->GetMarkerColor());
2121 TLegendEntry* ent_b22s=lege22s->AddEntry(hfsigx200_3,"#sigma res x < 250 #mum","PL"); // P: marker, L; line
2122 ent_b22s->SetTextColor(hfsigx200_3->GetMarkerColor());
2123 TLegendEntry* ent_c22s=lege22s->AddEntry(hfsigx300_3,"#sigma res x < 400 #mum","PL"); // P: marker, L; line
2124 ent_c22s->SetTextColor(hfsigx300_3->GetMarkerColor());
2125 lege22s->SetFillStyle(0);
2126 lege22s->Draw();
2127
2128 c10->cd(2);
2129 hfsigx100_4->SetTitle("fraction(#sigma x-residual < threshold) - layer 4");
2130 hfsigx100_4->SetLineColor(2);
2131 hfsigx100_4->SetMarkerStyle(31);
2132 hfsigx100_4->SetMarkerSize(2);
2133 hfsigx100_4->SetMarkerColor(2);
2134 hfsigx100_4->SetMinimum(0.);
2135 hfsigx100_4->SetMaximum(1.2);
2136 hfsigx100_4->GetXaxis()->SetTitle("run number");
2137 hfsigx100_4->GetYaxis()->SetTitleOffset(1.2);
2138 // hfsigx100_4->GetYaxis()->SetTitle("%");
2139 hfsigx100_4->Draw("P");
2140 hfsigx200_4->SetLineColor(4);
2141 hfsigx200_4->SetMarkerStyle(21);
2142 hfsigx200_4->SetMarkerColor(4);
2143 hfsigx200_4->Draw("same P");
2144 hfsigx300_4->SetLineColor(6);
2145 hfsigx300_4->SetMarkerStyle(22);
2146 hfsigx300_4->SetMarkerColor(6);
2147 hfsigx300_4->Draw("same P");
2148// l4->Draw();
2149 TLegend* leggi23s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2150 TLegendEntry* en_ta23s=leggi23s->AddEntry(hfsigx100_4,"#sigma res x < 200 #mum","PL"); // P: marker, L; line
2151 en_ta23s->SetTextColor(hfsigx100_4->GetMarkerColor());
2152 TLegendEntry* en_tb23s=leggi23s->AddEntry(hfsigx200_4,"#sigma res x < 250 #mum","PL"); // P: marker, L; line
2153 en_tb23s->SetTextColor(hfsigx200_4->GetMarkerColor());
2154 TLegendEntry* en_tc23s=leggi23s->AddEntry(hfsigx300_4,"#sigma res x < 400 #mum","PL"); // P: marker, L; line
2155 en_tc23s->SetTextColor(hfsigx300_4->GetMarkerColor());
2156 leggi23s->SetFillStyle(0);
2157 leggi23s->Draw();
2158 c10->SaveAs("SDD_sigx");
2159// pdfFileNames+=" SDD_sigx";
2160 c10->Update();
2161
2162 // -------------------------- SDD z -------------------------------------
2163
2164 TCanvas* c11=new TCanvas("c11","SDD z res width", 1200,800);
2165 c11->Divide(1,2);
2166 c11->cd(1);
2167 hfsigz100_3->SetTitle("fraction(#sigma z-residual < threshold) - layer 3");
2168 hfsigz100_3->SetLineColor(2);
2169 hfsigz100_3->SetMarkerStyle(31);
2170 hfsigz100_3->SetMarkerSize(2);
2171 hfsigz100_3->SetMarkerColor(2);
2172 hfsigz100_3->SetMinimum(0.);
2173 hfsigz100_3->SetMaximum(1.2);
2174 hfsigz100_3->GetXaxis()->SetTitle("run number");
2175 hfsigz100_3->GetYaxis()->SetTitleOffset(1.2);
2176 // hfsigz100_3->GetYaxis()->SetTitle("%");
2177 hfsigz100_3->Draw("P");
2178 hfsigz300_3->SetLineColor(4);
2179 hfsigz300_3->SetMarkerStyle(21);
2180 hfsigz300_3->SetMarkerColor(4);
2181 hfsigz300_3->Draw("same P");
2182 hfsigz500_3->SetLineColor(6);
2183 hfsigz500_3->SetMarkerStyle(22);
2184 hfsigz500_3->SetMarkerColor(6);
2185 hfsigz500_3->Draw("same P");
2186// l3->Draw();
2187 TLegend* lege24s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2188 TLegendEntry* ent_a24s=lege24s->AddEntry(hfsigz100_3,"#sigma res z < 150 #mum","PL"); // P: marker, L; line
2189 ent_a24s->SetTextColor(hfsigz100_3->GetMarkerColor());
2190 TLegendEntry* ent_b24s=lege24s->AddEntry(hfsigz300_3,"#sigma res z < 200 #mum","PL"); // P: marker, L; line
2191 ent_b24s->SetTextColor(hfsigz300_3->GetMarkerColor());
2192 TLegendEntry* ent_c24s=lege24s->AddEntry(hfsigz500_3,"#sigma res z < 350 #mum","PL"); // P: marker, L; line
2193 ent_c24s->SetTextColor(hfsigz500_3->GetMarkerColor());
2194 lege24s->SetFillStyle(0);
2195 lege24s->Draw();
2196
2197 c11->cd(2);
2198 hfsigz100_4->SetTitle("fraction(#sigma z-residual < threshold) - layer 4");
2199 hfsigz100_4->SetLineColor(2);
2200 hfsigz100_4->SetMarkerStyle(31);
2201 hfsigz100_4->SetMarkerSize(2);
2202 hfsigz100_4->SetMarkerColor(2);
2203 hfsigz100_4->SetMinimum(0.);
2204 hfsigz100_4->SetMaximum(1.2);
2205 hfsigz100_4->GetXaxis()->SetTitle("run number");
2206 hfsigz100_4->GetYaxis()->SetTitleOffset(1.2);
2207 // hfsigz100_4->GetYaxis()->SetTitle("%");
2208 hfsigz100_4->Draw("P");
2209 hfsigz300_4->SetLineColor(4);
2210 hfsigz300_4->SetMarkerStyle(21);
2211 hfsigz300_4->SetMarkerColor(4);
2212 hfsigz300_4->Draw("same P");
2213 hfsigz500_4->SetLineColor(6);
2214 hfsigz500_4->SetMarkerStyle(22);
2215 hfsigz500_4->SetMarkerColor(6);
2216 hfsigz500_4->Draw("same P");
2217// l4->Draw();
2218 TLegend* leggi25s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2219 TLegendEntry* en_ta25s=leggi25s->AddEntry(hfsigz100_4,"#sigma res z < 250 #mum","PL"); // P: marker, L; line
2220 en_ta25s->SetTextColor(hfsigz100_4->GetMarkerColor());
2221 TLegendEntry* en_tb25s=leggi25s->AddEntry(hfsigz300_4,"#sigma res z < 300 #mum","PL"); // P: marker, L; line
2222 en_tb25s->SetTextColor(hfsigz300_4->GetMarkerColor());
2223 TLegendEntry* en_tc25s=leggi25s->AddEntry(hfsigz500_4,"#sigma res | < 450 #mum","PL"); // P: marker, L; line
2224 en_tc25s->SetTextColor(hfsigz500_4->GetMarkerColor());
2225 leggi25s->SetFillStyle(0);
2226 leggi25s->Draw();
2227 c11->SaveAs("SDD_sigz");
2228// pdfFileNames+=" SDD_sigz";
2229 c11->Update();
2230
2231 //---------------------------- SSD x ----------------------------------------
2232
2233 TCanvas* c12=new TCanvas("c12","SSD x res width", 1200,800);
2234 c12->Divide(1,2);
2235 c12->cd(1);
2236 hfsigx100_5->SetTitle("fraction(#sigma x-residual < threshold) - layer 5");
2237 hfsigx100_5->SetLineColor(2);
2238 hfsigx100_5->SetMarkerStyle(31);
2239 hfsigx100_5->SetMarkerSize(2);
2240 hfsigx100_5->SetMarkerColor(2);
2241 hfsigx100_5->SetMinimum(0.);
2242 hfsigx100_5->SetMaximum(1.2);
2243 hfsigx100_5->GetXaxis()->SetTitle("run number");
2244 hfsigx100_5->GetYaxis()->SetTitleOffset(1.2);
2245 // hfsigx100_5->GetYaxis()->SetTitle("%");
2246 hfsigx100_5->Draw("P");
2247 hfsigx200_5->SetLineColor(4);
2248 hfsigx200_5->SetMarkerStyle(21);
2249 hfsigx200_5->SetMarkerColor(4);
2250 hfsigx200_5->Draw("same P");
2251 hfsigx300_5->SetLineColor(6);
2252 hfsigx300_5->SetMarkerStyle(22);
2253 hfsigx300_5->SetMarkerColor(6);
2254 hfsigx300_5->Draw("same P");
2255// l5->Draw();
2256 TLegend* lege26s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2257 TLegendEntry* ent_a26s=lege26s->AddEntry(hfsigx100_5,"#sigma res x < 70 #mum","PL"); // P: marker, L; line
2258 ent_a26s->SetTextColor(hfsigx100_5->GetMarkerColor());
2259 TLegendEntry* ent_b26s=lege26s->AddEntry(hfsigx200_5,"#sigma res x < 100 #mum","PL"); // P: marker, L; line
2260 ent_b26s->SetTextColor(hfsigx200_5->GetMarkerColor());
2261 TLegendEntry* ent_c26s=lege26s->AddEntry(hfsigx300_5,"#sigma res x < 300 #mum","PL"); // P: marker, L; line
2262 ent_c26s->SetTextColor(hfsigx300_5->GetMarkerColor());
2263 lege26s->SetFillStyle(0);
2264 lege26s->Draw();
2265
2266 c12->cd(2);
2267 hfsigx100_6->SetTitle("fraction(#sigma x-residual < threshold) - layer 6");
2268 hfsigx100_6->SetLineColor(2);
2269 hfsigx100_6->SetMarkerStyle(31);
2270 hfsigx100_6->SetMarkerSize(2);
2271 hfsigx100_6->SetMarkerColor(2);
2272 hfsigx100_6->SetMinimum(0.);
2273 hfsigx100_6->SetMaximum(1.2);
2274 hfsigx100_6->GetXaxis()->SetTitle("run number");
2275 hfsigx100_6->GetYaxis()->SetTitleOffset(1.2);
2276 // hfsigx100_6->GetYaxis()->SetTitle("%");
2277 hfsigx100_6->Draw("P");
2278 hfsigx200_6->SetLineColor(4);
2279 hfsigx200_6->SetMarkerStyle(21);
2280 hfsigx200_6->SetMarkerColor(4);
2281 hfsigx200_6->Draw("same P");
2282 hfsigx300_6->SetLineColor(6);
2283 hfsigx300_6->SetMarkerStyle(22);
2284 hfsigx300_6->SetMarkerColor(6);
2285 hfsigx300_6->Draw("same P");
2286// l6->Draw();
2287 TLegend* leggi27s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2288 TLegendEntry* en_ta27s=leggi27s->AddEntry(hfsigx100_6,"#sigma res x < 200 #mum","PL"); // P: marker, L; line
2289 en_ta27s->SetTextColor(hfsigx100_6->GetMarkerColor());
2290 TLegendEntry* en_tb27s=leggi27s->AddEntry(hfsigx200_6,"#sigma res x < 300 #mum","PL"); // P: marker, L; line
2291 en_tb27s->SetTextColor(hfsigx200_6->GetMarkerColor());
2292 TLegendEntry* en_tc27s=leggi27s->AddEntry(hfsigx300_6,"#sigma res x < 400 #mum","PL"); // P: marker, L; line
2293 en_tc27s->SetTextColor(hfsigx300_6->GetMarkerColor());
2294 leggi27s->SetFillStyle(0);
2295 leggi27s->Draw();
2296 c12->SaveAs("SSD_sigx");
2297// pdfFileNames+=" SSD_sigx";
2298 c12->Update();
2299
2300 // -------------------------- SSD z -------------------------------------
2301
2302 TCanvas* c13=new TCanvas("c13","SSD z res width", 1200,800);
2303 c13->Divide(1,2);
2304 c13->cd(1);
2305 hfsigz100_5->SetTitle("fraction(#sigma z-residual < threshold) - layer 5");
2306 hfsigz100_5->SetLineColor(2);
2307 hfsigz100_5->SetMarkerStyle(31);
2308 hfsigz100_5->SetMarkerSize(2);
2309 hfsigz100_5->SetMarkerColor(2);
2310 hfsigz100_5->SetMinimum(0.);
2311 hfsigz100_5->SetMaximum(1.2);
2312 hfsigz100_5->GetXaxis()->SetTitle("run number");
2313 hfsigz100_5->GetYaxis()->SetTitleOffset(1.2);
2314 // hfsigz100_5->GetYaxis()->SetTitle("%");
2315 hfsigz100_5->Draw("P");
2316 hfsigz300_5->SetLineColor(4);
2317 hfsigz300_5->SetMarkerStyle(21);
2318 hfsigz300_5->SetMarkerColor(4);
2319 hfsigz300_5->Draw("same P");
2320 hfsigz500_5->SetLineColor(6);
2321 hfsigz500_5->SetMarkerStyle(22);
2322 hfsigz500_5->SetMarkerColor(6);
2323 hfsigz500_5->Draw("same P");
2324// l5->Draw();
2325 TLegend* lege28s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2326 TLegendEntry* ent_a28s=lege28s->AddEntry(hfsigz100_5,"#sigma res z < 1000 #mum","PL"); // P: marker, L; line
2327 ent_a28s->SetTextColor(hfsigz100_5->GetMarkerColor());
2328 TLegendEntry* ent_b28s=lege28s->AddEntry(hfsigz300_5,"#sigma res z < 1100 #mum","PL"); // P: marker, L; line
2329 ent_b28s->SetTextColor(hfsigz300_5->GetMarkerColor());
2330 TLegendEntry* ent_c28s=lege28s->AddEntry(hfsigz500_5,"#sigma res z < 1500 #mum","PL"); // P: marker, L; line
2331 ent_c28s->SetTextColor(hfsigz500_5->GetMarkerColor());
2332 lege28s->SetFillStyle(0);
2333 lege28s->Draw();
2334
2335 c13->cd(2);
2336 hfsigz100_6->SetTitle("fraction(#sigma z-residual < threshold) - layer 6");
2337 hfsigz100_6->SetLineColor(2);
2338 hfsigz100_6->SetMarkerStyle(31);
2339 hfsigz100_6->SetMarkerSize(2);
2340 hfsigz100_6->SetMarkerColor(2);
2341 hfsigz100_6->SetMinimum(0.);
2342 hfsigz100_6->SetMaximum(1.2);
2343 hfsigz100_6->GetXaxis()->SetTitle("run number");
2344 hfsigz100_6->GetYaxis()->SetTitleOffset(1.2);
2345 // hfsigz100_6->GetYaxis()->SetTitle("%");
2346 hfsigz100_6->Draw("P");
2347 hfsigz300_6->SetLineColor(4);
2348 hfsigz300_6->SetMarkerStyle(21);
2349 hfsigz300_6->SetMarkerColor(4);
2350 hfsigz300_6->Draw("same P");
2351 hfsigz500_6->SetLineColor(6);
2352 hfsigz500_6->SetMarkerStyle(22);
2353 hfsigz500_6->SetMarkerColor(6);
2354 hfsigz500_6->Draw("same P");
2355// l6->Draw();
2356 TLegend* leggi29s=new TLegend(0.5,0.80,0.88,1.00); // x1, y1, x2, y2
2357 TLegendEntry* en_ta29s=leggi29s->AddEntry(hfsigz100_6,"#sigma res z < 1000 #mum","PL"); // P: marker, L; line
2358 en_ta29s->SetTextColor(hfsigz100_6->GetMarkerColor());
2359 TLegendEntry* en_tb29s=leggi29s->AddEntry(hfsigz300_6,"#sigma res z < 1100 #mum","PL"); // P: marker, L; line
2360 en_tb29s->SetTextColor(hfsigz300_6->GetMarkerColor());
2361 TLegendEntry* en_tc29s=leggi29s->AddEntry(hfsigz500_6,"#sigma res z < 1500 #mum","PL"); // P: marker, L; line
2362 en_tc29s->SetTextColor(hfsigz500_6->GetMarkerColor());
2363 leggi29s->SetFillStyle(0);
2364 leggi29s->Draw();
2365 c13->SaveAs("SSD_sigz");
2366// pdfFileNames+=" SSD_sigz";
2367 c13->Update();
2368
2369 // order pdf files
2370 pdfFileNames+=" SPD_resx";
2371 pdfFileNames+=" SPD_sigx";
2372 pdfFileNames+=" SPD_resz";
2373 pdfFileNames+=" SPD_sigz";
2374 pdfFileNames+=" SDD_resx_L3";
2375 pdfFileNames+=" SDD_resx_L4";
2376 pdfFileNames+=" SDD_sigx";
2377 pdfFileNames+=" SDD_resz";
2378 pdfFileNames+=" SDD_sigz";
2379 pdfFileNames+=" SSD_resx";
2380 pdfFileNames+=" SSD_sigx";
2381 pdfFileNames+=" SSD_resz";
2382 pdfFileNames+=" SSD_sigz";
2383
2384 // merge the pdf files
2385 TString command("gs -q -dNOPAUSE -dBATCH -sDEVICE=pdfwrite -sOutputFile=merged");
2386 command=command+"ITS_Align_trend_norm.pdf "+pdfFileNames;
2387 gSystem->Exec(command.Data());
2388 printf(" Merging the pdf file: %s \n",command.Data());
2389 delete [] myIndex;
2390 delete [] noRuns;
2391
2392}