Changeset 8636 for trunk/MagicSoft/Mars/mpointing

Timestamp:

07/15/07 19:25:45 (18 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mpointing

Files:

• MSrcPosCalc.cc (modified) (12 diffs)
• MSrcPosCalc.h (modified) (3 diffs)

trunk/MagicSoft/Mars/mpointing/MSrcPosCalc.cc

@@ -42 +42 @@
 // is taken from MGeomCam. The time is taken from MTime, and the
 // coordinates of the observatory from MObservatory.
-// 
+//
+// FIXME: Add here some information about the more-cycle calculation
+//
 // Input Container:
 //   MPointingPos
@@ -65 +67 @@
 #include "MSrcPosCalc.h"
 
+#include <TRandom.h>
 #include <TVector2.h>
 
@@ -96 +99 @@
     : fObservatory(NULL), fPointPos(NULL), fSourcePos(NULL), fDeviation(NULL),
     fSrcPosCam(NULL), fSrcPosAnti(NULL), fGeom(NULL), fTime(NULL), fCallback(NULL),
-    fMode(kDefault)
+    fMode(kDefault), fNumRandomOffPositions(0)
 {
     fName  = name  ? name  : "MSrcPosCalc";
@@ -157 +160 @@
 Int_t MSrcPosCalc::PreProcess(MParList *pList)
 {
+    // Reset fixed position
+    fFixedPos = TVector2();
+
     fSrcPosCam = (MSrcPosCam*)pList->FindCreateObj("MSrcPosCam");
     if (!fSrcPosCam)
@@ -172 +178 @@
             *fLog << inf;
             *fLog << "MSourcePos [MPointPos] not found... The source position" << endl;
-            *fLog << "set in MSrcPosCam will be set to  (0/0)  or in the case" << endl;
+            *fLog << "set in MSrcPosCam will be set to (0/0)  or in the case" << endl;
             *fLog << "of  Monte Carlo set to the appropriate wobble position." << endl;
             return kTRUE;
@@ -199 +205 @@
     *fLog << "Source Position: " << GetRaDec(*fSourcePos) << endl;
     if (fCallback)
-        *fLog << "Using " << fCallback->GetNumPasses() << " off-regions." << endl;
+        *fLog << "Calculating " << fCallback->GetNumPasses() << " off-regions." << endl;
+    if (fNumRandomOffPositions)
+        *fLog << "Calculating " << fNumRandomOffPositions << " random off-regions." << endl;
 
     // For the case ReInit is never called we try:
@@ -211 +219 @@
 // --------------------------------------------------------------------------
 //
+// Divide the 360deg into num+1 positions. Rotate the given vector
+// by pass+1 steps and return the result.
+//
+TVector2 MSrcPosCalc::Rotate(TVector2 v, Int_t pass, Int_t num) const
+{
+    const Double_t step = TMath::TwoPi()/(num+1);
+    return v.Rotate(step*(pass+1));
+}
+
+// --------------------------------------------------------------------------
+//
 // If fIsWobbleMode==kFALSE set source position to v and anto-source
 // position to -v, if fIsWobbleMode==kTRUE vice versa.
@@ -218 +237 @@
     if (fMode==kWobble)
     {
+        // The trick here is that the anti-source position in case
+        // of the off-source regions is always the on-source positon
+        // thus a proper anti-source cut is possible.
         fSrcPosAnti->SetXY(v);
-
         if (fCallback)
-        {
-            const Double_t step = TMath::TwoPi()/(fCallback->GetNumPasses()+1);
-            v = v.Rotate(step*(fCallback->GetNumPass()+1));
-        }
+            v = Rotate(v, fCallback->GetNumPass(), fCallback->GetNumPasses());
         else
             v *= -1;
-
         fSrcPosCam->SetXY(v);
     }
     else
     {
+        // Because we don't process this three times like in the
+        // wobble case we have to find another way to determine which
+        // off-source region is the right anti-source position
+        // We do it randomly.
         fSrcPosCam->SetXY(v);
-        v *= -1;
+        if (fNumRandomOffPositions>1)
+            v = Rotate(v, gRandom->Integer(fNumRandomOffPositions), fNumRandomOffPositions);
+        else
+            v *= -1;
         fSrcPosAnti->SetXY(v);
     }
@@ -247 +271 @@
     if (fMode==kOffData)
     {
-        SetSrcPos(TVector2());
+        // Set fixed position to camera center
+        fFixedPos = TVector2();
+
+        // It is set here for the cases in which Process is not called at all
+        fSrcPosCam->SetXY(fFixedPos);
         return kTRUE;
     }
@@ -284 +312 @@
     }
 
+    // Determine Monte Carlo position from Monte Carlo header
     TVector2 v(0, 0);
     if (h->GetWobbleMode()>0.5)
@@ -290 +319 @@
         v.Set(-120., 0.);
 
-    SetSrcPos(v);
-
+    // Set fixed position
+    fFixedPos = v;
+
+    // It is set here for the cases in which Process is not called at all
+    fSrcPosCam->SetXY(fFixedPos);
+
+    // Informal message output
     *fLog << inf;
-    *fLog << "Source Position set to x=" << fSrcPosCam->GetX() << "mm ";
-    *fLog << "y=" << fSrcPosCam->GetY() << "mm" << endl;
+    *fLog << "Source Position set to x=" << fFixedPos.X() << "mm ";
+    *fLog << "y=" << fFixedPos.Y() << "mm" << endl;
 
     return kTRUE;
@@ -350 +384 @@
     if (fRunType==MRawRunHeader::kRTMonteCarlo || !fSourcePos || !fTime || !fObservatory || fMode==kOffData)
     {
-        // If this is MC data do not change source position
-        if (fRunType==MRawRunHeader::kRTMonteCarlo)
-            return kTRUE;
-
-        // For real data reset source position to 0/0
-        SetSrcPos();
+        SetSrcPos(fFixedPos);
         return kTRUE;
     }

trunk/MagicSoft/Mars/mpointing/MSrcPosCalc.h

@@ -23 +23 @@
 public:
     enum Mode_t {
-        kDefault = 0,
-        kOffData = 1,
-        kWobble  = 2
+        kDefault = 0,   // Set source position to on-position
+        kOffData = 1,   // The source position is fixed to (0/0)
+        kWobble  = 2    // The source position is set to the wobble aka. anti-source position depending on the cycle number
     };
 private:
@@ -32 +32 @@
     };
 
-    MObservatory *fObservatory;
-    MPointingPos *fPointPos;
-    MPointingPos *fSourcePos;
-    MPointingDev *fDeviation;
-    MSrcPosCam   *fSrcPosCam;
-    MSrcPosCam   *fSrcPosAnti;
-    MGeomCam     *fGeom;
-    MTime        *fTime;
-    MTaskList    *fCallback;
+    MObservatory *fObservatory;   //! Observatory location
+    MPointingPos *fPointPos;      //! Present pointing position of the telescope in Zd/Az
+    MPointingPos *fSourcePos;     //! Source Postion in sky coordinates
+    MPointingDev *fDeviation;     //! Deviation calculated from starguider data
+    MSrcPosCam   *fSrcPosCam;     //! Output: Source position in the camera
+    MSrcPosCam   *fSrcPosAnti;    //! Output: Anti Source position in the camera
+    MGeomCam     *fGeom;          //! Camera geomety
+    MTime        *fTime;          //! Time of the current event
+    MTaskList    *fCallback;      //! Callback function to get the number of the cycle
 
     UShort_t fRunType;            //! Run Type to decide where to get pointing position from
 
-    Int_t fMode;
+    TVector2 fFixedPos;           //! Fixed source position
+
+    Int_t fMode;                  // Mode how the source position is calculated
+
+    Int_t fNumRandomOffPositions; // Number of possible random off-sourcr position
 
     // MSrcPosCalc
     void     SetSrcPos(TVector2 v=TVector2()) const;
+    TVector2 Rotate(TVector2 v, Int_t pass, Int_t num) const;
     TVector2 CalcXYinCamera(const MVector3 &pos0, const MVector3 &pos) const;
     TString  GetRaDec(const MPointingPos &pos) const;
@@ -71 +76 @@
     void SetMode(Mode_t m=kDefault) { fMode = m; }
     void SetCallback(MTaskList *list) { fCallback=list; }
+    void SetNumRandomOffPositions(Int_t num=0) { fNumRandomOffPositions=num; }
 
     ClassDef(MSrcPosCalc, 0) // Calculates the source position in the camera