/[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.8 by pam-rm2, Thu Aug 31 14:28:32 2006 UTC revision 1.14 by pam-rm2, Fri Mar 9 10:43:40 2007 UTC
# Line 59  Line 59 
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  * 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    
# Line 70  Line 75 
75  #include "TMultiGraph.h"  #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 87  Line 94 
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){  
 {  
98    
     
99    //------------------- Variables initilization -------------------------//    //------------------- Variables initilization -------------------------//
100    Long64_t nevents=0, runnevents;    Long64_t nevents=0, runnevents;
101    ULong_t lastime, firstime, primotempo, ultimotempo, primoffset=500000000, ultimoffset;    ULong_t lastime, firstime, primotempo, ultimotempo, primoffset=500000000, ultimoffset;
102    double obmin=0., time=0.;    double obmin=0.;
103    double obmax=0.;    double obmax=0.;
104    stringstream oss, oss1, oss2, oss3, noentries, stringa;    stringstream oss, oss1, oss2, oss3, noentries, stringa;
105    //------- load root file --------------  
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      
113    TTree *PhysicsTr = (TTree*)file->Get("Physics");    TTree *PhysicsTr = (TTree*)file->Get("Physics");
114    TBranch *headBr = PhysicsTr->GetBranch("Header");    TBranch *headBr = PhysicsTr->GetBranch("Header");
115    
116    pamela::EventHeader  *eh   = 0;    pamela::EventHeader  *eh   = new pamela::EventHeader;
117    pamela::PscuHeader   *ph   = 0;    pamela::PscuHeader   *ph   = new pamela::PscuHeader;
118    
119    PhysicsTr->SetBranchAddress("Header", &eh);    PhysicsTr->SetBranchAddress("Header", &eh);
120    nevents = headBr->GetEntries();    nevents = PhysicsTr->GetEntries();
121        const Int_t size = nevents;
122    
123    TTree *RunHeadTr = (TTree*)file->Get("RunHeader");  ///run header tree    TTree *RunHeadTr = (TTree*)file->Get("RunHeader");  ///run header tree
124    pamela::EventHeader *eH=0;    pamela::EventHeader *eH= new pamela::EventHeader;
125    pamela::RunHeaderEvent *reh=0;    pamela::RunHeaderEvent *reh=new pamela::RunHeaderEvent;
126        
127    RunHeadTr->SetBranchAddress("Header",&eH);    RunHeadTr->SetBranchAddress("Header",&eH);
128    RunHeadTr->SetBranchAddress("RunHeader",&reh);    RunHeadTr->SetBranchAddress("RunHeader",&reh);
129    runnevents = RunHeadTr->GetEntries();    runnevents = RunHeadTr->GetEntries();
   
130    
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();
# Line 129  void HeaderScan(TString base, TString ou Line 135  void HeaderScan(TString base, TString ou
135    if (nevents<=0) {    if (nevents<=0) {
136      printf("nevents = %i \n", nevents);      printf("nevents = %i \n", nevents);
137      printf(" \n");      printf(" \n");
   
138      TCanvas *canv = new TCanvas("No entries", "No entries ", 400, 200);      TCanvas *canv = new TCanvas("No entries", "No entries ", 400, 200);
139      canv->SetFillColor(10);      canv->SetFillColor(10);
140      canv->cd();      canv->cd();
   
141      TLatex *l = new TLatex();      TLatex *l = new TLatex();
142      l->SetTextAlign(12);      l->SetTextAlign(12);
143      l->SetTextSize(0.15);      l->SetTextSize(0.15);
# Line 144  void HeaderScan(TString base, TString ou Line 148  void HeaderScan(TString base, TString ou
148      noentries.str("");      noentries.str("");
149      noentries<< "No Physics entries for this files";      noentries<< "No Physics entries for this files";
150      l->DrawLatex(0.05, 0.5, noentries.str().c_str());      l->DrawLatex(0.05, 0.5, noentries.str().c_str());
   
151      if (outDir == "./") {      if (outDir == "./") {
152        oss.str("");        oss.str("");
153        oss << filename.Data() << "_HeaderScan_QL." << format.Data();        oss << filename.Data() << "_HeaderScan_QL." << format.Data();
154      } else {      } else {
155        oss.str("");        oss.str("");
156        oss << outDir.Data() << filename.Data() << "_HeaderScan_QL." << format.Data();        oss << outDir.Data() << filename.Data() << "_HeaderScan_QL." << format.Data();
157      }      }  
     
158      canv->Update();      canv->Update();
159      canv->SaveAs(oss.str().c_str());      canv->SaveAs(oss.str().c_str());
   
160      return;      return;
161    }    }
   
162    //-------------- to know the max and min OBT ----------------------------//    //-------------- to know the max and min OBT ----------------------------//
163     headBr->GetEntry(0);     headBr->GetEntry(0);
164     ph = eh->GetPscuHeader();     ph = eh->GetPscuHeader();
165     firstime = ph->GetOrbitalTime();     firstime = ph->GetOrbitalTime();
166       headBr->GetEntry(nevents);
167     int iii =0;     ph = eh->GetPscuHeader();
168     while(iii < nevents){     lastime = ph->GetOrbitalTime();
169       headBr->GetEntry(iii);     int i =0;
170       while(i < nevents){
171         headBr->GetEntry(i);
172       ph = eh->GetPscuHeader();       ph = eh->GetPscuHeader();
173       if((ph->GetOrbitalTime()) <= firstime) firstime=ph->GetOrbitalTime();       if((ph->GetOrbitalTime()) <= firstime) firstime=ph->GetOrbitalTime();
174       if((ph->GetOrbitalTime()) >= lastime) lastime=ph->GetOrbitalTime();       if((ph->GetOrbitalTime()) >= lastime) lastime=ph->GetOrbitalTime();
175       iii++;       i++;
176       }       }
177        
178    //------------------------ First histogram -----------------------------------//    //------------------------ First histogram -----------------------------------//
     
179    obmin=firstime;    obmin=firstime;
180    obmax=lastime;      obmax=lastime;
     
   oss1.str("");  
   oss1 << filename.Data() <<": Physics Packet per minute.  Start time = " << obmin << ", End time = "<< obmax <<" ms";  
181    Int_t  nbin   = (lastime-firstime)/60000;    Int_t  nbin   = (lastime-firstime)/60000;
   //TH1F *h1 = new TH1F ("histo1", oss1.str().c_str(), nbin, obmin, obmax);  
182    TH1F *h1 = new TH1F ("histo1", "" , nbin, obmin, obmax);    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      }
208    
209      //----------- Graph and MultiGraph -----------------------------------------------//
210      TMultiGraph *rate = new TMultiGraph();
211      TMultiGraph *packetLength = new TMultiGraph();
212      TMultiGraph *packeCounter = new TMultiGraph();
213    
214      oss1.str("");
215      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      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("");    oss2.str("");
224    oss2 << filename.Data() <<": Lenght of Physic packets";    oss2 <<"Lenght of Physic packets";
225    Int_t  nbin2  =(lastime-firstime);    packetLength1->SetMarkerColor(2);
226    //TH1F *packetLength = new TH1F ("packetLength", oss2.str().c_str(), nbin2, obmin, obmax);    packetLength1->SetMarkerSize(.3);
227    TH1F *packetLength = new TH1F ("packetLength", "", nbin2, obmin, obmax);    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("");    oss3.str("");
232    oss3 << filename.Data() <<": Physics Counter vs. OBT";    oss3 <<"Physics Counter vs. OBT";
233    //  TH1F *packeCounter = new TH1F ("packeCounter", oss3.str().c_str(), nbin2, obmin, obmax);    packeCounter1->SetMarkerColor(4);
234    TH1F *packeCounter = new TH1F ("packeCounter", "", nbin2, obmin, obmax);    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);    TCanvas *finalCanv = new TCanvas("Header", base, 1200, 1600);
240    finalCanv->Divide(1,5);    finalCanv->Divide(1,6);
241    finalCanv->SetFillColor(10);    finalCanv->SetFillColor(10);
242      
243    TPad *all2= new TPad ("","", 0, 0, 1, 1);    TPad *all2= new TPad ("","", 0, 0, 1, 1);
244    all2->SetFillColor(10);    all2->SetFillColor(10);
245    TPad *all3= new TPad ("","", 0, 0, 1, 1);    TPad *all3= new TPad ("","", 0, 0, 1, 1);
# Line 209  void HeaderScan(TString base, TString ou Line 250  void HeaderScan(TString base, TString ou
250    all->SetFillColor(10);    all->SetFillColor(10);
251    TPad *all1= new TPad ("","", 0, 0, 1, 1);    TPad *all1= new TPad ("","", 0, 0, 1, 1);
252    all1->SetFillColor(10);    all1->SetFillColor(10);
   TPad *pad = new TPad("pad","pad", .80,.50,.90,.80);  
   pad->SetFillColor(10);  
   //----------------------------------------------------  
   
   TMultiGraph *mg1 = new TMultiGraph();  
   
   
   TMultiGraph *mg2 = new TMultiGraph();  
