Changeset 9565 for trunk/MagicSoft/Mars/msimreflector

Timestamp:

03/30/10 14:47:27 (15 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/msimreflector

Files:

• MMirror.cc (modified) (1 diff)
• MMirror.h (modified) (3 diffs)
• MReflector.cc (modified) (3 diffs)
• MSimReflector.cc (modified) (9 diffs)
• MSimReflector.h (modified) (1 diff)

trunk/MagicSoft/Mars/msimreflector/MMirror.cc

@@ -99 +99 @@
 using namespace std;
 
+void MMirror::SetShape(Char_t c)
+{
+    switch (toupper(c))
+    {
+    case 'S':
+        fShape = 0;
+        break;
+
+    case 'P':
+        fShape = 1;
+        break;
+
+    default:
+        fShape = 0;
+    }
+}
+
 // --------------------------------------------------------------------------
 //

trunk/MagicSoft/Mars/msimreflector/MMirror.h

@@ -25 +25 @@
     Double_t  fSigmaPSF;
 
+    Int_t fShape; // Spherical=0, Parabolic=1
+
     //    MMirror *fNeighbors[964];
 
 public:
-    MMirror() : fSigmaPSF(-1)
+    MMirror() : fSigmaPSF(-1), fShape(0)
     {
     }
@@ -48 +50 @@
         fTilt.Rotate(-n.Theta(), TVector3(-n.Y(), n.X(), 0));
     }
+    void SetShape(Char_t c);
 
     void SetZ(Double_t z) { fPos.SetZ(z); }
@@ -95 +98 @@
     }
     */
-    ClassDef(MMirror, 1) // Base class to describe a mirror
+    ClassDef(MMirror, 2) // Base class to describe a mirror
 };
 

trunk/MagicSoft/Mars/msimreflector/MReflector.cc

@@ -162 +162 @@
                         atof(arr[5]->GetName()));
 
-    const Double_t F = atof(arr[6]->GetName());
-
-    const TString val = arr[7]->GetName();
+    UInt_t n = 6;
+
+    TString six = arr[n]->GetName();
+
+    Char_t shape = 'S';
+    if (!six.IsFloat())
+    {
+        shape = six[0];
+        n++;
+    }
+
+    const Double_t F = atof(arr[n++]->GetName());
+
+    const TString val = arr[n++]->GetName();
 
     const Double_t psf = val.IsFloat() ? val.Atof() : -1;
 
-    const UInt_t n = val.IsFloat() ? 9 : 8;
+    n += val.IsFloat() ? 1 : 0;
 
     TString type = arr[n-1]->GetName();
@@ -201 +212 @@
     }
 
-    m->SetFocalLength(F);
     m->SetPosition(pos);
     m->SetNorm(norm);
+    m->SetShape(shape);
+    m->SetFocalLength(F);
     m->SetSigmaPSF(psf>=0 ? psf : defpsf);
 
@@ -223 +235 @@
 //     x y z nx ny nz F [psf] Type ...
 //
-//  x:      x-coordinate of the mirror's center
-//  y:      y-coordinate of the mirror's center
-//  z:      z-coordinate of the mirror's center
-//  nx:     x-component of the normal vecor in the mirror center
-//  ny:     y-component of the normal vecor in the mirror center
-//  nz:     z-component of the normal vecor in the mirror center
-//  F:      Focal distance of a spherical mirror
-//  [psf]:  This number is the psf given in the units of x,y,z and
-//          defined at the focal distance F. It can be used to overwrite
-//          the second argument given in ReadFile for individual mirrors.
-//  Type:   A instance of a mirrot of the class Type MMirrorType is created
-//          (Type can be, for example, Hex for for MMirrorHex).
-//  ...:    Additional arguments as defined in MMirrorType::ReadM
+//  x:         x-coordinate of the mirror's center
+//  y:         y-coordinate of the mirror's center
+//  z:         z-coordinate of the mirror's center
+//  nx:        x-component of the normal vecor in the mirror center
+//  ny:        y-component of the normal vecor in the mirror center
+//  nz:        z-component of the normal vecor in the mirror center
+//  [shape]:   S for spherical <default>, P for parabolic
+//  F:         Focal distance of a spherical mirror
+//  [psf]:     This number is the psf given in the units of x,y,z and
+//             defined at the focal distance F. It can be used to overwrite
+//             the second argument given in ReadFile for individual mirrors.
+//  Type:      A instance of a mirrot of the class Type MMirrorType is created
+//             (Type can be, for example, Hex for for MMirrorHex).
+//  ...:       Additional arguments as defined in MMirrorType::ReadM
 //
 //

trunk/MagicSoft/Mars/msimreflector/MSimReflector.cc

@@ -18 +18 @@
 !   Author(s): Thomas Bretz,  1/2009 <mailto:tbretz@astro.uni-wuerzburg.de>
 !
-!   Copyright: CheObs Software Development, 2000-2009
+!   Copyright: CheObs Software Development, 2000-2010
 !
 !
@@ -33 +33 @@
 //  obsolete except they can later be "moved" inside the detector, e.g.
 //  if you use MSimPSF to emulate a PSF by moving photons randomly
-//  on the focal plane.
+//  on the focal plane. To switch off this check set detector margin to -1.
 //
 //////////////////////////////////////////////////////////////////////////////
