1 |
/** |
/** |
2 |
* FTrkCalibQLookExpert.cpp |
* FTrkCalibQLookExpert.cxx |
3 |
* |
* |
4 |
* autor: D.Fedele |
* autor: D.Fedele |
5 |
* version 2.0 |
* version v1r16 |
6 |
* Parameters: |
* Parameters: |
7 |
* file - the data file to analyze |
* file - the data file to analyze |
8 |
* step - select =1 in order to analyze one event at time |
* step - select =1 in order to analyze one event at time |
18 |
#include <TPaveText.h> |
#include <TPaveText.h> |
19 |
#include <TLatex.h> |
#include <TLatex.h> |
20 |
#include <TCanvas.h> |
#include <TCanvas.h> |
21 |
|
#include <TFile.h> |
22 |
#include <TTree.h> |
#include <TTree.h> |
23 |
#include <TGraph.h> |
#include <TGraph.h> |
24 |
#include <TStyle.h> |
#include <TStyle.h> |
30 |
#include <CalibTrk2Event.h> |
#include <CalibTrk2Event.h> |
31 |
// |
// |
32 |
|
|
33 |
|
typedef struct caltrk_def{ |
|
void stringcopy(TString& s1, const TString& s2, Int_t from=0, Int_t to=0){ |
|
|
if ( to == 0 ){ |
|
|
Int_t t2length = s2.Length(); |
|
|
s1 = ""; |
|
|
to = t2length; |
|
|
}; |
|
|
for (Int_t i = from; i<to; i++){ |
|
|
s1.Append(s2[i],1); |
|
|
}; |
|
|
}; |
|
|
|
|
|
typedef struct caltrk_def{ |
|
34 |
Int_t good0[2]; |
Int_t good0[2]; |
35 |
Int_t daqmode[12]; |
Int_t daqmode[12]; |
36 |
Int_t dspnum[12]; |
Int_t dspnum[12]; |
58 |
{ |
{ |
59 |
// |
// |
60 |
// obtain information about the data file and select the output dir |
// obtain information about the data file and select the output dir |
61 |
const string filepath=file.Data(); |
Int_t dwpos = file.Last('/'); |
62 |
Int_t dwpos = filepath.rfind("DW_"); |
Int_t dwpos1 = file.Last('.'); |
63 |
Int_t dwpos1 = filepath.find(".root"); |
TString base,ffile ; |
64 |
TString fpath=(filepath.c_str()); |
ffile=file(dwpos+1,dwpos1-(dwpos+1)); |
65 |
TString base,ffile; |
if(dwpos>0) base=file(0,dwpos); |
|
stringcopy(base,fpath,0,dwpos); |
|
|
stringcopy(ffile,fpath,dwpos,dwpos1); |
|
66 |
|
|
67 |
TString out; |
TString out; |
68 |
if(outdir.Length()==0){ |
if(outdir.Length()==0){ |
70 |
}else{ |
}else{ |
71 |
out = outdir; |
out = outdir; |
72 |
} |
} |
73 |
|
if(out.Last('/')+1<out.Length()) out+="/"; |
74 |
|
|
75 |
// |
// |
76 |
// inizialise the variables and open the file |
// inizialise the variables and open the file |
77 |
Int_t nevents=0; |
Int_t nevents=0; |
78 |
Int_t minevent = 0; |
Int_t minevent = 0; |
79 |
Int_t maxevent = 0; |
Int_t maxevent = 0; |
80 |
ULong64_t FOBT[2], HOBT=0, TOBT=0; |
ULong64_t FOBT[2]; |
81 |
|
|
82 |
FOBT[0]=0; |
FOBT[0]=0; |
83 |
FOBT[1]=0; |
FOBT[1]=0; |
121 |
pamela::EventHeader *eh1=0,*eh2=0,*eh4=0,*eh3=0; |
pamela::EventHeader *eh1=0,*eh2=0,*eh4=0,*eh3=0; |
122 |
pamela::PscuHeader *ph=0; |
pamela::PscuHeader *ph=0; |
123 |
pamela::CalibTrk1Event *trk1 = 0; |
pamela::CalibTrk1Event *trk1 = 0; |
124 |
pamela::CalibTrk2Event *trk2 = 0; |
pamela::CalibTrk2Event *trk2 = 0; |
125 |
|
pamela::EventCounter *cod=0; |
126 |
|
|
127 |
|
pamela::PacketType *pctp=0; |
128 |
|
|
129 |
hotr = (TTree*)calibFile->Get("CalibHeader"); |
hotr = (TTree*)calibFile->Get("CalibHeader"); |
130 |
hotr->SetBranchAddress("Header", &eh4); |
hotr->SetBranchAddress("Header", &eh4); |
176 |
} else if (toevent > nevents) { |
} else if (toevent > nevents) { |
177 |
maxevent = nevents; |
maxevent = nevents; |
178 |
} else { |
} else { |
179 |
maxevent = fromevent; |
maxevent = nevents; |
180 |
}; |
}; |
181 |
}; |
}; |
182 |
|
|
186 |
// |
// |
187 |
//********************************************************************** |
//********************************************************************** |
188 |
|
|
189 |
|
Int_t hcevent=hotr->GetEntries(); |
190 |
|
Int_t tcevent=totr->GetEntries(); |
191 |
|
ULong64_t HOBT[hcevent], TOBT[tcevent]; |
192 |
|
for (Int_t i = 0; i < hcevent; i++){ |
193 |
|
hotr->GetEntry(i); |
194 |
|
ph = eh4->GetPscuHeader(); |
195 |
|
HOBT[i]= ph->GetOrbitalTime(); |
196 |
|
} |
197 |
|
for (Int_t i = 0; i < tcevent; i++){ |
198 |
|
totr->GetEntry(i); |
199 |
|
ph = eh3->GetPscuHeader(); |
200 |
|
TOBT[i]= ph->GetOrbitalTime(); |
201 |
|
} |
202 |
|
|
203 |
printf("\n Scan of calibration packets from %i to %i ... \n\n",minevent+1,maxevent); |
printf("\n Scan of calibration packets from %i to %i ... \n\n",minevent+1,maxevent); |
204 |
for (Int_t i = minevent; i < maxevent; i++){ |
for (Int_t i = minevent; i < maxevent; i++){ |
205 |
|
|
206 |
otr1->GetEntry(i); |
otr1->GetEntry(i); |
207 |
otr2->GetEntry(i); |
otr2->GetEntry(i); |
208 |
totr->GetEntry(i); |
|
|
hotr->GetEntry(i); |
|
|
ph = eh4->GetPscuHeader(); |
|
|
HOBT= ph->GetOrbitalTime(); |
|
|
ph = eh3->GetPscuHeader(); |
|
|
TOBT= ph->GetOrbitalTime(); |
|
|
|
|
209 |
ctrk.good0[0]=trk1->good0; |
ctrk.good0[0]=trk1->good0; |
210 |
ctrk.good0[1]=trk2->good0; |
ctrk.good0[1]=trk2->good0; |
211 |
for (Int_t m = 0; m < 6; m++){ |
for (Int_t m = 0; m < 6; m++){ |
212 |
ph = eh1->GetPscuHeader(); |
ph = eh1->GetPscuHeader(); |
213 |
|
cod = eh1->GetCounter(); |
214 |
FOBT[0]= ph->GetOrbitalTime(); |
FOBT[0]= ph->GetOrbitalTime(); |
215 |
ctrk.daqmode[trk1->DSPnumber[m]-1]=trk1->DAQmode[m]; |
if(trk1->DSPnumber[m]>0){ |
216 |
ctrk.dspnum[trk1->DSPnumber[m]-1]=trk1->DSPnumber[m]; |
if(trk1->DSPnumber[m]<13){ |
217 |
ctrk.calibnum[trk1->DSPnumber[m]-1]=trk1->calibnumber[m]; |
ctrk.daqmode[trk1->DSPnumber[m]-1]=trk1->DAQmode[m]; |
218 |
ctrk.ncalev[trk1->DSPnumber[m]-1]=trk1->ncalib_event[m]; |
ctrk.dspnum[trk1->DSPnumber[m]-1]=trk1->DSPnumber[m]; |
219 |
ctrk.ped1[trk1->DSPnumber[m]-1]=trk1->ped_l1[m]; |
ctrk.calibnum[trk1->DSPnumber[m]-1]=trk1->calibnumber[m]; |
220 |
ctrk.ped2[trk1->DSPnumber[m]-1]=trk1->ped_l2[m]; |
ctrk.ncalev[trk1->DSPnumber[m]-1]=trk1->ncalib_event[m]; |
221 |
ctrk.ped3[trk1->DSPnumber[m]-1]=trk1->ped_l3[m]; |
ctrk.ped1[trk1->DSPnumber[m]-1]=trk1->ped_l1[m]; |
222 |
ctrk.sig1[trk1->DSPnumber[m]-1]=trk1->sig_l1[m]; |
ctrk.ped2[trk1->DSPnumber[m]-1]=trk1->ped_l2[m]; |
223 |
ctrk.sig2[trk1->DSPnumber[m]-1]=trk1->sig_l2[m]; |
ctrk.ped3[trk1->DSPnumber[m]-1]=trk1->ped_l3[m]; |
224 |
ctrk.sig3[trk1->DSPnumber[m]-1]=trk1->sig_l3[m]; |
ctrk.sig1[trk1->DSPnumber[m]-1]=trk1->sig_l1[m]; |
225 |
ctrk.nbad1[trk1->DSPnumber[m]-1]=trk1->nbad_l1[m]; |
ctrk.sig2[trk1->DSPnumber[m]-1]=trk1->sig_l2[m]; |
226 |
ctrk.nbad2[trk1->DSPnumber[m]-1]=trk1->nbad_l2[m]; |
ctrk.sig3[trk1->DSPnumber[m]-1]=trk1->sig_l3[m]; |
227 |
ctrk.nbad3[trk1->DSPnumber[m]-1]=trk1->nbad_l3[m]; |
ctrk.nbad1[trk1->DSPnumber[m]-1]=trk1->nbad_l1[m]; |
228 |
ctrk.calfl[trk1->DSPnumber[m]-1]=trk1->cal_flag[m]; |
ctrk.nbad2[trk1->DSPnumber[m]-1]=trk1->nbad_l2[m]; |
229 |
ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk1->crc_cal[m][0]; |
ctrk.nbad3[trk1->DSPnumber[m]-1]=trk1->nbad_l3[m]; |
230 |
ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk1->crc_cal[m][1]; |
ctrk.calfl[trk1->DSPnumber[m]-1]=trk1->cal_flag[m]; |
231 |
ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk1->crc_cal[m][2]; |
ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk1->crc_cal[m][0]; |
232 |
ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk1->crc_hcal[m]; |
ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk1->crc_cal[m][1]; |
233 |
for (Int_t j = 0; j < 3072; j++){ |
ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk1->crc_cal[m][2]; |
234 |
ctrk.dspped[trk1->DSPnumber[m]-1][j]=trk1->DSPped_par[m][j]; |
ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk1->crc_hcal[m]; |
235 |
ctrk.dspsig[trk1->DSPnumber[m]-1][j]=trk1->DSPsig_par[m][j]; |
for (Int_t j = 0; j < 3072; j++){ |
236 |
ctrk.dspbad[trk1->DSPnumber[m]-1][j]=trk1->DSPbad_par[m][j]; |
ctrk.dspped[trk1->DSPnumber[m]-1][j]=trk1->DSPped_par[m][j]; |
237 |
}; |
ctrk.dspsig[trk1->DSPnumber[m]-1][j]=trk1->DSPsig_par[m][j]; |
238 |
|
ctrk.dspbad[trk1->DSPnumber[m]-1][j]=trk1->DSPbad_par[m][j]; |
239 |
|
}; |
240 |
|
} |
241 |
|
} |
242 |
ph = eh2->GetPscuHeader(); |
ph = eh2->GetPscuHeader(); |
243 |
FOBT[1]= ph->GetOrbitalTime(); |
FOBT[1]= ph->GetOrbitalTime(); |
244 |
ctrk.daqmode[trk2->DSPnumber[m]-1]=trk2->DAQmode[m]; |
if(trk2->DSPnumber[m]>0){ |
245 |
ctrk.dspnum[trk2->DSPnumber[m]-1]=trk2->DSPnumber[m]; |
if(trk2->DSPnumber[m]<13){ |
246 |
ctrk.calibnum[trk2->DSPnumber[m]-1]=trk2->calibnumber[m]; |
ctrk.daqmode[trk2->DSPnumber[m]-1]=trk2->DAQmode[m]; |
247 |
ctrk.ncalev[trk2->DSPnumber[m]-1]=trk2->ncalib_event[m]; |
ctrk.dspnum[trk2->DSPnumber[m]-1]=trk2->DSPnumber[m]; |
248 |
ctrk.ped1[trk2->DSPnumber[m]-1]=trk2->ped_l1[m]; |
ctrk.calibnum[trk2->DSPnumber[m]-1]=trk2->calibnumber[m]; |
249 |
ctrk.ped2[trk2->DSPnumber[m]-1]=trk2->ped_l2[m]; |
ctrk.ncalev[trk2->DSPnumber[m]-1]=trk2->ncalib_event[m]; |
250 |
ctrk.ped3[trk2->DSPnumber[m]-1]=trk2->ped_l3[m]; |
ctrk.ped1[trk2->DSPnumber[m]-1]=trk2->ped_l1[m]; |
251 |
ctrk.sig1[trk2->DSPnumber[m]-1]=trk2->sig_l1[m]; |
ctrk.ped2[trk2->DSPnumber[m]-1]=trk2->ped_l2[m]; |
252 |
ctrk.sig2[trk2->DSPnumber[m]-1]=trk2->sig_l2[m]; |
ctrk.ped3[trk2->DSPnumber[m]-1]=trk2->ped_l3[m]; |
253 |
ctrk.sig3[trk2->DSPnumber[m]-1]=trk2->sig_l3[m]; |
ctrk.sig1[trk2->DSPnumber[m]-1]=trk2->sig_l1[m]; |
254 |
ctrk.nbad1[trk2->DSPnumber[m]-1]=trk2->nbad_l1[m]; |
ctrk.sig2[trk2->DSPnumber[m]-1]=trk2->sig_l2[m]; |
255 |
ctrk.nbad2[trk2->DSPnumber[m]-1]=trk2->nbad_l2[m]; |
ctrk.sig3[trk2->DSPnumber[m]-1]=trk2->sig_l3[m]; |
256 |
ctrk.nbad3[trk2->DSPnumber[m]-1]=trk2->nbad_l3[m]; |
ctrk.nbad1[trk2->DSPnumber[m]-1]=trk2->nbad_l1[m]; |
257 |
ctrk.calfl[trk2->DSPnumber[m]-1]=trk2->cal_flag[m]; |
ctrk.nbad2[trk2->DSPnumber[m]-1]=trk2->nbad_l2[m]; |
258 |
ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk2->crc_cal[m][0]; |
ctrk.nbad3[trk2->DSPnumber[m]-1]=trk2->nbad_l3[m]; |
259 |
ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk2->crc_cal[m][1]; |
ctrk.calfl[trk2->DSPnumber[m]-1]=trk2->cal_flag[m]; |
260 |
ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk2->crc_cal[m][2]; |
ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk2->crc_cal[m][0]; |
261 |
ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk2->crc_hcal[m]; |
ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk2->crc_cal[m][1]; |
262 |
for (Int_t j = 0; j < 3072; j++){ |
ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk2->crc_cal[m][2]; |
263 |
ctrk.dspped[trk2->DSPnumber[m]-1][j]=trk2->DSPped_par[m][j]; |
ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk2->crc_hcal[m]; |
264 |
ctrk.dspsig[trk2->DSPnumber[m]-1][j]=trk2->DSPsig_par[m][j]; |
for (Int_t j = 0; j < 3072; j++){ |
265 |
ctrk.dspbad[trk2->DSPnumber[m]-1][j]=trk2->DSPbad_par[m][j]; |
ctrk.dspped[trk2->DSPnumber[m]-1][j]=trk2->DSPped_par[m][j]; |
266 |
}; |
ctrk.dspsig[trk2->DSPnumber[m]-1][j]=trk2->DSPsig_par[m][j]; |
267 |
}; |
ctrk.dspbad[trk2->DSPnumber[m]-1][j]=trk2->DSPbad_par[m][j]; |
268 |
|
} |
269 |
|
} |
270 |
|
} |
271 |
|
} |
272 |
|
|
273 |
|
|
274 |
// |
// |
276 |
Int_t risposta=0; |
Int_t risposta=0; |
277 |
stringstream fromfile; |
stringstream fromfile; |
278 |
|
|
279 |
fromfile<<"FTrkCalibQLook_EXPERT File: "<<ffile<<" -- CalibHeader OBT= "<<HOBT<<" -- Calib pkt OBT= "<<FOBT[0]<<" -- CalibTrailer OBT= "<<TOBT<<" --"; |
fromfile<<"FTrkCalibQLook_EXPERT File: "<<ffile<<" -- CalibHeader OBT= "<<HOBT[(cod->Get(pctp->CalibHeader))-1]<<" -- Calib pkt OBT= "<<FOBT[0]<<" -- CalibTrailer OBT= "<<TOBT[(cod->Get(pctp->CalibTrailer))]<<" --"; |
280 |
|
|
281 |
gStyle->SetLabelSize(0.08,"x"); |
gStyle->SetLabelSize(0.07,"x"); |
282 |
gStyle->SetLabelSize(0.08,"y"); |
gStyle->SetLabelSize(0.07,"y"); |
283 |
gStyle->SetTitleFillColor(10); |
gStyle->SetTitleFillColor(10); |
284 |
gStyle->SetTitleFontSize(0.1); |
gStyle->SetTitleFontSize(0.08); |
285 |
gStyle->SetTitleOffset(0.8,"y"); |
gStyle->SetTitleOffset(0.8,"y"); |
286 |
gStyle->SetTitleOffset(1.,"x"); |
gStyle->SetTitleOffset(0.9,"x"); |
287 |
gStyle->SetTitleSize(0.06,"y"); |
gStyle->SetTitleSize(0.06,"y"); |
288 |
gStyle->SetTitleSize(0.06,"x"); |
gStyle->SetTitleSize(0.06,"x"); |
289 |
gStyle->SetOptStat(0); |
gStyle->SetOptStat(101110); |
290 |
|
gStyle->SetStatX(0.9); |
291 |
|
gStyle->SetStatW(0.4); |
292 |
|
gStyle->SetStatColor(10); |
293 |
|
gStyle->SetStatFontSize(0.1); |
294 |
|
|
295 |
// |
// |
296 |
// draw display area |
// draw display area |
297 |
|
|
298 |
TLatex *tzz=new TLatex(); |
TLatex *tzz=new TLatex(); |
299 |
tzz->SetTextFont(32); |
tzz->SetTextFont(32); |
300 |
tzz->SetTextColor(1); |
tzz->SetTextColor(1); |
311 |
tzz->DrawLatex(.01,0.98,fromfile.str().c_str()); |
tzz->DrawLatex(.01,0.98,fromfile.str().c_str()); |
312 |
tzz->DrawLatex(.90,0.98,"SIGMA"); |
tzz->DrawLatex(.90,0.98,"SIGMA"); |
313 |
|
|
314 |
|
|
315 |
|
TCanvas *sig=new TCanvas("sig","FTrkCalibQLook_EXPERT_histosig",canvasx,canvasy); |
316 |
|
sig->SetFillColor(10); |
317 |
|
tzz->DrawLatex(.01,0.98,fromfile.str().c_str()); |
318 |
|
tzz->DrawLatex(.85,0.97,"Histograms of the sigmas"); |
319 |
|
|
320 |
|
|
321 |
|
|
322 |
// draw pads |
// draw pads |
323 |
TPad *trkpad1[12],*trkpad2[12]; //pad for histos |
TPad *trkpad1[12],*trkpad2[12],*trkpad3[36]; //pad for histos |
324 |
TPaveText *trkpadtext[12]; //pad for header |
TPaveText *trkpadtext[12]; //pad for header |
325 |
TH1F *histosig[12]; //histos of sigma |
TH1F *histosig[12]; //histos of sigma |
326 |
TH1F *histoped[12]; //histos of pedestals |
TH1F *histoped[12]; //histos of pedestals |
327 |
TH1F *histoasig[12]; //histos of sigma |
TH1F *histoasig[12]; //histos of sigma |
328 |
TH1F *histoaped[12]; //histos of pedestals |
TH1F *histoaped[12]; //histos of pedestals |
329 |
|
|
330 |
Double_t posy = 0.95; // up y-coord - top pads |
TH1F *histosiglad[12][3]; //histos of sigma |
|
Double_t hpad = 0.15; // pad height |
|
|
Double_t posx1=0; // left x-coord - pad column |
|
|
Double_t posx2=0; // right x-coord - pad olumn |
|
|
Double_t posx0=0; // x-coord - column division |
|
|
Double_t wrel = 0.6; // relative x size of first sub-column |
|
|
Double_t marg = 0.004; // margin among pads |
|
331 |
stringstream title; |
stringstream title; |
332 |
stringstream hid; |
stringstream hid; |
333 |
|
|
334 |
|
Float_t posy = 0.95; // up y-coord - top pads |
335 |
|
Float_t hpad = 0.15; // pad height |
336 |
|
Float_t posx1=0; // left x-coord - pad column |
337 |
|
Float_t posx2=0; // right x-coord - pad olumn |
338 |
|
Float_t posx0=0; // x-coord - column division |
339 |
|
Float_t wrel = 0.6; // relative x size of first sub-column |
340 |
|
Float_t marg = 0.004; // margin among pads |
341 |
|
|
342 |
|
|
343 |
for(Int_t n = 0; n<12; n++){ |
for(Int_t n = 0; n<12; n++){ |
344 |
if ( (n+1)%2 ) { |
if ( (n+1)%2 ) { |
345 |
if(n>1)posy = posy-(marg*2+hpad); |
if(n>1)posy = posy-(marg*2+hpad); |
346 |
posx1 = marg; |
posx1 = marg; |
347 |
posx2 = 0.5 - marg; |
posx2 = 0.5 - marg; |
348 |
posx0 = 0.5*wrel; |
posx0 = 0.5*wrel; |
349 |
|
|
350 |
} else { |
} else { |
351 |
posx1 = posx1 + 0.5; |
posx1 = posx1 + 0.5; |
352 |
posx2 = posx2 + 0.5; |
posx2 = posx2 + 0.5; |
353 |
posx0 = posx0 + 0.5; |
posx0 = posx0 + 0.5; |
354 |
|
|
355 |
}; |
}; |
356 |
/* -----------> pad for histograms */ |
/* -----------> pad for histograms */ |
357 |
trkpad1[n] = new TPad("pad1"," ",posx1,posy-hpad,posx0-marg,posy,18,0,0); |
trkpad1[n] = new TPad("pad1"," ",posx1,posy-hpad,posx0-marg,posy,18,0,0); |
361 |
/* -----------> HISTOGRAMS */ |
/* -----------> HISTOGRAMS */ |
362 |
/* calibration parameters */ |
/* calibration parameters */ |
363 |
title<<"DSP "<<n+1; |
title<<"DSP "<<n+1; |
364 |
hid<<"h"<<n; |
hid<<"h"<<n<<"i"<<i; |
365 |
histosig[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3072,0.5,3072.5); |
histosig[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3072,0.5,3072.5); |
366 |
hid.str(""); |
hid.str(""); |
367 |
hid<<"hh"<<n; |
hid<<"hh"<<n<<"i"<<i; |
368 |
histoped[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3072,0.5,3072.5); |
histoped[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3072,0.5,3072.5); |
369 |
hid.str(""); |
hid.str(""); |
370 |
hid<<"hhh"<<n; |
hid<<"hhh"<<n<<"i"<<i; |
|
title.str(""); |
|
371 |
hid.str(""); |
hid.str(""); |
372 |
/* AVERAGE calibration parameters */ |
/* AVERAGE calibration parameters */ |
373 |
hid<<"ah"<<n; |
hid<<"ah"<<n<<"i"<<i; |
374 |
histoasig[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3,0.5,3072.5); |
histoasig[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3,0.5,3072.5); |
375 |
hid.str(""); |
hid.str(""); |
376 |
hid<<"ahh"<<n; |
hid<<"ahh"<<n<<"i"<<i; |
377 |
histoaped[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3,0.5,3072.5); |
histoaped[n] = new TH1F(hid.str().c_str(),title.str().c_str(),3,0.5,3072.5); |
378 |
hid.str(""); |
hid.str(""); |
379 |
|
for(int ii=0;ii<3;ii++){ |
380 |
|
title.str(""); |
381 |
|
title<<"DSP "<<n+1<<" / Lad "<<ii+1; |
382 |
|
hid<<"hhhh"<<n<<"i"<<i<<"ii"<<ii; |
383 |
|
histosiglad[n][ii] = new TH1F(hid.str().c_str(),title.str().c_str(),32,-0.5,30.5); |
384 |
|
hid.str(""); |
385 |
|
} |
386 |
|
title.str(""); |
387 |
}; //end loop on views |
}; //end loop on views |
388 |
|
|
389 |
|
Float_t tposy = 0.95; // up y-coord - top pads |
390 |
|
Float_t thpad = 0.; // pad height |
391 |
|
Float_t tposx1=0; // left x-coord - pad column |
392 |
|
Float_t tposx0=0; // x-coord - column division |
393 |
|
Float_t twrel = 0.; // relative x size of first sub-column |
394 |
|
Float_t tmarg = 0.002; // margin among pads |
395 |
|
thpad = (tposy-tmarg*11)/6; |
396 |
|
twrel = (1-tmarg*12)/6; |
397 |
|
|
398 |
|
for(Int_t n = 0; n<36; n++){ |
399 |
|
if ( (n+1)%6==1 ) { |
400 |
|
if(n>1) tposy = tposy-(tmarg*2+thpad); |
401 |
|
tposx1 = tmarg; |
402 |
|
tposx0 = tposx1 + twrel; |
403 |
|
} else { |
404 |
|
tposx1 = tposx0 + 2*tmarg; |
405 |
|
tposx0 = tposx1 + twrel; |
406 |
|
} |
407 |
|
trkpad3[n]= new TPad("pad3"," ",tposx1,tposy-thpad,tposx0,tposy,18,0,0); |
408 |
|
} |
409 |
|
|
410 |
|
|
411 |
stringstream message; |
stringstream message; |
412 |
|
|
413 |
//-------------------------------- |
//-------------------------------- |
414 |
//CHECK CALIBRATION procedure |
//CHECK CALIBRATION procedure |
415 |
//-------------------------------- |
//-------------------------------- |
417 |
Int_t calok = 0;//BAD |
Int_t calok = 0;//BAD |
418 |
for(Int_t n = 0; n<12; n++){ |
for(Int_t n = 0; n<12; n++){ |
419 |
if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD |
if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD |
420 |
|
|
421 |
|
if(ctrk.dspnum[n]==0) |
422 |
|
continue; |
423 |
|
|
424 |
nn=ctrk.dspnum[n]-1; |
nn=ctrk.dspnum[n]-1; |
425 |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
/*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.* |
426 |
* |
* |
493 |
/******************************************************/ |
/******************************************************/ |
494 |
/* fill histos */ |
/* fill histos */ |
495 |
for(Int_t j = 0; j < 3072; j++){ |
for(Int_t j = 0; j < 3072; j++){ |
496 |
histosig[nn]->Fill((Float_t)j,ctrk.dspsig[nn][j]); |
histosig[nn]->Fill(j,ctrk.dspsig[n][j]); |
497 |
histoped[nn]->Fill((Float_t)j,ctrk.dspped[nn][j]); |
histoped[nn]->Fill(j,ctrk.dspped[n][j]); |
498 |
|
if(j<1024) histosiglad[nn][0]->Fill(ctrk.dspsig[n][j]); |
499 |
|
if(j>=1024 && j<2048) histosiglad[nn][1]->Fill(ctrk.dspsig[n][j]); |
500 |
|
if(j>=2048 && j<3072) histosiglad[nn][2]->Fill(ctrk.dspsig[n][j]); |
501 |
}; |
}; |
502 |
histoasig[nn]->Fill(1,ctrk.sig1[nn]); |
histoasig[nn]->Fill(1,ctrk.sig1[n]); |
503 |
histoasig[nn]->Fill(1025,ctrk.sig2[nn]); |
histoasig[nn]->Fill(1025,ctrk.sig2[n]); |
504 |
histoasig[nn]->Fill(2049,ctrk.sig3[nn]); |
histoasig[nn]->Fill(2049,ctrk.sig3[n]); |
505 |
histoaped[nn]->Fill(1,ctrk.ped1[nn]); |
histoaped[nn]->Fill(1,ctrk.ped1[nn]); |
506 |
histoaped[nn]->Fill(1025,ctrk.ped2[nn]); |
histoaped[nn]->Fill(1025,ctrk.ped2[n]); |
507 |
histoaped[nn]->Fill(2049,ctrk.ped3[nn]); |
histoaped[nn]->Fill(2049,ctrk.ped3[n]); |
508 |
/******************************************************/ |
/******************************************************/ |
509 |
|
|
510 |
TLine li; |
TLine li; |
512 |
li.SetLineStyle(3); |
li.SetLineStyle(3); |
513 |
li.SetLineWidth(2); |
li.SetLineWidth(2); |
514 |
|
|
515 |
Double_t maxhist=0; |
Float_t maxhist=0; |
516 |
TBox b; |
TBox b; |
|
b.SetFillColor(6); |
|
|
b.SetFillStyle(3945); |
|
517 |
/* plot PEDESTAL */ |
/* plot PEDESTAL */ |
518 |
c1->cd(); |
c1->cd(); |
519 |
trkpadtext[nn]->Draw(); |
trkpadtext[nn]->Draw(); |
521 |
trkpad1[nn]->cd(); |
trkpad1[nn]->cd(); |
522 |
trkpad1[nn]->SetFillColor(10); |
trkpad1[nn]->SetFillColor(10); |
523 |
trkpad1[nn]->SetFrameFillColor(10); |
trkpad1[nn]->SetFrameFillColor(10); |
524 |
|
histoped[nn]->SetStats(kFALSE); |
525 |
histoped[nn]->SetLineColor(1); |
histoped[nn]->SetLineColor(1); |
526 |
histoped[nn]->SetFillColor(12); |
histoped[nn]->SetFillColor(12); |
527 |
histoped[nn]->SetLineWidth(1); |
histoped[nn]->SetLineWidth(1); |
531 |
if((nn+1)%2==0) histoped[nn]->GetYaxis()->SetRangeUser(700,1700); |
if((nn+1)%2==0) histoped[nn]->GetYaxis()->SetRangeUser(700,1700); |
532 |
histoaped[nn]->SetLineColor(5); |
histoaped[nn]->SetLineColor(5); |
533 |
histoaped[nn]->SetLineWidth(1); |
histoaped[nn]->SetLineWidth(1); |
534 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
535 |
histoped[nn]->Draw("b"); |
histoped[nn]->Draw("b"); |
536 |
if(nn==1){ |
if(nn==0){ |
537 |
maxhist=histoped[nn]->GetMaximum(); |
maxhist=histoped[nn]->GetMaximum(); |
538 |
b.DrawBox(2816.,700.,3060.,maxhist); |
b.SetFillColor(107); |
539 |
|
b.SetFillStyle(3945); |
540 |
|
b.DrawBox(768.,700.,2047.,maxhist); |
541 |
|
} |
542 |
|
else if(nn==1){ |
543 |
|
maxhist=histoped[nn]->GetMaximum(); |
544 |
|
b.SetFillColor(6); |
545 |
|
b.SetFillStyle(3945); |
546 |
|
b.DrawBox(2944.,700.,3060.,maxhist); |
547 |
|
|
548 |
|
b.SetFillColor(107); |
549 |
|
b.SetFillStyle(3954); |
550 |
|
b.DrawBox(2816.,700.,2944.,maxhist); |
551 |
|
b.DrawBox(2048.,700.,2176.,maxhist); |
552 |
|
} |
553 |
|
else if(nn==4){ |
554 |
|
maxhist=histoped[nn]->GetMaximum(); |
555 |
|
b.SetFillColor(107); |
556 |
|
b.SetFillStyle(3954); |
557 |
|
b.DrawBox(384.,2200.,512.,maxhist); |
558 |
} |
} |
559 |
else if(nn==6){ |
else if(nn==6){ |
560 |
maxhist=histoped[nn]->GetMaximum(); |
maxhist=histoped[nn]->GetMaximum(); |
561 |
|
b.SetFillColor(6); |
562 |
|
b.SetFillStyle(3945); |
563 |
b.DrawBox(2560.,2200.,2816.,maxhist); |
b.DrawBox(2560.,2200.,2816.,maxhist); |
564 |
|
b.DrawBox(1024.,2200.,1535.,maxhist); |
565 |
|
|
566 |
|
b.SetFillColor(107); |
567 |
|
b.SetFillStyle(3954); |
568 |
b.DrawBox(512.,2200.,768.,maxhist); |
b.DrawBox(512.,2200.,768.,maxhist); |
569 |
b.DrawBox(1024.,2200.,1792.,maxhist); |
b.DrawBox(1536.,2200.,1792.,maxhist); |
570 |
|
} |
571 |
|
else if(nn==7){ |
572 |
|
maxhist=histoped[nn]->GetMaximum(); |
573 |
|
b.SetFillColor(107); |
574 |
|
b.SetFillStyle(3954); |
575 |
|
b.DrawBox(512.,700.,768.,maxhist); |
576 |
|
} |
577 |
|
else if(nn==9){ |
578 |
|
maxhist=histoped[nn]->GetMaximum(); |
579 |
|
b.SetFillColor(107); |
580 |
|
b.SetFillStyle(3954); |
581 |
|
b.DrawBox(256.,700.,384.,maxhist); |
582 |
} |
} |
583 |
else if(nn==11){ |
else if(nn==11){ |
584 |
maxhist=histoped[nn]->GetMaximum(); |
maxhist=histoped[nn]->GetMaximum(); |
585 |
|
b.SetFillColor(6); |
586 |
|
b.SetFillStyle(3945); |
587 |
b.DrawBox(768.,700.,1024.,maxhist); |
b.DrawBox(768.,700.,1024.,maxhist); |
588 |
|
|
589 |
|
b.SetFillColor(107); |
590 |
|
b.SetFillStyle(3954); |
591 |
|
b.DrawBox(0.,700.,512.,maxhist); |
592 |
|
b.DrawBox(1920.,700.,2560.,maxhist); |
593 |
} |
} |
594 |
} |
// } |
595 |
else histoped[nn]->Draw("axis"); |
// else histoped[nn]->Draw("axis"); |
596 |
histoaped[nn]->Draw("same"); |
histoaped[nn]->Draw("same"); |
597 |
if((nn+1)%2==1) { |
if((nn+1)%2==1) { |
598 |
li.DrawLine(1024.5,2200,1024.5,3200); |
li.DrawLine(1024.5,2200,1024.5,3200); |
602 |
li.DrawLine(1024.5,700,1024.5,1700); |
li.DrawLine(1024.5,700,1024.5,1700); |
603 |
li.DrawLine(2048.5,700,2048.5,1700); |
li.DrawLine(2048.5,700,2048.5,1700); |
604 |
} |
} |
605 |
|
c1->Update();//draw pads in canvas |
606 |
|
|
607 |
|
|
608 |
|
|
609 |
/* plot SIGMA */ |
/* plot SIGMA */ |
610 |
Double_t max=500.; |
Float_t max=500.; |
611 |
c2->cd(); |
c2->cd(); |
612 |
trkpadtext[nn]->Draw(); |
trkpadtext[nn]->Draw(); |
613 |
trkpad2[nn]->SetLogy(); |
trkpad2[nn]->SetLogy(); |
615 |
trkpad2[nn]->cd(); |
trkpad2[nn]->cd(); |
616 |
trkpad2[nn]->SetFillColor(10); |
trkpad2[nn]->SetFillColor(10); |
617 |
trkpad2[nn]->SetFrameFillColor(10); |
trkpad2[nn]->SetFrameFillColor(10); |
618 |
|
histosig[nn]->SetStats(kFALSE); |
619 |
histosig[nn]->SetLineColor(1); |
histosig[nn]->SetLineColor(1); |
620 |
histosig[nn]->SetFillColor(12); |
histosig[nn]->SetFillColor(12); |
621 |
histosig[nn]->SetLineWidth(1); |
histosig[nn]->SetLineWidth(1); |
625 |
histosig[nn]->GetYaxis()->CenterTitle(); |
histosig[nn]->GetYaxis()->CenterTitle(); |
626 |
histoasig[nn]->SetLineColor(5); |
histoasig[nn]->SetLineColor(5); |
627 |
histoasig[nn]->SetLineWidth(1); |
histoasig[nn]->SetLineWidth(1); |
628 |
if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
// if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ |
629 |
histosig[nn]->Draw("b"); |
histosig[nn]->Draw("b"); |
630 |
if(nn==1){ |
if(nn==0){ |
631 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
632 |
b.DrawBox(2816.,0.,3060.,maxhist); |
b.SetFillColor(107); |
633 |
|
b.SetFillStyle(3945); |
634 |
|
b.DrawBox(768.,0.,2047.,maxhist); |
635 |
|
} |
636 |
|
else if(nn==1){ |
637 |
|
maxhist=histosig[nn]->GetMaximum(); |
638 |
|
b.SetFillColor(6); |
639 |
|
b.SetFillStyle(3945); |
640 |
|
b.DrawBox(2944.,0.,3060.,maxhist); |
641 |
|
|
642 |
|
b.SetFillColor(107); |
643 |
|
b.SetFillStyle(3954); |
644 |
|
b.DrawBox(2816.,0.,2944.,maxhist); |
645 |
|
b.DrawBox(2048.,0.,2176.,maxhist); |
646 |
|
} |
647 |
|
else if(nn==4){ |
648 |
|
maxhist=histosig[nn]->GetMaximum(); |
649 |
|
b.SetFillColor(107); |
650 |
|
b.SetFillStyle(3954); |
651 |
|
b.DrawBox(384.,0.,512.,maxhist); |
652 |
} |
} |
653 |
else if(nn==6){ |
else if(nn==6){ |
654 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
655 |
|
b.SetFillColor(6); |
656 |
|
b.SetFillStyle(3945); |
657 |
b.DrawBox(2560.,0.,2816.,maxhist); |
b.DrawBox(2560.,0.,2816.,maxhist); |
658 |
|
b.DrawBox(1024.,0.,1535.,maxhist); |
659 |
|
|
660 |
|
b.SetFillColor(107); |
661 |
|
b.SetFillStyle(3954); |
662 |
|
b.DrawBox(512.,0.,768.,maxhist); |
663 |
|
b.DrawBox(1536.,0.,1792.,maxhist); |
664 |
|
} |
665 |
|
else if(nn==7){ |
666 |
|
maxhist=histosig[nn]->GetMaximum(); |
667 |
|
b.SetFillColor(107); |
668 |
|
b.SetFillStyle(3954); |
669 |
b.DrawBox(512.,0.,768.,maxhist); |
b.DrawBox(512.,0.,768.,maxhist); |
670 |
b.DrawBox(1024.,0.,1792.,maxhist); |
} |
671 |
|
else if(nn==9){ |
672 |
|
maxhist=histosig[nn]->GetMaximum(); |
673 |
|
b.SetFillColor(107); |
674 |
|
b.SetFillStyle(3954); |
675 |
|
b.DrawBox(256.,0.,384.,maxhist); |
676 |
} |
} |
677 |
else if(nn==11){ |
else if(nn==11){ |
678 |
maxhist=histosig[nn]->GetMaximum(); |
maxhist=histosig[nn]->GetMaximum(); |
679 |
|
b.SetFillColor(6); |
680 |
|
b.SetFillStyle(3945); |
681 |
b.DrawBox(768.,0.,1024.,maxhist); |
b.DrawBox(768.,0.,1024.,maxhist); |
682 |
|
|
683 |
|
b.SetFillColor(107); |
684 |
|
b.SetFillStyle(3954); |
685 |
|
b.DrawBox(0.,0.,512.,maxhist); |
686 |
|
b.DrawBox(1920.,0.,2560.,maxhist); |
687 |
} |
} |
688 |
} |
// } |
689 |
else histosig[nn]->Draw("axis"); |
// else histosig[nn]->Draw("axis"); |
690 |
histoasig[nn]->Draw("same"); |
histoasig[nn]->Draw("same"); |
691 |
li.DrawLine(1024.5,0,1024.5,max); |
li.DrawLine(1024.5,0,1024.5,max); |
692 |
li.DrawLine(2048.5,0,2048.5,max); |
li.DrawLine(2048.5,0,2048.5,max); |
693 |
|
|
694 |
|
for(int ii=0;ii<3;ii++){ |
695 |
|
sig->cd(); |
696 |
|
trkpad3[nn*3+ii]->Draw(); |
697 |
|
trkpad3[nn*3+ii]->cd(); |
698 |
|
trkpad3[nn*3+ii]->SetFillColor(10); |
699 |
|
trkpad3[nn*3+ii]->SetFrameFillColor(10); |
700 |
|
trkpad3[nn*3+ii]->SetLogy(); |
701 |
|
histosiglad[nn][ii]->SetLineColor(1); |
702 |
|
histosiglad[nn][ii]->SetFillColor(1); |
703 |
|
histosiglad[nn][ii]->SetLineWidth(1); |
704 |
|
histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)"); |
705 |
|
histosiglad[nn][ii]->GetXaxis()->CenterTitle(); |
706 |
|
histosiglad[nn][ii]->Draw(""); |
707 |
|
} |
708 |
|
|
709 |
};//end loop on views |
};//end loop on views |
|
|
|
710 |
c1->Update();//draw pads in canvas |
c1->Update();//draw pads in canvas |
711 |
c2->Update();//draw pads in canvas |
c2->Update();//draw pads in canvas |
712 |
|
sig->Update();//draw pads in canvas |
713 |
stringstream nom1; |
|
714 |
stringstream nom2; |
stringstream nom1,nom2,nom3; |
715 |
|
|
716 |
if(!strcmp(outfile.Data(),"ps")||!strcmp(outfile.Data(),"pdf")){ |
if(!strcmp(outfile.Data(),"ps")||!strcmp(outfile.Data(),"pdf")){ |
717 |
nom1.str(""); |
nom1.str(""); |
718 |
nom2.str(""); |
nom2.str(""); |
719 |
|
nom3.str(""); |
720 |
nom1<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps("; |
nom1<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps("; |
721 |
nom2<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps)"; |
nom2<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps"; |
722 |
|
nom3<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps)"; |
723 |
c1->Print(nom1.str().c_str(),"Landscape"); |
c1->Print(nom1.str().c_str(),"Landscape"); |
724 |
c2->Print(nom2.str().c_str(),"Landscape"); |
c2->Print(nom2.str().c_str(),"Landscape"); |
725 |
|
sig->Print(nom3.str().c_str(),"Landscape"); |
726 |
|
|
727 |
if(!strcmp(outfile.Data(),"pdf")){ |
if(!strcmp(outfile.Data(),"pdf")){ |
728 |
stringstream com; |
stringstream com; |
729 |
com<<"ps2pdf13 "<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps "<<out<<ffile<<"_FTrkCalibQlook_EXPERT-pkt"<<i+1<<".pdf"; |
com<<"ps2pdf13 "<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps "<<out<<ffile<<"_FTrkCalibQlook_EXPERT-pkt"<<i+1<<".pdf"; |
739 |
else{ |
else{ |
740 |
nom1.str(""); |
nom1.str(""); |
741 |
nom2.str(""); |
nom2.str(""); |
742 |
|
nom3.str(""); |
743 |
nom1<<out<<ffile<<"_FTrkCalibQLook_EXPERT-ped-pkt"<<i+1<<"."<<outfile.Data(); |
nom1<<out<<ffile<<"_FTrkCalibQLook_EXPERT-ped-pkt"<<i+1<<"."<<outfile.Data(); |
744 |
nom2<<out<<ffile<<"_FTrkCalibQLook_EXPERT-sig-pkt"<<i+1<<"."<<outfile.Data(); |
nom2<<out<<ffile<<"_FTrkCalibQLook_EXPERT-sig-pkt"<<i+1<<"."<<outfile.Data(); |
745 |
|
nom3<<out<<ffile<<"_FTrkCalibQLook_EXPERT-histosig-pkt"<<i+1<<"."<<outfile.Data(); |
746 |
c1->Print(nom1.str().c_str()); |
c1->Print(nom1.str().c_str()); |
747 |
c2->Print(nom2.str().c_str()); |
c2->Print(nom2.str().c_str()); |
748 |
|
sig->Print(nom3.str().c_str()); |
749 |
} |
} |
750 |
|
|
751 |
if(step==1 && i!=maxevent-1 ){ |
if(step==1 && i!=maxevent-1 ){ |
752 |
printf("\n Press 1<enter> to continue, 2<enter> to quit.\n"); |
printf("\n Press 1<enter> to continue, 2<enter> to quit.\n"); |
753 |
cin>>risposta; |
cin>>risposta; |