Skip to content
Snippets Groups Projects
Commit 022122b8 authored by adrien-matta's avatar adrien-matta
Browse files

* Fixing nps compilation

        - missing using naespace std in NPIonIon
        - Fixing warning in neutron wall
parent cca502be
No related branches found
No related tags found
No related merge requests found
...@@ -64,6 +64,7 @@ ...@@ -64,6 +64,7 @@
#include "G4IonInelasticProcess.hh" #include "G4IonInelasticProcess.hh"
#include "G4GeneralSpaceNNCrossSection.hh" #include "G4GeneralSpaceNNCrossSection.hh"
using namespace std;
///////////////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////////////////////
NPIonIonInelasticPhysic::NPIonIonInelasticPhysic(const G4String& name): NPIonIonInelasticPhysic::NPIonIonInelasticPhysic(const G4String& name):
G4VPhysicsConstructor(name) G4VPhysicsConstructor(name)
...@@ -99,7 +100,7 @@ void NPIonIonInelasticPhysic::ConstructProcess() ...@@ -99,7 +100,7 @@ void NPIonIonInelasticPhysic::ConstructProcess()
G4IonsShenCrossSection* ShenCrossSections = new G4IonsShenCrossSection; G4IonsShenCrossSection* ShenCrossSections = new G4IonsShenCrossSection;
G4ComponentGGHadronNucleusXsc* GlauberGribovCrossSection = new G4ComponentGGHadronNucleusXsc; G4ComponentGGHadronNucleusXsc* GlauberGribovCrossSection = new G4ComponentGGHadronNucleusXsc;
G4CrossSectionInelastic* GlauberGribovDataSet = new G4CrossSectionInelastic(GlauberGribovCrossSection); //G4CrossSectionInelastic* GlauberGribovDataSet = new G4CrossSectionInelastic(GlauberGribovCrossSection);
// ****************** // ******************
// **** Elastic **** // **** Elastic ****
......
...@@ -60,9 +60,9 @@ using namespace CLHEP; ...@@ -60,9 +60,9 @@ using namespace CLHEP;
namespace NeutronWall_NS{ namespace NeutronWall_NS{
// Energy and time Resolution // Energy and time Resolution
const double EnergyThreshold = 0.1*MeV; const double EnergyThreshold = 0.1*MeV;
const double ResoTime = 0*ns ; //const double ResoTime = 0*ns ;
const double ResoEnergy = 0*MeV ; //const double ResoEnergy = 0*MeV ;
const double ResoPosition = 0*cm; //const double ResoPosition = 0*cm;
//The size of NS should depend on the distance between NeutronWall and plastic Bar right now //The size of NS should depend on the distance between NeutronWall and plastic Bar right now
double NS_X = 2020.0*mm; double NS_X = 2020.0*mm;
double NS_Y = 2020.0*mm; double NS_Y = 2020.0*mm;
......
...@@ -283,7 +283,7 @@ void Tiara::ReadSensitive(const G4Event* event){ ...@@ -283,7 +283,7 @@ void Tiara::ReadSensitive(const G4Event* event){
double EF = RandGauss::shoot(Info[0],ResoEnergyHyball); double EF = RandGauss::shoot(Info[0],ResoEnergyHyball);
if(EF>EnergyThreshold){ if(EF>EnergyThreshold){
int RingNumber=Info[8]; int RingNumber=Info[8];
RingNumber=abs(RingNumber-17); // RingNumber=abs(RingNumber-17);
Info[8]=RingNumber; Info[8]=RingNumber;
m_EventHyball->SetRingE(Info[7],Info[8],EF); m_EventHyball->SetRingE(Info[7],Info[8],EF);
m_EventHyball->SetRingT(Info[7],Info[8],Info[1]); m_EventHyball->SetRingT(Info[7],Info[8],Info[1]);
......
...@@ -48,9 +48,9 @@ using namespace CLHEP; ...@@ -48,9 +48,9 @@ using namespace CLHEP;
namespace TIARA{ namespace TIARA{
// Energy and time Resolution // Energy and time Resolution
const G4double ResoTime = 0 ; const G4double ResoTime = 0 ;
const G4double ResoEnergyInnerBarrel = 0.058*MeV ;// = 136keV of Resolution // Unit is MeV/2.35 const G4double ResoEnergyInnerBarrel = 0.058*MeV ;// = 136keV FWHM
const G4double ResoEnergyOuterBarrel = 0.058*MeV ;// = 136keV of Resolution // Unit is MeV/2.35 const G4double ResoEnergyOuterBarrel = 0.058*MeV ;// = 136keV FWHM
const G4double ResoEnergyHyball = 0.029*MeV ;// = 70keV of Resolution // Unit is MeV/2.35 const G4double ResoEnergyHyball = 0.029*MeV ;// = 70keV FWHM
const G4double EnergyThreshold = 200*keV; const G4double EnergyThreshold = 200*keV;
......
...@@ -5,7 +5,7 @@ Beam ...@@ -5,7 +5,7 @@ Beam
Particle= 22Na Particle= 22Na
ExcitationEnergy= 0 ExcitationEnergy= 0
Energy= 220 Energy= 220
SigmaEnergy= 1 SigmaEnergy= 0
SigmaThetaX= 0.01 SigmaThetaX= 0.01
SigmaPhiY= 0.01 SigmaPhiY= 0.01
SigmaX= 0.5 SigmaX= 0.5
......
...@@ -119,7 +119,8 @@ void Analysis::TreatEvent(){ ...@@ -119,7 +119,8 @@ void Analysis::TreatEvent(){
Energy = Si_E_TH; Energy = Si_E_TH;
// Evaluate energy using the thickness // Evaluate energy using the thickness
ELab = LightAl.EvaluateInitialEnergy( Energy ,0.4*micrometer , ThetaTHSurface); // ELab = LightAl.EvaluateInitialEnergy( Energy ,0.4*micrometer , ThetaTHSurface);
ELab = Energy;
// Target Correction // Target Correction
ELab = LightTarget.EvaluateInitialEnergy( ELab ,TargetThickness/2., ThetaNormalTarget); ELab = LightTarget.EvaluateInitialEnergy( ELab ,TargetThickness/2., ThetaNormalTarget);
......
/run/beamOn 100000
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