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