253    
   
   //--------------------------------------------------------------------------  
   for (Int_t l = 0; l < nevents; l++){  
     headBr->GetEntry(l);  
     ph = eh->GetPscuHeader();    
     h1->Fill(ph->GetOrbitalTime());  
     packetLength->Fill(ph->GetOrbitalTime(),ph->GetPacketLenght());  
     packeCounter->Fill(ph->GetOrbitalTime(),ph->GetCounter());  
   }    
   
   //------------ First Canvas ---------------------//  
254    TLine li;    TLine li;
255    li.SetLineStyle(4);    li.SetLineStyle(4);
256    li.SetLineWidth(2);    li.SetLineWidth(2);
257      
258      //----------- First PAD -------------------------------//
259    finalCanv->cd(1);    finalCanv->cd(1);
260    all2->Draw();    all2->Draw();
261    all2->cd();    all2->cd();
262      rate->Draw("ALP");
263    h1->SetStats(kFALSE);    rate->GetXaxis()->SetTitle("OBT   (ms)");
264    h1->GetXaxis()->SetTitle("OBT   (ms)");    rate->GetXaxis()->SetTitleSize(0.05);
265    h1->GetXaxis()->CenterTitle();    rate->GetXaxis()->CenterTitle();
266    h1->GetXaxis()->SetLabelSize(0.04);    rate->GetXaxis()->SetLabelSize(0.05);
267    h1->GetYaxis()->SetTitle("Number of events ");    rate->GetYaxis()->SetTitle("Number of events ");
268    h1->GetYaxis()->CenterTitle();    rate->GetYaxis()->CenterTitle();
269    h1->GetYaxis()->SetLabelSize(0.06);    rate->GetYaxis()->SetLabelSize(0.05);
270    h1->GetYaxis()->SetTitleSize(0.06);    rate->GetYaxis()->SetTitleSize(0.06);
271    h1->GetYaxis()->SetTitleOffset(0.6);    rate->GetYaxis()->SetTitleOffset(0.6);
   h1->Draw();  
