--- quicklook/tracker/flight/macros/FTrkQLook_EXPERT.cxx 2006/05/10 10:06:55 1.1 +++ quicklook/tracker/flight/macros/FTrkQLook_EXPERT.cxx 2006/05/18 13:06:08 1.2 @@ -137,7 +137,7 @@ const Int_t size=nevent; Int_t dsp=0; Double_t yd[size*12], yc[size*12], x[size], eventint[size]; - 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]; TGraph *comprtime,*temp[5]; TGraph *eventnumb; @@ -563,17 +563,17 @@ oss3<<"DSP h "<GetXaxis()->SetTitle("max signal"); landau1[n]->GetXaxis()->CenterTitle(); - 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); landau2[n]->GetXaxis()->SetTitle("max signal"); landau2[n]->GetXaxis()->CenterTitle(); - 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); landau3[n]->GetXaxis()->SetTitle("max signal"); landau3[n]->GetXaxis()->CenterTitle(); - 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); illuminazione[n]->GetXaxis()->SetTitle("strip with max signal"); illuminazione[n]->GetXaxis()->CenterTitle(); @@ -587,17 +587,24 @@ //********************************************************************************** // Fill Graphs and Histos //********************************************************************************** + + Double_t maxhist=0; + TBox b; + b.SetFillColor(5); for (Int_t i=0; i<12 ; i++){ for (Int_t ev=minevent; evFill(land[12*(ev-minevent)+i][1]); - landau2[i]->Fill(land[12*(ev-minevent)+i][2]); - landau3[i]->Fill(land[12*(ev-minevent)+i][3]); + // 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]); + if( land[12*(ev-minevent)+i][0]!=0) landau1[i]->Fill(land[12*(ev-minevent)+i][0]); + if( land[12*(ev-minevent)+i][1]!=0) landau2[i]->Fill(land[12*(ev-minevent)+i][1]); + if( land[12*(ev-minevent)+i][2]!=0) landau3[i]->Fill(land[12*(ev-minevent)+i][2]); for(Int_t j=0;j<3;j++){ - hi=(j*1024)+illum[12*(ev-minevent)+i][j]; - illuminazione[i]->Fill(hi); + if(land[12*(ev-minevent)+i][j]!=0){ + hi=(j*1024)+illum[12*(ev-minevent)+i][j]; + illuminazione[i]->Fill(hi); + } } } @@ -626,10 +633,33 @@ pad2[i]->Draw(); pad2[i]->cd(); pad2[i]->SetLogy(); - illuminazione [i]->SetStats(kFALSE); - illuminazione[i]->Draw(""); + // illuminazione[i]->SetStats(KFALSE); + if(i==1){ + illuminazione[i]->Draw("ah"); + maxhist= illuminazione[i]->GetMaximum(); + b.DrawBox(2816.,0.,3060.,maxhist); + illuminazione[i]->Draw("axis same"); + illuminazione[i]->Draw("same ]["); + } + else if(i==6){ + illuminazione[i]->Draw("ah"); + maxhist= illuminazione[i]->GetMaximum(); + b.DrawBox(2560.,0.,2816.,maxhist); + b.DrawBox(512.,0.,768.,maxhist); + b.DrawBox(1024.,0.,1792.,maxhist); + illuminazione[i]->Draw("axis same"); + illuminazione[i]->Draw("same]["); + } + else if(i==11){ + illuminazione[i]->Draw("ah"); + maxhist= illuminazione[i]->GetMaximum(); + b.DrawBox(768.,0.,1024.,maxhist); + illuminazione[i]->Draw("axis same"); + illuminazione[i]->Draw("same]["); + } + else illuminazione[i]->Draw(""); - CompTimeCanv->cd(); + CompTimeCanv->cd(); pad3[i]->SetFillColor(10); pad3[i]->SetFrameFillColor(10); pad3[i]->Draw();