ACORDE/macros/ReadACORDEDigits.C

   1
   2 void ReadACORDEDigits()
   3
   4 {
   5
   6   // get run loader loader
   7    AliRunLoader* rl =
   8      AliRunLoader::Open("galice.root",
   9                         AliConfig::GetDefaultEventFolderName(),"read");
  10    if (rl == 0x0) {
  11      cerr<<"Can not open session for file galice.root\n";
  12      return;
  13    }
  14
  15    rl->LoadgAlice();
  16    gAlice = rl->GetAliRun();
  17
  18    // get acorde loader
  19    AliACORDE* acorde = (AliACORDE*)gAlice->GetDetector("ACORDE");
  20    AliLoader* aloader =rl->GetLoader("ACORDELoader");
  21
  22    // loop over events
  23    Int_t nevt = rl->GetNumberOfEvents();
  24    cout << " There are " << nevt << " evts " << endl;
  25    for (Int_t ievt=0; ievt<nevt;ievt++) {
  26      rl->GetEvent(ievt);
  27      aloader->LoadDigits("READ");
  28      TTree* treeD = aloader->TreeD();
  29      TClonesArray *adigits = new TClonesArray ("AliACORDEdigit", 1000);
  30      treeD->GetBranch("ACORDEdigit")->SetAddress(&adigits);
  31
  32      // loop over entries
  33      Int_t nent = treeD->GetEntries();
  34      cout << " There are " << nent << " entries in event " << ievt << endl;
  35
  36      for (Int_t ient=0;ient<nent;ient++) {
  37        acorde->ResetDigits();
  38        treeD->GetEvent(ient);
  39        Int_t ndig = adigits->GetEntriesFast();
  40        cout << " There are " << ndig << " digits in entry " << ient << endl;
  41        for (Int_t idig=0;idig<ndig;idig++) {
  42          AliACORDEdigit* digit = (AliACORDEdigit*) adigits->At(idig);
  43          Int_t mod = digit->GetModule();
  44          Float_t time = digit->GetTime();
  45          cout << " Digit " << idig
  46               <<" : module="<< mod << " time=" << time << endl;
  47          for (Int_t i=0;i<3;i++)
  48            cout << " track " << i << " is " << digit->GetTrack(i) << endl;
  49        } // end loop over digits
  50      } // end loop over entries
  51      aloader->UnloadDigits();
  52    } // end loop over events
  53 }