[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

RE: Help in MCNP5



Yes, put an "e" card with your "f" card when you perform your tally. The

"e" card (Tally Energy card) can score your results in energy bins. I'm

looking at page 3-72, but that's of the MCNP4B manual. 



Mike



Michael G. Stabin, PhD, CHP

Assistant Professor of Radiology and Radiological Sciences 

Department of Radiology and Radiological Sciences 

Vanderbilt University 

1161 21st Avenue South

Nashville, TN 37232-2675 

Phone (615) 343-0068

Fax   (615) 322-3764

Pager (615) 835-5153

e-mail     michael.g.stabin@vanderbilt.edu 

internet   www.doseinfo-radar.com



 



-----Original Message-----

From: owner-radsafe@list.vanderbilt.edu

[mailto:owner-radsafe@list.vanderbilt.edu] On Behalf Of Pramoth CM

Sent: Tuesday, November 23, 2004 11:07 AM

To: radsafe@list.vanderbilt.edu

Subject: Help in MCNP5





Hi Group,

	I am a student working with MCNP5 and using a very small code to

obtain a 

tally of gamma ray spectrum of a particular gamma ray energy for a given



number of events or live time. I want to modify this code so that the

code 

would generate spectrums at increments at10 kev interval for 0 - 1000

kev. The only option of doing this which I can think of is running the

original 

code over and over at 10 kev increments - which would be time consuming.

Could anybody suggest the required change in the code, so that I could

run 

the code only once and then obtain a single OUTPUT file with a 100

columns 

( each column representing the spectrum at 10 kev intervals) . I am very



new to MCNP5 and read the required parts of the MCNP5 manual and could

not 

figure out how to do this. I will be really grateful if anybody could

help. 

I have provided the original code below.





c ************************** Cell definition cards

************************** c RightCircularCylinder Detector Assembly c

Sensor Cell Definition Cards c

     1   117  -3.67 (-1 2 -3 ) imp:p 1 $ NaI Sensor

c

c  End Cap

c

    20   213  -2.70   (84 -85 -82)                   imp:p 1 $ Al EndCap

Window

    21   444  -8.00   (82 -83 -85 81 ):(-80 81 -82 ) imp:p 1 $ SS EndCap

c

c  Vacuum

c

    32     0 (-84 1 -82)                            imp:p 1 $ Vacuum

c

c Everything else...

c

   988 220 -0.0012 (-900 -901 902 )(-81: 83: 85) &

                    &

                   imp:p 1 $ Detector environment

   999  0  (900: 901: -902)   imp:p 0 $ Outside world



c ************************* Surface definition cards

************************* c c Sensor Surface Definition Cards c

     1    pz      0.00 $ Top Face of sensor

     2    pz     -5.08 $ Bottom Face of sensor

     3    cz      2.54 $ Outer Radius of sensor

c

c End Cap and End Cap Window

c

    80    pz    -5.08 $ Inside Lower bound of detector end cap

    81    pz    -5.3975 $ Outside Lower bound of detector end cap

    82    cz     2.54 $ Inside radius of end cap

    83    cz     2.8575 $ Outside radius of end cap

    84    pz     0.5 $ Inside detector Face

    85    pz     0.55 $ Outside detector face

c

c Bounding surfaces

   900    cz  3.8575 $ Bounding cylinder

   901    pz  12 $ Bounding top Surface

   902    pz  -6.3975 $ Bounding bottom Surface



c ******************************** Data cards

******************************** c Mode P c c Material Cards c

m220     6000.01p    -0.000124 $ C   Air, rho = -0.001205 g/cc

          7000.01p    -0.755267 $ N

          8000.01p    -0.231781 $ O

         18000.01p    -0.012827 $ Ar

c

m117    11000.01p   1    $ NaI,  rho = -3.67 g/cc

         53000.01p   1

c

m444    6000.01p   -0.0024    $ Stainless Steel,  rho = -8 g/cc

         14000.01p   -0.064

         15000.01p   -0.0028

         16000.01p   -0.0016

         22000.01p   -0.032

         24000.01p   -1.44

         25000.01p   -0.16

         26000.01p   -5.5

         27000.01p   -0.002

         28000.01p   -0.8

c

m213    13000.01p   1    $ Aluminum,  rho = -2.7 g/cc

c

c Tally defination card(s)

c

f8:p  (1) $ Pulse Height tally on cell 1

c  Gaussian Energy Broadening Term

ft8 geb 0  0.05086  0.30486  $ for NAI

c

c Source definition cards

c

sdef  pos= 0  0  10  erg=d3

c

c Remember! The calculated efficiency is per photon emitted, so if you

have 

more

c than one photon energy, and they are all equally probable, multiply by

N c

  #   si3    sp3

      L       D

    1.0000   1.565E+15

c

c  Energy Range and bin width [in MeV!]

c

e0 0 1e-5 8e-3 512i 4.112 $ 8 KeV bins up to 4096 KeV

c

c Number of particles to follow

c

nps 67900000

c

c .OR. Minutes to run

c

c ctme   1

c

c Print & Dump Cycle card to create MCTAL output file

c

prdmp 2j 1



Thanks,

Pramoth





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

You are currently subscribed to the Radsafe mailing list. To

unsubscribe, send an e-mail to Majordomo@list.vanderbilt.edu  Put the

text "unsubscribe radsafe" (no quote marks) in the body of the e-mail,

with no subject line. You can view the Radsafe archives at

http://www.vanderbilt.edu/radsafe/



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

You are currently subscribed to the Radsafe mailing list. To

unsubscribe, send an e-mail to Majordomo@list.vanderbilt.edu  Put the

text "unsubscribe radsafe" (no quote marks) in the body of the e-mail,

with no subject line. You can view the Radsafe archives at

http://www.vanderbilt.edu/radsafe/