source: Sophya/trunk/Poubelle/archTOI.old/galcrosslocator.cc

// galcrosslocator.cc
// Eric Aubourg         CEA/DAPNIA/SPP   octobre 1999

#include "galcrosslocator.h"
#include "archexc.h"
#include "archparam.h"
#include "templocator.h"

#define alphaZenith     "alphaZenith"
#define deltaZenith     "deltaZenith"
#define rotSpeed        "rotSpeed0"
#define rotSpeedSample1 "rotSpeedSample1"
#define rotSpeedSample2 "rotSpeedSample2"
#define lastCrossSample "lastCrossSample"
#define nextCrossSample "nextCrossSample"
#define azimuthFPC      "azimuthFPC"
#define elvFPC          "elvFPC"
#define alphaFPC        "alphaFPC"
#define deltaFPC        "deltaFPC"
#define azimuthBolo     "azimuthBolo"
#define elvBolo         "elvBolo"
#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, "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"));  
}

string GalCrossLocator::getName() {
  return("GalCrossLocator 1.0");
}

bool GalCrossLocator::canGetValue(long sampleNum, TOI const&) {
  int dum;
  return (tempLocator.getCrossSamples(sampleNum, dum, dum) == 0);
}


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];
  double lat, lon, ts;
  
  for (map<TOI, TOIProducer*>::iterator i = m.begin(); i != m.end(); i++) {
    TOI const& inToi = (*i).first;
    TOIProducer* prod =  (*i).second;
    if (inToi.name == "latitude")  lat = prod->getValue(sampleNum, inToi);
    if (inToi.name == "longitude") lon = prod->getValue(sampleNum, inToi);
    if (inToi.name == "tsid")       ts = prod->getValue(sampleNum, inToi);
  }
    
  tempLocator.setEarthPos(lon,lat);
  tempLocator.setTSid(ts);

  if (toi.name == alphaZenith)     return tempLocator.getAlphaZenith();
  if (toi.name == deltaZenith)     return tempLocator.getDeltaZenith();
  if (toi.name == rotSpeed)        return tempLocator.getRotSpeed(sampleNum);

  if (toi.name == azimuthFPC)      return tempLocator.getAzimutCenter(sampleNum);
  if (toi.name == elvFPC)          return tempLocator.getElvCenter(sampleNum);
  if (toi.name == alphaFPC)        return tempLocator.getAlphaCenter(sampleNum);
  if (toi.name == deltaFPC)        return tempLocator.getDeltaCenter(sampleNum);
  
  if (toi.name == azimuthBolo)     return tempLocator.getAzimutBolo(sampleNum, toi.index);
  if (toi.name == elvBolo)         return tempLocator.getElvBolo(sampleNum, toi.index);
  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) {
     int SN1, SN2;
     tempLocator.getCrossSamples(sampleNum, SN1, SN2);
     return SN1;
  } 
  if (toi.name == nextCrossSample || toi.name == rotSpeedSample2) {
     int SN1, SN2;
     tempLocator.getCrossSamples(sampleNum, SN1, SN2);
     return SN2;
  } 
  
  return -1;
}


set<TOI> GalCrossLocator::reqTOIFor(TOI const&) {
  set<TOI> t; 
  t.insert(TOI("latitude",  TOI::unspec, "interp"));
  t.insert(TOI("longitude", TOI::unspec, "interp"));
  t.insert(TOI("tsid",      TOI::unspec));
  return t;
}