Changeset 39777 for branches/czw_branch/20160809/psastro/src/psastroModelAdjust.c

Timestamp:

Oct 18, 2016, 2:59:30 PM (10 years ago)

Author:

watersc1

Message:

Add rotator parity information to psastro models.

File:

• branches/czw_branch/20160809/psastro/src/psastroModelAdjust.c (modified) (3 diffs)

branches/czw_branch/20160809/psastro/src/psastroModelAdjust.c

@@ -58 +58 @@
     } 
 
+    int rotatorParity = psMetadataLookupS32(&status, recipe, "PSASTRO.MODEL.ROT.PARITY");
+    if (!status) psAbort ("Can't find recipe option PSASTRO.MODEL.ROT.PARITY");
+    psMetadataAddS32 (output->fpa->concepts, PS_LIST_TAIL, "FPA.ROT.PARITY", PS_META_REPLACE, "rotator parity parameter", rotatorParity);
+    
     // get reference chip from name
     pmChip *refChip = pmConceptsChipFromName (output->fpa, refChipName);
@@ -86 +90 @@
         psFree (PT);
     }
-
+    
     // rotate the chip-to-FPA transforms to have 0.0 posangle for refChip; 
     // compensate by rotating fpa to TPA transform
@@ -95 +99 @@
 
     // chipAngle should be refChipAngleNominal @ POSANGLE = 0.0
-    float posAngleOffset = chipAngle - refChipAngleNominal;
-
-    fprintf (stderr, "chipAngle is: %f, at PA = 0.0, it should be %f, rotating model by %f\n", 
-             chipAngle*PS_DEG_RAD, refChipAngleNominal*PS_DEG_RAD, posAngleOffset*PS_DEG_RAD);
+    float posAngleOffset = rotatorParity * (chipAngle - refChipAngleNominal);
+
+    fprintf (stderr, "chipAngle is: %f, at PA = 0.0, it should be %f, rotating model by %f (parity %d)\n", 
+             chipAngle*PS_DEG_RAD, refChipAngleNominal*PS_DEG_RAD, posAngleOffset*PS_DEG_RAD, rotatorParity);
 
     // psMetadataAddF32 (output->fpa->concepts, PS_LIST_TAIL, "FPA.POSANGLE", PS_META_REPLACE, "boresite parameter", posangle);