Hi Vladimir
Did you consider what is written in the Manual Note 1 of the USRYIELD
card?
1) While option USRBDX calculates angular distributions WITH RESPECT TO
THE NORMAL to the boundary at the point of crossing, USRYIELD's
distributions are calculated WITH RESPECT TO A FIXED DIRECTION (the
beam direction, or a different direction specified by the user with
SDUM = BEAMDEF).
The NORMAL and the beam direction coincide only at the back face of the
target, not for scoring of particles exiting laterally. You could test it
by defining a very thin region at the back end where the lateral
contribution is negligible. Here both scorings should give the same
result.
Cheers
Stefan
On Fri, 11 Jan 2019, babintsev wrote:
>
> Hello
>
> could you please help me understand the output of different
> scoring cards ?
>
> I am interested in the number of particles from target.
> I use the example from tutorial course.
>
> example : ... let us calculate the output of CHARGE Hadrons
>
> 1) from USRBDX card (unit 50) the output of Charge Had into AIR :
> -----------
> Sp3ChH :
> Tot. resp. (Part/cmq/pr) 4.5073396E-03 +/- 1.239388
> %
> ( --> (Part/pr) 1.486836 +/- 1.239388
> % )
>
> and I can get cumulative value :
>
> 5.5792654E-05 = 5.4125518E-05 * (5.011875 - 3.981074 )
>
> ... ok ...
>
> 2) from URSYIELD card (unit 57) the output of CHARGE Hadrons :
> -------------
>
> YieAng2 :
>
> Detector n: 2 YieAng2
> (User norm: 1.00000000
> sigma: 1.00000000 mb
> distr. scored: 218 ,
> from reg. 4 to reg. 5 )
> linear 1st variable (x1) binning from 0.0000E+00 to 3.1416E+00
> 18 bins ( 1.7453E-01 wide)
> 2nd variable (x2) ranges from 0.0000E+00 to 1.0000E+01
> 1st variable (x1) is: Laboratory Angle (radians)
> 2nd variable (x2) is: Laboratory Kinetic Energy
> The scored double differential yield is (normalized per primary
> particle):
> plain d2 N / d x1 d x2 where x1, x2 are
> the first and second variables
>
>
> 1st variable (x1) differential is a solid angle one
>
> Tot. response (integrated over x1) 0.1308400 +/- 1.208114
> %
> ( --> Tot. response without normaliz. 0.1308400 +/-
> 1.208114 % )
>
> ... the total number is 0.1308400
> for energy interval - 10 GeV :
>
> total number = 0.1308400 * 10 = 1.308400
>
> so, we have difference with initial value: 1.486836
>
> ... and I cannot get cumulative distribution here :
>
> for the 1-st value :
>
> 4.3209166E-05 * 0.1745329 = 7.54142104856139949e-06
> but in the output I see :
>
> 2.0000010E-05
>
> ---------------------
>
> so, I have the following questions:
>
> *) integral for charge hadrons are different for
>
> USRBDX and URSYIELD cards : 1.486836 and 1.308400
>
> *) and for the distribution I cannot get cumulative values in
>
> URSYIELD card output
>
> *) how to connect total number of charged hadrons from
>
> USRBDX card : 1.486836
>
> with Integral from USRBIN card : 1.68 (from plot)
>
> ---------------------------
>
> thank you,
> Vladimir
>
>
>
>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Mon Jan 14 2019 - 14:17:59 CET