272    for (Int_t l = 0; l < runnevents; l++){    for (Int_t l = 0; l < runnevents; l++){
273      RunHeadTr->GetEntry(l);      RunHeadTr->GetEntry(l);
274      ph = eH->GetPscuHeader();      ph = eH->GetPscuHeader();
275      int ws= reh->RM_ACQ_SETTING_MODE;      int ws= reh->RM_ACQ_SETTING_MODE;
276      int id = ph->GetPacketId1();      int id = ph->GetPacketId1();
277      Int_t obt = ph->GetOrbitalTime();      ULong_t utile1=ph->GetOrbitalTime();
278        Double_t obt = utile1;
279      if (ws==1){      if (ws==1){
280        li.SetLineColor(3);        li.SetLineColor(3);
281        li.DrawLine(obt,0,obt,h1->GetMaximum());        li.DrawLine(obt,0,obt,max);
282      }      }
283      else if (ws==2){      else if (ws==2){
284        li.SetLineColor(4);        li.SetLineColor(4);
285        li.DrawLine(obt,0,obt,h1->GetMaximum());        li.DrawLine(obt,0,obt,max);
286      }      }  
287    }      }
288      
289      RunHeadTr->GetEntry(0);
290      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);    finalCanv->cd(1);
297    stringstream ws1, ws2;    stringstream ws1, ws2;
298    ws1.str("");    ws1.str("");
# Line 283  void HeaderScan(TString base, TString ou Line 313  void HeaderScan(TString base, TString ou
313    pt1->SetFillColor(10);    pt1->SetFillColor(10);
314    pt1->SetBorderSize(0);    pt1->SetBorderSize(0);
315    pt1->Draw();    pt1->Draw();
   //  TPaveText *pt1=0;  
