START/AliSTARTHist.cxx

   1 #include "AliSTARTRawReader.h"
   2 #include "AliSTARTRawData.h"
   3 #include "/home/alla/AliRoot/verynew/RAW/AliRawReaderFile.h"
   4
   5
   6 #include <Riostream.h>
   7 #include "TMath.h"
   8 #include "TH1F.h"
   9 #include "TH2F.h"
  10 #include "TArrayI.h"
  11 #include "AliSTARTHist.h"
  12 #include "AliRunLoader.h"
  13 #include "AliLoader.h"
  14 #include "/home/alla/AliRoot/verynew/RAW/AliRawReader.h"
  15 #include "/home/alla/AliRoot/verynew/RAW/AliRawReaderFile.h"
  16 #include "AliSTARTRawReader.h"
  17 ClassImp(AliSTARTHist)
  18
  19   AliSTARTHist::AliSTARTHist() : TObject()
  20 {
  21   fhADCR = new TH1F("hADCR","ADC right",100,1.,500);
  22   fhADCL = new TH1F("hADCL","ADC right",100,1.,500);
  23   fhTimeR = new TH1F("hTimeR","Time Right",100,400,1000);
  24   fhTimeL = new TH1F("hTimeL","Time Left",100,10,500);
  25   fhBestTimeR = new TH1F("hBestTimeR","Time Right",100,400,1000);
  26   fhBestTimeL = new TH1F("hBestTimeL","Time Left",100,0,500);
  27   fhADCDet = new TH1F("hADCDet","ADC vs Ndet",30,0.,30);
  28   fhTimeDet = new TH1F("hTimeDet","Time vs Ndet",30,0,30);
  29   fhTimeDiff = new TH1F("hTimeDiff"," Time diff",100,350,450);
  30   fhMeanTime = new TH1F("hMeanTime"," Mean Time ",100,250,350);
  31   fhT0detR = new TH1F("hT0detR"," T0 vs Ndet right ",15,0,15);
  32   fhT0detL = new TH1F("hT0detL"," T0 vs Ndet left ",15,0,15);
  33   fhTimevsADCL = new TH2F("hTimevsADCL","time vs ADC left",100,0,500,100,1,500);
  34   fhTimevsADCR = new TH2F("hTimevsADCR","time vs ADC right",100,0,500,100,400,1000);
  35 }
  36
  37
  38   AliSTARTHist::~AliSTARTHist()
  39 {
  40   /*
  41   delete fhADCR;
  42   delete fhADCL;
  43   delete fhTimeR;
  44   delete fhTimeL;
  45   delete fhBestTimeR;
  46   delete fhBestTimeL;
  47   delete fhADCDet;
  48   delete fhTimeDet;
  49   delete fhTimeDiff;
  50   delete fhMeanTime;
  51   delete fhT0detR;
  52   delete fhT0detL;
  53   delete fhTimevsADCL;
  54   */
  55 }
  56
  57 void AliSTARTHist::FillHist(AliRunLoader* runLoader,
  58                             AliRawReader* rawReader) const
  59 {
  60
  61    TArrayI *fTime= new TArrayI(24);
  62    TArrayI *fADC = new TArrayI(24);
  63    AliSTARTRawReader myrawreader;
  64
  65
  66
  67   Int_t iEvent = 0;
  68   while (rawReader->NextEvent()) {
  69       runLoader->GetEvent(iEvent++);
  70
  71       myrawreader.NextThing( rawReader);
  72       myrawreader.GetTime(*fTime);
  73       myrawreader.GetADC(*fADC);
  74
  75       Int_t besttimeleft=99999;
  76       Int_t besttimeright=99999;
  77       Int_t ibestR=999; Int_t ibestL=999;
  78       for (Int_t i=0; i<12; i++ )
  79         {
  80           Int_t timel=fTime->At(i);
  81           if(timel<besttimeleft && timel!=0) { besttimeleft=timel; ibestL=i;}
  82           fhTimeL->Fill(timel);
  83           fhTimeDet->Fill(i,timel);
  84           Int_t adcl=fADC->At(i);
  85           fhADCL->Fill(adcl);
  86           fhADCDet->Fill(i,adcl);
  87           Int_t timer=fTime->At(i+12);
  88           if(timer<besttimeright && timer!=0){ besttimeright=timer; ibestR=i;}
  89           fhTimeDet->Fill(i+12,timer);
  90           fhTimeR->Fill(timer);
  91           Int_t adcr=fADC->At(i+12);
  92           fhADCR->Fill(adcr);
  93           fhADCDet->Fill(i+12,adcr);
  94           fhTimevsADCL->Fill(adcl,timel);
  95           fhTimevsADCR->Fill(adcr,timer);
  96
  97         }
  98
  99       fhBestTimeR->Fill(besttimeright);
 100       fhBestTimeL->Fill(besttimeleft);
 101       Int_t timeDiff1=besttimeright-besttimeleft;
 102       Float_t meanTime=(besttimeright+besttimeleft)/2.;
 103       fhTimeDiff->Fill(timeDiff1);
 104       fhMeanTime->Fill(meanTime);
 105       //     Float_t  timeDiff2=((besttimeright-besttimeleft)-(350.-69.7))/2;
 106       fhT0detR->Fill(ibestR); fhT0detL->Fill(ibestL);
 107   }
 108 }
 109
 110
 111