/[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.2 by pam-rm2, Sat Jun 17 11:02:09 2006 UTC revision 1.13 by pam-rm2, Thu Mar 8 09:53:44 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 49  Line 49 
49  * Date 16 June 2006 - Malvezzi  * Date 16 June 2006 - Malvezzi
50  *  *
51  * Description of changes:  * Description of changes:
52  *                         Implementation of case: numebr of events are <= 0.  *                         Implementation of the case: numebr of events <= 0.
53  *                         Remove of the graph "grPcktId1"; see PacketScan for the same information.  *                         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 all graphs have been converted in 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.
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    
# Line 64  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 76  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){  
 {  
98    
99    //------------------- Variables initilization -------------------------//    //------------------- Variables initilization -------------------------//
100      Long64_t nevents=0, runnevents;
101    TList *list = new TList;    ULong_t lastime, firstime, primotempo, ultimotempo, primoffset=500000000, ultimoffset;
   Int_t numkey;  
   TObject *key = new TObject;  
   const char *name;  
   //TTree* tr = new TTree;  
   //TTree* tr1 = new TTree;  
   Long64_t nevents=0; // ev=0, events=0;  
   ULong_t lastime, firstime;  
102    double obmin=0.;    double obmin=0.;
103    double obmax=0.;    double obmax=0.;
104    double obt;    stringstream oss, oss1, oss2, oss3, noentries, stringa;
105    ULong_t maxoffset, minoffset;  
106    Float_t  pcktLenght =0., pcktCounter=0., offset=0., allCounter=0.;    //-------------- Load root file, tree and branch  -------------------//
   stringstream oss, oss1, oss2, oss3, noentries;  
   //------- load root file --------------  
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 = ".";  
   list = file->GetListOfKeys(); //get list of trees in the file  
   numkey = file->GetNkeys();   //get number of trees in the file  
     
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 = PhysicsTr->GetEntries();
121      const Int_t size = nevents;
122    
123    nevents = headBr->GetEntries();    TTree *RunHeadTr = (TTree*)file->Get("RunHeader");  ///run header tree
124    const Int_t sizetot = nevents;    pamela::EventHeader *eH= new pamela::EventHeader;
125      pamela::RunHeaderEvent *reh=new pamela::RunHeaderEvent;
126      
127      RunHeadTr->SetBranchAddress("Header",&eH);
128      RunHeadTr->SetBranchAddress("RunHeader",&reh);
129      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 126  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 141  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 ----------------------------//
163       headBr->GetEntry(0);
164    ///-------------- to know the max and min File Offset ----------------------------//     ph = eh->GetPscuHeader();
165    /*for (Int_t i=0; i<numkey; i++){     firstime = ph->GetOrbitalTime();
166      key = list->At(i);     headBr->GetEntry(nevents);
167      name=(key->GetName());     ph = eh->GetPscuHeader();
168      tr = (TTree*)file->Get(name);     lastime = ph->GetOrbitalTime();
169      if (tr->IsZombie()) continue;     int i =0;
170      events = tr->GetEntries();     while(i < nevents){
171      tr->SetBranchAddress("Header", &eh);         headBr->GetEntry(i);
172      for (Int_t j = 0; j < events; j++){       ph = eh->GetPscuHeader();
173        tr->GetEntry(j);       if((ph->GetOrbitalTime()) <= firstime) firstime=ph->GetOrbitalTime();
174        ph = eh->GetPscuHeader();       if((ph->GetOrbitalTime()) >= lastime) lastime=ph->GetOrbitalTime();
175        ph->GetFileOffset();       i++;
176        if(ph->GetFileOffset() <= minoffset) minoffset=ph->GetFileOffset();       }
177        if(ph->GetFileOffset() >= maxoffset) maxoffset=ph->GetFileOffset();    
     }  
     }*/  
   //*************************** Histograms ************************************************************//  
178    //------------------------ First histogram -----------------------------------//    //------------------------ First histogram -----------------------------------//
   headBr->GetEntry(0);  
   ph = eh->GetPscuHeader();  
   firstime = ph->GetOrbitalTime();  
   headBr->GetEntry(nevents-1);  
   ph = eh->GetPscuHeader();  
   lastime = ph->GetOrbitalTime();  
179    obmin=firstime;    obmin=firstime;
180    obmax=lastime;      obmax=lastime;
   oss1.str("");  
   oss1 << filename.Data() << ": Physics Packet per minute;" <<" start @ " << firstime << ", end @ "<< lastime <<"ms";  
181    Int_t  nbin   = (lastime-firstime)/60000;    Int_t  nbin   = (lastime-firstime)/60000;
182    TH1F *h1 = new TH1F ("histo1", oss1.str().c_str(), nbin, obmin, obmax);    TH1F *h1 = new TH1F ("histo1", "" , nbin, obmin, obmax);
183    
184    //------------------------ Second histogram -----------------------------------//    //------------------- 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      Int_t max=0;
192      for (Int_t k = 0; k < nevents; k++){
193        headBr->GetEntry(k);
194        ph = eh->GetPscuHeader();  
195        h1->Fill(ph->GetOrbitalTime());
196        PscuCounter[k]= ph->GetCounter();
197        PacketLenght[k]=ph->GetPacketLenght();
198        utile= ph->GetOrbitalTime();
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 =(Int_t)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 nint = (lastime-firstime);    packetLength1->SetMarkerColor(2);
226    TH1F *PcktLenght = new TH1F ("histo2", oss2.str().c_str(), nint, obmin, obmax);    packetLength1->SetMarkerSize(.3);
227      packetLength1->SetMarkerStyle(21);
228      packetLength->Add(packetLength1);
229    
230   //------------------------ Third histogram -----------------------------------//    TGraph *packeCounter1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PscuCounter);
231    oss3.str("");    oss3.str("");
232    oss3 << filename.Data() << ": Physics Counter";    oss3 <<"Physics Counter vs. OBT";
233    TH1F *PcktCounter = new TH1F ("histo3", oss3.str().c_str(), nint, obmin, obmax);    packeCounter1->SetMarkerColor(4);
234      packeCounter1->SetMarkerSize(.2);
235   //------------------------ Fourth histogram -----------------------------------//    packeCounter1->SetMarkerStyle(21);
236    /*oss3.str("");    packeCounter->Add(packeCounter1);
237    oss3 << filename.Data() << ": Packet Counter";  
238    Int_t nintoffset = (maxoffset-minoffset)/10;    //------------ Create and Draw Canvas ---------------------//
239    TH1F *FileOffset = new TH1F ("histo4", oss3.str().c_str(), nintoffset, minoffset, maxoffset);    TCanvas *finalCanv = new TCanvas("Header", base, 1200, 1600);
240    */    finalCanv->Divide(1,6);
241    //**************************************************************************************************//    finalCanv->SetFillColor(10);
242    //------- fill histograms ---------//  
243      TPad *all2= new TPad ("","", 0, 0, 1, 1);
244    for (Int_t i = 0; i < nevents; i++){    all2->SetFillColor(10);
245      headBr->GetEntry(i);    TPad *all3= new TPad ("","", 0, 0, 1, 1);
246      ph = eh->GetPscuHeader();    all3->SetFillColor(10);
247      obt = ph->GetOrbitalTime();      TPad *all4= new TPad ("","", 0, 0, 1, 1);
248      pcktLenght=ph->GetPacketLenght();    all4->SetFillColor(10);
249      pcktCounter=ph->GetCounter();    TPad *all= new TPad ("","", 0, 0, 1, 1);
250      h1->Fill(obt);      all->SetFillColor(10);
251      PcktLenght->Fill(obt,pcktLenght);      TPad *all1= new TPad ("","", 0, 0, 1, 1);
252      PcktCounter->Fill(obt,pcktCounter);    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        Int_t obt = ph->GetOrbitalTime();
278        if (ws==1){
279          li.SetLineColor(3);
280          li.DrawLine(obt,0,obt,max);
281        }
282        else if (ws==2){
283          li.SetLineColor(4);
284          li.DrawLine(obt,0,obt,max);
285        }  
286    }    }
287      
288      RunHeadTr->GetEntry(0);
289      ph = eH->GetPscuHeader();
290      ULong_t TimeSync = reh->LAST_TIME_SYNC_INFO;
291      ULong_t ObtSync = reh->OBT_TIME_SYNC;
292      //cout<<"TimeSync "<<reh->LAST_TIME_SYNC_INFO<<"\n";
293      //cout<<"ObtSync "<<reh->OBT_TIME_SYNC<<"\n";
294    
295      finalCanv->cd(1);
296      stringstream ws1, ws2;
297      ws1.str("");
298      ws2.str("");
299      ws1<<"ACQ_SETTING_MODE = 1";
300      ws2<<"ACQ_SETTING_MODE = 2";
301      TPaveText *pt=0;
302      pt = new TPaveText (.60,.92,.76,.98);
303      pt->AddText(ws1.str().c_str());
304      pt->SetTextColor(3);
305      pt->SetFillColor(10);
306      pt->SetBorderSize(0);
307      pt->Draw();
308      TPaveText *pt1=0;
309      pt1 = new TPaveText (.76,.92,.92,.98);
310      pt1->AddText(ws2.str().c_str());
311      pt1->SetTextColor(4);
312      pt1->SetFillColor(10);
313      pt1->SetBorderSize(0);
314      pt1->Draw();
315      pt1 = new TPaveText (.05,.91,.6,1);
316      pt1->AddText(oss1.str().c_str());
317      pt1->SetTextColor(1);
318      pt1->SetFillColor(10);
319      pt1->SetBorderSize(0);
320      pt1->Draw();
321    
322      //----------- Second PAD -------------------------------//
323      finalCanv->cd(2);
324      all3->Draw();
325      all3->cd();
326      packetLength->Draw("AP");
327      packetLength->GetXaxis()->SetTitle("OBT    (ms)");
328      packetLength->GetXaxis()->SetTitleSize(0.05);
329      packetLength->GetXaxis()->CenterTitle();
330      packetLength->GetXaxis()->SetLabelSize(0.05);
331      packetLength->GetYaxis()->SetTitle("Lenght  (byte)");
332      packetLength->GetYaxis()->CenterTitle();
333      packetLength->GetYaxis()->SetLabelSize(0.05);
334      packetLength->GetYaxis()->SetTitleSize(0.06);
335      packetLength->GetYaxis()->SetTitleOffset(0.6);
336    
337    
338      finalCanv->cd(2);
339      pt = new TPaveText (.6,.91,.90,1);
340      pt->AddText(oss2.str().c_str());
341      pt->SetTextColor(2);
342      pt->SetFillColor(10);
343      pt->SetBorderSize(0);
344      pt->Draw();
345    
346      //----------- Third PAD -------------------------------//
347      finalCanv->cd(3);
348      all4->Draw();
349      all4->cd();
350      packeCounter->Draw("AP");
351      packeCounter->GetXaxis()->SetTitle("OBT    (ms)");
352      packeCounter->GetXaxis()->SetTitleSize(0.05);
353      packeCounter->GetXaxis()->CenterTitle();
354      packeCounter->GetXaxis()->SetLabelSize(0.05);
355      packeCounter->GetYaxis()->SetTitle("Counter");
356      packeCounter->GetYaxis()->SetTitleSize(0.05);
357      packeCounter->GetYaxis()->CenterTitle();
358      packeCounter->GetYaxis()->SetLabelSize(0.05);
359      packeCounter->GetYaxis()->SetTitleSize(0.06);
360      packeCounter->GetYaxis()->SetTitleOffset(0.6);
361    
362    
363      finalCanv->cd(3);
364      TPaveText *pt2=0;
365      pt2 = new TPaveText (.6,.91,.90,1);
366      pt2->AddText(oss3.str().c_str());
367      pt2->SetTextColor(4);
368      pt2->SetFillColor(10);
369      pt2->SetBorderSize(0);
370      pt2->Draw();  
371    
372      /**********************************************************************************************/
373    
374    /*for (Int_t i=0; i<numkey; i++){    TMultiGraph *mg1 = new TMultiGraph();
375      key = list->At(i);    TMultiGraph *mg2 = new TMultiGraph();
376      //---------------------- fill vectors and histogram --------------------------------------------------//
377      TList *list = new TList;
378      Int_t numkey;
379      TObject *key = new TObject;
380      const char *name;
381      char *SoftInfo="SoftInfo";
382      TTree* tr = new TTree;
383      Long64_t nevntskey=0;
384      list = file->GetListOfKeys();
385      numkey = file->GetNkeys();
386      ULong_t salto;
387      for (Int_t m=0; m<numkey; m++){
388        key = list->At(m);
389      name=(char *)(key->GetName());      name=(char *)(key->GetName());
390      tr1 = (TTree*)file->Get(name);      if(strcmp(name,SoftInfo)==0)continue;
391      if (tr1->IsZombie()) continue;      tr = (TTree*)file->Get(name);
392      tr1->SetBranchAddress("Header", &eh);      if (tr->IsZombie()) continue;
     ev = tr1->GetEntries();  
     for (Int_t j = 0; j < ev; j++){  
       tr1->GetEntry(j);    
       ph = eh->GetPscuHeader();  
       offset=ph->GetFileOffset();  
       allCounter= ph->GetCounter();  
       FileOffset->Fill(offset,allCounter);  
     }  
     }*/  
   
   //****************************** Canvases *******************************//  
   //TGaxis::SetMaxDigits(4);  
   //------------------- First Canvas --------------------------------//  
   TCanvas *finalCanv1 = new TCanvas("Header_1", base, 1280, 1024);  
   finalCanv1->Divide(1,3);  
   finalCanv1->SetFillColor(10);  
   
   finalCanv1->cd(1);  
   h1->SetStats(kFALSE);  
   h1->GetXaxis()->SetTitle("OBT   (ms)");  
   h1->GetXaxis()->CenterTitle();  
   h1->GetXaxis()->SetLabelSize(0.04);  
   h1->GetYaxis()->SetTitle("Number of events ");  
   h1->GetYaxis()->CenterTitle();  
   h1->GetYaxis()->SetTitleSize(0.06);  
   h1->GetYaxis()->SetTitleOffset(0.8);  
   h1->Draw();  
     
   finalCanv1->cd(2);  
   PcktLenght->SetStats(kFALSE);  
   PcktLenght->GetXaxis()->SetTitle("OBT    (ms)");  
   PcktLenght->GetXaxis()->CenterTitle();  
   PcktLenght->GetXaxis()->SetLabelSize(0.04);  
   PcktLenght->GetYaxis()->SetTitle("Lenght  (byte)");  
   PcktLenght->GetYaxis()->CenterTitle();  
   PcktLenght->GetYaxis()->SetLabelSize(0.04);  
   PcktLenght->GetYaxis()->SetTitleSize(0.06);  
   PcktLenght->GetYaxis()->SetTitleOffset(0.8);  
   PcktLenght->SetMarkerColor(2);  
   PcktLenght->SetMarkerSize(.5);  
   PcktLenght->SetMarkerStyle(21);  
   PcktLenght->Draw("9p");  
   
   finalCanv1->cd(3);  
   PcktCounter->SetStats(kFALSE);  
   PcktCounter->GetXaxis()->SetTitle("OBT    (ms)");  
   PcktCounter->GetXaxis()->SetTitleSize(0.05);  
   PcktCounter->GetXaxis()->CenterTitle();  
   PcktCounter->GetXaxis()->SetLabelSize(0.04);  
   PcktCounter->GetYaxis()->SetTitle("Counter");  
   PcktCounter->GetYaxis()->SetTitleSize(0.05);  
   PcktCounter->GetYaxis()->CenterTitle();  
   PcktCounter->GetYaxis()->SetLabelSize(0.04);  
   PcktCounter->GetYaxis()->SetTitleSize(0.06);  
   PcktCounter->GetYaxis()->SetTitleOffset(0.8);  
   PcktCounter->SetMarkerColor(4);  
   PcktCounter->SetMarkerSize(.5);  
   PcktCounter->SetMarkerStyle(21);  
   PcktCounter->Draw("9p");  
   //--------- Second Canvas -----------------------------//  
   /*TCanvas *finalCanv2 = new TCanvas("Header_2", base, 1280, 1024);  
   finalCanv2->Divide(1,2);  
   finalCanv2->SetFillColor(10);  
   
   finalCanv2->cd(1);  
   PcktCounter->SetStats(kFALSE);  
   PcktCounter->GetXaxis()->SetTitle("OBT    (ms)");  
   PcktCounter->GetXaxis()->SetTitleSize(0.05);  
   PcktCounter->GetXaxis()->CenterTitle();  
   PcktCounter->GetXaxis()->SetLabelSize(0.04);  
   PcktCounter->GetYaxis()->SetTitle("Counter");  
   PcktCounter->GetYaxis()->SetTitleSize(0.05);  
   PcktCounter->GetYaxis()->CenterTitle();  
   PcktCounter->GetYaxis()->SetLabelSize(0.04);  
   PcktCounter->SetMarkerColor(4);  
   PcktCounter->SetMarkerSize(.5);  
   PcktCounter->SetMarkerStyle(21);  
   PcktCounter->Draw("9p");  
     
   finalCanv2->cd(2);  
   FileOffset->SetStats(kFALSE);  
   FileOffset->GetXaxis()->SetTitle("File Offset");  
   FileOffset->GetXaxis()->CenterTitle();  
   FileOffset->GetXaxis()->SetTitleSize(0.05);  
   FileOffset->GetXaxis()->SetLabelSize(0.04);  
   FileOffset->GetYaxis()->SetTitle("Counter");  
   FileOffset->GetYaxis()->CenterTitle();  
   FileOffset->GetYaxis()->SetTitleSize(0.05);  
   FileOffset->GetYaxis()->SetLabelSize(0.04);  
   FileOffset->SetMarkerColor(3);  
   FileOffset->SetMarkerSize(.5);  
   FileOffset->SetMarkerStyle(21);  
   FileOffset->Draw("9p");*/  
     
   //-------to save---------------------------///  
   oss.str("");  
   oss1.str("");  
     
   oss << outDir.Data() << filename.Data();  
   oss1 << outDir.Data() << filename.Data();  
393    
394    oss  << "_HeaderScan_1." << format.Data();      tr->SetBranchAddress("Header", &eh);
395    oss1 << "_HeaderScan_2." << format.Data();      TBranch *Br = tr->GetBranch("Header");
396          nevntskey = tr->GetEntries();
397    finalCanv1->SaveAs(oss.str().c_str());      if(nevntskey !=0){
398    //finalCanv2->SaveAs(oss1.str().c_str());        Int_t size1=nevntskey;
399            Double_t *PscuCounter1 = new Double_t[size1];
400    file->Close();        Double_t *FileOffset1 = new Double_t[size1];
401            Double_t *tempo1 = new Double_t[size1];
402  }  
403          int n=0;
404          while(n<nevntskey){
405            Br->GetEntry(n);
406            ph = eh->GetPscuHeader();  
407            PscuCounter1[n]= ph->GetCounter();
408            FileOffset1[n]=ph->GetFileOffset();
409            utile=ph->GetOrbitalTime();
410            tempo1[n]=utile;
411            if((m==0) && (n==0)){
412              primotempo=ph->GetOrbitalTime();
413              salto=ph->GetOrbitalTime();
414            }
415            if(salto > ph->GetOrbitalTime())salto=ph->GetOrbitalTime();
416            if(ph->GetFileOffset()<= primoffset){
417              primoffset=ph->GetFileOffset();
418              primotempo=ph->GetOrbitalTime();
419            }
420            if(ph->GetFileOffset()>=ultimoffset){
421              ultimotempo=ph->GetOrbitalTime();
422              ultimoffset=ph->GetFileOffset();
423            }
424            n++;
425          }
426    
427          TGraph *graph3= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)PscuCounter1);
428          graph3->SetMarkerColor(3);
429          graph3->SetMarkerSize(.2);
430          graph3->SetMarkerStyle(21);
431          mg1->Add(graph3);
432          
433          TGraph *graph4= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)tempo1);
434          graph4->SetMarkerColor(kBlue);
435          graph4->SetMarkerSize(.2);
436          graph4->SetMarkerStyle(21);
437          mg2->Add(graph4);
438        }
439      }
440    
441      TLatex *lat = new TLatex();
442      lat->SetTextAlign(12);
443      lat->SetTextSize(0.15);
444      lat->SetTextColor(kBlue);
445    
446     //------------ Fourth PAD  ---------------------//
447       finalCanv->cd(4);
448       all1->Draw();
449       all1->cd();
450    
451       oss1.str("");
452       oss1 <<"PscuCounter vs FileOffset.";
453       mg1->Draw("AP");
454       mg1->GetXaxis()->SetTitle("File Offset");
455       mg1->GetXaxis()->CenterTitle();
456       mg1->GetXaxis()->SetTitleOffset(0.8);
457       mg1->GetXaxis()->SetTitleSize(0.05);
458       mg1->GetXaxis()->SetLabelSize(0.05);
459       mg1->GetYaxis()->SetTitle("Counter");
460       mg1->GetYaxis()->CenterTitle();
461       mg1->GetYaxis()->SetTitleSize(0.06);
462       mg1->GetYaxis()->SetLabelSize(0.06);
463       mg1->GetYaxis()->SetTitleOffset(0.6);
464       finalCanv->cd(4);
465       TPaveText *pt3=0;
466       pt3 = new TPaveText (.60,.91,.90,1);
467       pt3->AddText(oss1.str().c_str());
468       pt3->SetTextColor(3);
469       pt3->SetFillColor(10);
470       pt3->SetBorderSize(0);
471       pt3->Draw();
472    
473     //------------ Fifth PAD  ---------------------//
474       finalCanv->cd(5);
475       all->Draw();
476       all->cd();
477       oss3.str("");
478       oss3 << "OBT vs FileOffset";
479       mg2->Draw("AP");
480       mg2->GetXaxis()->SetTitle("File Offset");
481       mg2->GetXaxis()->CenterTitle();
482       mg2->GetXaxis()->SetTitleSize(0.05);
483       mg2->GetXaxis()->SetLabelSize(0.05);
484       mg2->GetYaxis()->SetTitle("OBT");
485       mg2->GetYaxis()->CenterTitle();
486       mg2->GetYaxis()->SetTitleSize(0.06);
487       mg2->GetYaxis()->SetLabelSize(0.05);
488       mg2->GetYaxis()->SetTitleOffset(0.6);
489    
490       finalCanv->cd(5);
491       TPaveText *pt4=0;
492       pt4 = new TPaveText (.70,.91,.90,1);
493       pt4->AddText(oss3.str().c_str());
494       pt4->SetTextColor(kBlue);
495       pt4->SetFillColor(10);
496       pt4->SetBorderSize(0);
497       pt4->Draw();
498    
499       finalCanv->cd(6);
500       ULong_t primotempoABS=TimeSync+((primotempo/1000)-ObtSync);
501       ULong_t obmaxABS=TimeSync+((lastime/1000)-ObtSync);
502       ULong_t saltoABS=TimeSync+((salto/1000)-ObtSync);
503       ULong_t ultimotempoABS=TimeSync+((ultimotempo/1000)-ObtSync);
504    
505       TPaveText *pt5=0;
506       pt5 = new TPaveText (0,0,1,1);
507       stringa.str("");
508       stringa << "  Filename:   "<<filename.Data()<<"\n";
509       TText *t1=pt5->AddText(0.25,0.95,stringa.str().c_str());
510       t1->SetTextSize(0.1);
511       stringa.str("");
512       stringa << "  OBT (ms)            ABS TIME (s)";
513       TText *t2=pt5->AddText(0.32,0.75,stringa.str().c_str());
514       t2->SetTextSize(0.07);
515       stringa.str("");
516       stringa << "New data start at:   "<<primotempo<<"         "<<primotempoABS;
517       TText *t3=pt5->AddText(0.25,0.60,stringa.str().c_str());
518       t3->SetTextSize(0.08);
519       stringa.str("");
520       stringa << "New data end at:     "<<lastime<<"         "<<obmaxABS;
521       TText *t4=pt5->AddText(0.25,0.50,stringa.str().c_str());
522       t4->SetTextSize(0.08);
523       if(primotempo!=salto || lastime!=ultimotempo){
524         stringa.str("");
525         stringa << "Old data start at:   "<<salto<<"        "<<saltoABS;
526         TText *t5=pt5->AddText(0.65,0.60,stringa.str().c_str());
527         t5->SetTextSize(0.08);
528         stringa.str("");
529         stringa << "Old data end at:     "<<ultimotempo<<"        "<<ultimotempoABS;
530         TText *t6=pt5->AddText(0.65,0.50,stringa.str().c_str());
531         t6->SetTextSize(0.08);
532         stringa.str("");
533         stringa << "  OBT (ms)          ABS TIME (s)";
534         TText *t2=pt5->AddText(0.72,0.75,stringa.str().c_str());
535         t2->SetTextSize(0.07);
536       }
537       pt5->SetTextColor(kBlack);
538       pt5->SetFillColor(10);
539       pt5->SetBorderSize(0);
540       pt5->Draw();
541    
542       finalCanv->Update();
543    
544       oss1.str("");
545       oss1 << outDir.Data() << filename.Data();
546       oss1 << "_HeaderScan"<<"." << format.Data();
547    
548       finalCanv->SaveAs(oss1.str().c_str());
549    
550       file->Close();
551      
552    }
553    
554  int main(int argc, char* argv[]){  int main(int argc, char* argv[]){
555    TString path;    TString path;
# Line 381  int main(int argc, char* argv[]){ Line 594  int main(int argc, char* argv[]){
594    }    }
595    HeaderScan(argv[1], outDir, format);    HeaderScan(argv[1], outDir, format);
596  }  }
597    
Legend:
Removed from v.1.1.1.2  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.23