Dear Adrian,
Do I understand correctly, that you want to score the energy deposition of photo- and Compton-electrons (and their secondaries) produced by the primary photon before it was subject to Rayleigh scattering, but that you want to exclude the energy deposition for photo- and Compton-electrons (and their secondaries) produced by a photon after it was subject Rayleigh scattering?
If yes, then this cannot be done with the simple if statement below. For the moment you only score the energy deposition of primary photons as well as all particles which have a generation number 3 (why?). Photons subject to Rayleigh scattering change the generation number, so if their second process is photo-electric effect or Compton scattering then the energy deposited by this third generation electrons will be counted.
If you only want to score the energy deposition by second generation electrons, it can be easily done with some similar IF statement (using JTRACK to identify that it is an electron, and LLO to identify that it is second generation). However, if you also want to score the energy deposited by the descendants of second generation electrons, then you need to use dedicated variables to keep track of particle histories (see page 49++ in
https://indico.cern.ch/event/489973/contributions/2000440/attachments/1272042/1972478/09_AdvancedUserRoutines2016.pdf). With this variables you can make particles "remember" if they are descendants of second generation electrons. Alternatively, you could just kill primary photons after they were subject to Rayleigh scattering. In both cases, it requires however a bit of programming and you would need to get more familiar with the user routines (please look at the course slides).
Cheers, Anton
________________________________________
From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of RADIOLOGY [a.knyziak_at_gum.gov.pl]
Sent: 08 December 2016 15:22
To: fluka-discuss_at_fluka.org
Subject: Re: [fluka-discuss]: the first interaction with a primary photon
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
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id¬c_info
Received on Fri Dec 09 2016 - 12:21:07 CET