Forum: Space Applications Not Logged In (login)
Show subscribers

Space Applications

The email gateway for this forum is: space_app-g4hn@slac.stanford.edu


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

Question Simulating Co60 source environment in GRAS  by Leonardo Ghizoni <Leonardo Ghizoni>,   Mar 10, 03:20
Hello,

In the context of Total Ionizing Dose, I wonder if it is possible in GRAS to simulate a Cobalt 60 source, so I can better understand the amount of needed shielding for a given component in space.

Is it just by declaring the Ion Cobalt 60 and decay physics lists? Then I would wish to have the incident 1.17 MeV and 1.33 MeV gamma-rays.

I know I can just shoot the gamma-rays, but I wonder if nuclear decay with the products is possible in GRAS.

Thanks a lot in advance.

/Leonardo.

None mesh scoring in GRAS  by Svetlana Shasharina <Svetlana Shasharina>,   Dec 06, 14:47
Hi, Is it possible to have mesh scoring in GRAS? Thank you, Sveta

Question Extreme Ultraviolet simulation  by Gianluca De Zanet <Gianluca De Zanet>,   27 Sep, 2017
Hi everyone, I'm a newbie with Geant4 and the environments it could simulate, so better asking directly: is it possible to simulate the UV environment or its energy is too low for Geant4? Thank you,

Gianluca

1 Note: Re: Extreme Ultraviolet simulation   (Tim Meehan - 28 Sep, 2017)
2 None: Re: Extreme Ultraviolet simulation   (php1ic - 29 Sep, 2017)
Question non-spherical sources for estimating geometric factor  by Elias Roussos <Elias Roussos>,   20 Sep, 2017
Hi,

I am trying to estimate the geometry factor of a particle telescope. &#913; macro file used for the input particles is:

/gps/particle e- /gps/pos/type Surface /gps/pos/shape Sphere /gps/centre 0 0 0 cm /gps/pos/radius 20 cm /gps/ang/type cos /gps/ene/mono 50. keV

If I use this, I can calculate the GF of the sensitive detector as:

GF = 4*PI^2*(Radius)^2*(N_detected/N_source)

and I get a result close to the expected value.

When I make a run with low energy particles that cannot penetrate the instrument shielding, it is enough to inject the paticles only from a surface in front of the collimator/telescope entrance, just to improve statistics. I defined a spherical cap with an opening e.g. of 60 deg like this:

/gps/hist/type biaspt

/gps/hist/point 0.0 0.0

/gps/hist/point 0.333 1.0

/gps/hist/point 1.0 0.0

The area of the spherical cap is: A=2*!PI*(1-cos(60))

How is the GF estimated then? unless I am doing something wron with the positioning biasing, I would have expected "N_detected" to be reduced by approx. A/(4*PI*R^2), but results do not seem to add up.

Thanks a lot in advance.

1 Note: Re: non-spherical sources for estimating geometric factor   (Elias Roussos - 20 Sep, 2017)
Question how can i use geant4 outcome to draw the dose-depth figure of mars soil.  by zhiyi han <zhiyi han>,   17 Sep, 2017
how I draw dose-depth profile in soil:

1. fitting the GCR spectra from 2Gev-100Gev, i get y1=x*10^((0.18*(log10(x))^3)-(2.6*(log10(x))^2)+9.8*(log10(x))-11);y1(particles/m2-s-sr-MeV/nuc) ;   x(MeV/nuc)


2,cut GCR spectra into 3 sections: 1-10Gev, 10-100Gev, 100-1000Gev.  
Take samples (for every sample, beam on 1000) from each section and run the simulation and get the dose data.  
(total soil depth is 2m, score mesh is 40,so dose is recorded every 5cm in soil)


3,in first soil layer(0-5cm),plot the Dose-Energy image for three sections. Fitting them. y(Gy); x(Mev),i get

1-10Gev		y2=(-(3.5*10^-32)*x^4+(7.8*10^-28)*x^3-(5.5*10^-24)*x^2+(1.6*10^-20)*x-1.5*10^-17);

10-100Gev 		y3=(-(3.1*10^-40)*x^5+(7.2*10^-35)*x^4-(5.5*10^-30)*x^3+(1.4*10^-25)*x^2+(2.3*10^-21)*x-1*10^-17);

100-1000Gev   y4=((1.8*10^-33)*x^3-(3.4*10^-27)*x^2+(1.8*10^-21)*x+3.6*10^-17);


4,dose that is produced by one particle * number of the particle = total dose.

base on my calculation, the number of particle is y1/x , dose is y2, y3, y4 for each section

so F1=(y1/x)*y2,integrate F1 from 1-10Gev,  
F2=(y1/x)*y3,  integrate F2 from 10-100Gev,
F3=(y1/x)*y4,  integrate F3 from 100-1000Gev.  sum them up, i get F4, F4=F3+F2+F1.

score box is 20m * 20m *2 m, so F4*400*(60*60*24*365)=F5, F5 is the total dose that is produced by all GCRs in this soil layer during the whole year.


but the outcome i get is about 10e-7 smaller than the reference. what is wrong? any help? THANKS for your attention and help in advance.



*****************************************************************

details:


Geant4 simulation of  radiation in Mars atmosphere and soil

version
GEANT4.10.02.p02

application
Based on B4a, changed it’s geometry and physics_list, added scorerwriter in B4a Main().     (.../geant4.10.02.p02/examples/basic/B4/B4a)

physics_list  
QGSP_BERT_HP 3.0

geometry 
world       			30*m X 30*m X 140000*m  (X Y  Z)
20 layers of air.(upon soil)   		20*m X 20*m X 70000*m  (X Y  Z)
Soil layer(under the air layers) 			20*m X 20*m X 20*m (X Y  Z)

primary generator(particle gun)
particle:proton
position(0, 0, -70000*m) (X Y  Z)
direction(0, 0, 1*m)(X Y  Z)


dose, energy deposition, secondaries data collected from run1.mac boxMesh. Scorer size is 20*m X 20*m X 2*m (X Y  Z), every layers is 5*cm along Z axies.

/score/create/boxMesh boxMesh_1
/score/mesh/boxSize 10 10 1 m
/score/mesh/nBin 1 1 40
/score/dumpQuantityToFile boxMesh_1 edep scoreenergydeposit.txt
/score/dumpQuantityToFile boxMesh_1 ndar scorenumberofsecondaries.txt
/score/dumpQuantityToFile boxMesh_1 dpo  scoredosedeposit.txt
None simplest example of reverse monte carlo for GRAS  by Svetlana Shasharina <Svetlana Shasharina>,   25 Aug, 2017
Hello, I have a setup using gdml and gps file. How do I modify it for an adjoint simulation in GRAS? I found the simple_geo example in GRAS, but was hoping to use the gps file in the adjoint simulation, if this is possible. Best regards, Sveta

1 None: Re: simplest example of reverse monte carlo for GRAS   (M.Tevfik Kaplanoglu - 25 Aug, 2017)
(_ None: Re: simplest example of reverse monte carlo for GRAS   (Svetlana Shasharina - 28 Aug, 2017)
Question capture high precision neutron  by jacques <jacques>,   23 Aug, 2017
hello, im working on measuring nbr of hits collected by few crystals located in a box, till now im able to record the total nbr of hits including the scattered ones...but i need to measure just the captured neutrons by each crystal. any suggestion how to add it? thanks in advance

Question a question about gps/pos for gaussian beam geometry  by <Claire>,   10 May, 2017
Hello I have a question about : gps/pos/ I would like to model a monochromatic beam with a Gaussian geometry (Gaussian symmetric shape in a XY plane perpendicular to the direction of propagation Z). The FWHM is the same in both directions X and Y, of 1 micrometer. So sigma = 1/2.355 =0.4246 micrometers. The beam is purely Gaussian, there is no central circular plateau.

Example 31 of the Gant4 toolkit is very close to what I would like to do:

# the beam spot is centered at the origin and is

# of 1d gaussian shape with a 1 mm central plateau

/gps/pos/shape Circle

/gps/pos/centre 0. 0. 0. mm

# Should I comment this line? /gps/pos/radius 1. mm

/gps/pos/sigma_r .4246 um

So my question is: Should I comment the line with radius, since our beam does not have a central plateau? Or should I take a radius with a very small value compared to sigma_r ? (otherwise we would define a circle shape without radius, which would maybe be not possible ?)Or is there another shape than "circle" for a pure Gaussian shape ? Thank you

Question GRAS returns error when using Heavy Ions primaries on composite shieldings  by Leonardo Ghizoni <Leonardo Ghizoni>,   30 Apr, 2017
Hello,

I'm simulating different mixtures for shielding materials, and I've got the following error when shooting Au197 ions:

-------- EEEE ------- G4Exception-START -------- EEEE -------
*** G4Exception : had006
      issued by : G4HadronicProcess::PostStepDoIt
In /home/lghizoni/Geant4/geant4.10.01/source/processes/hadronic/util/src/G4HadSecondary.cc, line 36: 
===> ATTEMPTING TO CREATE A SECONDARY WITH NEGATIVE KINETIC ENERGY.
Call for Binary Light Ion Cascade
Target element Re  Z= 75  A= 187
Unrecoverable error in the method ApplyYourself of ionInelastic
TrackID= 1  ParentID= 0  Au197
Ekin(GeV)= 10.2981;  direction= (0.00233913,-0.00273888,-0.999994)
Position(mm)= (3.05801,-13.2243,11.7177);  material NanoShield4
PhysicalVolume  <Shield4>
 ApplyYourself failed

*** Fatal Exception *** core dump *** -------- EEEE -------- G4Exception-END --------- EEEE -------

*** G4Exception: Aborting execution *** Aborted (core dumped)

I'm using GRAS 3.4 with GEANT4 10.01 patch3. The shielding material contains the following elements: Cr, Fe, Re and B.

Do anyone know know how to solve this?

Thanks in advance.

Leonardo.

1 Question: Re: GRAS returns error when using Heavy Ions primaries on composite shieldings   (Vladimir Ivanchenko - 01 May, 2017)
(_ Ok: Re: GRAS returns error when using Heavy Ions primaries on composite shieldings   (Leonardo Ghizoni - 01 May, 2017)
Question Understanding GRAS RMC Results  by Leonardo Ghizoni <Leonardo Ghizoni>,   20 Sep, 2016
Hi experts,

I'm having simulations on GRAS for a space system and I'm applying the RMC method. I'm a little confused about the results. I ask the simulation to calculate the total on a given device of my system:

/gras/analysis/dose/addModule TID_{TargetName}

/gras/analysis/dose/TID_{TargetName}/bookTuples false

/gras/analysis/dose/TID_{TargetName}/addVolume *{TargetName}*

/gras/analysis/dose/TID_{TargetName}/setUnit rad

/gras/analysis/dose/TID_{TargetName}/initialise

Then I get two different histograms:

TID Dose TID total dose vs primary kine

Summing the results of the second histogram, I should get the same value of dose as the first one, right? But they are very different. At the same time, the dose value of the first histogram gives extremely high numbers. Am I missing something here?

Thanks in advance.

None Different ionic states of elements for materials  by <MWestphal>,   02 Sep, 2016
Hello,

Is there a way to create materials out of elements in different ionic states (i.e. Fe +2, Fe +3, etc.) in Geant4?

Thanks in advance.

Question How to use and normalise an energy spectrum as input  by Germán Ros <Germán Ros>,   27 Jul, 2016
Hello!

I am trying to use an energy spectrum as an input and I have seen the example called "radioprotection" in the G4 advanced examples where this is done using the configuration file "primary.mac". Unfortunately I have several problems and things I do not really understand and I have not found the answer in Geant4 forum or elsewhere.

I would be very grateful if someone could help me.

I have the differential and integral spectrum of galactic cosmic rays obtained from SPENVIS. I have found an excel file that is used by several people in this forum to create the G4 configuration file given an energy spectrum, exactly what I need.

http://geant4.esa.int/tl_files/geant4/resources/CREME_M1_GCR_Flux_Spectra.xlt

Using it, I have created my configuration file that it is attached (GCRs_Carbon.mac). But I have several problems:

1) Following the .xls file, in the .mac file I have defined two histograms, one with the differential spectrum and a second one with the probabilities and weights for each energy. However in primary.mac only the first histo is defined, so ¿is the second one really needed?

2) Again, following the .xls file, I have calculated a normalization factor (at the end of GCRs_Carbon.mac). However, again in primary.mac a normalization factor is not used. I think it is needed since I have to simulate several nuclei and compare the different components, but ¿am I correct or this factor is not needed?

