bb.h

Go to the documentation of this file.

#ifndef GENBB_DECAY0_BB_H_
#define GENBB_DECAY0_BB_H_ 1

// Standard library:
#include <string>
#include <fstream>

namespace mygsl{
  class rng;
}

namespace genbb {

  class primary_event;

  namespace decay0 {

    struct enrange
    {
      double      ebb1; 
      double      ebb2; 
      double      toallevents; 
      int         levelE; // keV
      int         itrans02; // Flag for spin
      std::string chdspin; // Spin string

      enrange();
      void reset();
    };

    struct genbbpar
    {
      int nevents;
      int ievstart;
      int irndmst;
      int iwrfile;
      std::string chfile;
      int icurrent; // From common/currentev/
      std::ofstream fout;

      genbbpar();
      void reset();
    };

    // [SI - spectral index]
    enum modebb_type
      {
        MODEBB_UNDEF            = 0,  
        MODEBB_0NUBB_0_2N       = 1,  
        MODEBB_1                = MODEBB_0NUBB_0_2N,
        MODEBB_0NUBB_RHCL_0_2N  = 2,  
        MODEBB_2                = MODEBB_0NUBB_RHCL_0_2N,
        MODEBB_0NUBB_RHCL_02_N  = 3,  
        MODEBB_3                = MODEBB_0NUBB_RHCL_02_N,
        MODEBB_2NUBB_0_2N       = 4,  
        MODEBB_4                = MODEBB_2NUBB_0_2N,
        MODEBB_0NUBB_MAJ1_0_2N  = 5,  
        MODEBB_5                = MODEBB_0NUBB_MAJ1_0_2N,
        MODEBB_0NUBB_MAJ3_0_2N  = 6,  
        MODEBB_6                = MODEBB_0NUBB_MAJ3_0_2N,
        MODEBB_0NUBB_RHCL_2_2N  = 7,  
        MODEBB_7                = MODEBB_0NUBB_RHCL_2_2N,
        MODEBB_2NUBB_2_2N       = 8,  
        MODEBB_8                = MODEBB_2NUBB_2_2N,
        MODEBB_0NUKB_02         = 9,  
        MODEBB_9                = MODEBB_0NUKB_02,
        MODEBB_2NUKB_02         = 10, 
        MODEBB_10               = MODEBB_2NUKB_02,
        MODEBB_0NUKK_02         = 11, 
        MODEBB_11               = MODEBB_0NUKK_02,
        MODEBB_2NUKK_02         = 12, 
        MODEBB_12               = MODEBB_2NUKK_02,
        MODEBB_0NUBB_MAJ7_0_2N  = 13, 
        MODEBB_13               = MODEBB_0NUBB_MAJ7_0_2N,
        MODEBB_0NUBB_MAJ2_0_2N  = 14, 
        MODEBB_14               = MODEBB_0NUBB_MAJ2_0_2N,
        MODEBB_2NUBB_BOSE_0     = 15, 
        MODEBB_15               = MODEBB_2NUBB_BOSE_0,
        MODEBB_2NUBB_BOSE_2     = 16, 
        MODEBB_16               = MODEBB_2NUBB_BOSE_2,
        MODEBB_0NUBB_RHCE_0     = 17, 
        MODEBB_17               = MODEBB_0NUBB_RHCE_0,
        MODEBB_0NUBB_RHCE_0_NME = 18, 
        MODEBB_18               = MODEBB_0NUBB_RHCE_0_NME,
        MODEBB_MIN              = MODEBB_1,
        MODEBB_MAX              = MODEBB_18
      };

    std::string decay0_dbd_mode_label(int);

    bool decay0_supports_esum_limits(int modebb_);

    struct denrange
    {
      double dens;
      double denf;
      int    mode;
      denrange();
      void reset();
    };

    struct helpbb
    {
      double Zd;
      double Ad;
      double e0;
      double e1;
      helpbb();
      void reset();
    };

    struct eta_nme
    {
      double chi_GTw;
      double chi_Fw;
      double chip_GT;
      double chip_F;
      double chip_T;
      double chip_P;
      double chip_R;
      eta_nme();
      void reset();
    };

    struct bbpars : public enrange , public denrange, public helpbb, public eta_nme
    {
      static const unsigned int SPSIZE = 4300;
      int    modebb;   
      double Qbb;      
      double Edlevel;  
      double EK;       
      double Zdbb;     
      double Adbb;     
      int    istartbb; 
      double spthe1[SPSIZE];
      double spthe2[SPSIZE];
      double spmax;
      bbpars();
      void reset();
    };

    //
    // Output: - time of event's start (sec) (as supposed, =0 here);
    //         - time shift from previous time to calculate
    //           when particle was emitted (sec)
    //           (as supposed, =0 here);
    //         - toallevents
    // V.I.Tretyak, 26.11.1995; 17.02.2009.
    void decay0_bb(mygsl::rng & prng, primary_event & event, void *params);

  } // end of namespace decay0
} // end of namespace genbb

#endif // GENBB_DECAY0_BB_H_
//
// Local Variables: --
// mode: c++ --
// End: --