316    pt1 = new TPaveText (.05,.91,.6,1);    pt1 = new TPaveText (.05,.91,.6,1);
317    pt1->AddText(oss1.str().c_str());    pt1->AddText(oss1.str().c_str());
318    pt1->SetTextColor(1);    pt1->SetTextColor(1);
# Line 291  void HeaderScan(TString base, TString ou Line 320  void HeaderScan(TString base, TString ou
320    pt1->SetBorderSize(0);    pt1->SetBorderSize(0);
321    pt1->Draw();    pt1->Draw();
322    
323      //----------- Second PAD -------------------------------//
324    finalCanv->cd(2);    finalCanv->cd(2);
325    all3->Draw();    all3->Draw();
326    all3->cd();    all3->cd();
327    packetLength->SetStats(kFALSE);    packetLength->Draw("AP");
   packetLength->SetMarkerColor(2);  
   packetLength->SetMarkerSize(.5);  
   packetLength->SetMarkerStyle(21);  
328    packetLength->GetXaxis()->SetTitle("OBT    (ms)");    packetLength->GetXaxis()->SetTitle("OBT    (ms)");
329      packetLength->GetXaxis()->SetTitleSize(0.05);
330    packetLength->GetXaxis()->CenterTitle();    packetLength->GetXaxis()->CenterTitle();
331    packetLength->GetXaxis()->SetLabelSize(0.05);    packetLength->GetXaxis()->SetLabelSize(0.05);
332    packetLength->GetYaxis()->SetTitle("Lenght  (byte)");    packetLength->GetYaxis()->SetTitle("Lenght  (byte)");
# Line 306  void HeaderScan(TString base, TString ou Line 334  void HeaderScan(TString base, TString ou
334    packetLength->GetYaxis()->SetLabelSize(0.05);    packetLength->GetYaxis()->SetLabelSize(0.05);
335    packetLength->GetYaxis()->SetTitleSize(0.06);    packetLength->GetYaxis()->SetTitleSize(0.06);
336    packetLength->GetYaxis()->SetTitleOffset(0.6);    packetLength->GetYaxis()->SetTitleOffset(0.6);
337    packetLength->Draw("9p");  
338    
339    finalCanv->cd(2);    finalCanv->cd(2);
   //TPaveText *pt=0;  
340    pt = new TPaveText (.6,.91,.90,1);    pt = new TPaveText (.6,.91,.90,1);
341    pt->AddText(oss2.str().c_str());    pt->AddText(oss2.str().c_str());
342    pt->SetTextColor(2);    pt->SetTextColor(2);
# Line 316  void HeaderScan(TString base, TString ou Line 344  void HeaderScan(TString base, TString ou
344    pt->SetBorderSize(0);    pt->SetBorderSize(0);
345    pt->Draw();    pt->Draw();
346    
347      //----------- Third PAD -------------------------------//
348    finalCanv->cd(3);    finalCanv->cd(3);
349    all4->Draw();    all4->Draw();
350    all4->cd();    all4->cd();
351    packeCounter->GetTitle();    packeCounter->Draw("AP");
   packeCounter->SetStats(kFALSE);  
   packeCounter->SetMarkerColor(4);  
   packeCounter->SetMarkerSize(.2);  
   packeCounter->SetMarkerStyle(21);  
