Changeset 15887

Timestamp:

Dec 19, 2007, 3:33:48 PM (18 years ago)

Author:

eugene

Message:

update boresite model to fit boresite path

Location:

trunk/psastro/src

Files:

• Makefile.am (modified) (1 diff)
• psastroModelAnalysis.c (modified) (4 diffs)
• psastroModelBoresite.c (added)
• psastroModelFitBoresite.c (added)
• psastroStandAlone.h (modified) (1 diff)

trunk/psastro/src/Makefile.am

@@ -27 +27 @@
         psastroModelDataLoad.c      \
         psastroModelAnalysis.c      \
+        psastroModelBoresite.c      \
+        psastroModelFitBoresite.c   \
         psastroCleanup.c
 

trunk/psastro/src/psastroModelAnalysis.c

@@ -45 +45 @@
 
     // temp data vectors
-    psVector *posZero = psVectorAlloc (files->n, PS_TYPE_F32);
-    psVector *Lo      = psVectorAlloc (files->n, PS_TYPE_F32);
-    psVector *Mo      = psVectorAlloc (files->n, PS_TYPE_F32);
+    psVector *posZero  = psVectorAlloc (files->n, PS_TYPE_F32);
+    psVector *Po       = psVectorAlloc (files->n, PS_TYPE_F32);
+    psVector *Lo       = psVectorAlloc (files->n, PS_TYPE_F32);
+    psVector *Mo       = psVectorAlloc (files->n, PS_TYPE_F32);
 
     int n = 0;
@@ -71 +72 @@
         while (posZero->data.F32[n] <   0.0) posZero->data.F32[n] += 360.0;
 
+        Po->data.F32[n] = POSANGLE * PM_RAD_DEG;
         Lo->data.F32[n] = chip->toFPA->x->coeff[0][0];
         Mo->data.F32[n] = chip->toFPA->y->coeff[0][0];
@@ -79 +81 @@
     Lo->n = n;
     Mo->n = n;
+    Po->n = n;
     posZero->n = n;
 
@@ -86 +89 @@
     fprintf (stderr, "pos zero %f +/- %f\n", stats->sampleMedian, stats->sampleStdev);
 
+    psastroModelFitBoresite (Lo, Mo, Po);
+
     return true;
 }

trunk/psastro/src/psastroStandAlone.h

@@ -22 +22 @@
 bool psastroModelAnalysis (pmConfig *config);
 
+bool psastroModelFitBoresite (psVector *Lo, psVector *Mo, psVector *Po);
+psF32 psastroModelBoresite (psVector *deriv, const psVector *params, const psVector *coord);
+
+psF32 psastroModelBoresiteL (psVector *deriv, const psVector *params, const psVector *coord);
+psF32 psastroModelBoresiteM (psVector *deriv, const psVector *params, const psVector *coord);
+
+// these are used to define the boresite model parameters
+# define PAR_L0  0  // Lo = params[0] 
+# define PAR_M0  1  // Mo = params[1] 
+# define PAR_RL  2  // RL = params[2] 
+# define PAR_RM  3  // RM = params[3] 
+# define PAR_P0  4  // P0 = params[4]
+# define PAR_T0  5  // phi = params[4]
+
 #endif