Changeset 555 in Sophya for trunk/Poubelle/archTOI.old/galcrosslocator.cc

Timestamp:

Nov 9, 1999, 3:04:05 PM (26 years ago)

Author:

ansari

Message:

portage cxx en cours

File:

• trunk/Poubelle/archTOI.old/galcrosslocator.cc (modified) (4 diffs)

trunk/Poubelle/archTOI.old/galcrosslocator.cc

@@ -9 +9 @@
 #define alphaZenith     "alphaZenith"
 #define deltaZenith     "deltaZenith"
-#define rotSpeed        "rotSpeed"
+#define rotSpeed        "rotSpeed0"
 #define rotSpeedSample1 "rotSpeedSample1"
 #define rotSpeedSample2 "rotSpeedSample2"
@@ -22 +22 @@
 #define alphaBolo       "alphaBolo"
 #define deltaBolo       "deltaBolo"
+#define azimuthSST      "azimuthSST"
+#define elvSST          "elvSST"
+#define alphaSST        "alphaSST"
+#define deltaSST        "deltaSST"
 
+// Note : rotSpeed is an estimate, does not take into account balloon drift.
+// accurate value between two galactic crossings has to be determined afterwards
+// as a derived TOI.
 
 GalCrossLocator::GalCrossLocator() {
-  possibleTOIs.insert(TOI(alphaZenith,      TOI::unspec, "", "hours",   "galcross0"));  
-  possibleTOIs.insert(TOI(deltaZenith,      TOI::unspec, "", "degrees", "galcross0"));  
-  possibleTOIs.insert(TOI(rotSpeed,         TOI::unspec, "", "deg/s",   "galcross0"));  
-  possibleTOIs.insert(TOI(rotSpeedSample1,  TOI::unspec, "", "integer", "galcross0"));  
-  possibleTOIs.insert(TOI(rotSpeedSample2,  TOI::unspec, "", "integer", "galcross0"));  
-  possibleTOIs.insert(TOI(lastCrossSample,  TOI::unspec, "", "integer", "galcross0"));  
-  possibleTOIs.insert(TOI(lastCrossSample,  TOI::unspec, "", "integer", "galcross0"));  
-  possibleTOIs.insert(TOI(nextCrossSample,  TOI::unspec, "", "integer", "galcross0"));  
-  possibleTOIs.insert(TOI(azimuthFPC,       TOI::unspec, "", "degrees", "galcross0"));  
-  possibleTOIs.insert(TOI(elvFPC,           TOI::unspec, "", "degrees", "galcross0"));  
-  possibleTOIs.insert(TOI(alphaFPC,         TOI::unspec, "", "hours",   "galcross0"));  
-  possibleTOIs.insert(TOI(deltaFPC,         TOI::unspec, "", "degrees", "galcross0"));  
-  possibleTOIs.insert(TOI(azimuthBolo,      TOI::all,    "", "degrees", "galcross0"));  
-  possibleTOIs.insert(TOI(elvBolo,          TOI::all,    "", "degrees", "galcross0"));  
-  possibleTOIs.insert(TOI(alphaBolo,        TOI::all,    "", "hours",   "galcross0"));  
-  possibleTOIs.insert(TOI(deltaBolo,        TOI::all,    "", "degrees", "galcross0"));  
+  possibleTOIs.insert(TOI(alphaZenith,      TOI::unspec, "interp flag", "hours",   ""));  
+  possibleTOIs.insert(TOI(deltaZenith,      TOI::unspec, "interp flag", "degrees", ""));  
+  possibleTOIs.insert(TOI(rotSpeed,         TOI::unspec, "interp flag", "deg/s",   "galcross0"));  
+  possibleTOIs.insert(TOI(rotSpeedSample1,  TOI::unspec, "interp flag", "integer", "galcross0"));  
+  possibleTOIs.insert(TOI(rotSpeedSample2,  TOI::unspec, "interp flag", "integer", "galcross0"));  
+  possibleTOIs.insert(TOI(lastCrossSample,  TOI::unspec, "interp flag", "integer", "galcross0"));  
+  possibleTOIs.insert(TOI(nextCrossSample,  TOI::unspec, "interp flag", "integer", "galcross0"));  
+  possibleTOIs.insert(TOI(azimuthFPC,       TOI::unspec, "interp flag", "degrees", "galcross0"));  
+  possibleTOIs.insert(TOI(elvFPC,           TOI::unspec, "interp flag", "degrees", "galcross0"));  
+  possibleTOIs.insert(TOI(alphaFPC,         TOI::unspec, "interp flag", "hours",   "galcross0"));  
+  possibleTOIs.insert(TOI(deltaFPC,         TOI::unspec, "interp flag", "degrees", "galcross0"));  
+  possibleTOIs.insert(TOI(azimuthBolo,      TOI::all,    "interp flag", "degrees", "galcross0"));  
+  possibleTOIs.insert(TOI(elvBolo,          TOI::all,    "interp flag", "degrees", "galcross0"));  
+  possibleTOIs.insert(TOI(alphaBolo,        TOI::all,    "interp flag", "hours",   "galcross0"));  
+  possibleTOIs.insert(TOI(deltaBolo,        TOI::all,    "interp flag", "degrees", "galcross0"));  
+  possibleTOIs.insert(TOI(azimuthSST,       TOI::unspec, "interp flag", "degrees", "galcross0"));  
+  possibleTOIs.insert(TOI(elvSST,           TOI::unspec, "interp flag", "degrees", "galcross0"));  
+  possibleTOIs.insert(TOI(alphaSST,         TOI::unspec, "interp flag", "hours",   "galcross0"));  
+  possibleTOIs.insert(TOI(deltaSST,         TOI::unspec, "interp flag", "degrees", "galcross0"));  
 }
 
@@ -55 +65 @@
 double GalCrossLocator::getValue(long sampleNum, TOI const& toi) {
   if (!canGetValue(sampleNum, toi)) return -1;
+  
+  if (toi.options.find("flag") != toi.options.end()) {
+     int SN1, SN2;
+     tempLocator.getCrossSamples(sampleNum, SN1, SN2);
+     if (sampleNum == SN1 || sampleNum == SN2) return 1;
+     return 0;
+  }
 
-  map<TOI, TOIProducer*> m = neededTOIs[toi];
+  map<TOI, TOIProducer*>  & m = neededTOIs[toi];
   double lat, lon, ts;
   
   for (map<TOI, TOIProducer*>::iterator i = m.begin(); i != m.end(); i++) {
-    TOI inToi = (*i).first;
+    TOI const& inToi = (*i).first;
     TOIProducer* prod =  (*i).second;
     if (inToi.name == "latitude")  lat = prod->getValue(sampleNum, inToi);
@@ -83 +100 @@
   if (toi.name == alphaBolo)       return tempLocator.getAlphaBolo(sampleNum, toi.index);
   if (toi.name == deltaBolo)       return tempLocator.getDeltaBolo(sampleNum, toi.index);
+
+  if (toi.name == azimuthSST)      return tempLocator.getAzimutSST(sampleNum);
+  if (toi.name == elvSST)          return tempLocator.getElvSST(sampleNum);
+  if (toi.name == alphaSST)        return tempLocator.getAlphaSST(sampleNum);
+  if (toi.name == deltaSST)        return tempLocator.getDeltaSST(sampleNum);
 
   if (toi.name == lastCrossSample || toi.name == rotSpeedSample1) {