STARLIGHT/starlight/include/twophotonluminosity.h

   1 ///////////////////////////////////////////////////////////////////////////
   2 //
   3 //    Copyright 2010
   4 //
   5 //    This file is part of starlight.
   6 //
   7 //    starlight is free software: you can redistribute it and/or modify
   8 //    it under the terms of the GNU General Public License as published by
   9 //    the Free Software Foundation, either version 3 of the License, or
  10 //    (at your option) any later version.
  11 //
  12 //    starlight is distributed in the hope that it will be useful,
  13 //    but WITHOUT ANY WARRANTY; without even the implied warranty of
  14 //    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
  15 //    GNU General Public License for more details.
  16 //
  17 //    You should have received a copy of the GNU General Public License
  18 //    along with starlight. If not, see <http://www.gnu.org/licenses/>.
  19 //
  20 ///////////////////////////////////////////////////////////////////////////
  21 //
  22 // File and Version Information:
  23 // $Rev:: 164                         $: revision of last commit
  24 // $Author:: odjuvsla                 $: author of last commit
  25 // $Date:: 2013-10-06 16:18:08 +0200 #$: date of last commit
  26 //
  27 // Description:
  28 //
  29 //
  30 //
  31 ///////////////////////////////////////////////////////////////////////////
  32
  33
  34 #ifndef TWOPHOTONLUMINOSITY_H
  35 #define TWOPHOTONLUMINOSITY_H
  36
  37 #include "nucleus.h"
  38 #include "beam.h"
  39 #include "beambeamsystem.h"
  40 #include "starlightlimits.h"
  41
  42
  43 class twoPhotonLuminosity : public beamBeamSystem
  44 {
  45 public:
  46     twoPhotonLuminosity(beam beam_1, beam beam_2);
  47     ~twoPhotonLuminosity();
  48
  49 protected:
  50
  51
  52
  53 private:
  54     struct difflumiargs
  55     {
  56         twoPhotonLuminosity *self;
  57         double m;
  58         double y;
  59         double res;
  60     };
  61     void twoPhotonDifferentialLuminosity();
  62     double D2LDMDY(double M,double Y,double &Normalize);
  63     double D2LDMDY(double M,double Y) const;
  64     static void * D2LDMDY_Threaded(void *a);
  65
  66     double integral(double Normalize);
  67     double radmul(int N,double *Lower,double *Upper,int NIterMin,int NIterMax,double EPS,double *WK,int NIter,double &Result,double &ResErr,double &NFNEVL,double &Summary);
  68     double integrand(double N,double X[15]);
  69     double Nphoton(double W,double gamma,double Rho);
  70
  71     double _W1; //Energy of photon #1
  72     double _W2; //Energy of photon #2
  73     double _gamma; //Gamma of the system
  74
  75     const unsigned int _nWbins;
  76     const unsigned int _nYbins;
  77
  78     const double _wMin;
  79     const double _yMin;
  80     const double _wMax;
  81     const double _yMax;
  82
  83 };
  84
  85
  86 #endif  // TWOPHOTONLUMINOSITY_H