/[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.4 by pam-rm2, Tue Jul 11 14:27:33 2006 UTC revision 1.9 by pam-rm2, Fri Sep 8 08:26:47 2006 UTC
# 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;// offset, lastoffset, firstoffset;    ULong_t lastime, firstime, primotempo, ultimotempo, primoffset=500000000, ultimoffset;
95    double obmin=0.;    double obmin=0.;
96    double obmax=0.;    double obmax=0.;
97    stringstream oss, oss1, oss2, oss3, noentries;    stringstream oss, oss1, oss2, oss3, noentries, stringa;
98    //------- load root file --------------  
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      
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 103  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 113  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 128  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 ----------------------------//
    for (Int_t i = 0; i < nevents; i++){  
     headBr->GetEntry(i);  
     ph = eh->GetPscuHeader();  
     obt = ph->GetOrbitalTime();  
     //if(obt <= firstime) firstime=obt;  
     if(obt >= lastime) lastime=obt;  
     }  
     
   //------------------------ First histogram -----------------------------------//  
156     headBr->GetEntry(0);     headBr->GetEntry(0);
157     ph = eh->GetPscuHeader();     ph = eh->GetPscuHeader();
158     firstime = ph->GetOrbitalTime();     firstime = ph->GetOrbitalTime();
159     obmin=firstime;     headBr->GetEntry(nevents);
160     obmax=lastime;       ph = eh->GetPscuHeader();
161       lastime = ph->GetOrbitalTime();
162       int i =0;
163       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    oss1.str("");    //------------------------ First histogram -----------------------------------//
172    oss1 <<"File name= "<< filename.Data() <<". Start time= " << obmin << ", end time= "<< obmax <<" ms"<<". Physics Packet per minute";    obmin=firstime;
173      obmax=lastime;
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);
     
   //------------- Graph PscuCounter vs FileOffset ---------------------------------  
   const Int_t size = 10000;  
   Int_t ev=0;  
   Double_t PscuCounter[size], FileOffset[size], tempo[size], PcktLenght[size];  
   Int_t m=0;  
   Int_t pag= (nevents/10000)+1;  
176    
177    TMultiGraph *mg1 = new TMultiGraph();    //------------------- fill vectors and histogram -----------------------------//
178    oss1.str("");    Double_t *PscuCounter = new Double_t[size];
179    oss1 << filename.Data() <<": PscuCounter vs FileOffset.";    Double_t *PacketLenght = new Double_t[size];
180    mg1->SetTitle(oss1.str().c_str());    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      }
199    
200      //----------- Graph and MultiGraph -----------------------------------------------//
201      TMultiGraph *rate = new TMultiGraph();
202      TMultiGraph *packetLength = new TMultiGraph();
203      TMultiGraph *packeCounter = new TMultiGraph();
204    
   TMultiGraph *mg2 = new TMultiGraph();  
