/[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 by pam-rm2, Fri Jun 16 10:11:54 2006 UTC revision 1.9 by pam-rm2, Fri Sep 8 08:26:47 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 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  */  */
63    
64    
65    #include <fstream>
66    #include <math.h>
67    #include "TLatex.h"
68    #include "TF1.h"
69    #include "TPaveText.h"
70    #include "TMultiGraph.h"
71  #include <sstream>  #include <sstream>
72  #include <iostream>  #include <iostream>
73  #include "TString.h"  #include "TString.h"
# Line 61  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){  
 {  
   
   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;  
91    
92    //------- load root file --------------    //------------------- Variables initilization -------------------------//
93      Long64_t nevents=0, runnevents;
94      ULong_t lastime, firstime, primotempo, ultimotempo, primoffset=500000000, ultimoffset;
95      double obmin=0.;
96      double obmax=0.;
97      stringstream oss, oss1, oss2, oss3, noentries, stringa;
98    
99      //-------------- 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    if (outDir == "" ) outDir = ".";  
106        TTree *PhysicsTr = (TTree*)file->Get("Physics");
107    list = file->GetListOfKeys();    TBranch *headBr = PhysicsTr->GetBranch("Header");
108      
109    numkey = file->GetNkeys();    pamela::EventHeader  *eh   = 0;
110        pamela::PscuHeader   *ph   = 0;
111    pamela::EventHeader *eh=0;  
112    pamela::PscuHeader *ph=0;    PhysicsTr->SetBranchAddress("Header", &eh);
113        nevents = PhysicsTr->GetEntries();
114    ///-----to know the total number f events end of physics events----//////    const Int_t size = nevents;
115        
116    for (Int_t i=0; i<numkey; i++){    TTree *RunHeadTr = (TTree*)file->Get("RunHeader");  ///run header tree
117      key = list->At(i);    pamela::EventHeader *eH=0;
118      name=(key->GetName());    pamela::RunHeaderEvent *reh=0;
119      //cout<<name<<"\n";    
120      tr = (TTree*)file->Get(name);    RunHeadTr->SetBranchAddress("Header",&eH);
121      if (tr->IsZombie()) continue;    RunHeadTr->SetBranchAddress("RunHeader",&reh);
122      nevents = tr->GetEntries();    runnevents = RunHeadTr->GetEntries();
123      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;  
     }  
       
   }  
   
   const Long64_t totalevents=totevents;  
   const Long64_t totalphysevents=totphysevents;  
     
   /////////////////////////////////////////////////  
     
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();
   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;  
     
126    
127    totphysevents=0;    //----------- If nevents < = 0 ---------------------------------/
128    totevents=0;    if (nevents<=0) {
129    for (Int_t i=0; i<numkey; i++){      printf("nevents = %i \n", nevents);
130      key = list->At(i);      printf(" \n");
131      name=(char *)(key->GetName());      TCanvas *canv = new TCanvas("No entries", "No entries ", 400, 200);
132      tr = (TTree*)file->Get(name);      canv->SetFillColor(10);
133      if (tr->IsZombie()) continue;      canv->cd();
134      tr->SetBranchAddress("Header", &eh);      TLatex *l = new TLatex();
135      nevents = tr->GetEntries();      l->SetTextAlign(12);
136      for (Int_t j = 0; j < nevents; j++){      l->SetTextSize(0.15);
137        tr->GetEntry(j);        l->SetTextColor(2);
138        ph = eh->GetPscuHeader();      noentries.str("");
139        obt[j+totevents]=ph->GetOrbitalTime() ;      noentries<< "HeaderScan_QL:";
140        pcktId[j+totevents]=ph->GetPacketId1();      l->DrawLatex(0.05, 0.7, noentries.str().c_str());
141        offset[j+totevents]=ph->GetFileOffset();      noentries.str("");
142        if(pcktId[j+totevents]==16){      noentries<< "No Physics entries for this files";
143          pcktLenght[totphysevents]=ph->GetPacketLenght();      l->DrawLatex(0.05, 0.5, noentries.str().c_str());
144          pcktCounter[totphysevents]=ph->GetCounter();      if (outDir == "./") {
145          totphysevents=totphysevents+1;        oss.str("");
146        }        oss << filename.Data() << "_HeaderScan_QL." << format.Data();
147        offset[j+totevents]=ph->GetFileOffset();      } else {
148        allCounter[j+totevents]= ph->GetCounter();        oss.str("");
149      }        oss << outDir.Data() << filename.Data() << "_HeaderScan_QL." << format.Data();
150      totevents=totevents+nevents;      }  
151        canv->Update();
152        canv->SaveAs(oss.str().c_str());
153        return;
154    }    }
155        //-------------- to know the max and min OBT ----------------------------//
156         headBr->GetEntry(0);
157    Float_t  mintime=obt[0], maxtime=obt[0], minlen=pcktLenght[0], maxlen=pcktLenght[0], mincount=pcktCounter[0], maxcount=pcktCounter[0];       ph = eh->GetPscuHeader();
158         firstime = ph->GetOrbitalTime();
159         headBr->GetEntry(nevents);
160    for (Int_t t=0; t<totalevents; t++){     ph = eh->GetPscuHeader();
161      if(obt[t]<mintime) mintime=obt[t];     lastime = ph->GetOrbitalTime();
162      if(obt[t]>maxtime) maxtime=obt[t];     int i =0;
163    }     while(i < nevents){
164    for (Int_t t=0; t<totalphysevents; t++){       headBr->GetEntry(i);
165      if(pcktLenght[t]<=minlen) minlen=pcktLenght[t];       ph = eh->GetPscuHeader();
166      if(pcktLenght[t]>=maxlen) maxlen=pcktLenght[t];       if((ph->GetOrbitalTime()) <= firstime) firstime=ph->GetOrbitalTime();
167      if(pcktCounter[t]<=mincount) mincount=pcktCounter[t];       if((ph->GetOrbitalTime()) >= lastime) lastime=ph->GetOrbitalTime();
168      if(pcktCounter[t]>=maxcount) maxcount=pcktCounter[t];       i++;
169         }
170      
171      //------------------------ First histogram -----------------------------------//
172      obmin=firstime;
173      obmax=lastime;
174      Int_t  nbin   = (lastime-firstime)/60000;
175      TH1F *h1 = new TH1F ("histo1", "" , nbin, obmin, obmax);
176    
177      //------------------- fill vectors and histogram -----------------------------//
178      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();  
187        h1->Fill(ph->GetOrbitalTime());
188        PscuCounter[k]= ph->GetCounter();
189        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    }    }
     
     
   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);  
