Status of Mokka Fluka Integration

1
Status of Mokka / Fluka Integration

• Aim
• Method
• Current Status

2
Aim

• Systematic comparison between Mokka and Fluka
  implementation of test beam configuration
• Particularly interesting for hadronic
  interactions
• See DRWs Geant studies
• Want flexibility to
• Adapt to new mokka detector models
• Investigate full TDR type geometry
• Original problems
• Fluka geometry defined by data cards
• Only limited geometrical structures supported
• Repeated structures at 1 level only

3
Method

• Geomety and physics decoupled in G4 and Fluka
• Wrappers for f77/C
• Fluka authors tests
• Simple detectors
• Identical results Fluka,FlukaG4
• T36 calorimeter 81 layers Pb (10mm)-scint.(2.5mm)
  
• Consistent results
• My first test
• Use T36 calorimeter as above

From ATL-SOFT-98-039
4
Initial T36 calorimeter tests

• Normal Fluka use calculate particle fluence
  through defined regions (x,y,z) (radiation
  damage, dosage, etc.)
• Want complete flexibility in information output
• Extract hits with user
• supplied scoring on
• boundary crossing
• in mokka format
• e.g. Lateral response of
• T36 to 10 GeV ?-

5
Boundary Crossing Detection

• User routine callable when boundary between
  volumes detected during tracking a particle
• G4 replica or parametrised volumes Fluka
  lattice volumes
• Not seen during tracking
• Repeated placement of identical G4 physical
  volumes are seen
• Solution (from Flugg author)
• Use rudimentary drawing routine, gives user
  control at each step in tracking
• For Mokka study itself, can use Mokka CGA to
  determine boundaries within drawing routine

(Missing!) replicated volumes
6
Current status

• Mokka running within flugg/Fluka framework
• Using Mokka-01-00
• Flugg05 (Jan. 2003)
• Fluka 2002.1b (Dec. 2002)
• Procedure start from Mokka release and remove
• all classes except for detector construction,
  detector parametrisation, magnetic field
  construction
• corresponding include, variable and class
  definitions in .cc/.hh
• anything related to G4RunManager,
  DetectorMessenger
• code where SensitiveDetector is set
• interactive code, visualisation, etc.
• Some difficulties, e.g. with Mokka SD which
  contains sensitive detector definitions and
  rotation matrices

7
Operation
/home/nkw/lc/calice/flugg05/examples/fluggEx/emp
tyEx/src wildcard .cc -rw-r--r-- 1 nkw
opal 23998 Mar 13 1211 Control.cc
-rw-r--r-- 1 nkw opal 4275 Mar 6
1613 CGAInit.cc -rw-r--r-- 1 nkw opal
10139 Mar 6 0946 TPC00.cc -rw-r--r--
1 nkw opal 13878 Mar 5 1146 SD.cc
-rw-r--r-- 1 nkw opal 27482 Mar 5
1144 Proto01.cc -rw-r--r-- 1 nkw opal
20384 Mar 5 1140 Proto00.cc -rw-r--r--
1 nkw opal 17489 Mar 5 1046
CGAGeometryManager.cc -rw-r--r-- 1 nkw
opal 24553 Mar 5 1046 Ecal02.cc
-rw-r--r-- 1 nkw opal 20648 Mar 5
1046 Hcal02.cc -rw-r--r-- 1 nkw opal
30122 Mar 5 1045 Hcal03.cc -rw-r--r--
1 nkw opal 20115 Mar 5 1045
HcalWMod.cc -rw-r--r-- 1 nkw opal
5052 Mar 5 1036 Hodoscope00.cc -rw-r--r--
1 nkw opal 10789 Mar 5 0940
HECSD.cc -rw-r--r-- 1 nkw opal
6664 Mar 4 1519 ECSD.cc -rw-r--r-- 1 nkw
opal 5277 Mar 4 1244
VSubDetectorDriver.cc -rw-r--r-- 1 nkw
opal 8014 Aug 1 2002 SET00.cc
-rw-r--r-- 1 nkw opal 2197 Aug 1
2002 Coil00.cc -rw-r--r-- 1 nkw opal
1529 Jul 22 2002 CGAUtils.cc -rw-r--r--
1 nkw opal 4762 Jul 19 2002
Yoke00.cc -rw-r--r-- 1 nkw opal
3329 Jul 19 2002 Mask00.cc -rw-r--r-- 1 nkw
opal 7795 Feb 20 2002
MySQLWrapper.cc -rw-r--r-- 1 nkw opal
29959 Feb 20 2002 MyPlacement.cc
-rw-r--r-- 1 nkw opal 3333 Feb 20
2002 Field00.cc

