Current Fit for Ba-133 Runs

These are preliminary results of fitting the Ba-133 spectra using multi-mutau model.

Fit on 23 Jan 2004

The past model had four energy lines (31, 35, 53 and 81 keV). But exactly saying, there are seven lines in this energy region as listed below: 
      Energy[0] = 0.030625 *MeV ;
      Energy[1] = 0.030973 *MeV ;
      Energy[2] = 0.034967 *MeV ;
      Energy[3] = 0.036006 *MeV ;
      Energy[4] = 0.053148 *MeV ;
      Energy[5] = 0.079612 *MeV ;
      Energy[6] = 0.080989 *MeV ;
In the new model, I splited 31 keV into 30.625 and 30.973, 35 keV into 34.967 and 36.006, and so on. This improves the fit (reduced chi^2 = 6.8) especially arround 35 keV line. 
  EXT PARAMETER               CURRENT GUESS      PHYSICAL LIMITS       
  NO.   NAME      VALUE            ERROR       NEGATIVE      POSITIVE  
   1  Voltage      2.00000e+02   constant   
   2  Sigma        5.28268e+00   1.00000e+00   1.00000e+00   1.50000e+01
   3  Gain1        3.88910e+00   1.00000e-02   3.50000e+00   4.50000e+00
   4  Gain2        1.00865e+00   1.00000e-03   9.00000e-01   1.10000e+00
   5  Norm_31      9.71843e+00   5.28645e+00   1.05729e+00   1.05729e+02
   6  Norm_35      2.66071e+00   1.16886e+00   2.33773e-01   2.33773e+01
   7  Norm_53      2.45018e-01   1.16886e-01   2.33773e-02   2.33773e+00
   8  Norm_81      4.79439e+00   1.94881e+00   3.89763e-01   3.89763e+01
   9  Tail_length   2.65553e+01   1.00000e+01   1.00000e+01   1.00000e+03
  10  Tail_ratio   2.07810e-01   5.00000e-02   0.00000e+00   1.00000e+01

*******************************************************
reducedchisq = 6.8481 (DOF=292)
d_h->Integral() = 2.14281,   d_h->Integral(rxmin, rxmax) = 2.11458
m_h->Integral() = 2.1361,   m_h->Integral(rxmin, rxmax) = 2.10222
*******************************************************
The plot of normalization parameter (devided by the branching ratio). The fitted curve is a function written as a*exp(b*pow(x,c)). 
    a = 1.10364e+03
    b = -5.94087e+00
    c = -6.54569e-02
The plot of normalization parameter (devided by the branching ratio). The fitted curve is a function written as a*exp(b*pow(x,c)). 
    a is fixed to 15.0
    b = -1.44823e+01
    c = -1.00927e+00
Return to main BAT response page