]>
Commit | Line | Data |
---|---|---|
19f796ed | 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 | } |