• Two stage process
• First pass, to generate material definitions,
  fluka to G4 correspondance
• Second pass, to run fluka with G4 geometry

8
First pass, G4 ? Fluka conversion
B dens. 2.34g/cm3, nElem 1 Stored
as RON B2O3 dens. 2.34g/cm3, nElem 2
Stored as B2O3 Soda dens. 2.3g/cm3, nElem
2 Stored as SODA Pyrex dens.
2.23g/cm3, nElem 4 Stored as PYREX
TetraFluoroEthane dens. 0.00455g/cm3, nElem
3 Stored as TETRAFLU Isobutane dens.
0.00267g/cm3, nElem 2 Stored as
ISOBUTAN RPCGAS1 dens. 0.004032g/cm3,
nElem 4 Stored as RPCGAS1 Mix dens.
2.15747g/cm3, nElem 9 Stored as MIX
Printing FLUKA materials... Printing FLUKA
compounds... G4PhysicalVolumeStore
(0x40213a60) has 561 volumes. Printing
ASSIGNMAT... Printing Magnetic Field... No
field found... NavHistWithCountNavHistWithCount
createdNavHistWithCountNavHistWithCount
createdNavHistWithCountNavHistWithCount
createdgt Flugg FluggNavigatorUpdateNavigatorH
istory(0xcbb6f18) fHistoryHistory
depth0Level0 Phys NameWorldPhysical
TypeN
Connecting to the database models00Building
sub_detector P55, geometry db P55, driver
proto00Ecal prototype with 30 layers (10 x 2.8
mm, 10 x 2.8 mm and 10 x 4.2 mm) in vertical, 5 x
5 cells per waferConnecting to the database
P55Sub_detector P55 DONE!Fe-gtGetRadlen()
17.5854 mmAl-gtGetRadlen() 88.9253
mmAr-gtGetRadlen() 109831 mmAir-gtGetRadlen()
285161 mmSiO2-gtGetRadlen() 122.948
mmEpoxy-gtGetRadlen() 315.979
mmG10-gtGetRadlen() 162.017 mmKapton-gtGetRadlen
() 285.597 mmPolystyrene-gtGetRadlen()
425.434 mmB-gtGetRadlen() 224.928
mmB2O3-gtGetRadlen() 164.123 mmNa2O-gtGetRadlen(
) 126.811 mmPyrex-gtGetRadlen() 126.518
mmTFE-gtGetRadlen() 77265.4 mmIsobutane-gtGetRad
len() 169446 mmRPCGAS1-gtGetRadlen() 89867.3
mmCGAGeometryManager starting the detector
construction Asking for the model
P55Building Proto...total_W_plates
30MixDensite 2.15747 g/cm3Mix-gtGetRadlen()
75.0202 mmHalfDetectorX 362.4, HalfDetectorY
83.775, HalfDetectorZ 171.4Proto
done.Detector construction done.
G4PhysicalVolumeStore (0x40213a60) has 561
volumes. Iterating... Storing information...
Connecting to the database models00Detector
construction done. G4PhysicalVolumeStore
(0x40213a60) has 561 volumes. Iterating...
Storing information... Tungsten dens.
19.3g/cm3, nElem 1 Stored as TUNGSTEN
TungstenModified dens. 11g/cm3, nElem 1
Stored as TUNGST02 Copper dens. 8.96g/cm3,
nElem 1 Stored as COPPER Silicium
dens. 2.33g/cm3, nElem 1 Stored as
SILICIUM Iron dens. 7.87g/cm3, nElem
1 Stored as IRON Printing FLUKA
materials... Printing FLUKA compounds...
G4PhysicalVolumeStore (0x40213a60) has 561
volumes. Printing ASSIGNMAT... Printing
Magnetic Field... No field found...
NavHistWithCountNavHistWithCount
createdNavHistWithCountNavHistWithCount
createdNavHistWithCountNavHistWithCount
createdgt Flugg FluggNavigatorUpdateNavigatorH
istory(0xcbb6f18) fHistoryHistory
depth0Level0 Phys NameWorldPhysical
TypeN
9
Operation
MATERIAL 74.0 183.840 1.930e01
3.0 TUNGSTENMATERIAL
74.0 183.840 1.100e01 4.0
TUNGST02LOW-MAT 4.0

