--- quicklook/tracker/flight/macros/FTrkCalibQLook_BASIC.cxx 2006/06/29 11:11:59 1.4 +++ quicklook/tracker/flight/macros/FTrkCalibQLook_BASIC.cxx 2007/03/08 01:41:02 1.10 @@ -2,7 +2,7 @@ * FTrkCalibQLook_BASIC.cxx * * autor: D.Fedele - * version v1r04 + * version v1r11 * Parameters: * file - the data file to analyze * fromevent - first event to analyze @@ -54,7 +54,6 @@ { // // obtain information about the data file and select the output dir - const string filepath=file.Data(); Int_t dwpos = file.Last('/'); Int_t dwpos1 = file.Last('.'); TString base,ffile ; @@ -75,7 +74,7 @@ Int_t nevents=0; Int_t minevent = 0; Int_t maxevent = 0; - ULong64_t OBT[2]; + ULong_t OBT[2]; OBT[0]=0; OBT[1]=0; @@ -164,7 +163,7 @@ } else if (toevent > nevents) { maxevent = nevents; } else { - maxevent = fromevent; + maxevent = nevents; }; }; @@ -193,7 +192,7 @@ Int_t canvasx=900; Int_t canvasy=1200; - Int_t ndsp =0,alarm=0; + Int_t alarm=0; Float_t pedav[12][12],pedavtemp[12][12],sigav[12][12],sigavtemp[12][12]; Int_t flpedav[12][12],flsigav[12][12]; Float_t siglimsup[12][12],sigliminf[12][12],pedlimsup[12][12],pedliminf[12][12]; @@ -273,53 +272,61 @@ for (Int_t m = 0; m < 6; m++){ ph1 = eh1->GetPscuHeader(); OBT[0]= ph1->GetOrbitalTime(); - ctrk.daqmode[trk1->DSPnumber[m]-1]=trk1->DAQmode[m]; - ctrk.dspnum[trk1->DSPnumber[m]-1]=trk1->DSPnumber[m]; - ctrk.calibnum[trk1->DSPnumber[m]-1]=trk1->calibnumber[m]; - ctrk.ncalev[trk1->DSPnumber[m]-1]=trk1->ncalib_event[m]; - ctrk.ped1[trk1->DSPnumber[m]-1]=trk1->ped_l1[m]; - ctrk.ped2[trk1->DSPnumber[m]-1]=trk1->ped_l2[m]; - ctrk.ped3[trk1->DSPnumber[m]-1]=trk1->ped_l3[m]; - ctrk.sig1[trk1->DSPnumber[m]-1]=trk1->sig_l1[m]; - ctrk.sig2[trk1->DSPnumber[m]-1]=trk1->sig_l2[m]; - ctrk.sig3[trk1->DSPnumber[m]-1]=trk1->sig_l3[m]; - ctrk.nbad1[trk1->DSPnumber[m]-1]=trk1->nbad_l1[m]; - ctrk.nbad2[trk1->DSPnumber[m]-1]=trk1->nbad_l2[m]; - ctrk.nbad3[trk1->DSPnumber[m]-1]=trk1->nbad_l3[m]; - ctrk.calfl[trk1->DSPnumber[m]-1]=trk1->cal_flag[m]; - ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk1->crc_cal[m][0]; - ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk1->crc_cal[m][1]; - ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk1->crc_cal[m][2]; - ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk1->crc_hcal[m]; - for (Int_t j = 0; j < 3072; j++){ - ctrk.dspped[trk1->DSPnumber[m]-1][j]=trk1->DSPped_par[m][j]; - ctrk.dspsig[trk1->DSPnumber[m]-1][j]=trk1->DSPsig_par[m][j]; - ctrk.dspbad[trk1->DSPnumber[m]-1][j]=trk1->DSPbad_par[m][j]; + if(trk1->DSPnumber[m]>0){ + if(trk1->DSPnumber[m]<13){ + ctrk.daqmode[trk1->DSPnumber[m]-1]=trk1->DAQmode[m]; + ctrk.dspnum[trk1->DSPnumber[m]-1]=trk1->DSPnumber[m]; + ctrk.calibnum[trk1->DSPnumber[m]-1]=trk1->calibnumber[m]; + ctrk.ncalev[trk1->DSPnumber[m]-1]=trk1->ncalib_event[m]; + ctrk.ped1[trk1->DSPnumber[m]-1]=trk1->ped_l1[m]; + ctrk.ped2[trk1->DSPnumber[m]-1]=trk1->ped_l2[m]; + ctrk.ped3[trk1->DSPnumber[m]-1]=trk1->ped_l3[m]; + ctrk.sig1[trk1->DSPnumber[m]-1]=trk1->sig_l1[m]; + ctrk.sig2[trk1->DSPnumber[m]-1]=trk1->sig_l2[m]; + ctrk.sig3[trk1->DSPnumber[m]-1]=trk1->sig_l3[m]; + ctrk.nbad1[trk1->DSPnumber[m]-1]=trk1->nbad_l1[m]; + ctrk.nbad2[trk1->DSPnumber[m]-1]=trk1->nbad_l2[m]; + ctrk.nbad3[trk1->DSPnumber[m]-1]=trk1->nbad_l3[m]; + ctrk.calfl[trk1->DSPnumber[m]-1]=trk1->cal_flag[m]; + ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk1->crc_cal[m][0]; + ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk1->crc_cal[m][1]; + ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk1->crc_cal[m][2]; + ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk1->crc_hcal[m]; + for (Int_t j = 0; j < 3072; j++){ + ctrk.dspped[trk1->DSPnumber[m]-1][j]=trk1->DSPped_par[m][j]; + ctrk.dspsig[trk1->DSPnumber[m]-1][j]=trk1->DSPsig_par[m][j]; + ctrk.dspbad[trk1->DSPnumber[m]-1][j]=trk1->DSPbad_par[m][j]; + } + } } ph2 = eh2->GetPscuHeader(); OBT[1]= ph2->GetOrbitalTime(); - ctrk.daqmode[trk2->DSPnumber[m]-1]=trk2->DAQmode[m]; - ctrk.dspnum[trk2->DSPnumber[m]-1]=trk2->DSPnumber[m]; - ctrk.calibnum[trk2->DSPnumber[m]-1]=trk2->calibnumber[m]; - ctrk.ncalev[trk2->DSPnumber[m]-1]=trk2->ncalib_event[m]; - ctrk.ped1[trk2->DSPnumber[m]-1]=trk2->ped_l1[m]; - ctrk.ped2[trk2->DSPnumber[m]-1]=trk2->ped_l2[m]; - ctrk.ped3[trk2->DSPnumber[m]-1]=trk2->ped_l3[m]; - ctrk.sig1[trk2->DSPnumber[m]-1]=trk2->sig_l1[m]; - ctrk.sig2[trk2->DSPnumber[m]-1]=trk2->sig_l2[m]; - ctrk.sig3[trk2->DSPnumber[m]-1]=trk2->sig_l3[m]; - ctrk.nbad1[trk2->DSPnumber[m]-1]=trk2->nbad_l1[m]; - ctrk.nbad2[trk2->DSPnumber[m]-1]=trk2->nbad_l2[m]; - ctrk.nbad3[trk2->DSPnumber[m]-1]=trk2->nbad_l3[m]; - ctrk.calfl[trk2->DSPnumber[m]-1]=trk2->cal_flag[m]; - ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk2->crc_cal[m][0]; - ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk2->crc_cal[m][1]; - ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk2->crc_cal[m][2]; - ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk2->crc_hcal[m]; - for (Int_t j = 0; j < 3072; j++){ - ctrk.dspped[trk2->DSPnumber[m]-1][j]=trk2->DSPped_par[m][j]; - ctrk.dspsig[trk2->DSPnumber[m]-1][j]=trk2->DSPsig_par[m][j]; - ctrk.dspbad[trk2->DSPnumber[m]-1][j]=trk2->DSPbad_par[m][j]; + if(trk2->DSPnumber[m]>0){ + if(trk2->DSPnumber[m]<13){ + ctrk.daqmode[trk2->DSPnumber[m]-1]=trk2->DAQmode[m]; + ctrk.dspnum[trk2->DSPnumber[m]-1]=trk2->DSPnumber[m]; + ctrk.calibnum[trk2->DSPnumber[m]-1]=trk2->calibnumber[m]; + ctrk.ncalev[trk2->DSPnumber[m]-1]=trk2->ncalib_event[m]; + ctrk.ped1[trk2->DSPnumber[m]-1]=trk2->ped_l1[m]; + ctrk.ped2[trk2->DSPnumber[m]-1]=trk2->ped_l2[m]; + ctrk.ped3[trk2->DSPnumber[m]-1]=trk2->ped_l3[m]; + ctrk.sig1[trk2->DSPnumber[m]-1]=trk2->sig_l1[m]; + ctrk.sig2[trk2->DSPnumber[m]-1]=trk2->sig_l2[m]; + ctrk.sig3[trk2->DSPnumber[m]-1]=trk2->sig_l3[m]; + ctrk.nbad1[trk2->DSPnumber[m]-1]=trk2->nbad_l1[m]; + ctrk.nbad2[trk2->DSPnumber[m]-1]=trk2->nbad_l2[m]; + ctrk.nbad3[trk2->DSPnumber[m]-1]=trk2->nbad_l3[m]; + ctrk.calfl[trk2->DSPnumber[m]-1]=trk2->cal_flag[m]; + ctrk.crc_c[trk1->DSPnumber[m]-1][0]=trk2->crc_cal[m][0]; + ctrk.crc_c[trk1->DSPnumber[m]-1][1]=trk2->crc_cal[m][1]; + ctrk.crc_c[trk1->DSPnumber[m]-1][2]=trk2->crc_cal[m][2]; + ctrk.crc_hc[trk1->DSPnumber[m]-1]=trk2->crc_hcal[m]; + for (Int_t j = 0; j < 3072; j++){ + ctrk.dspped[trk2->DSPnumber[m]-1][j]=trk2->DSPped_par[m][j]; + ctrk.dspsig[trk2->DSPnumber[m]-1][j]=trk2->DSPsig_par[m][j]; + ctrk.dspbad[trk2->DSPnumber[m]-1][j]=trk2->DSPbad_par[m][j]; + } + } } } @@ -334,14 +341,15 @@ } } - Int_t nn,ok=0; + Int_t nn=0,ok=0; // // write warning if it occur for(Int_t n = 0; n<12; n++){ - ndsp = ctrk.dspnum[n]; - nn = ndsp-1; + if(ctrk.dspnum[n]==0) + continue; + nn = ctrk.dspnum[n]-1; for(Int_t iii=0;iii<3;iii++){ if(ctrk.crc_c[nn][iii]!=0){ @@ -408,7 +416,7 @@ space[wc]-=spacep; } - if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ + // if(ctrk.good0[0]==1 && ctrk.good0[1]==1){ // // evaluate the mean value of the sigma and pedestal for(Int_t j = 0; j < 3072; j++){ @@ -469,9 +477,10 @@ } } } - } + // } } - if(ok==0 && ctrk.good0[0]==1 && ctrk.good0[1]==1){ + if(ok==0){ +// if(ok==0 && ctrk.good0[0]==1 && ctrk.good0[1]==1){ if(space[wc]<=10){ wc+=1; flcanvas+=1;