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

[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 Cellules30 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 1e6Thank you.

 

[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

 

[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.
Thank you for your help.

 

[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

 

[PSRB191921]

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

 

[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

 

[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 :
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

 

[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

 

[PSRB191921]

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

 

[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?

I can't buy Moritz.

Thank you.

 

[PSRB191921]

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:

and for all the tracks it comes :

that's all folks
PSR

 

[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.
Thank you your help.

 

[PSRB191921]

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

 

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

 

[PSRB191921]

My results,
script1

script 2

script 3

script 4

script 5

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

 

[DePaco]

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