/[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.1.1.1 by pam-rm2, Fri Jun 16 10:11:54 2006 UTC revision 1.14 by pam-rm2, Fri Mar 9 10:43:40 2007 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 43  Line 43 
43  *                                      There is no default value, without this input the program will not run  *                                      There is no default value, without this input the program will not run
44  *              TString outDir        - the path where to save the output image (Default = ./)  *              TString outDir        - the path where to save the output image (Default = ./)
45  *              TString format        - the format which will be used fo rsave the produced images (Default = "jpg")  *              TString format        - the format which will be used fo rsave the produced images (Default = "jpg")
46    *
47    *
48    * Version 1.7
49    * Date 16 June 2006 - Malvezzi
50    *
51    * Description of changes:
52    *                         Implementation of the case: numebr of events <= 0.
53    *                         Remove graph "grPcktId1"; see PacketScan for the same information.
54    *             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.
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    * Version 1.9
63    * Date 7 March 2007 - Malvezzi
64    *
65    * Description: to read correctly the OBT changed the type of time variables from signed to unsigned
66    *
67  */  */
68    
69    
70    #include <fstream>
71    #include <math.h>
72    #include "TLatex.h"
73    #include "TF1.h"
74    #include "TPaveText.h"
75    #include "TMultiGraph.h"
76  #include <sstream>  #include <sstream>
77  #include <iostream>  #include <iostream>
78    #include <stdio.h>
79    #include <string.h>
80  #include "TString.h"  #include "TString.h"
81  #include "TStyle.h"  #include "TStyle.h"
82  #include "TFile.h"  #include "TFile.h"
# Line 61  Line 89 
89  #include "TGaxis.h"  #include "TGaxis.h"
90  #include "EventHeader.h"  #include "EventHeader.h"
91  #include "PscuHeader.h"  #include "PscuHeader.h"
92    #include "RunHeaderEvent.h"
93    #include "TPaveText.h"
94    
95  using namespace std;  using namespace std;
96    
97  void HeaderScan(TString base, TString outDir, TString format)  void HeaderScan(TString base, TString outDir, TString format){  
 {  
   
   TList *list = new TList;  
   Int_t numkey;  
   TObject *key = new TObject;  
   const char *name;  
   TTree* tr = new TTree;  
   Long64_t totevents=0, totphysevents=0;  
   Float_t id;  
   Long64_t nevents=0;  
98    
99    //------- load root file --------------    //------------------- Variables initilization -------------------------//
100      Long64_t nevents=0, runnevents;
101      ULong_t lastime, firstime, primotempo, ultimotempo, primoffset=500000000, ultimoffset;
102      double obmin=0.;
103      double obmax=0.;
104      stringstream oss, oss1, oss2, oss3, noentries, stringa;
105    
106      //-------------- Load root file, tree and branch  -------------------//
107    TFile *file = new TFile(base.Data());    TFile *file = new TFile(base.Data());
     
108    if (!file){    if (!file){
109      printf("No such file in the directory has been found");      printf("No such file in the directory has been found");
110      return;      return;
111    }    }
112    if (outDir == "" ) outDir = ".";  
113        TTree *PhysicsTr = (TTree*)file->Get("Physics");
114    list = file->GetListOfKeys();    TBranch *headBr = PhysicsTr->GetBranch("Header");
115      
116    numkey = file->GetNkeys();    pamela::EventHeader  *eh   = new pamela::EventHeader;
117        pamela::PscuHeader   *ph   = new pamela::PscuHeader;
118    pamela::EventHeader *eh=0;  
119    pamela::PscuHeader *ph=0;    PhysicsTr->SetBranchAddress("Header", &eh);
120        nevents = PhysicsTr->GetEntries();
121    ///-----to know the total number f events end of physics events----//////    const Int_t size = nevents;
122      
123    for (Int_t i=0; i<numkey; i++){    TTree *RunHeadTr = (TTree*)file->Get("RunHeader");  ///run header tree
124      key = list->At(i);    pamela::EventHeader *eH= new pamela::EventHeader;
125      name=(key->GetName());    pamela::RunHeaderEvent *reh=new pamela::RunHeaderEvent;
126      //cout<<name<<"\n";    
127      tr = (TTree*)file->Get(name);    RunHeadTr->SetBranchAddress("Header",&eH);
128      if (tr->IsZombie()) continue;    RunHeadTr->SetBranchAddress("RunHeader",&reh);
129      nevents = tr->GetEntries();    runnevents = RunHeadTr->GetEntries();
     totevents+=nevents;  
     tr->SetBranchAddress("Header", &eh);    
     for (Int_t j = 0; j < nevents; j++){  
       tr->GetEntry(j);  
       ph = eh->GetPscuHeader();  
       if(ph->GetPacketId1() == 0x10) ++totphysevents;  
     }  
       
   }  
130    
   const Long64_t totalevents=totevents;  
   const Long64_t totalphysevents=totphysevents;  
     
   /////////////////////////////////////////////////  
     
131    TString filename = ((TObjString*)base.Tokenize('/')->Last())->GetString();    TString filename = ((TObjString*)base.Tokenize('/')->Last())->GetString();
132    filename = ((TObjString*)filename.Tokenize('.')->First())->GetString();    filename = ((TObjString*)filename.Tokenize('.')->First())->GetString();
   stringstream oss;  
   oss.str("");  
   oss << filename.Data();  
     
   Float_t obt[totalevents], pcktId[totalevents], pcktLenght[totalphysevents], pcktCounter[totalphysevents], offset[totalevents], allCounter[totalevents];  
   Float_t obtphysevents=0;  
     
133    
134    totphysevents=0;    //----------- If nevents < = 0 ---------------------------------/
135    totevents=0;    if (nevents<=0) {
136    for (Int_t i=0; i<numkey; i++){      printf("nevents = %i \n", nevents);
137      key = list->At(i);      printf(" \n");
138      name=(char *)(key->GetName());      TCanvas *canv = new TCanvas("No entries", "No entries ", 400, 200);
139      tr = (TTree*)file->Get(name);      canv->SetFillColor(10);
140      if (tr->IsZombie()) continue;      canv->cd();
141      tr->SetBranchAddress("Header", &eh);      TLatex *l = new TLatex();
142      nevents = tr->GetEntries();      l->SetTextAlign(12);
143      for (Int_t j = 0; j < nevents; j++){      l->SetTextSize(0.15);
144        tr->GetEntry(j);        l->SetTextColor(2);
145        ph = eh->GetPscuHeader();      noentries.str("");
146        obt[j+totevents]=ph->GetOrbitalTime() ;      noentries<< "HeaderScan_QL:";
147        pcktId[j+totevents]=ph->GetPacketId1();      l->DrawLatex(0.05, 0.7, noentries.str().c_str());
148        offset[j+totevents]=ph->GetFileOffset();      noentries.str("");
149        if(pcktId[j+totevents]==16){      noentries<< "No Physics entries for this files";
150          pcktLenght[totphysevents]=ph->GetPacketLenght();      l->DrawLatex(0.05, 0.5, noentries.str().c_str());
151          pcktCounter[totphysevents]=ph->GetCounter();      if (outDir == "./") {
152          totphysevents=totphysevents+1;        oss.str("");
153        }        oss << filename.Data() << "_HeaderScan_QL." << format.Data();
154        offset[j+totevents]=ph->GetFileOffset();      } else {
155        allCounter[j+totevents]= ph->GetCounter();        oss.str("");
156      }        oss << outDir.Data() << filename.Data() << "_HeaderScan_QL." << format.Data();
157      totevents=totevents+nevents;      }  
158        canv->Update();
159        canv->SaveAs(oss.str().c_str());
160        return;
161    }    }
162        //-------------- to know the max and min OBT ----------------------------//
163         headBr->GetEntry(0);
164    Float_t  mintime=obt[0], maxtime=obt[0], minlen=pcktLenght[0], maxlen=pcktLenght[0], mincount=pcktCounter[0], maxcount=pcktCounter[0];       ph = eh->GetPscuHeader();
165         firstime = ph->GetOrbitalTime();
166         headBr->GetEntry(nevents);
167    for (Int_t t=0; t<totalevents; t++){     ph = eh->GetPscuHeader();
168      if(obt[t]<mintime) mintime=obt[t];     lastime = ph->GetOrbitalTime();
169      if(obt[t]>maxtime) maxtime=obt[t];     int i =0;
170    }     while(i < nevents){
171    for (Int_t t=0; t<totalphysevents; t++){       headBr->GetEntry(i);
172      if(pcktLenght[t]<=minlen) minlen=pcktLenght[t];       ph = eh->GetPscuHeader();
173      if(pcktLenght[t]>=maxlen) maxlen=pcktLenght[t];       if((ph->GetOrbitalTime()) <= firstime) firstime=ph->GetOrbitalTime();
174      if(pcktCounter[t]<=mincount) mincount=pcktCounter[t];       if((ph->GetOrbitalTime()) >= lastime) lastime=ph->GetOrbitalTime();
175      if(pcktCounter[t]>=maxcount) maxcount=pcktCounter[t];       i++;
176         }
177      
178      //------------------------ First histogram -----------------------------------//
179      obmin=firstime;
180      obmax=lastime;
181      Int_t  nbin   = (lastime-firstime)/60000;
182      TH1F *h1 = new TH1F ("histo1", "" , nbin, obmin, obmax);
183    
184      //------------------- fill vectors and histogram -----------------------------//
185      Double_t *PscuCounter = new Double_t[size];
186      Double_t *PacketLenght = new Double_t[size];
187      ULong_t utile=0;
188      double *OBTime = new Double_t[size];
189      Double_t *Eventsperminute= new Double_t[nbin];
190      Double_t *Minute= new Double_t[nbin];
191      Double_t max=0;
192      for (Int_t k = 0; k < nevents; k++){
193        headBr->GetEntry(k);
194        ph = eh->GetPscuHeader();
195        utile= ph->GetOrbitalTime();
196        h1->Fill(utile);
197        PscuCounter[k]= ph->GetCounter();
198        PacketLenght[k]=ph->GetPacketLenght();
199        OBTime[k]=utile;
200      }  
201      int l=0;
202      while(l<nbin){
203        Eventsperminute[l]=h1->GetBinContent(l);
204        Minute[l]=firstime+l*60000;
205        if(h1->GetBinContent(l) >= max)max =h1->GetBinContent(l);
206        l++;
207    }    }
     
     
   Float_t    nbintime   = (maxtime-mintime)/1000;  
   Float_t    nbinlen    = (maxlen-minlen)/100;  
   Float_t    nbincount  = (maxcount-mincount)/100;  
     
     
   ///---------------------------TO GRAPH---------------------------------------------///  
   TCanvas *finalCanv1 = new TCanvas("Header_1", base, 1280, 1024);  
   TCanvas *finalCanv2 = new TCanvas("Header_2", base, 1280, 1024);  
   finalCanv1->Divide(3);  
   finalCanv1->SetFillColor(10);  
   finalCanv2->Divide(2);  
   finalCanv2->SetFillColor(10);  
208    
209        //----------- Graph and MultiGraph -----------------------------------------------//
210    TGraph *grPcktId1 = new TGraph(totalevents, obt, pcktId);    TMultiGraph *rate = new TMultiGraph();
211    TGraph *grPcktLenght = new TGraph(totalevents, obt, pcktLenght);    TMultiGraph *packetLength = new TMultiGraph();
212    TGraph *grPcktCounter = new TGraph(totalevents, obt, pcktCounter);    TMultiGraph *packeCounter = new TMultiGraph();
   TGraph *grFileOffset = new TGraph(totalevents, offset, allCounter);  
     
   TGaxis::SetMaxDigits(3);  
213    
   //-----canvas 1-------------------------------///  
   finalCanv1->cd(1);  
     
   stringstream oss1;  
   oss1.str("");  
   oss1 << "OBT (ms) with t0 = " << mintime << "ms";  
     
   gStyle->SetTitleH(0.06);  
   grPcktId1->SetTitle(oss.str().c_str());  
   grPcktId1->GetXaxis()->SetTitle(oss1.str().c_str());  
   grPcktId1->GetXaxis()->CenterTitle();  
   grPcktId1->GetXaxis()->SetLabelSize(0.03);  
   grPcktId1->GetYaxis()->SetTitle("Packet ID");  
   grPcktId1->GetYaxis()->CenterTitle();  
   grPcktId1->SetMarkerSize(4);  
   grPcktId1->Draw("AP");  
     
     
   finalCanv1->cd(2);  
214    oss1.str("");    oss1.str("");
215    oss1 << "OBT (min) with t0 = " << mintime << "ms";    oss1 << "Physics Packet per minute.  Start time = " << obmin << ", End time = "<< obmax <<" ms";
216        TGraph *rate1= new TGraph(nbin, (const Double_t*)Minute, (const Double_t*)Eventsperminute);
217    TH1F *h1 = new TH1F ("h1", oss.str().c_str(), (int)(nbintime/60), mintime, maxtime);    rate1->SetMarkerColor(kBlack);
218      rate1->SetMarkerSize(.1);
219      rate1->SetMarkerStyle(21);
220      rate->Add(rate1);
221    
222      TGraph *packetLength1= new TGraph(nevents, (const double*)OBTime, (const double*)PacketLenght);
223      oss2.str("");
224      oss2 <<"Lenght of Physic packets";
225      packetLength1->SetMarkerColor(2);
226      packetLength1->SetMarkerSize(.3);
227      packetLength1->SetMarkerStyle(21);
228      packetLength->Add(packetLength1);
229    
230      TGraph *packeCounter1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PscuCounter);
231      oss3.str("");
232      oss3 <<"Physics Counter vs. OBT";
233      packeCounter1->SetMarkerColor(4);
234      packeCounter1->SetMarkerSize(.2);
235      packeCounter1->SetMarkerStyle(21);
236      packeCounter->Add(packeCounter1);
237    
238      //------------ Create and Draw Canvas ---------------------//
239      TCanvas *finalCanv = new TCanvas("Header", base, 1200, 1600);
240      finalCanv->Divide(1,6);
241      finalCanv->SetFillColor(10);
242    
243      TPad *all2= new TPad ("","", 0, 0, 1, 1);
244      all2->SetFillColor(10);
245      TPad *all3= new TPad ("","", 0, 0, 1, 1);
246      all3->SetFillColor(10);
247      TPad *all4= new TPad ("","", 0, 0, 1, 1);
248      all4->SetFillColor(10);
249      TPad *all= new TPad ("","", 0, 0, 1, 1);
250      all->SetFillColor(10);
251      TPad *all1= new TPad ("","", 0, 0, 1, 1);
252      all1->SetFillColor(10);
253    
254      TLine li;
255      li.SetLineStyle(4);
256      li.SetLineWidth(2);
257    
258      //----------- First PAD -------------------------------//
259      finalCanv->cd(1);
260      all2->Draw();
261      all2->cd();
262      rate->Draw("ALP");
263      rate->GetXaxis()->SetTitle("OBT   (ms)");
264      rate->GetXaxis()->SetTitleSize(0.05);
265      rate->GetXaxis()->CenterTitle();
266      rate->GetXaxis()->SetLabelSize(0.05);
267      rate->GetYaxis()->SetTitle("Number of events ");
268      rate->GetYaxis()->CenterTitle();
269      rate->GetYaxis()->SetLabelSize(0.05);
270      rate->GetYaxis()->SetTitleSize(0.06);
271      rate->GetYaxis()->SetTitleOffset(0.6);
272      for (Int_t l = 0; l < runnevents; l++){
273        RunHeadTr->GetEntry(l);
274        ph = eH->GetPscuHeader();
275        int ws= reh->RM_ACQ_SETTING_MODE;
276        int id = ph->GetPacketId1();
277        ULong_t utile1=ph->GetOrbitalTime();
278        Double_t obt = utile1;
279        if (ws==1){
280          li.SetLineColor(3);
281          li.DrawLine(obt,0,obt,max);
282        }
283        else if (ws==2){
284          li.SetLineColor(4);
285          li.DrawLine(obt,0,obt,max);
286        }  
287      }
288        
289    for (Int_t i=0; i<numkey; i++){    RunHeadTr->GetEntry(0);
290      key = list->At(i);    ph = eH->GetPscuHeader();
291      ULong_t TimeSync = reh->LAST_TIME_SYNC_INFO;
292      ULong_t ObtSync = reh->OBT_TIME_SYNC;
293      //cout<<"TimeSync "<<reh->LAST_TIME_SYNC_INFO<<"\n";
294      //cout<<"ObtSync "<<reh->OBT_TIME_SYNC<<"\n";
295    
296      finalCanv->cd(1);
297      stringstream ws1, ws2;
298      ws1.str("");
299      ws2.str("");
300      ws1<<"ACQ_SETTING_MODE = 1";
301      ws2<<"ACQ_SETTING_MODE = 2";
302      TPaveText *pt=0;
303      pt = new TPaveText (.60,.92,.76,.98);
304      pt->AddText(ws1.str().c_str());
305      pt->SetTextColor(3);
306      pt->SetFillColor(10);
307      pt->SetBorderSize(0);
308      pt->Draw();
309      TPaveText *pt1=0;
310      pt1 = new TPaveText (.76,.92,.92,.98);
311      pt1->AddText(ws2.str().c_str());
312      pt1->SetTextColor(4);
313      pt1->SetFillColor(10);
314      pt1->SetBorderSize(0);
315      pt1->Draw();
316      pt1 = new TPaveText (.05,.91,.6,1);
317      pt1->AddText(oss1.str().c_str());
318      pt1->SetTextColor(1);
319      pt1->SetFillColor(10);
320      pt1->SetBorderSize(0);
321      pt1->Draw();
322    
323      //----------- Second PAD -------------------------------//
324      finalCanv->cd(2);
325      all3->Draw();
326      all3->cd();
327      packetLength->Draw("AP");
328      packetLength->GetXaxis()->SetTitle("OBT    (ms)");
329      packetLength->GetXaxis()->SetTitleSize(0.05);
330      packetLength->GetXaxis()->CenterTitle();
331      packetLength->GetXaxis()->SetLabelSize(0.05);
332      packetLength->GetYaxis()->SetTitle("Lenght  (byte)");
333      packetLength->GetYaxis()->CenterTitle();
334      packetLength->GetYaxis()->SetLabelSize(0.05);
335      packetLength->GetYaxis()->SetTitleSize(0.06);
336      packetLength->GetYaxis()->SetTitleOffset(0.6);
337    
338    
339      finalCanv->cd(2);
340      pt = new TPaveText (.6,.91,.90,1);
341      pt->AddText(oss2.str().c_str());
342      pt->SetTextColor(2);
343      pt->SetFillColor(10);
344      pt->SetBorderSize(0);
345      pt->Draw();
346    
347      //----------- Third PAD -------------------------------//
348      finalCanv->cd(3);
349      all4->Draw();
350      all4->cd();
351      packeCounter->Draw("AP");
352      packeCounter->GetXaxis()->SetTitle("OBT    (ms)");
353      packeCounter->GetXaxis()->SetTitleSize(0.05);
354      packeCounter->GetXaxis()->CenterTitle();
355      packeCounter->GetXaxis()->SetLabelSize(0.05);
356      packeCounter->GetYaxis()->SetTitle("Counter");
357      packeCounter->GetYaxis()->SetTitleSize(0.05);
358      packeCounter->GetYaxis()->CenterTitle();
359      packeCounter->GetYaxis()->SetLabelSize(0.05);
360      packeCounter->GetYaxis()->SetTitleSize(0.06);
361      packeCounter->GetYaxis()->SetTitleOffset(0.6);
362    
363    
364      finalCanv->cd(3);
365      TPaveText *pt2=0;
366      pt2 = new TPaveText (.6,.91,.90,1);
367      pt2->AddText(oss3.str().c_str());
368      pt2->SetTextColor(4);
369      pt2->SetFillColor(10);
370      pt2->SetBorderSize(0);
371      pt2->Draw();  
372    
373      /**********************************************************************************************/
374    
375      TMultiGraph *mg1 = new TMultiGraph();
376      TMultiGraph *mg2 = new TMultiGraph();
377      //---------------------- fill vectors and histogram --------------------------------------------------//
378      TList *list = new TList;
379      Int_t numkey;
380      TObject *key = new TObject;
381      const char *name;
382      char *SoftInfo="SoftInfo";
383      TTree* tr = new TTree;
384      Long64_t nevntskey=0;
385      list = file->GetListOfKeys();
386      numkey = file->GetNkeys();
387      ULong_t salto;
388      for (Int_t m=0; m<numkey; m++){
389        key = list->At(m);
390      name=(char *)(key->GetName());      name=(char *)(key->GetName());
391      tr = (TTree*)file->Get(name);      if(strcmp(name,SoftInfo)==0)continue;
392        tr = (TTree*)file->Get(name);
393      if (tr->IsZombie()) continue;      if (tr->IsZombie()) continue;
394    
395      tr->SetBranchAddress("Header", &eh);      tr->SetBranchAddress("Header", &eh);
396      nevents = tr->GetEntries();      TBranch *Br = tr->GetBranch("Header");
397      for (Int_t j = 0; j < nevents; j++){      nevntskey = tr->GetEntries();
398        tr->GetEntry(j);        if(nevntskey !=0){
399        ph = eh->GetPscuHeader();        Int_t size1=nevntskey;
400        if((ph->GetPacketId1()) == 16){        Double_t *PscuCounter1 = new Double_t[size1];
401          obtphysevents=ph->GetOrbitalTime();        Double_t *FileOffset1 = new Double_t[size1];
402          h1->Fill(obtphysevents);        Double_t *tempo1 = new Double_t[size1];
403        }  
404      }        int n=0;
405          while(n<nevntskey){
406            Br->GetEntry(n);
407            ph = eh->GetPscuHeader();  
408            PscuCounter1[n]= ph->GetCounter();
409            FileOffset1[n]=ph->GetFileOffset();
410            utile=ph->GetOrbitalTime();
411            tempo1[n]=utile;
412            if((m==0) && (n==0)){
413              primotempo=ph->GetOrbitalTime();
414              salto=ph->GetOrbitalTime();
415            }
416            if(salto > ph->GetOrbitalTime())salto=ph->GetOrbitalTime();
417            if(ph->GetFileOffset()<= primoffset){
418              primoffset=ph->GetFileOffset();
419              primotempo=ph->GetOrbitalTime();
420            }
421            if(ph->GetFileOffset()>=ultimoffset){
422              ultimotempo=ph->GetOrbitalTime();
423              ultimoffset=ph->GetFileOffset();
424            }
425            n++;
426          }
427    
428          TGraph *graph3= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)PscuCounter1);
429          graph3->SetMarkerColor(3);
430          graph3->SetMarkerSize(.2);
431          graph3->SetMarkerStyle(21);
432          mg1->Add(graph3);
433          
434          TGraph *graph4= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)tempo1);
435          graph4->SetMarkerColor(kBlue);
436          graph4->SetMarkerSize(.2);
437          graph4->SetMarkerStyle(21);
438          mg2->Add(graph4);
439        }
440    }    }
     
   h1->SetMarkerStyle(6);  
   h1->GetXaxis()->SetTitle(oss1.str().c_str());  
   h1->GetXaxis()->CenterTitle();  
   h1->GetXaxis()->SetLabelSize(0.03);  
   h1->GetYaxis()->SetTitle("number of Physic packets");  
   h1->GetYaxis()->CenterTitle();  
   h1->Draw();  
     
     
   finalCanv1->cd(3);  
   oss1.str("");  
   oss1 << "OBT (ms) with t0 = " << mintime << "ms";  
     
   grPcktLenght->SetTitle(oss.str().c_str());  
   grPcktLenght->GetXaxis()->SetTitle(oss1.str().c_str());  
   grPcktLenght->GetXaxis()->CenterTitle();  
   grPcktLenght->GetXaxis()->SetLabelSize(0.03);  
   grPcktLenght->GetYaxis()->SetTitle("Lenght of Physic packets   (byte)");  
   grPcktLenght->GetYaxis()->CenterTitle();  
   grPcktLenght->GetYaxis()->SetLabelSize(0.03);  
   grPcktLenght->SetMarkerSize(4);  
   grPcktLenght->Draw("AP");  
     
     
   finalCanv1->cd(2);  
   h1->Draw();  
     
   finalCanv1->Update();  
     
     
   ///---------canvas 2-----------------------------//  
   finalCanv2->cd(1);  
     
   grPcktCounter->SetTitle(oss.str().c_str());  
   grPcktCounter->GetXaxis()->SetTitle(oss1.str().c_str());  
   grPcktCounter->GetXaxis()->SetTitleSize(0.04);  
   grPcktCounter->GetXaxis()->CenterTitle();  
   grPcktCounter->GetXaxis()->SetLabelSize(0.03);  
   grPcktCounter->GetYaxis()->SetTitle("PacketCounter of Physic packets");  
   grPcktCounter->GetYaxis()->SetTitleSize(0.04);  
   grPcktCounter->GetYaxis()->CenterTitle();  
   grPcktCounter->GetYaxis()->SetLabelSize(0.03);  
   grPcktCounter->SetMarkerSize(4);  
   grPcktCounter->SetMinimum(mincount);  
   grPcktCounter->SetMaximum(maxcount);  
   grPcktCounter->Draw("AP");  
     
   finalCanv2->cd(2);  
     
   grFileOffset->SetTitle(oss.str().c_str());  
   grFileOffset->GetXaxis()->SetTitle("File Offset");  
   grFileOffset->GetXaxis()->CenterTitle();  
   grFileOffset->GetXaxis()->SetTitleSize(0.04);  
   grFileOffset->GetXaxis()->SetLabelSize(0.03);  
   grFileOffset->GetYaxis()->SetTitle("Packet counter");  
   grFileOffset->GetYaxis()->CenterTitle();  
   grFileOffset->GetYaxis()->SetTitleSize(0.04);  
   grFileOffset->GetYaxis()->SetLabelSize(0.03);  
   grFileOffset->SetMarkerSize(4);  
   grFileOffset->Draw("AP");  
     
   finalCanv1->Update();  
     
     
   //-------to save---------------------------///  
   oss.str("");  
   oss1.str("");  
     
   oss << outDir.Data() << filename.Data();  
   oss1 << outDir.Data() << filename.Data();  
