Changeset 9362 for trunk/MagicSoft/Mars/msim/MSimPointingPos.cc

Timestamp:

02/22/09 23:36:01 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/msim/MSimPointingPos.cc (modified) (8 diffs)

trunk/MagicSoft/Mars/msim/MSimPointingPos.cc

@@ -25 +25 @@
 //////////////////////////////////////////////////////////////////////////////
 //
-//  MSimPointingPos
+// MSimPointingPos
+//
 //
 // This task is meant to simulate the pointing position (mirror orientation).
@@ -31 +32 @@
 // on the user request (e.g. Wobble mode).
 //
-// WARNING: Currently the telescope orientation is just fixed to the
-//          direction in the run-header (if a view cone was given) or
-//          the direction in the evt-header (if no view cone given)
+//
+// If fOffTragetDistance==0 the telescope is oriented depending on the
+// view cone option. If a view cone was given the orientation is fixed to
+// the main direction around the view cone was produced. If no view
+// cone was given the telescope is oriented parallel to the shower axis.
+//
+// If no view cone option was given and off-target observations are switched
+// on by setting fOffTargetDistance!=0 the poitnting position is calculated:
+//
+//  1) fOffTargetDistance < 0:
+//     The pointing position is randomly distributed in a disk of radius
+//     -fOffTragetDistance. fOffTargetDistance is silently ignored.
+//
+//  2) fOffTargetDistance > 0:
+//     The pointing position is set to a position in the given distance
+//     away from the shower axis. If fOffTargetPhi>=0 it is fixed at
+//     this phi value. For phi<0 it is randomly distributed at distances
+//     fOffTargetDistance. (phi==0 is the direction of positive theta)
+//
 //
 //  Input Containers:
@@ -41 +58 @@
 //  Output Containers:
 //   MPointingPos
-//   PointingCorsika [MPointingPos]
 //
 //////////////////////////////////////////////////////////////////////////////
@@ -68 +84 @@
 //
 MSimPointingPos::MSimPointingPos(const char* name, const char *title)
