added method to say if hit is in EMC or not

diff --git a/PHOS/AliPHOSGeometry.cxx b/PHOS/AliPHOSGeometry.cxx
index 718901736f10697d3e310b4a1adf903d42d35747..76f907bd7de0a7dfec3050884b4a6b47823fdc7f 100644 (file)
--- a/PHOS/AliPHOSGeometry.cxx
+++ b/PHOS/AliPHOSGeometry.cxx
@@ -121,6 +121,17 @@ void AliPHOSGeometry::Init(void)
   }
 }
 
+
+//____________________________________________________________________________
+Bool_t AliPHOSGeometry::IsInEMC(const Int_t id) const { 
+    Int_t relid[4] ; 
+    AbsToRelNumbering(id, relid) ; 
+    if ( relid[1] == 0 ) 
+      return kTRUE ; 
+    else 
+      return kFALSE ; 
+  }
+
 //____________________________________________________________________________
 Float_t AliPHOSGeometry::GetCPVBoxSize(Int_t index)  const
 {

diff --git a/PHOS/AliPHOSGeometry.h b/PHOS/AliPHOSGeometry.h
index 862d902b8826682512ad1ec3d7fc03ad9663d0c9..e92195f22bdbf228221c0ae45d006314500d90e0 100644 (file)
--- a/PHOS/AliPHOSGeometry.h
+++ b/PHOS/AliPHOSGeometry.h
@@ -84,7 +84,7 @@ public:
   void ImpactOnEmc(const Double_t theta, const Double_t phi, Int_t & ModuleNumber, Double_t & z, Double_t & x) const ; 
                                                          // calculates the impact coordinates of a neutral particle  
                                                          // emitted in direction theta and phi in ALICE
- 
+  Bool_t IsInEMC(const Int_t id) const ;
   void   RelPosInModule(const Int_t * RelId, Float_t & y, Float_t & z) const ; // gets the position of element (pad or Xtal) relative to 
                                                                          // center of PHOS module  
   void   RelPosInAlice(const Int_t AbsId, TVector3 &  pos) const ;             // gets the position of element (pad or Xtal) relative to