/[PAMELA software]/quicklook/tracker/flight/macros/FTrkCalibQLook_EXPERT.cxx
ViewVC logotype

Diff of /quicklook/tracker/flight/macros/FTrkCalibQLook_EXPERT.cxx

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.8 by pam-fi, Thu Jul 13 10:13:37 2006 UTC revision 1.9 by pam-fi, Fri Aug 11 10:24:09 2006 UTC
# Line 2  Line 2 
2   * FTrkCalibQLookExpert.cxx   * FTrkCalibQLookExpert.cxx
3   *   *
4   * autor: D.Fedele   * autor: D.Fedele
5   * version v1r06   * version v1r08
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
# Line 187  void FTrkCalibQLook_EXPERT(TString file, Line 187  void FTrkCalibQLook_EXPERT(TString file,
187    //**********************************************************************    //**********************************************************************
188            
189    Int_t hcevent=hotr->GetEntries();    Int_t hcevent=hotr->GetEntries();
   printf("\n%d\n",hcevent);  
190    ULong64_t HOBT[hcevent], TOBT[hcevent];    ULong64_t HOBT[hcevent], TOBT[hcevent];
191    for (Int_t i = 0; i < hcevent; i++){    for (Int_t i = 0; i < hcevent; i++){
192      totr->GetEntry(i);      totr->GetEntry(i);
# Line 266  void FTrkCalibQLook_EXPERT(TString file, Line 265  void FTrkCalibQLook_EXPERT(TString file,
265      Int_t risposta=0;      Int_t risposta=0;
266      stringstream fromfile;      stringstream fromfile;
267    
     printf("\n\t%d\n",cod->Get(pctp->CalibHeader));  
268      fromfile<<"FTrkCalibQLook_EXPERT      File: "<<ffile<<"        -- CalibHeader OBT= "<<HOBT[(cod->Get(pctp->CalibHeader))-1]<<" -- Calib pkt OBT= "<<FOBT[0]<<" -- CalibTrailer OBT= "<<TOBT[(cod->Get(pctp->CalibHeader))-1]<<" --";      fromfile<<"FTrkCalibQLook_EXPERT      File: "<<ffile<<"        -- CalibHeader OBT= "<<HOBT[(cod->Get(pctp->CalibHeader))-1]<<" -- Calib pkt OBT= "<<FOBT[0]<<" -- CalibTrailer OBT= "<<TOBT[(cod->Get(pctp->CalibHeader))-1]<<" --";
269            
270      gStyle->SetLabelSize(0.08,"x");      gStyle->SetLabelSize(0.07,"x");
271      gStyle->SetLabelSize(0.08,"y");      gStyle->SetLabelSize(0.07,"y");
272      gStyle->SetTitleFillColor(10);      gStyle->SetTitleFillColor(10);
273      gStyle->SetTitleFontSize(0.1);      gStyle->SetTitleFontSize(0.08);
274      gStyle->SetTitleOffset(0.8,"y");      gStyle->SetTitleOffset(0.8,"y");
275      gStyle->SetTitleOffset(1.,"x");      gStyle->SetTitleOffset(0.9,"x");
276      gStyle->SetTitleSize(0.06,"y");      gStyle->SetTitleSize(0.06,"y");
277      gStyle->SetTitleSize(0.06,"x");      gStyle->SetTitleSize(0.06,"x");
278      gStyle->SetOptStat(0);      gStyle->SetOptStat(101110);  
279        gStyle->SetStatX(0.9);  
280        gStyle->SetStatW(0.4);    
281        gStyle->SetStatColor(10);
282        gStyle->SetStatFontSize(0.1);
283            
284      //      //
285      //  draw display area      //  draw display area
286    
287      TLatex *tzz=new TLatex();      TLatex *tzz=new TLatex();
288      tzz->SetTextFont(32);      tzz->SetTextFont(32);
289      tzz->SetTextColor(1);      tzz->SetTextColor(1);
# Line 297  void FTrkCalibQLook_EXPERT(TString file, Line 300  void FTrkCalibQLook_EXPERT(TString file,
300      tzz->DrawLatex(.01,0.98,fromfile.str().c_str());      tzz->DrawLatex(.01,0.98,fromfile.str().c_str());
301      tzz->DrawLatex(.90,0.98,"SIGMA");      tzz->DrawLatex(.90,0.98,"SIGMA");
302    
303      
304        TCanvas *sig=new TCanvas("sig","FTrkCalibQLook_EXPERT_histosig",canvasx,canvasy);
305        sig->SetFillColor(10);
306        tzz->DrawLatex(.01,0.98,fromfile.str().c_str());
307        tzz->DrawLatex(.85,0.97,"Histograms of the sigmas");
308        
309    
310    
311      //  draw pads      //  draw pads
312      TPad *trkpad1[12],*trkpad2[12];        //pad for histos      TPad *trkpad1[12],*trkpad2[12],*trkpad3[36];        //pad for histos
313      TPaveText *trkpadtext[12]; //pad for header      TPaveText *trkpadtext[12]; //pad for header
314      TH1F *histosig[12];           //histos of sigma      TH1F *histosig[12];           //histos of sigma
315      TH1F *histoped[12];          //histos of pedestals        TH1F *histoped[12];          //histos of pedestals  
316      TH1F *histoasig[12];           //histos of sigma      TH1F *histoasig[12];           //histos of sigma
317      TH1F *histoaped[12];          //histos of pedestals      TH1F *histoaped[12];          //histos of pedestals
318    
319        TH1F *histosiglad[12][3];           //histos of sigma
320        stringstream title;
321        stringstream hid;
322            
323      Float_t posy = 0.95;    // up y-coord - top pads      Float_t posy = 0.95;    // up y-coord - top pads
324      Float_t hpad = 0.15;   // pad height      Float_t hpad = 0.15;   // pad height
# Line 312  void FTrkCalibQLook_EXPERT(TString file, Line 327  void FTrkCalibQLook_EXPERT(TString file,
327      Float_t posx0=0;          //       x-coord - column division      Float_t posx0=0;          //       x-coord - column division
328      Float_t wrel = 0.6;     // relative x size of first sub-column      Float_t wrel = 0.6;     // relative x size of first sub-column
329      Float_t marg = 0.004;   // margin among pads      Float_t marg = 0.004;   // margin among pads
     stringstream title;  
     stringstream hid;  
330    
331        
332      for(Int_t n = 0; n<12; n++){      for(Int_t n = 0; n<12; n++){
333        if ( (n+1)%2 ) {        if ( (n+1)%2 ) {
334          if(n>1)posy = posy-(marg*2+hpad);          if(n>1)posy = posy-(marg*2+hpad);
335          posx1 = marg;          posx1 = marg;
336          posx2 = 0.5 - marg;          posx2 = 0.5 - marg;
337          posx0 = 0.5*wrel;          posx0 = 0.5*wrel;
338    
339        } else {        } else {
340          posx1 = posx1 + 0.5;          posx1 = posx1 + 0.5;
341          posx2 = posx2 + 0.5;          posx2 = posx2 + 0.5;
342          posx0 = posx0 + 0.5;          posx0 = posx0 + 0.5;
343    
344        };        };
345        /* -----------> pad for histograms  */        /* -----------> pad for histograms  */
346        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);
# Line 341  void FTrkCalibQLook_EXPERT(TString file, Line 357  void FTrkCalibQLook_EXPERT(TString file,
357        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);
358        hid.str("");        hid.str("");
359        hid<<"hhh"<<n<<"i"<<i;        hid<<"hhh"<<n<<"i"<<i;
       title.str("");  
360        hid.str("");        hid.str("");
361        /* AVERAGE calibration parameters */            /* AVERAGE calibration parameters */    
362        hid<<"ah"<<n<<"i"<<i;        hid<<"ah"<<n<<"i"<<i;
# Line 350  void FTrkCalibQLook_EXPERT(TString file, Line 365  void FTrkCalibQLook_EXPERT(TString file,
365        hid<<"ahh"<<n<<"i"<<i;        hid<<"ahh"<<n<<"i"<<i;
366        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);
367        hid.str("");        hid.str("");
368          for(int ii=0;ii<3;ii++){
369            title.str("");
370            title<<"DSP "<<n+1<<" / Lad "<<ii+1;
371            hid<<"hhhh"<<n<<"i"<<i<<"ii"<<ii;
372            histosiglad[n][ii] = new TH1F(hid.str().c_str(),title.str().c_str(),32,-0.5,30.5);
373            hid.str("");
374          }
375          title.str("");
376      };                     //end loop on views      };                     //end loop on views
377        
378        Float_t tposy = 0.95;    // up y-coord - top pads
379        Float_t thpad = 0.;   // pad height
380        Float_t tposx1=0;          // left  x-coord - pad column
381        Float_t tposx0=0;          //       x-coord - column division
382        Float_t twrel = 0.;     // relative x size of first sub-column
383        Float_t tmarg = 0.002;   // margin among pads
384        thpad = (tposy-tmarg*11)/6;
385        twrel = (1-tmarg*12)/6;
386    
387        for(Int_t n = 0; n<36; n++){
388          if ( (n+1)%6==1 ) {
389            if(n>1) tposy = tposy-(tmarg*2+thpad);
390            tposx1 = tmarg;
391            tposx0 = tposx1 + twrel;
392          } else {
393            tposx1 = tposx0 + 2*tmarg;
394            tposx0 = tposx1 + twrel;
395          }
396          trkpad3[n]= new TPad("pad3"," ",tposx1,tposy-thpad,tposx0,tposy,18,0,0);
397        }
398        
399      
400      stringstream message;      stringstream message;
401      
402      //--------------------------------      //--------------------------------
403      //CHECK CALIBRATION procedure      //CHECK CALIBRATION procedure
404      //--------------------------------      //--------------------------------
# Line 362  void FTrkCalibQLook_EXPERT(TString file, Line 406  void FTrkCalibQLook_EXPERT(TString file,
406      Int_t calok = 0;//BAD        Int_t calok = 0;//BAD  
407      for(Int_t n = 0; n<12; n++){      for(Int_t n = 0; n<12; n++){
408        if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD        if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD
409          
410                
411        nn=ctrk.dspnum[n]-1;        nn=ctrk.dspnum[n]-1;
412        /*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*        /*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
# Line 438  void FTrkCalibQLook_EXPERT(TString file, Line 482  void FTrkCalibQLook_EXPERT(TString file,
482        for(Int_t j = 0; j < 3072; j++){        for(Int_t j = 0; j < 3072; j++){
483          histosig[nn]->Fill((Float_t)j,ctrk.dspsig[nn][j]);          histosig[nn]->Fill((Float_t)j,ctrk.dspsig[nn][j]);
484          histoped[nn]->Fill((Float_t)j,ctrk.dspped[nn][j]);          histoped[nn]->Fill((Float_t)j,ctrk.dspped[nn][j]);
485            if(j<1024) histosiglad[nn][0]->Fill(ctrk.dspsig[nn][j]);
486            if(j>=1024 && j<2048) histosiglad[nn][1]->Fill(ctrk.dspsig[nn][j]);
487            if(j>=2048 && j<3072) histosiglad[nn][2]->Fill(ctrk.dspsig[nn][j]);
488        };        };
489        histoasig[nn]->Fill(1,ctrk.sig1[nn]);        histoasig[nn]->Fill(1,ctrk.sig1[nn]);
490        histoasig[nn]->Fill(1025,ctrk.sig2[nn]);        histoasig[nn]->Fill(1025,ctrk.sig2[nn]);
# Line 454  void FTrkCalibQLook_EXPERT(TString file, Line 501  void FTrkCalibQLook_EXPERT(TString file,
501                
502        Float_t maxhist=0;        Float_t maxhist=0;
503        TBox b;        TBox b;
       b.SetFillColor(6);  
       b.SetFillStyle(3945);  
504        /* plot PEDESTAL */        /* plot PEDESTAL */
505        c1->cd();        c1->cd();
506        trkpadtext[nn]->Draw();        trkpadtext[nn]->Draw();
# Line 463  void FTrkCalibQLook_EXPERT(TString file, Line 508  void FTrkCalibQLook_EXPERT(TString file,
508        trkpad1[nn]->cd();        trkpad1[nn]->cd();
509        trkpad1[nn]->SetFillColor(10);        trkpad1[nn]->SetFillColor(10);
510        trkpad1[nn]->SetFrameFillColor(10);        trkpad1[nn]->SetFrameFillColor(10);
511          histoped[nn]->SetStats(kFALSE);
512        histoped[nn]->SetLineColor(1);        histoped[nn]->SetLineColor(1);
513        histoped[nn]->SetFillColor(12);        histoped[nn]->SetFillColor(12);
514        histoped[nn]->SetLineWidth(1);        histoped[nn]->SetLineWidth(1);
# Line 476  void FTrkCalibQLook_EXPERT(TString file, Line 522  void FTrkCalibQLook_EXPERT(TString file,
522          histoped[nn]->Draw("b");          histoped[nn]->Draw("b");
523          if(nn==1){          if(nn==1){
524            maxhist=histoped[nn]->GetMaximum();            maxhist=histoped[nn]->GetMaximum();
525            b.DrawBox(2816.,700.,3060.,maxhist);            b.SetFillColor(6);
526              b.SetFillStyle(3945);
527              b.DrawBox(2944.,700.,3060.,maxhist);
528    
529              b.SetFillColor(107);
530              b.SetFillStyle(3954);
531              b.DrawBox(2816.,700.,2944.,maxhist);
532              b.DrawBox(2048.,700.,2176.,maxhist);
533          }          }
534          else if(nn==6){          else if(nn==6){
535            maxhist=histoped[nn]->GetMaximum();            maxhist=histoped[nn]->GetMaximum();
536              b.SetFillColor(6);
537              b.SetFillStyle(3945);
538            b.DrawBox(2560.,2200.,2816.,maxhist);            b.DrawBox(2560.,2200.,2816.,maxhist);
           b.DrawBox(512.,2200.,768.,maxhist);  
539            b.DrawBox(1024.,2200.,1792.,maxhist);            b.DrawBox(1024.,2200.,1792.,maxhist);
540              
541              b.SetFillColor(107);
542              b.SetFillStyle(3954);
543              b.DrawBox(512.,2200.,768.,maxhist);
544            }
545            else if(nn==7){
546              maxhist=histoped[nn]->GetMaximum();
547              b.SetFillColor(107);
548              b.SetFillStyle(3954);
549              b.DrawBox(512.,700.,768.,maxhist);
550          }          }
551          else if(nn==11){          else if(nn==11){
552            maxhist=histoped[nn]->GetMaximum();            maxhist=histoped[nn]->GetMaximum();
553              b.SetFillColor(6);
554              b.SetFillStyle(3945);
555            b.DrawBox(768.,700.,1024.,maxhist);            b.DrawBox(768.,700.,1024.,maxhist);
556              
557              b.SetFillColor(107);
558              b.SetFillStyle(3954);
559              b.DrawBox(0.,700.,512.,maxhist);
560              b.DrawBox(1920.,700.,2048.,maxhist);
561          }          }
562        }        }
563        else histoped[nn]->Draw("axis");        else histoped[nn]->Draw("axis");
# Line 511  void FTrkCalibQLook_EXPERT(TString file, Line 582  void FTrkCalibQLook_EXPERT(TString file,
582        trkpad2[nn]->cd();        trkpad2[nn]->cd();
583        trkpad2[nn]->SetFillColor(10);        trkpad2[nn]->SetFillColor(10);
584        trkpad2[nn]->SetFrameFillColor(10);        trkpad2[nn]->SetFrameFillColor(10);
585          histosig[nn]->SetStats(kFALSE);
586        histosig[nn]->SetLineColor(1);        histosig[nn]->SetLineColor(1);
587        histosig[nn]->SetFillColor(12);        histosig[nn]->SetFillColor(12);
588        histosig[nn]->SetLineWidth(1);        histosig[nn]->SetLineWidth(1);
# Line 524  void FTrkCalibQLook_EXPERT(TString file, Line 596  void FTrkCalibQLook_EXPERT(TString file,
596          histosig[nn]->Draw("b");          histosig[nn]->Draw("b");
597          if(nn==1){          if(nn==1){
598            maxhist=histosig[nn]->GetMaximum();            maxhist=histosig[nn]->GetMaximum();
599            b.DrawBox(2816.,0.,3060.,maxhist);            b.SetFillColor(6);
600              b.SetFillStyle(3945);
601              b.DrawBox(2944.,0.,3060.,maxhist);
602    
603              b.SetFillColor(107);
604              b.SetFillStyle(3954);
605              b.DrawBox(2816.,0.,2944.,maxhist);
606              b.DrawBox(2048.,0.,2176.,maxhist);
607          }          }
608          else if(nn==6){          else if(nn==6){
609            maxhist=histosig[nn]->GetMaximum();            maxhist=histosig[nn]->GetMaximum();
610              b.SetFillColor(6);
611              b.SetFillStyle(3945);
612            b.DrawBox(2560.,0.,2816.,maxhist);            b.DrawBox(2560.,0.,2816.,maxhist);
           b.DrawBox(512.,0.,768.,maxhist);  
613            b.DrawBox(1024.,0.,1792.,maxhist);            b.DrawBox(1024.,0.,1792.,maxhist);
614              
615              b.SetFillColor(107);
616              b.SetFillStyle(3954);
617              b.DrawBox(512.,0.,768.,maxhist);
618            }
619            else if(nn==7){
620              maxhist=histosig[nn]->GetMaximum();
621              b.SetFillColor(107);
622              b.SetFillStyle(3954);
623              b.DrawBox(512.,0.,768.,maxhist);
624          }          }
625          else if(nn==11){          else if(nn==11){
626            maxhist=histosig[nn]->GetMaximum();            maxhist=histosig[nn]->GetMaximum();
627              b.SetFillColor(6);
628              b.SetFillStyle(3945);
629            b.DrawBox(768.,0.,1024.,maxhist);            b.DrawBox(768.,0.,1024.,maxhist);
630              
631              b.SetFillColor(107);
632              b.SetFillStyle(3954);
633              b.DrawBox(0.,0.,512.,maxhist);
634              b.DrawBox(1920.,0.,2048.,maxhist);
635          }          }
636        }        }
637        else histosig[nn]->Draw("axis");        else histosig[nn]->Draw("axis");
# Line 542  void FTrkCalibQLook_EXPERT(TString file, Line 639  void FTrkCalibQLook_EXPERT(TString file,
639        li.DrawLine(1024.5,0,1024.5,max);        li.DrawLine(1024.5,0,1024.5,max);
640        li.DrawLine(2048.5,0,2048.5,max);        li.DrawLine(2048.5,0,2048.5,max);
641                
642          for(int ii=0;ii<3;ii++){
643            sig->cd();
644      };//end loop on views          trkpad3[nn*3+ii]->Draw();
645            trkpad3[nn*3+ii]->cd();
646            trkpad3[nn*3+ii]->SetFillColor(10);
647            trkpad3[nn*3+ii]->SetFrameFillColor(10);
648            trkpad3[nn*3+ii]->SetLogy();
649            histosiglad[nn][ii]->SetLineColor(1);
650            histosiglad[nn][ii]->SetFillColor(1);
651            histosiglad[nn][ii]->SetLineWidth(1);
652            histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)");
653            histosiglad[nn][ii]->GetXaxis()->CenterTitle();
654            histosiglad[nn][ii]->Draw("");  
655          }
656            
657        };//end loop on views
658      c1->Update();//draw pads in canvas      c1->Update();//draw pads in canvas
659      c2->Update();//draw pads in canvas      c2->Update();//draw pads in canvas
660        sig->Update();//draw pads in canvas
661      stringstream nom1;      
662      stringstream nom2;      stringstream nom1,nom2,nom3;
663          
664      if(!strcmp(outfile.Data(),"ps")||!strcmp(outfile.Data(),"pdf")){      if(!strcmp(outfile.Data(),"ps")||!strcmp(outfile.Data(),"pdf")){
665        nom1.str("");        nom1.str("");
666        nom2.str("");        nom2.str("");
667          nom3.str("");
668        nom1<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps(";        nom1<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps(";
669        nom2<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps)";        nom2<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps";
670          nom3<<out<<ffile<<"_FTrkCalibQLook_EXPERT-pkt"<<i+1<<".ps)";
671        c1->Print(nom1.str().c_str(),"Landscape");        c1->Print(nom1.str().c_str(),"Landscape");
672        c2->Print(nom2.str().c_str(),"Landscape");        c2->Print(nom2.str().c_str(),"Landscape");
673          sig->Print(nom3.str().c_str(),"Landscape");
674          
675        if(!strcmp(outfile.Data(),"pdf")){        if(!strcmp(outfile.Data(),"pdf")){
676          stringstream com;          stringstream com;
677          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";
# Line 574  void FTrkCalibQLook_EXPERT(TString file, Line 687  void FTrkCalibQLook_EXPERT(TString file,
687      else{      else{
688        nom1.str("");        nom1.str("");
689        nom2.str("");        nom2.str("");
690          nom3.str("");
691        nom1<<out<<ffile<<"_FTrkCalibQLook_EXPERT-ped-pkt"<<i+1<<"."<<outfile.Data();        nom1<<out<<ffile<<"_FTrkCalibQLook_EXPERT-ped-pkt"<<i+1<<"."<<outfile.Data();
692        nom2<<out<<ffile<<"_FTrkCalibQLook_EXPERT-sig-pkt"<<i+1<<"."<<outfile.Data();        nom2<<out<<ffile<<"_FTrkCalibQLook_EXPERT-sig-pkt"<<i+1<<"."<<outfile.Data();
693          nom3<<out<<ffile<<"_FTrkCalibQLook_EXPERT-histosig-pkt"<<i+1<<"."<<outfile.Data();
694        c1->Print(nom1.str().c_str());        c1->Print(nom1.str().c_str());
695        c2->Print(nom2.str().c_str());        c2->Print(nom2.str().c_str());
696          sig->Print(nom3.str().c_str());
697      }      }
698              
699      if(step==1 && i!=maxevent-1 ){        if(step==1 && i!=maxevent-1 ){  
700        printf("\n Press 1<enter> to continue, 2<enter> to quit.\n");        printf("\n Press 1<enter> to continue, 2<enter> to quit.\n");
701        cin>>risposta;        cin>>risposta;
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.23