Do not include from subdirectories

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

diff --git a/HMPID/HESDfromKin.C b/HMPID/HESDfromKin.C
index b26a1efd78ae11e146edda34474bb1d13ee66663..042f4ca2c77ce4fb6e4bd56724c21882395a711a 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,11 +89,12 @@ 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
+      trk.SetLabel(i);
       pEsd->AddTrack(&trk);
     
       Int_t status;
-      if(!htaCheck) status = pTracker.Recon         (pEsd,pH->CluLst(),pNmean,pQthre);
+      if(!htaCheck) status = pTracker.ReconFromKin  (pEsd,pH->CluLst(),pNmean,pQthre);
       else          status = pTracker.ReconHiddenTrk(pEsd,pH->CluLst(),pNmean,pQthre);
 
 //      Printf("status %i",status);
@@ -114,7 +121,9 @@ 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->SetSpecificStorage("GRP/GRP/Data", Form("local://%s",gSystem->pwd()));
+  pCDB->SetRun(0);
   AliCDBEntry *pQthreEnt=pCDB->Get("HMPID/Calib/Qthre",0);
   AliCDBEntry *pNmeanEnt=pCDB->Get("HMPID/Calib/Nmean",0);
   
@@ -126,4 +135,3 @@ Bool_t OpenCalib()
   if(!pQthre || !pNmean) return kFALSE;  
   return kTRUE;
 }
-//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++