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();
pHL->TreeR()->GetEntry(0);
AliStack *pStack=gAL->Stack();
Int_t nTrkHMPID=0;
+
for(Int_t i=0;i<pStack->GetNtrack();i++){
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
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);
- // if(!(iEvt%50)) Printf("Number of events processed: %i with tracks %i in HMPID",iEvt,nTrkHMPID);
- Printf("Number of events processed: %i with tracks %i in HMPID",iEvt,nTrkHMPID);
-
Int_t status;
if(!htaCheck) status = pTracker.Recon (pEsd,pH->CluLst(),pNmean,pQthre);
else status = pTracker.ReconHiddenTrk(pEsd,pH->CluLst(),pNmean,pQthre);
- Printf("status %i",status);
+// Printf("status %i",status);
if(status !=0) continue;
}// track loop
+
+ if(!(iEvt%50)) Printf("Number of events processed: %i with tracks %i in HMPID",iEvt,nTrkHMPID);
+// Printf("Number of events processed: %i with tracks %i in HMPID",iEvt,nTrkHMPID);
+
gEsdTr->Fill();
pEsd->Reset();
}// event loop
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);
if(!pQthre || !pNmean) return kFALSE;
return kTRUE;
}
-//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++