Use access methods to hit data.

diff --git a/MUON/MUONacc.C b/MUON/MUONacc.C
index 66f4efeba633c6abad6aa846b5068786a05a9af6..439f6659feab184784bf5b1fdf1875bb500d5b78 100644 (file)
--- a/MUON/MUONacc.C
+++ b/MUON/MUONacc.C
@@ -142,10 +142,7 @@ void MUONacc (Float_t d=3., Int_t evNumber1=0, Int_t evNumber2=0)
                        
                        if (Particle != kMuonPlus && Particle != kMuonMinus) continue;
                        
-                       Float_t P    =
-                           TMath::Sqrt(mHit->Cx()*mHit->Cx()+
-                                       mHit->Cy()*mHit->Cy()+
-                                       mHit->Cz()*mHit->Cz());
+                       Float_t P   = mhit->Momentum();
                        Float_t R   = TMath::Sqrt(x*x+y*y);
                        TParticlePDG* Part = DataBase->GetParticle(Particle);
                        Double_t mass = Part->Mass();

diff --git a/MUON/MUONhits.C b/MUON/MUONhits.C
index cae5ce67537ba35a29c8a44e545e0bbf72f990f6..55bb20c9c5c225de5202444715b5d76f89579bda 100644 (file)
--- a/MUON/MUONhits.C
+++ b/MUON/MUONhits.C
@@ -135,10 +135,7 @@ void MUONhits (Int_t evNumber1=0,Int_t evNumber2=0, Int_t ic=1)
                    Float_t Eloss =    mHit->Eloss();
                    Float_t Theta =    mHit->Theta();
                    Float_t Particle = mHit->Particle();
-                   Float_t P    =     
-                       TMath::Sqrt(mHit->Cx()*mHit->Cx()+
-                                   mHit->Cy()*mHit->Cy()+
-                                   mHit->Cz()*mHit->Cz());
+                   Float_t P     =    mHit->Momentum(); 
                    TParticlePDG* Part = DataBase->GetParticle(Particle);
                    Double_t mass = Part->Mass();
                    

diff --git a/MUON/MUONstraggling.C b/MUON/MUONstraggling.C
index 1674453f456a5ecca3b210953319de6a97d18b26..b6ee680f0c816591520d4aa090878c0edb62d89d 100644 (file)
--- a/MUON/MUONstraggling.C
+++ b/MUON/MUONstraggling.C
@@ -78,22 +78,21 @@ void MUONstraggling (Int_t evNumber1=0,Int_t evNumber2=0)
                   mHit;
                   mHit=(AliMUONHit*)MUON->NextHit()) 
               {
-                  Int_t   nch   = mHit->fChamber;  // chamber number
-                  Float_t x     = mHit->fX;        // x-pos of hit
-                  Float_t y     = mHit->fY;        // y-pos
-                  Float_t z     = mHit->fZ;        // y-pos
-                  Float_t p=mHit->fPTot;
-                  Float_t px=mHit->fCxHit;
-                  Float_t py=mHit->fCyHit;
-                  Float_t pz=mHit->fCzHit;
+                  Int_t   nch   = mHit->Chamber();  // chamber number
+                  Float_t x     = mHit->X();        // x-pos of hit
+                  Float_t y     = mHit->Y();        // y-pos
+                  Float_t z     = mHit->Z();        // y-pos
+                  Float_t p=mHit->Momentum();
+                  Float_t px=mHit->Px();
+                  Float_t py=mHit->Py();
+                  Float_t pz=mHit->Pz();
                   
                   if (nch != 1) continue;
                   
-                  Int_t ipart = mHit->fParticle;
-                  TClonesArray *fPartArray = gAlice->Particles();
+                  Int_t ipart = mHit->Particle();
                   TParticle *Part;
-                  Int_t ftrack = mHit->fTrack;
-                  Part = (TParticle*) fPartArray->UncheckedAt(ftrack);
+                  Int_t ftrack = mHit->Track();
+                  Part = gAlice->Particle(ftrack);
                   Int_t ipart = Part->GetPdgCode();
                   TParticle *Mother;
                   Float_t px0=Part->Px();

diff --git a/MUON/MUONtestabso.C b/MUON/MUONtestabso.C
index 8a8feaffcde44baa75a75bbf05d20ae1f125853d..8ea69da93f00990401a790e60563faaf60bf35be 100644 (file)
--- a/MUON/MUONtestabso.C
+++ b/MUON/MUONtestabso.C
@@ -137,8 +137,8 @@ void MUONtestabso (Int_t evNumber1=0,Int_t evNumber2=0)
                       if (dE<5.5) hde2->Fill(dE, (float) 1);
                       xmuon=mHit->fX;
                       ymuon=mHit->fY;
-                      Float_t thetax=TMath::ATan2(mHit->fCxHit, mHit->fPTot);
-                      Float_t thetay=TMath::ATan2(mHit->fCyHit, mHit->fPTot);
+                      Float_t thetax=TMath::ATan2(mHit->Px(), mHit->Momentum());
+                      Float_t thetay=TMath::ATan2(mHit->Py(), mHit->Momentum());
                       dtheta->Fill((thetax-thetax0)*1000., 1.);
                       dtheta->Fill((thetay-thetay0)*1000., 1.);                       
                   }