From fb78d9ada2e076aa551be15af7f3844ea379483f Mon Sep 17 00:00:00 2001
From: matta <matta@npt>
Date: Mon, 7 Dec 2009 05:43:28 +0000
Subject: [PATCH] * Fixing bug in TransfertToReasonance and Transfert Event
generator * Heavy particle was emmitted with the wrong Phi angle,
problem is fixed
---
Inputs/EventGenerator/10He.reaction | 2 +-
NPSimulation/src/EventGeneratorTransfert.cc | 4 ++--
NPSimulation/src/EventGeneratorTransfertToResonance.cc | 4 ++--
3 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/Inputs/EventGenerator/10He.reaction b/Inputs/EventGenerator/10He.reaction
index 1987df736..8567e32e3 100644
--- a/Inputs/EventGenerator/10He.reaction
+++ b/Inputs/EventGenerator/10He.reaction
@@ -19,7 +19,7 @@ TransfertToResonance
ResonanceDecayA= 8
CrossSectionPath= 11Li(d,3He)10He.txt
ShootLight= 1
- ShootHeavy= 0
+ ShootHeavy= 1
ShootDecayProduct= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
diff --git a/NPSimulation/src/EventGeneratorTransfert.cc b/NPSimulation/src/EventGeneratorTransfert.cc
index 98cf44b25..d4bb00cc4 100644
--- a/NPSimulation/src/EventGeneratorTransfert.cc
+++ b/NPSimulation/src/EventGeneratorTransfert.cc
@@ -450,8 +450,8 @@ void EventGeneratorTransfert::GenerateEvent(G4Event* anEvent , G4ParticleGun* pa
sin(ThetaLight) * sin(phi),
cos(ThetaLight));
// Momentum in beam frame for heavy particle
- G4ThreeVector momentum_kineHeavy_beam(sin(ThetaHeavy) * cos(phi),
- sin(ThetaHeavy) * sin(phi),
+ G4ThreeVector momentum_kineHeavy_beam(sin(ThetaHeavy) * cos(phi+pi),
+ sin(ThetaHeavy) * sin(phi+pi),
cos(ThetaHeavy));
// write angles/energy to ROOT file
diff --git a/NPSimulation/src/EventGeneratorTransfertToResonance.cc b/NPSimulation/src/EventGeneratorTransfertToResonance.cc
index 270f2eada..dab91dec0 100644
--- a/NPSimulation/src/EventGeneratorTransfertToResonance.cc
+++ b/NPSimulation/src/EventGeneratorTransfertToResonance.cc
@@ -510,8 +510,8 @@ void EventGeneratorTransfertToResonance::GenerateEvent(G4Event* anEvent , G4Part
sin(ThetaLight) * sin(phi),
cos(ThetaLight));
// Momentum in beam frame for heavy particle
- G4ThreeVector momentum_kineHeavy_beam(sin(ThetaHeavy) * cos(phi),
- sin(ThetaHeavy) * sin(phi),
+ G4ThreeVector momentum_kineHeavy_beam(sin(ThetaHeavy) * cos(phi+pi),
+ sin(ThetaHeavy) * sin(phi+pi),
cos(ThetaHeavy));
//////////////////////////////////////////////////
--
GitLab