Commit f2239925 authored by Adrien Matta's avatar Adrien Matta
Browse files

* fixing issue in Miniball simulation and analysis

        - The Interraction scorer was not cleared, making the angle
          wrong in the analysis
parent 043a5cd1
Pipeline #116665 failed with stages
in 186 minutes and 44 seconds
......@@ -82,4 +82,5 @@ void PS_Interactions::EndOfEvent(G4HCofThisEvent*){
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
void PS_Interactions::clear(){
m_DataVector.clear();
m_InterractionCoordinates->Clear();
}
%%%%%%%%%%%%%%%%%%%%%% S1107 at Triumf %%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Beam
Particle= 79Zn
ExcitationEnergy= 0 MeV
Energy= 435 MeV
SigmaEnergy= 0.448 MeV
SigmaThetaX= 0.01 deg
SigmaPhiY= 0.01 deg
SigmaX= 0.5 mm
SigmaY= 0.5 mm
MeanThetaX= 0 deg
MeanPhiY= 0 deg
MeanX= 0 mm
MeanY= 0 mm
%EnergyProfilePath=
%XThetaXProfilePath=
%YPhiYProfilePath=
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
TwoBodyReaction
Beam= 79Zn
Target= 2H
Light= 1H
Heavy= 80Zn
ExcitationEnergyLight= 0.0 MeV
ExcitationEnergyHeavy= 3.7 MeV
CrossSectionPath= 79Zndp_d5half_3.7MeV.dat CSR
ShootLight= 1
ShootHeavy= 1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Decay 80Zn
Threshold= 3.7
Daughter= 80Zn gamma gamma gamma gamma
BranchingRatio= 0.5
ExcitationEnergy= 0 0.879 0.842 0.482 1.497
Shoot= 1 1 1 1 1
LifeTime= 0 ns
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Decay 80Zn
Threshold= 3.7
Daughter= 80Zn gamma gamma gamma
BranchingRatio= 0.5
ExcitationEnergy= 0 1.497 0.482 1.721
Shoot= 1 1 1 1
LifeTime= 0 ns
......@@ -47,7 +47,7 @@ void Analysis::Init(){
else
ligth = "deuteron";
heavy = "Mg29";
heavy = "Zn79";
Sharc = (TSharcPhysics*) m_DetectorManager -> GetDetector("Sharc");
......@@ -145,8 +145,8 @@ void Analysis::TreatEvent(){
X_Trifoil = P.X();
Y_Trifoil = P.Y();
myReaction->SetThetaCM(0);
myReaction->SetExcitationHeavy(0);
// myReaction->SetThetaCM(0);
// myReaction->SetExcitationHeavy(0);
/************************************************/
......@@ -200,8 +200,8 @@ void Analysis::TreatEvent(){
X_Trifoil = P.X();
Y_Trifoil = P.Y();
myReaction->SetThetaCM(0);
myReaction->SetExcitationHeavy(0);
// myReaction->SetThetaCM(0);
// myReaction->SetExcitationHeavy(0);
/************************************************/
......
EmPhysicsList Option4
DefaultCutOff 1000000
IonBinaryCascadePhysics 0
NPIonInelasticPhysics 0
EmExtraPhysics 0
HadronElasticPhysics 0
StoppingPhysics 0
OpticalPhysics 0
HadronPhysicsINCLXX 0
HadronPhysicsQGSP_BIC_HP 0
Decay 0
Menate_R 0
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment