quicklook/QLflightTmtc_Header/HeaderScan.cpp

/**
* Header_Scan
* Author  Nagni
* version 1.0
*
* Version 1.1 - 28 December 2004
* If outList does not exist the function exit to prompt
* If nevents = 0 the function exit to promt
*
* Version 1.2 - 3 January 2005
* Two canvases are created to see the graphs better
*
* Version 1.25 - 13 January 2005
* Changed Int_t to Float because of variable range size
* (UInt_t has been excluded beacuse of uncompatibility with TGraph)
*
* version 1.3  - 22 February 2005 - Nagni
* For compatibility with batch mode excution:
* 1) Added "include <iostream>" and "using namespace std" 
* 2) Removed gROOT->Reset()
*
* Version 1.4
* Date 08 March 2005
* Added "format" parameter for saving the produced image in various formats
* (for a complete list of types refer to TPad::SaveAs method)
*
* Version 1.5
* Date 09 February 2006 - Marcelli
* Update to work with new Yoda output
*
* Version 1.6
* Date 27 February 2006 - Marcelli
* For compilation: 
*     Added function "int main(int argc, char* argv[])"
*
*
* Description: This script creates two canvases with five pads. The first pad shows packetID variable (for all packets) vs. OBT.
*              The second pad shows the number of physic packets vs. OBT. The third pad shows the lenght of Physic packets (byte) vs. OBT.
*              The fourth pad shows the packetcounter of physic packets vs. OBT. The fifth pad shows PacketCounter vs. File Offset.
*
* Parameters:
*              TSTring base          - the path to the root directory for the specific Pamela unpack session
*                                      There is no default value, without this input the program will not run
*              TString outDir        - the path where to save the output image (Default = ./) 
*              TString format        - the format which will be used fo rsave the produced images (Default = "jpg") 
*
*
* Version 1.7
* Date 16 June 2006 - Malvezzi
*
* Description of changes:
*                         Implementation of the case: numebr of events <= 0.
*                         Remove graph "grPcktId1"; see PacketScan for the same information.
*             Fixed bugs: for a large namber of events is not possible to have vectors, so I have subsituted graphs with histograms
*                         or divided the graphs in more than one canvas.
*
* Version 1.8
* Date 8 August 2006 - Malvezzi
*
* Description: changed the scale in the second and third graph of the first canvas; added a pad of text in the second canvas 
*
*/


#include <fstream>
#include <math.h>
#include "TLatex.h"
#include "TF1.h"
#include "TPaveText.h"
#include "TMultiGraph.h"
#include <sstream>
#include <iostream>
#include <stdio.h>
#include <string.h>
#include "TString.h"
#include "TStyle.h"
#include "TFile.h"
#include "TList.h"
#include "TTree.h"
#include "TObjString.h"
#include "TCanvas.h"
#include "TGraph.h"
#include "TH1F.h"
#include "TGaxis.h"
#include "EventHeader.h"
#include "PscuHeader.h"
#include "RunHeaderEvent.h"
#include "TPaveText.h"

using namespace std;

