Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

MCNP VisEd Plot tracks :Error->File xsdir not available

  1. Jul 28, 2017 #1
    Hello,
    I use free VisEd and I want to plot the collision, the particle's transport and the particle in tally but I have an error. For the source the plot is ok. How can I solve this problem freely?
    Someone can help me?

    My input :


    C Cellules


    30 50 -1 (-1 3 19 ):(-2 3 ):(-4 3 -6 5 -7 1 2 ) imp:p=1

    31 50 -1 -19 3 #(-4 3 -6 5 -7 1 2 ) imp:p=1 $Nodule eau

    32 51 -1.19 -8 -9 #((-1 3 ):(-2 3 ):(-4 3 -6 5 -7 1 2 ): $cou PMMA

    (-19 3 )) imp:p=1

    33 52 -2.7 -10 11 imp:p=1

    34 53 -0.0012 (-11 -13 ):-16 :-18 imp:p=1

    35 51 -1.19 -11 13 -14 12 imp:p=1

    36 54 -3.58 -11 14 imp:p=1

    37 55 -3.67 -12 imp:p=1

    38 56 -11.35 -15 16 imp:p=1

    39 57 -7.87 -17 18 imp:p=1

    40 53 -0.0012 -20 #(-8 -9 )10 15 17 imp:p=1

    80 0 20 imp:p=0

    c Structures

    c Thyroide

    1 sq 0.982 0.982 0.171 0 0 0 -1 1.4122 0 0

    2 sq 0.982 0.982 0.171 0 0 0 -1 -1.4122 0 0

    3 cz 1.009

    4 cz 1.412

    5 pz -2.017

    6 pz 0

    7 py 0

    8 rcc 0 3.528 -3 0 0 6 5.4

    9 py 1.412

    c Detecteur NaI

    10 41 rcc 0 0 0 0 3 0 2.85

    11 41 rcc 0 0.04 0 0 2.92 0 2.81

    12 41 rcc 0 0.18 0 0 2.54 0 2.54

    13 41 py 0.18

    14 41 py 2.72

    15 41 rcc 0 3 0 0 1 0 2.85

    16 41 rcc 0 3 0 0 1 0 2.35

    17 41 rcc 0 4 0 0 0.32 0 2.85

    18 41 rcc 0 4 0 0 0.32 0 2.64

    19 sx 1.77 0.6

    20 so 80

    c tr41 0 -21.192 0 $DDC=15cm

    c tr41 0 -36.192 0 $DDC=30cm

    mode p

    c eau

    m50 1000. 2

    8000. 1

    c PMMA

    m51 1000. -0.080541

    6000. -0.599846 8000. -0.319613

    c Aluminium

    m52 13000. 1

    c Air

    m53 6000. -0.000124

    7000. -0.755268 8000. -0.231781 18000. -0.012827

    c MgO

    m54 12000. -0.603036

    8000. -0.396964

    c NaI

    m55 11000. -0.153373

    53000. -0.846627

    c Plomb

    m56 82000. 1

    c Inox (Fer)

    m57 26000. 1

    c Reste

    c deplacement detecteur

    tr41 0 -6.192 0

    c Definition des paramètres uptake

    sdef PAR=2 CEL=d1 POS=FCEL d2 AXS=0 0 1 RAD=FCEL d5 EXT=FCEL d8 ERG=FCEL d20

    c Choix des cellules et leur uptake

    si1 L 30 31

    sp1 0 1

    c SP1 0.1 0.9

    c SP1 0.2 0.8

    c SP1 0.3 0.7

    c SP1 0.5 0.5

    c SP1 0.7 0.3

    c SP1 1 0

    c Position centrale du cylindre englobant chaque source

    ds2 S 3 4

    si3 L 0 0 0

    sp3 1

    si4 L 1.77 0 0

    sp4 1

    c Rayon du cylindre englobant chaque source

    ds5 S 6 7

    si6 0 2.5

    sp6 -21 1

    si7 0 0.7

    sp7 -21 1

    c Limite inférieure et supérieure de chaque cylindre

    ds8 S 9 10

    si9 -3 3

    sp9 -21 0

    si10 -0.7 0.7

    sp10 -21 0

    c Energie des photons de chaque source

    ds20 S 21 22

    si21 L 0.356

    sp21 D 1

    si22 L 0.356

    sp22 D 1

    f8:p 37

    e8 0.0 1e-4 0.320 100i 0.370

    nps 1e6


    Thank you.
     
  2. jcsd
  3. Jul 28, 2017 #2
    Hi,
    I am not sure to understand your problem. If you want to see the interaction and the particles transport you must have a "ptrac"
    PSR
     
  4. Aug 1, 2017 #3
    If you use Vised, it's possible to display all particles in source. I am trying to display particle transport from source to target (NaI) and It's not possible because "file xsdir" is not available. Do you have a solution to solve this error? For the precision, I use a free version so I need a free solution.
    Thank you for your help.
     
  5. Aug 1, 2017 #4
    Hi,

    I am surprised because xsdir is essential for MCNP. Can you run your file with no problem ?
    I am ok that you can see the source emission with vised, but if you want to "see" the particles transport and the interaction from the source to your INa detector you must put something like:
    PTRAC write=all file=asc max=1000 cell=37
    I tried your file with this line and moritz and it seams to be ok,
    PSR
     
  6. Aug 1, 2017 #5
    pfXY3D.jpg


    With your file and the PTRAC line It gives this type of image
     
  7. Aug 2, 2017 #6
    Hi,
    thank you for your answer but how do you display this image. Do you use MCNP or VisEd? I have tried in MCNP5 by adding your line at the end of my input (last part of code) and I have just the output ptrac. With VisEd when I display the source I obtain this image.

    Thank you.

    PO
     

    Attached Files:

  8. Aug 2, 2017 #7
    Hi,
    yes, first you run MCNP with the line PTRAC I give you. A file "ptrac" is created. I read this file with "moritz". I think that it is possible to read a ptrac file with "vised". You must have something like :
    read a ptrac file,
    read a track file,
    ...
    I don't know vised but with moritz it is:
    file-> Read-> read particle tracks
    hope this helps
    PSR
     
  9. Aug 3, 2017 #8
    Hi, with mcnp5 I obtained these files. Try with my ptrac file to display your image. It's not ok in VisEd. Moritz is a free software?

    Thanks
     

    Attached Files:

  10. Aug 3, 2017 #9
    Hi,
    it gives that (ptrad.txt), with Moritz (not free)
    but I think it is possible with vised
     

    Attached Files:

  11. Aug 4, 2017 #10
    Hi,
    thank you for your answer but your plot shows that the source is out off the nodule whereas the source is only on the nodule. In the distribution, we have 100% in the nodule and nothing out off (see TestPlot). Why the particle doesn't exist in y direction with your plot? If you change the ration uptake in nodule (70%) and 30% in out off (see MCNPCodeD70), how the transport become?

    I can't buy Moritz.

    Thank you.
     

    Attached Files:

  12. Aug 4, 2017 #11
    ok, I understand my mistake,
    I wrote
    SI1 L 30 31
    SP1 0.2 0.8
    instead of
    SI1 L 30 31
    SP1 0 1
    For y direction it's because I calculated only the tracks with a contribution in the tally

    "with a contribution in the tally" it comes:

    upload_2017-8-4_22-44-55.png


    and for all the tracks it comes :
    upload_2017-8-4_22-53-48.png

    that's all folks
    PSR
     
  13. Aug 7, 2017 #12
    Hi,
    I have to explain why the counting efficiency decreases when the ratio of uptake in the nodule increases. So can you show me for all tracks these images :
    -70% at contact (script1)
    - 100% in nodule at 15 cm (script2) and 70% at 15 cm (script3)
    - 100% in nodule at 30 cm (script4) and 70% at 30 cm (script5)

    I am sorry to ask you many questions.
    Thank you your help.
     

    Attached Files:

  14. Aug 7, 2017 #13
    Sorry,
    But I don't have time to run 5 MCNP's files
    PSR
     
  15. Aug 9, 2017 #14
    Hi,
    I know you have no time but please it's my last questions before to close this post. Sorry again for all my questions.


    PS:"You have already two beers in Paris. When you want!"
     
  16. Aug 9, 2017 #15
    My results,
    script1
    upload_2017-8-9_23-18-27.png

    script 2
    upload_2017-8-9_23-18-55.png

    script 3
    upload_2017-8-9_23-19-20.png

    script 4
    upload_2017-8-9_23-19-55.png

    script 5

    upload_2017-8-9_23-20-26.png

    So as you can say "Champagne au Moulin rouge" ...
     
  17. Aug 10, 2017 #16
    Thank you very match. No problem for your Champagne. Send me a message if you are in Paris.


    Best regards.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted



Similar Discussions: MCNP VisEd Plot tracks :Error->File xsdir not available
  1. Batch file (Replies: 2)

  2. Deleting these files? (Replies: 1)

Loading...