X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=EMCAL%2FAliEMCALv1.cxx;h=52c48c1dd60a5c9ff93012286a704b3460c2bffb;hb=ccc5bf38dcae4c7351ba2b97390f10e79cfccef0;hp=02aa79562fcdd64a197810ba831c589860149955;hpb=e939a978a98d29901b86bdb02cdcaf5157e717de;p=u%2Fmrichter%2FAliRoot.git diff --git a/EMCAL/AliEMCALv1.cxx b/EMCAL/AliEMCALv1.cxx index 02aa79562fc..52c48c1dd60 100644 --- a/EMCAL/AliEMCALv1.cxx +++ b/EMCAL/AliEMCALv1.cxx @@ -30,6 +30,7 @@ // 2. Timing signal is collected and added to hit // --- ROOT system --- +#include #include #include @@ -41,6 +42,7 @@ #include "AliEMCALGeometry.h" #include "AliRun.h" #include "AliMC.h" +#include "AliStack.h" ClassImp(AliEMCALv1) @@ -132,6 +134,9 @@ void AliEMCALv1::StepManager(void){ AliEMCALGeometry * geom = GetGeometry() ; + TParticle *part; + Int_t parent; + static Int_t idXPHI = gMC->VolId("XPHI"); if(gMC->CurrentVolID(copy) == idXPHI ) { // We are in a Scintillator Layer Float_t depositedEnergy ; @@ -142,16 +147,20 @@ void AliEMCALv1::StepManager(void){ if (fCurParent==-1 || tracknumber != fCurTrack) { // Check parentage - Int_t parent=tracknumber; + //Int_t parent=tracknumber; + parent=tracknumber; if (fCurParent != -1) { while (parent != fCurParent && parent != -1) { - TParticle *part=gAlice->GetMCApp()->Particle(parent); + //TParticle *part=gAlice->GetMCApp()->Particle(parent); + part=gAlice->GetMCApp()->Particle(parent); parent=part->GetFirstMother(); } } if (fCurParent==-1 || parent==-1) { - Int_t parent=tracknumber; - TParticle *part=gAlice->GetMCApp()->Particle(parent); + //Int_t parent=tracknumber; + //TParticle *part=gAlice->GetMCApp()->Particle(parent); + parent=tracknumber; + part=gAlice->GetMCApp()->Particle(parent); while (parent != -1 && geom->IsInEMCAL(part->Vx(),part->Vy(),part->Vz())) { parent=part->GetFirstMother(); if (parent!=-1) @@ -161,7 +170,8 @@ void AliEMCALv1::StepManager(void){ if (fCurParent==-1) Error("StepManager","Cannot find parent"); else { - TParticle *part=gAlice->GetMCApp()->Particle(fCurParent); + //TParticle *part=gAlice->GetMCApp()->Particle(fCurParent); + part=gAlice->GetMCApp()->Particle(fCurParent); ienergy = part->Energy(); } while (parent != -1) { @@ -224,6 +234,7 @@ void AliEMCALv1::StepManager(void){ } } +//___________________________________________________________ void AliEMCALv1::RemapTrackHitIDs(Int_t *map) { // remap track index numbers for primary and parent indices // (Called by AliStack::PurifyKine) @@ -242,6 +253,7 @@ void AliEMCALv1::RemapTrackHitIDs(Int_t *map) { } } +//___________________________________________________________ void AliEMCALv1::FinishPrimary() { // finish primary fCurPrimary=-1;