Transcript TFluka::InitPhysics()

TFluka Status
Alberto Fasso
Endre Futo
Andreas Morsch
ALICE Offline Week 1/9/2003
Outline
●
Updates
●
Physics Configuration
●
Optical Properties for Cerenkov Photons
●
User produced secondaries
Updates
●
Code updated to work with ...
–
Latest Fluka version
●
●
–
Cannot use DPMJET as primary particle generator
and TFluka for tracking
Geant4 5.2
●
–
Clash of names for subroutines between DPMJET
and FLUKA
Loading sequence of libraries
New IO
●
FlukaConfig.C
Physics configuration
●
Strategy
–
Collect information from all SetProcess and
SetCuts calls from Config.C
–
Generate consistent set of FLUKA input cards in
TFluka::InitPhysics()
Processes and Cuts
●
Not always one to one correspondence with
G3
–
LOSS = 2 (complete ionisation loss fluctuations)
–
Particle cuts are not by region (except electron,
positron and gamma via EMFCUT)
Processes and Cuts
●
To do
–
FLUKA input cards need number of materials
and number of regions which are only available
after geometry initialization
–
Wait for TGeo
Cerenkov Optical Properties
●
Two possibilities
–
Use FLUKA input cards which define optical
properties (refraction index, absorption length,
quantum efficiency) as 3rd order polynomials
●
–
Precision too bad to describe ALICE RICH response
Use user routines
●
queffc.cxx
●
rflctv.cxx
●
rfrndx.cxx
●
abscff.cxx
Cerenkov Optical Properties
●
Where to store optical properties
–
In TFluka object as list of materials in which
Cerenkov photons should be produced and/or
transported
●
Prototype by SuSt Martina Bielikova
–
–
–
TGeoCerenkovMaterial
TClonesArray + Look-Up Table
In TGeo ?
Secondaries
●
●
FLUKA produced secondaries are managed
by FLUKA
–
TVirtualMCStack object only for recording
(passive)
–
Flagged as secondaries that have been already
transported
User produced secondaries are put on
TVirtualMCStack object
–
Example: Feedback Photons
–
Flagged as secondaries to be transported
Secondaries
●
●
When source is called by feeder for the
transport of the next particle there are three
possibilities
–
primary particle
–
tracked secondary
–
untracked secondary
Code works correctly if
–
PostTrack() and FinishPrimary() as well as
PreTrack() are only called if next particle is
primary
Tests
●
●
●
Code has been tested for TPC and RICH hit
production including feedback photon
production
However tests limited by FLUGG because of
crashes especially for transport of electrons,
gammas and Cerenkov photons
Wait for TGeo
–
On the critical path
–
Should be highest priority now !