352    packeCounter->GetXaxis()->SetTitle("OBT    (ms)");    packeCounter->GetXaxis()->SetTitle("OBT    (ms)");
353    packeCounter->GetXaxis()->SetTitleSize(0.05);    packeCounter->GetXaxis()->SetTitleSize(0.05);
354    packeCounter->GetXaxis()->CenterTitle();    packeCounter->GetXaxis()->CenterTitle();
# Line 334  void HeaderScan(TString base, TString ou Line 359  void HeaderScan(TString base, TString ou
359    packeCounter->GetYaxis()->SetLabelSize(0.05);    packeCounter->GetYaxis()->SetLabelSize(0.05);
360    packeCounter->GetYaxis()->SetTitleSize(0.06);    packeCounter->GetYaxis()->SetTitleSize(0.06);
361    packeCounter->GetYaxis()->SetTitleOffset(0.6);    packeCounter->GetYaxis()->SetTitleOffset(0.6);
362     Double_t min = 500000000.;  
363     for (Int_t l = 0; l < nevents; l++){  
      if((packeCounter->GetBinContent(l))<=  min && (packeCounter->GetBinContent(l))!= 0) min = packeCounter->GetBinContent(l);  
      }    
   packeCounter->SetMinimum(min);  
   packeCounter->Draw("9p");  
364    finalCanv->cd(3);    finalCanv->cd(3);
365    TPaveText *pt2=0;    TPaveText *pt2=0;
366    pt2 = new TPaveText (.6,.91,.90,1);    pt2 = new TPaveText (.6,.91,.90,1);
# Line 347  void HeaderScan(TString base, TString ou Line 368  void HeaderScan(TString base, TString ou
368    pt2->SetTextColor(4);    pt2->SetTextColor(4);
369    pt2->SetFillColor(10);    pt2->SetFillColor(10);
370    pt2->SetBorderSize(0);    pt2->SetBorderSize(0);
371    pt2->Draw();    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;    TList *list = new TList;
379    Int_t numkey;    Int_t numkey;
380    TObject *key = new TObject;    TObject *key = new TObject;
381    const char *name;    const char *name;
382      char *SoftInfo="SoftInfo";
383    TTree* tr = new TTree;    TTree* tr = new TTree;
384    Long64_t nevntskey=0;    Long64_t nevntskey=0;
385    list = file->GetListOfKeys();    list = file->GetListOfKeys();
386    numkey = file->GetNkeys();    numkey = file->GetNkeys();
387    int ll=0;    ULong_t salto;
388    int kk=0;    for (Int_t m=0; m<numkey; m++){
389        key = list->At(m);
   for (Int_t i=0; i<numkey; i++){  
     key = list->At(i);  
390      name=(char *)(key->GetName());      name=(char *)(key->GetName());
391        if(strcmp(name,SoftInfo)==0)continue;
392      tr = (TTree*)file->Get(name);      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      TBranch *Br = tr->GetBranch("Header");      TBranch *Br = tr->GetBranch("Header");
397      nevntskey = tr->GetEntries();      nevntskey = tr->GetEntries();
   
398      if(nevntskey !=0){      if(nevntskey !=0){
399        Int_t size1=nevntskey;        Int_t size1=nevntskey;
400        Double_t *PscuCounter1 = new Double_t[size1];        Double_t *PscuCounter1 = new Double_t[size1];
401        Double_t *FileOffset1 = new Double_t[size1];        Double_t *FileOffset1 = new Double_t[size1];
402        Double_t *tempo1 = new Double_t[size1];        Double_t *tempo1 = new Double_t[size1];
403        //Double_t *salto= new Double_t[100];  
404        //Double_t *salto1= new Double_t[100];        int n=0;
405                while(n<nevntskey){
406        int l=0;          Br->GetEntry(n);
       while(l<nevntskey){  
         Br->GetEntry(l);  
407          ph = eh->GetPscuHeader();            ph = eh->GetPscuHeader();  
408          PscuCounter1[l]= ph->GetCounter();          PscuCounter1[n]= ph->GetCounter();
409          FileOffset1[l]=ph->GetFileOffset();          FileOffset1[n]=ph->GetFileOffset();
410          tempo1[l]=ph->GetOrbitalTime();          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){          if(ph->GetFileOffset()<= primoffset){
418            primoffset =ph->GetFileOffset();            primoffset=ph->GetFileOffset();
419            primotempo=ph->GetOrbitalTime();            primotempo=ph->GetOrbitalTime();
420          }          }
421          if(ph->GetFileOffset()>= ultimoffset){          if(ph->GetFileOffset()>=ultimoffset){
422            ultimoffset = ph->GetFileOffset();            ultimotempo=ph->GetOrbitalTime();
423            ultimotempo = ph->GetOrbitalTime();            ultimoffset=ph->GetFileOffset();
424            }          }
425          /*if(l>0){          n++;
           if(tempo1[l] < tempo1[l-1]){  
             salto1[kk] =(const Double_t*)tempo1[l-1];  
             salto[ll]= ph->GetOrbitalTime();  
             ll=ll+1;  
             kk=kk+1;  
           }  
           }*/  
         l++;  
426        }        }
427          
428        TGraph *graph3= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)PscuCounter1);        TGraph *graph3= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)PscuCounter1);
429        graph3->SetMarkerColor(3);        graph3->SetMarkerColor(3);
430        graph3->SetMarkerSize(.2);        graph3->SetMarkerSize(.2);
# Line 418  void HeaderScan(TString base, TString ou Line 438  void HeaderScan(TString base, TString ou
438        mg2->Add(graph4);        mg2->Add(graph4);
439      }      }
440    }    }
   
   kk=kk-1;  
   ll=ll-1;  