3) Why in the .xls file the normalization factor has /cm2/s as units but when translated to G4 code is set as cm2?

4) When trying to run geant4 with my GCRs_Carbon.mac the line /analysis/normalise 7.305E+08 give me a core dump with no other information. ¿Do you know what am I doing wrong?

5) A more general question. When you use a spectrum as input, how are particles simulated in time? I mean, all of them are thrown at the same time, uniformly in a time interval that should be defined (how?) or something else?

Thank you very much in advanced

Best Regards,

Germán Ros

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/07/27/09.02-52585-GCRs_Carbon.mac

1 None: Re: How to use and normalise an energy spectrum as input   (Sveta - 28 Jul, 2016)
2 None: Re: How to use and normalise an energy spectrum as input   (Germán Ros - 07 Sep, 2016)
1 Question: Re: How to use and normalise an energy spectrum as input   (panda - Jan 22, 17:56)
Question Isotropic flux with Planetocosmics: GPS comands not working properly   by Michael Milde <Michael Milde>,   12 Jul, 2016
Hello everyone,

I need to generate an isotropic particle flux on Earth´s surface through the Planetocosmics framework. For the source geometry I use the following Geant4 /gps commands:

/gps/pos/type Surface

/gps/pos/shape Sphere

/gps/pos/centre 0. 0. 0. km

/gps/pos/radius 7870 km

/gps/ang/type cos

/gps/ang/mintheta 0.00 deg

/gps/ang/maxtheta 90.00 deg

