Index: trunk/MagicSoft/Mars/Changelog
===================================================================
--- trunk/MagicSoft/Mars/Changelog	(revision 2576)
+++ trunk/MagicSoft/Mars/Changelog	(revision 2577)
@@ -5,5 +5,10 @@
                                                  -*-*- END OF LINE -*-*-
 
-
+  2003/12/1: Abelardo Moralejo
+
+  * mimage/MStereoPar.cc
+    - Added calculation of fCT1Impact2, fCT2Impact2, the estimates of
+      impact parameter using the a priori assumption that the events
+      come from a point source located at the center of the camera.
 
   2003/11/27: Wolfgang Wittek
Index: trunk/MagicSoft/Mars/mimage/MStereoPar.cc
===================================================================
--- trunk/MagicSoft/Mars/mimage/MStereoPar.cc	(revision 2576)
+++ trunk/MagicSoft/Mars/mimage/MStereoPar.cc	(revision 2577)
@@ -234,4 +234,25 @@
 		       scalar * scalar );
 
+    //
+    // Now calculate i.p. assuming source is point-like and placed in
+    // the center of the camera.
+    //
+    scalar = (fCoreX2-ct1_x)*(-sin(mcevt1.GetTelescopeTheta())*
+			     cos(mcevt1.GetTelescopePhi()))  +
+      (fCoreY2-ct1_y)*(-sin(mcevt1.GetTelescopeTheta())*
+		      sin(mcevt1.GetTelescopePhi()));
+
+    fCT1Impact2 = sqrt( (fCoreX2-ct1_x)*(fCoreX2-ct1_x) +
+		       (fCoreY2-ct1_y)*(fCoreY2-ct1_y) -
+		       scalar * scalar );
+
+    scalar = (fCoreX2-ct2_x)*(-sin(mcevt2.GetTelescopeTheta())*
+			     cos(mcevt2.GetTelescopePhi()))  +
+      (fCoreY2-ct2_y)*(-sin(mcevt2.GetTelescopeTheta())*
+		      sin(mcevt2.GetTelescopePhi()));
+
+    fCT2Impact2 = sqrt( (fCoreX2-ct2_x)*(fCoreX2-ct2_x) +
+		       (fCoreY2-ct2_y)*(fCoreY2-ct2_y) -
+		       scalar * scalar );
 
  
@@ -287,14 +308,2 @@
 
 }
-
-
-
-
-
-
-
-
-
-
-
-