---

Spectrum.cxx

Go to the documentation of this file.

// Spectrum.cxx: implementation of the Spectrum class.
//

#include "Spectrum.h"
#include <cmath>

// CLHEP
#include "CLHEP/Random/RandomEngine.h"
#include "CLHEP/Random/RandGeneral.h"
#include "CLHEP/Random/RandExponential.h"
#include "CLHEP/Random/RandFlat.h"

// Construction/Destruction

Spectrum::~Spectrum(){}


double Spectrum::flux (double time ) const {
    return 0.; // flag that we don't have a flux
}

double Spectrum::solidAngle( )const
{
    return 1.0; // flag that doesn't calculate
}

std::pair<float,float> Spectrum::dir(float energy)const
{
    // Purpose: return solid angle pair (costh, phi) for the given energy
    // Input:: the given energy.
    // default: random except for Earth occultation
    //here's an attempt at random default distributions as above:
    return std::make_pair(((RandFlat::shoot(1.0))*1.4)-0.4,(RandFlat::shoot(1.0))*2*M_PI);
    
}


const char * Spectrum::particleName()const{
    static const char x='p';
    return &x;
}

double Spectrum::energySrc(HepRandomEngine* engine, double time)
{
    // default implementation, which works for other Spectrum objects
    return (*this)(engine->flat());
}


std::pair<double,double> Spectrum::dir(double energy, HepRandomEngine* engine)
{
    // default that needs fixing!
    return dir(energy);
}

void Spectrum::parseParamList(std::string input, std::vector<float>& output) const
{
    
    int i=0;
    for(;!input.empty() && i!=std::string::npos;){
        float f = ::atof( input.c_str() );
        output.push_back(f);
        i=input.find_first_of(",");
        input= input.substr(i+1);
    } 
}


double Spectrum::interval (double time)
{
    double  r = (solidAngle()*flux(time)*6.);
    
    if (r == 0){ return -1.;
    }else{  
        double p = RandFlat::shoot(1.);
        return (-1.)*(log(1.-p))/r;
    }
}

---