Re: FLUKA: usreou.f
Dear Ludwik
FLUKA has a built-in DETECTOR card , that allows coincidences and
anti-coincideneces, but probably it is not sufficient for your problem.
Of course you can use mgdraw.f to write on file energy depositions etc, and
then
process them; this can be quite space-consuming. What I did in the past for
calorimeter study, was to implement some set of
EVENTBIN cards , and re-read them offline producing data files or n-tuples
with
the combined info for each event. Note that in EVENTBIN you can define
user-variables, with the 8 type, or "fake" the standard spatial grid by
changing
the scoring coordinates in the comscw.f or fluscw.f routines ( see card
EXTRAWEI
, but be careful, they are normally there only to change yhe scored
quantity...).
In principle, one could also access the binning data with a user routine in
FLUKA, but they are stored in the blank common with dynamic allocation, I
think
it is too complex.
Sorry, I'll not look at your input, not to be unkind, but if I start with
one....... but feel free to ask if you meet problems.
Paola
Ludwik Pienkowski wrote:
> Dear Paola,
>
> Thanks for explanations. I'll try to use mgdraw.f subroutine. However,
> I'm not sure that this is a solution for my problem.
>
> We've collected large set of experimental data, neutron multiplicity
> distributions of 1-5 GeV/c hadrons (mainly protons) spallation reactions
> on thin (few mm) and thick (up to 35 cm) targets. Part of the results are
> already published, for example: Phys.Rev C56 1909 (1997), Nucl.Instr. and
> Meth. A 414 117 (1998), and most recently Nucl.Instr. and Meth. B 170 299
> (200) including the comparison with HERMES simulations.
> The data were collected at CERN (PS208 collaboration, spoke-person
> Dr.D.Hilscher) and in Juelich (NESSI collaboration, Dr.U.Jahnke)
> using high efficiency, 4-pi, Gd-loaded, liquid scintillator detector
> (1500l tank).
>
> I hope that the neutron kinetic energy spectra I can obtain using the
> cards USRYIELD. However, our neutron detector was able to detect
> neutrons only for events with a defined trigger condition, and
> not for all beam particles.
>
> Is it possible to write a user subroutine which can simulate the
> experimental trigger conditions and accept or reject a given event?
>
> I'm a beginner using FLUKA program. My FLUKA input file is a very
> short now, only few lines. Please have a look at this input file
> which is below.
> Of course our aim is not only to make a simple simulation, but
> to simulate whole neutrons detection process in our neutron
> ball.
>
> Best regards
> Ludwik Pienkowski
> Heavy Ion Laboratory Warsaw University
>
> My FLUKA input file:
>
> TITLE
> Test of 2 GeV proton on Pb
> DEFAULTS
>NEW-DEFA
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
>12345678
> BEAM -2. PROTON
> BEAMPOS 0.0 0.0 -20.0
> GEOBEGIN 0.0
>COMBINAT
> BNB thin/thick target 10cm diameter, 2cm thick
> SPH 1 0.0 0.0 0.0 +100000.0
> SPH 2 0.0 0.0 0.0 25.0
> ZCC 3 0.0 0.0 5.0
> XYP 4 -1.0
> XYP 5 1.0
> END
> * black hole
> bh1 5 +1 -2
> * vacuum
> va2 5OR +2 -3OR +2 +3 +4OR +2 +3 -5
> * target
> pb3 5 +3 -4 +5
> END
> GEOEND
> *
> * Lead material n. 16, normal temperature, overrides default mat. 16
> MATERIAL 82.0 207.19 11.35 16.0 LEAD
> LOW-MAT 16.0 82.0 -2.0 293.0 LEAD
> LOW-NEUT
> *
> * target
> ASSIGNMAT 16.0 3.0
> * External Black Hole
> ASSIGNMAT 1.0 1.0
> * Vacuum
> ASSIGNMAT 2.0 2.0
> *
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
>12345678
> *
> USRYIELD 1401. 8. 48. 2. 1. 1. BNB
> USRYIELD 2. 0. 100. 3.1415926 0. 163. &
> USRYIELD
>BEAMDEF
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
>12345678
> START 10000. 99999999. 1.0 1.0
> STOP
>
> and PEMF input file for lead:
>
> ELEM
> &INP IRAYL=1, RHO=11.35, ISSB=1, IEV=823.0, CBAR=6.2018, X0=0.3776,
> X1=3.8073, AFACT=0.09359, SK=3.1608, DELTA0=0.14 &END
> LEAD
> PB
> ENER
> &INP AE=1.511, UE=100000001., AP=.5, UP=100000000. &END
> PWLF
> &INP NALE=300, NALG=400, NALR=100 &END
> DECK
> &INP &END
Partial thread listing: