From: Alberto Fasso' (fasso@SLAC.Stanford.EDU)
Date: Fri May 26 2006 - 20:03:02 CEST
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
This archive was generated by hypermail 2.1.6 : Sat May 27 2006 - 23:26:29 CEST