Re: Ionisation losses fluctuations

From: Antonello Carloni <antonello.carloni_at_iss.infn.it>
Date: Thu, 26 Aug 2010 01:21:55 +0200

Dear Andrea,
thanks for the suggestions and also, emilia Basile suggested me to run at
least 5 times the same job:
$FLUPRO/flutil/rfluka -N0 -M5 ppp.
to obtain a good fluctuation result and so i did.
Then I joined all the files together with usbsuw.f and otained the final
result with command usbrea.f.
Sometimes it seems so easy but you do not notice the "lines" between the
lines in the manual.
I posted this reply as a future reference for beginners like me! Thanks a
lot for the help, and do not worry, you will hear again from me, i'm sure
.....
Thanks again to Emilia Basile too!
Antonello Carloni

----- Original Message -----
From: <Andrea.Mairani_at_mi.infn.it>
To: <fluka-discuss_at_fluka.org>
Sent: Monday, August 23, 2010 12:21 PM
Subject: Re: Ionisation losses fluctuations

> Dear Antonello,
> I run your input (attached below) changing the output unit as I mentioned
> in the previous answer and the USRBIN results look reasonable for me.
> Andrea
>
>
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> TITLE
> Beam divergence and energy on pad
> DEFAULTS 0.0 0.0 0.0 0.0 0.0 0.0
> HADROTHE
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * 65 MeV proton beam, gaussian profile 5mm in x and y axes
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> BEAM -0.065 0.0
> 0.0 -0.5 -0.5 -1.PROTON
> * proton source is set at a distance of 5 cm from ion chamber window
> BEAMPOS -1.487 +1.487 -5.0 0.0 0.0 0.
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> GEOBEGIN 0.0 0.0 0.0 0.0 0.0
> 0.COMBNAME
> 0 0 proton beam on ion chamber
> * outer limits of external blackhole
> RPP hol -1000.0 +1000.0 -1000.0 +1000.0 -1000.0 +1000.0
> * outer limits of vacuum
> RPP vac -500.0 +500.0 -500.0 +500.0 -500.0 +500.0
> * boundaries of ion chamber
> RPP ich -1.53 +1.53 -1.53 +1.53 -0.106 +0.106
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * vacuum detector
> RPP bec -10.0 +10.0 -10.0 +10.0 +25.0 +25.1
> * brick in the center of pad
> RPP bri -1.4925 -1.4825 +1.475 +1.4995 +0.0955 +0.096
> *strip on back of the ouw
> RPP str -1.53 +1.53 +1.482 +1.492 +0.106 +0.1065
> * ionization chamber inlet window plane
> XYP BP1 -0.104
> * ionization chamber outlet window plane
> XYP BP2 +0.096
> * plane dividing mylar from aluminium in inlet win
> XYP BP3 -0.1045
> * plane dividing upstream and downstream detector
> XYP PLA +25.05
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * phantom
> * outer water phantom
> ZCC OCO -1.487 +1.487 +25.0
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * upper border cylinder plane
> XYP BOL +100.0
> * BOTTOM border cylinder plane
> XYP BOR +50.0
> * Sketch of a pad
> * First vertical plane along the x axis
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> PLA P1 -1.000 0.000 -0.000 -1.4999 +0.000 +0.000
> * Second plane parallel to the x axis
> PLA P2 -1.000 +0.000 -0.000 -1.53 +0.000 +0.000
> * plane on top
> XYP U1 +0.096
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * plane on bottom
> PLA U2 +0.000 +0.000 -1.000 +0.000 +0.000 +0.0955
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * upper plane at 45 degrees on the xy plane
> PLA UP1 +0.707 +0.707 0.0 -1.53 +1.5225 +0.096
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * bottom plane at 45 degrees on the xy plane
> PLA UP2 +0.707 -0.707 -0.0 -1.53 +1.4525 +0.096
> * sketch of the second pad
> * First vertical plane along the y axis
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> PLA P3 0.000 1.000 0.000 0.000 +1.475 0.000
> * second vertical plane along the y axis
> PLA P4 -0.000 -1.000 0.000 0.000 +1.445 0.000
> * plane on top
> XYP U3 +0.096
> * plane on bottom
> PLA U4 +0.000 +0.000 -1.000 +0.000 +0.000 +0.0955
> * left plane at 45 degrees on the pad
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> PLA UP3 -0.707 +0.707 0.000 -1.4925 +1.475 +0.096
> * left plane at 45 degrees on pad
> PLA UP4 +0.707 +0.707 0.000 -1.4825 +1.475 +0.096
> * sketch of the third pad
> * plane on top
> XYP U5 +0.096
> * plane on bottom
> PLA U6 +0.000 +0.000 -1.000 +0.000 +0.000 +0.0955
> * First vertical plane along the y axis
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> PLA P5 0.000 -1.000 -0.000 0.000 +1.4995 0.000
> * second vertical plane along the y axis
> PLA P6 0.000 1.000 0.000 0.000 +1.53 0.000
> * right plane at 45 degrees on the pad
> PLA UP5 0.707 -0.707 -0.000 -1.4825 +1.4995 +0.096
> * left plane at 45 degrees on the pad
> PLA UP6 -0.707 -0.707 0.000 -1.4925 +1.4995 +0.096
> * sketch of the fourth pad
> * plane on top
> XYP U7 +0.096
> * plane on bottom
> PLA U8 +0.000 +0.000 -1.000 +0.000 +0.000 +0.0955
> * First vertical plane along the x axis
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> YZP P7 -1.445
> * Second vertical plane along the x axis
> PLA P8 -1.000 0.000 -0.000 -1.475 0.000 0.000
> * bottom plane at 45 degrees on the pad
> PLA UP7 -0.707 -0.707 0.000 -1.445 +1.4525 +0.096
> * upper plane at 45 degrees on pad
> PLA UP8 -0.707 0.707 0.000 -1.445 +1.5225 +0.096
> END
> * regions
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * outer black hole
> BLK 7 +hol -vac
> * vacuum around
> VAK 7 +vac -ich -bec -str -( +OCO +BOL -BOR)
> * inlet mylar window
> MYL 7 +ich +BP3
> * inlet aluminium window
> ALI 7 +ich +BP1 -BP3
> * out window
> OUW 7 +ich -BP2
> * ion chamber
> DEC 7 +ich +BP2 -BP1 -bri -(| -P1 +P2 +U1 +U2 +UP1 +UP2)
> - (| +P3 +P4 +U1 +U2 +UP3 +UP4) - (| +P5 +P6 +U5 +U6 +UP5 +UP6)
> - (| +P7 +P8 +U7 +U8 +UP7 +UP8)
> * upstream detector
> UPS 7 +bec +PLA
> * downstream detector
> DWN 7 +bec -PLA
> * pad region number one
> PA1 7 | -P1 +P2 +U1 +U2 +UP1 +UP2
> * pad region number two
> PA2 7 | +P3 +P4 +U1 +U2 +UP3 +UP4
> * pad region number three
> PA3 7 | +P5 +P6 +U5 +U6 +UP5 +UP6
> * pad region number four
> PA4 7 | +P7 +P8 +U7 +U8 +UP7 +UP8
> *brick region
> BRA 7 +bri
> * strip region
> STA 7 +str
> * outer cylinder
> CEX 7 +OCO +BOL -BOR
> END
> GEOEND
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * materials
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> MATERIAL 7.0 14.007 0.001251 7.0
> NITROGEN
> MATERIAL 1.0 1.00794 0.0000837 3.0
> HYDROGEN
> MATERIAL 6.0 12.0107 2.0 6.0
> CARBON
> MATERIAL 8.0 15.9994 0.00133 8.0
> OXYGEN
> MATERIAL 18.0 39.948 0.00166 20.0
> ARGON
> MATERIAL 13.0 26.9815 2.699 10.0
> ALUMINUM
> MATERIAL 29.0 63.546 8.960 12.0
> COPPER
> MATERIAL 0.00195 27.0
> DIOXIDE
> MATERIAL 0.00170 28.0 GAS
> MATERIAL 1.4 29.0
> MYLAR
> MATERIAL 1.42 30.0
> KAPTON
> MATERIAL 1. 1.0079 .0000899 3.0
> 1.HYDROGEN
> MATERIAL 8. 15.9994 .001429 8.0
> OXYGEN
> MATERIAL 1.0 31.0
> WATER
> COMPOUND 2. HYDROGEN 1. OXYGEN
> WATER
> COMPOUND 1.0 CARBON 2.0 OXYGEN
> DIOXIDE
> COMPOUND -30.0 DIOXIDE -70.0 ARGON GAS
> COMPOUND -4.2 HYDROGEN -62.5 CARBON -33.3 OXYGEN
> MYLAR
> COMPOUND -2.6 HYDROGEN -69.1 CARBON -7.32 NITROGEN
> KAPTON
> COMPOUND -20.98 OXYGEN
> KAPTON
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * Material assignment to regions
> ASSIGNMAT BLCKHOLE BLK
> ASSIGNMAT VACUUM VAK
> ASSIGNMAT VACUUM UPS UPS 1.0
> ASSIGNMAT VACUUM DWN DWN 1.0
> ASSIGNMAT ALUMINUM ALI ALI 1.0
> ASSIGNMAT MYLAR MYL MYL 1.0
> ASSIGNMAT GAS DEC DEC 1.0
> ASSIGNMAT KAPTON OUW OUW 1.0
> ASSIGNMAT COPPER PA1 1.0
> ASSIGNMAT COPPER PA2 1.0
> ASSIGNMAT COPPER PA3 1.0
> ASSIGNMAT COPPER PA4 1.0
> ASSIGNMAT COPPER BRA 1.0
> ASSIGNMAT COPPER STA 1.0
> ASSIGNMAT WATER CEX 1.0
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> *Fluence and divergence of proton beam through vacuum detector
> USRBDX 101.0 208.0 -32.0 UPS DWN 400
> USRBDX 0.065 0.0643 100.0 0.24 0.0 50.0&
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> * Energy released by proton beam into OCO water cylinder (GeV/cm3)
> USRBIN 11.0 208.0 -19.0 +25.0
> +100.0cylinder
> USRBIN +0.0 +50.0 1.0 1.0 1.0&
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
> RANDOMIZE 1.0
> START 10000000.
> STOP
>
>
>
>
>
>
>
>
>
>> Dear users, it is me again,
>> in my program I'm asking fluka a simple Energy (particle code 208)
>> released output on a cilindrical body with the card USRBIN with an
>> unformatted output unit.
>> When it comes to results, i process them using usbsuw and then usbrea
>> and
>> i receive my data matrix but when it comes to the fluctuations, in the
>> second part of the unit, i get always the same result 9.999 e-10 when
>> the
>> correspondin energy bin is different from zero otherwise i get zero as
>> well.
>> the card is :
>> * Energy released by proton beam into OCO water cylinder (GeV/cm3)
>> USRBIN 11.0 208.0 -19.0 +25.0
>> +100.0cylinder
>> USRBIN +0.0 +50.0 100.0 100.0 200.0&
>>
>> What do you thin the problem might be!?
>> I mean, fluctuations should be different from on value to the other and
>> not all identical....
>> Thank you all!
>> Antonello Carloni
>>
>>
>
>
Received on Thu Aug 26 2010 - 15:39:07 CEST

This archive was generated by hypermail 2.2.0 : Thu Aug 26 2010 - 15:39:14 CEST