PointingInfo.cxx

Go to the documentation of this file.

class MsgStream; // needed for Exposure.
#include "Event/MonteCarlo/Exposure.h"
#include "FluxSvc/PointingInfo.h"
#include "astro/GPS.h"
#include "astro/SkyDir.h"
#include "astro/JulianDate.h"
#include "astro/EarthCoordinate.h"
#include "ntupleWriterSvc/INTupleWriterSvc.h"



void PointingInfo::set()
{
    // Purpose: save the current status until the next tick
    using namespace astro;
 
    // The GPS singleton has current time and orientation
    GPS* gps = GPS::instance();
    start = gps->time(); 
    CLHEP::Hep3Vector pos_km = gps->position();
    CLHEP::Hep3Vector location = 1.e3* pos_km; 
    
    // cartesian location of the LAT (in m)
    sc_position[0] = location.x();
    sc_position[1] = location.y();
    sc_position[2] = location.z(); 

    ra_zenith = gps->zenithDir().ra();
    dec_zenith= gps->zenithDir().dec();

    ra_scx =    gps->xAxisDir().ra();
    dec_scx =   gps->xAxisDir().dec();

    ra_scz =    gps->zAxisDir().ra();
    dec_scz =   gps->zAxisDir().dec();

    lat_geo = gps->lat(); 
    lon_geo = gps->lon(); 
    
    // override altitude by using shape of earth; access magnetic stuff
    EarthCoordinate loc = gps->earthpos();
    rad_geo = loc.altitude(); 
    L=loc.L();
    B=loc.B();
    lat_mag = loc.geolat();
    in_saa= loc.insideSAA()? 1:0;
    zenith_scz = 180/M_PI* gps->zenithDir().difference(gps->zAxisDir());

}
//------------------------------------------------------------------------
Event::Exposure* PointingInfo::forTDS()const
{

    Event::Exposure* entry = new Event::Exposure;
    entry->init(start,lat_geo,lon_geo,rad_geo,
        sc_position[0],sc_position[1],sc_position[2],ra_scx,dec_scx,ra_scz,dec_scz);
    return entry;

}
//------------------------------------------------------------------------
void PointingInfo::finish(double stop_time, double live)
{
    stop = stop_time;
    livetime = live;

}
//------------------------------------------------------------------------
void PointingInfo::setFT2Tuple(INTupleWriterSvc* tuple, const std::string& tname)
{
    if( tuple==0 ) return;
    tuple->addItem(tname, "start",     &start);
    tuple->addItem(tname, "stop",      &stop);
    tuple->addItem(tname, "sc_position[3]", sc_position);
    tuple->addItem(tname, "lat_geo",   &lat_geo);
    tuple->addItem(tname, "lat_mag",   &lat_mag);
    tuple->addItem(tname, "lon_geo",   &lon_geo);
    tuple->addItem(tname, "rad_geo",   &rad_geo);
    tuple->addItem(tname, "ra_zenith", &ra_zenith);
    tuple->addItem(tname, "dec_zenith",&dec_zenith);
    tuple->addItem(tname, "B_McIlwain",&B);
    tuple->addItem(tname, "L_McIlwain",&L);
    tuple->addItem(tname, "ra_scz",    &ra_scz);
    tuple->addItem(tname, "dec_scz",   &dec_scz);
    tuple->addItem(tname, "ra_scx",    &ra_scx);
    tuple->addItem(tname, "dec_scx",   &dec_scx);
    tuple->addItem(tname, "in_saa",   &in_saa);
    tuple->addItem(tname, "livetime",  &livetime);
}
//------------------------------------------------------------------------
void PointingInfo::setPtTuple(INTupleWriterSvc* tuple, const std::string& tname)
{

    if( tuple==0 ) return;
    tuple->addItem(tname, "PtTime",   &start);
    tuple->addItem(tname, "PtPos[3]", sc_position);
    tuple->addItem(tname, "PtLat",    &lat_geo);
    tuple->addItem(tname, "PtLon",    &lon_geo);
    tuple->addItem(tname, "PtMagLat", &lat_mag);
    tuple->addItem(tname, "PtAlt",    &rad_geo);
    tuple->addItem(tname, "PtRaz",    &ra_scz);
    tuple->addItem(tname, "PtDecz",   &dec_scz);
    tuple->addItem(tname, "PtRax",    &ra_scx);
    tuple->addItem(tname, "PtDecx",   &dec_scx);
    tuple->addItem(tname, "PtSCzenith", &zenith_scz);
    tuple->addItem(tname, "PtMcIlwainB", &B);
    tuple->addItem(tname, "PtMcIlwainL", &L);

}

---