Message: Re: ./rdecay02 proton-beam.mac output file Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: ./rdecay02 proton-beam.mac output file 

Forum: Documentation and Examples
Re: Question ./rdecay02 proton-beam.mac output file (Cristian Bungau)
Date: 19 Feb, 2016
From: Laurent Desorgher <Laurent Desorgher>

Dear Christian,

As you probably have noticed this example is to illustrate how the activation of a target from a high energy proton beam can be computed with Geant4.

For some obscure reason the proton nuclear physics was switched off (uncommented in the cc code) in the physics lits of the example. Leading to no activation at all, and therefore an empty activity results file.

If I understand well your message you have found this yourself and have activated the nuclear proton physics in the cc code? If you indeed do such and rerun the code you should get activation results like the one produced below:

Radioactivities in decay window no. 0 Z A E Activity (decays/window) Error (decays/window) 21 48 0 1.43956e-09 1.43956e-09 32 71 0 1.23432e-10 1.23432e-10 34 75 0 6.68183e-12 6.68183e-12 33 77 0 7.41712e-10 7.41712e-10 .........

Radioactivities in decay window no. 1 Z A E Activity (decays/window) Error (decays/window) 2 6 0 0 0 13 31 0 0 0 14 31 0 2.72891e-08 2.72891e-08 26 59 0 7.21078e-11 7.21078e-11 26 60 0 3.35337e-18 3.35337e-18 27 60 0 1.64548e-23 1.64548e-23 33 76 0 2.91363e-09 2.91363e-09 37 84 0 9.77523e-11 9.77523e-11 39 88 0 3.00937e-11 3.00937e-11

Radioactivities in decay window no. 2 Z A E Activity (decays/window) Error (decays/window) 2 6 0 0 0 21 50 0 1.15434e-35 1.15434e-35 .....

This file gives the number of decays per time windows for each radioactive products per time window . To get the activity you need to divide this number by the width of the time windows.

Best regards

Laurent

On Sat, 06 Feb 2016 18:54:18 GMT, Cristian Bungau wrote:

> Hello,
> 
> I am trying to run the rdecay02 example and I get an output file with no
> useful entries:
> 
> ./rdecay02 proton-beam.mac
> 
> produces the proton-beam.activity file which contains
> 
> Cristians-iMac:rdecay02-build cristian$ more proton-beam.activity
> 
> Radioactivities in decay window no. 0
> 
> Z A E Activity (decays/window) Error (decays/window)
> 
> Radioactivities in decay window no. 1
> 
> Z A E Activity (decays/window) Error (decays/window)
> 
> Radioactivities in decay window no. 2
> 
> Z A E Activity (decays/window) Error (decays/window)
> 
> Radioactivities in decay window no. 3
> 
> Z A E Activity (decays/window) Error (decays/window)
> 
> Yes, I have activated the hadronic processes in the PhysicsList and I do
> get out messages during the run:
> 
>  Mean energy deposit in target,   in time window = 18.939 MeV;  rms = 30.027 MeV
>  Mean energy deposit in detector, in time window = 4.7488 MeV;  rms = 24.593 MeV
> 
>  Process calls frequency in target :
>                 Decay=    233      RadioactiveDecay=   2011        Transportation=  24155
>               annihil=    545                 compt=   1458                  conv=    200
>                 eIoni=  13422                 hIoni=    813            hadElastic=   1084
>               ionIoni=   2660                   msc=   1053                muIoni=     11
>      neutronInelastic=    207                  phot=   4333          pi+Inelastic=      8
>          pi-Inelastic=      4       protonInelastic=    647 
> 
>  Process calls frequency in detector:
>                 Decay=    110      RadioactiveDecay=   3476        Transportation=  16642
>               annihil=    158                 compt=   4302                  conv=    153
>                 eIoni=   5639                 hIoni=    460            hadElastic=   2438
>               ionIoni=   3543                   msc=   1774                muIoni=     54
>              nCapture=     11      neutronInelastic=    953                  phot=   1038
>          pi+Inelastic=     22          pi-Inelastic=     14       protonInelastic=     38
> 
>  List of generated particles in target:
>           Ag105:       1  Emean =  67.869 MeV   ( 67.869 MeV --> 67.869 MeV)
>           Ag110:       1  Emean =  68.943 MeV   ( 68.943 MeV --> 68.943 MeV)
>           Ag111:       1  Emean =  37.268 MeV   ( 37.268 MeV --> 37.268 MeV)
>            As75:       1  Emean =   80.28 MeV   ( 80.28 MeV --> 80.28 MeV)
>            As78:       1  Emean =  84.192 MeV   ( 84.192 MeV --> 84.192 MeV)
>           Au182:       1  Emean =  1.0584 MeV   ( 1.0584 MeV --> 1.0584 MeV)
>           Au184:       1  Emean =  1.3786 MeV   ( 1.3786 MeV --> 1.3786 MeV)
>    Au185[8.900]:       2  Emean =  2.6959 MeV   ( 240.87 keV --> 5.1509 MeV)
>           Au186:       3  Emean =  3.3054 MeV   ( 1.745 MeV --> 5.3904 MeV)
>           Au187:       2  Emean =    2.43 MeV   ( 1.3988 MeV --> 3.4611 MeV)
>           Au188:       4  Emean =  5.1658 MeV   ( 2.381 MeV --> 7.8083 MeV)
>           Au189:       4  Emean =  2.7413 MeV   ( 1.1048 MeV --> 4.0091 MeV)
>    Au189[9.920]:       1  Emean =  351.02 keV   ( 351.02 keV --> 351.02 keV)
>           Au190:       7  Emean =  1.1151 MeV   ( 162.01 keV --> 3.417 MeV)
>           Au191:       3  Emean =  2.3402 MeV   ( 555.03 keV --> 4.9121 MeV)
>           Au192:       3  Emean =  1.5772 MeV   ( 549.85 keV --> 2.1798 MeV)
>           Au193:       2  Emean =  5.4543 MeV   ( 1.0481 MeV --> 9.8604 MeV)
>           Au194:       1  Emean =  7.5994 MeV   ( 7.5994 MeV --> 7.5994 MeV)
>           Au195:       1  Emean =  662.04 keV   ( 662.04 keV --> 662.04 keV)
>           Au201:       1  Emean =  836.96 keV   ( 836.96 keV --> 836.96 keV)
>           Au203:       1  Emean =  546.55 keV   ( 546.55 keV --> 546.55 keV)
>    Au204[0.161]:       1  Emean =  2.0761 MeV   ( 2.0761 MeV --> 2.0761 MeV)
>           Bi195:       2  Emean =  326.79 keV   ( 199.17 keV --> 454.4 keV)
>           Bi196:       1  Emean =  678.88 keV   ( 678.88 keV --> 678.88 keV)
>           Bi197:       5  Emean =  1.1179 MeV   ( 164.21 keV --> 3.4634 MeV)
>           Bi198:       1  Emean =  2.2176 MeV   ( 2.2176 MeV --> 2.2176 MeV)
>           Bi199:       3  Emean =  1.3622 MeV   ( 1.179 MeV --> 1.5187 MeV)
>           Bi200:       2  Emean =  1.4773 MeV   ( 1.4337 MeV --> 1.5208 MeV)
>           Bi201:       1  Emean =    81.4 keV   ( 81.4 keV --> 81.4 keV)
>           Bi202:       6  Emean =  302.61 keV   ( 136.3 keV --> 827.05 keV)
>           Bi203:       3  Emean =   401.1 keV   ( 27.951 keV --> 672.28 keV)
>           Bi204:       4  Emean =  304.73 keV   ( 102.35 keV --> 542.45 keV)
>           Bi205:       4  Emean =  202.85 keV   ( 77.354 keV --> 362.59 keV)
>           Bi206:       4  Emean =  235.44 keV   ( 30.038 keV --> 423.81 keV)
>           Bi207:       2  Emean =  129.55 keV   ( 24.068 keV --> 235.02 keV)
>           Bi208:       1  Emean =  35.615 keV   ( 35.615 keV --> 35.615 keV)
>           Bi209:       1  Emean =  22.786 keV   ( 22.786 keV --> 22.786 keV)
>            Br77:       2  Emean =  59.013 MeV   ( 54.322 MeV --> 63.703 MeV)
>            Br78:       1  Emean =  94.333 MeV   ( 94.333 MeV --> 94.333 MeV)
>            Br80:       1  Emean =  84.936 MeV   ( 84.936 MeV --> 84.936 MeV)
>            Br82:       1  Emean =  62.702 MeV   ( 62.702 MeV --> 62.702 MeV)
>            Br83:       1  Emean =  93.001 MeV   ( 93.001 MeV --> 93.001 MeV)
>             C11:       1  Emean =   110.5 MeV   ( 110.5 MeV --> 110.5 MeV)
>             C15:       1  Emean =  38.224 MeV   ( 38.224 MeV --> 38.224 MeV)
>            Ca43:       1  Emean =  142.93 MeV   ( 142.93 MeV --> 142.93 MeV)
>           Cd108:       1  Emean =  81.354 MeV   ( 81.354 MeV --> 81.354 MeV)
>           Cd111:       1  Emean =  59.986 MeV   ( 59.986 MeV --> 59.986 MeV)
>           Cd113:       2  Emean =  42.783 MeV   ( 24.531 MeV --> 61.035 MeV)
>           Cd114:       2  Emean =  42.121 MeV   ( 39.357 MeV --> 44.886 MeV)
>            Co62:       1  Emean =   85.48 MeV   ( 85.48 MeV --> 85.48 MeV)
>           Cs125:       1  Emean =  12.885 MeV   ( 12.885 MeV --> 12.885 MeV)
>           Cs126:       1  Emean =  9.0905 MeV   ( 9.0905 MeV --> 9.0905 MeV)
>            Cu63:       1  Emean =  101.07 MeV   ( 101.07 MeV --> 101.07 MeV)
>            Ga70:       2  Emean =  93.824 MeV   ( 86.966 MeV --> 100.68 MeV)
>            Ga75:       1  Emean =  96.096 MeV   ( 96.096 MeV --> 96.096 MeV)
>            Ge71:       1  Emean =  85.701 MeV   ( 85.701 MeV --> 85.701 MeV)
>            Ge72:       2  Emean =  71.497 MeV   ( 68.078 MeV --> 74.917 MeV)
>            Ge74:       2  Emean =  77.607 MeV   ( 65.839 MeV --> 89.374 MeV)
>            Ge76:       2  Emean =  86.874 MeV   ( 80.179 MeV --> 93.569 MeV)
>             He3:      21  Emean =  51.976 MeV   ( 15.06 MeV --> 112.52 MeV)
>           Hf170:       2  Emean =  2.9443 MeV   ( 2.8676 MeV --> 3.0211 MeV)
>           Hf174:       1  Emean =  1.0796 MeV   ( 1.0796 MeV --> 1.0796 MeV)
>           Hg184:       2  Emean =  4.1154 MeV   ( 3.5344 MeV --> 4.6963 MeV)
>           Hg185:       1  Emean =  1.5178 MeV   ( 1.5178 MeV --> 1.5178 MeV)
>           Hg186:       1  Emean =  1.1794 MeV   ( 1.1794 MeV --> 1.1794 MeV)
>           Hg187:       5  Emean =  2.1247 MeV   ( 574.51 keV --> 3.3553 MeV)
>           Hg188:       1  Emean =  2.8868 MeV   ( 2.8868 MeV --> 2.8868 MeV)
>           Hg189:       5  Emean =  4.0627 MeV   ( 2.7552 MeV --> 5.5772 MeV)
>           Hg190:      16  Emean =  2.7563 MeV   ( 833.07 keV --> 6.8775 MeV)
>           Hg191:       4  Emean =  2.3036 MeV   ( 503 keV --> 4.4017 MeV)
>           Hg192:      11  Emean =   2.045 MeV   ( 370.18 keV --> 5.1467 MeV)
>           Hg193:      12  Emean =   2.374 MeV   ( 288.7 keV --> 8.196 MeV)
>           Hg194:       7  Emean =  1.7537 MeV   ( 71.81 keV --> 4.4129 MeV)
>           Hg195:       6  Emean =   1.684 MeV   ( 327.67 keV --> 3.4164 MeV)
>           Hg196:       5  Emean =   1.551 MeV   ( 110.63 keV --> 4.1166 MeV)
>           Hg197:       7  Emean =  1.6293 MeV   ( 372.31 keV --> 3.0024 MeV)
>           Hg198:       4  Emean =  685.76 keV   ( 271.65 keV --> 1.4871 MeV)
>           Hg199:       4  Emean =  1.2919 MeV   ( 503.85 keV --> 2.2238 MeV)
>           Hg200:       1  Emean =  698.96 keV   ( 698.96 keV --> 698.96 keV)
>           Hg201:       1  Emean =  1.0988 MeV   ( 1.0988 MeV --> 1.0988 MeV)
>           Hg202:       1  Emean =  255.44 keV   ( 255.44 keV --> 255.44 keV)
>           Hg203:       1  Emean =  404.76 keV   ( 404.76 keV --> 404.76 keV)
>           Hg204:       1  Emean =  978.58 keV   ( 978.58 keV --> 978.58 keV)
>           Hg205:       1  Emean =  76.767 keV   ( 76.767 keV --> 76.767 keV)
>           Hg206:       1  Emean =  1.3796 MeV   ( 1.3796 MeV --> 1.3796 MeV)
>           In113:       1  Emean =  48.077 MeV   ( 48.077 MeV --> 48.077 MeV)
>           In118:       1  Emean =  48.659 MeV   ( 48.659 MeV --> 48.659 MeV)
>           Ir177:       1  Emean =  986.36 keV   ( 986.36 keV --> 986.36 keV)
>           Ir179:       2  Emean =  4.0903 MeV   ( 811.62 keV --> 7.369 MeV)
>           Ir180:       1  Emean =  2.7305 MeV   ( 2.7305 MeV --> 2.7305 MeV)
>           Ir181:       5  Emean =  4.7688 MeV   ( 1.7941 MeV --> 7.6972 MeV)
>           Ir182:       2  Emean =  4.9992 MeV   ( 4.3202 MeV --> 5.6782 MeV)
>           Ir183:       1  Emean =  2.4486 MeV   ( 2.4486 MeV --> 2.4486 MeV)
>           Ir184:       5  Emean =  1.2979 MeV   ( 355.41 keV --> 2.6711 MeV)
>           Ir185:       2  Emean =  3.5668 MeV   ( 195.76 keV --> 6.9378 MeV)
>    Ir185[5.800]:       1  Emean =  1.0252 MeV   ( 1.0252 MeV --> 1.0252 MeV)
>           Ir186:       1  Emean =  5.4325 MeV   ( 5.4325 MeV --> 5.4325 MeV)
>           Ir187:       1  Emean =  2.1736 MeV   ( 2.1736 MeV --> 2.1736 MeV)
>           Ir188:       1  Emean =  2.5296 MeV   ( 2.5296 MeV --> 2.5296 MeV)
>            Kr80:       2  Emean =  43.091 MeV   ( 36.17 MeV --> 50.011 MeV)
>            Kr82:       1  Emean =  75.971 MeV   ( 75.971 MeV --> 75.971 MeV)
>            Kr83:       1  Emean =  70.643 MeV   ( 70.643 MeV --> 70.643 MeV)
>            Kr85:       1  Emean =  48.998 MeV   ( 48.998 MeV --> 48.998 MeV)
>            Kr86:       1  Emean =  97.198 MeV   ( 97.198 MeV --> 97.198 MeV)
>           La130:       1  Emean =  30.026 MeV   ( 30.026 MeV --> 30.026 MeV)
>             Li6:       1  Emean =  26.759 MeV   ( 26.759 MeV --> 26.759 MeV)
>           Lu163:       1  Emean =  22.074 MeV   ( 22.074 MeV --> 22.074 MeV)
>           Lu167:       1  Emean =  2.6118 MeV   ( 2.6118 MeV --> 2.6118 MeV)
>            Mn54:       1  Emean =  52.155 MeV   ( 52.155 MeV --> 52.155 MeV)
>            Mn56:       1  Emean =  74.788 MeV   ( 74.788 MeV --> 74.788 MeV)
>            Mo95:       1  Emean =  71.223 MeV   ( 71.223 MeV --> 71.223 MeV)
>            Mo96:       1  Emean =  41.846 MeV   ( 41.846 MeV --> 41.846 MeV)
>             N14:       1  Emean =  67.112 MeV   ( 67.112 MeV --> 67.112 MeV)
>            Nb93:       3  Emean =  64.866 MeV   ( 30.511 MeV --> 99.553 MeV)
>            Nb94:       2  Emean =  63.055 MeV   ( 57.247 MeV --> 68.862 MeV)
>            Nb95:       1  Emean =  48.249 MeV   ( 48.249 MeV --> 48.249 MeV)
>            Nb96:       3  Emean =  79.922 MeV   ( 65.09 MeV --> 91.402 MeV)
>            Nb98:       1  Emean =  41.265 MeV   ( 41.265 MeV --> 41.265 MeV)
>            Nb99:       1  Emean =  51.021 MeV   ( 51.021 MeV --> 51.021 MeV)
>            Ni61:       1  Emean =  89.637 MeV   ( 89.637 MeV --> 89.637 MeV)
>           Os173:       1  Emean =  10.031 MeV   ( 10.031 MeV --> 10.031 MeV)
>           Os174:       1  Emean =  448.84 keV   ( 448.84 keV --> 448.84 keV)
>           Os175:       2  Emean =  1.3112 MeV   ( 783.64 keV --> 1.8388 MeV)
>           Os176:       3  Emean =  3.0598 MeV   ( 1.8345 MeV --> 4.0593 MeV)
>           Os178:       7  Emean =   5.546 MeV   ( 283.6 keV --> 14.288 MeV)
>           Os179:       1  Emean =  3.2551 MeV   ( 3.2551 MeV --> 3.2551 MeV)
>           Os180:       4  Emean =   2.816 MeV   ( 585.47 keV --> 6.2229 MeV)
>           Os181:       3  Emean =  5.3332 MeV   ( 2.2455 MeV --> 8.8751 MeV)
>           Os182:       2  Emean =   4.773 MeV   ( 4.1747 MeV --> 5.3712 MeV)
>           Os185:       1  Emean =  398.31 keV   ( 398.31 keV --> 398.31 keV)
>           Pb188:       1  Emean =  2.1818 MeV   ( 2.1818 MeV --> 2.1818 MeV)
>           Pb190:       1  Emean =  735.12 keV   ( 735.12 keV --> 735.12 keV)
>           Pb191:       2  Emean =  2.4688 MeV   ( 1.6298 MeV --> 3.3077 MeV)
>           Pb192:       5  Emean =  1.0334 MeV   ( 124.37 keV --> 3.1213 MeV)
>           Pb193:       5  Emean =  2.4228 MeV   ( 1.1147 MeV --> 4.7738 MeV)
>           Pb194:       7  Emean =  1.6385 MeV   ( 498.16 keV --> 2.2943 MeV)
>           Pb195:       8  Emean =   1.916 MeV   ( 192.07 keV --> 2.8607 MeV)
>           Pb196:      10  Emean =  787.91 keV   ( 36.979 keV --> 1.5127 MeV)
>           Pb197:       7  Emean =  1.4064 MeV   ( 262.9 keV --> 6.1721 MeV)
>           Pb198:       9  Emean =   934.6 keV   ( 283.42 keV --> 2.5077 MeV)
>           Pb199:       9  Emean =  1.0646 MeV   ( 264.57 keV --> 2.6094 MeV)
>           Pb200:      13  Emean =  955.35 keV   ( 85.273 keV --> 2.2468 MeV)
>           Pb201:       9  Emean =  1.2695 MeV   ( 197.11 keV --> 2.7812 MeV)
>           Pb202:       9  Emean =  867.93 keV   ( 158.27 keV --> 3.3965 MeV)
>           Pb203:      14  Emean =  508.75 keV   ( 50.536 keV --> 2.0848 MeV)
>           Pb204:      44  Emean =   265.3 keV   ( 176.8 eV  --> 1.5147 MeV)
>           Pb205:      28  Emean =  257.71 keV   ( 37.145 keV --> 980.33 keV)
>    Pb205[2.329]:       7  Emean =  327.01 keV   ( 119.86 keV --> 757.43 keV)
>           Pb206:     335  Emean =  25.674 keV   ( 51.866 eV  --> 1.181 MeV)
>           Pb207:     286  Emean =  31.462 keV   ( 15.756 eV  --> 1.1856 MeV)
>   Pb207[133.977]:       1  Emean =  102.95 keV  ( 102.95 keV --> 102.95 keV)
>   Pb207[136.844]:       1  Emean =  89.265 keV  ( 89.265 keV --> 89.265 keV)
>   Pb207[95.496]:       1  Emean =  368.26 keV   ( 368.26 keV --> 368.26 keV)
>           Pb208:     628  Emean =  14.279 keV   ( 6.5008 eV  --> 793.44 keV)
>   Pb208[1087.792]:       1  Emean =  259.72 keV ( 259.72 keV --> 259.72 keV)
>           Pd102:       1  Emean =  35.437 MeV   ( 35.437 MeV --> 35.437 MeV)
>           Pd104:       1  Emean =  64.997 MeV   ( 64.997 MeV --> 64.997 MeV)
>           Pd107:       1  Emean =  68.322 MeV   ( 68.322 MeV --> 68.322 MeV)
>           Pd110:       1  Emean =  50.322 MeV   ( 50.322 MeV --> 50.322 MeV)
>           Po198:       1  Emean =  2.6278 MeV   ( 2.6278 MeV --> 2.6278 MeV)
>           Po202:       1  Emean =  793.72 keV   ( 793.72 keV --> 793.72 keV)
>           Pt179:       1  Emean =  2.9353 MeV   ( 2.9353 MeV --> 2.9353 MeV)
>           Pt180:       2  Emean =  6.4751 MeV   ( 3.3984 MeV --> 9.5518 MeV)
>           Pt181:       2  Emean =  4.5481 MeV   ( 3.0388 MeV --> 6.0574 MeV)
>           Pt182:       7  Emean =  4.6463 MeV   ( 730.76 keV --> 11.211 MeV)
>           Pt183:      11  Emean =  2.4979 MeV   ( 471.58 keV --> 6.4522 MeV)
>           Pt184:       8  Emean =  2.8524 MeV   ( 360.05 keV --> 9.2713 MeV)
>           Pt185:       5  Emean =  3.2371 MeV   ( 1.3506 MeV --> 5.0815 MeV)
>           Pt186:       3  Emean =  3.3565 MeV   ( 2.2569 MeV --> 5.1326 MeV)
>           Pt187:       4  Emean =  4.5095 MeV   ( 2.6896 MeV --> 6.6558 MeV)
>           Pt188:       2  Emean =  1.9525 MeV   ( 1.1919 MeV --> 2.7132 MeV)
>           Pt189:       2  Emean =  3.0549 MeV   ( 1.2791 MeV --> 4.8308 MeV)
>           Pt191:       1  Emean =  58.934 keV   ( 58.934 keV --> 58.934 keV)
>           Pt192:       1  Emean =  826.19 keV   ( 826.19 keV --> 826.19 keV)
>           Pt195:       1  Emean =  491.56 keV   ( 491.56 keV --> 491.56 keV)
>            Rb82:       1  Emean =  67.278 MeV   ( 67.278 MeV --> 67.278 MeV)
>            Rb84:       1  Emean =  51.537 MeV   ( 51.537 MeV --> 51.537 MeV)
>            Rb85:       1  Emean =  76.358 MeV   ( 76.358 MeV --> 76.358 MeV)
>            Rb86:       1  Emean =  58.174 MeV   ( 58.174 MeV --> 58.174 MeV)
>            Rb87:       1  Emean =  84.035 MeV   ( 84.035 MeV --> 84.035 MeV)
>           Re173:       1  Emean =  1.1281 MeV   ( 1.1281 MeV --> 1.1281 MeV)
>           Re174:       1  Emean =  738.15 keV   ( 738.15 keV --> 738.15 keV)
>           Re175:       3  Emean =  5.1554 MeV   ( 734.78 keV --> 8.2775 MeV)
>           Re176:       1  Emean =  970.88 keV   ( 970.88 keV --> 970.88 keV)
>           Re177:       2  Emean =   3.557 MeV   ( 432.2 keV --> 6.6817 MeV)
>           Re178:       1  Emean =   2.442 MeV   ( 2.442 MeV --> 2.442 MeV)
>           Re179:       1  Emean =  5.0514 MeV   ( 5.0514 MeV --> 5.0514 MeV)
>           Re180:       1  Emean =  2.1975 MeV   ( 2.1975 MeV --> 2.1975 MeV)
>           Re181:       2  Emean =  11.713 MeV   ( 11.212 MeV --> 12.214 MeV)
>           Re182:       1  Emean =  8.3965 MeV   ( 8.3965 MeV --> 8.3965 MeV)
>           Re184:       2  Emean =  3.8791 MeV   ( 3.2335 MeV --> 4.5246 MeV)
>           Re187:       1  Emean =  2.4867 MeV   ( 2.4867 MeV --> 2.4867 MeV)
>           Rh102:       1  Emean =  45.646 MeV   ( 45.646 MeV --> 45.646 MeV)
>           Rh103:       1  Emean =   40.33 MeV   ( 40.33 MeV --> 40.33 MeV)
>           Rh107:       1  Emean =  47.543 MeV   ( 47.543 MeV --> 47.543 MeV)
>           Ru100:       2  Emean =  55.398 MeV   ( 53.422 MeV --> 57.374 MeV)
>           Ru101:       1  Emean =  79.006 MeV   ( 79.006 MeV --> 79.006 MeV)
>           Ru102:       1  Emean =  76.102 MeV   ( 76.102 MeV --> 76.102 MeV)
>            Ru99:       1  Emean =  85.573 MeV   ( 85.573 MeV --> 85.573 MeV)
>             S34:       1  Emean =  117.35 MeV   ( 117.35 MeV --> 117.35 MeV)
>           Sb118:       2  Emean =  30.321 MeV   ( 27.16 MeV --> 33.482 MeV)
>            Se76:       1  Emean =  69.124 MeV   ( 69.124 MeV --> 69.124 MeV)
>            Se77:       1  Emean =  90.462 MeV   ( 90.462 MeV --> 90.462 MeV)
>            Se80:       1  Emean =  90.173 MeV   ( 90.173 MeV --> 90.173 MeV)
>           Sn115:       1  Emean =  61.473 MeV   ( 61.473 MeV --> 61.473 MeV)
>           Sn118:       1  Emean =  43.838 MeV   ( 43.838 MeV --> 43.838 MeV)
>            Sr87:       3  Emean =  77.573 MeV   ( 60.058 MeV --> 103.14 MeV)
>            Sr88:       1  Emean =  77.836 MeV   ( 77.836 MeV --> 77.836 MeV)
>            Sr91:       1  Emean =  102.45 MeV   ( 102.45 MeV --> 102.45 MeV)
>            Sr92:       1  Emean =  79.242 MeV   ( 79.242 MeV --> 79.242 MeV)
>           Ta166:       1  Emean =  4.7107 MeV   ( 4.7107 MeV --> 4.7107 MeV)
>           Ta170:       2  Emean =   2.026 MeV   ( 933.42 keV --> 3.1187 MeV)
>           Ta174:       1  Emean =  8.2524 MeV   ( 8.2524 MeV --> 8.2524 MeV)
>           Ta175:       1  Emean =  3.4302 MeV   ( 3.4302 MeV --> 3.4302 MeV)
>           Tc100:       1  Emean =  72.671 MeV   ( 72.671 MeV --> 72.671 MeV)
>           Tc102:       2  Emean =   38.66 MeV   ( 36.383 MeV --> 40.937 MeV)
>            Tc99:       1  Emean =  54.775 MeV   ( 54.775 MeV --> 54.775 MeV)
>            Ti47:       1  Emean =  97.674 MeV   ( 97.674 MeV --> 97.674 MeV)
>           Tl186:       1  Emean =  1.6104 MeV   ( 1.6104 MeV --> 1.6104 MeV)
>           Tl187:       1  Emean =  3.3722 MeV   ( 3.3722 MeV --> 3.3722 MeV)
>           Tl188:       1  Emean =  4.2498 MeV   ( 4.2498 MeV --> 4.2498 MeV)
>           Tl189:       2  Emean =  3.3274 MeV   ( 2.6445 MeV --> 4.0104 MeV)
>           Tl190:       1  Emean =   6.497 MeV   ( 6.497 MeV --> 6.497 MeV)
>           Tl191:       9  Emean =  1.4788 MeV   ( 347.69 keV --> 3.2514 MeV)
>           Tl192:       7  Emean =  1.5407 MeV   ( 966.36 keV --> 2.3841 MeV)
>           Tl193:       5  Emean =  740.84 keV   ( 403.97 keV --> 1.4271 MeV)
>           Tl194:       7  Emean =  1.2713 MeV   ( 289.4 keV --> 3.3155 MeV)
>           Tl195:       8  Emean =  1.2208 MeV   ( 53.185 keV --> 2.0412 MeV)
>           Tl196:       8  Emean =  1.4907 MeV   ( 230.42 keV --> 2.8933 MeV)
>           Tl197:       8  Emean =  2.2203 MeV   ( 701.5 keV --> 4.9371 MeV)
>           Tl198:       5  Emean =   1.556 MeV   ( 554.22 keV --> 3.312 MeV)
>           Tl199:      11  Emean =  945.02 keV   ( 56.503 keV --> 1.8689 MeV)
>           Tl200:       3  Emean =   1.516 MeV   ( 797.72 keV --> 2.1402 MeV)
>           Tl201:      10  Emean =  1.0317 MeV   ( 420.37 keV --> 2.0071 MeV)
>           Tl202:       2  Emean =  905.63 keV   ( 287.28 keV --> 1.524 MeV)
>           Tl203:       8  Emean =  1.0918 MeV   ( 117.02 keV --> 2.7433 MeV)
>           Tl204:       6  Emean =  803.32 keV   ( 115.05 keV --> 1.3793 MeV)
>           Tl205:      13  Emean =  553.07 keV   ( 4.3422 keV --> 1.8611 MeV)
>           Tl206:       8  Emean =  697.33 keV   ( 66.916 keV --> 1.3412 MeV)
>           Tl207:       8  Emean =   494.2 keV   ( 35.076 keV --> 1.2712 MeV)
>   Tl207[99.095]:       1  Emean =  516.92 keV   ( 516.92 keV --> 516.92 keV)
>           Tl208:       1  Emean =  712.39 keV   ( 712.39 keV --> 712.39 keV)
>             V50:       1  Emean =  93.971 MeV   ( 93.971 MeV --> 93.971 MeV)
>            W169:       2  Emean =   819.3 keV   ( 387 keV --> 1.2516 MeV)
>            W172:       2  Emean =  3.5805 MeV   ( 3.0541 MeV --> 4.1069 MeV)
>            W173:       1  Emean =    12.2 MeV   ( 12.2 MeV --> 12.2 MeV)
>            W174:       4  Emean =  3.8323 MeV   ( 1.5817 MeV --> 9.0271 MeV)
>            W175:       1  Emean =  1.4302 MeV   ( 1.4302 MeV --> 1.4302 MeV)
>           Xe122:       1  Emean =  11.974 MeV   ( 11.974 MeV --> 11.974 MeV)
>           Xe126:       1  Emean =   55.08 MeV   ( 55.08 MeV --> 55.08 MeV)
>             Y87:       1  Emean =  71.522 MeV   ( 71.522 MeV --> 71.522 MeV)
>             Y90:       2  Emean =  68.396 MeV   ( 55.795 MeV --> 80.996 MeV)
>             Y91:       5  Emean =  65.102 MeV   ( 45.854 MeV --> 94.735 MeV)
>             Y92:       1  Emean =   46.69 MeV   ( 46.69 MeV --> 46.69 MeV)
>             Y93:       1  Emean =  75.443 MeV   ( 75.443 MeV --> 75.443 MeV)
>           Yb162:       1  Emean =  4.2822 MeV   ( 4.2822 MeV --> 4.2822 MeV)
>           Yb164:       1  Emean =  3.4843 MeV   ( 3.4843 MeV --> 3.4843 MeV)
>           Yb166:       2  Emean =  9.4699 MeV   ( 8.8152 MeV --> 10.125 MeV)
>           Yb167:       1  Emean =  12.696 MeV   ( 12.696 MeV --> 12.696 MeV)
>            Zr90:       2  Emean =  55.654 MeV   ( 55.305 MeV --> 56.002 MeV)
>            Zr91:       1  Emean =  80.042 MeV   ( 80.042 MeV --> 80.042 MeV)
>            Zr92:       3  Emean =   77.04 MeV   ( 69.381 MeV --> 81.284 MeV)
>            Zr93:       2  Emean =  77.971 MeV   ( 48.61 MeV --> 107.33 MeV)
>           alpha:     630  Emean =  8.7785 MeV   ( 1.6173 MeV --> 92.698 MeV)
>       anti_nu_e:      99  Emean =   3.455 MeV   ( 16.193 keV --> 45.127 MeV)
>      anti_nu_mu:      13  Emean =   34.05 MeV   ( 27.724 MeV --> 51.361 MeV)
>        deuteron:     170  Emean =  26.242 MeV   ( 397.02 keV --> 76.006 MeV)
>              e+:     660  Emean =  14.636 MeV   ( 103.73 eV  --> 500.8 MeV)
>              e-:   13055  Emean =  1.0045 MeV   ( 2.255 eV  --> 478.36 MeV)
>             eta:       1  Emean =  71.285 MeV   ( 71.285 MeV --> 71.285 MeV)
>           gamma:    7610  Emean =  6.5691 MeV   ( 208.36 eV  --> 716.51 MeV)
>             mu+:       3  Emean =  4.1198 MeV   ( 4.1198 MeV --> 4.1198 MeV)
>             mu-:      11  Emean =  10.309 MeV   ( 4.1198 MeV --> 45.755 MeV)
>         neutron:    8035  Emean =   26.96 MeV   ( 1.1718 keV --> 971.47 MeV)
>            nu_e:    2118  Emean =  1.0872 MeV   ( 536.8 eV  --> 39.132 MeV)
>           nu_mu:      12  Emean =  31.336 MeV   ( 13.34 MeV --> 49.681 MeV)
>             pi+:     192  Emean =  141.79 MeV   ( 8.3805 MeV --> 724.92 MeV)
>             pi-:     125  Emean =  90.981 MeV   ( 406.61 keV --> 493.74 MeV)
>             pi0:     152  Emean =  139.31 MeV   ( 7.8327 MeV --> 598.12 MeV)
>          proton:    1263  Emean =  159.92 MeV   ( 364.44 keV --> 992.46 MeV)
>          triton:      55  Emean =  22.371 MeV   ( 150.28 keV --> 74.495 MeV)
> 
>  List of generated particles in detector:
>            As72:       2  Emean =  437.31 keV   ( 303.19 keV --> 571.43 keV)
>            As74:       1  Emean =  95.262 keV   ( 95.262 keV --> 95.262 keV)
>            Co56:       1  Emean =  307.03 keV   ( 307.03 keV --> 307.03 keV)
>            Co58:       1  Emean =  569.66 keV   ( 569.66 keV --> 569.66 keV)
>            Co59:       2  Emean =  3.5193 MeV   ( 3.3366 MeV --> 3.7021 MeV)
>            Co60:       1  Emean =  1.4135 MeV   ( 1.4135 MeV --> 1.4135 MeV)
>            Cr49:       1  Emean =   11.75 MeV   ( 11.75 MeV --> 11.75 MeV)
>            Cu60:       1  Emean =   352.1 keV   ( 352.1 keV --> 352.1 keV)
>            Cu61:       2  Emean =  5.9689 MeV   ( 4.5164 MeV --> 7.4215 MeV)
>            Cu62:       1  Emean =  1.1162 MeV   ( 1.1162 MeV --> 1.1162 MeV)
>            Cu63:       3  Emean =  1.7098 MeV   ( 1.2782 MeV --> 2.0463 MeV)
>            Cu64:       2  Emean =  2.6956 MeV   ( 1.951 MeV --> 3.4403 MeV)
>            Cu65:       4  Emean =  782.11 keV   ( 197.17 keV --> 1.9061 MeV)
>            Cu66:       3  Emean =  1.5395 MeV   ( 1.1418 MeV --> 1.972 MeV)
>            Ga66:       1  Emean =  910.93 keV   ( 910.93 keV --> 910.93 keV)
>            Ga67:       3  Emean =  1.9948 MeV   ( 807.22 keV --> 3.1097 MeV)
>            Ga68:       5  Emean =   631.6 keV   ( 54.49 keV --> 1.674 MeV)
>            Ga69:       9  Emean =  864.65 keV   ( 170.97 keV --> 1.8592 MeV)
>            Ga70:      12  Emean =  644.52 keV   ( 38.042 keV --> 1.5066 MeV)
>            Ga71:       8  Emean =  425.49 keV   ( 111.73 keV --> 944.73 keV)
>            Ga72:      10  Emean =  344.29 keV   ( 12.267 keV --> 913.26 keV)
>            Ga73:       4  Emean =  146.49 keV   ( 25.209 keV --> 301.79 keV)
>            Ga74:       1  Emean =  273.71 keV   ( 273.71 keV --> 273.71 keV)
>            Ge68:       3  Emean =  1.1881 MeV   ( 510.2 keV --> 1.9486 MeV)
>            Ge69:      18  Emean =  680.44 keV   ( 59.684 keV --> 2.9284 MeV)
>            Ge70:     697  Emean =  54.745 keV   ( 24.793 eV  --> 4.2886 MeV)
>            Ge71:      38  Emean =   484.1 keV   ( 8.4236 keV --> 2.1524 MeV)
>            Ge72:     868  Emean =  59.624 keV   ( 11.753 eV  --> 1.1921 MeV)
>            Ge73:     285  Emean =  72.376 keV   ( 15.326 eV  --> 1.5638 MeV)
>            Ge74:    1151  Emean =  40.306 keV   ( 6.3175 eV  --> 992.62 keV)
>            Ge75:      10  Emean =  340.27 keV   ( 2.7609 keV --> 1.6506 MeV)
>            Ge76:     257  Emean =  40.983 keV   ( 36.585 eV  --> 761.83 keV)
>            Ge77:       1  Emean =  658.81 eV    ( 658.81 eV  --> 658.81 eV )
>             He3:       3  Emean =    28.4 MeV   ( 9.2128 MeV --> 63.002 MeV)
>            Ni57:       1  Emean =  185.79 keV   ( 185.79 keV --> 185.79 keV)
>            Ni60:       2  Emean =  1.4061 MeV   ( 1.145 MeV --> 1.6672 MeV)
>            Ni61:       4  Emean =  2.9325 MeV   ( 1.6072 MeV --> 3.8197 MeV)
>            Ni62:       6  Emean =  2.4641 MeV   ( 595.37 keV --> 6.3399 MeV)
>            Ni64:       1  Emean =  4.7604 MeV   ( 4.7604 MeV --> 4.7604 MeV)
>            Zn62:       1  Emean =  1.5052 MeV   ( 1.5052 MeV --> 1.5052 MeV)
>            Zn63:       1  Emean =  533.74 keV   ( 533.74 keV --> 533.74 keV)
>            Zn64:       7  Emean =  1.1173 MeV   ( 383.94 keV --> 3.2059 MeV)
>            Zn65:       2  Emean =  2.1855 MeV   ( 1.8343 MeV --> 2.5367 MeV)
>            Zn66:      23  Emean =  1.2464 MeV   ( 176.18 keV --> 3.6293 MeV)
>            Zn67:      10  Emean =  1.1542 MeV   ( 281.77 keV --> 3.1521 MeV)
>            Zn68:       4  Emean =  2.2507 MeV   ( 943.69 keV --> 4.0866 MeV)
>            Zn69:       4  Emean =  843.94 keV   ( 367.5 keV --> 1.3913 MeV)
>            Zn70:       3  Emean =  1.0988 MeV   ( 45.314 keV --> 2.8362 MeV)
>            Zn71:       1  Emean =  317.02 keV   ( 317.02 keV --> 317.02 keV)
>           alpha:      63  Emean =  9.3174 MeV   ( 4.2484 MeV --> 27.709 MeV)
>       anti_nu_e:      32  Emean =   31.41 MeV   ( 12.183 MeV --> 48.573 MeV)
>      anti_nu_mu:      52  Emean =  33.206 MeV   ( 20.598 MeV --> 52.14 MeV)
>        deuteron:      29  Emean =  17.275 MeV   ( 624.66 keV --> 52.332 MeV)
>              e+:     175  Emean =  34.158 MeV   ( 105.78 keV --> 255.27 MeV)
>              e-:    5522  Emean =  1.5054 MeV   ( 111.79 eV  --> 390.45 MeV)
>           gamma:    2616  Emean =  1.6044 MeV   ( 13.271 keV --> 177.85 MeV)
>             mu+:      21  Emean =  4.1198 MeV   ( 4.1198 MeV --> 4.1198 MeV)
>             mu-:      30  Emean =  4.1198 MeV   ( 4.1198 MeV --> 4.1198 MeV)
>         neutron:    1546  Emean =  7.9857 MeV   ( 13.836 keV --> 346.94 MeV)
>            nu_e:      22  Emean =  26.434 MeV   ( 5.7158 MeV --> 47.176 MeV)
>           nu_mu:      53  Emean =  33.726 MeV   ( 13.328 MeV --> 52.19 MeV)
>             pi+:       4  Emean =  75.005 MeV   ( 31.618 MeV --> 106.37 MeV)
>             pi-:       6  Emean =     134 MeV   ( 7.5546 MeV --> 374.17 MeV)
>             pi0:       2  Emean =  73.865 MeV   ( 43.724 MeV --> 104.01 MeV)
>          proton:     183  Emean =  47.601 MeV   ( 1.8391 MeV --> 774.1 MeV)
>          triton:       3  Emean =  12.304 MeV   ( 9.3208 MeV --> 15.114 MeV)
> 
> --------- Ranecu engine status ---------
>  Initial seed (index) = 0
>  Current couple of seeds = 1982454223, 1571101230
> ----------------------------------------
> G4WT1 > ================== Deleting memory pools ===================
> G4WT1 > Number of memory pools allocated: 11; of which, static: 0
> G4WT1 > Dynamic pools deleted: 11 / Total memory freed: 0.21 MB
> G4WT1 > ============================================================
> G4WT1 > Thread-local UImanager is to be deleted.
> G4WT1 > There should not be any thread-local G4cout/G4cerr hereafter.
> G4WT0 > ================== Deleting memory pools ===================
> G4WT0 > Number of memory pools allocated: 11; of which, static: 0
> G4WT0 > Dynamic pools deleted: 11 / Total memory freed: 0.28 MB
> G4WT0 > ============================================================
> G4WT0 > Thread-local UImanager is to be deleted.
> G4WT0 > There should not be any thread-local G4cout/G4cerr hereafter.
> G4WT2 > ================== Deleting memory pools ===================
> G4WT2 > Number of memory pools allocated: 11; of which, static: 0
> G4WT2 > Dynamic pools deleted: 11 / Total memory freed: 0.25 MB
> G4WT2 > ============================================================
> G4WT2 > Thread-local UImanager is to be deleted.
> G4WT2 > There should not be any thread-local G4cout/G4cerr hereafter.
> G4WT3 > ================== Deleting memory pools ===================
> G4WT3 > Number of memory pools allocated: 11; of which, static: 0
> G4WT3 > Dynamic pools deleted: 11 / Total memory freed: 0.27 MB
> G4WT3 > ============================================================
> G4WT3 > Thread-local UImanager is to be deleted.
> G4WT3 > There should not be any thread-local G4cout/G4cerr hereafter.
> G4 kernel has come to Quit state.
> ================== Deleting memory pools ===================
> Number of memory pools allocated: 4; of which, static: 0
> Dynamic pools deleted: 4 / Total memory freed: 0.0058 MB
> ============================================================
> RunManagerKernel is deleted. Good bye :)
> Cristians-iMac:rdecay02-build cristian$ ls -lrt
> total 792
> -rw-r--r--   1 cristian  staff    2578  6 Feb 18:28 vis.mac
> -rw-r--r--   1 cristian  staff     360  6 Feb 18:28 rdecay02.in
> -rw-r--r--   1 cristian  staff     694  6 Feb 18:28 plotHisto.C
> -rw-r--r--   1 cristian  staff     217  6 Feb 18:28 debug.mac
> -rw-r--r--   1 cristian  staff    2177  6 Feb 18:28 cmake_install.cmake
> -rw-r--r--   1 cristian  staff   19718  6 Feb 18:28 Makefile
> -rw-r--r--   1 cristian  staff   14011  6 Feb 18:28 CMakeCache.txt
> -rw-r--r--@  1 cristian  staff     793  6 Feb 18:31 proton-beam.mac
> -rw-r--r--@  1 cristian  staff      10  6 Feb 18:31 beam.data
> -rw-r--r--@  1 cristian  staff      52  6 Feb 18:31 measures.data
> -rw-r--r--@  1 cristian  staff     656  6 Feb 18:33 proton.mac
> -rwxr-xr-x   1 cristian  staff  323964  6 Feb 18:37 rdecay02
> drwxr-xr-x  15 cristian  staff     510  6 Feb 18:39 CMakeFiles
> -rw-r--r--   1 cristian  staff     388  6 Feb 18:39 proton-beam.activity
> Cristians-iMac:rdecay02-build cristian$ more proton-beam.activity 
> Radioactivities in decay window no. 0
> Z       A       E       Activity (decays/window)        Error (decays/window) 
> 
> Radioactivities in decay window no. 1 Z A E Activity (decays/window)
> Error (decays/window)
> 
> Radioactivities in decay window no. 2 Z A E Activity (decays/window)
> Error (decays/window)
> 
> Radioactivities in decay window no. 3 Z A E Activity (decays/window)
> Error (decays/window)
> 
> Also I fail to understand (sorry) the file
> 
> $ more measures.data
> 
> 0. 1
> 
> 1. 0
> 
> 1000 1
> 
> 1001 0
> 
> 10000 1
> 
> 10001 0
> 
> 1e7 1
> 
> 1e9 0
> 
> Is this suppose to mean the isotopes will decay in 4 times window (0 - 1
> seconds) (1000 - 1001 seconds) (10000 - 10001 seconds) and (1e7-1e9
> seconds) ? Are the units really seconds?
> 
> Thank you,
> 
> Cristian
> 

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

1 More: Re: ./rdecay02 proton-beam.mac output file   (Cristian Bungau - 24 Feb, 2016)
(_ None: Re: ./rdecay02 proton-beam.mac output file   (michel maire - 24 Feb, 2016)
(_ More: Re: ./rdecay02 proton-beam.mac output file   (Cristian Bungau - 24 Feb, 2016)
(_ Idea: Re: ./rdecay02 proton-beam.mac output file   (Michael H. Kelsey - 25 Feb, 2016)
(_ More: Re: ./rdecay02 proton-beam.mac output file   (Cristian Bungau - 26 Feb, 2016)
(_ Idea: Re: ./rdecay02 proton-beam.mac output file   (Michael H. Kelsey - 26 Feb, 2016)
(_ More: Re: ./rdecay02 proton-beam.mac output file   (Cristian Bungau - 26 Feb, 2016)
(_ Ok: Re: ./rdecay02 proton-beam.mac output file   (Michael H. Kelsey - 26 Feb, 2016)
(_ Ok: Re: ./rdecay02 proton-beam.mac output file   (Cristian Bungau - 26 Feb, 2016)
 Add Message Add Message
to: "Re: ./rdecay02 proton-beam.mac output file"

 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 ]