441    
442    TLatex *lat = new TLatex();    TLatex *lat = new TLatex();
443    lat->SetTextAlign(12);    lat->SetTextAlign(12);
444    lat->SetTextSize(0.15);    lat->SetTextSize(0.15);
445    lat->SetTextColor(kBlue);    lat->SetTextColor(kBlue);
446   //------------ Second Canvas ---------------------//  
447     //------------ Fourth PAD  ---------------------//
448     finalCanv->cd(4);     finalCanv->cd(4);
449     all1->Draw();     all1->Draw();
450     all1->cd();     all1->cd();
451    
452     oss1.str("");     oss1.str("");
453     oss1 << filename.Data() <<": PscuCounter vs FileOffset.";     oss1 <<"PscuCounter vs FileOffset.";
    //mg1->SetTitle(oss1.str().c_str());  
454     mg1->Draw("AP");     mg1->Draw("AP");
455     mg1->GetXaxis()->SetTitle("File Offset");     mg1->GetXaxis()->SetTitle("File Offset");
456     mg1->GetXaxis()->CenterTitle();     mg1->GetXaxis()->CenterTitle();
457     mg1->GetXaxis()->SetTitleSize(0.045);     mg1->GetXaxis()->SetTitleOffset(0.8);
458       mg1->GetXaxis()->SetTitleSize(0.05);
459     mg1->GetXaxis()->SetLabelSize(0.05);     mg1->GetXaxis()->SetLabelSize(0.05);
460     mg1->GetYaxis()->SetTitle("Counter");     mg1->GetYaxis()->SetTitle("Counter");
461     mg1->GetYaxis()->CenterTitle();     mg1->GetYaxis()->CenterTitle();
462     mg1->GetYaxis()->SetTitleSize(0.05);     mg1->GetYaxis()->SetTitleSize(0.06);
463     mg1->GetYaxis()->SetLabelSize(0.06);     mg1->GetYaxis()->SetLabelSize(0.06);
464     mg1->GetYaxis()->SetTitleOffset(0.7);     mg1->GetYaxis()->SetTitleOffset(0.6);
465     finalCanv->cd(4);     finalCanv->cd(4);
466     TPaveText *pt3=0;     TPaveText *pt3=0;
467     pt3 = new TPaveText (.60,.91,.90,1);     pt3 = new TPaveText (.60,.91,.90,1);
# Line 453  void HeaderScan(TString base, TString ou Line 471  void HeaderScan(TString base, TString ou
471     pt3->SetBorderSize(0);     pt3->SetBorderSize(0);
472     pt3->Draw();     pt3->Draw();
473    
474     //------------ Fifth PAD  ---------------------//
475     finalCanv->cd(5);     finalCanv->cd(5);
476     all->Draw();     all->Draw();
477     all->cd();     all->cd();
     
