Changeset 3892 for trunk/MagicSoft/Mars/manalysis

Timestamp:

04/29/04 16:44:48 (21 years ago)

Author:

tonello

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• MCT1PointingCorrCalc.cc (modified) (8 diffs)
• MCT1PointingCorrCalc.h (modified) (3 diffs)

trunk/MagicSoft/Mars/manalysis/MCT1PointingCorrCalc.cc

@@ -17 +17 @@
 !
 !   Author(s): Wolfgang Wittek  03/2003 <mailto:wittek@mppmu.mpg.de>
-!
+!              Nadia Tonello    05/2003 <mailto:tonello@mppmu.mpg.de>
 !   Copyright: MAGIC Software Development, 2000-2003
 !
@@ -27 +27 @@
 //  MCT1PointingCorrCalc                                                   //
 //                                                                         //
-//  This is a task to do the CT1 pointing correction                       //
-//                                                                         //
+//  This is a task to do the CT1 pointing correction.                      // 
+//                                                                         //
+//  NT: You can correct the Mkn421 data (default setting), or the 1ES1959  //
+//      data.                                                              // 
+//      To change to the correction needed for 1ES1959,  you have to call  //
+//      the member funcion: SetPointedSource                               //
+//                                                                         //
+//  Example:                                                               //
+//      MCT1PointingCorrCalc correct;                                      //
+//      correct.SetPointedSource(MCT1PointingCorrectionCalc::K1959)        //
+//                                                                         // 
 /////////////////////////////////////////////////////////////////////////////
 
@@ -45 +54 @@
 
 using namespace std;
-
 // --------------------------------------------------------------------------
 //
@@ -52 +60 @@
 MCT1PointingCorrCalc::MCT1PointingCorrCalc(const char *srcname,
                                            const char *name, const char *title)
-    : fSrcName(srcname)
+  : fSrcName(srcname), fPoiSource(k421)
 {
     fName  = name  ? name  : "MCT1PointingCorrCalc";
@@ -70 +78 @@
     }
 
-    fHourAngle = (MParameterD*)pList->FindObject("HourAngle", "MParameterD");
+   fHourAngle = (MParameterD*)pList->FindObject("HourAngle", "MParameterD");
     if (!fHourAngle)
     {
@@ -89 +97 @@
 }
 
+
+// --------------------------------------------------------------------------
+//
+//Implemented Daniel Kranich's pointing correction for Mkn421 (2001 data)
+//
+
+void MCT1PointingCorrCalc::PointCorr421()
+{
+   //*fLog << "MCT1PointingCorrCalc::Process; fhourangle = " 
+   //      << fhourangle << endl;
+
+   Float_t fhourangle = fHourAngle->GetVal();
+
+   Float_t cx = -0.05132 - 0.001064 * fhourangle 
+                         - 3.530e-6 * fhourangle * fhourangle;
+   cx /= fMm2Deg;
+
+   Float_t cy = -0.04883 - 0.0003175* fhourangle
+                         - 2.165e-5 * fhourangle * fhourangle;
+   cy /= fMm2Deg;
+
+   fSrcPos->SetXY(cx, cy);
+
+   //*fLog << "MCT1PointingCorrCal::Process; fhourangle, cx, cy, fMm2Deg = "
+   //      << fhourangle << ",  " << cx << ",  " << cy << ",  " 
+   //      << fMm2Deg << endl;
+
+   fSrcPos->SetReadyToSave();
+   return;
+}
+
+
+// --------------------------------------------------------------------------
+//
+// NT :Implemente Daniel Kranich's pointing correction for 1ES1959 (2002 data)
+
+void MCT1PointingCorrCalc::PointCorr1959()
+{
+   //*fLog << "MCT1PointingCorrCalc::Process; fhourangle = " 
+   //      << fhourangle << endl;
+
+   Float_t fhourangle = fHourAngle->GetVal();
+
+   Float_t cx = -0.086 - 0.00091 * fhourangle ;
+   cx /= fMm2Deg;
+
+   Float_t cy = -0.083 - 0.001 * fhourangle ;
+   cy /= fMm2Deg;
+
+   fSrcPos->SetXY(cx, cy);
+
+   //*fLog << "MCT1PointingCorrCal::Process; fhourangle, cx, cy, fMm2Deg = "
+   //      << fhourangle << ",  " << cx << ",  " << cy << ",  " 
+   //      << fMm2Deg << endl;
+
+   fSrcPos->SetReadyToSave();
+   return;
+
+}
 // --------------------------------------------------------------------------
 //
 // Do the pointing correction
 //
-// the parametrization is for Mkn421 2001 data (Daniel Kranich)
-// 
 Int_t MCT1PointingCorrCalc::Process()
 {
@@ -100 +165 @@
    // (cx, cy) is the source position in the camera [mm]
    //
-   Float_t fhourangle = fHourAngle->GetVal();
-
-   //*fLog << "MCT1PointingCorrCalc::Process; fhourangle = " 
-   //      << fhourangle << endl;
-
-   Float_t cx = -0.05132 - 0.001064 * fhourangle 
-                         - 3.530e-6 * fhourangle * fhourangle;
-   cx /= fMm2Deg;
-
-   Float_t cy = -0.04883 - 0.0003175* fhourangle
-                         - 2.165e-5 * fhourangle * fhourangle;
-   cy /= fMm2Deg;
-
-   fSrcPos->SetXY(cx, cy);
-
-   //*fLog << "MCT1PointingCorrCal::Process; fhourangle, cx, cy, fMm2Deg = "
-   //      << fhourangle << ",  " << cx << ",  " << cy << ",  " 
-   //      << fMm2Deg << endl;
-
-   fSrcPos->SetReadyToSave();
-
+   switch (fPoiSource)
+    {
+    case k421:
+        PointCorr421();
+    case k1959:
+        PointCorr1959();
+    }  
    return kTRUE;
 }
@@ -126 +177 @@
 
 
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+

trunk/MagicSoft/Mars/manalysis/MCT1PointingCorrCalc.h

@@ -20 +20 @@
 class MCT1PointingCorrCalc : public MTask
 {
+public:
+    typedef enum {
+        k421,
+        k1959
+    } PointedSource_t;
+
 private:
 
@@ -26 +32 @@
     MParameterD  *fHourAngle;
 
+    PointedSource_t fPoiSource;
     Float_t       fMm2Deg;
 
+    void PointCorr421();
+    void PointCorr1959();
+
 public:
+
     MCT1PointingCorrCalc(const char *srcname="MSrcPosCam",
                          const char *name=NULL, const char *title=NULL);
@@ -35 +46 @@
     Int_t Process();
 
-    ClassDef(MCT1PointingCorrCalc, 0)   // Task to do the CT1 pointing correction for Mkn421 2001 data
+    void SetPointedSource(PointedSource_t s) { fPoiSource = s; }
+ 
+    ClassDef(MCT1PointingCorrCalc, 0)   // Task to do the CT1 pointing correction for Mkn421 2001 data or 1ES1959 2002 data
 };