199    
200        //----------- Graph and MultiGraph -----------------------------------------------//
201    TGraph *grPcktId1 = new TGraph(totalevents, obt, pcktId);    TMultiGraph *rate = new TMultiGraph();
202    TGraph *grPcktLenght = new TGraph(totalevents, obt, pcktLenght);    TMultiGraph *packetLength = new TMultiGraph();
203    TGraph *grPcktCounter = new TGraph(totalevents, obt, pcktCounter);    TMultiGraph *packeCounter = new TMultiGraph();
   TGraph *grFileOffset = new TGraph(totalevents, offset, allCounter);  
     
   TGaxis::SetMaxDigits(3);  
204    
   //-----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);  
205    oss1.str("");    oss1.str("");
206    oss1 << "OBT (min) with t0 = " << mintime << "ms";    oss1 << filename.Data() <<": 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    TH1F *h1 = new TH1F ("h1", oss.str().c_str(), (int)(nbintime/60), mintime, maxtime);    rate1->SetMarkerColor(kBlack);
209      rate1->SetMarkerSize(.1);
210      rate1->SetMarkerStyle(21);
211      rate->Add(rate1);
212    
213      TGraph *packetLength1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PacketLenght);
214      oss2.str("");
215      oss2 << filename.Data() <<": Lenght of Physic packets";
216      packetLength1->SetMarkerColor(2);
217      packetLength1->SetMarkerSize(.5);
218      packetLength1->SetMarkerStyle(21);
219      packetLength->Add(packetLength1);
220    
221      TGraph *packeCounter1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PscuCounter);
222      oss3.str("");
223      oss3 << filename.Data() <<": Physics Counter vs. OBT";
224      packeCounter1->SetMarkerColor(4);
225      packeCounter1->SetMarkerSize(.2);
226      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,5);
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      TPad *pad = new TPad("pad","pad", .80,.45,.90,.75);
245      pad->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    for (Int_t i=0; i<numkey; i++){    finalCanv->cd(1);
282      key = list->At(i);    stringstream ws1, ws2;
283      ws1.str("");
284      ws2.str("");
285      ws1<<"ACQ_SETTING_MODE = 1";
286      ws2<<"ACQ_SETTING_MODE = 2";
287      TPaveText *pt=0;
288      pt = new TPaveText (.60,.92,.76,.98);
289      pt->AddText(ws1.str().c_str());
290      pt->SetTextColor(3);
291      pt->SetFillColor(10);
292      pt->SetBorderSize(0);
293      pt->Draw();
294      TPaveText *pt1=0;
295      pt1 = new TPaveText (.76,.92,.92,.98);
296      pt1->AddText(ws2.str().c_str());
297      pt1->SetTextColor(4);
298      pt1->SetFillColor(10);
299      pt1->SetBorderSize(0);
300      pt1->Draw();
301      pt1 = new TPaveText (.05,.91,.6,1);
302      pt1->AddText(oss1.str().c_str());
303      pt1->SetTextColor(1);
304      pt1->SetFillColor(10);
305      pt1->SetBorderSize(0);
306      pt1->Draw();
307    
308      //----------- Second PAD -------------------------------//
309      finalCanv->cd(2);
310      all3->Draw();
311      all3->cd();
312      packetLength->Draw("AP");
313      packetLength->GetXaxis()->SetTitle("OBT    (ms)");
314      packetLength->GetXaxis()->SetTitleSize(0.05);
315      packetLength->GetXaxis()->CenterTitle();
316      packetLength->GetXaxis()->SetLabelSize(0.05);
317      packetLength->GetYaxis()->SetTitle("Lenght  (byte)");
318      packetLength->GetYaxis()->CenterTitle();
319      packetLength->GetYaxis()->SetLabelSize(0.05);
320      packetLength->GetYaxis()->SetTitleSize(0.06);
321      packetLength->GetYaxis()->SetTitleOffset(0.6);
322    
323    
324      finalCanv->cd(2);
325      pt = new TPaveText (.6,.91,.90,1);
326      pt->AddText(oss2.str().c_str());
327      pt->SetTextColor(2);
328      pt->SetFillColor(10);
329      pt->SetBorderSize(0);
330      pt->Draw();
331    
332      //----------- Third PAD -------------------------------//
333      finalCanv->cd(3);
334      all4->Draw();
335      all4->cd();
336      packeCounter->Draw("AP");
337      packeCounter->GetXaxis()->SetTitle("OBT    (ms)");
338      packeCounter->GetXaxis()->SetTitleSize(0.05);
339      packeCounter->GetXaxis()->CenterTitle();
340      packeCounter->GetXaxis()->SetLabelSize(0.05);
341      packeCounter->GetYaxis()->SetTitle("Counter");
342      packeCounter->GetYaxis()->SetTitleSize(0.05);
343      packeCounter->GetYaxis()->CenterTitle();
344      packeCounter->GetYaxis()->SetLabelSize(0.05);
345      packeCounter->GetYaxis()->SetTitleSize(0.06);
346      packeCounter->GetYaxis()->SetTitleOffset(0.6);
347    
348    
349      finalCanv->cd(3);
350      TPaveText *pt2=0;
351      pt2 = new TPaveText (.6,.91,.90,1);
352      pt2->AddText(oss3.str().c_str());
353      pt2->SetTextColor(4);
354      pt2->SetFillColor(10);
355      pt2->SetBorderSize(0);
356      pt2->Draw();  
357    
358      /**********************************************************************************************/
359    
360      TMultiGraph *mg1 = new TMultiGraph();
361      TMultiGraph *mg2 = new TMultiGraph();
362      //---------------------- fill vectors and histogram --------------------------------------------------//
363      TList *list = new TList;
364      Int_t numkey;
365      TObject *key = new TObject;
366      const char *name;
367      TTree* tr = new TTree;
368      Long64_t nevntskey=0;
369      list = file->GetListOfKeys();
370      numkey = file->GetNkeys();
371      Double_t salto;
372      for (Int_t m=0; m<numkey; m++){
373        key = list->At(m);
374      name=(char *)(key->GetName());      name=(char *)(key->GetName());
375      tr = (TTree*)file->Get(name);      tr = (TTree*)file->Get(name);
376      if (tr->IsZombie()) continue;      if (tr->IsZombie()) continue;
377      tr->SetBranchAddress("Header", &eh);      tr->SetBranchAddress("Header", &eh);
378      nevents = tr->GetEntries();      TBranch *Br = tr->GetBranch("Header");
379      for (Int_t j = 0; j < nevents; j++){      nevntskey = tr->GetEntries();
380        tr->GetEntry(j);    
381        ph = eh->GetPscuHeader();      if(nevntskey !=0){
382        if((ph->GetPacketId1()) == 16){        Int_t size1=nevntskey;
383          obtphysevents=ph->GetOrbitalTime();        Double_t *PscuCounter1 = new Double_t[size1];
384          h1->Fill(obtphysevents);        Double_t *FileOffset1 = new Double_t[size1];
385        }        Double_t *tempo1 = new Double_t[size1];
386      }  
387          int n=0;
388          while(n<nevntskey){
389            Br->GetEntry(n);
390            ph = eh->GetPscuHeader();  
391            PscuCounter1[n]= ph->GetCounter();
392            FileOffset1[n]=ph->GetFileOffset();
393            tempo1[n]=ph->GetOrbitalTime();
394            if((m==0) && (n==0)){
395              primotempo=ph->GetOrbitalTime();
396              salto=ph->GetOrbitalTime();
397            }
398            if(salto > ph->GetOrbitalTime())salto=ph->GetOrbitalTime();
399            if(ph->GetFileOffset()<= primoffset){
400              primoffset=ph->GetFileOffset();
401              primotempo=ph->GetOrbitalTime();
402            }
403            if(ph->GetFileOffset()>=ultimoffset){
404              ultimotempo=ph->GetOrbitalTime();
405              ultimoffset=ph->GetFileOffset();
406            }
407            n++;
408          }
409          TGraph *graph3= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)PscuCounter1);
410          graph3->SetMarkerColor(3);
411          graph3->SetMarkerSize(.2);
412          graph3->SetMarkerStyle(21);
413          mg1->Add(graph3);
414          
415          TGraph *graph4= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)tempo1);
416          graph4->SetMarkerColor(kBlue);
417          graph4->SetMarkerSize(.2);
418          graph4->SetMarkerStyle(21);
419          mg2->Add(graph4);
420        }
421    }    }
422      
423    h1->SetMarkerStyle(6);    TLatex *lat = new TLatex();
424    h1->GetXaxis()->SetTitle(oss1.str().c_str());    lat->SetTextAlign(12);
425    h1->GetXaxis()->CenterTitle();    lat->SetTextSize(0.15);
426    h1->GetXaxis()->SetLabelSize(0.03);    lat->SetTextColor(kBlue);
427    h1->GetYaxis()->SetTitle("number of Physic packets");  
428    h1->GetYaxis()->CenterTitle();   //------------ Fourth PAD  ---------------------//
429    h1->Draw();     finalCanv->cd(4);
430         all1->Draw();
431         all1->cd();
432    finalCanv1->cd(3);  
433    oss1.str("");     oss1.str("");
434    oss1 << "OBT (ms) with t0 = " << mintime << "ms";     oss1 << filename.Data() <<": PscuCounter vs FileOffset.";
435         mg1->Draw("AP");
436    grPcktLenght->SetTitle(oss.str().c_str());     mg1->GetXaxis()->SetTitle("File Offset");
437    grPcktLenght->GetXaxis()->SetTitle(oss1.str().c_str());     mg1->GetXaxis()->CenterTitle();
438    grPcktLenght->GetXaxis()->CenterTitle();     mg1->GetXaxis()->SetTitleOffset(0.8);
439    grPcktLenght->GetXaxis()->SetLabelSize(0.03);     mg1->GetXaxis()->SetTitleSize(0.05);
440    grPcktLenght->GetYaxis()->SetTitle("Lenght of Physic packets   (byte)");     mg1->GetXaxis()->SetLabelSize(0.05);
441    grPcktLenght->GetYaxis()->CenterTitle();     mg1->GetYaxis()->SetTitle("Counter");
442    grPcktLenght->GetYaxis()->SetLabelSize(0.03);     mg1->GetYaxis()->CenterTitle();
443    grPcktLenght->SetMarkerSize(4);     mg1->GetYaxis()->SetTitleSize(0.06);
444    grPcktLenght->Draw("AP");     mg1->GetYaxis()->SetLabelSize(0.06);
445         mg1->GetYaxis()->SetTitleOffset(0.6);
446         finalCanv->cd(4);
447    finalCanv1->cd(2);     TPaveText *pt3=0;
448    h1->Draw();     pt3 = new TPaveText (.60,.91,.90,1);
449         pt3->AddText(oss1.str().c_str());
450    finalCanv1->Update();     pt3->SetTextColor(3);
451         pt3->SetFillColor(10);
452         pt3->SetBorderSize(0);
453    ///---------canvas 2-----------------------------//     pt3->Draw();
454    finalCanv2->cd(1);  
455       //------------ Fifth PAD  ---------------------//
456    grPcktCounter->SetTitle(oss.str().c_str());     finalCanv->cd(5);
457    grPcktCounter->GetXaxis()->SetTitle(oss1.str().c_str());     all->Draw();
458    grPcktCounter->GetXaxis()->SetTitleSize(0.04);     all->cd();
   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");  
459        
460    finalCanv2->cd(2);     oss3.str("");
461         oss3 << filename.Data() <<" OBT vs FileOffset. First packet at "<<primotempo <<" ms, last packet at "<<ultimotempo<<" ms.";
462    grFileOffset->SetTitle(oss.str().c_str());     mg2->Draw("AP");
463    grFileOffset->GetXaxis()->SetTitle("File Offset");     mg2->GetXaxis()->SetTitle("File Offset");
464    grFileOffset->GetXaxis()->CenterTitle();     mg2->GetXaxis()->CenterTitle();
465    grFileOffset->GetXaxis()->SetTitleSize(0.04);     mg2->GetXaxis()->SetTitleSize(0.05);
466    grFileOffset->GetXaxis()->SetLabelSize(0.03);     mg2->GetXaxis()->SetLabelSize(0.05);
467    grFileOffset->GetYaxis()->SetTitle("Packet counter");     mg2->GetYaxis()->SetTitle("OBT");
468    grFileOffset->GetYaxis()->CenterTitle();     mg2->GetYaxis()->CenterTitle();
469    grFileOffset->GetYaxis()->SetTitleSize(0.04);     mg2->GetYaxis()->SetTitleSize(0.06);
470    grFileOffset->GetYaxis()->SetLabelSize(0.03);     mg2->GetYaxis()->SetLabelSize(0.05);
471    grFileOffset->SetMarkerSize(4);     mg2->GetYaxis()->SetTitleOffset(0.6);
472    grFileOffset->Draw("AP");     pad->Draw();
473         pad->cd();
474    finalCanv1->Update();     stringa.str("");
475         stringa << "jump at: "<<salto<<" ms \n";
476         if(salto != primotempo) lat->DrawLatex(0.08, 0.8,stringa.str().c_str());
477    //-------to save---------------------------///  
478    oss.str("");     finalCanv->cd(5);
479    oss1.str("");     TPaveText *pt4=0;
480         pt4 = new TPaveText (.38 ,.91,.92,1);
481    oss << outDir.Data() << filename.Data();     pt4->AddText(oss3.str().c_str());
482    oss1 << outDir.Data() << filename.Data();     pt4->SetTextColor(kBlue);
483       pt4->SetFillColor(10);
484       pt4->SetBorderSize(0);
485       pt4->Draw();
486    
487       finalCanv->Update();
488    
489       oss1.str("");
490       oss1 << outDir.Data() << filename.Data();
491       oss1 << "_HeaderScan"<<"." << format.Data();
492    
493    oss  << "_HeaderScan_1." << format.Data();     finalCanv->SaveAs(oss1.str().c_str());
494    oss1 << "_HeaderScan_2." << format.Data();  
495         file->Close();
   finalCanv1->SaveAs(oss.str().c_str());  
   finalCanv2->SaveAs(oss1.str().c_str());  
     
     
   file->Close();  
