/[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.15 by pam-fi, Mon Oct 2 13:05:39 2006 UTC revision 1.27 by pam-fi, Mon Feb 26 16:41:05 2007 UTC
# Line 2  Line 2 
2   * FTrkCalibQLookExpert.cxx   * FTrkCalibQLookExpert.cxx
3   *   *
4   * autor: D.Fedele   * autor: D.Fedele
5   * version v1r13   * version v1r25
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 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]==0)
423            continue;
424    
425        nn=ctrk.dspnum[n]-1;        nn=ctrk.dspnum[n]-1;
426        /*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*        /*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
427         *         *
# Line 491  void FTrkCalibQLook_EXPERT(TString file, Line 494  void FTrkCalibQLook_EXPERT(TString file,
494        /******************************************************/        /******************************************************/
495        /* fill histos */        /* fill histos */
496        for(Int_t j = 0; j < 3072; j++){        for(Int_t j = 0; j < 3072; j++){
497          histosig[nn]->Fill((Float_t)j,ctrk.dspsig[nn][j]);          histosig[nn]->Fill(j,ctrk.dspsig[n][j]);
498          histoped[nn]->Fill((Float_t)j,ctrk.dspped[nn][j]);          histoped[nn]->Fill(j,ctrk.dspped[n][j]);
499          if(j<1024) histosiglad[nn][0]->Fill(ctrk.dspsig[nn][j]);          if(j<1024) histosiglad[nn][0]->Fill(ctrk.dspsig[n][j]);
500          if(j>=1024 && j<2048) histosiglad[nn][1]->Fill(ctrk.dspsig[nn][j]);          if(j>=1024 && j<2048) histosiglad[nn][1]->Fill(ctrk.dspsig[n][j]);
501          if(j>=2048 && j<3072) histosiglad[nn][2]->Fill(ctrk.dspsig[nn][j]);          if(j>=2048 && j<3072) histosiglad[nn][2]->Fill(ctrk.dspsig[n][j]);
502        };        };
503        histoasig[nn]->Fill(1,ctrk.sig1[nn]);        histoasig[nn]->Fill(1,ctrk.sig1[n]);
504        histoasig[nn]->Fill(1025,ctrk.sig2[nn]);        histoasig[nn]->Fill(1025,ctrk.sig2[n]);
505        histoasig[nn]->Fill(2049,ctrk.sig3[nn]);        histoasig[nn]->Fill(2049,ctrk.sig3[n]);
506        histoaped[nn]->Fill(1,ctrk.ped1[nn]);        histoaped[nn]->Fill(1,ctrk.ped1[nn]);
507        histoaped[nn]->Fill(1025,ctrk.ped2[nn]);        histoaped[nn]->Fill(1025,ctrk.ped2[n]);
508        histoaped[nn]->Fill(2049,ctrk.ped3[nn]);        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 529  void FTrkCalibQLook_EXPERT(TString file, Line 535  void FTrkCalibQLook_EXPERT(TString file,
535        if((nn+1)%2==0)  histoped[nn]->GetYaxis()->SetRangeUser(700,1700);        if((nn+1)%2==0)  histoped[nn]->GetYaxis()->SetRangeUser(700,1700);
536        histoaped[nn]->SetLineColor(5);        histoaped[nn]->SetLineColor(5);
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();  
575            b.SetFillColor(6);            b.SetFillColor(6);
576            b.SetFillStyle(3945);            b.SetFillStyle(3945);
577            b.DrawBox(2560.,2200.,2816.,maxhist);            b.DrawBox(1024.,2200.,1280.,maxhist);
           b.DrawBox(1024.,2200.,1535.,maxhist);  
578                        
579            b.SetFillColor(107);            b.SetFillColor(107);
580            b.SetFillStyle(3954);            b.SetFillStyle(3954);
581            b.DrawBox(512.,2200.,768.,maxhist);            b.DrawBox(512.,2200.,768.,maxhist);
582            b.DrawBox(1536.,2200.,1792.,maxhist);            b.DrawBox(1280.,2200.,1792.,maxhist);
583              b.DrawBox(2560.,2200.,2816.,maxhist);
584          }          }
585          else if(nn==7){          else if(nn==7){
           maxhist=histoped[nn]->GetMaximum();  
586            b.SetFillColor(107);            b.SetFillColor(107);
587            b.SetFillStyle(3954);            b.SetFillStyle(3954);
588            b.DrawBox(512.,700.,768.,maxhist);            b.DrawBox(512.,700.,768.,maxhist);
589              b.DrawBox(896.,700.,1535.,maxhist);
590            }
591            else if(nn==8){
592              b.SetFillColor(107);
593              b.SetFillStyle(3954);
594              b.DrawBox(512.,2200.,768.,maxhist);
595          }          }
596          else if(nn==9){          else if(nn==9){
           maxhist=histoped[nn]->GetMaximum();  
597            b.SetFillColor(107);            b.SetFillColor(107);
598            b.SetFillStyle(3954);            b.SetFillStyle(3954);
599            b.DrawBox(256.,700.,384.,maxhist);            b.DrawBox(256.,700.,384.,maxhist);
600              b.DrawBox(1280.,700.,1535.,maxhist);
601              b.DrawBox(1792.,700.,1920.,maxhist);
602            }
603            else if(nn==10){
604              b.SetFillColor(107);
605              b.SetFillStyle(3954);
606              b.DrawBox(1024.,2200.,1152.,maxhist);
607              b.DrawBox(2048.,2200.,3070.,maxhist);
608          }          }
609          else if(nn==11){          else if(nn==11){
           maxhist=histoped[nn]->GetMaximum();  
610            b.SetFillColor(6);            b.SetFillColor(6);
611            b.SetFillStyle(3945);            b.SetFillStyle(3945);
612            b.DrawBox(768.,700.,1024.,maxhist);            b.DrawBox(768.,700.,1024.,maxhist);
# Line 581  void FTrkCalibQLook_EXPERT(TString file, Line 614  void FTrkCalibQLook_EXPERT(TString file,
614            b.SetFillColor(107);            b.SetFillColor(107);
615            b.SetFillStyle(3954);            b.SetFillStyle(3954);
616            b.DrawBox(0.,700.,512.,maxhist);            b.DrawBox(0.,700.,512.,maxhist);
617            b.DrawBox(1920.,700.,2048.,maxhist);            b.DrawBox(1920.,700.,2560.,maxhist);
           b.DrawBox(2176.,700.,2304.,maxhist);  
618          }          }
619        }    //     }
620        else histoped[nn]->Draw("axis");  //       else histoped[nn]->Draw("axis");
621        histoaped[nn]->Draw("same");        histoaped[nn]->Draw("same");
622        if((nn+1)%2==1) {        if((nn+1)%2==1) {
623          li.DrawLine(1024.5,2200,1024.5,3200);          for(int va=1; va<24; va++){
624          li.DrawLine(2048.5,2200,2048.5,3200);            va1x=128*va;
625              liva1.DrawLine(va1x,2200.,va1x,3200.);
626            }
627            li.DrawLine(1024.5,2200.,1024.5,3200.);
628            li.DrawLine(2048.5,2200.,2048.5,3200.);
629        }        }
630        if((nn+1)%2==0) {        if((nn+1)%2==0) {
631            for(int va=1; va<24; va++){
632              va1x=128*va;
633              liva1.DrawLine(va1x,700.,va1x,1700.);
634            }
635          li.DrawLine(1024.5,700,1024.5,1700);          li.DrawLine(1024.5,700,1024.5,1700);
636          li.DrawLine(2048.5,700,2048.5,1700);          li.DrawLine(2048.5,700,2048.5,1700);
637        }        }
638          c1->Update();//draw pads in canvas
639                
640    
641                
# Line 617  void FTrkCalibQLook_EXPERT(TString file, Line 658  void FTrkCalibQLook_EXPERT(TString file,
658        histosig[nn]->GetYaxis()->CenterTitle();        histosig[nn]->GetYaxis()->CenterTitle();
659        histoasig[nn]->SetLineColor(5);        histoasig[nn]->SetLineColor(5);
660        histoasig[nn]->SetLineWidth(1);        histoasig[nn]->SetLineWidth(1);
661        if(ctrk.good0[0]==1 && ctrk.good0[1]==1){  //      if(ctrk.good0[0]==1 && ctrk.good0[1]==1){
662          histosig[nn]->Draw("b");          histosig[nn]->Draw("b");
663          if(nn==1){          maxhist=histosig[nn]->GetMaximum();
664            maxhist=histosig[nn]->GetMaximum();          if(nn==0){
665              b.SetFillColor(107);
666              b.SetFillStyle(3945);
667              b.DrawBox(256.,0.,384.,maxhist);
668              b.DrawBox(768.,0.,2047.,maxhist);
669            }
670            else if(nn==1){
671            b.SetFillColor(6);            b.SetFillColor(6);
672            b.SetFillStyle(3945);            b.SetFillStyle(3945);
673            b.DrawBox(2944.,0.,3060.,maxhist);            b.DrawBox(2944.,0.,3070.,maxhist);
674    
675            b.SetFillColor(107);            b.SetFillColor(107);
676            b.SetFillStyle(3954);            b.SetFillStyle(3954);
677              b.DrawBox(384.,0.,512.,maxhist);
678              b.DrawBox(2048.,0.,2432.,maxhist);
679            b.DrawBox(2816.,0.,2944.,maxhist);            b.DrawBox(2816.,0.,2944.,maxhist);
680            b.DrawBox(2048.,0.,2176.,maxhist);          }
681            else if(nn==3){
682              b.SetFillColor(107);
683              b.SetFillStyle(3954);
684              b.DrawBox(2816.,0.,3070.,maxhist);
685          }          }
686          else if(nn==4){          else if(nn==4){
           maxhist=histosig[nn]->GetMaximum();  
687            b.SetFillColor(107);            b.SetFillColor(107);
688            b.SetFillStyle(3954);            b.SetFillStyle(3954);
689            b.DrawBox(384.,0.,512.,maxhist);            b.DrawBox(256.,0.,512.,maxhist);
690              b.DrawBox(2816.,0.,3070.,maxhist);
691            }
692            else if(nn==5){
693              b.SetFillColor(107);
694              b.SetFillStyle(3954);
695              b.DrawBox(896.,0.,1024.,maxhist);      
696          }          }
697          else if(nn==6){          else if(nn==6){
           maxhist=histosig[nn]->GetMaximum();  
698            b.SetFillColor(6);            b.SetFillColor(6);
699            b.SetFillStyle(3945);            b.SetFillStyle(3945);
700            b.DrawBox(2560.,0.,2816.,maxhist);            b.DrawBox(1024.,0.,1280.,maxhist);
           b.DrawBox(1024.,0.,1535.,maxhist);  
701                        
702            b.SetFillColor(107);            b.SetFillColor(107);
703            b.SetFillStyle(3954);            b.SetFillStyle(3954);
704            b.DrawBox(512.,0.,768.,maxhist);            b.DrawBox(512.,0.,768.,maxhist);
705            b.DrawBox(1536.,0.,1792.,maxhist);            b.DrawBox(1280.,0.,1792.,maxhist);
706              b.DrawBox(2560.,0.,2816.,maxhist);
707          }          }
708          else if(nn==7){          else if(nn==7){
709            maxhist=histosig[nn]->GetMaximum();            b.SetFillColor(107);
710              b.SetFillStyle(3954);
711              b.DrawBox(512.,0.,768.,maxhist);
712              b.DrawBox(896.,0.,1535.,maxhist);
713            }
714            else if(nn==8){
715            b.SetFillColor(107);            b.SetFillColor(107);
716            b.SetFillStyle(3954);            b.SetFillStyle(3954);
717            b.DrawBox(512.,0.,768.,maxhist);            b.DrawBox(512.,0.,768.,maxhist);
718          }          }
719          else if(nn==9){          else if(nn==9){
           maxhist=histosig[nn]->GetMaximum();  
720            b.SetFillColor(107);            b.SetFillColor(107);
721            b.SetFillStyle(3954);            b.SetFillStyle(3954);
722            b.DrawBox(256.,0.,384.,maxhist);            b.DrawBox(256.,0.,384.,maxhist);
723              b.DrawBox(1280.,0.,1535.,maxhist);
724              b.DrawBox(1792.,0.,1920.,maxhist);
725            }
726            else if(nn==10){
727              b.SetFillColor(107);
728              b.SetFillStyle(3954);
729              b.DrawBox(1024.,0.,1152.,maxhist);
730              b.DrawBox(2048.,0.,3070.,maxhist);
731          }          }
732          else if(nn==11){          else if(nn==11){
           maxhist=histosig[nn]->GetMaximum();  
733            b.SetFillColor(6);            b.SetFillColor(6);
734            b.SetFillStyle(3945);            b.SetFillStyle(3945);
735            b.DrawBox(768.,0.,1024.,maxhist);            b.DrawBox(768.,0.,1024.,maxhist);
# Line 669  void FTrkCalibQLook_EXPERT(TString file, Line 737  void FTrkCalibQLook_EXPERT(TString file,
737            b.SetFillColor(107);            b.SetFillColor(107);
738            b.SetFillStyle(3954);            b.SetFillStyle(3954);
739            b.DrawBox(0.,0.,512.,maxhist);            b.DrawBox(0.,0.,512.,maxhist);
740            b.DrawBox(1920.,0.,2048.,maxhist);            b.DrawBox(1920.,0.,2560.,maxhist);
           b.DrawBox(2176.,0.,2304.,maxhist);  
741          }          }
742        }  //       }
743        else histosig[nn]->Draw("axis");  //       else histosig[nn]->Draw("axis");
744        histoasig[nn]->Draw("same");        histoasig[nn]->Draw("same");
745          for(int va=1; va<24; va++){
746            va1x=128*va;
747            liva1.DrawLine(va1x,0.,va1x,max);
748          }
749        li.DrawLine(1024.5,0,1024.5,max);        li.DrawLine(1024.5,0,1024.5,max);
750        li.DrawLine(2048.5,0,2048.5,max);        li.DrawLine(2048.5,0,2048.5,max);
751          c2->Update();//draw pads in canvas
752                
753        for(int ii=0;ii<3;ii++){        for(int ii=0;ii<3;ii++){
754          sig->cd();          sig->cd();
# Line 691  void FTrkCalibQLook_EXPERT(TString file, Line 763  void FTrkCalibQLook_EXPERT(TString file,
763          histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)");          histosiglad[nn][ii]->GetXaxis()->SetTitle("SIG (ADC channels)");
764          histosiglad[nn][ii]->GetXaxis()->CenterTitle();          histosiglad[nn][ii]->GetXaxis()->CenterTitle();
765          histosiglad[nn][ii]->Draw("");            histosiglad[nn][ii]->Draw("");  
766        }        }
767              
768      };//end loop on views      };//end loop on views
769      c1->Update();//draw pads in canvas      c1->Update();//draw pads in canvas
770      c2->Update();//draw pads in canvas      c2->Update();//draw pads in canvas
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.23