**From:** Alberto Fasso' (*fasso@SLAC.Stanford.EDU*)

**Date:** Fri May 26 2006 - 20:03:02 CEST

**Previous message:**Dipanwita Dutta: "RE: LET from USRYIELD: units"**In reply to:**Dipanwita Dutta: "RE: LET from USRYIELD: units"**Next in thread:**Reddell, Brandon D: "Vacuum stopping error message with Alpha's"**Reply:**Reddell, Brandon D: "Vacuum stopping error message with Alpha's"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

Dear Dr. Dutta,

There is no way to calculate d2N/dz dE with the standard options.

However, this can be achieved by adding a command USERWEIG and

modifying the user routine fluscw. I suggest that you request as many USRYIELD

detectors as the number of z intervals you want, scoring in each of

them d2N/dE d(omega), with E first scoring variable (as many intervals

as you wish) and omega second variable (1 interval covering the whole

solid angle: 0 to 180 degrees). In fluscw, you can write a logic

of this kind:

...................

INCLUDE '(SCOHLP)'

*

LOGICAL LFIRST

SAVE LFIRST, ZMIN, DELTAZ

DATA LFIRST /.TRUE./

* Initialization, to be done only at first call

IF (LFIRST) THEN

LFIRST = .FALSE.

* It is always good to write a warning on standard output that a

* user routine has been used

WRITE(LUNOUT,*)

& "**** Using routine fluscw to get differential yield in z ****"

ZMIN = .... ! write here the limits of the z range

ZMAX = .... ! where you want to do the scoring

NINTRV = .... ! and the number of z intervals

DELTAZ = (ZMAX - ZMIN)/DBLE(NINTRV)

END IF

FLUSCW = ONEONE

LSCZER = .FALSE.

* The following applies only to yield detectors (there could be others)

IF(ISCRNG .EQ. 5) THEN

* find the z interval number

INTERV = INT((ZZ - ZMIN)/DELTAZ) + 1

* don't score if detector number is different from interval number

LSCZER = JSCRNG .NE. INTERV

* multiply the score by 4pi to get integral result in omega

IF (.NOT. LSCZER) FLUSCW = FOUFOU * PIPIPI

END IF

RETURN

END

Of course, this is only a possible suggestion (and there might be errors:

I have not checked it). You may think of possible variations

Alberto

On Tue, 23 May 2006, Dipanwita Dutta wrote:

*> Hi,
*

*> I want to calculate the spatial distribution of neutron yield
*

*> for an thick cylindrical target ( radius 10 cm and length 100cm
*

*> situated from z=0 to z-100)
*

*>
*

*> i.e. d2N/dz dE where z is the beam direction. From USRYIELD
*

*> one can calculate the d2N/d(omega) dE, which doesnot
*

*> correspond to equal z interval. How to get neutrons
*

*> in equal z interval ?
*

*> If somebody can help?
*

*>
*

*> regds
*

*> Dipanwita
*

*>
*

*> --------------------------------------------------------------------
*

*> | Dr. Dipanwita Dutta | Phone: 91-22-25592041 (Off.) |
*

*> | Scientific Officer | 91-22-25589546 (Res.) |
*

*> | Nuclear Physics Division | |
*

*> | Bhabha Atomic Research Centre | Email: ddutta@apsara.barc.ernet.in |
*

*> | Trombay, Mumbai-400085, India | dutta_dipa03@yahoo.com |
*

*> --------------------------------------------------------------------
*

*>
*

*>
*

*>
*

-- Alberto Fassò SLAC-RP, MS 48, 2575 Sand Hill Road, Menlo Park CA 94025 Phone: (1 650) 926 4762 Fax: (1 650) 926 3569 fasso@slac.stanford.edu

**Previous message:**Dipanwita Dutta: "RE: LET from USRYIELD: units"**In reply to:**Dipanwita Dutta: "RE: LET from USRYIELD: units"**Next in thread:**Reddell, Brandon D: "Vacuum stopping error message with Alpha's"**Reply:**Reddell, Brandon D: "Vacuum stopping error message with Alpha's"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

*
This archive was generated by hypermail 2.1.6
: Sat May 27 2006 - 23:26:29 CEST
*