Changeset 4072 for trunk/MagicSoft/Mars/mtemp/mifae/library/MSrcRotate.cc

Timestamp:

05/14/04 16:35:57 (21 years ago)

Author:

rico

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mtemp/mifae/library/MSrcRotate.cc (modified) (10 diffs)

trunk/MagicSoft/Mars/mtemp/mifae/library/MSrcRotate.cc

@@ -26 +26 @@
 // MSrcRotate
 //
-// Task to rotate the position of the source as a function of Azimuth and 
-// Zenith angles
+// Task to rotate the position of the source as a function of time
 //
 //  Input Containers:
@@ -45 +44 @@
 #include "MRawEvtHeader.h"
 #include "MSrcRotate.h"
-#include "MSrcPosCam.h"
-#include "MDCA.h"
 
 #include "MAstroSky2Local.h"
 #include "MObservatory.h"
+#include "MSrcPosCam.h"
+#include "MDCA.h"
 
 #include "MLog.h"
@@ -59 +58 @@
 
 static const TString gsDefName  = "MSrcRotate";
-static const TString gsDefTitle = "De-rotate position of the source";
+static const TString gsDefTitle = "Rotate position of the source";
 
 // -------------------------------------------------------------------------
@@ -68 +67 @@
 //
 MSrcRotate::MSrcRotate(const char* srcPos, const char* dca, const char *name, const char *title)
-  : fSrcPos(NULL), fDCA(NULL), fRA(0), fDEC(0), fRunNumber(0)
+  : fRA(0), fDEC(0), fRefMJD(0), fRunNumber(0)
 {
     fName  = name  ? name  : gsDefName.Data();
     fTitle = title ? title : gsDefTitle.Data();
 
-    fSrcPosName  = srcPos;
-    fDCAName     = dca;
+    SetSrcPosName(srcPos);
+    SetDCAName(dca);
 }
 
@@ -84 +83 @@
 Int_t MSrcRotate::PreProcess(MParList *pList)
 {
-    // look for/create MSrcPosCam
-  fSrcPos = (MSrcPosCam*)pList->FindObject(AddSerialNumber(fSrcPosName), "MSrcPosCam");
-  if (!fSrcPos)
-    {
-      *fLog << warn << AddSerialNumber(fSrcPosName) << " [MSrcPosCam] not found... creating default container." << endl;
-      fSrcPos = (MSrcPosCam*)pList->FindCreateObj("MSrcPosCam", AddSerialNumber(fSrcPosName));
-      if(!fSrcPos)
-        return kFALSE;
-    }
-
-  // look for/create MDCA
-  fDCA = (MDCA*)pList->FindObject(AddSerialNumber(fDCAName), "MDCA");
-  if (!fDCA)
-    {
-      *fLog << warn << AddSerialNumber(fDCAName) << " [MDCA] not found... creating default container." << endl;
-      fDCA = (MDCA*)pList->FindCreateObj("MDCA", AddSerialNumber(fDCAName));
-      if(!fDCA)
-        return kFALSE;
-    }
-
+  if(!SearchForSrcPos(pList))
+    return kFALSE;
+  
   // look for MRawRunHeader
   fRunHeader = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
@@ -141 +123 @@
 // FIXME: for the time being, this is computed by assuming constant event rate
 //
-Int_t MSrcRotate::Process()
+Int_t MSrcRotate::ComputeNewSrcPosition()
 {  
   if(fRunHeader->GetRunNumber()!=fRunNumber)
     {
       fRunNumber=fRunHeader->GetRunNumber();
-
+      
       // save the number of events, initial and final times
       fNEvts    = fRunHeader->GetNumEvents();
@@ -155 +137 @@
       if((ULong_t)TMath::Nint(fIniTime.GetMjd())!=(ULong_t)TMath::Nint(fFinTime.GetMjd()))
         {
-          *fLog << err << "MSrcRotate::Process Error: Inial and final MJDs are different ("<<fIniTime.GetMjd()<<"!="<<fFinTime.GetMjd()<<")" << endl;
+          *fLog << err << "MSrcRotate::ComputeNewSrcPosition Error: Inial and final MJDs are different ("<<fIniTime.GetMjd()<<"!="<<fFinTime.GetMjd()<<")" << endl;
           return kFALSE;
         }
-
+      
 #ifdef DEBUG
       cout << endl << "********************"  << endl;
@@ -169 +151 @@
     }
   
-
+  
   // Compute the event time
   // FIXME: for the time being, this is computed by assuming constant event rate
@@ -177 +159 @@
   MTime refTime;
   refTime.SetMjd(fRefMJD);
-
+  
   // de-rotate the source position
   const MAstroSky2Local Observation(eventTime, *fObservatory);
   const MAstroSky2Local RefObservation(refTime, *fObservatory);
-
+  
 #ifdef DEBUG
   printf("Run:%d, Event:%d, iniMJD=%15.5f, finMJD=%15.5f, fDeltaT=%15.5f, newMJD=%15.5f, fRefMJD=%15.5f, rotation=%15.5f, ref=%15.5f\n",
@@ -190 +172 @@
   //  cout << "newMJD=" << newMJD << ", fRefMJD="<<fRefMJD<<", rotation="<<Observation.RotationAngle(fRA,fDEC)<<", ref="<<RefObservation.RotationAngle(fRA,fDEC)<< endl;
 #endif
-
-  Double_t rotationAngle = Observation.RotationAngle(fRA,fDEC)-RefObservation.RotationAngle(fRA,fDEC);
+  
+  Double_t rotationAngle = fRefMJD ? 
+    Observation.RotationAngle(fRA,fDEC)-RefObservation.RotationAngle(fRA,fDEC) :
+    Observation.RotationAngle(fRA,fDEC) ;
+  MSrcPosCam* srcpos = GetSrcPosCam();
+  MDCA* dca = GetDCA();
   
   Float_t c = TMath::Cos(rotationAngle);
   Float_t s = TMath::Sin(rotationAngle);
   // perform a rotation of -rotationAngle to move the source back to the "initial" position
-  Float_t newX = c*fSrcPos->GetX()-s*fSrcPos->GetY();
-  Float_t newY = s*fSrcPos->GetX()+c*fSrcPos->GetY();
-
+  Float_t newX = c*srcpos->GetX()-s*srcpos->GetY();
+  Float_t newY = s*srcpos->GetX()+c*srcpos->GetY();
+  
 #ifdef DEBUG
   Double_t rotationAngleComp = fObservatory->RotationAngle(0.1256,2.63);
   cout << "Event " << fEvtHeader->GetDAQEvtNumber() << endl;
-  cout << "newMJD=" << newMJD <<", rotationAngle=" << rotationAngle <<", rotationAngleComp=" << rotationAngleComp << ", oldX="<<fIniSrcPos.GetX()<< ", oldY="<<fIniSrcPos.GetY()<< ", newX="<<newX<< ", newY="<<newY << endl,
+  cout << "newMJD=" << newMJD <<", rotationAngle=" << rotationAngle <<", rotationAngleComp=" << rotationAngleComp << ", oldX="<<fIniSrcPos.GetX()<< ", oldY="<<fIniSrcPos.GetY()<< ", newX="<<newX<< ", newY="<<newY << endl;
 #endif
-
-  fSrcPos->SetX(newX);
-  fSrcPos->SetY(newY);
-  fDCA->SetRefPoint(newX,newY);
+  
+  srcpos->SetX(newX);
+  srcpos->SetY(newY);
+  dca->SetRefPoint(newX,newY);
 
   return kTRUE;