---

CompositeDiffuse.h

Go to the documentation of this file.

// CompositeDiffuse.h: interface for the CompositeDiffuse class.
#include "FluxSvc/EventSource.h"
#include "CompositeSource.h"
#include <vector>

#include "dom/DOM_Document.hpp"
#include "dom/DOM_Element.hpp"
#include "xml/Dom.h"
#include "xml/IFile.h"
#include "xml/XmlParser.h"
#include <string>
#include <typeinfo>

//#include <cmath>
//#include <algorithm>
//#include <functional>
#include <fstream>
//#include <iostream>

class FluxSource;

class CompositeDiffuse: public CompositeSource { 
public:
    
    typedef struct{
        double l;
        double b;
        double flux;
        double energyIndex;
    }PointSourceData;

    CompositeDiffuse(double totalFlux):
      m_totalFlux(totalFlux),m_unclaimedFlux(totalFlux)
      {
          fillTable();
          setFileFlux();
      }
      
      ~CompositeDiffuse(){}
      
      
      FluxSource* event (double time);
      
      void CompositeDiffuse::addSource (EventSource* aSource);
      
      void addNewSource();
      
      double remainingFluxInterval();
      
      //sets the total flux integrated over the input file, for later use
      void setFileFlux();
      
      void fillTable();
      
      //stolen from FluxMgr, this collects xml and dtd and prepares an input
      //file for parsing by the DOM parser.
      std::string writeXmlFile(const std::vector<std::string>& fileList);

      char* writeLogHistogram();

      void writeSourceCharacteristic(std::ostream& out);
      
private:
    double m_totalFlux; // The total flux from the entire sky, particles/sec/m^2/steradian.
    double m_unclaimedFlux; // The amount of the flux "unaccounted for" by the known sources
    double getRandomFlux();
    //long double pofi(long double intensity);
    long double logNlogS(long double flux);
    std::vector<std::pair<double,double> > m_logNLogS; // inputted logN.logS graph, to be used for finding random fluxes.
    std::vector<PointSourceData> m_listOfDiffuseSources; //information on the sources which got added randomly.
    double m_minFlux; //the minimum flux to be considered in the spectrum
    double m_maxFlux; //the maximum flux to be considered.
    double m_totalIntegratedFlux; //integrated flux, using input file.
    std::string m_dtd; //the DTD document for the XML data to be read in.
    std::map<std::string, DOM_Element > m_sources;
};

---