Message: Proton physic list between 300 - 6600 MeV Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Proton physic list between 300 - 6600 MeV 

Forum: Physics List
Date: 29 Aug, 2011
From: Manuel Castro Avila <Manuel Castro Avila>

Hi Geant4 Users,

I am simulating the interaction the protons with a CZT detector. So, my proton physic list is:

 if (particleName == "proton") {

pmanager->AddProcess(new G4hMultipleScattering, -1, 1, 1);
pmanager->AddProcess(new G4hIonisation(),-1, 2, 2);

G4HadronElasticProcess* theElasticProcess = new G4HadronElasticProcess();
G4LElastic* theElasticModel = new G4LElastic();
theElasticProcess->RegisterMe(theElasticModel);
pmanager->AddDiscreteProcess(theElasticProcess);
      
G4ProtonInelasticProcess* theInelasticProcess = new G4ProtonInelasticProcess();
G4LEProtonInelastic* theLEInelasticModel = new G4LEProtonInelastic;
theInelasticProcess->RegisterMe(theLEInelasticModel);
pmanager->AddDiscreteProcess(theInelasticProcess);

}

I am using this file proton.mac to describe the source:

/gps/particle proton
/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/pos/radius 60. cm
/gps/pos/centre 0.0 0.0 0.0 cm
/gps/ang/type cos
/gps/ang/mintheta 0. deg
/gps/ang/maxtheta 90. deg
/gps/ang/minphi 0. deg
/gps/ang/maxphi 360. deg
/gps/ene/type Pow
/gps/ene/min 300. MeV
/gps/ene/max 6600. MeV
/gps/ene/alpha -1.6
/run/beamOn 26602349

The compilation is successful, but when I run the simulation with the proton.mac file, an error appears:

Segmentation fault (core dumped)

Any idea?

Thanks in advance,

Manuel

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Idea: Re: Proton physic list between 300 - 6600 MeV   (Vladimir Ivanchenko - 29 Aug, 2011)
 Add Message Add Message
to: "Proton physic list between 300 - 6600 MeV"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]