...

f1 = ROOT.TF1("f1","([2]/pow(x,[3]))/(1+exp(([0]-x)/[4]))+[1]",300,2000)

Results

the following plot plots show the evolution of the fitted parameter 0 as a function of time. This parameter is likely to described the mean FLE (FHE) threshold.

• the The red cross are the obtained values and the error bars are the errors given by MINUIT
• First plot is the fitted FLE threshold as a function of time and the second plot is the same but for FHE threshold
  

• FLE : The mean value is 98 Mev not 100MeV : due to the method?
• FLE FHE : mean is > 1000 MeV
  We clearly see an oscillation in both plot

TODO

• check the result with the CALCalib code
• Make a cron job and run it daily. The results will be stored in a accessible zone.