/[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.9 by pam-rm2, Fri Sep 8 08:26:47 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 106  void HeaderScan(TString base, TString ou Line 113  void HeaderScan(TString base, TString ou
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();    nevents = PhysicsTr->GetEntries();
121    const Int_t size = nevents;    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();
133    
# Line 177  void HeaderScan(TString base, TString ou Line 184  void HeaderScan(TString base, TString ou
184    //------------------- fill vectors and histogram -----------------------------//    //------------------- fill vectors and histogram -----------------------------//
185    Double_t *PscuCounter = new Double_t[size];    Double_t *PscuCounter = new Double_t[size];
186    Double_t *PacketLenght = new Double_t[size];    Double_t *PacketLenght = new Double_t[size];
187    Double_t *OBTime = new Double_t[size];    ULong_t utile=0;
188      double *OBTime = new Double_t[size];
189    Double_t *Eventsperminute= new Double_t[nbin];    Double_t *Eventsperminute= new Double_t[nbin];
190    Double_t *Minute= new Double_t[nbin];    Double_t *Minute= new Double_t[nbin];
191    Int_t max=0;    Double_t max=0;
192    for (Int_t k = 0; k < nevents; k++){    for (Int_t k = 0; k < nevents; k++){
193      headBr->GetEntry(k);      headBr->GetEntry(k);
194      ph = eh->GetPscuHeader();        ph = eh->GetPscuHeader();
195      h1->Fill(ph->GetOrbitalTime());      utile= ph->GetOrbitalTime();
196        h1->Fill(utile);
197      PscuCounter[k]= ph->GetCounter();      PscuCounter[k]= ph->GetCounter();
198      PacketLenght[k]=ph->GetPacketLenght();      PacketLenght[k]=ph->GetPacketLenght();
199      OBTime[k]=ph->GetOrbitalTime();      OBTime[k]=utile;
200    }      }  
201    int l=0;    int l=0;
202    while(l<nbin){    while(l<nbin){
203      Eventsperminute[l]=h1->GetBinContent(l);      Eventsperminute[l]=h1->GetBinContent(l);
204      Minute[l]=firstime+l*60000;      Minute[l]=firstime+l*60000;
205      if(h1->GetBinContent(l) >= max)max =(Int_t)h1->GetBinContent(l);      if(h1->GetBinContent(l) >= max)max =h1->GetBinContent(l);
206      l++;      l++;
207    }    }
208    
# Line 203  void HeaderScan(TString base, TString ou Line 212  void HeaderScan(TString base, TString ou
212    TMultiGraph *packeCounter = new TMultiGraph();    TMultiGraph *packeCounter = new TMultiGraph();
213    
214    oss1.str("");    oss1.str("");
215    oss1 << filename.Data() <<": Physics Packet per minute.  Start time = " << obmin << ", End time = "<< obmax <<" ms";    oss1 << "Physics Packet per minute.  Start time = " << obmin << ", End time = "<< obmax <<" ms";
216    TGraph *rate1= new TGraph(nbin, (const Double_t*)Minute, (const Double_t*)Eventsperminute);    TGraph *rate1= new TGraph(nbin, (const Double_t*)Minute, (const Double_t*)Eventsperminute);
217    rate1->SetMarkerColor(kBlack);    rate1->SetMarkerColor(kBlack);
218    rate1->SetMarkerSize(.1);    rate1->SetMarkerSize(.1);
219    rate1->SetMarkerStyle(21);    rate1->SetMarkerStyle(21);
220    rate->Add(rate1);    rate->Add(rate1);
221    
222    TGraph *packetLength1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PacketLenght);    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    packetLength1->SetMarkerColor(2);    packetLength1->SetMarkerColor(2);
226    packetLength1->SetMarkerSize(.5);    packetLength1->SetMarkerSize(.3);
227    packetLength1->SetMarkerStyle(21);    packetLength1->SetMarkerStyle(21);
228    packetLength->Add(packetLength1);    packetLength->Add(packetLength1);
229    
230    TGraph *packeCounter1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PscuCounter);    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    packeCounter1->SetMarkerColor(4);    packeCounter1->SetMarkerColor(4);
234    packeCounter1->SetMarkerSize(.2);    packeCounter1->SetMarkerSize(.2);
235    packeCounter1->SetMarkerStyle(21);    packeCounter1->SetMarkerStyle(21);
# Line 228  void HeaderScan(TString base, TString ou Line 237  void HeaderScan(TString base, TString ou
237    
238    //------------ Create and Draw Canvas ---------------------//    //------------ 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);
# Line 241  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,.45,.90,.75);  
   pad->SetFillColor(10);  
253    
254    TLine li;    TLine li;
255    li.SetLineStyle(4);    li.SetLineStyle(4);
# Line 267  void HeaderScan(TString base, TString ou Line 274  void HeaderScan(TString base, TString ou
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,max);        li.DrawLine(obt,0,obt,max);
# Line 278  void HeaderScan(TString base, TString ou Line 286  void HeaderScan(TString base, TString ou
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 364  void HeaderScan(TString base, TString ou Line 379  void HeaderScan(TString base, TString ou
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    Double_t salto;    ULong_t salto;
388    for (Int_t m=0; m<numkey; m++){    for (Int_t m=0; m<numkey; m++){
389      key = list->At(m);      key = list->At(m);
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];
# Line 390  void HeaderScan(TString base, TString ou Line 407  void HeaderScan(TString base, TString ou
407          ph = eh->GetPscuHeader();            ph = eh->GetPscuHeader();  
408          PscuCounter1[n]= ph->GetCounter();          PscuCounter1[n]= ph->GetCounter();
409          FileOffset1[n]=ph->GetFileOffset();          FileOffset1[n]=ph->GetFileOffset();
410          tempo1[n]=ph->GetOrbitalTime();          utile=ph->GetOrbitalTime();
411            tempo1[n]=utile;
412          if((m==0) && (n==0)){          if((m==0) && (n==0)){
413            primotempo=ph->GetOrbitalTime();            primotempo=ph->GetOrbitalTime();
414            salto=ph->GetOrbitalTime();            salto=ph->GetOrbitalTime();
# Line 406  void HeaderScan(TString base, TString ou Line 424  void HeaderScan(TString base, TString ou
424          }          }
425          n++;          n++;
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 431  void HeaderScan(TString base, TString ou Line 450  void HeaderScan(TString base, TString ou
450     all1->cd();     all1->cd();
451    
452     oss1.str("");     oss1.str("");
453     oss1 << filename.Data() <<": PscuCounter vs FileOffset.";     oss1 <<"PscuCounter vs FileOffset.";
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();
# Line 456  void HeaderScan(TString base, TString ou Line 475  void HeaderScan(TString base, TString ou
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 at "<<primotempo <<" ms, last packet at "<<ultimotempo<<" ms.";     oss3 << "OBT vs FileOffset";
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();
# Line 469  void HeaderScan(TString base, TString ou Line 487  void HeaderScan(TString base, TString ou
487     mg2->GetYaxis()->SetTitleSize(0.06);     mg2->GetYaxis()->SetTitleSize(0.06);
488     mg2->GetYaxis()->SetLabelSize(0.05);     mg2->GetYaxis()->SetLabelSize(0.05);
489     mg2->GetYaxis()->SetTitleOffset(0.6);     mg2->GetYaxis()->SetTitleOffset(0.6);
    pad->Draw();  
    pad->cd();  
    stringa.str("");  
    stringa << "jump at: "<<salto<<" ms \n";  
    if(salto != primotempo) lat->DrawLatex(0.08, 0.8,stringa.str().c_str());  
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->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("");
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.23