496        
497  }  }
498    
   
   
499  int main(int argc, char* argv[]){  int main(int argc, char* argv[]){
500    TString path;    TString path;
501    TString outDir     = "./";    TString outDir     = "./";
502    TString format     = "jpg";    TString format     = "jpg";
       
503   if (argc < 2){   if (argc < 2){
504      printf("You have to insert at least the file to analyze \n");      printf("You have to insert at least the file to analyze \n");
505      printf("Try '--help' for more information. \n");      printf("Try '--help' for more information. \n");
506      exit(1);      exit(1);
507    }      }  
   
508    if (!strcmp(argv[1], "--help")){    if (!strcmp(argv[1], "--help")){
509          printf( "Usage: HeaderScan FILE [OPTION] \n");          printf( "Usage: HeaderScan FILE [OPTION] \n");
510          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 512  int main(int argc, char* argv[]){
512          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");
513          exit(1);          exit(1);
514    }    }
   
     
515    path=argv[1];    path=argv[1];
516      for (int i = 2; i < argc; i++){
   for (int i = 2; i < argc; i++){  
     
517      if (!strcmp(argv[i], "-outDir")){      if (!strcmp(argv[i], "-outDir")){
518        if (++i >= argc){        if (++i >= argc){
519          printf( "-outDir needs arguments. \n");          printf( "-outDir needs arguments. \n");
# Line 350  int main(int argc, char* argv[]){ Line 524  int main(int argc, char* argv[]){
524          outDir = argv[i];          outDir = argv[i];
525          continue;          continue;
526        }        }
527      }      }  
       
528      if (!strcmp(argv[i], "-format")){      if (!strcmp(argv[i], "-format")){
529        if (++i >= argc){        if (++i >= argc){
530          printf( "-format needs arguments. \n");          printf( "-format needs arguments. \n");
# Line 363  int main(int argc, char* argv[]){ Line 536  int main(int argc, char* argv[]){
536          continue;          continue;
537        }        }
538      }      }
     
   
539    }    }
   
540    HeaderScan(argv[1], outDir, format);    HeaderScan(argv[1], outDir, format);
   
541  }  }
542    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.23