void HeaderScan(TString base, TString outDir, TString format){  

  //------------------- Variables initilization -------------------------//
  Long64_t nevents=0, runnevents; 
  ULong_t lastime, firstime, primotempo, ultimotempo, primoffset=500000000, ultimoffset;
  double obmin=0.;
  double obmax=0.; 
  stringstream oss, oss1, oss2, oss3, noentries, stringa;

  //-------------- Load root file, tree and branch  -------------------//
  TFile *file = new TFile(base.Data());
  if (!file){
    printf("No such file in the directory has been found");
    return;
  } 
 
  TTree *PhysicsTr = (TTree*)file->Get("Physics");
  TBranch *headBr = PhysicsTr->GetBranch("Header");

  pamela::EventHeader  *eh   = new pamela::EventHeader;
  pamela::PscuHeader   *ph   = new pamela::PscuHeader;

  PhysicsTr->SetBranchAddress("Header", &eh);
  nevents = PhysicsTr->GetEntries();
  const Int_t size = nevents;

  TTree *RunHeadTr = (TTree*)file->Get("RunHeader");  ///run header tree
  pamela::EventHeader *eH= new pamela::EventHeader;
  pamela::RunHeaderEvent *reh=new pamela::RunHeaderEvent;
  
  RunHeadTr->SetBranchAddress("Header",&eH);
  RunHeadTr->SetBranchAddress("RunHeader",&reh);
  runnevents = RunHeadTr->GetEntries();

  TString filename = ((TObjString*)base.Tokenize('/')->Last())->GetString(); 
  filename = ((TObjString*)filename.Tokenize('.')->First())->GetString();

  //----------- If nevents < = 0 ---------------------------------/
  if (nevents<=0) {
    printf("nevents = %i \n", nevents);
    printf(" \n");
    TCanvas *canv = new TCanvas("No entries", "No entries ", 400, 200);
    canv->SetFillColor(10);
    canv->cd();
    TLatex *l = new TLatex();
    l->SetTextAlign(12);
    l->SetTextSize(0.15);
    l->SetTextColor(2);
    noentries.str("");
    noentries<< "HeaderScan_QL:";
    l->DrawLatex(0.05, 0.7, noentries.str().c_str()); 
    noentries.str("");
    noentries<< "No Physics entries for this files";
    l->DrawLatex(0.05, 0.5, noentries.str().c_str()); 
    if (outDir == "./") {
      oss.str("");
      oss << filename.Data() << "_HeaderScan_QL." << format.Data();
    } else {
      oss.str("");
      oss << outDir.Data() << filename.Data() << "_HeaderScan_QL." << format.Data();
    }  
    canv->Update();
    canv->SaveAs(oss.str().c_str()); 
    return;
  }
  //-------------- to know the max and min OBT ----------------------------//
   headBr->GetEntry(0);
   ph = eh->GetPscuHeader();
   firstime = ph->GetOrbitalTime();
   headBr->GetEntry(nevents);
   ph = eh->GetPscuHeader();
   lastime = ph->GetOrbitalTime();
   int i =0;
   while(i < nevents){
     headBr->GetEntry(i);
     ph = eh->GetPscuHeader();
     if((ph->GetOrbitalTime()) <= firstime) firstime=ph->GetOrbitalTime();
     if((ph->GetOrbitalTime()) >= lastime) lastime=ph->GetOrbitalTime();
     i++;
     }
  
  //------------------------ First histogram -----------------------------------//
  obmin=firstime;
  obmax=lastime; 
  Int_t  nbin   = (lastime-firstime)/60000;
  TH1F *h1 = new TH1F ("histo1", "" , nbin, obmin, obmax);

  //------------------- fill vectors and histogram -----------------------------//
  Double_t *PscuCounter = new Double_t[size];
  Double_t *PacketLenght = new Double_t[size];
  Double_t *OBTime = new Double_t[size];
  Double_t *Eventsperminute= new Double_t[nbin];
  Double_t *Minute= new Double_t[nbin];
  Int_t max=0;
  for (Int_t k = 0; k < nevents; k++){
    headBr->GetEntry(k);
    ph = eh->GetPscuHeader();  
    h1->Fill(ph->GetOrbitalTime()); 
    PscuCounter[k]= ph->GetCounter();
    PacketLenght[k]=ph->GetPacketLenght();
    OBTime[k]=ph->GetOrbitalTime();
  }   
  int l=0;
  while(l<nbin){
    Eventsperminute[l]=h1->GetBinContent(l);
    Minute[l]=firstime+l*60000;
    if(h1->GetBinContent(l) >= max)max =(Int_t)h1->GetBinContent(l);
    l++;
  }

  //----------- Graph and MultiGraph -----------------------------------------------//
  TMultiGraph *rate = new TMultiGraph();
  TMultiGraph *packetLength = new TMultiGraph();
  TMultiGraph *packeCounter = new TMultiGraph();

  oss1.str("");
  oss1 << "Physics Packet per minute.  Start time = " << obmin << ", End time = "<< obmax <<" ms";
  TGraph *rate1= new TGraph(nbin, (const Double_t*)Minute, (const Double_t*)Eventsperminute);
  rate1->SetMarkerColor(kBlack);
  rate1->SetMarkerSize(.1);
  rate1->SetMarkerStyle(21);
  rate->Add(rate1);

  TGraph *packetLength1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PacketLenght);
  oss2.str("");
  oss2 <<"Lenght of Physic packets";
  packetLength1->SetMarkerColor(2);
  packetLength1->SetMarkerSize(.3);
  packetLength1->SetMarkerStyle(21);
  packetLength->Add(packetLength1);

  TGraph *packeCounter1= new TGraph(nevents, (const Double_t*)OBTime, (const Double_t*)PscuCounter);
  oss3.str("");
  oss3 <<"Physics Counter vs. OBT";
  packeCounter1->SetMarkerColor(4);
  packeCounter1->SetMarkerSize(.2);
  packeCounter1->SetMarkerStyle(21);
  packeCounter->Add(packeCounter1);

  //------------ Create and Draw Canvas ---------------------//
  TCanvas *finalCanv = new TCanvas("Header", base, 1200, 1600);
  finalCanv->Divide(1,6);
  finalCanv->SetFillColor(10);

  TPad *all2= new TPad ("","", 0, 0, 1, 1);
  all2->SetFillColor(10);
  TPad *all3= new TPad ("","", 0, 0, 1, 1);
  all3->SetFillColor(10);
  TPad *all4= new TPad ("","", 0, 0, 1, 1);
  all4->SetFillColor(10);
  TPad *all= new TPad ("","", 0, 0, 1, 1);
  all->SetFillColor(10);
  TPad *all1= new TPad ("","", 0, 0, 1, 1);
  all1->SetFillColor(10);

  TLine li;
  li.SetLineStyle(4);
  li.SetLineWidth(2);

  //----------- First PAD -------------------------------//
  finalCanv->cd(1);
  all2->Draw();
  all2->cd();
  rate->Draw("ALP");
  rate->GetXaxis()->SetTitle("OBT   (ms)");
  rate->GetXaxis()->SetTitleSize(0.05);
  rate->GetXaxis()->CenterTitle();
  rate->GetXaxis()->SetLabelSize(0.05);
  rate->GetYaxis()->SetTitle("Number of events ");
  rate->GetYaxis()->CenterTitle();
  rate->GetYaxis()->SetLabelSize(0.05);
  rate->GetYaxis()->SetTitleSize(0.06);
  rate->GetYaxis()->SetTitleOffset(0.6);
  for (Int_t l = 0; l < runnevents; l++){
    RunHeadTr->GetEntry(l);
    ph = eH->GetPscuHeader();
    int ws= reh->RM_ACQ_SETTING_MODE;
    int id = ph->GetPacketId1();
    Int_t obt = ph->GetOrbitalTime();
    if (ws==1){
      li.SetLineColor(3);
      li.DrawLine(obt,0,obt,max);
    }
    else if (ws==2){
      li.SetLineColor(4);
      li.DrawLine(obt,0,obt,max);
    }  
  }
  
  RunHeadTr->GetEntry(0);
  ph = eH->GetPscuHeader();
  ULong_t TimeSync = reh->LAST_TIME_SYNC_INFO; 
  ULong_t ObtSync = reh->OBT_TIME_SYNC; 
  //cout<<"TimeSync "<<reh->LAST_TIME_SYNC_INFO<<"\n";
  //cout<<"ObtSync "<<reh->OBT_TIME_SYNC<<"\n";

  finalCanv->cd(1);
  stringstream ws1, ws2;
  ws1.str("");
  ws2.str("");
  ws1<<"ACQ_SETTING_MODE = 1";
  ws2<<"ACQ_SETTING_MODE = 2";
  TPaveText *pt=0;
  pt = new TPaveText (.60,.92,.76,.98);
  pt->AddText(ws1.str().c_str());
  pt->SetTextColor(3);
  pt->SetFillColor(10);
  pt->SetBorderSize(0);
  pt->Draw();
  TPaveText *pt1=0;
  pt1 = new TPaveText (.76,.92,.92,.98);
  pt1->AddText(ws2.str().c_str());
  pt1->SetTextColor(4);
  pt1->SetFillColor(10);
  pt1->SetBorderSize(0);
  pt1->Draw();
  pt1 = new TPaveText (.05,.91,.6,1);
  pt1->AddText(oss1.str().c_str());
  pt1->SetTextColor(1);
  pt1->SetFillColor(10);
  pt1->SetBorderSize(0);
  pt1->Draw();

  //----------- Second PAD -------------------------------//
  finalCanv->cd(2);
  all3->Draw();
  all3->cd();
  packetLength->Draw("AP");
  packetLength->GetXaxis()->SetTitle("OBT    (ms)");
  packetLength->GetXaxis()->SetTitleSize(0.05);
  packetLength->GetXaxis()->CenterTitle();
  packetLength->GetXaxis()->SetLabelSize(0.05);
  packetLength->GetYaxis()->SetTitle("Lenght  (byte)");
  packetLength->GetYaxis()->CenterTitle();
  packetLength->GetYaxis()->SetLabelSize(0.05);
  packetLength->GetYaxis()->SetTitleSize(0.06);
  packetLength->GetYaxis()->SetTitleOffset(0.6);


  finalCanv->cd(2);
  pt = new TPaveText (.6,.91,.90,1);
  pt->AddText(oss2.str().c_str());
  pt->SetTextColor(2);
  pt->SetFillColor(10);
  pt->SetBorderSize(0);
  pt->Draw();

  //----------- Third PAD -------------------------------//
  finalCanv->cd(3);
  all4->Draw();
  all4->cd();
  packeCounter->Draw("AP");
  packeCounter->GetXaxis()->SetTitle("OBT    (ms)");
  packeCounter->GetXaxis()->SetTitleSize(0.05);
  packeCounter->GetXaxis()->CenterTitle();
  packeCounter->GetXaxis()->SetLabelSize(0.05);
  packeCounter->GetYaxis()->SetTitle("Counter");
  packeCounter->GetYaxis()->SetTitleSize(0.05);
  packeCounter->GetYaxis()->CenterTitle();
  packeCounter->GetYaxis()->SetLabelSize(0.05);
  packeCounter->GetYaxis()->SetTitleSize(0.06);
  packeCounter->GetYaxis()->SetTitleOffset(0.6);


  finalCanv->cd(3);
  TPaveText *pt2=0;
  pt2 = new TPaveText (.6,.91,.90,1);
  pt2->AddText(oss3.str().c_str());
  pt2->SetTextColor(4);
  pt2->SetFillColor(10);
  pt2->SetBorderSize(0);
  pt2->Draw();  
 
  /**********************************************************************************************/

  TMultiGraph *mg1 = new TMultiGraph();
  TMultiGraph *mg2 = new TMultiGraph();
  //---------------------- fill vectors and histogram --------------------------------------------------//
  TList *list = new TList;
  Int_t numkey;
  TObject *key = new TObject;
  const char *name;
  char *SoftInfo="SoftInfo";
  TTree* tr = new TTree;
  Long64_t nevntskey=0;
  list = file->GetListOfKeys(); 
  numkey = file->GetNkeys();
  ULong_t salto;
  for (Int_t m=0; m<numkey; m++){
    key = list->At(m);
    name=(char *)(key->GetName());
    if(strcmp(name,SoftInfo)==0)continue;
    tr = (TTree*)file->Get(name);
    if (tr->IsZombie()) continue;

    tr->SetBranchAddress("Header", &eh);
    TBranch *Br = tr->GetBranch("Header");
    nevntskey = tr->GetEntries();
    if(nevntskey !=0){
      Int_t size1=nevntskey;
      Double_t *PscuCounter1 = new Double_t[size1];
      Double_t *FileOffset1 = new Double_t[size1];
      Double_t *tempo1 = new Double_t[size1];

      int n=0;
      while(n<nevntskey){
        Br->GetEntry(n);
        ph = eh->GetPscuHeader();  
        PscuCounter1[n]= ph->GetCounter();
        FileOffset1[n]=ph->GetFileOffset();
        tempo1[n]=ph->GetOrbitalTime();
        if((m==0) && (n==0)){
          primotempo=ph->GetOrbitalTime();
          salto=ph->GetOrbitalTime();
        }
        if(salto > ph->GetOrbitalTime())salto=ph->GetOrbitalTime();
        if(ph->GetFileOffset()<= primoffset){
          primoffset=ph->GetFileOffset();
          primotempo=ph->GetOrbitalTime();
        }
        if(ph->GetFileOffset()>=ultimoffset){
          ultimotempo=ph->GetOrbitalTime();
          ultimoffset=ph->GetFileOffset();
        }
        n++;
      }

      TGraph *graph3= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)PscuCounter1);
      graph3->SetMarkerColor(3);
      graph3->SetMarkerSize(.2);
      graph3->SetMarkerStyle(21);
      mg1->Add(graph3);
      
      TGraph *graph4= new TGraph(nevntskey, (const Double_t*)FileOffset1, (const Double_t*)tempo1);
      graph4->SetMarkerColor(kBlue);
      graph4->SetMarkerSize(.2);
      graph4->SetMarkerStyle(21);
      mg2->Add(graph4);
    }
  }

  TLatex *lat = new TLatex();
  lat->SetTextAlign(12);
  lat->SetTextSize(0.15);
  lat->SetTextColor(kBlue);

 //------------ Fourth PAD  ---------------------//
   finalCanv->cd(4);
   all1->Draw();
   all1->cd();

   oss1.str("");
   oss1 <<"PscuCounter vs FileOffset.";
   mg1->Draw("AP");
   mg1->GetXaxis()->SetTitle("File Offset");
   mg1->GetXaxis()->CenterTitle();
   mg1->GetXaxis()->SetTitleOffset(0.8);
   mg1->GetXaxis()->SetTitleSize(0.05);
   mg1->GetXaxis()->SetLabelSize(0.05);
   mg1->GetYaxis()->SetTitle("Counter");
   mg1->GetYaxis()->CenterTitle();
   mg1->GetYaxis()->SetTitleSize(0.06);
   mg1->GetYaxis()->SetLabelSize(0.06);
   mg1->GetYaxis()->SetTitleOffset(0.6);
   finalCanv->cd(4);
   TPaveText *pt3=0;
   pt3 = new TPaveText (.60,.91,.90,1);
   pt3->AddText(oss1.str().c_str());
   pt3->SetTextColor(3);
   pt3->SetFillColor(10);
   pt3->SetBorderSize(0);
   pt3->Draw();

 //------------ Fifth PAD  ---------------------//
   finalCanv->cd(5);
   all->Draw();
   all->cd();
   oss3.str("");
   oss3 << "OBT vs FileOffset";
   mg2->Draw("AP");
   mg2->GetXaxis()->SetTitle("File Offset");
   mg2->GetXaxis()->CenterTitle();
   mg2->GetXaxis()->SetTitleSize(0.05);
   mg2->GetXaxis()->SetLabelSize(0.05);
   mg2->GetYaxis()->SetTitle("OBT");
   mg2->GetYaxis()->CenterTitle();
   mg2->GetYaxis()->SetTitleSize(0.06);
   mg2->GetYaxis()->SetLabelSize(0.05);
   mg2->GetYaxis()->SetTitleOffset(0.6);

   finalCanv->cd(5);
   TPaveText *pt4=0;
   pt4 = new TPaveText (.70,.91,.90,1);
   pt4->AddText(oss3.str().c_str());
   pt4->SetTextColor(kBlue);
   pt4->SetFillColor(10);
   pt4->SetBorderSize(0);
   pt4->Draw();

   finalCanv->cd(6);
   ULong_t primotempoABS=TimeSync+((primotempo/1000)-ObtSync);
   ULong_t obmaxABS=TimeSync+((lastime/1000)-ObtSync);
   ULong_t saltoABS=TimeSync+((salto/1000)-ObtSync);
   ULong_t ultimotempoABS=TimeSync+((ultimotempo/1000)-ObtSync);

   TPaveText *pt5=0;
   pt5 = new TPaveText (0,0,1,1);
   stringa.str("");
   stringa << "  Filename:   "<<filename.Data()<<"\n";
   TText *t1=pt5->AddText(0.25,0.95,stringa.str().c_str());
   t1->SetTextSize(0.1);
   stringa.str("");
   stringa << "  OBT (ms)            ABS TIME (s)";
   TText *t2=pt5->AddText(0.32,0.75,stringa.str().c_str());
   t2->SetTextSize(0.07);
   stringa.str("");
   stringa << "New data start at:   "<<primotempo<<"         "<<primotempoABS;
   TText *t3=pt5->AddText(0.25,0.60,stringa.str().c_str());
   t3->SetTextSize(0.08);
   stringa.str("");
   stringa << "New data end at:     "<<lastime<<"         "<<obmaxABS;
   TText *t4=pt5->AddText(0.25,0.50,stringa.str().c_str());
   t4->SetTextSize(0.08);
   if(primotempo!=salto || lastime!=ultimotempo){
     stringa.str("");
     stringa << "Old data start at:   "<<salto<<"        "<<saltoABS;
     TText *t5=pt5->AddText(0.65,0.60,stringa.str().c_str());
     t5->SetTextSize(0.08);
     stringa.str("");
     stringa << "Old data end at:     "<<ultimotempo<<"        "<<ultimotempoABS;
     TText *t6=pt5->AddText(0.65,0.50,stringa.str().c_str());
     t6->SetTextSize(0.08);
     stringa.str("");
     stringa << "  OBT (ms)          ABS TIME (s)";
     TText *t2=pt5->AddText(0.72,0.75,stringa.str().c_str());
     t2->SetTextSize(0.07);
   }
   pt5->SetTextColor(kBlack);
   pt5->SetFillColor(10);
   pt5->SetBorderSize(0);
   pt5->Draw();

   finalCanv->Update();

   oss1.str("");
   oss1 << outDir.Data() << filename.Data();
   oss1 << "_HeaderScan"<<"." << format.Data();

   finalCanv->SaveAs(oss1.str().c_str());

   file->Close();
  
}