-    : fRunHeader(0), fEvtHeader(0), fPointingCorsika(0), fPointingLocal(0),
-    fOffTargetDistance(-1), fOffTargetPhi(-1)
+    : fRunHeader(0), fEvtHeader(0), fPointing(0),
+    fOffTargetDistance(0), fOffTargetPhi(-1)
 
 {
@@ -83 +99 @@
 Int_t MSimPointingPos::PreProcess(MParList *pList)
 {
-//    fPointingCorsika = (MPointingPos*)pList->FindCreateObj("MPointingPos", "PointingCorsika");
-//    if (!fPointingCorsika)
-//        return kFALSE;
-
-    fPointingLocal = (MPointingPos*)pList->FindCreateObj("MPointingPos");
-    if (!fPointingLocal)
+    fPointing = (MPointingPos*)pList->FindCreateObj("MPointingPos");
+    if (!fPointing)
         return kFALSE;
 
@@ -105 +117 @@
     }
 
-    // FIXED offset
-    // Diffuse?  ( FOV of camera folded with mirror diameter as Corsika input? )
-    // Hour angle!
-    // Angle to magnetic field!
-
-    if (IsOffTargetObservation())
-    {
-        *fLog << inf;
-        *fLog << "Off target observations switched on with" << endl;
+    if (!IsOffTargetObservation())
+        return kTRUE;
+
+    *fLog << inf;
+    *fLog << "Off target observations switched on with" << endl;
+    if (fOffTargetDistance>0)
+    {
         *fLog <<"   a pointing distance of " << GetOffTargetDistance() << "deg ";
         if (fOffTargetPhi<0)
@@ -120 +130 @@
             *fLog << "and phi=" << GetOffTargetPhi() << "deg." << endl;
     }
+    else
+        *fLog << "   a homogenous distribution up to a distance of " << -GetOffTargetDistance() << "deg " << endl;
 
     return kTRUE;
 }
 
-/*
 Bool_t MSimPointingPos::ReInit(MParList *pList)
 {
+    if (fRunHeader->HasViewCone() && IsOffTargetObservation())
+    {
+        *fLog << warn;
+        *fLog << "WARNING - Combining the view cone option with off-target observations doesn't make sense." << endl;
+        *fLog << "          Option for off-target observations will be ignored." << endl;
+    }
     // FIXME: Check also the enlightened region on the ground!
     return kTRUE;
 }
-*/
+
+void MSimPointingPos::GetDelta(Double_t &dtheta, Double_t &dphi) const
+{
+    if (fOffTargetDistance>0)
+    {
+        dtheta = fOffTargetDistance;
+        dphi   = fOffTargetPhi>0 ? fOffTargetPhi : gRandom->Uniform(TMath::TwoPi());
+    }
+    else
+    {
+        dtheta = TMath::Sqrt(gRandom->Uniform(fOffTargetDistance));
+        dphi   = gRandom->Uniform(TMath::TwoPi());
+    }
+}
 
 // --------------------------------------------------------------------------
@@ -136 +166 @@
 Int_t MSimPointingPos::Process()
 {
-
     // If a view cone is given use the fixed telescope orientation
-    const Bool_t fixed = fRunHeader->HasViewCone();
+    const Bool_t viewcone = fRunHeader->HasViewCone();
 
     // Local sky coordinates (direction of telescope axis)
-    /*const*/ Double_t zd = fixed ? fRunHeader->GetZdMin() : fEvtHeader->GetZd()*TMath::RadToDeg();  // x==north
-    /*const*/ Double_t az = fixed ? fRunHeader->GetAzMin() : fEvtHeader->GetAz()*TMath::RadToDeg();  // y==west
-
-    if (!fixed && IsOffTargetObservation())
-    {
+    Double_t zd = viewcone ? fRunHeader->GetZdMin() : fEvtHeader->GetZd()*TMath::RadToDeg();  // x==north
+    Double_t az = viewcone ? fRunHeader->GetAzMin() : fEvtHeader->GetAz()*TMath::RadToDeg();  // y==west
+
+    if (!viewcone)
+    {
+        Double_t dtheta, dphi;
+        GetDelta(dtheta, dphi);
+
         const Double_t theta = zd*TMath::DegToRad();
         const Double_t phi   = az*TMath::DegToRad();
 
-        /*const*/ TVector3 source;
-        source.SetMagThetaPhi(1, theta, phi);
-
-        /*const*/ TVector3 point;
-        point.SetMagThetaPhi(1, theta+fOffTargetDistance, phi);
-
-        const Double_t delta = fOffTargetPhi>0 ? fOffTargetPhi :
-            gRandom->Uniform(TMath::TwoPi());
-
-        point.Rotate(delta, source);
-
-        zd = point.Theta()*TMath::RadToDeg();
-        az = point.Phi()*TMath::RadToDeg();
-    }
+        TVector3 src, pnt;
+        src.SetMagThetaPhi(1, theta,        phi);
+        pnt.SetMagThetaPhi(1, theta+dtheta, phi);
+
+        pnt.Rotate(dphi, src);
+
+        zd = pnt.Theta()*TMath::RadToDeg();
+        az = pnt.Phi()  *TMath::RadToDeg();
+    }
+
+    // Transform the corsika coordinate system (north is magnetic north)
+    // into the telescopes local coordinate system. Note, that all vectors
+    // are already correctly oriented.
+    az += fRunHeader->GetMagneticFieldAz()*TMath::RadToDeg();
 
     // Setup the pointing position
-//    fPointingCorsika->SetLocalPosition(zd, az/*+fRunHeader->GetMagneticFieldAz()*TMath::RadToDeg()*/);
-    fPointingLocal->SetLocalPosition(zd, az+fRunHeader->GetMagneticFieldAz()*TMath::RadToDeg());
+    fPointing->SetLocalPosition(zd, az);
 
     // Calculate incident angle between magnetic field direction