Message: Re: Other questions through TestEm7 Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Re: Other questions through TestEm7 

Keywords: TestEm7 absorbed dose
Forum: Medical Applications
Re: Sad the calculation of absorbed dose through TestEm7 (Shahrokh)
Re: None Re: the calculation of absorbed dose through TestEm7 (michel maire)
Re: Question Re: Other questions through TestEm7 (Shahrokh)
Re: None Re: Other questions through TestEm7 (michel maire)
Date: 29 Aug, 2004
From: Shahrokh <Shahrokh>

1) You amazed me with your answer since you could run dose program
in both conditions and got the same behavior. 
Perhaps it needs to tell you exactly whatever I do. 
Condition 1: I define n_particle=1 at PrimaryGeneratorAction.cc 
and UI command  /run/beamOn 100000 at macrofile and 
then compile and execute with the command:
#./dose rundose.mac. I take the following output:
Cumulated Doses: 	Edep      	Edep/Ebeam 	Dose  	
tally 0: 	 	18.0511 keV 	1.46757e-05 % 	2.8921e-12Gy 
tally 1: 	 	0 eV  	        0 % 	        0 Gy  	
tally 2: 	 	548.862 eV  	4.46229e-07 % 	8.79374e-14 Gy 
tally 3: 	 	0 eV  	        0 % 	 	0 Gy 	
tally 4: 	 	10.5051 keV 	8.54077e-06 % 	1.68311e-12 Gy 
tally 5: 	 	0 eV  	        0 % 		0 Gy  	
tally 6: 	 	0 eV  	        0 % 	 	0 Gy  	
tally 7: 	 	0 eV  	        0 % 	 	0 Gy
tally 8: 	 	307.583 keV 	0.000250068 % 	4.92803e-11 Gy 
tally 9: 	 	17.2944 keV 	1.40605e-05 % 	2.77087e-12Gy
tally 10: 	 	0 eV  	        0 %             0 Gy 	
tally 11: 	 	0 eV  	        0 %             0 Gy
tally 12: 	 	0 eV  	        0 %           	0 Gy  	
tally 13: 	 	40.3388 keV 	3.27958e-05 % 	6.463e-12 Gy
tally 14: 	 	0 eV  	        0 %           	0 Gy 	
tally 15: 	 	0 eV  	        0 %           	0 Gy 
tally 16: 	 	0 eV  	        0 %           	0 Gy  	
tally 17: 	 	0 eV     	0 %            	0 Gy  
tally 18: 	 	0 eV  	        0 %           	0 Gy 	
tally 19: 	 	0 eV     	0 %           	0 Gy  	

Condition 2: I define again n_particle=100000 and /run/beamOn 1 
and then compile and execute with previous command. 
I take the following output: 
Cumulated Doses: 	Edep    Edep/Ebeam 	Dose  	
tally 0: 	 	0 eV  	0 % 	 	0 Gy 	 
tally 1: 	 	0 eV  	0 % 	 	0 Gy 	 
tally 2: 	 	0 eV  	0 % 	 	0 Gy 	 
tally 3: 	 	0 eV  	0 % 	 	0 Gy 	 
tally 4: 	 	0 eV  	0 % 	 	0 Gy 	 
tally 5: 	 	0 eV  	0 % 	 	0 Gy 	 
tally 6: 	 	0 eV  	0 % 	 	0 Gy 	 
tally 7: 	 	0 eV  	0 % 	 	0 Gy 	 
tally 8: 	 	0 eV  	0 %             0 Gy 
tally 9: 	 	0 eV  	0 %             0 Gy 
tally 10: 	 	0 eV  	0 %             0 Gy 
tally 11: 	 	0 eV  	0 %           	0 Gy 
tally 12: 	 	0 eV  	0 %            	0 Gy
tally 13: 	 	0 eV  	0 %            	0 Gy 
tally 14: 	 	0 eV  	0 %            	0 Gy 
tally 15: 	 	0 eV  	0 %            	0 Gy 
tally 16: 	 	0 eV  	0 % 	  	0 Gy 
tally 17: 	 	0 eV  	0 % 	 	0 Gy 
tally 18: 	 	0 eV  	0 % 	 	0 Gy 
tally 19: 	 	0 eV  	0 % 	 	0 Gy 	 	
The tallies in the same positions are filled with the values that
differ with each other!!!
If the number of events processed is less than 10000, the energy deposited
in all tallies will be equal to zero. 
In other words, I think whatever the number of events processed increases;
the number of the tallies that their Edep are equal to zero decreases. 
Why is it?

