RE: particle distribution after the collimator

From: <Andrea.Mairani_at_mi.infn.it>
Date: Fri, 13 Nov 2009 18:30:08 +0100 (CET)

Dear Vitaly,
 I've compiled/linked in this way:
    (FROM THE MANUAL)
    * compile the modified routines (with the fff script on LINUX/UNIX):
$FLUPRO/fff mgdraw.f (produces a new file mgdraw.o)

    * link them (with the lfluka script on LINUX/UNIX) to the FLUKA library:

         $FLUPRO/lfluka -o myfluka -m fluka mgdraw.o

This will produce a new executable (indicated here as myfluka). To run the
new executable, launch the usual rfluka script with the option -e myfluka.

I've run it and I've looked at the file "distribv1001_results.txt", and it
is not empty:

1.0000 29.906 0.2026608E-02-0.9287661E-04 10.00000
0.88145287E-040.15333768E-030.99999998
1.0000 29.906 88.14643 153.3361 1000000.0
0.88145287E-040.15333768E-030.99999998
1.0000 1.8091 -0.8620296 0.2507047E-01 10.00000 -.30555668
 0.88358964E-020.95213289
1.0000 1.8091 -305554.5 8835.835 952133.6 -.30555668
 0.88358964E-020.95213289
1.0000 29.903 0.2368082E-02 0.1537777E-02 10.00000
0.11509366E-030.16178755E-030.99999998
1.0000 29.903 115.0949 161.7875 1000000.0
0.11509366E-030.16178755E-030.99999998
1.0000 29.907 -0.4543866E-02-0.1401436E-02 10.00000
-.83927268E-03-.50902004E-040.99999965

Ciao!
Andrea

> Dear FLUKA experts,
>
> I succeeded with the FLUKA compilation. But, unfortunately, another
> problem is occurred. There are no data in the file RESULTS.TXT which I
> created in order to write the needed data there. Could you say what is
> wrong? In the attachment you can find the modified version of MGDRAW.F and
> unchanged input file.
>
> Thank you for any suggestion.
>
> Best regards,
> Vitaliy
>
>
>
>> -----Original Message-----
>> From: owner-fluka-discuss_at_mi.infn.it
>> [mailto:owner-fluka-discuss_at_mi.infn.it] On Behalf Of
>> Andrea.Mairani_at_mi.infn.it
>> Sent: 11 2009 . 13:53
>> To: fluka-discuss_at_fluka.org
>> Subject: Re: particle distribution after the collimator
>> Importance: High
>>
>>
>> Dear Vitaly,
>> I looked at your mgdraw.f, some comments:
>> * you have to define
>> CHARACTER*8 MRGNAM, NRGNAM
>> * instead of
>> IF ((MRGNAM.EQ."TARG".OR.MRGNAM.EQ."VOID").and.
>> & JTRACK.EQ.1) THEN
>> you have to use
>> IF ((MRGNAM.EQ."TARG".AND.NRGNAM.EQ."VOID").and.
>> & JTRACK.EQ.1) THEN
>> if you want to score all particles leaving the target.
>> On the contrary if you are interested on the forward emission
>> after the target you have to define a tiny region at the end
>> of the RCC ( for example using XYP 9.995).
>>
>> Ciao
>> Andrea
>>
>>
>> > Dear FLUKA experts,
>> >
>> > Recently I asked the FLUKA community about the way to obtain the
>> > particle distribution after the collimators.
>> >
>> > ...I would kindly ask you to help me in solution of one problem. We
>> > are going to calculate the passage of the particles through the
>> > graphite collimators in details. This means that we want to
>> know the
>> > change of original particle distribution including particles going
>> > directly through the graphite. There is a certain particle
>> > distribution before the collimator. Then the particles pass mostly
>> > through the collimator slit but part of them pass directly
>> through the
>> > collimator material. Graphite is not so effective as iron
>> in sense of
>> > particle cutting. The radiation length for graphite is about 30 cm
>> > only. Therefore, we need longer collimators but, unfortunately, we
>> > have not enough space for them. We are interested in
>> particles which
>> > can come to the transfer beam acceptance together with antiprotons,
>> > for example, pions. Pions have similar rigidity in our case and
>> > theoretically can be transport to the ring and make some
>> harm there.
>> > Could you recommend some proper way to obtain such a particle
>> > distribution after collimator. Do we need some additional
>> code? The
>> > idea is to have a particle distribution before and after
>> collimator in
>> > order to understand how many particles can come to the
>> septum magnet
>> > of the following storage ring. It gives us an information about
>> > possible problems...
>> >
>> > There were several useful recommendations but,
>> unfortunately, I didn't
>> > solve the problem up to now. I tried to use the MGDRAW
>> subroutine and
>> > wrote some lines in mgdraw.f (namely, in entry BXDRAW). I'm
>> not sure
>> > that this is correct and the best variant. In addition, maybe
>> > something wrong with my input file. Nevertheless, there are some
>> > errors occurred during compilation. Maybe I'm using the
>> CALL GEOR2N in
>> > wrong way because there are some errors due to different type of
>> > variables. I didn't understand how it really works.
>> >
>> > I appreciate any help. Thank you in advance.
>> >
>> > Best regards,
>> > Vitaliy Gostishchev
>> >
>> > p.s. In the attachment you can find my input and the
>> modified MGDRAW.F
>> > files. <<distribv1.inp>> <<mgdraw.f>>
>> >
>>
>>
>
>
Received on Fri Nov 13 2009 - 20:34:35 CET

This archive was generated by hypermail 2.2.0 : Fri Nov 13 2009 - 20:34:35 CET