//_________________________________________________________________________
// Hits class for PHOS
// A hit in PHOS is the sum of all hits in a single crystal
-//
-//*-- Author: Maxime Volkov (RRC KI) & Yves Schutz (SUBATECH)
+//*--
+//*-- Author: Maxime Volkov (RRC KI) & Yves Schutz (SUBATECH) & Dmitri Peressounko (RRC KI & SUBATECH)
// --- ROOT system ---
// --- Standard library ---
-#include <stdio.h>
#include <string.h>
#include <stdlib.h>
-#include <strstream.h>
// --- AliRoot header files ---
#include "AliPHOSHit.h"
#include "AliRun.h"
#include "AliConst.h"
-
+#include "AliPHOSLoader.h"
+#include "AliPHOSGeometry.h"
+#include "AliPHOS.h"
ClassImp(AliPHOSHit)
-
-//____________________________________________________________________________
-AliPHOSHit::AliPHOSHit(const AliPHOSHit & hit)
+
+ //____________________________________________________________________________
+ AliPHOSHit::AliPHOSHit(const AliPHOSHit & hit)
{
- // copy ctor
-
+ // copy ctor
+ fX = hit.fX ;
+ fY = hit.fY ;
+ fZ = hit.fZ ;
fId = hit.fId ;
fELOS = hit.fELOS ;
fPrimary = hit.fPrimary ;
fTrack = hit.fTrack ;
- fX = hit.fX ;
- fY = hit.fY ;
- fZ = hit.fZ ;
- fPid = hit.fPid ;
-
+ fTime = hit.fTime ;
+
}
//____________________________________________________________________________
-AliPHOSHit::AliPHOSHit(Int_t Shunt, Int_t primary, Int_t Track, Int_t id, Float_t *hits, Int_t pid) : AliHit(Shunt, Track)
+AliPHOSHit::AliPHOSHit(Int_t shunt, Int_t primary, Int_t track, Int_t id, Float_t *hits): AliHit(shunt, track)
{
- // ctor
+ //
+ // Create a CPV hit object
+ //
- fId = id ;
- fTrack = Track;
- fX = hits[0] ;
- fY = hits[1] ;
- fZ = hits[2] ;
- fELOS = hits[3] ;
- fPrimary = primary ;
- fPid = pid ;
+ fX = hits[0] ;
+ fY = hits[1] ;
+ fZ = hits[2] ;
+ fTime = hits[3] ;
+ fId = id ;
+ fELOS = hits[4] ;
+ fPrimary = primary ;
+}
+//____________________________________________________________________________
+Float_t AliPHOSHit::X() const
+{
+ // if(fX < -1000.){
+ TVector3 pos ;
+ AliPHOS * phos = static_cast<AliPHOS*> (gAlice->GetDetector("PHOS")) ;
+ phos->GetGeometry() ->RelPosInAlice(GetId(), pos) ;
+ return pos.X() ;
+ // fX = pos.X() ;
+ // fY = pos.Y() ;
+ // fZ = pos.Z() ;
+ // }
+ // return fX;
+}
+//____________________________________________________________________________
+Float_t AliPHOSHit::Y() const
+{
+ // if(fY < -1000.){
+ TVector3 pos ;
+ AliPHOS * phos = static_cast<AliPHOS*> (gAlice->GetDetector("PHOS")) ;
+ phos->GetGeometry() ->RelPosInAlice(GetId(), pos) ;
+ return pos.Y();
+ // fX = pos.X() ;
+ // fY = pos.Y() ;
+ // fZ = pos.Z() ;
+ // }
+ // return fY;
+}
+//____________________________________________________________________________
+Float_t AliPHOSHit::Z() const
+{
+ // if(fY < -1000.){
+ TVector3 pos ;
+ AliPHOS * phos = static_cast<AliPHOS*> (gAlice->GetDetector("PHOS")) ;
+ phos->GetGeometry() ->RelPosInAlice(GetId(), pos) ;
+ return pos.Z() ;
+ // fX = pos.X() ;
+ // fY = pos.Y() ;
+ // fZ = pos.Z() ;
+ // }
+ // return fZ;
}
-
//____________________________________________________________________________
Bool_t AliPHOSHit::operator==(AliPHOSHit const &rValue) const
{
- // Two hits are identical if they have the same Id and originate from the same primary
+ // Two hits are identical if they have the same Id and originate from the same primary
Bool_t rv = kFALSE ;
- if ( fId == rValue.GetId() && fPrimary == rValue.GetPrimary() )
+ if ( (fId == rValue.GetId()) && ( fPrimary == rValue.GetPrimary() ) )
rv = kTRUE;
return rv;
}
//____________________________________________________________________________
-AliPHOSHit AliPHOSHit::operator+(const AliPHOSHit &rValue) const
+AliPHOSHit AliPHOSHit::operator+(const AliPHOSHit &rValue)
{
// Add the energy of the hit
- AliPHOSHit added(*this);
-
- // the accumulated hit position is the position of the first hi
- // added.fX = rValue.fX ;
- // added.fY = rValue.fY ;
- // added.fZ = rValue.fZ ;
+ fELOS += rValue.GetEnergy() ;
- added.fELOS += rValue.GetEnergy() ;
+ if(rValue.GetTime() < fTime)
+ fTime = rValue.GetTime() ;
- return added;
+ return *this;
}
{
// Print out Id and energy
- out << "AliPHOSHit = " << hit.GetId() << " " << hit.GetEnergy() << endl ;
+ //out << "AliPHOSHit = " << hit.GetId() << " " << hit.GetEnergy() << " " << hit.GetTime() << endl ;
+ Warning("operator <<", "Implement differently") ;
return out ;
}