GNGen.hh

Go to the documentation of this file.

#ifndef GNGEN_HH
#define GNGEN_HH

#include "TRandom3.h"
#include "TH3F.h"

class GNGen{
public:
   // all masses in units of solar mass
   GNGen(double mSMBH, double mmin, double mmax, double beta=2); 
   GNGen(const GNGen& x);
   ~GNGen();
   
   // 
   void setFreqCutoff(double f);
   
   // prints or saves in a file a list of m1, m2, rp, e values to be used by PN codes! 
   void generateEvents(int n, char* fn=0);  
   
   // generates m1, m2, rp, e (with evolution of orbit to higher frequency)
   // returns leading order merger time estimate (not very accurate)
   double generateEvent(double& m1, double& m2, double& rp, double& e);
   
private:
   double minM, maxM, beta, smbhM;
   TH3F* dGammadmdMdr;
   TRandom3 rnd;
   double freqCutoff;
   void EvolveRa(double m1, double m2, double& rp, double& ra);
};

#endif