...

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

Results

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

...

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

TODO

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