HMPID tracking classes updated according to the modifications in the central tracking...

[u/mrichter/AliRoot.git] / HMPID / HESDfromKin.C

diff --git a/HMPID/HESDfromKin.C b/HMPID/HESDfromKin.C
index effd7d04410d9995aa25fd49a128da3500098543..4e80cce8143dff82eec0e6532157f098750772c4 100644 (file)
--- a/HMPID/HESDfromKin.C
+++ b/HMPID/HESDfromKin.C
@@ -73,6 +73,7 @@ void SimEsd(AliLoader *pHL,AliESDEvent *pEsd,Bool_t htaCheck)
       if(!pStack->IsPhysicalPrimary(i)) continue;
       TParticle *pTrack=pStack->Particle(i); 
       if(pTrack->GetPDG()->Charge()==0) continue;
+      Printf("track n. %i",i);
       AliESDtrack trk(pTrack); 
       Float_t xPc,yPc,xRa,yRa,thRa,phRa;
       Int_t iCh=pTracker.IntTrkCha(&trk,xPc,yPc,xRa,yRa,thRa,phRa);         //get chamber intersected by this track 
@@ -87,7 +88,7 @@ void SimEsd(AliLoader *pHL,AliESDEvent *pEsd,Bool_t htaCheck)
       
       if(phRa<0) phRa += TMath::TwoPi(); // to be verified
       
-      trk.SetHMPIDtrk(xRa,yRa,thRa,phRa);                                                        //store initial infos
+      trk.SetHMPIDtrk(xPc,yPc,thRa,phRa);                                                        //store initial infos
       pEsd->AddTrack(&trk);
     
       Int_t status;
@@ -322,7 +323,9 @@ Bool_t InitGRP() {
   else {
     Printf(Form("Cannot determine the beam type from %s, assume no LHC magnet field",beamtype));
   }
-  
+  Printf("------------------------------");
+  Printf(" Summary for B: %s",s.Data());
+  Printf("------------------------------");
   AliMagF* fld = new AliMagF("MagneticFieldMap", s.Data(), 2, fcL3, fcDip, 10., map, path, 
                             btype,beamenergy);
   TGeoGlobalMagField::Instance()->SetField( fld );