int main(int argc, char* argv[]){
  TString path;
  TString outDir     = "./";
  TString format     = "jpg";
 if (argc < 2){
    printf("You have to insert at least the file to analyze \n");
    printf("Try '--help' for more information. \n");
    exit(1);
  }  
  if (!strcmp(argv[1], "--help")){
        printf( "Usage: HeaderScan FILE [OPTION] \n");
        printf( "\t --help                 Print this help and exit \n");
        printf( "\t -outDir[path]          Path where to put the output [default ./] \n");
        printf( "\t -format[jpg|ps|gif]    Format for output files [default 'jpg'] \n");
        exit(1);
  }
  path=argv[1];
  for (int i = 2; i < argc; i++){ 
    if (!strcmp(argv[i], "-outDir")){
      if (++i >= argc){
        printf( "-outDir needs arguments. \n");
        printf( "Try '--help' for more information. \n");
        exit(1);
      } 
      else{
        outDir = argv[i];
        continue;
      }
    }  
    if (!strcmp(argv[i], "-format")){
      if (++i >= argc){
        printf( "-format needs arguments. \n");
        printf( "Try '--help' for more information. \n");
        exit(1);
      }
      else{
        format = argv[i];
        continue;
      }
    }
  }
  HeaderScan(argv[1], outDir, format);
}

