/[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.16 by pam-fi, Thu Oct 12 16:58:34 2006 UTC revision 1.28 by pam-fi, Thu Mar 8 01:41:02 2007 UTC
# Line 2  Line 2 
2   * FTrkCalibQLookExpert.cxx   * FTrkCalibQLookExpert.cxx
3   *   *
4   * autor: D.Fedele   * autor: D.Fedele
5   * version v1r14   * version v1r26
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 77  void FTrkCalibQLook_EXPERT(TString file, Line 77  void FTrkCalibQLook_EXPERT(TString 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];    ULong_t FOBT[2];
81        
82    FOBT[0]=0;    FOBT[0]=0;
83    FOBT[1]=0;    FOBT[1]=0;
# Line 87  void FTrkCalibQLook_EXPERT(TString file, Line 87  void FTrkCalibQLook_EXPERT(TString file,
87    ctrk.good0[1]=0;    ctrk.good0[1]=0;
88    for(Int_t i=0;i<12;i++){    for(Int_t i=0;i<12;i++){
89      ctrk.daqmode[i]=0;      ctrk.daqmode[i]=0;
90      ctrk.dspnum[i]=100;      ctrk.dspnum[i]=0;
91      ctrk.calibnum[i]=0;      ctrk.calibnum[i]=0;
92      ctrk.ncalev[i]=0;      ctrk.ncalev[i]=0;
93      ctrk.calfl[i]=0;      ctrk.calfl[i]=0;
# Line 188  void FTrkCalibQLook_EXPERT(TString file, Line 188  void FTrkCalibQLook_EXPERT(TString file,
188            
189    Int_t hcevent=hotr->GetEntries();    Int_t hcevent=hotr->GetEntries();
190    Int_t tcevent=totr->GetEntries();    Int_t tcevent=totr->GetEntries();
191    ULong64_t HOBT[hcevent], TOBT[tcevent];    ULong_t HOBT[hcevent], TOBT[tcevent];
192    for (Int_t i = 0; i < hcevent; i++){    for (Int_t i = 0; i < hcevent; i++){
193      hotr->GetEntry(i);      hotr->GetEntry(i);
194      ph = eh4->GetPscuHeader();      ph = eh4->GetPscuHeader();
# Line 273  void FTrkCalibQLook_EXPERT(TString file, Line 273  void FTrkCalibQLook_EXPERT(TString file,
273    
274      //      //
275      // other variables definitions      // other variables definitions
276    
277      Int_t risposta=0;      Int_t risposta=0;
278      stringstream fromfile;      stringstream fromfile;
279    
# Line 418  void FTrkCalibQLook_EXPERT(TString file, Line 419  void FTrkCalibQLook_EXPERT(TString file,
419      for(Int_t n = 0; n<12; n++){      for(Int_t n = 0; n<12; n++){
420        if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD        if(ctrk.ncalev[n]==0 && ctrk.calfl[n]==0)calok = 1;//GOOD
421            
422        if(ctrk.dspnum[n]>13)        if(ctrk.dspnum[n]==0)
423          continue;          continue;
424    
425        nn=ctrk.dspnum[n]-1;        nn=ctrk.dspnum[n]-1;
# Line 506  void FTrkCalibQLook_EXPERT(TString file, Line 507  void FTrkCalibQLook_EXPERT(TString file,
507        histoaped[nn]->Fill(1025,ctrk.ped2[n]);        histoaped[nn]->Fill(1025,ctrk.ped2[n]);
508        histoaped[nn]->Fill(2049,ctrk.ped3[n]);        histoaped[nn]->Fill(2049,ctrk.ped3[n]);
509        /******************************************************/        /******************************************************/
510                    
511        TLine li;        TLine li,liva1;
512        li.SetLineColor(38);        li.SetLineColor(38);
513        li.SetLineStyle(3);        li.SetLineStyle(4);
514        li.SetLineWidth(2);        li.SetLineWidth(2);
515                liva1.SetLineColor(42);
516        Float_t maxhist=0;        liva1.SetLineStyle(3);
517          liva1.SetLineWidth(1);
518    
519          Float_t maxhist=0,va1x=0;
520        TBox b;        TBox b;
521        /* plot PEDESTAL */        /* plot PEDESTAL */
522        c1->cd();        c1->cd();
# Line 533  void FTrkCalibQLook_EXPERT(TString file, Line 537  void FTrkCalibQLook_EXPERT(TString file,
537        histoaped[nn]->SetLineWidth(1);        histoaped[nn]->SetLineWidth(1);
538  //       if(ctrk.good0[0]==1 && ctrk.good0[1]==1){  //       if(ctrk.good0[0]==1 && ctrk.good0[1]==1){
539          histoped[nn]->Draw("b");          histoped[nn]->Draw("b");
540          if(nn==1){          maxhist=histoped[nn]->GetMaximum();
541            maxhist=histoped[nn]->GetMaximum();          if(nn==0){
542              b.SetFillColor(107);
543              b.SetFillStyle(3945);
544              b.DrawBox(256.,2200.,384.,maxhist);
545              b.DrawBox(768.,2200.,2047.,maxhist);
546            }
547            else if(nn==1){
548            b.SetFillColor(6);            b.SetFillColor(6);
549            b.SetFillStyle(3945);            b.SetFillStyle(3945);
550            b.DrawBox(2944.,700.,3060.,maxhist);            b.DrawBox(2944.,700.,3060.,maxhist);
551    
552            b.SetFillColor(107);            b.SetFillColor(107);
553            b.SetFillStyle(3954);            b.SetFillStyle(3954);
554              b.DrawBox(384.,700.,512.,maxhist);
555              b.DrawBox(2048.,700.,2432.,maxhist);
556            b.DrawBox(2816.,700.,2944.,maxhist);            b.DrawBox(2816.,700.,2944.,maxhist);
557            b.DrawBox(2048.,700.,2176.,maxhist);          }
558            else if(nn==3){
559              b.SetFillColor(107);
560              b.SetFillStyle(3954);
561              b.DrawBox(2816.,700.,3070.,maxhist);
562          }          }
563          else if(nn==4){          else if(nn==4){
           maxhist=histoped[nn]->GetMaximum();  
564            b.SetFillColor(107);            b.SetFillColor(107);
565            b.SetFillStyle(3954);            b.SetFillStyle(3954);
566            b.DrawBox(384.,2200.,512.,maxhist);            b.DrawBox(256.,2200.,512.,maxhist);
567              b.DrawBox(2816.,2200.,3070.,maxhist);
568            }
569            else if(nn==5){
570              b.SetFillColor(107);
571              b.SetFillStyle(3954);
572              b.DrawBox(896.,700.,1024.,maxhist);      
573          }          }
574          else if(nn==6){          else if(nn==6){
           maxhist=histoped[nn]->GetMaximum();  
           b.SetFillColor(6);  
           b.SetFillStyle(3945);  
           b.DrawBox(2560.,2200.,2816.,maxhist);  
           b.DrawBox(1024.,2200.,1535.,maxhist);  
             
575            b.SetFillColor(107);            b.SetFillColor(107);
576            b.SetFillStyle(3954);            b.SetFillStyle(3954);
577            b.DrawBox(512.,2200.,768.,maxhist);            b.DrawBox(512.,2200.,768.,maxhist);
578            b.DrawBox(1536.,2200.,1792.,maxhist);            b.DrawBox(1024.,2200.,1280.,maxhist);
579              b.DrawBox(1280.,2200.,1792.,maxhist);
580              b.DrawBox(2560.,2200.,2816.,maxhist);
581          }          }
582          else if(nn==7){          else if(nn==7){
           maxhist=histoped[nn]->GetMaximum();  
583            b.SetFillColor(107);            b.SetFillColor(107);
584            b.SetFillStyle(3954);            b.SetFillStyle(3954);
585            b.DrawBox(512.,700.,768.,maxhist);            b.DrawBox(512.,700.,768.,maxhist);
586              b.DrawBox(896.,700.,1535.,maxhist);
587              b.DrawBox(2024.,700.,2280.,maxhist);
588            }
589            else if(nn==8){
590              b.SetFillColor(107);
591              b.SetFillStyle(3954);
592              b.DrawBox(512.,2200.,768.,maxhist);
593          }          }
594          else if(nn==9){          else if(nn==9){
           maxhist=histoped[nn]->GetMaximum();  
595            b.SetFillColor(107);            b.SetFillColor(107);
596            b.SetFillStyle(3954);            b.SetFillStyle(3954);
597            b.DrawBox(256.,700.,384.,maxhist);            b.DrawBox(256.,700.,384.,maxhist);
598              b.DrawBox(1280.,700.,1535.,maxhist);
599              b.DrawBox(1792.,700.,1920.,maxhist);
600            }
601            else if(nn==10){
602              b.SetFillColor(107);
603              b.SetFillStyle(3954);
604              b.DrawBox(1024.,2200.,1152.,maxhist);
605              b.DrawBox(2048.,2200.,3070.,maxhist);
606          }          }
607          else if(nn==11){          else if(nn==11){
           maxhist=histoped[nn]->GetMaximum();  
608            b.SetFillColor(6);            b.SetFillColor(6);
609            b.SetFillStyle(3945);            b.SetFillStyle(3945);
610            b.DrawBox(768.,700.,1024.,maxhist);            b.DrawBox(768.,700.,1024.,maxhist);
# Line 589  void FTrkCalibQLook_EXPERT(TString file, Line 618  void FTrkCalibQLook_EXPERT(TString file,
618  //       else histoped[nn]->Draw("axis");  //       else histoped[nn]->Draw("axis");
619        histoaped[nn]->Draw("same");        histoaped[nn]->Draw("same");
620        if((nn+1)%2==1) {        if((nn+1)%2==1) {
621          li.DrawLine(1024.5,2200,1024.5,3200);          for(int va=1; va<24; va++){
622          li.DrawLine(2048.5,2200,2048.5,3200);            va1x=128*va;
623              liva1.DrawLine(va1x,2200.,va1x,3200.);
624            }
625            li.DrawLine(1024.5,2200.,1024.5,3200.);
626            li.DrawLine(2048.5,2200.,2048.5,3200.);
627        }        }
628        if((nn+1)%2==0) {        if((nn+1)%2==0) {
629            for(int va=1; va<24; va++){
630              va1x=128*va;
631              liva1.DrawLine(va1x,700.,va1x,1700.);
632            }
633          li.DrawLine(1024.5,700,1024.5,1700);          li.DrawLine(1024.5,700,1024.5,1700);
634          li.DrawLine(2048.5,700,2048.5,1700);          li.DrawLine(2048.5,700,2048.5,1700);
635        }        }
# Line 620  void FTrkCalibQLook_EXPERT(TString file, Line 657  void FTrkCalibQLook_EXPERT(TString file,
657        histoasig[nn]->SetLineColor(5);        histoasig[nn]->SetLineColor(5);
658        histoasig[nn]->SetLineWidth(1);        histoasig[nn]->SetLineWidth(1);
659  //      if(ctrk.good0[0]==1 && ctrk.good0[1]==1){  //      if(ctrk.good0[0]==1 && ctrk.good0[1]==1){
660          histosig[nn]->Draw("b");          histosig[nn]->Draw("b");
661          if(nn==1){          maxhist=histosig[nn]->GetMaximum();
662            maxhist=histosig[nn]->GetMaximum();          if(nn==0){
663              b.SetFillColor(107);
664              b.SetFillStyle(3945);
665              b.DrawBox(256.,0.,384.,maxhist);
666              b.DrawBox(768.,0.,2047.,maxhist);
667            }
668            else if(nn==1){
669            b.SetFillColor(6);            b.SetFillColor(6);
670            b.SetFillStyle(3945);            b.SetFillStyle(3945);
671            b.DrawBox(2944.,0.,3060.,maxhist);            b.DrawBox(2944.,0.,3070.,maxhist);
672    
673            b.SetFillColor(107);            b.SetFillColor(107);
674            b.SetFillStyle(3954);            b.SetFillStyle(3954);
675              b.DrawBox(384.,0.,512.,maxhist);
676              b.DrawBox(2048.,0.,2432.,maxhist);
677            b.DrawBox(2816.,0.,2944.,maxhist);            b.DrawBox(2816.,0.,2944.,maxhist);
678            b.DrawBox(2048.,0.,2176.,maxhist);          }
679            else if(nn==3){
680              b.SetFillColor(107);
681              b.SetFillStyle(3954);
682              b.DrawBox(2816.,0.,3070.,maxhist);
683          }          }
684          else if(nn==4){          else if(nn==4){
           maxhist=histosig[nn]->GetMaximum();  
685            b.SetFillColor(107);            b.SetFillColor(107);
686            b.SetFillStyle(3954);            b.SetFillStyle(3954);
687            b.DrawBox(384.,0.,512.,maxhist);            b.DrawBox(256.,0.,512.,maxhist);
688              b.DrawBox(2816.,0.,3070.,maxhist);
689            }
690            else if(nn==5){
691              b.SetFillColor(107);
692              b.SetFillStyle(3954);
693              b.DrawBox(896.,0.,1024.,maxhist);      
694          }          }
695          else if(nn==6){          else if(nn==6){
           maxhist=histosig[nn]->GetMaximum();  
           b.SetFillColor(6);  
           b.SetFillStyle(3945);  
           b.DrawBox(2560.,0.,2816.,maxhist);  
           b.DrawBox(1024.,0.,1535.,maxhist);  
             
696            b.SetFillColor(107);            b.SetFillColor(107);
697            b.SetFillStyle(3954);            b.SetFillStyle(3954);
698            b.DrawBox(512.,0.,768.,maxhist);            b.DrawBox(512.,0.,768.,maxhist);
699            b.DrawBox(1536.,0.,1792.,maxhist);            b.DrawBox(1024.,0.,1280.,maxhist);
700              b.DrawBox(1280.,0.,1792.,maxhist);
701              b.DrawBox(2560.,0.,2816.,maxhist);
702          }          }
703          else if(nn==7){          else if(nn==7){
704            maxhist=histosig[nn]->GetMaximum();            b.SetFillColor(107);
705              b.SetFillStyle(3954);
706              b.DrawBox(512.,0.,768.,maxhist);
707              b.DrawBox(896.,0.,1535.,maxhist);
708              b.DrawBox(2024.,0.,2280.,maxhist);
709            }
710            else if(nn==8){
711            b.SetFillColor(107);            b.SetFillColor(107);
712            b.SetFillStyle(3954);            b.SetFillStyle(3954);
713            b.DrawBox(512.,0.,768.,maxhist);            b.DrawBox(512.,0.,768.,maxhist);
714          }          }
715          else if(nn==9){          else if(nn==9){
           maxhist=histosig[nn]->GetMaximum();  
716            b.SetFillColor(107);            b.SetFillColor(107);
717            b.SetFillStyle(3954);            b.SetFillStyle(3954);
718            b.DrawBox(256.,0.,384.,maxhist);            b.DrawBox(256.,0.,384.,maxhist);
719              b.DrawBox(1280.,0.,1535.,maxhist);
720              b.DrawBox(1792.,0.,1920.,maxhist);
721            }
722            else if(nn==10){
723              b.SetFillColor(107);
724              b.SetFillStyle(3954);
725              b.DrawBox(1024.,0.,1152.,maxhist);
726              b.DrawBox(2048.,0.,3070.,maxhist);
727          }          }
728          else if(nn==11){          else if(nn==11){
           maxhist=histosig[nn]->GetMaximum();  
729            b.SetFillColor(6);            b.SetFillColor(6);
730            b.SetFillStyle(3945);            b.SetFillStyle(3945);
731            b.DrawBox(768.,0.,1024.,maxhist);            b.DrawBox(768.,0.,1024.,maxhist);
# Line 676  void FTrkCalibQLook_EXPERT(TString file, Line 738  void FTrkCalibQLook_EXPERT(TString file,
738  //       }  //       }
739  //       else histosig[nn]->Draw("axis");  //       else histosig[nn]->Draw("axis");
740        histoasig[nn]->Draw("same");        histoasig[nn]->Draw("same");
741          for(int va=1; va<24; va++){
742            va1x=128*va;
743            liva1.DrawLine(va1x,0.,va1x,max);
744          }
745        li.DrawLine(1024.5,0,1024.5,max);        li.DrawLine(1024.5,0,1024.5,max);
746        li.DrawLine(2048.5,0,2048.5,max);        li.DrawLine(2048.5,0,2048.5,max);
747          c2->Update();//draw pads in canvas
748                
749        for(int ii=0;ii<3;ii++){        for(int ii=0;ii<3;ii++){
750          sig->cd();          sig->cd();
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.23