TUNGSTENMATERIAL 29.0 63.546 8.960e00
5.0 COPPERMATERIAL
14.0 28.090 2.330e00 6.0
SILICIUMMATERIAL 26.0 55.850
7.870e00 7.0
IRONMATERIAL 13.0 26.980 2.700e00
8.0 ALUMINUMMATERIAL
18.0 39.950 1.780e-03 9.0
ARGONMATERIAL
1.290e-03 10.0
AIRMATERIAL 7.0 14.010 9.990e-01
11.0 NITROGENMATERIAL
8.0 16.000 9.990e-01 12.0
OXIGENMATERIAL
2.200e00 13.0
QUARTZMATERIAL 14.0 28.090 9.990e-01
14.0 SILICONMATERIAL
1.300e00 15.0
EPOXYMATERIAL 1.0 1.010
9.990e-01 16.0
HYDROGENMATERIAL 6.0 12.010 9.990e-01
17.0 CARBONMATERIAL
1.700e00 18.0
G10MATERIAL 17.0 35.450
9.990e-01 19.0
CHLORINEMATERIAL 1.420e00
20.0 KAPTONMATERIAL
1.032e00 21.0
POLYSTYRMATERIAL
2.340e00 22.0
BORONMATERIAL 5.0 10.000 9.990e-01
23.0 10.0B10MATERIAL
5.0 11.000 9.990e-01 24.0
11.0B11MATERIAL 2.340e00
25.0 B2O3MATERIAL
2.300e00 26.0
SODAMATERIAL 11.0 22.990 9.990e-01
27.0 SODIUMMATERIAL
2.230e00 28.0
PYREXMATERIAL
4.550e-03 29.0
TETRAFLUMATERIAL 9.0 18.998 9.990e-01
30.0 FLUORINEMATERIAL
2.670e-03 31.0
ISOBUTANMATERIAL
4.032e-03 32.0
RPCGAS1MATERIAL 2.157e00
33.0 MIX
10
Operation
GEANT4 VOLUMES
.......1........2.....
...3........4........5........6........7...
1 WorldPhysical 102
WafferPhysical 103
CuPhys 104 CuPhys
2 SensWafferPhys
3 WSlabPhys
4 WafferPhysical 5
WafferPhysical 6
WafferPhysical 7
WafferPhysical 8
WafferPhysical 9
WafferPhysical 10
WafferPhysical 11
WafferPhysical 12
WafferPhysical 13
WafferPhysical 14
WafferPhysical 15
WafferPhysical etc...
546 LayerPhys 547
LayerPhys 548 LayerPhys
549 LayerPhys
550 LayerPhys
551 LayerPhys 552
LayerPhys 553
LayerPhys 554 LayerPhys
555 LayerPhys
556 LayerPhys
557 LayerPhys 558
LayerPhys 559
LayerPhys 560 LayerPhys
561 LayerPhys

11
Ongoing work

• Understand runtime warnings, e.g.
• Check again Fluka user provided inputs, e.g.
• low energy neutron cross-section
• PEMF, etc.,
• Extract information using CGA to Mokka format
• Compare with DRWs work

NavHistWithCountNavHistWithCount
createdNavHistWithCountNavHistWithCount
createdgt Flugg FluggNavigatorUpdateNavigatorH
istory(0xcbb7278) fHistoryHistory
depth0Level0 Phys NameWorldPhysical
TypeN) ...lt Flugg FluggNavigatorUpdateNav
igatorHistory(0xcbb7278)WARNING - Attempt to
delete the physical volume store while geometry
closed !WARNING - Attempt to delete the logical
volume store while geometry closed !WARNING -
Attempt to delete the solid store while geometry
closed !
12
Summary

• Identified easy way of comparing G4/Fluka
• Cannot use G-Fluka, completely obsolete
• Dont want to use standalone Fluka
• inefficient, for variations in geometry
• Improved understanding, moderate progres
• Integration with Mokka geometry OK with some work
• Debugging phase, to understand Fluka/Flugg