Re: [fluka-discuss]: the first interaction with a primary photon

From: Mikhail Polkovnikov <pmk_at_ihep.ru>
Date: Fri, 9 Dec 2016 12:01:30 +0300

On 08.12.2016 17:22, RADIOLOGY wrote:
> Dear FLUKA Experts and Users,
>
> I try to get the energy deposited in a volume by electrons (and their
> progeny) produced in the first interaction with a primary photon. For
> this purpose I wrote a simple condition in comscw.f:
>
> IF (ISCRNG .EQ. 1) THEN
> IF ((JI .NE. 210) .AND. (LLO .NE. 3)) LSCZER = .TRUE.
> END IF
>
> But I am not sure it is correct.
>
> Best regards,
>
> Adrian
>
> my input file:
>
> TITLE
> the free-air ionization chamber standard
> * Set the defaults for precision simulations
> DEFAULTS PRECISIO
> * Define the beam characteristics
> BEAM -5E-05 0.904 0.904 PHOTON
> * Define the beam position
> BEAMPOS 0.0 0.0 -100.0
> GEOBEGIN COMBNAME
> 0 0
> * Black body
> SPH blkbody 0.0 0.0 0.0 100000.0
> * Void sphere
> SPH void 0.0 0.0 0.0 10000.0
> RPP frontAL -20.05 20.05 -13.0 13.0 0.0 34.30
> RPP frontAIR -19.05 19.05 -12.0 12.0 0.0 34.30
> RPP colREG -14.05 14.05 -12.0 12.0 34.30 44.30
> RPP leftAIR -19.05 -14.05 -12.0 12.0 34.30 44.30
> RPP rightAIR 14.05 19.05 -12.0 12.0 34.30 44.30
> RPP leftAL -20.05 -19.05 -12.0 12.0 34.30 44.30
> RPP rightAL 19.05 20.05 -12.0 12.0 34.30 44.30
> RPP upAL -20.05 20.05 12.0 13.0 34.30 44.30
> RPP downAL -20.05 20.05 -13.0 -12.0 34.30 44.30
> RPP backAL -20.05 20.05 -13.0 13.0 44.30 78.60
> RPP backAIR -19.05 19.05 -12.0 12.0 44.30 78.60
> END
> * Black hole
> BLKBODY 5 +blkbody -void
> * Void around
> VOID 5 +void -frontAL -colREG -leftAIR -rightAIR -leftAL
> -rightAL -upAL -downAL -backAL
> FRONTAL 5 +frontAL -frontAIR
> FRONTAIR 5 +frontAIR
> COLREG 5 +colREG
> LEFTAIR 5 +leftAIR
> RIGHTAIR 5 +rightAIR
> LEFTAL 5 +leftAL
> RIGHTAL 5 +rightAL
> UPAL 5 +upAL
> DOWNAL 5 +downAL
> BACKAL 5 +backAL -backAIR
> BACKAIR 5 +backAIR
> END
> GEOEND
> * ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
> ASSIGNMA BLCKHOLE BLKBODY
> ASSIGNMA VACUUM VOID
> ASSIGNMA ALUMINUM FRONTAL
> ASSIGNMA AIR FRONTAIR
> ASSIGNMA AIR COLREG
> ASSIGNMA AIR LEFTAIR
> ASSIGNMA AIR RIGHTAIR
> ASSIGNMA ALUMINUM LEFTAL
> ASSIGNMA ALUMINUM RIGHTAL
> ASSIGNMA ALUMINUM UPAL
> ASSIGNMA ALUMINUM DOWNAL
> ASSIGNMA ALUMINUM BACKAL
> ASSIGNMA AIR BACKAIR
> EMF
> EMFCUT -1E-06 1E-6 FRONTAL BACKAIR 10.
> EMFRAY 2. FRONTAL BACKAIR 1.
> EMFFLUO -1. ALUMINUM AIR 1.
> USRBIN 10. ENERGY -21. 14.05 12.0
> 44.30ScorRegA
> USRBIN -14.05 -12.0 34.30 100.0 100.0 100.0 &
> USRBIN 10. ENERGY -22. 20.05 13.0
> 44.30ScorRegB
> USRBIN -20.05 -13.0 34.30 100.0 100.0 100.0 &
> AUXSCORE USRBIN ELECTRON ScorRegA ScorRegB
> USERWEIG 1.
> * Set the random number seed
> RANDOMIZ 1.0
> * Set the number of primary histories to be simulated in the run
> START 1E6 0.0
> STOP
>
> On 12/07/2016 08:05 AM, Mikhail Polkovnikov <Михаил Полковников> wrote:
>> On 06.12.2016 16:06, a.knyziak_at_gum.gov.pl a.knyziak_at_gum.gov.pl wrote:
>>> Dear Mikhail,
>>>
>>> Thank you for answer. Is there easier way? Is there any examples? I
>>> can not
>>> programming in fortran. I'm learning now.
>>>
>>> Best regards,
>>>
>>> Adrian
>>>
>>> __________________________________________________________________________
>>>
>>> You can manage unsubscription from this mailing list at
>>> https://www.fluka.org/fluka.php?id=acc_info
>>
>> Dear Adrian,
>>
>> I suppose the easier way (without FORTRAN programming) is to use
>> AUXSCORE card with different type of scoring particles. At least you
>> will get a bit more information about energy deposition of different
>> kinds of particles.
>>
>> BEAMPART - energy deposition of primary photons only.
>>
>> ELECTRONS - energy deposition of electrons.
>>
>> PHOTONS - energy deposition of primary and scattered photons.
>>
>>
>> __________________________________________________________________________
>>
>> You can manage unsubscription from this mailing list at
>> https://www.fluka.org/fluka.php?id=acc_info
>>
>
> __________________________________________________________________________
>
> You can manage unsubscription from this mailing list at
> https://www.fluka.org/fluka.php?id=acc_info

Dear Adrian,

1. There is a typo in your example "JI", must be "IJ".

2. You can use JTRACK to check FLUKA particle number, for electron JTRACK=3.

3. Primary particle has LLO = 1, so second generation electrons will
have LLO = 2.

C Energy deposition from second generation electrons.
      IF (ISCRNG .EQ. 1) THEN
          IF ((JTRACK .NE. 3) .AND. (LLO .NE. 2)) LSCZER = .TRUE.
       END IF


Best regards,

Mikhail



__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Fri Dec 09 2016 - 11:53:00 CET

This archive was generated by hypermail 2.3.0 : Fri Dec 09 2016 - 11:53:05 CET