From d7b8c2a401a756d065a7feb74e1b60ad6486a285 Mon Sep 17 00:00:00 2001 From: Morfouace <pierre.morfouace@gmail.com> Date: Sat, 5 Sep 2020 18:04:12 +0200 Subject: [PATCH] Small modification to ChiNu: --- NPSimulation/Detectors/ChiNu/ChiNu.cc | 5 ++++- Projects/ChiNu/Analysis.cxx | 4 ++-- 2 files changed, 6 insertions(+), 3 deletions(-) diff --git a/NPSimulation/Detectors/ChiNu/ChiNu.cc b/NPSimulation/Detectors/ChiNu/ChiNu.cc index c6f7e6e05..2f2eaa5b7 100644 --- a/NPSimulation/Detectors/ChiNu/ChiNu.cc +++ b/NPSimulation/Detectors/ChiNu/ChiNu.cc @@ -420,7 +420,10 @@ void ChiNu::ReadConfiguration(NPL::InputParser parser){ // Construct detector and inialise sensitive part. // Called After DetecorConstruction::AddDetector Method void ChiNu::ConstructDetector(G4LogicalVolume* world){ - + + G4Material* Air = MaterialManager::getInstance()->GetMaterialFromLibrary("Air"); + //world->SetMaterial(Air); + for (unsigned short i = 0 ; i < m_R.size() ; i++) { G4double wX = m_R[i] * sin(m_Theta[i] ) * cos(m_Phi[i] ) ; diff --git a/Projects/ChiNu/Analysis.cxx b/Projects/ChiNu/Analysis.cxx index c594da765..05863c377 100644 --- a/Projects/ChiNu/Analysis.cxx +++ b/Projects/ChiNu/Analysis.cxx @@ -44,7 +44,7 @@ void Analysis::Init(){ InitInputBranch(); InitOutputBranch(); - my_Reaction = new NPL::Reaction("1n(238U,1n)238U@1.5"); + my_Reaction = new NPL::Reaction("1n(238U,1n)238U@0.75"); neutron = new NPL::Nucleus("1n"); } @@ -83,7 +83,7 @@ void Analysis::TreatEvent(){ double DeltaTheta = atan(89.0/Rdet); double random_ThetaLab = ra.Uniform(init_ThetaLab-DeltaTheta, init_ThetaLab+DeltaTheta); - double dEx = my_Reaction->ReconstructRelativistic(Elab[i], init_ThetaLab); + double dEx = my_Reaction->ReconstructRelativistic(Elab[i], random_ThetaLab); ThetaLab.push_back(random_ThetaLab/deg); Ex.push_back(dEx); -- GitLab