updating to new Shuttle processing logic

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

diff --git a/HMPID/HESDfromKin.C b/HMPID/HESDfromKin.C
index b26a1efd78ae11e146edda34474bb1d13ee66663..f0694adff2a6daadddb9c413991de50ec647a715 100644 (file)
--- a/HMPID/HESDfromKin.C
+++ b/HMPID/HESDfromKin.C
@@ -54,7 +54,12 @@ void SimEsd(AliLoader *pHL,AliESDEvent *pEsd,Bool_t htaCheck)
     Printf("| SimESD: Utility to embed ESD from kinematics|");
     Printf("-----------------------------------------------");
 }
-  AliHMPIDTracker::SetFieldMap(gAL->GetAliRun()->Field(),kTRUE);
+
+  AliGRPManager *grpMan = new AliGRPManager();
+  
+  grpMan->ReadGRPEntry();
+  grpMan->SetMagField();
+
   AliHMPIDTracker pTracker;
   AliHMPID *pH=(AliHMPID*)gAL->GetAliRun()->GetDetector("HMPID");
   Int_t iNevt=gAL->GetNumberOfEvents();
@@ -69,6 +74,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 
@@ -83,7 +89,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;
@@ -114,7 +120,8 @@ void SimEsd(AliLoader *pHL,AliESDEvent *pEsd,Bool_t htaCheck)
 Bool_t OpenCalib()
 {
   AliCDBManager* pCDB = AliCDBManager::Instance();
-  pCDB->SetDefaultStorage("local://$HOME");
+  pCDB->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
+  pCDB->SetRun(0);
   AliCDBEntry *pQthreEnt=pCDB->Get("HMPID/Calib/Qthre",0);
   AliCDBEntry *pNmeanEnt=pCDB->Get("HMPID/Calib/Nmean",0);
   
@@ -126,4 +133,3 @@ Bool_t OpenCalib()
   if(!pQthre || !pNmean) return kFALSE;  
   return kTRUE;
 }
-//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++