HLT/PHOS/ana/energySpectrum/energySpectrum.C

   1
   2 Int_t energySpectrum(const char* runNb)
   3 {
   4   gROOT->ProcessLine(".L ../digits/AliHLTPHOSAltroConfig.cxx++");
   5   gROOT->ProcessLine(".L ../digits/AliHLTPHOSDigit.cxx++");
   6   TH1F *spectrumHist = new TH1F("spectrumHist", "Digit energy spectrum", 200, 0, 200);
   7   TH1F *spectrumHistSingle = new TH1F("spectrumHist", "Digit energy spectrum", 200, 0, 200);
   8   TClonesArray *digArray = new TClonesArray("AliHLTPHOSDebugRawDigit" , 100);
   9   TChain *tree= new TChain("digitTree");
  10   AliHLTPHOSDebugRawDigit *digit = 0;
  11
  12   char filepath [50];
  13   sprintf(filepath, "/tmp/phoshlt/analysis/data/run%s/*", runNb);
  14   tree->Add(filepath);
  15
  16   tree->SetBranchAddress("DebugRawDigit", &digArray);
  17   cout << endl << "Entries in tree: " << tree->GetEntries() << endl;
  18
  19   for(int k = 0; k < tree->GetEntries(); k++)
  20     {
  21       tree->GetEntry(k);
  22       for(int j = 0; j < digArray->GetEntriesFast(); j++)
  23         {
  24           digit = (AliHLTPHOSDebugRawDigit*)digArray->At(j);
  25           if(digit->GetAmplitude() > 20&& digit->GetCrazyness()==1 && digit->GetGain() == 1)
  26             {
  27             spectrumHist->Fill(digit->GetAmplitude());
  28           if(digit->fX = 30 &&digit->fZ = 30)
  29             {
  30               spectrumHistSingle->Fill(digit->GetAmplitude());
  31             }
  32             }
  33         }
  34     }
  35   spectrumHist->Draw();
  36
  37   return 0;
  38 }
  39
  40