@@ -206 +206 @@
 //    Focal length: F=R/2                        |  Focal length: F = 1/4p
 //                                               |
-//   r^2 + (z-2*F)^2 = (2*F)^2                   |            z = F/4*r^2
+//   r^2 + (z-2*F)^2 = (2*F)^2                   |           z = r^2/4F
 //                                               |
 //          z        = -sqrt(4*F*F - r*r) + 2*F  |
 //          z-2*F    = -sqrt(4*F*F - r*r)        |
 //         (z-2*F)^2 = 4*F*F - r*r               |
-//   z^2-4*F*z+4*F^2 = 4*F*F - r*r  (4F^2-r^2>0) |  z - F/4*r^2 = 0
+//   z^2-4*F*z+4*F^2 = 4*F*F - r*r  (4F^2-r^2>0) |  z - r^2/4F = 0
 //   z^2-4*F*z+r^2   = 0
 //
@@ -223 +223 @@
 // ------------------
 //
-//   z^2*(1+|v|^2) - 2*z*(2*F+q*v) + |q|^2 =  0  |  z^2*|v|^2 - 2*(2/F+q*v)*z + |q|^2 = 0
+//   z^2*(1+|v|^2) - 2*z*(2*F+q*v) + |q|^2 =  0  |  z^2*|v|^2 - 2*z*(2*F+q*v) + |q|^2 = 0
 //
 //                                   z = (-b +- sqrt(b*b - 4ac))/(2*a)
 //
 //              a =   1+|v|^2                    |             a = |v|^2
-//              b = - 2*(2*F+q*v)                |             b = - 2*(2/F+q*v)
+//              b = - 2*b'  with  b' = 2*F+q*v   |             b = - 2*b'  with  b' = 2*F+q*v
 //              c =   |q|^2                      |             c = |q|^2
 //                                               |
@@ -234 +234 @@
 //                                        substitute b := 2*b'
 //
-//                            z = (-2*b' +- 2*sqrt(b'*b' - ac))/(2*a)
-//                            z = (-  b' +-   sqrt(b'*b' - ac))/a
-//                            z = (-b'/a +-   sqrt(b'*b' - ac))/a
+//                            z = (2*b' +- 2*sqrt(b'*b' - ac))/(2*a)
+//                            z = (  b' +-   sqrt(b'*b' - ac))/a
+//                            z = (b'/a +-   sqrt(b'*b' - ac))/a
 //
 //                                        substitute f := b'/a
 //
-//                                      z = (-f +- sqrt(f^2 - c/a)
+//                                      z = f +- sqrt(f^2 - c/a)
 //
 // =======================================================================================
@@ -265 +265 @@
 
     // Find the incident point of the vector to the mirror
-    // u corresponds to downwaqrd going particles, thus we use -u here
+    // u corresponds to downward going particles, thus we use -u here
     const Double_t b = G - q*v;
     const Double_t a = v.Mod2();
@@ -271 +271 @@
 
     // Solution for q spherical (a+1) (parabolic mirror (a) instead of (a+1))
-    const Double_t f = b/(a+1);
-    const Double_t g = c/(a+1);
+    const Double_t A = fShape ? a : a+1;
+
+    const Double_t f = b/A;
+    const Double_t g = c/A;
 
     // Solution of second order polynomial (transformed: a>0)
     // (The second solution can be omitted, it is the intersection
     //  with the upper part of the sphere)
-    //    const Double_t dz = a==0 ? c/(2*b) : f - TMath::Sqrt(f*f - g);
-    const Double_t z = f - TMath::Sqrt(f*f - g);
+    const Double_t z = a==0 ? c/(2*b) : f - TMath::Sqrt(f*f - g);
 
     // Move the photon along its trajectory to the x/y plane of the
@@ -293 +294 @@
 
     // Get normal vector for reflection by calculating the derivatives
-    // of a spherical mirror along x and y
-    const Double_t d = TMath::Sqrt(G*G - p.R2());
-
-    // This is a normal vector at the incident point
+    // of a the mirror's surface along x and y
+    const Double_t d = fShape ? G : TMath::Sqrt(G*G - p.R2());
+
+    // The solution for the normal vector is
+    //    TVector3 n(-p.X()/d, -p.Y()/d, 1));
+    // Since the normal vector doesn't need to be of normal
+    // length we can avoid an obsolete division
     TVector3 n(p.X(), p.Y(), -d);
-    // This is the obvious solution for the normal vector
-    //  TVector3 n(-p.X()/d, -p.Y()/d, 1));
 
     if (fSigmaPSF>0)
@@ -540 +542 @@
         // It is detector specific not reflector specific
         // Discard all photons which definitly can not hit the detector surface
-        if (!fGeomCam->HitDetector(p, fDetectorMargin))
+        if (fDetectorMargin>=0 && !fGeomCam->HitDetector(p, fDetectorMargin))
             continue;
 

trunk/MagicSoft/Mars/msimreflector/MSimReflector.h

@@ -49 +49 @@
     void SetNameReflector(const char *name="MReflector") { fNameReflector = name; }
 
+    void SetDetectorMargin(Double_t m=0) { fDetectorMargin = m; }
+
     ClassDef(MSimReflector, 0) // Task to calculate reflection on a mirror
 };