--- PamCut/CaloCuts/CaloGeomCut/CaloGeomCut.cpp	2011/07/28 13:26:31	1.5
+++ PamCut/CaloCuts/CaloGeomCut/CaloGeomCut.cpp	2011/09/02 16:25:38	1.6
@@ -14,9 +14,9 @@
 const Float_t CaloGeomCut::_zCaloTop = -26.181;
 
 CaloGeomCut::CaloGeomCut(const char *cutName, CaloAxis *xCaloAxis, CaloAxis *yCaloAxis, UInt_t iMethod,
-    Float_t xTolCaloTrack, Float_t yTolCaloTrack, Float_t mass, Float_t resMax, Float_t qualCut, Float_t chi2Cut) :
+    Float_t xTolCaloTrack, Float_t yTolCaloTrack, Float_t mass, Float_t zeta, Float_t resMax, Float_t qualCut, Float_t chi2Cut) :
     PamCut(cutName), _xCaloAxis(xCaloAxis), _yCaloAxis(yCaloAxis), _iMethod(iMethod), _xTolCaloTrack(xTolCaloTrack), _yTolCaloTrack(
-        yTolCaloTrack), _mass(mass), _resMax(resMax), _qualCut(qualCut), _chi2Cut(chi2Cut), _traj(NULL), _zIn(_nPoint) {
+        yTolCaloTrack), _mass(mass),  _zeta(zeta), _resMax(resMax), _qualCut(qualCut), _chi2Cut(chi2Cut), _traj(NULL), _zIn(_nPoint) {
 
   if (_iMethod == 1) {
     _traj = (Trajectory*) new Trajectory(_nPoint);
@@ -56,7 +56,7 @@
     alCaloTop[3] = TMath::Pi() + atan2(tany, tanx); // [phi]
 
     beta = event->GetToFLevel2()->CalcBeta(0, _resMax, _qualCut, _chi2Cut); // TOF stand-alone beta
-    alCaloTop[4] = 1. / (_mass * beta / sqrt(1. - pow(beta, 2))); // [etaP] [rho = (m/Z)*(beta/(sqrt(1-beta^2)))]
+    alCaloTop[4] = 1. / ((_mass/_zeta) * beta / sqrt(1. - pow(beta, 2))); // [etaP] [rho = (m/Z)*(beta/(sqrt(1-beta^2)))]
 
     _traj->DoTrack(alCaloTop, _zCaloTop);