205    oss1.str("");    oss1.str("");
206    oss1 << filename.Data() <<": OBT vs FileOffset.";    oss1 << filename.Data() <<": Physics Packet per minute.  Start time = " << obmin << ", End time = "<< obmax <<" ms";
207    mg2->SetTitle(oss1.str().c_str());    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    TMultiGraph *mg3 = new TMultiGraph();    TGraph *packetLength1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PacketLenght);
214    oss2.str("");    oss2.str("");
215    oss2 << filename.Data() <<": Lenght of Physic packets";    oss2 << filename.Data() <<": Lenght of Physic packets";
216    mg3->SetTitle(oss2.str().c_str());    packetLength1->SetMarkerColor(2);
217      packetLength1->SetMarkerSize(.5);
218      packetLength1->SetMarkerStyle(21);
219      packetLength->Add(packetLength1);
220    
221    TMultiGraph *mg4 = new TMultiGraph();    TGraph *packeCounter1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PscuCounter);
222    oss3.str("");    oss3.str("");
223    oss3 << filename.Data() <<": Physics Counter";    oss3 << filename.Data() <<": Physics Counter vs. OBT";
224    mg4->SetTitle(oss3.str().c_str());    packeCounter1->SetMarkerColor(4);
225      packeCounter1->SetMarkerSize(.2);
226    TCanvas *finalCanv1 = new TCanvas("Header_1", base, 1280, 1024);    packeCounter1->SetMarkerStyle(21);
227    finalCanv1->Divide(1,3);    packeCounter->Add(packeCounter1);
228    finalCanv1->SetFillColor(10);  
229      //------------ Create and Draw Canvas ---------------------//
230    TCanvas *finalCanv2 = new TCanvas("Header_2", base, 1280, 1024);    TCanvas *finalCanv = new TCanvas("Header", base, 1200, 1600);
231    finalCanv2->Divide(1,2);    finalCanv->Divide(1,5);
232    finalCanv2->SetFillColor(10);    finalCanv->SetFillColor(10);
233    
234   /////ciclo sulle pagine////////////////////    TPad *all2= new TPad ("","", 0, 0, 1, 1);
235    for (Int_t interval=0; interval<pag; interval++){    all2->SetFillColor(10);
236      TPad *all3= new TPad ("","", 0, 0, 1, 1);
237     for (Int_t l = 0; l < size; l++){     //ciclo su un sottogruppo    all3->SetFillColor(10);
238      headBr->GetEntry(ev+l);    TPad *all4= new TPad ("","", 0, 0, 1, 1);
239      ph = eh->GetPscuHeader();      all4->SetFillColor(10);
240      PscuCounter[l]= ph->GetCounter(); // = PcktCounter[l]    TPad *all= new TPad ("","", 0, 0, 1, 1);
241      FileOffset[l]=ph->GetFileOffset();    all->SetFillColor(10);
242      tempo[l]=ph->GetOrbitalTime();    TPad *all1= new TPad ("","", 0, 0, 1, 1);
243      PcktLenght[l]=ph->GetPacketLenght();    all1->SetFillColor(10);
244      h1->Fill(ph->GetOrbitalTime());      TPad *pad = new TPad("pad","pad", .80,.45,.90,.75);
245      m=l;    pad->SetFillColor(10);
246     }    
247      TLine li;
248     TGraph *graph1= new TGraph(m, tempo, PcktLenght);    li.SetLineStyle(4);
249     graph1->SetMarkerColor(2);    li.SetLineWidth(2);
250     graph1->SetMarkerSize(.5);  
251     graph1->SetMarkerStyle(21);    //----------- First PAD -------------------------------//
252     mg3->Add(graph1);    finalCanv->cd(1);
253      all2->Draw();
254     TGraph *graph2= new TGraph(m, tempo, PscuCounter);    all2->cd();
255     graph2->SetMarkerColor(4);    rate->Draw("ALP");
256     graph2->SetMarkerSize(.5);    rate->GetXaxis()->SetTitle("OBT   (ms)");
257     graph2->SetMarkerStyle(21);    rate->GetXaxis()->SetTitleSize(0.05);
258     mg4->Add(graph2);    rate->GetXaxis()->CenterTitle();
259      rate->GetXaxis()->SetLabelSize(0.05);
260     TGraph *graph3= new TGraph(m, FileOffset, PscuCounter);    rate->GetYaxis()->SetTitle("Number of events ");
261     graph3->SetMarkerColor(3);    rate->GetYaxis()->CenterTitle();
262     graph3->SetMarkerSize(.5);    rate->GetYaxis()->SetLabelSize(0.05);
263     graph3->SetMarkerStyle(21);    rate->GetYaxis()->SetTitleSize(0.06);
264     mg1->Add(graph3);    rate->GetYaxis()->SetTitleOffset(0.6);
265      for (Int_t l = 0; l < runnevents; l++){
266     TGraph *graph4= new TGraph(m, FileOffset, tempo);      RunHeadTr->GetEntry(l);
267     graph4->SetMarkerColor(kBlue);      ph = eH->GetPscuHeader();
268     graph4->SetMarkerSize(.5);      int ws= reh->RM_ACQ_SETTING_MODE;
269     graph4->SetMarkerStyle(21);      int id = ph->GetPacketId1();
270     mg2->Add(graph4);      Int_t obt = ph->GetOrbitalTime();
271        if (ws==1){
272     ev=ev+10000;        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    }    }
   //------------ First Canvas ---------------------//  
   
   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();  
280        
281    finalCanv1->cd(2);    finalCanv->cd(1);
282    mg3->Draw("AP");    stringstream ws1, ws2;
283    mg3->GetXaxis()->SetTitle("OBT    (ms)");    ws1.str("");
284    mg3->GetXaxis()->CenterTitle();    ws2.str("");
285    mg3->GetXaxis()->SetLabelSize(0.04);    ws1<<"ACQ_SETTING_MODE = 1";
286    mg3->GetYaxis()->SetTitle("Lenght  (byte)");    ws2<<"ACQ_SETTING_MODE = 2";
287    mg3->GetYaxis()->CenterTitle();    TPaveText *pt=0;
288    mg3->GetYaxis()->SetLabelSize(0.04);    pt = new TPaveText (.60,.92,.76,.98);
289    mg3->GetYaxis()->SetTitleSize(0.06);    pt->AddText(ws1.str().c_str());
290    mg3->GetYaxis()->SetTitleOffset(0.8);    pt->SetTextColor(3);
291      pt->SetFillColor(10);
292      pt->SetBorderSize(0);
293    finalCanv1->cd(3);    pt->Draw();
294    mg4->Draw("AP");    TPaveText *pt1=0;
295    mg4->GetXaxis()->SetTitle("OBT    (ms)");    pt1 = new TPaveText (.76,.92,.92,.98);
296    mg4->GetXaxis()->SetTitleSize(0.05);    pt1->AddText(ws2.str().c_str());
297    mg4->GetXaxis()->CenterTitle();    pt1->SetTextColor(4);
298    mg4->GetXaxis()->SetLabelSize(0.04);    pt1->SetFillColor(10);
299    mg4->GetYaxis()->SetTitle("Counter");    pt1->SetBorderSize(0);
300    mg4->GetYaxis()->SetTitleSize(0.05);    pt1->Draw();
301    mg4->GetYaxis()->CenterTitle();    pt1 = new TPaveText (.05,.91,.6,1);
302    mg4->GetYaxis()->SetLabelSize(0.04);    pt1->AddText(oss1.str().c_str());
303    mg4->GetYaxis()->SetTitleSize(0.06);    pt1->SetTextColor(1);
304    mg4->GetYaxis()->SetTitleOffset(0.8);    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    finalCanv1->Update();    TMultiGraph *mg1 = new TMultiGraph();
361        TMultiGraph *mg2 = new TMultiGraph();
362    oss.str("");    //---------------------- fill vectors and histogram --------------------------------------------------//
363    oss << outDir.Data() << filename.Data();    TList *list = new TList;
364    oss  << "_HeaderScan_1." << format.Data();    Int_t numkey;
365    finalCanv1->SaveAs(oss.str().c_str());    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());
375        tr = (TTree*)file->Get(name);
376        if (tr->IsZombie()) continue;
377        tr->SetBranchAddress("Header", &eh);
378        TBranch *Br = tr->GetBranch("Header");
379        nevntskey = tr->GetEntries();
380    
381        if(nevntskey !=0){
382          Int_t size1=nevntskey;
383          Double_t *PscuCounter1 = new Double_t[size1];
384          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      TLatex *lat = new TLatex();
424      lat->SetTextAlign(12);
425      lat->SetTextSize(0.15);
426      lat->SetTextColor(kBlue);
427    
428     //------------ Fourth PAD  ---------------------//
429       finalCanv->cd(4);
430       all1->Draw();
431       all1->cd();
432    
433    //------------ Second Canvas ---------------------//     oss1.str("");
434     finalCanv2->cd(1);     oss1 << filename.Data() <<": PscuCounter vs FileOffset.";
435     mg1->Draw("AP");     mg1->Draw("AP");
436     mg1->GetXaxis()->SetTitle("File Offset");     mg1->GetXaxis()->SetTitle("File Offset");
437     mg1->GetXaxis()->CenterTitle();     mg1->GetXaxis()->CenterTitle();
438       mg1->GetXaxis()->SetTitleOffset(0.8);
439     mg1->GetXaxis()->SetTitleSize(0.05);     mg1->GetXaxis()->SetTitleSize(0.05);
440     mg1->GetXaxis()->SetLabelSize(0.04);     mg1->GetXaxis()->SetLabelSize(0.05);
441     mg1->GetYaxis()->SetTitle("Counter");     mg1->GetYaxis()->SetTitle("Counter");
442     mg1->GetYaxis()->CenterTitle();     mg1->GetYaxis()->CenterTitle();
443     mg1->GetYaxis()->SetTitleSize(0.05);     mg1->GetYaxis()->SetTitleSize(0.06);
444     mg1->GetYaxis()->SetLabelSize(0.04);     mg1->GetYaxis()->SetLabelSize(0.06);
445       mg1->GetYaxis()->SetTitleOffset(0.6);
446     finalCanv2->cd(2);     finalCanv->cd(4);
447       TPaveText *pt3=0;
448       pt3 = new TPaveText (.60,.91,.90,1);
449       pt3->AddText(oss1.str().c_str());
450       pt3->SetTextColor(3);
451       pt3->SetFillColor(10);
452       pt3->SetBorderSize(0);
453       pt3->Draw();
454    
455     //------------ Fifth PAD  ---------------------//
456       finalCanv->cd(5);
457       all->Draw();
458       all->cd();
459      
460       oss3.str("");
461       oss3 << filename.Data() <<" OBT vs FileOffset. First packet at "<<primotempo <<" ms, last packet at "<<ultimotempo<<" ms.";
462     mg2->Draw("AP");     mg2->Draw("AP");
463     mg2->GetXaxis()->SetTitle("File Offset");     mg2->GetXaxis()->SetTitle("File Offset");
464     mg2->GetXaxis()->CenterTitle();     mg2->GetXaxis()->CenterTitle();
465     mg2->GetXaxis()->SetTitleSize(0.05);     mg2->GetXaxis()->SetTitleSize(0.05);
466     mg2->GetXaxis()->SetLabelSize(0.04);     mg2->GetXaxis()->SetLabelSize(0.05);
467     mg2->GetYaxis()->SetTitle("OBT");     mg2->GetYaxis()->SetTitle("OBT");
468     mg2->GetYaxis()->CenterTitle();     mg2->GetYaxis()->CenterTitle();
469     mg2->GetYaxis()->SetTitleSize(0.05);     mg2->GetYaxis()->SetTitleSize(0.06);
470     mg2->GetYaxis()->SetLabelSize(0.04);     mg2->GetYaxis()->SetLabelSize(0.05);
471       mg2->GetYaxis()->SetTitleOffset(0.6);
472       pad->Draw();
473       pad->cd();
474       stringa.str("");
475       stringa << "jump at: "<<salto<<" ms \n";
476       if(salto != primotempo) lat->DrawLatex(0.08, 0.8,stringa.str().c_str());
477    
478       finalCanv->cd(5);
479       TPaveText *pt4=0;
480       pt4 = new TPaveText (.38 ,.91,.92,1);
481       pt4->AddText(oss3.str().c_str());
482       pt4->SetTextColor(kBlue);
483       pt4->SetFillColor(10);
484       pt4->SetBorderSize(0);
485       pt4->Draw();
486    
487     finalCanv2->Update();     finalCanv->Update();
488    
489     oss1.str("");     oss1.str("");
490     oss1 << outDir.Data() << filename.Data();     oss1 << outDir.Data() << filename.Data();
491     oss1 << "_HeaderScan_2"<<"." << format.Data();     oss1 << "_HeaderScan"<<"." << format.Data();
492    
493       finalCanv->SaveAs(oss1.str().c_str());
494    
    finalCanv2->SaveAs(oss1.str().c_str());  
     
495     file->Close();     file->Close();
496        
497  }  }
# Line 364  int main(int argc, char* argv[]){ Line 539  int main(int argc, char* argv[]){
539    }    }
540    HeaderScan(argv[1], outDir, format);    HeaderScan(argv[1], outDir, format);
541  }  }
542    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.23