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 Mon Aug 23 2010 - 14:48:17 CEST