RE: LET from USRYIELD: units

From: Alberto Fasso' (fasso@SLAC.Stanford.EDU)
Date: Fri May 26 2006 - 20:03:02 CEST

  • Next message: Sunil C.: "Read USRBDX binary output"

    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
    

  • Next message: Sunil C.: "Read USRBDX binary output"

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