Changeset 2577

Timestamp:

12/01/03 13:59:37 (21 years ago)

Author:

moralejo

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mimage/MStereoPar.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -5 +5 @@
                                                  -*-*- 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

trunk/MagicSoft/Mars/mimage/MStereoPar.cc

@@ -234 +234 @@
                        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 +308 @@
 
 }
-
-
-
-
-
-
-
-
-
-
-
-