Changeset 9234 for trunk/MagicSoft/Mars/mgeom

Timestamp:

01/21/09 14:26:15 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mgeom

Files:

• MGeomCam.h (modified) (2 diffs)
• MGeomCamDwarf.cc (modified) (1 diff)
• MGeomCamDwarf.h (modified) (1 diff)
• MGeomCamMagic.cc (modified) (1 diff)
• MGeomCamMagic.h (modified) (1 diff)

trunk/MagicSoft/Mars/mgeom/MGeomCam.h

@@ -13 +13 @@
 #ifndef ROOT_TArrayS
 #include <TArrayS.h>
+#endif
+#ifndef MARS_MQuaternion
+#include "MQuaternion.h"
 #endif
 
@@ -109 +112 @@
     Int_t GetNeighbor(UInt_t idx, Int_t dir) const;
 
+    virtual Bool_t HitFrame(MQuaternion p, const MQuaternion &u) const { return kFALSE; }
+    virtual Bool_t HitDetector(const MQuaternion &p) const { return kFALSE; }
+
     void Print(Option_t *opt=NULL) const;
 

trunk/MagicSoft/Mars/mgeom/MGeomCamDwarf.cc

@@ -112 +112 @@
 // --------------------------------------------------------------------------
 //
+// Check if the photon which is flying along the trajectory u has passed
+// (or will pass) the frame of the camera (and consequently get
+// absorbed). The position p and direction u must be in the
+// telescope coordinate frame, which is z parallel to the focal plane,
+// x to the right and y upwards, looking from the mirror towards the camera.
+//
+// The units are cm.
+//
+Bool_t MGeomCamDwarf::HitFrame(MQuaternion p, const MQuaternion &u) const
+{
+    // z is defined from the mirror (0) to the camera (z>0).
+    // Thus we just propagate to the focal plane (z=fDist)
+    //p -= 1700./u.Z()*u;
+    p.PropagateZ(u, GetCameraDist()*100); // m->cm
+
+    // Add 10% to the max radius and convert from mm to cm
+    return p.R()<GetMaxRadius()*0.11;//TMath::Abs(p.X())<65 && TMath::Abs(p.Y())<65;
+}
+
+// --------------------------------------------------------------------------
+//
+// Check if the position given in the focal plane (so z can be ignored)
+// is a position which might hit the detector. It is meant to be a rough
+// and fast estimate not a precise calculation. All positions dicarded
+// must not hit the detector. All positions accepted might still miss
+// the detector.
+//
+// The units are cm.
+//
+Bool_t MGeomCamDwarf::HitDetector(const MQuaternion &v) const
+{
+    // Add 10% to the max radius and convert from mm to cm
+    return v.R()<GetMaxRadius()*0.11;//60.2;
+}
+
+// --------------------------------------------------------------------------
+//
 // Calculate in the direction 0-5 (kind of sector) in the ring-th ring
 // the x and y coordinate of the i-th pixel. The unitx are unity,

trunk/MagicSoft/Mars/mgeom/MGeomCamDwarf.h

@@ -28 +28 @@
     TObject *Clone(const char *newname) const;
 
-    ClassDef(MGeomCamDwarf, 1)          // Geometry class for the Dwarf camera
+    Bool_t HitFrame(MQuaternion p, const MQuaternion &u) const;
+    Bool_t HitDetector(const MQuaternion &p) const;
+
+    ClassDef(MGeomCamDwarf, 1) // Geometry class for the Dwarf camera
 };
 
 #endif
-

trunk/MagicSoft/Mars/mgeom/MGeomCamMagic.cc

@@ -49 +49 @@
     CreateNN();
     InitGeometry();
+}
+
+// --------------------------------------------------------------------------
+//
+// Check if the photon which is flying along the trajectory u has passed
+// (or will pass) the frame of the camera (and consequently get
+// absorbed). The position p and direction u must be in the
+// telescope coordinate frame, which is z parallel to the focal plane,
+// x to the right and y upwards, looking from the mirror towards the camera.
+//
+Bool_t MGeomCamMagic::HitFrame(MQuaternion p, const MQuaternion &u) const
+{
+    // z is defined from the mirror (0) to the camera (z>0).
+    // Thus we just propagate to the focal plane (z=fDist)
+    //p -= 1700./u.Z()*u;
+    p.PropagateZ(u, GetCameraDist()*100);
+
+    return TMath::Abs(p.X())<65 && TMath::Abs(p.Y())<65;
+}
+
+// --------------------------------------------------------------------------
+//
+// Check if the position given in the focal plane (so z can be ignored)
+// is a position which might hit the detector. It is meant to be a rough
+// and fast estimate not a precise calculation. All positions dicarded
+// must not hit the detector. All positions accepted might still miss
+// the detector.
+//
+Bool_t MGeomCamMagic::HitDetector(const MQuaternion &v) const
+{
+    return v.R()<60.2;
 }
 

trunk/MagicSoft/Mars/mgeom/MGeomCamMagic.h

@@ -15 +15 @@
     MGeomCamMagic(const char *name=NULL);
 
+    Bool_t HitFrame(MQuaternion p, const MQuaternion &u) const;
+    Bool_t HitDetector(const MQuaternion &p) const;
+
     ClassDef(MGeomCamMagic, 1)          // Geometry class for the Magic camera
 };