START/readDigits.C

   1 void readDigits(Int_t evNumber=1)
   2 {
   3
   4   // Dynamically link some shared libs
   5   if (gClassTable->GetID("AliRun") < 0) {
   6     gROOT->LoadMacro("loadlibs.C");
   7     loadlibs();
   8   }
   9
  10   // Connect the Root Galice file containing Geometry, Kine and Hits
  11   TFile *file =  (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
  12   //TFile *file =  (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
  13   if (!file) file = new TFile("production.root","UPDATE");
  14
  15   // Get AliRun object from file or create it if not on file
  16   if (!gAlice) {
  17     gAlice = (AliRun*)file->Get("gAlice");
  18     if (gAlice) printf("AliRun object found on file\n");
  19     if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
  20   }
  21   char nameTD[8],nameTR[8];
  22
  23   TH1F *hTimediff = new TH1F("hTimediff","Time difference",100,0,200);
  24   TH1F *hTimePs = new TH1F("hTimePs","Time in Ps",100,2000,3000);
  25
  26   digits = new AliSTARTdigit();
  27
  28
  29  // Event ------------------------- LOOP
  30   for (j=0; j<evNumber; j++){
  31     gAlice->GetEvent(j);
  32     sprintf(nameTD,"START_D_%d",j);
  33     printf("%s\n",nameTD);
  34   TObject *td = (TObject*)gDirectory->Get(nameTD);
  35   digits->Read(nameTD);
  36   digits->Dump();
  37   printf("time %d\n",digits->GetTime());
  38
  39     if(digits->GetTime()!=999999){
  40       Int_t timediff = digits->GetMeanTime();
  41       //     Double_t timePs=(timediff-128)*10.; // time in Ps channel_width =10ps
  42       Double_t timePs=(timediff)*10.; // time in Ps channel_width =10ps
  43       cout<<"timediff "<<timediff<<" timePs "<<timePs<<endl;
  44       hTimediff->Fill(timediff);
  45       hTimePs->Fill(timePs);
  46     }
  47   }
  48
  49       Hfile = new TFile("figs.root","UPDATE","Histograms for STASRT digits");
  50     printf("Writting histograms to root file \n");
  51     Hfile->cd();
  52 //Create a canvas, set the view range, show histograms
  53     gStyle->SetOptStat(111111);
  54   //  TCanvas *c1 = new TCanvas("c1","Alice START Time ",400,10,600,600);
  55   hTimePs->SetXTitle("arriving time, ps");
  56   hTimePs->SetYTitle("number of events");
  57   hTimePs->Write();
  58   Hfile->Close();
  59
  60
  61 } // end of macro
  62
  63
  64
  65