SpecTime.h

Go to the documentation of this file.

#ifndef __SpecTime__
#define __SpecTime__
 
#include <TF1.h>
#include <TH2.h>
#include <TGraphAsymmErrors.h>
#include <TRandom3.h>
#include "Extinction.h"
 
using namespace std;
 
class SpecTime: public Extinction{
  
 public:
  
 
  SpecTime(const int aNt, const double *aTimeBins,
           const int aNe, const double *aEnergyBins);
  
 
  SpecTime(const int aNt, const double aTimeMin, const double aTimeMax, const bool aUseTimeLog,
       const int aNe, const double aEnergyMin, const double aEnergyMax, const bool aUseEnergyLog);
  
 
  SpecTime(TH2D* aFlux, TH2D* aRate=NULL);
  
 
  virtual ~SpecTime(void);  
 
  inline void SetPhotoElectricCrossSection(TGraph *aGraph){
    delete photoelectric_xsec;
    photoelectric_xsec = (TGraph*)aGraph->Clone("sp_xsec");
    photoelectric_xsec->SetBit(TGraph::kIsSortedX);
  };
  
 
  void SetFlux(TH2D *aFlux, TGraph *aArea=NULL);
  
 
  void ApplyExtinction(TGraph *aArea=NULL);
 
 
  void MakeRate(const double aArea);
 
 
  void MakeRate(TGraph *aArea);
 
 
  int GetParticleN(const ULong64_t aTime, const UInt_t aDuration, const double aCosTheta=1.0);
 
 
  inline double GetTimeMax(void){ return flux->GetXaxis()->GetBinUpEdge(flux->GetNbinsX()); };
 
 
  inline double GetTimeMin(void){ return flux->GetXaxis()->GetBinLowEdge(1); };
 
 
  TH1D* GetLightCurve(const double aEnergyMin=0.0, const double aEnergyMax=-1.0);
 
 
  TH1D* GetMagnitude(TGraph *aEnergyTransmission);
 
 
  double GetMagnitude(const double aTimeMin, const double aTimeMax, TGraph *aEnergyTransmission);
 
 
  TH1D* GetSpectrum(const double aTimeMin=0.0, const double aTimeMax=-1.0, const bool aObserved=false);
 
 
  void MakeUniform(const double aFlux, TGraph *aArea=NULL);
 
 
  void MakeBandLimited(const double aFlux,
                       const double aEnergyMin, const double aEnergyMax,
                       TGraph *aArea=NULL);
 
 
  void MakeExpDecayTime(const double aFlux, const double aTimeDecay,
            const double aEnergyMin, const double aEnergyMax,
            TGraph *aArea=NULL);
 
 
  void MakeExpDecayEnergy(const double aFlux, const double aEnergyDecay,
              TGraph *aArea=NULL);
  
 
  void MakeGRBModel0(const double aFlux,
                     const double aTime0, const double aAlpha,
                     const double aEnergy0, const double aBeta,
                     TGraph *aArea=NULL);
 
 
  void MakeGRBModel1(const double aFlux, const double aEnergy0,
              const double aAlpha, const double aBeta,
              TGraph *aArea=NULL);
 
 
  void MakeGRBData(TH1D *aLightCurveExp,
                   const double aE1, const double aE2,
                   const double aPhotonIndex,
                   const double aNint, const double aNgal,
                   const double aZ, TGraph *aAreaExp,
                   TGraph *aArea=NULL);
 
 
  double DrawEnergy(const double aTime);
 
 
  inline TH2D* GetFlux(void){ return flux_x; };
 
 
  inline TH2D* GetRate(void){ return rate; };
 
 protected:
 
  TH2D *flux;                  
  TH2D *flux_x;                
  TH2D *rate;                  
  TRandom3 *randgen;           
 
 private:
  TGraph *photoelectric_xsec;  
  double **rate_cum;           
 
  void MakeCum(void);          
 
  void ConstructHistos(const int aNt, const double *aTimeBins,
                       const int aNe, const double *aEnergyBins);
 
  ClassDef(SpecTime,0)
};
 
#endif