Dear experts,

I am dealing with a problem where alpha and triton particles released from a LiF target hit by thermal neutrons. The detector is a box 3 x 3 x 1 cm (1cm is the drift region).

With the EVENTBIN scoring estimator I could get the energy deposition as follows:

EVENTBIN -12.0 ENERGY -97. GAS TOT

EVENTBIN 1. 0.0 0.0 &

with a region binning, which is equivalent with the following option with Cartesian binning:

EVENTBIN -10.0 ENERGY -46. 3.0 3.0 1.000035GASREG5

EVENTBIN -3.0 -3.0 3.50E-5 1. 1. 1. &

The results are the same as expected.

In addition, I tried to calculate the energy deposition using mgdraw.f routine:

From the main entry I record continuous energy deposition in the detector region

```
IF (MREG.EQ.5) THEN
IF (MTRACK.GT.0 .AND. (JTRACK.EQ.-4 .OR. JTRACK.EQ.-6)) THEN
WRITE(67,*) NCASE,JTRACK,SNGL(EKIN),SNGL(CTRACK),
+ SNGL(XTRACK(1:NTRACK)),SNGL(YTRACK(1:NTRACK)),
+ SNGL(ZTRACK(1:NTRACK)),SNGL(CXTRCK),SNGL(CYTRCK),
+ SNGL(CZTRCK),SNGL(ANGLE2),SNGL(DTRACK(1:MTRACK))
```

If I sum the SNGL(DTRACK(1:MTRACK) per primary (or use the entry ENDRAW, which is equivalent since RULL = 0 in this case -local energy depositions-) I am getting the same energy deposition spectra per event as using the two EVENTBIN cards above.

My question is the following: Is it possible to get in the main mgdraw entry the Cartesian X,Y,Z inside in the same NTUPLE with Edep/hit and XTRACK, YTRACK, ZTRACK?

In other words, how this EVENTBIN card:

EVENTBIN -10.0 ENERGY -46. 3.0 3.0 1.000035GASREG5

EVENTBIN -3.0 -3.0 3.50E-5 1. 1. 1. &

can be calculated inside mgdraw?

I would like at the end to have a complete NTUPLE built inside mgdraw where I could apply conditions not only per track lenght (position pointed by X,Y,Ztrack) but also per X,Y,Zvolume inside the detector region.

Thanks a lot in advance

I am looking forward to hearing from you

Georgios