478     oss3.str("");     oss3.str("");
479     oss3 << filename.Data() <<" OBT vs FileOffset. First packet = "<<primotempo <<" ms,  Last packet = "<<ultimotempo<<" ms.";     oss3 << "OBT vs FileOffset";
    //mg2->SetTitle(oss3.str().c_str());  
480     mg2->Draw("AP");     mg2->Draw("AP");
481     mg2->GetXaxis()->SetTitle("File Offset");     mg2->GetXaxis()->SetTitle("File Offset");
482     mg2->GetXaxis()->CenterTitle();     mg2->GetXaxis()->CenterTitle();
483     mg2->GetXaxis()->SetTitleSize(0.045);     mg2->GetXaxis()->SetTitleSize(0.05);
484     mg2->GetXaxis()->SetLabelSize(0.05);     mg2->GetXaxis()->SetLabelSize(0.05);
485     mg2->GetYaxis()->SetTitle("OBT");     mg2->GetYaxis()->SetTitle("OBT");
486     mg2->GetYaxis()->CenterTitle();     mg2->GetYaxis()->CenterTitle();
487     mg2->GetYaxis()->SetTitleSize(0.05);     mg2->GetYaxis()->SetTitleSize(0.06);
488     mg2->GetYaxis()->SetLabelSize(0.055);     mg2->GetYaxis()->SetLabelSize(0.05);
489     mg2->GetYaxis()->SetTitleOffset(0.6);     mg2->GetYaxis()->SetTitleOffset(0.6);
    /*pad->Draw();  
    pad->cd();  
    stringa.str("");  
    stringa << "new data from "<< primotempo;//<<" to "<< salto1;  
    lat->DrawLatex(0.08, 0.8, stringa.str().c_str());  
    double min2 = 0.8;  
    if(ll > 0){  
    pad->Draw();  
    pad->cd();  
     for(Int_t k=0; k <(kk); k++){  
       stringa.str("");  
       stringa << "old data from "<< salto[kk]<< " to "<<salto1[kk+1];  
       min2=min2-0.1;  
       lat->DrawLatex(0.08, min2,stringa.str().c_str());  
       //cout<<salto[kk]<<";\n";  
     }  
     }*/  
490    
491     finalCanv->cd(5);     finalCanv->cd(5);
492     TPaveText *pt4=0;     TPaveText *pt4=0;
493     pt4 = new TPaveText (.38 ,.91,.92,1);     pt4 = new TPaveText (.70,.91,.90,1);
494     pt4->AddText(oss3.str().c_str());     pt4->AddText(oss3.str().c_str());
495     pt4->SetTextColor(kBlue);     pt4->SetTextColor(kBlue);
496     pt4->SetFillColor(10);     pt4->SetFillColor(10);
497     pt4->SetBorderSize(0);     pt4->SetBorderSize(0);
498     pt4->Draw();     pt4->Draw();
499    
500     // finalCanv->Draw();     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();     finalCanv->Update();
544    
545     oss1.str("");     oss1.str("");
# Line 507  void HeaderScan(TString base, TString ou Line 548  void HeaderScan(TString base, TString ou
548    
549     finalCanv->SaveAs(oss1.str().c_str());     finalCanv->SaveAs(oss1.str().c_str());
550    
     
   
551     file->Close();     file->Close();
552        
553  }  }
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.23