diff --git a/Inputs/EventGenerator/10He.reaction b/Inputs/EventGenerator/10He.reaction index 1987df73649f40bdf632faad34251b0c2a9d64b3..8567e32e3e2aad009e6f1f14e403149c23216e86 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 98cf44b25ad8ec2fac8977cebda7751989d56747..d4bb00cc459f78c8b1e4cb2726bc3cc529202095 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 270f2eadae742b4d751c5211cbee822406692ce2..dab91dec0c35eed7f27f43f8df34e9893adca48a 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)); //////////////////////////////////////////////////