Modifs en vue de l'ajustement des baselines, Reza 11/02/2019

cxbeam.h

@@ -15,23 +15,28 @@ using namespace SOPHYA;
 class CxBeam  {  //  Auto-correlation beam 
 public:
   CxBeam(ACBeam & a1, ACBeam & a2, Vector3d & baseline)
-    : a1_(a1), a2_(a2), baseline_(baseline)
+    : a1_(a1), a2_(a2), base_baseline_(baseline), baseline_(baseline)
   {
     if (fabs(a1.getLambda()-a2.getLambda())>1e-9)
       throw ParmError("CxBeam::CxBeam() - not same lambda in the two different beams !");
     facphase_=2.*M_PI/a1.getLambda();
   }
   CxBeam(CxBeam const & c)
-    : a1_(c.a1_), a2_(c.a2_), baseline_(c.baseline_), facphase_(c.facphase_)
+    : a1_(c.a1_), a2_(c.a2_), base_baseline_(c.base_baseline_), baseline_(c.baseline_), facphase_(c.facphase_)
   {
   }
     // copy operator 
   CxBeam& operator = (CxBeam const& c)
   {
-    a1_=c.a1_;   a2_=c.a2_;   baseline_=c.baseline_;   facphase_=c.facphase_;
+    a1_=c.a1_;   a2_=c.a2_;   base_baseline_=c.base_baseline_; baseline_=c.baseline_;   facphase_=c.facphase_;
     return (*this);
   }
 
+  inline void ShiftBaseline(Vector3d & shift_baseline)
+  {
+    baseline_ = base_baseline_;
+    baseline_.Add(shift_baseline);
+  }
     //--- return the Beam value for a given direction 
   inline complex<double> Value(LongitudeLatitude const& ll)
   {
@@ -50,7 +55,8 @@ public:
   }
 
   ACBeam a1_, a2_;
-  Vector3d baseline_;
+  Vector3d base_baseline_;   // baseline value before shift 
+  Vector3d baseline_;        // current baseline-value = base_baseline_+shift 
   double facphase_; 
 };
 

gcxfit.h

 /*  PAON4 analysis software 
-    Class for fitting auto-correlations (antenna geometry)  
+    Class for fitting cross-correlations : individual phases on each baseline  
     R. Ansari, December 2018                                             */
 
 #ifndef GCXFIT_H_SEEN