---

Flux.cxx

Go to the documentation of this file.

// $Header: /nfs/slac/g/glast/ground/cvs/FluxSvc/src/Flux.cxx,v 1.23 2002/10/07 23:42:20 srobinsn Exp $

// Original author: T. Burnett

#include "Flux.h"

#include "FluxSvc/FluxSource.h"
#include "FluxSvc/EventSource.h"
#include "FluxMgr.h"
#include "GPS.h"

Flux::Flux(std::string name) 
: m_time(0)
, m_flux(0)
{
    m_event = s_mgr->source(name);
}
Flux::~Flux() 
{
    delete m_flux;
}

FluxMgr* Flux::s_mgr=0;

void Flux::mgr(FluxMgr* m){ s_mgr=m;}
// name of the flux
std::string Flux::name() const
{
    return m_flux->name();
}

std::string Flux::title()const 
{
    return m_event->fullTitle();
}


void Flux::generate()
{
    // Purpose and Method: generate a new entry trajectory, set FluxSource, increment local time
    // Inputs  - none
    // Outputs - none
    
    m_flux = m_event->event(time());
    double timepass = m_event->interval(time());
    m_time+= timepass;
    pass(timepass);
}

// the particle generated 
std::string Flux::particleName()const{
    return std::string(m_flux->spectrum()->particleName());
}

// its kinetic energy
double Flux::energy()const
{
    return m_flux->energy();
}

// starting point 
HepPoint3D Flux::launchPoint()const
{
    return m_flux->launchPoint();
}

double Flux::time()const 
{
    return m_time ;
}


void Flux::pass ( double t){
    s_mgr->pass(t);
}

/*GPStime*//*int*/double Flux::gpsTime () const{
return s_mgr->time();
}


// direction
HepVector3D Flux::launchDir()const
{
    return m_flux->launchDir();
}

// rate ( /mm**2 /s)
double Flux::rate()const
{
    return m_event->rate(time());
}

void Flux::setTargetArea( double area)
{
    m_event->totalArea(area);
}

double Flux::targetArea()const
{
    return m_event->totalArea();
}


std::string Flux::findSource()const
{
    return m_event->findSource();
}

int Flux::numSource()const
{
    return m_event->numSource();
    
}

#if 0
void Flux::addFactory( const IFactory* factory ) {
    FactoryTable::instance()->addFactory( factory );
}
#endif

void Flux::addFactory(std::string name, const ISpectrumFactory* factory ) {
    SpectrumFactoryTable::instance()->addFactory(name,factory);
}


Rotation Flux::CELTransform(double time)const{
    return s_mgr->CELTransform(time);
}

Rotation Flux::orientTransform(double time)const{
    return s_mgr->orientTransform(time);
}

Rotation Flux::transformGlastToGalactic(double time)const{
    
    return s_mgr->transformGlastToGalactic(time);
}

void Flux::writeSourceCharacteristic(std::ostream& out){
    m_event->writeSourceCharacteristic(out);
}

---