Error importing ICRU-145

[RafaelOScatena]

Hello,

I get an error while loading ICRU-145 mesh phantom, independently of my geometry input block. I made the pegs4 file with the materials...

Thank you in advance,

Rafael Scatena

error log

RAYLEIGH OPTION REQUESTED FOR MEDIUM NUMBER187

END OF FILE ON UNIT 12

 PROGRAM STOPPED IN HATCH BECAUSE THE
FOLLOWING NAMES WERE NOT RECOGNIZED:

                                        100
                                       200                     
                                        300                    
                                        301                    

:start` geometry definition:
    :start geometry:
        name        = my_mesh
        library     = egs_mesh
        # or, using the TetGen format:
         file       = MRCP_AF.ele
        hownear method  = 1
        :start media input:
         media = icru145-female.pegs4dat

         set medium = 100 1
         set medium = 200 1
         set medium = 300 2
         set medium = 301 2
         set medium = 302 2
         set medium = 303 2
         set medium = 400 2
         set medium = 401 2
         set medium = 402 2
         set medium = 403 2
         set medium = 404 2
         set medium = 405 2
         set medium = 500 3
         set medium = 501 3
         set medium = 600 3
         set medium = 700 2
         set medium = 800 2
         set medium = 801 2
         set medium = 802 2
         set medium = 803 2
         set medium = 804 2
         set medium = 805 2
         set medium = 806 2
         set medium = 807 2
         set medium = 808 2
         set medium = 813 2
         set medium = 900 4
         set medium = 910 4
         set medium = 1000 4
         set medium = 1010 4
         set medium = 1100 4
         set medium = 1110 4
         set medium = 1200 4
         set medium = 1210 4
         set medium = 1300 5
         set medium = 1400 7
         set medium = 1500 6
         set medium = 1600 5
         set medium = 1700 8
         set medium = 1800 6
         set medium = 1900 5
         set medium = 2000 8
         set medium = 2100 6
         set medium = 2200 8
         set medium = 2300 8
         set medium = 2400 5
         set medium = 2500 9
         set medium = 2600 5
         set medium = 2700 10
         set medium = 2800 5
         set medium = 2900 11
         set medium = 3000 6
         set medium = 3100 5
         set medium = 3200 11
         set medium = 3300 6
         set medium = 3400 5
         set medium = 3500 8
         set medium = 3600 8
         set medium = 3700 5
         set medium = 3800 8
         set medium = 3900 5
         set medium = 4000 12
         set medium = 4100 5
         set medium = 4200 13
         set medium = 4300 5
         set medium = 4400 14
        set medium = 4500 5
         set medium = 4600 15
         set medium = 4700 5
         set medium = 4800 16
         set medium = 4900 5
         set medium = 5000 17
         set medium = 5100 5
         set medium = 5200 18
         set medium = 5300 5
         set medium = 5400 19
         set medium = 5500 5
         set medium = 5600 20
         set medium = 5700 21
         set medium = 5800 21
         set medium = 6100 22
         set medium = 6200 23
         set medium = 6300 24
         set medium = 6400 23
         set medium = 6500 24
         set medium = 6600 25
         set medium = 6601 25
         set medium = 6700 26
         set medium = 6701 27
         set medium = 6702 28
         set medium = 6800 25
         set medium = 6801 25
         set medium = 6900 26
         set medium = 6901 27
         set medium = 6902 28
         set medium = 7000 2
         set medium = 7100 29
         set medium = 7200 30
         set medium = 7201 30
         set medium = 7202 30
         set medium = 7203 30
         set medium = 7300 33
         set medium = 7400 31
         set medium = 7401 31
         set medium = 7402 31
         set medium = 7403 31
         set medium = 7500 31
         set medium = 7501 31
        set medium = 7600 32
         set medium = 7601 32
         set medium = 7602 32
         set medium = 7700 32
        set medium = 7800 32
         set medium = 7801 32
         set medium = 7802 32
        set medium = 7900 32
         set medium = 8000 32
         set medium = 8001 32
        set medium = 8002 32
         set medium = 8100 32
         set medium = 8200 32
         set medium = 8201 32
         set medium = 8202 32
         set medium = 8300 32
        set medium = 8400 32
         set medium = 8401 32
         set medium = 8402 32
         set medium = 8500 32
         set medium = 8600 32
         set medium = 8700 34
         set medium = 8800 4
         set medium = 8900 35
         set medium = 9000 35
         set medium = 9100 35
        set medium = 9200 35
         set medium = 9300 35
         set medium = 9400 35
         set medium = 9500 36
         set medium = 9700 37
         set medium = 9900 37
         set medium = 10000 38
         set medium = 10100 38
         set medium = 10200 38
         set medium = 10300 38
         set medium = 10400 38
         set medium = 10500 38
         set medium = 10600 39
         set medium = 10700 39
         set medium = 10800 39
         set medium = 10900 39
         set medium = 11000 40
        set medium = 11001 40
         set medium = 11002 40
         set medium = 11003 40
         set medium = 11100 41
         set medium = 11200 41
         set medium = 11300 42
         set medium = 11400 2
         set medium = 11500 43
         set medium = 11600 44
         set medium = 11700 44
         set medium = 11800 44
         set medium = 11900 44
         set medium = 12000 2
         set medium = 12100 2
         set medium = 12200 45
         set medium = 12201 45
         set medium = 12300 45
         set medium = 12301 45
        set medium = 12400 45
         set medium = 12401 45
        set medium = 12500 45
         set medium = 12501 45
        set medium = 12600 2
         set medium = 12700 46
         set medium = 12800 47
         set medium = 12801 47
         set medium = 12900 41
         set medium = 13000 41
         set medium = 13100 2
         set medium = 13200 48
         set medium = 13300 3
         set medium = 13301 3
        set medium = 13400 2
         set medium = 13500 2
         set medium = 13600 2
         set medium = 13700 49
         set medium = 13701 49
         set medium = 13800 50
         set medium = 13900 43
         set medium = 14000 51 
  
       :stop media input:
    :stop geometry:

    simulation geometry = my_mesh
:stop geometry definition:

[rtownson]

If you send me an email I can send you a zip file with egs++ input files ready for use with the MCRP phantoms, along with the materials as .density files, and instructions. This will be included as supplementary material in the 2026 release.

[cviniciusgomesf]

Hi Reid,

I'm trying to follow the mesh phantom instructions. Initially, when I navigated to HEN_HOUSE/egs++ and ran

git checkout add-mesh-dose-scoring

I received:
fatal: not a git repository (or any of the parent directories): .git

I realized my EGSnrc installation was not a Git clone, so I performed a fresh clone from the GitHub repository. After cloning, I was able to successfully check out the develop branch. However, when I tried again with

git checkout add-mesh-dose-scoring

I received:
error: pathspec 'add-mesh-dose-scoring' did not match any file(s) known to git

and git branch -a does not list that branch. Could you please advise on the recommended installation procedure for the mesh phantom examples with the current EGSnrc version? Is the develop branch sufficient, or are there additional steps, branches, or materials that I should use?

Best regards,
Carlos

[rtownson]

Hi Carlos,

I just answered your email, but just to have it here too - the add-mesh-dose-scoring branch doesn't exist because it was merged into the develop branch. So it's ok to just use that, and soon that will all end up in the official release.

To increase the number of materials, adjust $MXMED in HEN_HOUSE/src/egsnrc.macros (line 182), and then recompile whichever application you're using.

[rtownson]

I added the zip file with MCRP materials to the 2025 release assets:

https://github.com/nrc-cnrc/EGSnrc/releases/tag/v2025

Note that it suggests to use git checkout add-mesh-dose-scoring in the instructions. You may wish to use git checkout develop instead to be on the most recent version (before running the installer).

[lpfismed]

Hi @rtownson , thank you for sharing these files. I am trying to run a simulation. First, I changed the maximum number of media in the array_sizes.h file, then I recompiled the usercode. Now, I am getting this message:

EGSnrc SUCCESSFULLY 'HATCHED' FOR 187 MEDIA.

Error entering get_inputs: Asked for values from 187 to 1 with a max of 100

This implies a bug in the calling routine
Fix it up and try again. Stopping now.

I tried to change REPLACE {$NVALUE} WITH {100}; for 300 in HEN_HOUSE/src/egsnrc.macros, but it didn't work. Can you help me?
Cheers!
Lucas