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

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

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

revision 1.1 by pam-fi, Wed May 10 10:06:55 2006 UTC revision 1.2 by pam-fi, Thu May 18 13:06:08 2006 UTC
# Line 137  void FTrkQLook_EXPERT(TString file,Int_t Line 137  void FTrkQLook_EXPERT(TString file,Int_t
137    const Int_t size=nevent;    const Int_t size=nevent;
138    Int_t dsp=0;    Int_t dsp=0;
139    Double_t yd[size*12], yc[size*12], x[size], eventint[size];    Double_t yd[size*12], yc[size*12], x[size], eventint[size];
140    Double_t yyc[size], hi=0, land[size*12][3], illum[size*12][3] ;    Double_t yyc[size], hi=0, land[size*12][3], illum[size*12][3];
141    
142    TGraph *comprtime,*temp[5];    TGraph *comprtime,*temp[5];
143    TGraph *eventnumb;    TGraph *eventnumb;
# Line 563  void FTrkQLook_EXPERT(TString file,Int_t Line 563  void FTrkQLook_EXPERT(TString file,Int_t
563      oss3<<"DSP h "<<n+1;      oss3<<"DSP h "<<n+1;
564      oss4<<"DSP   h"<<n+1;      oss4<<"DSP   h"<<n+1;
565            
566      landau1[n]=new TH1F(oss1.str().c_str(),oss1.str().c_str(),1200,0.5,1200.5);      landau1[n]=new TH1F(oss1.str().c_str(),oss1.str().c_str(),1201,-0.5,1200.5);
567      landau1[n]->GetXaxis()->SetTitle("max signal");      landau1[n]->GetXaxis()->SetTitle("max signal");
568      landau1[n]->GetXaxis()->CenterTitle();      landau1[n]->GetXaxis()->CenterTitle();
569      landau2[n]=new TH1F(oss2.str().c_str(),oss1.str().c_str(),1200,0.5,1200.5);      landau2[n]=new TH1F(oss2.str().c_str(),oss1.str().c_str(),1201,-0.5,1200.5);
570      landau2[n]->GetXaxis()->SetTitle("max signal");      landau2[n]->GetXaxis()->SetTitle("max signal");
571      landau2[n]->GetXaxis()->CenterTitle();      landau2[n]->GetXaxis()->CenterTitle();
572      landau3[n]=new TH1F(oss3.str().c_str(),oss1.str().c_str(),1200,0.5,1200.5);      landau3[n]=new TH1F(oss3.str().c_str(),oss1.str().c_str(),1201,-0.5,1200.5);
573      landau3[n]->GetXaxis()->SetTitle("max signal");      landau3[n]->GetXaxis()->SetTitle("max signal");
574      landau3[n]->GetXaxis()->CenterTitle();      landau3[n]->GetXaxis()->CenterTitle();
575    
576      illuminazione[n]=new TH1F(oss4.str().c_str(),oss1.str().c_str(),3072,0.5,3072.5);      illuminazione[n]=new TH1F(oss4.str().c_str(),oss1.str().c_str(),3073,-0.5,3072.5);
577      illuminazione[n]->GetXaxis()->SetTitle("strip with max signal");      illuminazione[n]->GetXaxis()->SetTitle("strip with max signal");
578      illuminazione[n]->GetXaxis()->CenterTitle();      illuminazione[n]->GetXaxis()->CenterTitle();
579    
# Line 587  void FTrkQLook_EXPERT(TString file,Int_t Line 587  void FTrkQLook_EXPERT(TString file,Int_t
587    //**********************************************************************************    //**********************************************************************************
588    // Fill Graphs and Histos    // Fill Graphs and Histos
589    //**********************************************************************************    //**********************************************************************************
590      
591      Double_t maxhist=0;
592      TBox b;
593      b.SetFillColor(5);
594    
595    for (Int_t i=0; i<12 ; i++){    for (Int_t i=0; i<12 ; i++){
596    
597      for (Int_t ev=minevent; ev<maxevent; ev++){      for (Int_t ev=minevent; ev<maxevent; ev++){
598        yyc[(ev-minevent)]=yc[12*(ev-minevent)+i];        yyc[(ev-minevent)]=yc[12*(ev-minevent)+i];
599        landau1[i]->Fill(land[12*(ev-minevent)+i][1]);        // if( land[12*(ev-minevent)+i][0]!=0 && illum[12*(ev-minevent)+i][0]>510 && illum[12*(ev-minevent)+i][0]<515)  landau1[i]->Fill(land[12*(ev-minevent)+i][0]);
600        landau2[i]->Fill(land[12*(ev-minevent)+i][2]);        if( land[12*(ev-minevent)+i][0]!=0)  landau1[i]->Fill(land[12*(ev-minevent)+i][0]);
601        landau3[i]->Fill(land[12*(ev-minevent)+i][3]);        if( land[12*(ev-minevent)+i][1]!=0)  landau2[i]->Fill(land[12*(ev-minevent)+i][1]);
602          if( land[12*(ev-minevent)+i][2]!=0)  landau3[i]->Fill(land[12*(ev-minevent)+i][2]);
603        for(Int_t j=0;j<3;j++){        for(Int_t j=0;j<3;j++){
604          hi=(j*1024)+illum[12*(ev-minevent)+i][j];          if(land[12*(ev-minevent)+i][j]!=0){
605          illuminazione[i]->Fill(hi);            hi=(j*1024)+illum[12*(ev-minevent)+i][j];
606              illuminazione[i]->Fill(hi);
607            }
608        }        }
609      }      }
610            
# Line 626  void FTrkQLook_EXPERT(TString file,Int_t Line 633  void FTrkQLook_EXPERT(TString file,Int_t
633      pad2[i]->Draw();      pad2[i]->Draw();
634      pad2[i]->cd();      pad2[i]->cd();
635      pad2[i]->SetLogy();      pad2[i]->SetLogy();
636      illuminazione [i]->SetStats(kFALSE);      // illuminazione[i]->SetStats(KFALSE);
637      illuminazione[i]->Draw("");     if(i==1){
638          illuminazione[i]->Draw("ah");
639          maxhist= illuminazione[i]->GetMaximum();
640          b.DrawBox(2816.,0.,3060.,maxhist);
641          illuminazione[i]->Draw("axis same");
642          illuminazione[i]->Draw("same ][");
643        }
644        else if(i==6){
645          illuminazione[i]->Draw("ah");
646          maxhist= illuminazione[i]->GetMaximum();
647          b.DrawBox(2560.,0.,2816.,maxhist);
648          b.DrawBox(512.,0.,768.,maxhist);
649          b.DrawBox(1024.,0.,1792.,maxhist);
650          illuminazione[i]->Draw("axis same");
651          illuminazione[i]->Draw("same][");
652        }
653        else if(i==11){
654          illuminazione[i]->Draw("ah");
655          maxhist= illuminazione[i]->GetMaximum();
656          b.DrawBox(768.,0.,1024.,maxhist);
657          illuminazione[i]->Draw("axis same");
658          illuminazione[i]->Draw("same][");
659        }
660        else illuminazione[i]->Draw("");
661    
662      CompTimeCanv->cd();     CompTimeCanv->cd();
663      pad3[i]->SetFillColor(10);      pad3[i]->SetFillColor(10);
664      pad3[i]->SetFrameFillColor(10);      pad3[i]->SetFrameFillColor(10);
665      pad3[i]->Draw();      pad3[i]->Draw();
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.23