/[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.12 by pam-rm2, Tue Oct 24 08:23:17 2006 UTC
# Line 1  Line 1 
1  /**  /**
2  * Header Scan  * Header_Scan
3  * Author  Nagni  * Author  Nagni
4  * version 1.0  * version 1.0
5  *  *
# Line 54  Line 54 
54  *             Fixed bugs: for a large namber of events is not possible to have vectors, so I have subsituted graphs with histograms  *             Fixed bugs: for a large namber of events is not possible to have vectors, so I have subsituted graphs with histograms
55  *                         or divided the graphs in more than one canvas.  *                         or divided the graphs in more than one canvas.
56  *  *
57    * Version 1.8
58    * Date 8 August 2006 - Malvezzi
59    *
60    * Description: changed the scale in the second and third graph of the first canvas; added a pad of text in the second canvas
61    *
62  */  */
63    
64    
# Line 65  Line 70 
70  #include "TMultiGraph.h"  #include "TMultiGraph.h"
71  #include <sstream>  #include <sstream>
72  #include <iostream>  #include <iostream>
73    #include <stdio.h>
74    #include <string.h>
75  #include "TString.h"  #include "TString.h"
76  #include "TStyle.h"  #include "TStyle.h"
77  #include "TFile.h"  #include "TFile.h"
# Line 77  Line 84 
84  #include "TGaxis.h"  #include "TGaxis.h"
85  #include "EventHeader.h"  #include "EventHeader.h"
86  #include "PscuHeader.h"  #include "PscuHeader.h"
87    #include "RunHeaderEvent.h"
88    #include "TPaveText.h"
89    
90  using namespace std;  using namespace std;
91    
92  void HeaderScan(TString base, TString outDir, TString format)  void HeaderScan(TString base, TString outDir, TString format){  
 {  
93    
94    //------------------- Variables initilization -------------------------//    //------------------- Variables initilization -------------------------//
95    Long64_t nevents=0;    Long64_t nevents=0, runnevents;
96    ULong_t lastime, firstime,obt;    ULong_t lastime, firstime, primotempo, ultimotempo, primoffset=500000000, ultimoffset;
97    double obmin=0.;    double obmin=0.;
98    double obmax=0.;    double obmax=0.;
99    Float_t maxoffset, minoffset;    stringstream oss, oss1, oss2, oss3, noentries, stringa;
100    stringstream oss, oss1, oss2, oss3, noentries;  
101    //------- load root file --------------    //-------------- Load root file, tree and branch  -------------------//
102    TFile *file = new TFile(base.Data());    TFile *file = new TFile(base.Data());
103    if (!file){    if (!file){
104      printf("No such file in the directory has been found");      printf("No such file in the directory has been found");
105      return;      return;
106    }    }
107      
108    TTree *PhysicsTr = (TTree*)file->Get("Physics");    TTree *PhysicsTr = (TTree*)file->Get("Physics");
109    TBranch *headBr = PhysicsTr->GetBranch("Header");    TBranch *headBr = PhysicsTr->GetBranch("Header");
110    
111    pamela::EventHeader  *eh   = 0;    pamela::EventHeader  *eh   = new pamela::EventHeader;
112    pamela::PscuHeader   *ph   = 0;    pamela::PscuHeader   *ph   = new pamela::PscuHeader;
113    
114    PhysicsTr->SetBranchAddress("Header", &eh);    PhysicsTr->SetBranchAddress("Header", &eh);
115      nevents = PhysicsTr->GetEntries();
116      const Int_t size = nevents;
117    
118    nevents = headBr->GetEntries();    TTree *RunHeadTr = (TTree*)file->Get("RunHeader");  ///run header tree
119      pamela::EventHeader *eH= new pamela::EventHeader;
120      pamela::RunHeaderEvent *reh=new pamela::RunHeaderEvent;
121      
122      RunHeadTr->SetBranchAddress("Header",&eH);
123      RunHeadTr->SetBranchAddress("RunHeader",&reh);
124      runnevents = RunHeadTr->GetEntries();
125    
126    TString filename = ((TObjString*)base.Tokenize('/')->Last())->GetString();    TString filename = ((TObjString*)base.Tokenize('/')->Last())->GetString();
127    filename = ((TObjString*)filename.Tokenize('.')->First())->GetString();    filename = ((TObjString*)filename.Tokenize('.')->First())->GetString();
# Line 114  void HeaderScan(TString base, TString ou Line 130  void HeaderScan(TString base, TString ou
130    if (nevents<=0) {    if (nevents<=0) {
131      printf("nevents = %i \n", nevents);      printf("nevents = %i \n", nevents);
132      printf(" \n");      printf(" \n");
   
133      TCanvas *canv = new TCanvas("No entries", "No entries ", 400, 200);      TCanvas *canv = new TCanvas("No entries", "No entries ", 400, 200);
134      canv->SetFillColor(10);      canv->SetFillColor(10);
135      canv->cd();      canv->cd();
   
136      TLatex *l = new TLatex();      TLatex *l = new TLatex();
137      l->SetTextAlign(12);      l->SetTextAlign(12);
138      l->SetTextSize(0.15);      l->SetTextSize(0.15);
# Line 129  void HeaderScan(TString base, TString ou Line 143  void HeaderScan(TString base, TString ou
143      noentries.str("");      noentries.str("");
144      noentries<< "No Physics entries for this files";      noentries<< "No Physics entries for this files";
145      l->DrawLatex(0.05, 0.5, noentries.str().c_str());      l->DrawLatex(0.05, 0.5, noentries.str().c_str());
   
146      if (outDir == "./") {      if (outDir == "./") {
147        oss.str("");        oss.str("");
148        oss << filename.Data() << "_HeaderScan_QL." << format.Data();        oss << filename.Data() << "_HeaderScan_QL." << format.Data();
149      } else {      } else {
150        oss.str("");        oss.str("");
151        oss << outDir.Data() << filename.Data() << "_HeaderScan_QL." << format.Data();        oss << outDir.Data() << filename.Data() << "_HeaderScan_QL." << format.Data();
152      }      }  
     
153      canv->Update();      canv->Update();
154      canv->SaveAs(oss.str().c_str());      canv->SaveAs(oss.str().c_str());
   
155      return;      return;
156    }    }
   
157    //-------------- to know the max and min OBT ----------------------------//    //-------------- to know the max and min OBT ----------------------------//
158    for (Int_t i = 0; i < nevents; i++){     headBr->GetEntry(0);
159      headBr->GetEntry(i);     ph = eh->GetPscuHeader();
160      ph = eh->GetPscuHeader();     firstime = ph->GetOrbitalTime();
161      obt = ph->GetOrbitalTime();     headBr->GetEntry(nevents);
162      if(obt <= firstime) firstime=obt;     ph = eh->GetPscuHeader();
163      if(obt >= lastime) lastime=obt;     lastime = ph->GetOrbitalTime();
164      }     int i =0;
165    //*************************** Histograms ************************************************************//     while(i < nevents){
166         headBr->GetEntry(i);
167         ph = eh->GetPscuHeader();
168         if((ph->GetOrbitalTime()) <= firstime) firstime=ph->GetOrbitalTime();
169         if((ph->GetOrbitalTime()) >= lastime) lastime=ph->GetOrbitalTime();
170         i++;
171         }
172      
173    //------------------------ First histogram -----------------------------------//    //------------------------ First histogram -----------------------------------//
174    obmin=firstime;    obmin=firstime;
175    obmax=lastime;      obmax=lastime;
   oss1.str("");  
   oss1 <<"File name= "<< filename.Data() <<". Start time= " << obmin << ", end time= "<< obmax <<" ms"<<". Physics Packet per minute";  
176    Int_t  nbin   = (lastime-firstime)/60000;    Int_t  nbin   = (lastime-firstime)/60000;
177    TH1F *h1 = new TH1F ("histo1", oss1.str().c_str(), nbin, obmin, obmax);    TH1F *h1 = new TH1F ("histo1", "" , nbin, obmin, obmax);
   
   //------------------------ Second histogram -----------------------------------//  
   oss2.str("");  
   oss2 << " Lenght of Physic packets;";  
   Int_t nint = (lastime-firstime);  
   TH1F *PcktLenght = new TH1F ("histo2", oss2.str().c_str(), nint, obmin, obmax);  
178    
179   //------------------------ Third histogram -----------------------------------//    //------------------- fill vectors and histogram -----------------------------//
180    oss3.str("");    Double_t *PscuCounter = new Double_t[size];
181    oss3 << " Physics Counter";    Double_t *PacketLenght = new Double_t[size];
182    TH1F *PcktCounter = new TH1F ("histo3", oss3.str().c_str(), nint, obmin, obmax);    Double_t *OBTime = new Double_t[size];
183      Double_t *Eventsperminute= new Double_t[nbin];
184    //**************************************************************************************************//    Double_t *Minute= new Double_t[nbin];
185    //------- fill histograms ---------//    Int_t max=0;
186      for (Int_t k = 0; k < nevents; k++){
187    for (Int_t i = 0; i < nevents; i++){      headBr->GetEntry(k);
     headBr->GetEntry(i);  
188      ph = eh->GetPscuHeader();        ph = eh->GetPscuHeader();  
189      h1->Fill(ph->GetOrbitalTime());        h1->Fill(ph->GetOrbitalTime());
190      PcktLenght->Fill(ph->GetOrbitalTime(),ph->GetPacketLenght());        PscuCounter[k]= ph->GetCounter();
191      PcktCounter->Fill(ph->GetOrbitalTime(),ph->GetCounter());      PacketLenght[k]=ph->GetPacketLenght();
192        OBTime[k]=ph->GetOrbitalTime();
193      }  
194      int l=0;
195      while(l<nbin){
196        Eventsperminute[l]=h1->GetBinContent(l);
197        Minute[l]=firstime+l*60000;
198        if(h1->GetBinContent(l) >= max)max =(Int_t)h1->GetBinContent(l);
199        l++;
200    }    }
201    
202    //****************************** Canvases *******************************//    //----------- Graph and MultiGraph -----------------------------------------------//
203    //TGaxis::SetMaxDigits(4);    TMultiGraph *rate = new TMultiGraph();
204    //------------------- First Canvas --------------------------------//    TMultiGraph *packetLength = new TMultiGraph();
205    TCanvas *finalCanv1 = new TCanvas("Header_1", base, 1280, 1024);    TMultiGraph *packeCounter = new TMultiGraph();
206    finalCanv1->Divide(1,3);  
207    finalCanv1->SetFillColor(10);    oss1.str("");
208      oss1 << "Physics Packet per minute.  Start time = " << obmin << ", End time = "<< obmax <<" ms";
209    finalCanv1->cd(1);    TGraph *rate1= new TGraph(nbin, (const Double_t*)Minute, (const Double_t*)Eventsperminute);
210    h1->SetStats(kFALSE);    rate1->SetMarkerColor(kBlack);
211    h1->GetXaxis()->SetTitle("OBT   (ms)");    rate1->SetMarkerSize(.1);
212    h1->GetXaxis()->CenterTitle();    rate1->SetMarkerStyle(21);
213    h1->GetXaxis()->SetLabelSize(0.04);    rate->Add(rate1);
214    h1->GetYaxis()->SetTitle("Number of events ");  
215    h1->GetYaxis()->CenterTitle();    TGraph *packetLength1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PacketLenght);
216    h1->GetYaxis()->SetTitleSize(0.06);    oss2.str("");
217    h1->GetYaxis()->SetTitleOffset(0.8);    oss2 <<"Lenght of Physic packets";
218    h1->Draw();    packetLength1->SetMarkerColor(2);
219      packetLength1->SetMarkerSize(.3);
220      packetLength1->SetMarkerStyle(21);
221      packetLength->Add(packetLength1);
222    
223      TGraph *packeCounter1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PscuCounter);
224      oss3.str("");
225      oss3 <<"Physics Counter vs. OBT";
226      packeCounter1->SetMarkerColor(4);
227      packeCounter1->SetMarkerSize(.2);
228      packeCounter1->SetMarkerStyle(21);
229      packeCounter->Add(packeCounter1);
230    
231      //------------ Create and Draw Canvas ---------------------//
232      TCanvas *finalCanv = new TCanvas("Header", base, 1200, 1600);
233      finalCanv->Divide(1,6);
234      finalCanv->SetFillColor(10);
235    
236      TPad *all2= new TPad ("","", 0, 0, 1, 1);
237      all2->SetFillColor(10);
238      TPad *all3= new TPad ("","", 0, 0, 1, 1);
239      all3->SetFillColor(10);
240      TPad *all4= new TPad ("","", 0, 0, 1, 1);
241      all4->SetFillColor(10);
242      TPad *all= new TPad ("","", 0, 0, 1, 1);
243      all->SetFillColor(10);
244      TPad *all1= new TPad ("","", 0, 0, 1, 1);
245      all1->SetFillColor(10);
246    
247      TLine li;
248      li.SetLineStyle(4);
249      li.SetLineWidth(2);
250    
251      //----------- First PAD -------------------------------//
252      finalCanv->cd(1);
253      all2->Draw();
254      all2->cd();
255      rate->Draw("ALP");
256      rate->GetXaxis()->SetTitle("OBT   (ms)");
257      rate->GetXaxis()->SetTitleSize(0.05);
258      rate->GetXaxis()->CenterTitle();
259      rate->GetXaxis()->SetLabelSize(0.05);
260      rate->GetYaxis()->SetTitle("Number of events ");
261      rate->GetYaxis()->CenterTitle();
262      rate->GetYaxis()->SetLabelSize(0.05);
263      rate->GetYaxis()->SetTitleSize(0.06);
264      rate->GetYaxis()->SetTitleOffset(0.6);
265      for (Int_t l = 0; l < runnevents; l++){
266        RunHeadTr->GetEntry(l);
267        ph = eH->GetPscuHeader();
268        int ws= reh->RM_ACQ_SETTING_MODE;
269        int id = ph->GetPacketId1();
270        Int_t obt = ph->GetOrbitalTime();
271        if (ws==1){
272          li.SetLineColor(3);
273          li.DrawLine(obt,0,obt,max);
274        }
275        else if (ws==2){
276          li.SetLineColor(4);
277          li.DrawLine(obt,0,obt,max);
278        }  
279      }
280        
281    finalCanv1->cd(2);    RunHeadTr->GetEntry(0);
282    PcktLenght->SetStats(kFALSE);    ph = eH->GetPscuHeader();
283    PcktLenght->GetXaxis()->SetTitle("OBT    (ms)");    ULong_t TimeSync = reh->LAST_TIME_SYNC_INFO;
284    PcktLenght->GetXaxis()->CenterTitle();    ULong_t ObtSync = reh->OBT_TIME_SYNC;
285    PcktLenght->GetXaxis()->SetLabelSize(0.04);    //cout<<"TimeSync "<<reh->LAST_TIME_SYNC_INFO<<"\n";
286    PcktLenght->GetYaxis()->SetTitle("Lenght  (byte)");    //cout<<"ObtSync "<<reh->OBT_TIME_SYNC<<"\n";
287    PcktLenght->GetYaxis()->CenterTitle();  
288    PcktLenght->GetYaxis()->SetLabelSize(0.04);    finalCanv->cd(1);
289    PcktLenght->GetYaxis()->SetTitleSize(0.06);    stringstream ws1, ws2;
290    PcktLenght->GetYaxis()->SetTitleOffset(0.8);    ws1.str("");
291    PcktLenght->SetMarkerColor(2);    ws2.str("");
292    PcktLenght->SetMarkerSize(.5);    ws1<<"ACQ_SETTING_MODE = 1";
293    PcktLenght->SetMarkerStyle(21);    ws2<<"ACQ_SETTING_MODE = 2";
294    PcktLenght->Draw("9p");    TPaveText *pt=0;
295      pt = new TPaveText (.60,.92,.76,.98);
296    finalCanv1->cd(3);    pt->AddText(ws1.str().c_str());
297    PcktCounter->SetStats(kFALSE);    pt->SetTextColor(3);
298    PcktCounter->GetXaxis()->SetTitle("OBT    (ms)");    pt->SetFillColor(10);
299    PcktCounter->GetXaxis()->SetTitleSize(0.05);    pt->SetBorderSize(0);
300    PcktCounter->GetXaxis()->CenterTitle();    pt->Draw();
301    PcktCounter->GetXaxis()->SetLabelSize(0.04);    TPaveText *pt1=0;
302    PcktCounter->GetYaxis()->SetTitle("Counter");    pt1 = new TPaveText (.76,.92,.92,.98);
303    PcktCounter->GetYaxis()->SetTitleSize(0.05);    pt1->AddText(ws2.str().c_str());
304    PcktCounter->GetYaxis()->CenterTitle();    pt1->SetTextColor(4);
305    PcktCounter->GetYaxis()->SetLabelSize(0.04);    pt1->SetFillColor(10);
306    PcktCounter->GetYaxis()->SetTitleSize(0.06);    pt1->SetBorderSize(0);
307    PcktCounter->GetYaxis()->SetTitleOffset(0.8);    pt1->Draw();
308    PcktCounter->SetMarkerColor(4);    pt1 = new TPaveText (.05,.91,.6,1);
309    PcktCounter->SetMarkerSize(.5);    pt1->AddText(oss1.str().c_str());
310    PcktCounter->SetMarkerStyle(21);    pt1->SetTextColor(1);
311    PcktCounter->Draw("9p");    pt1->SetFillColor(10);
312      pt1->SetBorderSize(0);
313    oss.str("");    pt1->Draw();
314    oss << outDir.Data() << filename.Data();  
315    oss  << "_HeaderScan_1." << format.Data();    //----------- Second PAD -------------------------------//
316    finalCanv1->SaveAs(oss.str().c_str());    finalCanv->cd(2);
317    //------------- Graph PscuCounter vs FileOffset ---------------------------------    all3->Draw();
318    const Int_t size = 100000;    all3->cd();
319    Int_t ev=0;    packetLength->Draw("AP");
320    Double_t PscuCounter[size], FileOffset[size];    packetLength->GetXaxis()->SetTitle("OBT    (ms)");
321    Int_t m=0;    packetLength->GetXaxis()->SetTitleSize(0.05);
322    Int_t pag= (nevents/100000)+1;    packetLength->GetXaxis()->CenterTitle();
323      packetLength->GetXaxis()->SetLabelSize(0.05);
324      packetLength->GetYaxis()->SetTitle("Lenght  (byte)");
325      packetLength->GetYaxis()->CenterTitle();
326      packetLength->GetYaxis()->SetLabelSize(0.05);
327      packetLength->GetYaxis()->SetTitleSize(0.06);
328      packetLength->GetYaxis()->SetTitleOffset(0.6);
329    
330    
331      finalCanv->cd(2);
332      pt = new TPaveText (.6,.91,.90,1);
333      pt->AddText(oss2.str().c_str());
334      pt->SetTextColor(2);
335      pt->SetFillColor(10);
336      pt->SetBorderSize(0);
337      pt->Draw();
338    
339      //----------- Third PAD -------------------------------//
340      finalCanv->cd(3);
341      all4->Draw();
342      all4->cd();
343      packeCounter->Draw("AP");
344      packeCounter->GetXaxis()->SetTitle("OBT    (ms)");
345      packeCounter->GetXaxis()->SetTitleSize(0.05);
346      packeCounter->GetXaxis()->CenterTitle();
347      packeCounter->GetXaxis()->SetLabelSize(0.05);
348      packeCounter->GetYaxis()->SetTitle("Counter");
349      packeCounter->GetYaxis()->SetTitleSize(0.05);
350      packeCounter->GetYaxis()->CenterTitle();
351      packeCounter->GetYaxis()->SetLabelSize(0.05);
352      packeCounter->GetYaxis()->SetTitleSize(0.06);
353      packeCounter->GetYaxis()->SetTitleOffset(0.6);
354    
355    
356      finalCanv->cd(3);
357      TPaveText *pt2=0;
358      pt2 = new TPaveText (.6,.91,.90,1);
359      pt2->AddText(oss3.str().c_str());
360      pt2->SetTextColor(4);
361      pt2->SetFillColor(10);
362      pt2->SetBorderSize(0);
363      pt2->Draw();  
364    
365      /**********************************************************************************************/
366    
367    TCanvas *finalCanv2 = new TCanvas("Header_2", base, 1280, 1024);    TMultiGraph *mg1 = new TMultiGraph();
368    finalCanv2->SetFillColor(10);    TMultiGraph *mg2 = new TMultiGraph();
369      //---------------------- fill vectors and histogram --------------------------------------------------//
370      TList *list = new TList;
371      Int_t numkey;
372      TObject *key = new TObject;
373      const char *name;
374      char *SoftInfo="SoftInfo";
375      TTree* tr = new TTree;
376      Long64_t nevntskey=0;
377      list = file->GetListOfKeys();
378      numkey = file->GetNkeys();
379      ULong_t salto;
380      for (Int_t m=0; m<numkey; m++){
381        key = list->At(m);
382        name=(char *)(key->GetName());
383        if(strcmp(name,SoftInfo)==0)continue;
384        tr = (TTree*)file->Get(name);
385        if (tr->IsZombie()) continue;
386    
387        tr->SetBranchAddress("Header", &eh);
388        TBranch *Br = tr->GetBranch("Header");
389        nevntskey = tr->GetEntries();
390        if(nevntskey !=0){
391          Int_t size1=nevntskey;
392          Double_t *PscuCounter1 = new Double_t[size1];
393          Double_t *FileOffset1 = new Double_t[size1];
394          Double_t *tempo1 = new Double_t[size1];
395    
396          int n=0;
397          while(n<nevntskey){
398            Br->GetEntry(n);
399            ph = eh->GetPscuHeader();  
400            PscuCounter1[n]= ph->GetCounter();
401            FileOffset1[n]=ph->GetFileOffset();
402            tempo1[n]=ph->GetOrbitalTime();
403            if((m==0) && (n==0)){
404              primotempo=ph->GetOrbitalTime();
405              salto=ph->GetOrbitalTime();
406            }
407            if(salto > ph->GetOrbitalTime())salto=ph->GetOrbitalTime();
408            if(ph->GetFileOffset()<= primoffset){
409              primoffset=ph->GetFileOffset();
410              primotempo=ph->GetOrbitalTime();
411            }
412            if(ph->GetFileOffset()>=ultimoffset){
413              ultimotempo=ph->GetOrbitalTime();
414              ultimoffset=ph->GetFileOffset();
415            }
416            n++;
417          }
418    
419   /////ciclo sulle pagine////////////////////        TGraph *graph3= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)PscuCounter1);
420    for (Int_t interval=0; interval<pag; interval++){        graph3->SetMarkerColor(3);
421          graph3->SetMarkerSize(.2);
422          graph3->SetMarkerStyle(21);
423          mg1->Add(graph3);
424          
425          TGraph *graph4= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)tempo1);
426          graph4->SetMarkerColor(kBlue);
427          graph4->SetMarkerSize(.2);
428          graph4->SetMarkerStyle(21);
429          mg2->Add(graph4);
430        }
431      }
432    
433     for (Int_t l = 0; l < size; l++){     //ciclo su un sottogruppo    TLatex *lat = new TLatex();
434      headBr->GetEntry(ev+l);    lat->SetTextAlign(12);
435      ph = eh->GetPscuHeader();      lat->SetTextSize(0.15);
436      PscuCounter[l]= ph->GetCounter();    lat->SetTextColor(kBlue);
437      FileOffset[l]=ph->GetFileOffset();  
438      m=l;   //------------ Fourth PAD  ---------------------//
439     }       finalCanv->cd(4);
440       all1->Draw();
441       all1->cd();
442    
443     oss1.str("");     oss1.str("");
444     oss1 << filename.Data() <<": PscuCounter vs FileOffset. Canvas n."<< (interval+1);     oss1 <<"PscuCounter vs FileOffset.";
445     TGraph *graph= new TGraph(m, FileOffset, PscuCounter);     mg1->Draw("AP");
446       mg1->GetXaxis()->SetTitle("File Offset");
447       mg1->GetXaxis()->CenterTitle();
448       mg1->GetXaxis()->SetTitleOffset(0.8);
449       mg1->GetXaxis()->SetTitleSize(0.05);
450       mg1->GetXaxis()->SetLabelSize(0.05);
451       mg1->GetYaxis()->SetTitle("Counter");
452       mg1->GetYaxis()->CenterTitle();
453       mg1->GetYaxis()->SetTitleSize(0.06);
454       mg1->GetYaxis()->SetLabelSize(0.06);
455       mg1->GetYaxis()->SetTitleOffset(0.6);
456       finalCanv->cd(4);
457       TPaveText *pt3=0;
458       pt3 = new TPaveText (.60,.91,.90,1);
459       pt3->AddText(oss1.str().c_str());
460       pt3->SetTextColor(3);
461       pt3->SetFillColor(10);
462       pt3->SetBorderSize(0);
463       pt3->Draw();
464    
465     //------------ Fifth PAD  ---------------------//
466       finalCanv->cd(5);
467       all->Draw();
468       all->cd();
469       oss3.str("");
470       oss3 << "OBT vs FileOffset";
471       mg2->Draw("AP");
472       mg2->GetXaxis()->SetTitle("File Offset");
473       mg2->GetXaxis()->CenterTitle();
474       mg2->GetXaxis()->SetTitleSize(0.05);
475       mg2->GetXaxis()->SetLabelSize(0.05);
476       mg2->GetYaxis()->SetTitle("OBT");
477       mg2->GetYaxis()->CenterTitle();
478       mg2->GetYaxis()->SetTitleSize(0.06);
479       mg2->GetYaxis()->SetLabelSize(0.05);
480       mg2->GetYaxis()->SetTitleOffset(0.6);
481    
482       finalCanv->cd(5);
483       TPaveText *pt4=0;
484       pt4 = new TPaveText (.70,.91,.90,1);
485       pt4->AddText(oss3.str().c_str());
486       pt4->SetTextColor(kBlue);
487       pt4->SetFillColor(10);
488       pt4->SetBorderSize(0);
489       pt4->Draw();
490    
491       finalCanv->cd(6);
492       ULong_t primotempoABS=TimeSync+((primotempo/1000)-ObtSync);
493       ULong_t obmaxABS=TimeSync+((lastime/1000)-ObtSync);
494       ULong_t saltoABS=TimeSync+((salto/1000)-ObtSync);
495       ULong_t ultimotempoABS=TimeSync+((ultimotempo/1000)-ObtSync);
496    
497       TPaveText *pt5=0;
498       pt5 = new TPaveText (0,0,1,1);
499       stringa.str("");
500       stringa << "  Filename:   "<<filename.Data()<<"\n";
501       TText *t1=pt5->AddText(0.25,0.95,stringa.str().c_str());
502       t1->SetTextSize(0.1);
503       stringa.str("");
504       stringa << "  OBT (ms)            ABS TIME (s)";
505       TText *t2=pt5->AddText(0.32,0.75,stringa.str().c_str());
506       t2->SetTextSize(0.07);
507       stringa.str("");
508       stringa << "New data start at:   "<<primotempo<<"         "<<primotempoABS;
509       TText *t3=pt5->AddText(0.25,0.60,stringa.str().c_str());
510       t3->SetTextSize(0.08);
511       stringa.str("");
512       stringa << "New data end at:     "<<lastime<<"         "<<obmaxABS;
513       TText *t4=pt5->AddText(0.25,0.50,stringa.str().c_str());
514       t4->SetTextSize(0.08);
515       if(primotempo!=salto || lastime!=ultimotempo){
516         stringa.str("");
517         stringa << "Old data start at:   "<<salto<<"        "<<saltoABS;
518         TText *t5=pt5->AddText(0.65,0.60,stringa.str().c_str());
519         t5->SetTextSize(0.08);
520         stringa.str("");
521         stringa << "Old data end at:     "<<ultimotempo<<"        "<<ultimotempoABS;
522         TText *t6=pt5->AddText(0.65,0.50,stringa.str().c_str());
523         t6->SetTextSize(0.08);
524         stringa.str("");
525         stringa << "  OBT (ms)          ABS TIME (s)";
526         TText *t2=pt5->AddText(0.72,0.75,stringa.str().c_str());
527         t2->SetTextSize(0.07);
528       }
529       pt5->SetTextColor(kBlack);
530       pt5->SetFillColor(10);
531       pt5->SetBorderSize(0);
532       pt5->Draw();
533    
534     finalCanv2->cd();     finalCanv->Update();
    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);  
    graph->SetMarkerColor(3);  
    graph->SetMarkerSize(.5);  
    graph->SetMarkerStyle(21);  
    graph->Draw("AP");  
535    
536     oss1.str("");     oss1.str("");
537     oss1 << outDir.Data() << filename.Data();     oss1 << outDir.Data() << filename.Data();
538     oss1 << "_HeaderScan_2_"<<(interval+1)<<"." << format.Data();     oss1 << "_HeaderScan"<<"." << format.Data();
539    
540     finalCanv2->SaveAs(oss1.str().c_str());     finalCanv->SaveAs(oss1.str().c_str());
    finalCanv2->Clear();  
541    
542     ev=ev+100000;     file->Close();
   }  
   
   file->Close();  
543        
544  }  }
545    
# Line 334  int main(int argc, char* argv[]){ Line 586  int main(int argc, char* argv[]){
586    }    }
587    HeaderScan(argv[1], outDir, format);    HeaderScan(argv[1], outDir, format);
588  }  }
589    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.23