441    
442    oss  << "_HeaderScan_1." << format.Data();    TLatex *lat = new TLatex();
443    oss1 << "_HeaderScan_2." << format.Data();    lat->SetTextAlign(12);
444        lat->SetTextSize(0.15);
445    finalCanv1->SaveAs(oss.str().c_str());    lat->SetTextColor(kBlue);
446    finalCanv2->SaveAs(oss1.str().c_str());  
447       //------------ Fourth PAD  ---------------------//
448         finalCanv->cd(4);
449    file->Close();     all1->Draw();
450         all1->cd();
451  }  
452       oss1.str("");
453       oss1 <<"PscuCounter vs FileOffset.";
454       mg1->Draw("AP");
455       mg1->GetXaxis()->SetTitle("File Offset");
456       mg1->GetXaxis()->CenterTitle();
457       mg1->GetXaxis()->SetTitleOffset(0.8);
458       mg1->GetXaxis()->SetTitleSize(0.05);
459       mg1->GetXaxis()->SetLabelSize(0.05);
460       mg1->GetYaxis()->SetTitle("Counter");
461       mg1->GetYaxis()->CenterTitle();
462       mg1->GetYaxis()->SetTitleSize(0.06);
463       mg1->GetYaxis()->SetLabelSize(0.06);
464       mg1->GetYaxis()->SetTitleOffset(0.6);
465       finalCanv->cd(4);
466       TPaveText *pt3=0;
467       pt3 = new TPaveText (.60,.91,.90,1);
468       pt3->AddText(oss1.str().c_str());
469       pt3->SetTextColor(3);
470       pt3->SetFillColor(10);
471       pt3->SetBorderSize(0);
472       pt3->Draw();
473    
474     //------------ Fifth PAD  ---------------------//
475       finalCanv->cd(5);
476       all->Draw();
477       all->cd();
478       oss3.str("");
479       oss3 << "OBT vs FileOffset";
480       mg2->Draw("AP");
481       mg2->GetXaxis()->SetTitle("File Offset");
482       mg2->GetXaxis()->CenterTitle();
483       mg2->GetXaxis()->SetTitleSize(0.05);
484       mg2->GetXaxis()->SetLabelSize(0.05);
485       mg2->GetYaxis()->SetTitle("OBT");
486       mg2->GetYaxis()->CenterTitle();
487       mg2->GetYaxis()->SetTitleSize(0.06);
488       mg2->GetYaxis()->SetLabelSize(0.05);
489       mg2->GetYaxis()->SetTitleOffset(0.6);
490    
491       finalCanv->cd(5);
492       TPaveText *pt4=0;
493       pt4 = new TPaveText (.70,.91,.90,1);
494       pt4->AddText(oss3.str().c_str());
495       pt4->SetTextColor(kBlue);
496       pt4->SetFillColor(10);
497       pt4->SetBorderSize(0);
498       pt4->Draw();
499    
500       finalCanv->cd(6);
501       ULong_t primotempoABS=TimeSync+((primotempo/1000)-ObtSync);
502       ULong_t obmaxABS=TimeSync+((lastime/1000)-ObtSync);
503       ULong_t saltoABS=TimeSync+((salto/1000)-ObtSync);
504       ULong_t ultimotempoABS=TimeSync+((ultimotempo/1000)-ObtSync);
505    
506       TPaveText *pt5=0;
507       pt5 = new TPaveText (0,0,1,1);
508       stringa.str("");
509       stringa << "  Filename:   "<<filename.Data()<<"\n";
510       TText *t1=pt5->AddText(0.25,0.95,stringa.str().c_str());
511       t1->SetTextSize(0.1);
512       stringa.str("");
513       stringa << "  OBT (ms)            ABS TIME (s)";
514       TText *t2=pt5->AddText(0.32,0.75,stringa.str().c_str());
515       t2->SetTextSize(0.07);
516       stringa.str("");
517       stringa << "New data start at:   "<<primotempo<<"         "<<primotempoABS;
518       TText *t3=pt5->AddText(0.25,0.60,stringa.str().c_str());
519       t3->SetTextSize(0.08);
520       stringa.str("");
521       stringa << "New data end at:     "<<lastime<<"         "<<obmaxABS;
522       TText *t4=pt5->AddText(0.25,0.50,stringa.str().c_str());
523       t4->SetTextSize(0.08);
524       if(primotempo!=salto || lastime!=ultimotempo){
525         stringa.str("");
526         stringa << "Old data start at:   "<<salto<<"        "<<saltoABS;
527         TText *t5=pt5->AddText(0.65,0.60,stringa.str().c_str());
528         t5->SetTextSize(0.08);
529         stringa.str("");
530         stringa << "Old data end at:     "<<ultimotempo<<"        "<<ultimotempoABS;
531         TText *t6=pt5->AddText(0.65,0.50,stringa.str().c_str());
532         t6->SetTextSize(0.08);
533         stringa.str("");
534         stringa << "  OBT (ms)          ABS TIME (s)";
535         TText *t2=pt5->AddText(0.72,0.75,stringa.str().c_str());
536         t2->SetTextSize(0.07);
537       }
538       pt5->SetTextColor(kBlack);
539       pt5->SetFillColor(10);
540       pt5->SetBorderSize(0);
541       pt5->Draw();
542    
543       finalCanv->Update();
544    
545       oss1.str("");
546       oss1 << outDir.Data() << filename.Data();
547       oss1 << "_HeaderScan"<<"." << format.Data();
548    
549       finalCanv->SaveAs(oss1.str().c_str());
550    
551       file->Close();
552      
553    }
554    
555  int main(int argc, char* argv[]){  int main(int argc, char* argv[]){
556    TString path;    TString path;
557    TString outDir     = "./";    TString outDir     = "./";
558    TString format     = "jpg";    TString format     = "jpg";
       
559   if (argc < 2){   if (argc < 2){
560      printf("You have to insert at least the file to analyze \n");      printf("You have to insert at least the file to analyze \n");
561      printf("Try '--help' for more information. \n");      printf("Try '--help' for more information. \n");
562      exit(1);      exit(1);
563    }      }  
   
564    if (!strcmp(argv[1], "--help")){    if (!strcmp(argv[1], "--help")){
565          printf( "Usage: HeaderScan FILE [OPTION] \n");          printf( "Usage: HeaderScan FILE [OPTION] \n");
566          printf( "\t --help                 Print this help and exit \n");          printf( "\t --help                 Print this help and exit \n");
# Line 334  int main(int argc, char* argv[]){ Line 568  int main(int argc, char* argv[]){
568          printf( "\t -format[jpg|ps|gif]    Format for output files [default 'jpg'] \n");          printf( "\t -format[jpg|ps|gif]    Format for output files [default 'jpg'] \n");
569          exit(1);          exit(1);
570    }    }
   
     
571    path=argv[1];    path=argv[1];
572      for (int i = 2; i < argc; i++){
   for (int i = 2; i < argc; i++){  
     
573      if (!strcmp(argv[i], "-outDir")){      if (!strcmp(argv[i], "-outDir")){
574        if (++i >= argc){        if (++i >= argc){
575          printf( "-outDir needs arguments. \n");          printf( "-outDir needs arguments. \n");
# Line 350  int main(int argc, char* argv[]){ Line 580  int main(int argc, char* argv[]){
580          outDir = argv[i];          outDir = argv[i];
581          continue;          continue;
582        }        }
583      }      }  
       
584      if (!strcmp(argv[i], "-format")){      if (!strcmp(argv[i], "-format")){
585        if (++i >= argc){        if (++i >= argc){
586          printf( "-format needs arguments. \n");          printf( "-format needs arguments. \n");
# Line 363  int main(int argc, char* argv[]){ Line 592  int main(int argc, char* argv[]){
592          continue;          continue;
593        }        }
594      }      }
     
   
595    }    }
   
596    HeaderScan(argv[1], outDir, format);    HeaderScan(argv[1], outDir, format);
   
597  }  }
598    
Legend:
Removed from v.1.1.1.1  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.23