Skip to content
Snippets Groups Projects
Commit fb78d9ad authored by matta's avatar matta
Browse files

* Fixing bug in TransfertToReasonance and Transfert Event generator

	* Heavy particle was emmitted with the wrong Phi angle, problem is fixed
parent 59b26a71
No related branches found
No related tags found
No related merge requests found
...@@ -19,7 +19,7 @@ TransfertToResonance ...@@ -19,7 +19,7 @@ TransfertToResonance
ResonanceDecayA= 8 ResonanceDecayA= 8
CrossSectionPath= 11Li(d,3He)10He.txt CrossSectionPath= 11Li(d,3He)10He.txt
ShootLight= 1 ShootLight= 1
ShootHeavy= 0 ShootHeavy= 1
ShootDecayProduct= 0 ShootDecayProduct= 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
...@@ -450,8 +450,8 @@ void EventGeneratorTransfert::GenerateEvent(G4Event* anEvent , G4ParticleGun* pa ...@@ -450,8 +450,8 @@ void EventGeneratorTransfert::GenerateEvent(G4Event* anEvent , G4ParticleGun* pa
sin(ThetaLight) * sin(phi), sin(ThetaLight) * sin(phi),
cos(ThetaLight)); cos(ThetaLight));
// Momentum in beam frame for heavy particle // Momentum in beam frame for heavy particle
G4ThreeVector momentum_kineHeavy_beam(sin(ThetaHeavy) * cos(phi), G4ThreeVector momentum_kineHeavy_beam(sin(ThetaHeavy) * cos(phi+pi),
sin(ThetaHeavy) * sin(phi), sin(ThetaHeavy) * sin(phi+pi),
cos(ThetaHeavy)); cos(ThetaHeavy));
// write angles/energy to ROOT file // write angles/energy to ROOT file
......
...@@ -510,8 +510,8 @@ void EventGeneratorTransfertToResonance::GenerateEvent(G4Event* anEvent , G4Part ...@@ -510,8 +510,8 @@ void EventGeneratorTransfertToResonance::GenerateEvent(G4Event* anEvent , G4Part
sin(ThetaLight) * sin(phi), sin(ThetaLight) * sin(phi),
cos(ThetaLight)); cos(ThetaLight));
// Momentum in beam frame for heavy particle // Momentum in beam frame for heavy particle
G4ThreeVector momentum_kineHeavy_beam(sin(ThetaHeavy) * cos(phi), G4ThreeVector momentum_kineHeavy_beam(sin(ThetaHeavy) * cos(phi+pi),
sin(ThetaHeavy) * sin(phi), sin(ThetaHeavy) * sin(phi+pi),
cos(ThetaHeavy)); cos(ThetaHeavy));
////////////////////////////////////////////////// //////////////////////////////////////////////////
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment