Changeset 4056 for trunk/MagicSoft/Mars/mtemp/mifae/library

Timestamp:

05/12/04 18:50:10 (21 years ago)

Author:

rico

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mtemp/mifae/library

Files:

• MSrcRotate.cc (modified) (4 diffs)
• MSrcRotate.h (modified) (2 diffs)

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

@@ -149 +149 @@
       // save the number of events, initial and final times
       fNEvts    = fRunHeader->GetNumEvents();
-      fFirstEvt = fEvtHeader->GetDAQEvtNumber();
       fIniTime  = fRunHeader->GetRunStart();
       fFinTime  = fRunHeader->GetRunEnd();
@@ -164 +163 @@
       cout << "Run number: " << fRunHeader->GetRunNumber() << endl;
       cout << "Number of events: " << fNEvts << endl;
-      cout << "First event: " << fFirstEvt << endl;
       cout << "Initial MJD date: " << fIniTime.GetMjd() << endl;
       cout << "Final MJD   date: " << fFinTime.GetMjd() << endl;
@@ -175 +173 @@
   // FIXME: for the time being, this is computed by assuming constant event rate
   MTime eventTime;
-  Double_t newMJD = fIniTime.GetMjd() + (fFinTime.GetMjd()-fIniTime.GetMjd())*(fEvtHeader->GetDAQEvtNumber()-fFirstEvt)/fNEvts;  
-  eventTime.SetMjd(newMJD);  
-  
+  Double_t newMJD = fIniTime.GetMjd() + (fFinTime.GetMjd()-fIniTime.GetMjd())*fEvtHeader->GetDAQEvtNumber()/fNEvts;  
+  eventTime.SetMjd(newMJD);    
+  MTime refTime;
+  refTime.SetMjd(fRefMJD);
+
   // de-rotate the source position
   const MAstroSky2Local Observation(eventTime, *fObservatory);
-  const MAstroSky2Local RefObservation(fRefMJD, *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",
+         fRunHeader->GetRunNumber(),fEvtHeader->GetDAQEvtNumber(),
+         fIniTime.GetMjd(),fFinTime.GetMjd(),fDeltaT,
+         newMJD,fRefMJD,Observation.RotationAngle(fRA,fDEC),
+         RefObservation.RotationAngle(fRA,fDEC));
+  //  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);
   
@@ -186 +196 @@
   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*fSrcPos->GetX()-s*fSrcPos->GetY();
+  Float_t newY = s*fSrcPos->GetX()+c*fSrcPos->GetY();
 
 #ifdef DEBUG

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

@@ -18 +18 @@
 {
 private:
-    MSrcPosCam*    fSrcPos;      //!  Pointer to the source position
-    MDCA*          fDCA;         //!  Pointer to the MDCA object
-    MRawEvtHeader* fEvtHeader;   //!  Pointer to the event header
-    MRawRunHeader* fRunHeader;   //!  Pointer to the run header
-    MObservatory*  fObservatory; //!  Pointer to the MObservatory    
+    MSrcPosCam*    fSrcPos;      //  Pointer to the source position
+    MDCA*          fDCA;         //  Pointer to the MDCA object
+    MRawEvtHeader* fEvtHeader;   //  Pointer to the event header
+    MRawRunHeader* fRunHeader;   //  Pointer to the run header
+    MObservatory*  fObservatory; //  Pointer to the MObservatory    
 
     TString     fSrcPosName;
@@ -31 +31 @@
     Double_t   fRefMJD;      //  [MJ date] reference time for rotation
     UInt_t     fNEvts;       //  Number of events in file
-    UInt_t     fFirstEvt;    //  Event index for first event in file
     MTime      fIniTime;     //  Run initial time 
     MTime      fFinTime;     //  Run final time