Also, I didn't understand your mind from this sentence that 
“for the computation of the ratio, be care of EbeamCumul”.
Could you explain me a bit more?

2)I think that the energy deposited at a tally must decrease 
with increase in distance between particle generator and the tally. 
But in practice, there was a little variation between the simulation 
results and the theory. 
Is this affected by circular source instead of point source?
For example: n_particle=1,  /run/beamOn 10000000

tally number and position   Edep      	 Edep/Ebeam 	Dose  	
tally 0   (-10 0 0 cm) 	  7.66043 MeV    6.22799e-05 % 	1.22734e-09Gy 
tally 1    ( -9 0 0 cm)	  4.91947 MeV    3.99957e-05 % 	7.88186e-10Gy 
tally 2    ( -8 0 0 cm)   7.17745 MeV    5.83532e-05 % 	1.14995e-09Gy 
tally 3    ( -7 0 0 cm)   10.9176 MeV    8.87606e-05 % 	1.74919e-09Gy 
tally 4    ( -6 0 0 cm)   4.31018 MeV    3.50421e-05 % 	6.90567e-10Gy 
tally 5    ( -5 0 0 cm)   7.63259 MeV    6.20536e-05 % 	1.22288e-09Gy 
tally 6    ( -4 0 0 cm)   3.39132 MeV    2.75717e-05 % 	5.43349e-10Gy 
tally 7    ( -3 0 0 cm)   4.88245 MeV    3.96947e-05 % 	7.82254e-10Gy 
tally 8    ( -2 0 0 cm)   4.13689 MeV    3.36332e-05 % 	6.62803e-10Gy 
tally 9    ( -1 0 0 cm)   5.45894 MeV    4.43817e-05 % 	8.7462e-10Gy 
tally 10   ( 0 0 0 cm) 	  3.21238 MeV    2.6117e-05 % 	5.14681e-10Gy 
tally 11   ( 1 0 0 cm) 	  5.12507 MeV    4.16672e-05 % 	8.21127e-10Gy 
tally 12   ( 2 0 0 cm) 	  2.29093 MeV    1.86254e-05 % 	3.67047e-10Gy 
tally 13   ( 3 0 0 cm) 	  4.49537 MeV    3.65478e-05 % 	7.20239e-10Gy 
tally 14   ( 4 0 0 cm) 	  909.778 keV    7.39657e-06 % 	1.45763e-10Gy 
tally 15   ( 5 0 0 cm) 	  3.61464 MeV    2.93873e-05 % 	5.79129e-10Gy 
tally 16   ( 6 0 0 cm) 	  3.28858 MeV    2.67364e-05 % 	5.26889e-10Gy 
tally 17   ( 7 0 0 cm) 	  2.45891 MeV    1.99912e-05 % 	3.93962e-10Gy 
tally 18   ( 8 0 0 cm) 	  2.70589 MeV    2.19991e-05 % 	4.33531e-10Gy 
tally 19   ( 9 0 0 cm) 	  2.26599 MeV    1.84227e-05 % 	3.63051e-10Gy 
3) I expected to occur build up dose in within phantom and in the near
distant from phantom surface.
   

4) Can particle generator be defined in out of world? 
In definition of tally position in rundose.mac, 
if its position was in out of phantom, the executation of program 
would abort. 
This result wouldn’t occur if particle generator was in out of world, 
why is it?
5) How can several particle generators be defined that bombard phantom, simultaneously?

Thank you

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

 Add Message Add Message
to: "Re: Other questions through TestEm7"

 Subscribe Subscribe

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


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