The expected flux from the starting positions on the source sphere, as I´ve read in many threads here, should have a cos modulation centred around the perpendicular direction at the source sphere (and thus directed to the Earth's centre, since my source sphere is centred in 0 0 0 in a planetocentric coordinate frame).

Now my problem: if I run the simulation, the particles are correctly generated from the source sphere centred around the Earth, but their direction is centred towards a point somewhere in space, and not at the Earth's centre, as I would like. In the attachments you can find pictures (1-2) of the tracks. In picture 3 and 4 you can see the result if I set /gps/ang/mintheta 0.00 deg and /gps/ang/maxtheta 0.00 deg, so just the (expected) "zenith" direction is allowed. Here, all the tracks are pointed towards some point in space.

I already tried to add commands like /gps/ang/surfnorm, /gps/ang/user_coor, but nothing seems to help. Also the circle source has the same behaviour.

Did someone encounter this or a similar problem already and was able to solve it? Do you have some advices?

Thanks! Best Regards,

Michael

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/07/12/03.34-42312-PCSimulationSphere1.png
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/07/12/03.34-7371-PCSimulationSphere2.png
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/07/12/03.34-56069-PCSimulationSphere3.png
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/07/12/03.34-7442-PCSimulationSphere4.png

1 None: Re: Isotropic flux with Planetocosmics: GPS comands not working properly   (sveta - 12 Jul, 2016)
(_ Question: Re: Isotropic flux with Planetocosmics: GPS comands not working properly   (Michael Milde - 12 Jul, 2016)
(_ None: Re: Isotropic flux with Planetocosmics: GPS comands not working properly   (sveta - 12 Jul, 2016)
Question GRAS license  by Svetlana Shasharina <Svetlana Shasharina>,   23 May, 2016
I have been using GRAS for a while now and have been very happy with this application. Some files of GRAS refer to G4 license, but some state that the implementation is the intellectual property of ESA. What do these statements mean for the GRAS be distribution? Best regards, Sveta

1 None: Re: GRAS license   (Vladimir Ivanchenko - 25 Jul, 2016)
(_ None: Re: GRAS license   (sveta - 25 Jul, 2016)
(_ None: Re: GRAS license   (Vladimir Ivanchenko - 25 Jul, 2016)
Question Getting SSAT tool  by Leonardo Ghizoni <Leonardo Ghizoni>,   27 Mar, 2016
Hello,

I'm trying to get a local copy of SSAT tool, but I'm unable to find the proper place to get it.

I've been using the tool from SPENVIS, but I know it has some limitations regarding the processing and geometry.

Thank you in advance. Leonardo.

Question Unexpected results for fluence and TID using Spenvis/GRAS  by Svetlana Shasharina <Svetlana Shasharina>,   08 Mar, 2016
I did the following calculations using Spenvis/GRAS(gdml):

1. Independent runs for two vacuum spheres with r = 19cm and r = 15 cm. Source: Monoenergetic (100MeV) electrons with fluence 100 (1/cm^2 s), sphere with r = 20cm.

SPENVIS calculated fluences 28390 and 17644 for the first and the second spheres. The ratio of fluences/surface_area is constant for these cases (6.24 and 6.26) and is very close to 2pi :-) but is very different from the incident fluence (100).

First confusion: I thought the calculated fluence is normalized to the target's areas. But looks like it is not. Second confusion: after dividing by the surface areas I still have a number very different to the incident environment number (100).

Maybe I do not understand what the incident fluence really mean? How do I get a number that would relate to it?

2. Ta sphere r = 15cm, calculating TID in MeV with the source described. I get 1.6564e+06MeV/s.

If all e were to absorb, I would expect 100(1/cm^2 s)*100MeV*SurfaceArea = 1e+4 MeV/s *2826cm^2 = 2.8e+7 MeV/s.

So about 1/16 of the electrons gets absorbed, I guess? If I change normalization to CURRENT instead of FLUX, I will get 1/4 difference, I assume. Still a big difference. Any intuitive reaction to this difference?

Best regards, Sveta

1 None: Re: Unexpected results for fluence and TID using Spenvis/GRAS   (Svetlana G. Shasharina - 08 Mar, 2016)
1 None: Re: Unexpected results for fluence and TID using Spenvis/GRAS   (Svetlana Shasharina - 10 Mar, 2016)
2 More: Re: Unexpected results for fluence and TID using Spenvis/GRAS   (Roger Byrd - 02 Apr, 2016)
1 None: Re: Unexpected results for fluence and TID using Spenvis/GRAS   (sveta - 09 Apr, 2016)
None Changing initial seed in GRAS  by Leonardo Ghizoni <Leonardo Ghizoni>,   07 Mar, 2016
Hello,

I wander what is the best way of changing the initial seed in GRAS.

I see that I could add in the main file (assuming I'm using GRAS as an external application):

#include "Randomize.hh"
#include "time.h"

CLHEP::HepRandom::setTheEngine(new CLHEP::RanecuEngine());

G4long seed = time(NULL);

CLHEP::HepRandom::setTheSeed(seed);

But I would like to know if there's a way of changing this via macro, as that's the GRAS "way".

My initial seed is always the same:

Initial seed = 9876
 u[] = 0.174685 0.883674 0.958512 0.690013 0.0250537 0.142294 0.783649 0.356081

0.214774 0.69287 0.700302 0.465109 0.22215 0.768752 0.461089 0.397602 0.522215

0.121546 0.960227 0.510892 0.0116568 0.0286447 0.238777 0.585219 0.780612

0.709182 0.082106 0.828601 0.61446 0.0798209 0.688002 0.00978529 0.370831

0.376435 0.710184 0.864017 0.167338 0.21114 0.349235 0.522479 0.943713

0.0702708 0.778742 0.639458 0.506555 0.151751 0.990123 0.564688 0.430497

0.889832 0.995398 0.725368 0.710512 0.0550886 0.0540335 0.0159176 0.728379

0.947141 0.34922 0.20551 0.351029 0.672539 0.118976 0.655784 0.195375 0.0380308

0.293556 0.237497 0.79723 0.158375 0.029761 0.602215 0.616657 0.995638 0.734827

0.943425 0.130007 0.277064 0.306344 0.973337 0.183926 0.547536 0.914344

0.666572 0.210074 0.300351 0.83602 0.267525 0.248549 0.533332 0.1754 0.299529

0.538308 0.555472 0.283894 0.485254 0.222273

c = 0.0216029, cd = 0.456233, cm = 1

i97 = 96, u[i97] = 0.222273

j97 = 32, u[j97] = 0.370831

As a consequence, I always get exactly the same results.

Thank a lot. Leonardo.

1 None: Re: Changing initial seed in GRAS   (Hugh Evans - 07 Mar, 2016)
(_ None: Re: Changing initial seed in GRAS   (Leonardo Ghizoni - 07 Mar, 2016)
Question Non-Ionizing energy loss (NIEL)  by aimsphere <aimsphere>,   05 Feb, 2016
Hello G4 developers and users,

I would like to estimate displacement damage (DD) in a silicon detector in space. I understand that it is calculated as a product of the Non-Ionizing energy loss (NIEL) and the particle fluence. My question here is about NIEL calculation.

I found that there is a method called GetNonIonizingEnergyDeposit() which can be accessed from the G4step class. Is this method appropriate to estimate NIEL and thus the DD dose? Has anyone tried it?

Thank you.

1 None: Re: Non-Ionizing energy loss (NIEL)   (Hugh Evans - 11 Feb, 2016)
(_ None: Re: Non-Ionizing energy loss (NIEL)   (xue - 28 Mar, 2016)
None tessellated in gdml  by andrea lucaroni <andrea lucaroni>,   23 Dec, 2015
Dear all

i have a geometry in gdml format.
I created this geometry with Instep, the code is loaded correctly, but it
seems that the volumes are empty, when i run with gras, i'm reporting under
the line of analysis of geometry by gras.

[GRASGeometryUtil::ListLogicalVolumes]: 471 volumes in the store
   Volume <CSP_1_vol>   Material: Al   Mass: [kg] 1.62854e-12   ([kg]
1.62854e-12  including sub-tree)
   Volume <CSP_2_vol>   Material: Si   Mass: [kg] 8.45983e-13   ([kg]
8.45983e-13  including sub-tree)
   Volume <CSP_3_vol>   Material: Al   Mass: [kg] 1.16744e-14   ([kg]
1.16744e-14  including sub-tree)
   Volume <CSP_4_vol>   Material: Al   Mass: [kg] 1.1665e-14   ([kg]
1.1665e-14  including sub-tree)
   Volume <CSP_5_vol>   Material: Al   Mass: [kg] 5.7097e-14   ([kg]
5.7097e-14  including sub-tree)
   Volume <CSP_6_vol>   Material: Si   Mass: [kg] 5.25532e-13   ([kg]
5.25532e-13  including sub-tree)
   Volume <CSP_7_vol>   Material: Al   Mass: [kg] 4.4183e-14   ([kg]
4.4183e-14  including sub-tree)


for example:
mass of CSP_2_vol must be about 4e-3 kg not   8.45983e-13 kg

Are there commands from put in the gdml file to consider the solid volumes ?

thanks
best regardes
Andrea


Andrea

 [ MIME part of type text/html without a name stripped ]

1 None: Re: tessellated in gdml   (sveta - 23 Dec, 2015)
Question Dose-depth curves from SSAT  by Michael Harig <Michael Harig>,   15 Dec, 2015
My apologies from the beginning. I am a new user to SPENVIS/GEANT4.

I currently have Geant4 working without any problems.

I am trying to evaluate radiation transport through various shielding materials in space environments such as on the ISS and/or Mars using SSAT tool in SPENVIS.

Will someone please help point me in the right direction so that I can get the dose-depth curves plotted?

1 None: Re: Dose-depth curves from SSAT   (Hugh Evans - 16 Dec, 2015)
None Fwd: install gras  by andrea lucaroni <andrea lucaroni>,   12 Dec, 2015
Hi,

I have just installed geant4.9.6.p04 and GRAS-03-03 on my linux computer. I
have enabled GDML, i try to run the exemple:

automatic_normalisation/

gras perevent_normalisation.g4mac


I get this error:

4GDML: Reading 'geo1D.gdml' done!

Stripping off GDML names of materials, solids and volumes ...

[PhysicsList::ConstructProcess]

unusable number for repFlag: repFlag=1071881836

In
/home/dicrescenzo/ARTE/soft/geant4.9.6.p03/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPElasticFS.cc,
line 178:

===> G4NeutronHPElasticFS::Init -- unusable number for repFlag

terminate called after throwing an instance of 'G4HadronicException'

  what():  std::exception

Aborted


Thank you in advance for any help,

Andrea

 [ MIME part of type text/html without a name stripped ]

1 Idea: Re: Fwd: install gras   (Vladimir Ivanchenko - 13 Dec, 2015)
(_ None: Re: Fwd: install gras   (andrea lucaroni - 30 Dec, 2015)
(_ None: Re: Fwd: install gras   (andrea lucaroni - 30 Dec, 2015)
(_ Idea: Re: Fwd: install gras   (Vladimir Ivanchenko - 04 Jan, 2016)
(_ None: Re: Fwd: install gras   (andrea lucaroni - 04 Jan, 2016)
(_ Idea: Re: Fwd: install gras   (Vladimir Ivanchenko - 05 Jan, 2016)
Question GRAS crashes for Z >= 2  by Leonardo Ghizoni <Leonardo Ghizoni>,   12 Nov, 2015
Hello,

I'm using GRAS 3.3 as external application with GEANT4.9.6.p04. I'm building a satellite for LEO mission and taking the macro files from SPENVIS as input. I use QBBC as my physics list.

The simulation runs well for trapped electrons, trapped protons and solar protons. When running with GCR Z = 1 source, I start to get some warning messages:

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : had012
      issued by : G4HadronicProcess:CheckResult()
Warning: Bad energy non-conservation detected, will re-sample the interaction
 Process / Model: ionInelastic / Binary Light Ion Cascade
 Primary: Li7[0.0] (1000030070), E= 15557.6, target nucleus (13,27)
 E(initial - final) = 6668.06 MeV.

*** This is just a warning message. *** -------- WWWW -------- G4Exception-END --------- WWWW -------

But the simulation ends with no further problems. When I run with GCR Z >= 2, the warning messages become very frequent, till the simulation crashes:

In /home/lghizoni/Geant4/geant4.9.6/source/processes/hadronic/models/parton_string/management/src/G4VPartonStringModel.cc, line 81: ===> G4VPartonStringModel::Scatter(): fails to generate strings

-------- EEEE ------- G4Exception-START -------- EEEE -------
*** G4Exception : had006
      issued by : G4HadronicProcess::PostStepDoIt
In /home/lghizoni/Geant4/geant4.9.6/source/processes/hadronic/models/parton_string/management/src/G4VPartonStringModel.cc, line 81:
===> G4VPartonStringModel::Scatter(): fails to generate strings
Call for FTFP
Target element H  Z= 1  A= 1
Unrecoverable error in the method ApplyYourself of ionInelastic
TrackID= 1  ParentID= 0  Li7[0.0]
Ekin(GeV)= 50.4311;  direction= (-0.30308,-0.102135,0.947476)
Position(mm)= (-59.0097,-284.055,-1586.32);  material Hydrazine
PhysicalVolume  <Fuel>
 ApplyYourself failed

*** Fatal Exception *** core dump *** -------- EEEE -------- G4Exception-END --------- EEEE -------

*** G4Exception: Aborting execution *** Aborted

I tried to change the physics list to FTFP_BERT_HP and Shielding, but the simulation still crashes.

I thank very much for help, since I cannot find a solution for the problem.

Leonardo.

1 None: Re: GRAS crashes for Z >= 2   (Hugh Evans - 13 Nov, 2015)
(_ Idea: Re: GRAS crashes for Z >= 2   (Vladimir Ivanchenko - 13 Nov, 2015)
(_ Question: Re: GRAS crashes for Z >= 2   (Leonardo Ghizoni - 13 Nov, 2015)
1 Idea: Re: GRAS crashes for Z >= 2   (Vladimir Ivanchenko - 14 Nov, 2015)
1 Sad: Re: GRAS crashes for Z >= 2   (Leonardo Ghizoni - 14 Nov, 2015)
... 4 Message(s)
2 News: Re: GRAS crashes for Z >= 2   (Hugh Evans - 16 Nov, 2015)
1 Ok: Re: GRAS crashes for Z >= 2   (Leonardo Ghizoni - 16 Nov, 2015)
Question Simulation of anisotropic particle fluxes  by Michael Milde <Michael Milde>,   23 Sep, 2015
Hello,

I want to simulate the response of my detector in the anisotropic flux of trapped protons in the South Atlantic Anomaly. The angular distribution for the flux is calculated by the SPENVIS webtool in form of a histogram. At the moment I am not sure how to build the particle gun in a way that I have the same flux at every position in the world volume.

Maybe someone has experience in creating such a gun. Thanks in advance.

Greetings,

Michael

1 None: Re: Simulation of anisotropic particle fluxes   (Hugh Evans - 23 Sep, 2015)
(_ None: Re: Simulation of anisotropic particle fluxes   (Hugh Evans - 24 Sep, 2015)
(_ Ok: Re: Simulation of anisotropic particle fluxes   (Michael Milde - 01 Oct, 2015)
Question how to find the cross section data of the nuclear reactions  by liyong <liyong>,   16 Mar, 2015
Dear all,

I am trying to simulate the generation of cosmogenic nuclei 10Be and 26Al in lunar soil. I didn't find the cross section data of the reaction 'O(n,x)10Be' in G4NDL4.5/IsotopeProduction/CrossSection, however, 10Be is produced in the simulation of 100-200MeV neutron interact with the Oxygen. Could you tell me which cross section data is used in this reaction. And could you tell how can I find the cross section data of the nuclear reactions, such as 'O(p,x)10Be'.

Thank you in advance. Liyong

1 None: Re: how to find the cross section data of the nuclear reactions   (Dennis H. Wright - 16 Mar, 2015)
Question Initial run of GRAS-03-03 with geant4.9.6.p04  by Richard Gillham Darnley <Richard Gillham Darnley>,   04 Mar, 2015
Hi,

I have just installed geant4.9.6.p04 and GRAS-03-03 on my linux computer. I have enabled GDML options and ran example G01 to test my basic GDML and XercesC install. This example worked fine.

When running the configure_and_install.sh script for GRAS-03-03 I had to add the following cmake options so that it could find XercesC:

-DXERCESC_LIBRARY=/home/richard/xercesc/xerces-c-src_2_8_0/lib/libxerces-c.so -DXERCESC_INCLUDE_DIR=/home/richard/xercesc/xerces-c-src_2_8_0/include

Now when I attempt to run gras by typing 'gras' in the terminal I get the following message:

gras: error while loading shared libraries: libxerces-c.so.28: cannot open shared object file: No such file or directory

This file does exist (its in the same location as libxerces-c.so). How should I tackle this problem? Should I change the library cmake option to look for 'libxerces-c.so.28' instead or is there another step in the process that I am missing?

Thank you in advance for any help,

Richard

1 None: Re: Initial run of GRAS-03-03 with geant4.9.6.p04   (sveta - 04 Mar, 2015)
1 Feedback: Re: Initial run of GRAS-03-03 with geant4.9.6.p04   (Richard Gillham Darnley - 04 Mar, 2015)
2 Idea: Re: Initial run of GRAS-03-03 with geant4.9.6.p04   (Michael H. Kelsey - 04 Mar, 2015)
1 Feedback: Re: Initial run of GRAS-03-03 with geant4.9.6.p04   (Richard Gillham Darnley - 04 Mar, 2015)
None I need to study the depth dose-distributions and partial fragmentation cross sections or yields of charged nuclear fragm  by Wu Zheng <Wu Zheng>,   07 Nov, 2014
For example ,I want to Simulated fragmentation distribution as a function of the fragment Z after 1 g/cm of Al hit by 1 GeV/nucleon Fe ions by geant4. Who can give me some suggestions?

Question Dose Normalization  by Leonardo Ghizoni <Leonardo Ghizoni>,   06 Oct, 2014
Hello Experts,

I'm having some difficulties when trying to normalize my outputs, like the Total Dose.

First I'm running with a simple geometry: a spherical target sorrounded by a spherical shell with 2mm thickness.

To compute the dose, I'm using the primitive scorer G4PSDoseDeposit and using GPS as the source.

I took a look at a ppt document of G. Santin about normalization, and as I understood, the "real world" Dose is calculated as:

Dr = Nr*Ds/Ns,

where Nr is the number of particles in the real world, Ds is the dose calculated in the simulation and Ns is the number of events simulated, right?

Now, the "real world" number of particles Nr is calculated by the flux. That's where I'm a little confused. I do have the fluxes as functions of the energy, in a table calculated by Spenvis and the respective macro created. As I understad, the number Nr is calculated by:

Nr = flux * 4*pi^2*r^2.

I don't know which flux to use, since I have a table based on energies. And the radius is defined by me on the source, right? I defined the source like this:

/gps/particle proton

/gps/ene/type Arb

/gps/hist/type arb

/gps/hist/point 1.000000E-01 2.407800E+01

/gps/hist/point 1.500000E-01 2.397600E+01

/gps/hist/point 2.000000E-01 1.974300E+01

/gps/hist/point 3.000000E-01 1.135400E+01

/gps/hist/point 4.000000E-01 1.059100E+01

/gps/hist/point 5.000000E-01 9.840100E+00

/gps/hist/point 6.000000E-01 1.154300E+01

/gps/hist/point 7.000000E-01 1.277300E+01

/gps/hist/point 1.000000E+00 8.546000E+00

/gps/hist/point 1.500000E+00 3.924800E+00

/gps/hist/point 2.000000E+00 3.838800E+00

/gps/hist/point 3.000000E+00 3.667500E+00

/gps/hist/point 4.000000E+00 3.831800E+00

/gps/hist/point 5.000000E+00 3.993600E+00

/gps/hist/point 6.000000E+00 4.341400E+00

/gps/hist/point 7.000000E+00 4.587900E+00

/gps/hist/point 1.000000E+01 3.820900E+00

/gps/hist/point 1.500000E+01 2.932800E+00

/gps/hist/point 2.000000E+01 2.538000E+00

/gps/hist/point 3.000000E+01 2.066400E+00

/gps/hist/point 4.000000E+01 1.757500E+00

/gps/hist/point 5.000000E+01 1.592800E+00

/gps/hist/point 6.000000E+01 1.462600E+00

/gps/hist/point 7.000000E+01 1.370200E+00

/gps/hist/point 1.000000E+02 1.118000E+00

/gps/hist/point 1.500000E+02 7.635400E-01

/gps/hist/point 2.000000E+02 4.962900E-01

/gps/hist/point 3.000000E+02 2.196200E-01

/gps/hist/point 4.000000E+02 2.586700E-02

/gps/hist/inter Spline

/gps/pos/type Surface

/gps/pos/shape Sphere

/gps/pos/centre 0. 0. 0. mm

/gps/pos/radius 9 mm

/gps/ang/type cos

/gps/ang/mintheta 0.000E+00 deg

/gps/ang/maxtheta 9.000E+01 deg

/gps/source/list

So, my radius will be 9mm.

Finally, my normalized dose changes values when I change the number of primary events...this sounds incorret, since I'm normalizing it.

I know I'm doing it wrong, just don't know where.

Please help? I'm really having a hard time on this.

Thank you in advance.

Best Regards, Leonardo.

Question modeling space environment in testing facilities  by Svetlana Shasharina <Svetlana Shasharina>,   30 Jul, 2014
Sorry, if I am posting the second time (do not see my first post).

I am given integral fluences for electrons and protons (Europa environment). The testing facility will give us particular monoenergetic fluxes in microsec pulses in particular energies. What is the accepted practice to reproduce the TID effects in such facilities?

Thank you, Sveta

Question calculating secondary electrons from electron source (using GRAS)  by Svetlana Shasharina <Svetlana Shasharina>,   30 Jul, 2014
Hi everybody,

I need to calculate the flux of "extra" electrons produced by enviroment electrons in Europa environment. I have beed using GRAS to do my simulations and can calculate fluxes "to", "from" and "between" volumes but I cannot directly calculate the flux of the new electrons that might appear due to e radiation.

Questions: 1. Can I ignore the secondary electrons from electron sources and consider electrons produced by protons dominating (I know how to calculate secondary electrons when I have ion source:-)?

2. Can I estimate the secondary flux by subtracting "in" flux from the "from" flux? I understand that the logic is not flawless as some incident electrons might not leave the volume.

Thank you, Sveta

Question How to compute Dose with GPS?  by Leonardo Ghizoni <Leonardo Ghizoni>,   22 Jul, 2014
Hello,

I'm using the GPS to simulate the space radiation environment...

What would be the right angular distribution to use, assuming I want an omnidirectional distribution?

Other question: how do I compute the dose on a specified part of my geometry usin the GPS? And finally, is it possible to compute the dose on more than one object on the same run?

Thank you in advance.

Best regards. Leonardo

1 None: Re: How to compute Dose with GPS?   (Hugh Evans - 23 Jul, 2014)
(_ Question: Re: How to compute Dose with GPS?   (Leonardo Ghizoni - 23 Jul, 2014)
1 None: Re: How to compute Dose with GPS?   (sveta - 23 Jul, 2014)
2 None: Re: How to compute Dose with GPS?   (Hugh Evans - 23 Jul, 2014)
1 Question: Re: How to compute Dose with GPS?   (Leonardo Ghizoni - 23 Jul, 2014)
... 2 Message(s)
Question Using defualt GRAS histograms instead of AIDA  by Svetlana Shasharina <Svetlana Shasharina>,   10 Jun, 2014
I see the following lines in a g4mac generated by Spenvis:

/gras/histo/fileName spenvis_gras31

/gras/histo/AIDAFileType root

Does this imply that the output uses AIDA? How do I change these lines in order not to depend on AIDA? Should I just ommit the second line?

Thanks, Sveta

1 Idea: Re: Using defualt GRAS histograms instead of AIDA   (Vladimir Ivanchenko - 01 Jul, 2014)
None Model a crossing muon in a ccd cell  by Pierre <Pierre>,   14 May, 2014
Dear Geant4 users,

I am really new in the use of G4 and as its mentionned in the title, I would like to create a box representing my ccd cell and measure the energy deposited by crossing muons which come from isotropic (from a sphere).

I have well written the DetectorConstruction.cc and the PhysicsList.cc files but I dont know exactly how to code the run.cc and primarygenerator.cc files. I attached these files.

Precisely I don't know how to code the incident muons coming from the sphere.

I am a little bit lost in the use og Geant4, everything is still not clear for me. I hope I could get some help for that task.

Cheers

Pierre

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/05/14/04.39-88-DetectorConstruction.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/05/14/04.39-34516-PhysicsList.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/05/14/04.39-46915-PrimaryGeneratorAction.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/05/14/04.39-83636-Run.cc

1 Question: Re: Model a crossing muon in a ccd cell   (Pierre - 16 May, 2014)
(_ None: Re: Model a crossing muon in a ccd cell   (Pierre - 19 May, 2014)
(_ None: Re: Model a crossing muon in a ccd cell   (michel maire - 19 May, 2014)
(_ None: Re: Model a crossing muon in a ccd cell   (John Allison - 19 May, 2014)
Question extra parts add to TID: unexpected result  by Svetlana Shasharina <Svetlana Shasharina>,   05 May, 2014
Hi I am using SPENVIS/GRAS. I calculated fluxes and TIDs in two configurations: (1) cylindrical shield with an opening on top (2) the same shield plus a hollow cone sitting on the opening and having the same size of opening. I used the same fluxes and spherical source in both cases. To my surprise, the TID from electrons of the "detector" inside the cylindrical shields and a proton flux on them was higher in the second case. I then added calculations of fluxes from the cone and it was zero. Flux between the cone and the detector was zero too. I am puzzled by the result. Any ideas? Thanks, Sveta

1 Note: Re: extra parts add to TID: unexpected result   (Svetlana Shasharina - 01 Jul, 2014)
Question local GRAS generating empty histograms  by Svetlana Shasharina <Svetlana Shasharina>,   05 Apr, 2014
I used spenvis and ran TID analysis. I tried to run the same mac and gdml files locally and am getting validation errors:

G4GDML: VALIDATION ERROR! ID value 'v01' is not unique at line: 99
G4GDML: VALIDATION ERROR! ID value 'v02' is not unique at line: 104
G4GDML: VALIDATION ERROR! ID value 'v03' is not unique at line: 109
G4GDML: VALIDATION ERROR! ID value 'v04' is not unique at line: 114
G4GDML: VALIDATION ERROR! ID value 'v05' is not unique at line: 119
G4GDML: VALIDATION ERROR! ID attribute 'G4_Ta' is referenced but was never declared at line: 131
 etc

and getting empty histograms.

Once I follow the suggestion from https://www.spenvis.oma.be/forum/viewtopic.php?f=27&t=228 (although I do not this additional line in any gdml examples) and am still getting the above errors and empty histograms (see attached).

Should I install AIDA?

Thanks, Sveta

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/04/05/06.57-30254-nvis_gras31_e_local.g4mac
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/04/05/06.57-66185-spenvis_suda.gdml
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/04/05/06.57-33506-spenvis_gras31.csv

Question Normalization way for a generating sphere  by Manuel Castro Avila <Manuel Castro Avila>,   20 Mar, 2014
Hello everyone,

I am trying to normalize a counting rate in my simulation. I want to put it in units of counts cm^{-2} s^{-1} kev{-1}

In my simulation I am considering the input spectrum from the SPENVIS system (SPENVIS give me the spectrum in units of particles cm{-2} s^{-1} MeV{-1},
 I am considering the background produced by trapped protons and implemented into GEANT4 via /gps/hist/point).

So, using the input spectrum, I am generating a number N of events from a sphere's surface that hit a circular area of radius r. This circular area represents my detector.

Now, to normalize I proceed on this way:

-From the simulation I get the energy deposited for each event in my detector. Using these values I build a histogram in the energy range I am considering and I get a spectrum in Counts keV^{-1}. Am I right?
-Now, integrating (numerically) the spectrum from SPENVIS in the energy range, I have a flux in particles cm{-2} s{-1} in that energy range and using the geometrical factor
  (for a generating sphere is 2*pi*Acirc, where Acirc is the detector's area), I can calculate what is the integration time corresponding with the number N of events I am considering. Now dividied by
  a integration time I have the spectrum in counts kev{-1} s{-1}.
-Finally, as I know the area of my detector, dividing the previous result by the detector's area I have a spectrum in counts cm{-1} s{-1} kev{-1}.

My question is, am I proceeding on the right way?

Any idea..

Thanks

Manuel

Question CREME_M1_GCR_Flux_Spectra.xlt  by Ivan Brilkov <Ivan Brilkov>,   14 Mar, 2014
I've tried to use that file to to get fluxes in gps format: http://geant4.esa.int/tl_files/geant4/resources/CREME_M1_GCR_Flux_Spectra.xlt But macros seems not working in my version of Excel, in what version it works correctly?

PS: is still available Qinetiq sites about gps, mulassis, and others: http://reat.space.qinetiq.com/mulassis/

1 None: Re: CREME_M1_GCR_Flux_Spectra.xlt   (Hugh Evans - 17 Mar, 2014)
(_ None: Re: CREME_M1_GCR_Flux_Spectra.xlt   (Ivan Brilkov - 26 Sep, 2014)
None The 10th Geant4 Space Users Workshop  by Svetlana Shasharina <Svetlana Shasharina>,   19 Feb, 2014
If this is an open event, could somebody post the URL please? Thanks, Sveta

1 None: Re: The 10th Geant4 Space Users Workshop   (Asai, Makoto - 19 Feb, 2014)
Question differences in TIDs produced by gras and mulassis  by Svetlana Shasharina <Svetlana Shasharina>,   12 Feb, 2014
I was using gras with mulassis type of geometry (hollow sphere with a solid sphere inside) and am getting huge differences in TID (the ratios are not constant, though):

TID Al (gras) TID Cu (gras) 1.07E+15 1.01E+13

TID Al (mulassis) TID Cu (mulassis) 5.70E+06 1.40E+05

The inputs were generated using spenvis portal but I also ran gras 2.1 locally and confirmed that the resulting TIDs are huge.

Anybody can help me to interpret these differences?

Thanks, Sveta

1 None: Re: differences in TIDs produced by gras and mulassis   (Sveta - 12 Feb, 2014)
2 None: Re: differences in TIDs produced by gras and mulassis   (Laurent Desorgher - 24 Feb, 2014)
1 None: Re: differences in TIDs produced by gras and mulassis   (sveta - 25 Feb, 2014)
(_ More: Re: differences in TIDs produced by gras and mulassis   (Roger Byrd - 05 Aug, 2015)
(_ None: Re: differences in TIDs produced by gras and mulassis   (Sveta - 05 Aug, 2015)
(_ More: Re: differences in TIDs produced by gras and mulassis   (Roger Byrd - 09 Aug, 2015)
(_ None: Re: differences in TIDs produced by gras and mulassis   (sveta - 09 Aug, 2015)
(_ None: Re: differences in TIDs produced by gras and mulassis   (John Allison - 10 Aug, 2015)
(_ None: Re: differences in TIDs produced by gras and mulassis   (sveta - 10 Aug, 2015)
(_ None: Re: differences in TIDs produced by gras and mulassis   (Hugh Evans - 12 Aug, 2015)
1 None: Re: differences in TIDs produced by gras and mulassis   (Hugh Evans - 12 Aug, 2015)
... 1 Message(s)
2 None: Re: differences in TIDs produced by gras and mulassis   (Sveta - 18 Sep, 2015)
Question use integral particle flux as input  by Svetlana Shasharina <Svetlana Shasharina>,   28 Jan, 2014
Hi, I see that I need to set up my source using the shape of radiation, angular distribution and the energy spectrum. Quite often I see that only integral omnidirectional flux is available as data. How one woudl create a G4 particle source from this? I can assume isotropy.

Thanks, Sveta

Question Can we simulate the imaging process of Wolter-I type telescope?  by <zhaodh@bao.ac.cn>,   26 Aug, 2013
Hi everyone,

Can we simulate the imaging process of Wolter-I type X-ray ( a few keV) telescope and save the image using Geant4?

If we can, how to do it? Since there is no total internal reflection for X-ray in Geant4.

Looking forward to your reply. Thank you and best regards,

D.H Z

Question GPS interpolation  by Wolfgang Mai <Wolfgang Mai>,   25 Jun, 2013
Dear colleagues,

do you know how the general particle source is interpolating arbitrary spectra?

I was using GRAS 3.1 with an SPENVIS generated macro with an diff. electron spectrum (AE-8MAX, GEO), which looked like this:

/gps/particle e-
/gps/ene/type Arb
/gps/hist/type arb
/gps/hist/point 4.000000E-02 1.138400E+10
/gps/hist/point 1.000000E-01 5.653800E+09
/gps/hist/point 2.000000E-01 5.623300E+08
/gps/hist/point 3.000000E-01 1.717600E+08
/gps/hist/point 4.000000E-01 8.587400E+07
/gps/hist/point 5.000000E-01 5.386000E+07
/gps/hist/point 6.000000E-01 3.011800E+07
/gps/hist/point 7.000000E-01 2.467300E+07
/gps/hist/point 8.000000E-01 2.042700E+07
/gps/hist/point 1.000000E+00 1.383600E+07
/gps/hist/point 1.250000E+00 8.481900E+06
/gps/hist/point 1.500000E+00 4.907300E+06
/gps/hist/point 1.750000E+00 2.598900E+06
/gps/hist/point 2.000000E+00 1.649900E+06
/gps/hist/point 2.250000E+00 9.743300E+05
/gps/hist/point 2.500000E+00 4.355700E+05
/gps/hist/point 2.750000E+00 1.797200E+05
/gps/hist/point 3.000000E+00 8.610400E+04
/gps/hist/point 3.250000E+00 4.061000E+04
/gps/hist/point 3.500000E+00 2.145000E+04
/gps/hist/point 3.750000E+00 1.119400E+04
/gps/hist/point 4.000000E+00 4.665900E+03
/gps/hist/point 4.250000E+00 1.943400E+03
/gps/hist/point 4.500000E+00 8.365100E+02
/gps/hist/point 4.750000E+00 3.363000E+02
/gps/hist/point 5.000000E+00 5.920500E+01
/gps/hist/point 5.500000E+00 3.613200E+00
/gps/hist/point 6.000000E+00 0.000000E+00
#/gps/hist/point 6.500000E+00 0.000000E+00
#/gps/hist/point 7.000000E+00 0.000000E+00
/gps/hist/inter Lin
/gps/ang/type cos
/gps/ang/maxtheta 90.00 deg
/gps/pos/type Surface
/gps/shape Sphere
/gps/pos/centre 0.0000E+00 0.0000E+00 0.0000E+00 mm
/gps/pos/radius 5.100000E+00 mm

#
#Normalisation
#
/gras/analysis/setNormalisationType FLUENCE
/gras/analysis/setSourceFluence 2.0767E+17 cm-2

Now what I did was cutting the first value of the spectrum and changing the fluence value in order to get a better statistic, as these particles don't contribute remarkable to my dose measurement. The macro then looked like this:

...
/gps/hist/type arb
#/gps/hist/point 4.000000E-02 1.138400E+10
/gps/hist/point 1.000000E-01 5.653800E+09
...
/gras/analysis/setSourceFluence 4.6473E+16 cm-2

The strange thing is now, that the measured dose is only half the value of the first one, although these low energy electrons don't contribute.
After that I measured the fluence of these sources I saw, that there are much more electrons over 0.1 MeV in the first case than in the second. It seems, that gps distributes the particles in a different way as they really are distributed.

Does anyone know about this problem?

Thank you in advance
Wolfgang

Question Problem building GRAS code  by Leonardo Ghizoni <Leonardo Ghizoni>,   24 May, 2013
Hello,

I've downloaded the GRAS code to perform Geant4 space applications simulations. When building the code, I ran the command with cmake:

lghizoni@EQSAPhysics:~/Geant4/Gras-build$ cmake -DCMAKE_INSTALL_PREFIX=/home/lghizoni/Geant4/geant4.9.6-install/ -DWITH_ROOT=ON -DSTATIC_BUILD=ON -DGRAS_INSTALL_PREFIX=/home/lghizoni/Geant4/Gras-build/ /home/lghizoni/Geant4/Gras

That went OK, no errors! But then I ran 'make', and got the following error:

lghizoni@EQSAPhysics:~/Geant4/Gras-build$ make
Scanning dependencies of target gras_lib
[  0%] Building CXX object CMakeFiles/gras_lib.dir/source/analysis/src/GRASNIELAnalysisModuleMessenger.cc.o
[  1%] Building CXX object CMakeFiles/gras_lib.dir/source/analysis/src/GRASAnalysisManager.cc.o
[  1%] Building CXX object CMakeFiles/gras_lib.dir/source/analysis/src/GRASFluenceAnalysisModule.cc.o
[  1%] Building CXX object CMakeFiles/gras_lib.dir/source/analysis/src/GRASNIELCurveSetGaAsSummers.cc.o
[  2%] Building CXX object CMakeFiles/gras_lib.dir/source/analysis/src/GRASChargeCollectionAnalysisModule.cc.o
[  2%] Building CXX object CMakeFiles/gras_lib.dir/source/analysis/src/GRASEqDoseAnalysisModule.cc.o
[  3%] Building CXX object CMakeFiles/gras_lib.dir/source/analysis/src/GRASNIELCurveSet_JPL_NRL_NASA_InP.cc.o
[  3%] Building CXX object CMakeFiles/gras_lib.dir/source/analysis/src/GRASActivationAnalysisModule.cc.o
/home/lghizoni/Geant4/Gras/source/analysis/src/GRASActivationAnalysisModule.cc: In member function ‘virtual void GRASActivationAnalysisModule::EndOfRun(const GRASInfo*)’:
/home/lghizoni/Geant4/Gras/source/analysis/src/GRASActivationAnalysisModule.cc:127:73: error: cannot convert ‘std::map<CLHEP::Hep3Vector, CLHEP::Hep2Vector>*’ to ‘std::map<CLHEP::Hep3Vector, double>*’ in initialization
/home/lghizoni/Geant4/Gras/source/analysis/src/GRASActivationAnalysisModule.cc: In member function ‘virtual void GRASActivationAnalysisModule::PrintCSVFile(const GRASInfo*, SpenvisCSVCollection*)’:
/home/lghizoni/Geant4/Gras/source/analysis/src/GRASActivationAnalysisModule.cc:223:73: error: cannot convert ‘std::map<CLHEP::Hep3Vector, CLHEP::Hep2Vector>*’ to ‘std::map<CLHEP::Hep3Vector, double>*’ in initialization
In file included from /home/lghizoni/Geant4/geant4.9.6-install/include/Geant4/G4PhysicalConstants.hh:27:0,
                 from /home/lghizoni/Geant4/geant4.9.6-install/include/Geant4/globals.hh:70,
                 from /home/lghizoni/Geant4/Gras/source/analysis/include/GRASActivationAnalysisModule.hh:38,
                 from /home/lghizoni/Geant4/Gras/source/analysis/src/GRASActivationAnalysisModule.cc:35:
/home/lghizoni/Geant4/geant4.9.6-install/include/Geant4/CLHEP/Units/PhysicalConstants.h: At global scope:
/home/lghizoni/Geant4/geant4.9.6-install/include/Geant4/CLHEP/Units/PhysicalConstants.h:79:21: warning: ‘CLHEP::electron_charge’ defined but not used [-Wunused-variable]
In file included from /home/lghizoni/Geant4/geant4.9.6-install/include/Geant4/Randomize.hh:32:0,
                 from /home/lghizoni/Geant4/geant4.9.6-install/include/Geant4/G4SteppingManager.hh:62,
                 from /home/lghizoni/Geant4/geant4.9.6-install/include/Geant4/G4TrackingManager.hh:56,
                 from /home/lghizoni/Geant4/Gras/source/global/include/GRASInfo.hh:42,
                 from /home/lghizoni/Geant4/Gras/source/analysis/src/GRASActivationAnalysisModule.cc:39:
/home/lghizoni/Geant4/geant4.9.6-install/include/Geant4/CLHEP/Random/Randomize.h:59:12: warning: ‘CLHEP::HepRandomGenActive’ defined but not used [-Wunused-variable]
make[2]: *** [CMakeFiles/gras_lib.dir/source/analysis/src/GRASActivationAnalysisModule.cc.o] Error 1
make[1]: *** [CMakeFiles/gras_lib.dir/all] Error 2
make: *** [all] Error 2

Can someone help me?? I just don't know what to fix.

Thanks a lot

Leonardo.

1 None: Re: Problem building GRAS code   (Colin Paul Gloster - 24 May, 2013)
1 None: Re: Problem building GRAS code   (Leonardo Ghizoni - 24 May, 2013)
2 Idea: Re: Problem building GRAS code   (Vladimir Ivanchenko - 27 May, 2013)
Question Please, tell me any example   by Batom <Batom>,   03 May, 2013
Dear

Please, advise me any example for how to use the function ComputeSolid in the Parameterisation class?

Regards

Question Simulate Total Ionizing Dose and Single Event Effect  by Leonardo Ghizoni <Leonardo Ghizoni>,   09 Apr, 2013
Hello,

I'm new with Geant4, and I'm trying to build a simulation to calculate the total ionizing dose and some single event effects on components of a satellite.

I studied a lot the GPS Class, but I still don't know where to begin with...

I also already have my geometry complete and was able to build it with GPS...

Can someone help how to begin with and maybe have some plots too??

Thanks a lot.

Best regards.

1 None: RE: Simulate Total Ionizing Dose and Single Event Effect   (Susanna Guatelli - 09 Apr, 2013)
1 None: Re: RE: Simulate Total Ionizing Dose and Single Event Effect   (Leonardo Ghizoni - 10 Apr, 2013)
2 None: Re: Simulate Total Ionizing Dose and Single Event Effect   (Hugh Evans - 10 Apr, 2013)
1 None: Re: Simulate Total Ionizing Dose and Single Event Effect   (Leonardo Ghizoni - 10 Apr, 2013)
None Sector Shilding Analysis Tool installation  by T Ergin <T Ergin>,   10 Feb, 2013
Hi, 

I've recently downloaded the latest source code of SSAT from
     http://reat.space.qinetiq.com/ssat/
I have also installed Geant4 with the GDML option set to ON. 

But the ssat_env.csh script is not working, because the variable GDMLTOP is not set and it is not clear what directory it should point to. 

When I look into the what the script is doing with GDMLTOP is that is goes into the subdirectories of this package that GDMLTOP is supposed to point at. But it is not clear where exactly or what package exactly this variable must point at. In the directory structure of the installed Geant4, it doesn’t fit anywhere. Even the latest version of GDML from the Geant4 web site does not have such a subdirectory structure. 

Would you please let me know where exactly I should set this variable? 


 Thanks.

    T.
None Trouble with MULASSIS install  by Brian Evans <Brian Evans>,   18 Sep, 2012
Hi All,

I'll admit I'm a complete novice when it comes to building applications with Geant4, so please forgive me if I'm missing something obvious.

I trying to build MULASSIS 1.23, with Geant4 9.5 on 64-bit RedHat Enterprise Linux 6, but I keep running into errors. For starters, I've successfully built CLHEP 2.1.1.0, XercesC 3.1.1, Geant4 9.5, and GRAS 3.1, and installed each in their own directory in /opt/local.

I've unpacked the MULASSIS tar file in /tmp/swbuild, which is also where the source for each of the previous packages are located as well (also each in their own directory).

Now, I followed instructions in some other posts for setting up environment variables, and it seems the most important for MULASSIS is G4INSTALL, which I've set to /opt/local/geant4.

So, the errors: When I run "gmake" from the MULASSIS source directory, I get an error saying it can't find files in the "config" dir, which indeed does not exist in /opt/local/geant4. But it does exist in the geant4 source dir, so the error above seems to indicate it's looking for the geant4 source files, not the install directory, which seems at odds with the instructions I got to set G4INSTALL to point to the install dir.

Any idea what I'm doing wrong? I'd be glad to provide any further information if needed.

Thanks in advance for any assistance!

Question Simulating deep space radiation environment and normalization of the results  by Ali Ebrahimi <Ali Ebrahimi>,   08 Aug, 2012
Hello,

I would to simulate deep space radiation environment (GCR+SPE) using GPS for my application. I have got GCR proton flux from CREME and wrote the following macro for GPS and want to know if it is correct:

/gps/particle proton /gps/ene/type Arb /gps/hist/type arb /gps/hist/point 1.0023 0.0055595 /gps/hist/point 1.01623 0.023156 ... /gps/hist/point 100000 5.8474E-05

/gps/hist/inter Log

/gps/pos/type Surface /gps/pos/shape Sphere /gps/pos/centre 0. 0. 0. mm /gps/pos/radius 20. m

/gps/ang/type cos /gps/ang/mintheta 0. deg /gps/ang/maxtheta 90. deg

If it is correct, who can I normalize it? I have seen a couple of of other questions regarding normalizations and saw the files on http://geant4.in2p3.fr/2007/prog/prog.htm but it is still not clear to me in particular since I'm using arbitrary type . I would appreciate it if you could please elaborate on this or give a reference so I can read myself.

Also, from the previous topics in the space category and from GPS examples I learned that it is possible to have several particle source and simulate SPE at the same time. To simulate SPE, say the worst case, should I do it like GCR with arbitrary types? If yes is there a source for flux data?

Thank you very much in advance, Ali

None How to get Primary Particle Energy ?  by Mehmet <Mehmet>,   02 Jul, 2012
Dear all,

It is probably a very simple question, but I could not figure it out although I have been searching for several hours. Basically, I would like to get the primary particle's (say proton) energy, generated using GPS, at the time of injection and print it out with the information of the secondary (say muon). I would appreciate if anyone could tell me how to do or at least direct me to location with example.

Thank you, Mehmet

1 None: Re: How to get Primary Particle Energy ?   (John Allison - 02 Jul, 2012)
1 Feedback: Re: How to get Primary Particle Energy ?   (Mehmet - 03 Jul, 2012)
(_ Feedback: Re: How to get Primary Particle Energy ?   (Michael H. Kelsey - 03 Jul, 2012)
2 Feedback: Re: How to get Primary Particle Energy ?   (Neil Fazel - 02 Jul, 2012)
Question Simulating and calculating Geometric Factors in GEANT4  by Bruno Morgado <Bruno Morgado>,   26 Jun, 2012
Hello,
I'm trying to calculate the "real" Geometric Factor (GF) of an instrument in Geant, but I'm having some problems with the results.

My macro file is:
# Set GPS spherical source
/gps/particle e-			# particle type
/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/centre 0 0 0 m			# source centre
/gps/pos/radius 10 cm
/gps/ang/type cos
/gps/ene/mono 52. keV
/run/beamOn 200000000

And I finish calculating the GF with:
double geometryFactor(double numberDetectedParticles, double numberEmitedParticles, const double sourceRadius) {
    return ((8 * numberDetectedParticles * pow(TMath::Pi(), 2) * pow(sourceRadius, 2)) / numberEmitedParticles);
}

Still, like I said the results don't seem to add up. Can you direct me to some place where I can check if these are the correct calculations or if there is some mistake?
1 Idea: Re: Simulating and calculating Geometric Factors in GEANT4   (Laurent Desorgher - 02 Oct, 2012)
(_ None: Re: Simulating and calculating Geometric Factors in GEANT4   (Bruno Morgado - 08 Oct, 2012)
(_ None: Re: Simulating and calculating Geometric Factors in GEANT4   (Laurent Desorgher - 08 Oct, 2012)
Question Dose equivalent  by Geng <Geng>,   23 Jun, 2012
Dear experts:

I want to know something detail about the dose equivalent part.

I only have the method to calculate the absorb dose but not the dose equivalent.

Could you tell me the detail method in geant4 to realize the dose equivalent calculation with quality factor?

Or you may only send an example code of this part,I would study it by myself.

Thanks and best regards!

1 None: Re: Dose equivalent   (Youming Yang - 24 Jun, 2012)
(_ Question: Re: Dose equivalent   (Geng - 24 Jun, 2012)
Question How to get secondary particle (e.g. muon) flux?  by Halil <Halil>,   15 May, 2012
Hello Geant4 users;

I am working on the simulation of the cosmic muons using Geant4. I am trying to get flux at sea level using the primary particle flux which is given in unit (m2 sr s GeV)^-1 from balloon measurements. I shoot primary particles on the top of the atmosphere using GPS (general particle source) and record all the muons’ energy, angle etc. at sea level. The problem I am having is that I do not know how to express the muon's flux with the unit (cm2 sr s GeV)^-1.

Because of the possible scatterings and decays, muons spread a wide area at sea level with different angles. (For instance, when I shoot the protons vertically, muons reach the earth's surface with an angle in a range from 0 to 25 degree.) Therefore, I don’t know wheter I should count all the muons reaching the sea level regardless of their angle or I sould do correction by dividing the total number by the detector area etc. Namely, how can I convert the number of muon events I obtained into the muon flux in proper units for a comparison with the measurements?

Thanks in advance

Question Define "differential energy spectrum"  by Geng <Geng>,   11 Mar, 2012
Hi Sir/madam

   Thanks in advance!
   I am a fresh man in Geant4.I want to simulate a source with a "differential energy spectrum" , just like a GCR proton spectrum.
   Could you give me some advice to solute it?
   I need your suggestions urgently!
Thanks very much!
Geng

1 None: RE: Define "differential energy spectrum"   (Susanna Guatelli - 11 Mar, 2012)
2 None: Re: Define   (michel maire - 12 Mar, 2012)
Question 2011 Geant4 Space User Workshop  by wangym <wangym>,   17 Sep, 2011
hollo, I am a Geant4 user. I want to ask that Do 2011 Geant4 Space User Workshop call for paper ? I want to submit a paper, my teacher tell me if the paper has be received , I can attend the meeting.

1 Question: Re: 2011 Geant4 Space User Workshop   (wangym - 20 Sep, 2011)
Question Defining the integration time using a primary event number  by Manuel Castro Avila <Manuel Castro Avila>,   14 Oct, 2010
Hi everyone,

I like your help.


I have to simulate the behavior of an instrument during a exposition time to the radiation, in this case, to gamma photons. The instrument has a detector with size of 169 cm^{2} and 50 cm away from the detector, there is an aperture of 1332.25 cm^{2}.

So, I have the power law that represent the energy distribution of the photon spectrum that I go to use. The power law is the following:

5.16x10^{-2} E^{-1.81} ph cm^{-2} Sr^{-1} s^{-1} MeV^{-1}

This power law is valid between 0.04 MeV to 10 Mev.

I need to calculate the number of primary events that represent a exposition time, for instance, 4 hours. To calculate this number, I calculate the integral among the energy range and I get:

0.854 ph cm^{-2} Sr^{-1} s^{-1} MeV{-1}

-The integration time is: 4h=14400 s
-I want to integrate over 4*PI Sr (I considere a isotropic flux into my simulation using GPS from GEANT 4)
-The detector area is of 169 cm^{2}

Now, to get the total number photons, I multiply the number above (0.854 ph ...) by 14400*4*PI*169 and I get the total primary event number equals 26116617.68 photons.

My macro file is:

/gps/particle gamma
/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/ene/type Pow
/gps/ene/min 1 MeV
/gps/ene/max 10 MeV
/gps/ene/alpha -1.8

/run/beamOn 26116617

My question is: Is this procedure above correct?? Am I using the correct way to calculate the integration time??

By the other hand, the power law represent the photon background measured in the atmosphere to a atmospheric depth of 3.5 g cm^{-2}, is there a way to adapt this power law to the case where the atmospheric depth is of 2.8 g cm^{-2}???

Thanks,

Manuel

Question How to simulate GCR spectra   by aimsphere <aimsphere>,   15 Jun, 2010
Hello,

I have a question regarding the use of GPS for simulating the entire GCR. I wanted to get some suggestion regarding how to do this. I have made simulations for individual particle specie using mac files for each individually. Now since I would like to consider ions from Atomic no. 1 to 28 I was wondering if the only way to do is to make simulations for each ion using their respective individual mac files or is there some other way to do this?

Thank you.
1 Question: Re: How to simulate GCR spectra   (Geng - 24 Jun, 2012)
2 None: Re: How to simulate GCR spectra   (Hugh Evans - 25 Jun, 2012)
Question Normalization of an isotropic flux in space  by valentina <valentina>,   03 Mar, 2009
Dear all,

I am simulating the interaction of a space telescope with the cosmic ray environment and, after a long discussion with my colleagues, I am a bit confused on the normalization process that translates the number of counts in my detector in unity of counts/cm^2 sec.

So I really appreciate if somebody could help me to clarify this topic!

I show here what is my normalization way.

I have an isotropic proton flux F:

F = prot /cm^2 sec sr

I run the protons from a sphere of radius R within an half cone angle q to increase the statistics. My spacecraft is placed within the emission angle, that is the spacecraft is placed within the sphere of radius r, with

r = tan(q) x R

The GPS macro is:

/gps/particle proton

/gps/pos/type Surface

/gps/pos/shape Sphere

/gps/pos/radius 50.0 m ---> R

/gps/pos/centre 0.0 0.0 0.0 mm

/gps/ang/type iso

/gps/ang/mintheta 0 rad

/gps/ang/maxtheta 0.01 rad ---> q

while the spectral shape is made via /gps/hist .

At the end of the run, I detect a number of counts C on the focal plane. If N is the number of emitted protons, we call the fraction of detected/emitted particles Eff:

Eff = C/N

In order to trasform the number C to a rate P in counts/sec, I apply this formula:

P = Eff x F x (4 x PI) x PI x r^2 counts/sec

Then I divide the count rate P by the detecting area Adet, and the final flux B is given by:

B = P/Adet counts/cm^2 sec

The same normalization can be achieved using the simulation time:

D = F x A_sphere x omega prot/sec

where:

- A_sphere = 4 x PI x R^2 (the source sphere)

- omega = 2 x PI x (1 - cos(q)) (solide angle of the emission cone)

Then we calculate the simulated time:

Time = N/D (sec)

The flux B on the focal plane is given by:

B = C/(Time x Adet) counts/cm^2 sec

My normalization way is based on what shown by Sullivan (1971).

Now my question is: is this the right way?

Thank you all for the help!

Best regards, Valentina

1 Feedback: Re: Normalization of an isotropic flux in space   (Giovanni Santin - 03 Mar, 2009)
1 Question: Re: Normalization of an isotropic flux in space   (valentina - 03 Mar, 2009)
(_ Note: Re: Normalization of an isotropic flux in space   (Juan - 04 Mar, 2009)
(_ Ok: Re: Normalization of an isotropic flux in space   (valentina - 04 Mar, 2009)
(_ None: Re: Normalization of an isotropic flux in space   (Juan - 04 Mar, 2009)
2 None: Re: Normalization of an isotropic flux in space   (Paul Nicholas Colin Gloster - 03 Mar, 2009)
Question gps sphere source   by Daniel Braun <Daniel Braun>,   09 Oct, 2008
Hallo,

I need to use a sphere as a particle source with the gps class.

and I want the particles to shoot just in the inside of the sphere maximally tangential to the sphere. They also should be distributed even on the surface of the sphere.

What is the best way to realize this?

Thanks for any help.

Regards Daniel

1 Feedback: Re: gps sphere source   (Giovanni Santin - 10 Oct, 2008)
2 Feedback: Re: gps sphere source   (Giovanni Santin - 10 Oct, 2008)
1 None: Re: gps sphere source   (Daniel Braun - 10 Oct, 2008)
1 Feedback: Re: gps sphere source   (Giovanni Santin - 16 Oct, 2008)
2 None: Re: gps sphere source   (Daniel Braun - 30 Oct, 2008)
Question dose depth curve  by Daniel Braun <Daniel Braun>,   12 Sep, 2008
Hello everybody,

I tried to get a Dose depth curve in an Aluminium Sphere and compared it with the one I get from www.Spenvis.oma.be.

But for a thin shield I get not enough Energy deposited and for a thick shield I get to much Energy deposited.

Does anybody know what the Problem is

My Physiks are

// // // $Id: ExN03PhysicsList.cc,v 1.21 2007/07/02 13:22:08 vnivanch Exp $ // GEANT4 tag $Name: geant4-09-01-patch-02 $ // //

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

#include "ExN03PhysicsList.hh"

#include "G4ProcessManager.hh"
#include "G4ParticleTypes.hh"

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

ExN03PhysicsList::ExN03PhysicsList():  G4VUserPhysicsList()
{
  defaultCutValue = 1*mm;
  SetVerboseLevel(1);
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

ExN03PhysicsList::~ExN03PhysicsList() {}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void ExN03PhysicsList::ConstructParticle()
{
  // In this method, static member functions should be called
  // for all particles which you want to use.
  // This ensures that objects of these particle types will be
  // created in the program. 

  ConstructBosons();
  ConstructLeptons();
  ConstructMesons();
  ConstructBaryons();
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void ExN03PhysicsList::ConstructBosons()
{
  // pseudo-particles
  G4Geantino::GeantinoDefinition();
  G4ChargedGeantino::ChargedGeantinoDefinition();

  // gamma
  G4Gamma::GammaDefinition();

  // optical photon
  G4OpticalPhoton::OpticalPhotonDefinition();
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void ExN03PhysicsList::ConstructLeptons()
{
  // leptons
  G4Electron::ElectronDefinition();
  G4Positron::PositronDefinition();
  G4MuonPlus::MuonPlusDefinition();
  G4MuonMinus::MuonMinusDefinition();

  G4NeutrinoE::NeutrinoEDefinition();
  G4AntiNeutrinoE::AntiNeutrinoEDefinition();
  G4NeutrinoMu::NeutrinoMuDefinition();
  G4AntiNeutrinoMu::AntiNeutrinoMuDefinition();
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void ExN03PhysicsList::ConstructMesons()
{
 //  mesons
  G4PionPlus::PionPlusDefinition();
  G4PionMinus::PionMinusDefinition();
  G4PionZero::PionZeroDefinition();
  G4Eta::EtaDefinition();
  G4EtaPrime::EtaPrimeDefinition();
  G4KaonPlus::KaonPlusDefinition();
  G4KaonMinus::KaonMinusDefinition();
  G4KaonZero::KaonZeroDefinition();
  G4AntiKaonZero::AntiKaonZeroDefinition();
  G4KaonZeroLong::KaonZeroLongDefinition();
  G4KaonZeroShort::KaonZeroShortDefinition();
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void ExN03PhysicsList::ConstructBaryons()
{
  //  barions
  G4Proton::ProtonDefinition();
  G4AntiProton::AntiProtonDefinition();
  G4Neutron::NeutronDefinition();
  G4AntiNeutron::AntiNeutronDefinition();
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void ExN03PhysicsList::ConstructProcess()
{
  AddTransportation();
  ConstructEM();
  ConstructDecay();
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

#include "G4ComptonScattering.hh"
#include "G4GammaConversion.hh"
#include "G4PhotoElectricEffect.hh"

#include "G4MultipleScattering.hh"

#include "G4eIonisation.hh"
#include "G4eBremsstrahlung.hh"
#include "G4eplusAnnihilation.hh"

#include "G4MuIonisation.hh"
#include "G4MuBremsstrahlung.hh"
#include "G4MuPairProduction.hh"

#include "G4hIonisation.hh"

#include "G4LowEnergyIonisation.hh" 
#include "G4LowEnergyBremsstrahlung.hh" 

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void ExN03PhysicsList::ConstructEM()
{
  theParticleIterator->reset();
  while( (*theParticleIterator)() ){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    G4String particleName = particle->GetParticleName();

    if (particleName == "gamma") {
      // gamma         
      pmanager->AddDiscreteProcess(new G4PhotoElectricEffect);
      pmanager->AddDiscreteProcess(new G4ComptonScattering);
      pmanager->AddDiscreteProcess(new G4GammaConversion);
      //pmanager->AddProcess(new G4eBremsstrahlung,   -1, 3,3);

    } else if (particleName == "e-") {
      //electron
      pmanager->AddProcess(new G4MultipleScattering,-1, 1,1);
      pmanager->AddProcess(new G4eIonisation,       -1, 2,2);
      pmanager->AddProcess(new G4eBremsstrahlung,   -1, 3,3);      

    } else if (particleName == "e+") {
      //positron
      pmanager->AddProcess(new G4MultipleScattering,-1, 1,1);
      pmanager->AddProcess(new G4eIonisation,       -1, 2,2);
      pmanager->AddProcess(new G4eBremsstrahlung,   -1, 3,3);
      pmanager->AddProcess(new G4eplusAnnihilation,  0,-1,4);

    } else if( particleName == "mu+" || 
               particleName == "mu-"    ) {
      //muon  
      pmanager->AddProcess(new G4MultipleScattering,-1, 1,1);
      pmanager->AddProcess(new G4MuIonisation,      -1, 2,2);
      pmanager->AddProcess(new G4MuBremsstrahlung,  -1, 3,3);
      pmanager->AddProcess(new G4MuPairProduction,  -1, 4,4);       

    } else if ((!particle->IsShortLived()) &&
	       (particle->GetPDGCharge() != 0.0) && 
	       (particle->GetParticleName() != "chargedgeantino")) {
      //all others charged particles except geantino
      pmanager->AddProcess(new G4MultipleScattering,-1, 1,1);
      pmanager->AddProcess(new G4hIonisation,       -1, 2,2); 

    }
  }
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

#include "G4Decay.hh"

void ExN03PhysicsList::ConstructDecay()
{
  // Add Decay Process
  G4Decay* theDecayProcess = new G4Decay();
  theParticleIterator->reset();
  while( (*theParticleIterator)() ){
    G4ParticleDefinition* particle = theParticleIterator->value();
    G4ProcessManager* pmanager = particle->GetProcessManager();
    if (theDecayProcess->IsApplicable(*particle)) { 
      pmanager ->AddProcess(theDecayProcess);
      // set ordering for PostStepDoIt and AtRestDoIt
      pmanager ->SetProcessOrdering(theDecayProcess, idxPostStep);
      pmanager ->SetProcessOrdering(theDecayProcess, idxAtRest);
    }
  }
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void ExN03PhysicsList::SetCuts()
{
  if (verboseLevel >0){
    G4cout << "ExN03PhysicsList::SetCuts:";
    G4cout << "CutLength : " << G4BestUnit(defaultCutValue,"Length") << G4endl;
  }

  // set cut values for gamma at first and for e- second and next for e+,
  // because some processes for e+/e- need cut values for gamma
  //

G4double lowlimit=250*eV;
  G4ProductionCutsTable::GetProductionCutsTable()
    ->SetEnergyRange(lowlimit, 100.*GeV);

  // default cuts for world volume
  G4double cutValue = 1*mm;
  SetCutValue(cutValue,"gamma");
  SetCutValue(cutValue,"e-");
  SetCutValue(cutValue,"e+");

  if (verboseLevel>0) DumpCutValuesTable();
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

From my Detector Construction

I use aluminium as shield and Silicon as detector Material defined like this.

new G4Material("Aluminium", z=13., a=26.98*g/mole, density=2.70*g/cm3); new G4Material("Silicon" , z=14., a= 28.09*g/mole, density= 2.33*g/cm3);

I'm new with geant4 and don't have a clue where the mistake could be.

Thanks for any help

Kind regards

Daniel

1 None: Re: dose depth curve   (michel maire - 12 Sep, 2008)
(_ None: Re: dose depth curve   (Daniel Braun - 12 Sep, 2008)
(_ None: Re: dose depth curve   (Hugh Evans - 15 Sep, 2008)
(_ None: Re: dose depth curve   (Daniel Braun - 16 Sep, 2008)
(_ None: Re: dose depth curve   (Hugh Evans - 17 Sep, 2008)
(_ None: Re: dose depth curve   (Daniel Braun - 22 Sep, 2008)
None SPE Primary Particles   by Steven Avery <Steven Avery>,   21 Apr, 2008

I'm using the radioprotection example but I would like to only use solar particle events (SPE) as my primary particle. The example only has ascii files for galactic cosmic rays (GCR) data. How do I set the example to only generate SPE events?

1 Feedback: Re: SPE Primary Particles   (Giovanni Santin - 21 Apr, 2008)
Question Grazing angle proton scattering  by valentina <valentina>,   13 Feb, 2008

Hi, I have a question about the simulation of grazing angle proton scattering.

In the article:

"Low Angle Scattering of Protons on the XMM-Newton Optics and Effects on the On-Board CCD detectors"

Nartallo et al. 2001

Nartallo uses the multiple scattering to simulate the interaction of protons with a mirror reflecting surface layer (in this case XMM-Newton optics).

As a conseguence of:

"Grazing Angle Proton Scattering: effects on Chandra and XMM-Newton X-ray telescopes"

Dichter and Woolf, 2003

, where the authors show that the multiple scattering under predicts the expected proton flux at the focal plane,

Nartallo in the article:

"Update on the use of Geant4 for the simulation of low energy protons scattering off X-ray mirrors at grazing incidence angles"

implements a new physics process, G4FirsovScattering, in the G4 Monte Carlo code, to treat the scattering of low energy protons at grazing incidence angles.

Now the question is: someone knows where is G4FirsovScattering? I can't find it in G4 library!

Thanks!

Valentina

None Generating a Pulse height spectra through GEANT  Keywords: Getting a Pulse height spectra from GEANT
by rajani <rajani>,   02 Feb, 2007

Hi All!!

In a normal space experiment, the particles ( suppose photons here) interact with the detector, giving rise to the absorption or transmission of the particle. In GEANT, we can trace the physics processes it undergoes, the energy it carries after that & the step at which it does that. In a Lab- experiment, what we get out is a Pulse Height Spectra that tells the variation of no of counts with the voltage, My concern is, Can GEANT in any way tell me about the Pulse height spectra? If yes, Can you tell me how & the reference? Thanking You

Rajani

1 None: Re: Generating a Pulse height spectra through GEANT   (Maria Grazia Pia - 02 Feb, 2007)
Question Sunyaev-Zeldowich effect and GEANT4  Keywords: Cosmology, Astroparticle Physics
by Forough <Forough>,   22 Nov, 2006
   Hi,
   We are going to use GEANT4 in scattering of electrons
   from ionized atoms, so-called Sunyaev-Zeldowich effect
   in cosmology and astroparticle physics.
   Up to the best of our knowledge, this topic has not been
   worked until now. We are pleased to receive your comments
   about the application of GEANT4 in Sunyaev-Zeldowoch effect.
   Could you please help us? 
   Thanks if you reply.
   Sincerely,
   Forough Nasseri
Question hemisphere generation using General Particle Source Module  Keywords: hemisphere in gps
by Manju Sudhakar <Manju Sudhakar>,   28 Aug, 2006

Hi,

I would like to know how to generate a hemisphere using General Particle Source Module. It is easy to do so in the G4UserPrimaryGeneratorAction.cc, but I want to use gps because my energy is in the form of a user defined histogram. I checked the queries posted in this forum regarding generation of a hemisphere using gps, and I found out that "biasing" should be used and the extended example "exgps" has a macro "test32.g4mac, which gives the methodology to do so. But, I donot know how the biasing works and the gps manual doesnot give any details. Could someone please help me out in this.. It is really important, as I want to model the space-craft background due to photons generated on the inner surface of a hemisphere that is placed behind the space-craft, ie, the photons are incident from the hemisphere inner surface and impinging on the back of the detector.

Please could someone tell me how to go about this ??

Thanks, Manju Sudhakar.

None Physics Process to include to take care of x-ray polarization  by Manju Sudhakar <Manju Sudhakar>,   30 May, 2006

Hi,

I am working on a first step x-ray polarimeter design.I want to measure the polarization of low energy x-rays after scattering from a target. So, inorder to do that, I need to include, in the physics list a process that takes care of polarization. There is a physics process "G4LowEnergyPolarizedComptonScattering", but it is not used in any of the examples ..

And another thing is that in the http://geant4.esa.int/events/g4suw03/ site, there are a number of presentations that mention polarized photoelectric, Rayleigh scattering, pair production etc. And in the photoelectric presentation , there is a page where they show the actual track of the photon, and it is mentioned "G4LowEnergyPolarizedPhotoelectric".

I am using GEANT version 4.7.0.p01.
In the directory /geant4/geant4.7.0.p01/source/processes/electromagnetic/lowenergy/include, 
the only low energy polarization physics process is G4LowEnergyPolarizedComptonScattering
is there. Also, there is no PolarizationManager. In the paper by Ralph Dollan, Karim Laihem 
and Andreas Schalicke "Monte Carlo based studies of a polarized positron 
source for International Linear Collider (ILC)" (arXiv:physics/0512192 v1), they talk about 
polarization framework in GEANT4 .. They talk about "proposed implementation of polarized processes" 
 ... so are these yet to be implemented ?? In the latest version of GEANT, 
are these implemented ??

So, in the version of GEANT that I am using, 4.7.0.p01, the only physics process that is available is "G4LowEnergyPolarizedComptonScattering".

Could you please give me more information about the proposed implementation of polarization ? Is it already implemented ?

Thanks in advance, Manju Sudhakar.

1 None: Re: Physics Process to include to take care of x-ray polarization   (Maria Grazia Pia - 02 Jun, 2006)
None Geant4 applications at IEEE NSS 2006  by Maria Grazia Pia <Maria Grazia Pia>,   29 May, 2006

Geant4 users in space applications may be interested to submit an abstract on their ongoing projects to the IEEE Nuclear Science Symposium (http://www.nss-mic.org/2006). IEEE NSS is a major conference in experimental technologies, including a Software & Computing session. This year it will take place in San Diego, the week before the SPENVIS-Geant4 Space User Workshop also hosted in California in 2006.

The official deadline for abstract submission is over, but direct submission to the conveners is still accepted.

If you are interested, please send your abstract to the Software and Computing convener (Maria.Grazia.Pia@cern.ch, with cc to nss2006@ornl.gov).

Looking forward to many exciting Geant4 results,

Maria Grazia Pia

NSS Software & Computing co-Convener

None regarding polarization in gps  by Manju Sudhakar <Manju Sudhakar>,   25 May, 2006

Hi,

This is the first time I will be simulating polarization and I'll be doing it using gps.

The command /gps/polarizarion Px Py Pz where Px, Py and Pz are the doubles that sets the polarization of the source .. Are these the Stokes parameters ?? The parameters for linear polarization, circular polarization and elliptical polarization ..So if my source photons are linearly polarized only, then the values for Px, Py and Pz would be 1,0,0 .. Is my understanding correct or is this something else ...??? These are not vectors ..but doubles .. Could you help me out ??

Regards and thanks, Manju Sudhakar

 

Question How to define the flux of particles?  Keywords: flux
by Lee Xuesheng <chaosl@pku.edu.cn>,   06 Mar, 2006

When using the geant4 in application for Space Radiation,how can I define the flux of particles(e- or gamma)? thanks

1 None: Re: How to define the flux of particles?   (Hugh Evans - 07 Mar, 2006)
1 None: Re: How to define the flux of particles?   (Andres - 25 Mar, 2006)
3 None: Re: How to define the flux of particles?   (Andres - 25 Mar, 2006)
1 None: Re: How to define the flux of particles?   (Hugh Evans - 27 Mar, 2006)
...
2 None: Re: How to define the flux of particles?   (Lee Xuesheng - 05 Oct, 2006)
...
None Untitled  by Andres <Andres.Russu@uv.es>,   06 Feb, 2006

I have two questions about spacecraft geant simulation.

1.- The number of particles. I am using gps and sphere source with spenvis histogram but I am not sure how many particles per second I have to simulate.

2.- Long term radiation. I would like to simulate long term secondary particles coming from the desactivation of the ISS materials. Do you have any idea about how to do this?

Thanks

1 None: Re: Untitled   (Hugh Evans - 06 Feb, 2006)
(_ None: Re: Untitled   (Andres.Russu@uv.es - 22 Feb, 2006)
None Spacecraft activation from bkg environment  by <Andres.Russu@uv.es>,   31 Jan, 2006

Dear friends,

I am simulating a gamma instrument for the ISS. I would like to simulate how the spacecraft is getting radioactive in a long term simulation. I mean I want to take into account the particles coming form the decays of the materials around the detector plane.

Thanks in advance, Andres

Question doubt concerning the radioactive decay  Keywords: energy deposit by radioactive decay as a function of time
by Manju Sudhakar <Manju Sudhakar>,   02 Jan, 2006
Hi, 

   I would like to simulate the spallation reactions that occur CsI detectors 
due to cosmic ray protons.First I want to do an experiment on the energy deposited 
in CsI due to irradiation by a beam of monoenergetic protons, as a function of time.
I want to get the energy deposit spectrum in CsI after irradiation, as a function of time.

    I have looked at the examples/extended/radioactivity/exrdm, where the concept of 
theStep->GetPreStepPoint()->GetGlobalTime() is used.

Please, could someone help me in understanding this concept and how it is used in 
radioactive decay ????
How can I use this to get the energy deposit as a function of time since irradiation ???
Could someone please answer my query ...

With thanks, 
Manju Sudhakar.
Question doubts concerning grdm  Keywords: doubts about grdm commands
by Manju Sudhakar <Manju Sudhakar>,   29 Dec, 2005
Hi, 
  I would like to simulate the decay of induced radio-isotopes when protons 
are incident on a space-borne gamma-ray detector, cosmic-ray protons.I want 
to model the detector background.To simulate radioactive decay, we need to 
use grdm.There are some aspects about it I don't understand.

1) the command /grdm/sourceTimeProfile <filename>
    
    as given in the manual, it is for defining the "source time profile"
I donot understand ... what does this mean ?? Why is it used ?? Is it 
mandatory ?? In the extended example "radioactivity/exrdm" in geant 4.7.0, 
in one of the macros "proton-b.mac" /grdm/sourceTimeProfile beam.data is a 
command.When I opened beam.data, the file looked like this 
0  1e6
60  0

Donot understand head or tail of this ... please help me get clarified on 
this point.


2) the command /grdm/decayBiasProfile <filename>

       What about this command ?In the extended example "radioactivity/exrdm" 
in geant 4.7.0, in the macro "proton-b.mac", there is a command 
/grdm/decayBiasProfile measures.data.And, when I opened it
1000     1
1100     0
3000     1
3100     0
1e5      1
1.1e5    0

    No clue to what this file means.

Could someone please explain all this to me clearly and in detail ?

With lots and lots of advance thanks 
Manju Sudhakar.
None Re: Doubt regarding General Particle Source Module  by Hugh Evans <Hugh Evans>,   30 Aug, 2005
Hello Manju,

The energy units in the GPS histogram that you're building are in MeV. The
absolute values of the differential "flux" are not important as they are all
normalised to a probability distribution [0..1], therefore the relative
difference between bins in the histogram are of the greatest importance, not
the absolute value.

While cosmic ray fluxes aren't yet available, the trapped proton, trapped
electron and solar proton flux/fluence spectra from SPENVIS
(http://www.spenvis.oma.be/spenvis/) can be converted to GPS macros via its
Mulassis application.  By playing with this, you can easily determine how the
environmental spectra are converted to the GPS macro commands.

Alternatively, I can also send you an Excel spreadsheet within which the
CREME-86 spectra are implemented and the GPS macros are formed.

Regards,
Hugh

-----
H. Evans
ESA/ESTEC/TEC-EES
Postbus 299
2200 AG Noordwijk
The Netherlands

Phone:+31 (0)71 565-5109
Fax:     +31 (0)71 565-4999



|---------+--------------------------->
|         |           Manju Sudhakar  |
|         |           <manjus@isac.ern|
|         |           et.in>          |
|         |                           |
|         |           30/08/2005 09:25|
|         |                           |
|---------+--------------------------->
  >--------------------------------------------------------------------------------------------------------------------------|
  |                                                                                                                          |
  |        To:      space_app-g4hn@slac.stanford.edu                                                                         |
  |        cc:                                                                                                               |
  |        Subject: Doubt regarding General Particle Source Module                                                           |
  >--------------------------------------------------------------------------------------------------------------------------|




*** Discussion title: Space Applications
Email replies to space_app-g4hn@slac.stanford.edu must include:
  In-Reply-To: <"/space_app/3"@geant4-hn.slac.stanford.edu>
  Subject: ...change this to be about your reply.

Dear GEANT users ,

I have a doubt regarding one of the commands in General Particle Source
Module.

I have constructed a detector system mounted on a space platform, and I
wish to irradiate it with cosmic-ray protons, at 1AU.I have the spectrum
of the protons with me and I have to enter the energy and flux as a
user-defined arbitrary point-wise function of the form, /gps/ene/type
Arb /gps/hist/type arb /gps/hist/point ...(energy) ... (flux) .... ....
.... .... /gps/hist/inter Lin /run/beamOn 100000

        This is fine, but my doubts are as follows ,
1) What is the default value of the energy points ? Is it MeV ?
2) What is the default value of the differential flux to be entered as in the

   second column ?

Could you please answer my doubts. Regards, Manju Sudhakar.





1 None: Re: Doubt regarding General Particle Source Module   (Andres - 28 Sep, 2005)
Question Doubt regarding General Particle Source Module  Keywords: units of differential spectrum in gps.
by Manju Sudhakar <manjus@isac.ernet.in>,   30 Aug, 2005

Dear GEANT users ,

I have a doubt regarding one of the commands in General Particle Source Module.

I have constructed a detector system mounted on a space platform, and I wish to irradiate it with cosmic-ray protons, at 1AU.I have the spectrum of the protons with me and I have to enter the energy and flux as a user-defined arbitrary point-wise function of the form, /gps/ene/type Arb /gps/hist/type arb /gps/hist/point ...(energy) ... (flux) .... .... .... .... /gps/hist/inter Lin /run/beamOn 100000

        This is fine, but my doubts are as follows ,
1) What is the default value of the energy points ? Is it MeV ?
2) What is the default value of the differential flux to be entered as in the 
   second column ? 

Could you please answer my doubts. Regards, Manju Sudhakar.

None Geant4 Space Users Workshop  by Makoto Asai <Makoto Asai>,   02 Jun, 2005
Dear colleague,

On 3-7 October 2005, a combined workshop for the SPENVIS and 
GEANT4 Space applications communities will be held in the 
historic university town of Leuven, Belgium. The workshop venue 
is the Faculty Club of the Catholic University, situated in the 
beautifully restored Grand Beguinage, a UNESCO World Heritage 
Site.

The SPENVIS workshop will take place on 3-4 Oct, the Geant4 
workshop will follow on 6-7 Oct, with a joint session in the 
morning of 5 Oct.

The SPENVIS workshop is the second in a series that started at 
ESTEC in 2002. The emphasis during the SPENVIS workshop will 
be on user experience and requirements for future developments. 
We aim for an informal workshop with a limited number of 
presentations (and a poster session) which can serve as starting 
points for round table discussions addressing specific
developments and applications. The programme is deliberately 
kept open to allow users to actively participate in the 
organisation. If you have inputs or suggestions for the workshop 
programme, please send a note to D.Heynderickx@oma.be .

Several GEANT4 based radiation effects packages are now available 
or under development, which are or will be implemented in the 
SPENVIS user interface. In view of these developments, we opted 
for a combined workshop where users and developers of both 
communities can meet and interact.

The GEANT4 Space Users' workshop follows the events organised 
at ESTEC in 2003 and Vanderbilt University in 2004. The focus 
will be on new results for component, sensor, and shielding 
analyses, as well as on developments for Exploration and Science 
missions. Users are invited to highlight new requirements based 
on their applications.

Registration is now open at:
 http://www.spenvis.oma.be/spenvis/workshop/workshop.html 
Workshop participants are invited to submit abstracts before 
1st September, via e-mail attachment in pdf format to 
D.Heynderickx@oma.be . Please indicate for which part of the 
workshop you are submitting, and your preference for oral or 
poster presentation (subject to the final decision of the 
scientific committee). Abstracts and a more detailed programme
will be put on-line when available.

More information on SPENVIS and GEANT4 space applications can 
be found at
 http://www.spenvis.oma.be/ 
and
 http://geant4.esa.int/ 
respectively.


 Add Message Add Message
to: "Space Applications"

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 ]