calibration_parameters.C

Go to the documentation of this file.

void Write() {
  using namespace Diana;
  QGlobalDataManager dm;
  dm.SetOwner("calibration_parameters");
 
  QCalibrationParameters calpar;
  TF1 f("pol2","[0]+[1]*x+[2]*x*x");
  f.SetParameter(0,0);
  f.SetParameter(1,1.5);
  f.SetParameter(2,0.3);
  calpar.fFunction = f;
  calpar.fSource = "Th";
 
  QCalibrationResiduals res2615;
  QCalibrationResiduals res911;
  QCalibrationResiduals res511;
 
  res2615.PeakEnergy = 2614.533;
  res2615.Residual = 0.5;
  res2615.ErrorResidual = 0.1;
  calpar.fResiduals.push_back(res2615);
 
  res911.PeakEnergy = 911.5;
  res911.Residual = 0.5;
  res911.ErrorResidual = 0.1;
  calpar.fResiduals.push_back(res911);
 
  res511.PeakEnergy = 511.0;
  res511.Residual = 0.5;
  res511.ErrorResidual = 0.1;
  calpar.fResiduals.push_back(res511);
 
  QCalibrationParametersHandle paramhandle(1,6034);
  paramhandle.Set(calpar);
  cout<<calpar<<endl;
 
  dm.Set(&paramhandle,"calibration_parameters.root");
}
 
void Read(const char* input) {
  using namespace Diana;
  QGlobalDataManager dm;
  dm.SetOwner("calibration_parameters");
  QCalibrationParametersHandle paramhandle(1,6034);
  dm.Get("calibration_parameters",&paramhandle,input);
 
 
  if(!paramhandle.IsValid()) {
    cout<<paramhandle.GetError()<<endl;
    return;
  }
 const QCalibrationParameters& params = paramhandle.Get(); 
 cout<<params<<endl;
}