/[PAMELA software]/quicklook/QLflightTmtc_Header/HeaderScan.cpp
ViewVC logotype

Diff of /quicklook/QLflightTmtc_Header/HeaderScan.cpp

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

revision 1.2 by pam-rm2, Thu Jun 22 12:16:18 2006 UTC revision 1.3 by pam-rm2, Mon Jun 26 09:42:17 2006 UTC
# Line 85  void HeaderScan(TString base, TString ou Line 85  void HeaderScan(TString base, TString ou
85    
86    //------------------- Variables initilization -------------------------//    //------------------- Variables initilization -------------------------//
87    Long64_t nevents=0;    Long64_t nevents=0;
88    ULong_t lastime, firstime,obt;    ULong_t lastime, firstime,obt, offset, lastoffset, firstoffset;
89    double obmin=0.;    double obmin=0.;
90    double obmax=0.;    double obmax=0.;
   Float_t maxoffset, minoffset;  
91    stringstream oss, oss1, oss2, oss3, noentries;    stringstream oss, oss1, oss2, oss3, noentries;
92    //------- load root file --------------    //------- load root file --------------
93    TFile *file = new TFile(base.Data());    TFile *file = new TFile(base.Data());
# Line 149  void HeaderScan(TString base, TString ou Line 148  void HeaderScan(TString base, TString ou
148      headBr->GetEntry(i);      headBr->GetEntry(i);
149      ph = eh->GetPscuHeader();      ph = eh->GetPscuHeader();
150      obt = ph->GetOrbitalTime();      obt = ph->GetOrbitalTime();
151        offset=ph->GetFileOffset();
152      if(obt <= firstime) firstime=obt;      if(obt <= firstime) firstime=obt;
153      if(obt >= lastime) lastime=obt;      if(obt >= lastime) lastime=obt;
154      }      }
155    
156    //*************************** Histograms ************************************************************//    //*************************** Histograms ************************************************************//
157    //------------------------ First histogram -----------------------------------//    //------------------------ First histogram -----------------------------------//
158    obmin=firstime;    obmin=firstime;
# Line 174  void HeaderScan(TString base, TString ou Line 175  void HeaderScan(TString base, TString ou
175    
176    //**************************************************************************************************//    //**************************************************************************************************//
177    //------- fill histograms ---------//    //------- fill histograms ---------//
178      Int_t l=0;
179    for (Int_t i = 0; i < nevents; i++){    for (Int_t i = 0; i < nevents; i++){
180      headBr->GetEntry(i);      headBr->GetEntry(i);
181      ph = eh->GetPscuHeader();        ph = eh->GetPscuHeader();  
# Line 237  void HeaderScan(TString base, TString ou Line 238  void HeaderScan(TString base, TString ou
238    oss << outDir.Data() << filename.Data();    oss << outDir.Data() << filename.Data();
239    oss  << "_HeaderScan_1." << format.Data();    oss  << "_HeaderScan_1." << format.Data();
240    finalCanv1->SaveAs(oss.str().c_str());    finalCanv1->SaveAs(oss.str().c_str());
241    
242    //------------- Graph PscuCounter vs FileOffset ---------------------------------    //------------- Graph PscuCounter vs FileOffset ---------------------------------
243    const Int_t size = 100000;    const Int_t size = 100000;
244    Int_t ev=0;    Int_t ev=0;
# Line 244  void HeaderScan(TString base, TString ou Line 246  void HeaderScan(TString base, TString ou
246    Int_t m=0;    Int_t m=0;
247    Int_t pag= (nevents/100000)+1;    Int_t pag= (nevents/100000)+1;
248    
249      TMultiGraph *mg1 = new TMultiGraph();
250      oss1.str("");
251      oss1 << filename.Data() <<": PscuCounter vs FileOffset.";
252      mg1->SetTitle(oss1.str().c_str());
253    
254    TCanvas *finalCanv2 = new TCanvas("Header_2", base, 1280, 1024);    TCanvas *finalCanv2 = new TCanvas("Header_2", base, 1280, 1024);
255    finalCanv2->SetFillColor(10);    finalCanv2->SetFillColor(10);
256    
# Line 258  void HeaderScan(TString base, TString ou Line 265  void HeaderScan(TString base, TString ou
265      m=l;      m=l;
266     }       }  
267    
    oss1.str("");  
    oss1 << filename.Data() <<": PscuCounter vs FileOffset. Canvas n."<< (interval+1);  
268     TGraph *graph= new TGraph(m, FileOffset, PscuCounter);     TGraph *graph= new TGraph(m, FileOffset, PscuCounter);
   
    finalCanv2->cd();  
    graph->SetTitle(oss1.str().c_str());  
    graph->GetXaxis()->SetTitle("File Offset");  
    graph->GetXaxis()->CenterTitle();  
    graph->GetXaxis()->SetTitleSize(0.04);  
    graph->GetXaxis()->SetLabelSize(0.02);  
    graph->GetYaxis()->SetTitle("Counter");  
    graph->GetYaxis()->CenterTitle();  
    graph->GetYaxis()->SetTitleSize(0.04);  
    graph->GetYaxis()->SetLabelSize(0.02);  
269     graph->SetMarkerColor(3);     graph->SetMarkerColor(3);
270     graph->SetMarkerSize(.5);     graph->SetMarkerSize(.5);
271     graph->SetMarkerStyle(21);     graph->SetMarkerStyle(21);
272     graph->Draw("AP");     mg1->Add(graph);
273    
274       ev=ev+100000;
275      }
276       finalCanv2->cd();
277       mg1->Draw("AP");
278       mg1->SetTitle(oss1.str().c_str());
279       mg1->GetXaxis()->SetTitle("File Offset");
280       mg1->GetXaxis()->CenterTitle();
281       mg1->GetXaxis()->SetTitleSize(0.04);
282       mg1->GetXaxis()->SetLabelSize(0.02);
283       mg1->GetYaxis()->SetTitle("Counter");
284       mg1->GetYaxis()->CenterTitle();
285       mg1->GetYaxis()->SetTitleSize(0.04);
286       mg1->GetYaxis()->SetLabelSize(0.02);
287    
288    
289       finalCanv2->Update();
290    
291     oss1.str("");     oss1.str("");
292     oss1 << outDir.Data() << filename.Data();     oss1 << outDir.Data() << filename.Data();
293     oss1 << "_HeaderScan_2_"<<(interval+1)<<"." << format.Data();     oss1 << "_HeaderScan_2"<<"." << format.Data();
294    
295     finalCanv2->SaveAs(oss1.str().c_str());     finalCanv2->SaveAs(oss1.str().c_str());
296     finalCanv2->Clear();    
297       file->Close();
    ev=ev+100000;  
   }  
   
   file->Close();  
298        
299  }  }
300    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23