PMD/AliPMDHits2SDigits.C

   1 // ----------------------------------------------------//
   2 //                                                     //
   3 //       This macro does Hits to SDigits               //
   4 //                                                     //
   5 // ----------------------------------------------------//
   6
   7 #include <stdlib.h>
   8 #include "Riostream.h"
   9 #include "TROOT.h"
  10 #include "TFile.h"
  11 #include "TNetFile.h"
  12 #include "TRandom.h"
  13 #include "TTree.h"
  14 #include "TBranch.h"
  15 #include "TClonesArray.h"
  16 #include "TStopwatch.h"
  17
  18 void AliPMDHits2SDigits(Int_t nevt=1)
  19 {
  20   TStopwatch timer;
  21   timer.Start();
  22   Float_t zpos = 361.5;
  23
  24   // Input (and output) file name
  25   Char_t *alifile = "galice.root";
  26
  27   // Create the PMD digitzer
  28   AliPMDDigitizer *digitizer = new AliPMDDigitizer();
  29
  30   // Open the AliRoot file
  31   digitizer->OpengAliceFile(alifile,"HS");
  32
  33   digitizer->SetZPosition(zpos);
  34
  35   // Create the sdigits
  36
  37   for (Int_t ievt=0; ievt<nevt; ievt++)
  38     {
  39       digitizer->Hits2SDigits(ievt);
  40     }
  41
  42   digitizer->UnLoad("S");
  43   timer.Stop();
  44   timer.Print();
  45
  46 }
  47