ReactorConstants.hh Source File

00001 
00006 #ifndef REACTORCONSTANTS_HH
00007 #define REACTORCONSTANTS_HH
00008 
00013   //
00014   // Common set of constants for simulation.
00015   //
00016   // For neutrino event generation
00017   // XCorder can be changed to 0 to use neglect 1/M effects in the nubar 
00018   // cross section.
00019   // Emin,Emax set the range of thrown energy.
00020   // Rmaxgen sets the largest production radius
00021   //
00022   static const int XCorder = 1;  
00023   static const double Emin=1.85;
00024   static const double Emax=8.15;
00025   static const double RMAXGEN = 265.;
00026 
00027   //
00028   // Default isotope number
00029   //
00030   static const int ISOTOPE = 252;
00031   //
00032   // Number of all particles we can handle
00033   //
00034   static const int Ndim = 200;
00035   //
00036   // throw neutron via Gaussian or track
00037   //
00038   static const double SlowNeutronStep=0.1;
00039   //
00040   // units = MeV/Kelvin
00041   //
00042   static const double KBOLTZMAN = 8.617342e-11;
00043   //
00044   // basics
00045   //
00046   static const double PI = 3.1415928;
00047   static const double MPROTON = 938.27200;
00048   //
00049   // Mneutron-Mproton
00050   //
00051   static const double DELTA = 1.2933318;
00052   static const double MELECTRON = 0.510998902;
00053   static const double MMUON = 105.6583692;
00054   //
00055   // 1/alpha
00056   //
00057   static const double ALPHAINV = 137.03599976;
00058   //
00059   // Gfermi in MeV**(-2)
00060   // 
00061   static const double GFERMI = 1.16639e-11;
00062   static const double NAVOGADRO = 6.02214199e23;
00063   //
00064   //classical electron radius
00065   //
00066   static const double RELECTRON = 2.817940285e-13;  
00067   //
00068   // conversion for cross sections from (Mev)**-2 to cm**2.
00069   //
00070   static const double XcMeVtoCmsqrd = 0.389379292e-21;
00071   //
00072   //  c in cm/ns
00073   //
00074   static const double c = 29.9792458;
00075   //
00076   // Scintillator and mineral oil constants taken from Chooz.
00077   // Only H are C are taken into account.
00078   //
00079   static const double density = 0.85;
00080   static const double A0=1.;
00081   static const double A1=12.;
00082   static const double A2=155.;
00083   static const double A3=157.;
00084   static const double Z0=1.;
00085   static const double Z1=6.;
00086   static const double Z2=64.;
00087   static const double Z3=64.;
00088   //
00089   // H and C mass fractions for scintillator
00090   //
00091   static const double f0=0.122;
00092   static const double f1=0.844;
00093   //
00094   // note these are the abundances of Gd155,Gd157
00095   //
00096   static const double f2=0.148;
00097   static const double f3=0.156;
00098   //
00099   // H and C mass fractions for mineral oil
00100   //
00101   static const double f4=0.133;
00102   static const double f5=0.855;
00103   //
00104   // Ionization potentials for dE/dx
00105   //
00106   static const double I0=19.2e-6;
00107   static const double I1=78.0e-6;  
00108   //
00109   // Gd fraction(%)
00110   //
00111   // NOTE: GdConcentration should not be blindly changed because the
00112   // values for mfpGd, tGd, GdCaptureFraction below are based on a 
00113   // GdConcentration of 0.1
00114   // Instead use ReactorDetector::SetParam_GdConcentration which 
00115   // handles resetting all these values.
00116   //
00117   static const double GdConcentrationRef = 0.1;
00118   //
00119   //
00120   // properties of capture in Gd region
00121   // GdCaptureFraction is the fraction of neutrons captured on Gd
00122   // The program assumes capture on Gd155 or 157.  Gd155frac is the
00123   // fraction pf Gd155 capture, etc.
00124   // Gd155peak is the total energy release from Gd155 capture.
00125   //
00126   static const double Gd155frac = 0.1848;
00127   static const double Gd155peak = 8.536;
00128   static const double Gd157frac = 0.8152;
00129   static const double Gd157peak = 7.937;
00130   //
00131   // This is an ad-hoc guess for the average number of gammas/per Gd.
00132   // The program in practice generates at least two gammas per capture
00133   // and an additional Poisson fluctuated GammasPerGd-2.
00134   //
00135   static const double GammasPerGd = 3.;
00136   //
00137   // From capture on H
00138   // This is the gamma energy from np->D gamma.
00139   //
00140   static const double Hpeak = 2.2246;
00141   //
00142   // For capture on C12
00143   //
00144   static const double Cpeak = 4.9461;
00145   //
00146   // Positron annihilation times from KAMLAND 
00147   //    Y. Kino et al.
00148   //    J. Nucl. Radiochem. Science 1, 63, 200
00149   // Not yet used.
00150   //
00151   static const double ShortPosFrac0 = 0.30;
00152   static const double ShortPosFrac1 = 0.21;
00153   static const double LongPosFrac = 0.49;
00154   static const double ShortPosTime0 = 0.19;
00155   static const double ShortPosTime1 = 0.48;
00156   static const double LongPosTime = 3.41;
00157   //
00158 #endif