1	pam-rm2	1.1	/**
2	pam-rm2	1.4	* Header_Scan
3	pam-rm2	1.1	* Author Nagni
4			* version 1.0
5			*
6			* Version 1.1 - 28 December 2004
7			* If outList does not exist the function exit to prompt
8			* If nevents = 0 the function exit to promt
9			*
10			* Version 1.2 - 3 January 2005
11			* Two canvases are created to see the graphs better
12			*
13			* Version 1.25 - 13 January 2005
14			* Changed Int_t to Float because of variable range size
15			* (UInt_t has been excluded beacuse of uncompatibility with TGraph)
16			*
17			* version 1.3 - 22 February 2005 - Nagni
18			* For compatibility with batch mode excution:
19			* 1) Added "include <iostream>" and "using namespace std"
20			* 2) Removed gROOT->Reset()
21			*
22			* Version 1.4
23			* Date 08 March 2005
24			* Added "format" parameter for saving the produced image in various formats
25			* (for a complete list of types refer to TPad::SaveAs method)
26			*
27			* Version 1.5
28			* Date 09 February 2006 - Marcelli
29			* Update to work with new Yoda output
30			*
31			* Version 1.6
32			* Date 27 February 2006 - Marcelli
33			* For compilation:
34			* Added function "int main(int argc, char* argv[])"
35			*
36			*
37			* Description: This script creates two canvases with five pads. The first pad shows packetID variable (for all packets) vs. OBT.
38			* The second pad shows the number of physic packets vs. OBT. The third pad shows the lenght of Physic packets (byte) vs. OBT.
39			* The fourth pad shows the packetcounter of physic packets vs. OBT. The fifth pad shows PacketCounter vs. File Offset.
40			*
41			* Parameters:
42			* TSTring base - the path to the root directory for the specific Pamela unpack session
43			* There is no default value, without this input the program will not run
44			* TString outDir - the path where to save the output image (Default = ./)
45			* TString format - the format which will be used fo rsave the produced images (Default = "jpg")
46	pam-rm2	1.2	*
47			*
48			* Version 1.7
49			* Date 16 June 2006 - Malvezzi
50			*
51			* Description of changes:
52			* Implementation of the case: numebr of events <= 0.
53			* Remove graph "grPcktId1"; see PacketScan for the same information.
54			* Fixed bugs: for a large namber of events is not possible to have vectors, so I have subsituted graphs with histograms
55			* or divided the graphs in more than one canvas.
56			*
57	pam-rm2	1.6	* 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	pam-rm2	1.1	*/
63
64
65	pam-rm2	1.2	#include <fstream>
66			#include <math.h>
67			#include "TLatex.h"
68			#include "TF1.h"
69			#include "TPaveText.h"
70			#include "TMultiGraph.h"
71	pam-rm2	1.1	#include <sstream>
72			#include <iostream>
73	pam-rm2	1.11	#include <stdio.h>
74			#include <string.h>
75	pam-rm2	1.1	#include "TString.h"
76			#include "TStyle.h"
77			#include "TFile.h"
78			#include "TList.h"
79			#include "TTree.h"
80			#include "TObjString.h"
81			#include "TCanvas.h"
82			#include "TGraph.h"
83			#include "TH1F.h"
84			#include "TGaxis.h"
85			#include "EventHeader.h"
86			#include "PscuHeader.h"
87	pam-rm2	1.8	#include "RunHeaderEvent.h"
88			#include "TPaveText.h"
89	pam-rm2	1.1
90			using namespace std;
91
92	pam-rm2	1.9	void HeaderScan(TString base, TString outDir, TString format){
93	pam-rm2	1.1
94	pam-rm2	1.2	//------------------- Variables initilization -------------------------//
95	pam-rm2	1.8	Long64_t nevents=0, runnevents;
96	pam-rm2	1.7	ULong_t lastime, firstime, primotempo, ultimotempo, primoffset=500000000, ultimoffset;
97	pam-rm2	1.9	double obmin=0.;
98	pam-rm2	1.2	double obmax=0.;
99	pam-rm2	1.6	stringstream oss, oss1, oss2, oss3, noentries, stringa;
100	pam-rm2	1.9
101			//-------------- Load root file, tree and branch -------------------//
102	pam-rm2	1.1	TFile *file = new TFile(base.Data());
103			if (!file){
104			printf("No such file in the directory has been found");
105			return;
106	pam-rm2	1.9	}
107
108	pam-rm2	1.2	TTree PhysicsTr = (TTree)file->Get("Physics");
109			TBranch *headBr = PhysicsTr->GetBranch("Header");
110
111	pam-rm2	1.11	pamela::EventHeader *eh = new pamela::EventHeader;
112			pamela::PscuHeader *ph = new pamela::PscuHeader;
113	pam-rm2	1.2
114			PhysicsTr->SetBranchAddress("Header", &eh);
115	pam-rm2	1.9	nevents = PhysicsTr->GetEntries();
116			const Int_t size = nevents;
117	pam-rm2	1.11
118	pam-rm2	1.8	TTree RunHeadTr = (TTree)file->Get("RunHeader"); ///run header tree
119	pam-rm2	1.11	pamela::EventHeader *eH= new pamela::EventHeader;
120			pamela::RunHeaderEvent *reh=new pamela::RunHeaderEvent;
121	pam-rm2	1.8
122			RunHeadTr->SetBranchAddress("Header",&eH);
123			RunHeadTr->SetBranchAddress("RunHeader",&reh);
124			runnevents = RunHeadTr->GetEntries();
125	pam-rm2	1.11
126	pam-rm2	1.1	TString filename = ((TObjString*)base.Tokenize('/')->Last())->GetString();
127			filename = ((TObjString*)filename.Tokenize('.')->First())->GetString();
128	pam-rm2	1.2
129			//----------- If nevents < = 0 ---------------------------------/
130			if (nevents<=0) {
131			printf("nevents = %i \n", nevents);
132			printf(" \n");
133			TCanvas *canv = new TCanvas("No entries", "No entries ", 400, 200);
134			canv->SetFillColor(10);
135			canv->cd();
136			TLatex *l = new TLatex();
137			l->SetTextAlign(12);
138			l->SetTextSize(0.15);
139			l->SetTextColor(2);
140			noentries.str("");
141			noentries<< "HeaderScan_QL:";
142			l->DrawLatex(0.05, 0.7, noentries.str().c_str());
143			noentries.str("");
144			noentries<< "No Physics entries for this files";
145			l->DrawLatex(0.05, 0.5, noentries.str().c_str());
146			if (outDir == "./") {
147			oss.str("");
148			oss << filename.Data() << "_HeaderScan_QL." << format.Data();
149			} else {
150			oss.str("");
151			oss << outDir.Data() << filename.Data() << "_HeaderScan_QL." << format.Data();
152	pam-rm2	1.9	}
153	pam-rm2	1.2	canv->Update();
154			canv->SaveAs(oss.str().c_str());
155			return;
156	pam-rm2	1.1	}
157	pam-rm2	1.2	//-------------- to know the max and min OBT ----------------------------//
158	pam-rm2	1.5	headBr->GetEntry(0);
159			ph = eh->GetPscuHeader();
160			firstime = ph->GetOrbitalTime();
161	pam-rm2	1.9	headBr->GetEntry(nevents);
162			ph = eh->GetPscuHeader();
163			lastime = ph->GetOrbitalTime();
164			int i =0;
165			while(i < nevents){
166			headBr->GetEntry(i);
167	pam-rm2	1.8	ph = eh->GetPscuHeader();
168			if((ph->GetOrbitalTime()) <= firstime) firstime=ph->GetOrbitalTime();
169			if((ph->GetOrbitalTime()) >= lastime) lastime=ph->GetOrbitalTime();
170	pam-rm2	1.9	i++;
171	pam-rm2	1.8	}
172	pam-rm2	1.4
173	pam-rm2	1.2	//------------------------ First histogram -----------------------------------//
174	pam-rm2	1.8	obmin=firstime;
175	pam-rm2	1.9	obmax=lastime;
176			Int_t nbin = (lastime-firstime)/60000;
177			TH1F *h1 = new TH1F ("histo1", "" , nbin, obmin, obmax);
178
179			//------------------- fill vectors and histogram -----------------------------//
180			Double_t *PscuCounter = new Double_t[size];
181			Double_t *PacketLenght = new Double_t[size];
182			Double_t *OBTime = new Double_t[size];
183			Double_t *Eventsperminute= new Double_t[nbin];
184			Double_t *Minute= new Double_t[nbin];
185			Int_t max=0;
186			for (Int_t k = 0; k < nevents; k++){
187			headBr->GetEntry(k);
188			ph = eh->GetPscuHeader();
189			h1->Fill(ph->GetOrbitalTime());
190			PscuCounter[k]= ph->GetCounter();
191			PacketLenght[k]=ph->GetPacketLenght();
192			OBTime[k]=ph->GetOrbitalTime();
193			}
194			int l=0;
195			while(l<nbin){
196			Eventsperminute[l]=h1->GetBinContent(l);
197			Minute[l]=firstime+l*60000;
198			if(h1->GetBinContent(l) >= max)max =(Int_t)h1->GetBinContent(l);
199			l++;
200			}
201
202			//----------- Graph and MultiGraph -----------------------------------------------//
203			TMultiGraph *rate = new TMultiGraph();
204			TMultiGraph *packetLength = new TMultiGraph();
205			TMultiGraph *packeCounter = new TMultiGraph();
206
207	pam-rm2	1.2	oss1.str("");
208	pam-rm2	1.10	oss1 << "Physics Packet per minute. Start time = " << obmin << ", End time = "<< obmax <<" ms";
209	pam-rm2	1.9	TGraph rate1= new TGraph(nbin, (const Double_t)Minute, (const Double_t*)Eventsperminute);
210			rate1->SetMarkerColor(kBlack);
211			rate1->SetMarkerSize(.1);
212			rate1->SetMarkerStyle(21);
213			rate->Add(rate1);
214
215			TGraph packetLength1= new TGraph(nevents, (const Double_t)OBTime, (const Double_t*)PacketLenght);
216	pam-rm2	1.2	oss2.str("");
217	pam-rm2	1.10	oss2 <<"Lenght of Physic packets";
218	pam-rm2	1.9	packetLength1->SetMarkerColor(2);
219	pam-rm2	1.10	packetLength1->SetMarkerSize(.3);
220	pam-rm2	1.9	packetLength1->SetMarkerStyle(21);
221			packetLength->Add(packetLength1);
222	pam-rm2	1.2
223	pam-rm2	1.9	TGraph packeCounter1= new TGraph(nevents, (const Double_t)OBTime, (const Double_t*)PscuCounter);
224	pam-rm2	1.2	oss3.str("");
225	pam-rm2	1.10	oss3 <<"Physics Counter vs. OBT";
226	pam-rm2	1.9	packeCounter1->SetMarkerColor(4);
227			packeCounter1->SetMarkerSize(.2);
228			packeCounter1->SetMarkerStyle(21);
229			packeCounter->Add(packeCounter1);
230
231			//------------ Create and Draw Canvas ---------------------//
232	pam-rm2	1.8	TCanvas *finalCanv = new TCanvas("Header", base, 1200, 1600);
233	pam-rm2	1.10	finalCanv->Divide(1,6);
234	pam-rm2	1.8	finalCanv->SetFillColor(10);
235	pam-rm2	1.9
236	pam-rm2	1.6	TPad *all2= new TPad ("","", 0, 0, 1, 1);
237			all2->SetFillColor(10);
238			TPad *all3= new TPad ("","", 0, 0, 1, 1);
239			all3->SetFillColor(10);
240			TPad *all4= new TPad ("","", 0, 0, 1, 1);
241			all4->SetFillColor(10);
242	pam-rm2	1.8	TPad *all= new TPad ("","", 0, 0, 1, 1);
243			all->SetFillColor(10);
244			TPad *all1= new TPad ("","", 0, 0, 1, 1);
245			all1->SetFillColor(10);
246	pam-rm2	1.6
247	pam-rm2	1.8	TLine li;
248			li.SetLineStyle(4);
249			li.SetLineWidth(2);
250	pam-rm2	1.9
251			//----------- First PAD -------------------------------//
252	pam-rm2	1.8	finalCanv->cd(1);
253	pam-rm2	1.6	all2->Draw();
254			all2->cd();
255	pam-rm2	1.9	rate->Draw("ALP");
256			rate->GetXaxis()->SetTitle("OBT (ms)");
257			rate->GetXaxis()->SetTitleSize(0.05);
258			rate->GetXaxis()->CenterTitle();
259			rate->GetXaxis()->SetLabelSize(0.05);
260			rate->GetYaxis()->SetTitle("Number of events ");
261			rate->GetYaxis()->CenterTitle();
262			rate->GetYaxis()->SetLabelSize(0.05);
263			rate->GetYaxis()->SetTitleSize(0.06);
264			rate->GetYaxis()->SetTitleOffset(0.6);
265	pam-rm2	1.8	for (Int_t l = 0; l < runnevents; l++){
266			RunHeadTr->GetEntry(l);
267			ph = eH->GetPscuHeader();
268			int ws= reh->RM_ACQ_SETTING_MODE;
269			int id = ph->GetPacketId1();
270			Int_t obt = ph->GetOrbitalTime();
271			if (ws==1){
272			li.SetLineColor(3);
273	pam-rm2	1.9	li.DrawLine(obt,0,obt,max);
274	pam-rm2	1.8	}
275			else if (ws==2){
276			li.SetLineColor(4);
277	pam-rm2	1.9	li.DrawLine(obt,0,obt,max);
278			}
279			}
280
281	pam-rm2	1.10	RunHeadTr->GetEntry(0);
282			ph = eH->GetPscuHeader();
283			ULong_t TimeSync = reh->LAST_TIME_SYNC_INFO;
284			ULong_t ObtSync = reh->OBT_TIME_SYNC;
285			//cout<<"TimeSync "<<reh->LAST_TIME_SYNC_INFO<<"\n";
286			//cout<<"ObtSync "<<reh->OBT_TIME_SYNC<<"\n";
287
288	pam-rm2	1.8	finalCanv->cd(1);
289			stringstream ws1, ws2;
290			ws1.str("");
291			ws2.str("");
292			ws1<<"ACQ_SETTING_MODE = 1";
293			ws2<<"ACQ_SETTING_MODE = 2";
294			TPaveText *pt=0;
295			pt = new TPaveText (.60,.92,.76,.98);
296			pt->AddText(ws1.str().c_str());
297			pt->SetTextColor(3);
298			pt->SetFillColor(10);
299			pt->SetBorderSize(0);
300			pt->Draw();
301			TPaveText *pt1=0;
302			pt1 = new TPaveText (.76,.92,.92,.98);
303			pt1->AddText(ws2.str().c_str());
304			pt1->SetTextColor(4);
305			pt1->SetFillColor(10);
306			pt1->SetBorderSize(0);
307			pt1->Draw();
308			pt1 = new TPaveText (.05,.91,.6,1);
309			pt1->AddText(oss1.str().c_str());
310			pt1->SetTextColor(1);
311			pt1->SetFillColor(10);
312			pt1->SetBorderSize(0);
313			pt1->Draw();
314	pam-rm2	1.6
315	pam-rm2	1.9	//----------- Second PAD -------------------------------//
316	pam-rm2	1.8	finalCanv->cd(2);
317	pam-rm2	1.6	all3->Draw();
318			all3->cd();
319	pam-rm2	1.9	packetLength->Draw("AP");
320	pam-rm2	1.6	packetLength->GetXaxis()->SetTitle("OBT (ms)");
321	pam-rm2	1.9	packetLength->GetXaxis()->SetTitleSize(0.05);
322	pam-rm2	1.6	packetLength->GetXaxis()->CenterTitle();
323	pam-rm2	1.8	packetLength->GetXaxis()->SetLabelSize(0.05);
324	pam-rm2	1.6	packetLength->GetYaxis()->SetTitle("Lenght (byte)");
325			packetLength->GetYaxis()->CenterTitle();
326	pam-rm2	1.8	packetLength->GetYaxis()->SetLabelSize(0.05);
327	pam-rm2	1.6	packetLength->GetYaxis()->SetTitleSize(0.06);
328	pam-rm2	1.8	packetLength->GetYaxis()->SetTitleOffset(0.6);
329	pam-rm2	1.9
330
331	pam-rm2	1.8	finalCanv->cd(2);
332			pt = new TPaveText (.6,.91,.90,1);
333			pt->AddText(oss2.str().c_str());
334			pt->SetTextColor(2);
335			pt->SetFillColor(10);
336			pt->SetBorderSize(0);
337			pt->Draw();
338	pam-rm2	1.2
339	pam-rm2	1.9	//----------- Third PAD -------------------------------//
340	pam-rm2	1.8	finalCanv->cd(3);
341	pam-rm2	1.6	all4->Draw();
342			all4->cd();
343	pam-rm2	1.9	packeCounter->Draw("AP");
344	pam-rm2	1.6	packeCounter->GetXaxis()->SetTitle("OBT (ms)");
345			packeCounter->GetXaxis()->SetTitleSize(0.05);
346			packeCounter->GetXaxis()->CenterTitle();
347	pam-rm2	1.8	packeCounter->GetXaxis()->SetLabelSize(0.05);
348	pam-rm2	1.6	packeCounter->GetYaxis()->SetTitle("Counter");
349			packeCounter->GetYaxis()->SetTitleSize(0.05);
350			packeCounter->GetYaxis()->CenterTitle();
351	pam-rm2	1.8	packeCounter->GetYaxis()->SetLabelSize(0.05);
352	pam-rm2	1.6	packeCounter->GetYaxis()->SetTitleSize(0.06);
353	pam-rm2	1.8	packeCounter->GetYaxis()->SetTitleOffset(0.6);
354	pam-rm2	1.9
355
356	pam-rm2	1.8	finalCanv->cd(3);
357			TPaveText *pt2=0;
358			pt2 = new TPaveText (.6,.91,.90,1);
359			pt2->AddText(oss3.str().c_str());
360			pt2->SetTextColor(4);
361			pt2->SetFillColor(10);
362			pt2->SetBorderSize(0);
363	pam-rm2	1.9	pt2->Draw();
364
365			/**********************************************************************************************/
366
367			TMultiGraph *mg1 = new TMultiGraph();
368			TMultiGraph *mg2 = new TMultiGraph();
369			//---------------------- fill vectors and histogram --------------------------------------------------//
370	pam-rm2	1.6	TList *list = new TList;
371			Int_t numkey;
372			TObject *key = new TObject;
373			const char *name;
374	pam-rm2	1.11	char *SoftInfo="SoftInfo";
375	pam-rm2	1.6	TTree* tr = new TTree;
376			Long64_t nevntskey=0;
377			list = file->GetListOfKeys();
378			numkey = file->GetNkeys();
379	pam-rm2	1.10	ULong_t salto;
380	pam-rm2	1.9	for (Int_t m=0; m<numkey; m++){
381			key = list->At(m);
382	pam-rm2	1.6	name=(char *)(key->GetName());
383	pam-rm2	1.11	if(strcmp(name,SoftInfo)==0)continue;
384	pam-rm2	1.6	tr = (TTree*)file->Get(name);
385			if (tr->IsZombie()) continue;
386	pam-rm2	1.11
387	pam-rm2	1.6	tr->SetBranchAddress("Header", &eh);
388			TBranch *Br = tr->GetBranch("Header");
389			nevntskey = tr->GetEntries();
390	pam-rm2	1.8	if(nevntskey !=0){
391			Int_t size1=nevntskey;
392			Double_t *PscuCounter1 = new Double_t[size1];
393			Double_t *FileOffset1 = new Double_t[size1];
394			Double_t *tempo1 = new Double_t[size1];
395	pam-rm2	1.9
396			int n=0;
397			while(n<nevntskey){
398			Br->GetEntry(n);
399	pam-rm2	1.8	ph = eh->GetPscuHeader();
400	pam-rm2	1.9	PscuCounter1[n]= ph->GetCounter();
401			FileOffset1[n]=ph->GetFileOffset();
402			tempo1[n]=ph->GetOrbitalTime();
403			if((m==0) && (n==0)){
404			primotempo=ph->GetOrbitalTime();
405			salto=ph->GetOrbitalTime();
406			}
407			if(salto > ph->GetOrbitalTime())salto=ph->GetOrbitalTime();
408	pam-rm2	1.8	if(ph->GetFileOffset()<= primoffset){
409	pam-rm2	1.9	primoffset=ph->GetFileOffset();
410	pam-rm2	1.8	primotempo=ph->GetOrbitalTime();
411			}
412	pam-rm2	1.9	if(ph->GetFileOffset()>=ultimoffset){
413			ultimotempo=ph->GetOrbitalTime();
414			ultimoffset=ph->GetFileOffset();
415			}
416			n++;
417	pam-rm2	1.8	}
418	pam-rm2	1.11
419	pam-rm2	1.8	TGraph graph3= new TGraph(nevntskey, (const Double_t)FileOffset1, (const Double_t*)PscuCounter1);
420			graph3->SetMarkerColor(3);
421			graph3->SetMarkerSize(.2);
422			graph3->SetMarkerStyle(21);
423			mg1->Add(graph3);
424
425			TGraph graph4= new TGraph(nevntskey, (const Double_t)FileOffset1, (const Double_t*)tempo1);
426			graph4->SetMarkerColor(kBlue);
427			graph4->SetMarkerSize(.2);
428			graph4->SetMarkerStyle(21);
429			mg2->Add(graph4);
430	pam-rm2	1.6	}
431			}
432	pam-rm2	1.8
433			TLatex *lat = new TLatex();
434			lat->SetTextAlign(12);
435			lat->SetTextSize(0.15);
436			lat->SetTextColor(kBlue);
437	pam-rm2	1.9
438			//------------ Fourth PAD ---------------------//
439	pam-rm2	1.8	finalCanv->cd(4);
440	pam-rm2	1.6	all1->Draw();
441			all1->cd();
442	pam-rm2	1.8
443			oss1.str("");
444	pam-rm2	1.10	oss1 <<"PscuCounter vs FileOffset.";
445	pam-rm2	1.3	mg1->Draw("AP");
446			mg1->GetXaxis()->SetTitle("File Offset");
447			mg1->GetXaxis()->CenterTitle();
448	pam-rm2	1.9	mg1->GetXaxis()->SetTitleOffset(0.8);
449			mg1->GetXaxis()->SetTitleSize(0.05);
450	pam-rm2	1.8	mg1->GetXaxis()->SetLabelSize(0.05);
451	pam-rm2	1.3	mg1->GetYaxis()->SetTitle("Counter");
452			mg1->GetYaxis()->CenterTitle();
453	pam-rm2	1.9	mg1->GetYaxis()->SetTitleSize(0.06);
454	pam-rm2	1.8	mg1->GetYaxis()->SetLabelSize(0.06);
455	pam-rm2	1.9	mg1->GetYaxis()->SetTitleOffset(0.6);
456	pam-rm2	1.8	finalCanv->cd(4);
457			TPaveText *pt3=0;
458			pt3 = new TPaveText (.60,.91,.90,1);
459			pt3->AddText(oss1.str().c_str());
460			pt3->SetTextColor(3);
461			pt3->SetFillColor(10);
462			pt3->SetBorderSize(0);
463			pt3->Draw();
464
465	pam-rm2	1.9	//------------ Fifth PAD ---------------------//
466	pam-rm2	1.8	finalCanv->cd(5);
467	pam-rm2	1.6	all->Draw();
468			all->cd();
469			oss3.str("");
470	pam-rm2	1.10	oss3 << "OBT vs FileOffset";
471	pam-rm2	1.4	mg2->Draw("AP");
472			mg2->GetXaxis()->SetTitle("File Offset");
473			mg2->GetXaxis()->CenterTitle();
474	pam-rm2	1.9	mg2->GetXaxis()->SetTitleSize(0.05);
475	pam-rm2	1.8	mg2->GetXaxis()->SetLabelSize(0.05);
476	pam-rm2	1.4	mg2->GetYaxis()->SetTitle("OBT");
477			mg2->GetYaxis()->CenterTitle();
478	pam-rm2	1.9	mg2->GetYaxis()->SetTitleSize(0.06);
479			mg2->GetYaxis()->SetLabelSize(0.05);
480	pam-rm2	1.8	mg2->GetYaxis()->SetTitleOffset(0.6);
481	pam-rm2	1.6
482	pam-rm2	1.8	finalCanv->cd(5);
483			TPaveText *pt4=0;
484	pam-rm2	1.10	pt4 = new TPaveText (.70,.91,.90,1);
485	pam-rm2	1.8	pt4->AddText(oss3.str().c_str());
486			pt4->SetTextColor(kBlue);
487			pt4->SetFillColor(10);
488			pt4->SetBorderSize(0);
489			pt4->Draw();
490	pam-rm2	1.6
491	pam-rm2	1.10	finalCanv->cd(6);
492			ULong_t primotempoABS=TimeSync+((primotempo/1000)-ObtSync);
493			ULong_t obmaxABS=TimeSync+((lastime/1000)-ObtSync);
494			ULong_t saltoABS=TimeSync+((salto/1000)-ObtSync);
495			ULong_t ultimotempoABS=TimeSync+((ultimotempo/1000)-ObtSync);
496
497			TPaveText *pt5=0;
498			pt5 = new TPaveText (0,0,1,1);
499			stringa.str("");
500			stringa << " Filename: "<<filename.Data()<<"\n";
501			TText *t1=pt5->AddText(0.25,0.95,stringa.str().c_str());
502			t1->SetTextSize(0.1);
503			stringa.str("");
504			stringa << " OBT (ms) ABS TIME (s)";
505			TText *t2=pt5->AddText(0.32,0.75,stringa.str().c_str());
506			t2->SetTextSize(0.07);
507			stringa.str("");
508			stringa << "New data start at: "<<primotempo<<" "<<primotempoABS;
509			TText *t3=pt5->AddText(0.25,0.60,stringa.str().c_str());
510			t3->SetTextSize(0.08);
511			stringa.str("");
512			stringa << "New data end at: "<<lastime<<" "<<obmaxABS;
513			TText *t4=pt5->AddText(0.25,0.50,stringa.str().c_str());
514			t4->SetTextSize(0.08);
515			if(primotempo!=salto \|\| lastime!=ultimotempo){
516			stringa.str("");
517			stringa << "Old data start at: "<<salto<<" "<<saltoABS;
518			TText *t5=pt5->AddText(0.65,0.60,stringa.str().c_str());
519			t5->SetTextSize(0.08);
520			stringa.str("");
521			stringa << "Old data end at: "<<ultimotempo<<" "<<ultimotempoABS;
522			TText *t6=pt5->AddText(0.65,0.50,stringa.str().c_str());
523			t6->SetTextSize(0.08);
524			stringa.str("");
525			stringa << " OBT (ms) ABS TIME (s)";
526			TText *t2=pt5->AddText(0.72,0.75,stringa.str().c_str());
527			t2->SetTextSize(0.07);
528			}
529			pt5->SetTextColor(kBlack);
530			pt5->SetFillColor(10);
531			pt5->SetBorderSize(0);
532			pt5->Draw();
533
534	pam-rm2	1.8	finalCanv->Update();
535	pam-rm2	1.6
536	pam-rm2	1.2	oss1.str("");
537			oss1 << outDir.Data() << filename.Data();
538	pam-rm2	1.8	oss1 << "_HeaderScan"<<"." << format.Data();
539
540			finalCanv->SaveAs(oss1.str().c_str());
541	pam-rm2	1.2
542	pam-rm2	1.3	file->Close();
543	pam-rm2	1.1
544			}
545
546			int main(int argc, char* argv[]){
547			TString path;
548			TString outDir = "./";
549			TString format = "jpg";
550			if (argc < 2){
551			printf("You have to insert at least the file to analyze \n");
552			printf("Try '--help' for more information. \n");
553			exit(1);
554			}
555			if (!strcmp(argv[1], "--help")){
556			printf( "Usage: HeaderScan FILE [OPTION] \n");
557			printf( "\t --help Print this help and exit \n");
558			printf( "\t -outDir[path] Path where to put the output [default ./] \n");
559			printf( "\t -format[jpg\|ps\|gif] Format for output files [default 'jpg'] \n");
560			exit(1);
561			}
562			path=argv[1];
563	pam-rm2	1.2	for (int i = 2; i < argc; i++){
564	pam-rm2	1.1	if (!strcmp(argv[i], "-outDir")){
565			if (++i >= argc){
566			printf( "-outDir needs arguments. \n");
567			printf( "Try '--help' for more information. \n");
568			exit(1);
569			}
570			else{
571			outDir = argv[i];
572			continue;
573			}
574	pam-rm2	1.2	}
575	pam-rm2	1.1	if (!strcmp(argv[i], "-format")){
576			if (++i >= argc){
577			printf( "-format needs arguments. \n");
578			printf( "Try '--help' for more information. \n");
579			exit(1);
580			}
581			else{
582			format = argv[i];
583			continue;
584			}
585			}
586			}
587			HeaderScan(argv[1], outDir, format);
588			}
589	pam-rm2	1.7