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

Tags:
1. Jul 28, 2017

### DePaco

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

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. Jul 28, 2017

### PSRB191921

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

3. Aug 1, 2017

### DePaco

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.

4. Aug 1, 2017

### PSRB191921

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

5. Aug 1, 2017

### PSRB191921

With your file and the PTRAC line It gives this type of image

6. Aug 2, 2017

### DePaco

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

File size:
48.8 KB
Views:
51
7. Aug 2, 2017

### PSRB191921

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 :
...
I don't know vised but with moritz it is:
hope this helps
PSR

8. Aug 3, 2017

### DePaco

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

File size:
817 bytes
Views:
20
File size:
190.9 KB
Views:
21
9. Aug 3, 2017

### PSRB191921

Hi,
it gives that (ptrad.txt), with Moritz (not free)
but I think it is possible with vised

#### Attached Files:

File size:
29.4 KB
Views:
35
• ###### pf.pdf
File size:
98.9 KB
Views:
33
10. Aug 4, 2017

### DePaco

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?

Thank you.

#### Attached Files:

File size:
2.8 KB
Views:
28
File size:
45.7 KB
Views:
41
• ###### MCNPCodeD70.txt
File size:
2.8 KB
Views:
22
11. Aug 4, 2017

### PSRB191921

ok, I understand my mistake,
I wrote
SI1 L 30 31
SP1 0.2 0.8
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:

and for all the tracks it comes :

that's all folks
PSR

12. Aug 7, 2017

### DePaco

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.

#### Attached Files:

File size:
2.8 KB
Views:
31
File size:
2.8 KB
Views:
28
File size:
2.8 KB
Views:
23
File size:
2.8 KB
Views:
20
• ###### Script5.txt
File size:
2.8 KB
Views:
23
13. Aug 7, 2017

### PSRB191921

Sorry,
But I don't have time to run 5 MCNP's files
PSR

14. Aug 9, 2017

### DePaco

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!"

15. Aug 9, 2017

### PSRB191921

My results,
script1

script 2

script 3

script 4

script 5

So as you can say "Champagne au Moulin rouge" ...

16. Aug 10, 2017

### DePaco

Thank you very match. No problem for your Champagne. Send me a message if you are in Paris.

Best regards.