-void AliEMCALv3::FinishEvent()
-{ // 26-may-05
- static double de=0.;
- fHNhits->Fill(double(fHits->GetEntries()));
- de = GetDepositEnergy(0);
- if(fHDe) fHDe->Fill(de);
-}
-
-Double_t AliEMCALv3::GetDepositEnergy(int print)
-{ // 23-mar-05 - for testing
- // cout<<"AliEMCALv3::GetDepositEnergy() : fHits "<<fHits<<endl;
- if(fHits == 0) return 0.;
- static AliEMCALHitv1 *hit=0;
- static Double_t de=0., deHit=0., zhl=0.0, zShift = 0.16*77;
- for(int ih=0; ih<fHits->GetEntries(); ih++) {
- hit = (AliEMCALHitv1*)fHits->UncheckedAt(ih);
- deHit = (double)hit->GetEnergy();
- zhl = (double)hit->Z() + zShift;
- if(zhl>2.*zShift) zhl = 2*zShift - 0.002;
- de += deHit;
- if(fHDeDz) fHDeDz->Fill(zhl, deHit);
- }
- if(print>0) {
- printf(" #hits %i de %f \n", fHits->GetEntries(), de);
- if(print>1) {
- printf(" #primary particles %i\n", gAlice->GetHeader()->GetNprimary());
- }
- }
- return de;
-}
-