}
//____________________________________________________________________________
-TVector3 AliEMCALPIDv1::GetMomentumDirection(AliEMCALTowerRecPoint * emc)const
+TVector3 AliEMCALPIDv1::GetMomentumDirection(AliEMCALRecPoint * emc)const
{
// Calculates the momentum direction:
// direction is given by IP and this RecPoint
TVector3 dir(0,0,0) ;
TVector3 emcglobalpos ;
- TMatrix dummy ;
+ // TMatrix dummy ;
- emc->GetGlobalPosition(emcglobalpos, dummy) ;
+ emc->GetGlobalPosition(emcglobalpos) ;
dir = emcglobalpos ;
- dir.SetZ( -dir.Z() ) ; // why ?
// dir.SetMag(1.) ; Removed to avoid warings !!!!!!!!!!!!!! TO BE REVISED
//account correction to the position of IP
Print("") ;
return ;
}
- AliEMCALGetter * gime = AliEMCALGetter::Instance(GetTitle()) ;
-
+ AliEMCALGetter * gime = AliEMCALGetter::Instance() ;
+
if (fLastEvent == -1)
fLastEvent = gime->MaxEvent() - 1 ;
else
Int_t ievent ;
-
for (ievent = fFirstEvent; ievent <= fLastEvent; ievent++) {
gime->Event(ievent,"R") ;
MakeRecParticles() ;
if(strstr(option,"deb"))
PrintRecParticles(option) ;
//increment the total number of rec particles per run
- fRecParticlesInRun += gime->RecParticles()->GetEntriesFast() ;
-
+ fRecParticlesInRun += gime->RecParticles()->GetEntriesFast() ;
}
if(strstr(option,"tim")){
gBenchmark->Stop("EMCALPID");
recParticles->Clear();
TIter next(aECARecPoints) ;
- AliEMCALTowerRecPoint * eca ;
+ AliEMCALRecPoint * eca ;
Int_t index = 0 ;
AliEMCALRecParticle * rp ;
- while ( (eca = (AliEMCALTowerRecPoint *)next()) ) {
+ while ( (eca = (AliEMCALRecPoint *)next()) ) {
new( (*recParticles)[index] ) AliEMCALRecParticle() ;
rp = (AliEMCALRecParticle *)recParticles->At(index) ;
rp->SetFirstDaughter(-1);
rp->SetLastDaughter(-1);
rp->SetPolarisation(0,0,0);
+ //Set the position in global coordinate system from the RecPoint
+ //AliEMCALGeometry * geom = gime->EMCALGeometry() ;
+ //AliEMCALTowerRecPoint * erp = gime->ECARecPoint(rp->GetEMCALRPIndex()) ;
+ TVector3 pos ;
+ //geom->GetGlobal(erp, pos) ; !!!!!!!!!! to check
+ rp->SetPos(pos);
+
+
index++ ;
}