From 39e5ce41c2b6b2bcfb8bc84dcf8a518bbb26aa96 Mon Sep 17 00:00:00 2001 From: tiwillia <tiwillia> Date: Tue, 1 Mar 2016 17:18:52 +0000 Subject: [PATCH] Bug fix to correctly calculate uncertainties on theta0 parameter git-svn-id: svn+ssh://svn.cern.ch/reps/lhcb/Geant4/trunk@202377 4525493e-7705-40b1-a816-d608a930855b --- .../extended/electromagnetic/G4TestEm5/doc/release.notes | 2 ++ .../electromagnetic/G4TestEm5/srcnew/src/RMSPlotter.cc | 3 ++- 2 files changed, 4 insertions(+), 1 deletion(-) diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/doc/release.notes b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/doc/release.notes index d2ac698..9055213 100644 --- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/doc/release.notes +++ b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/doc/release.notes @@ -3,6 +3,8 @@ ! Responsible : Timothy Williams ! Purpose : Standalone Geant4 example used to test and validate multiple cattering !----------------------------------------------------------------------------- +! 2016-03-01 - Timothy Williams + Bug fix to correctly calculate uncertainy on theta0 parameter ! 2015-11-02 - Marco Clemencic - Added .gitignore file diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/srcnew/src/RMSPlotter.cc b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/srcnew/src/RMSPlotter.cc index bdbcf12..22abcf9 100644 --- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/srcnew/src/RMSPlotter.cc +++ b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/srcnew/src/RMSPlotter.cc @@ -46,6 +46,7 @@ TGraphErrors* RMSPlotter::MakePlots(){ std::ofstream outt("Results_Table.txt"); outt<<"Energy/MeV "<<"RMS of Scattering Angle/mRad "<<" Uncertainty/mRad" <<std::endl; int D= DataPoints.size(); + double dD=D; Double_t* E=new Double_t[D]; Double_t* R=new Double_t[D]; Double_t* RU=new Double_t[D]; @@ -55,7 +56,7 @@ TGraphErrors* RMSPlotter::MakePlots(){ for(;Di!=DataPoints.end();++Di,++a){ E[a]=(*Di).Energy; R[a]=(*Di).rms=(*Di).Distribution->GetMean(); - RU[a]=(*Di).rmsUncert=(*Di).Distribution->GetRMS(); + RU[a]=(*Di).rmsUncert=((*Di).Distribution->GetRMS())/dD; EU[a]=0; outt<<E[a]<<" , "<<R[a]<<" , "<<RU[a]<<std::endl; } -- GitLab