Status of Mokka Fluka Integration Aim Method Current

Status of Mokka / Fluka Integration <Aim <Method <Current Status Nigel Watson / RAL PPD-Birmingham 1 Imperial, 18/03/2003

Aim < Systematic comparison between Mokka and Fluka implementation of test beam configuration < Particularly interesting for hadronic interactions } See DRW’s 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 Nigel Watson / RAL PPD-Birmingham 2 Imperial, 18/03/2003

Method < Geomety and physics decoupled in G 4 and Fluka < Wrappers for f 77/C++ < Fluka authors’ tests } Simple detectors Þ Identical results Fluka, Fluka+G 4 } T 36 calorimeter: 81 layers Pb (10 mm)-scint. (2. 5 mm) Þ Consistent results < My first test } Use T 36 calorimeter as above [From ATL-SOFT-98 -039] Nigel Watson / RAL PPD-Birmingham 3 Imperial, 18/03/2003

Initial T 36 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 T 36 to 10 Ge. V Nigel Watson / RAL PPD-Birmingham 4 Imperial, 18/03/2003

Boundary Crossing Detection < User routine callable when boundary between volumes detected during tracking a particle < G 4 replica or parametrised volumes Fluka “lattice volumes” } Not seen during tracking < Repeated placement of identical G 4 physical volumes are seen < Solution (from Flugg author) } Use rudimentary drawing routine, gives user control at each step in (Missing!) tracking replicated } For Mokka study itself, can use Mokka CGA to determine volumes… boundaries within drawing routine Nigel Watson / RAL PPD-Birmingham 5 Imperial, 18/03/2003

Current status < Mokka running within flugg/Fluka framework } Using Mokka-01 -00 } Flugg 05 (Jan. 2003) } Fluka 2002. 1 b (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 G 4 Run. Manager, Detector. Messenger } code where Sensitive. Detector is set } interactive code, visualisation, etc. < Some difficulties, e. g. with Mokka SD which contains sensitive detector definitions and rotation matrices Nigel Watson / RAL PPD-Birmingham 6 Imperial, 18/03/2003

Operation lice/flugg 05/examples/flugg. Ex/empty. Ex/src: nkw nkw nkw nkw nkw nkw opal opal opal opal opal opal 23998 4275 10139 13878 27482 20384 17489 24553 20648 30122 20115 5052 10789 6664 5277 8014 2197 1529 4762 3329 7795 29959 3333 Mar Mar Mar Mar Aug Jul Jul Feb Feb 13 6 6 5 5 5 5 5 4 4 1 1 22 19 19 20 20 20 12: 11 16: 13 09: 46 11: 44 11: 40 10: 46 10: 45 10: 36 09: 40 15: 19 12: 44 2002 2002 < Two stage process Control. cc CGAInit. cc TPC 00. cc SD. cc Proto 01. cc Proto 00. cc CGAGeometry. Manager. cc Ecal 02. cc Hcal 03. cc Hcal. WMod. cc Hodoscope 00. cc HECSD. cc VSub. Detector. Driver. cc SET 00. cc Coil 00. cc CGAUtils. cc Yoke 00. cc Mask 00. cc My. SQLWrapper. cc My. Placement. cc Field 00. cc Nigel Watson / RAL PPD-Birmingham < First pass, to generate material definitions, fluka to G 4 correspondance < Second pass, to run fluka with G 4 geometry 7 Imperial, 18/03/2003

First pass, G 4 Fluka conversion Connecting to the database models 00 Building sub_detector P 55, geometry db P 55, driver proto 00: Ecal 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 wafer Connecting to the database P 55 Connecting to the database models 00 Sub_detector P 55 DONE! Detector construction done. * G 4 Physical. Volume. Store (0 x 40213 a 60) has 5 * Storing information. . . + Tungsten: dens. = 19. 3 g/cm 3, Stored as TUNGSTEN + Tungsten. Modified: dens. = 11 Stored as TUNGST 02 + Copper: dens. = 8. 96 g/cm 3, n Stored as COPPER + Silicium: dens. = 2. 33 g/cm 3, Stored as SILICIUM + Iron: dens. = 7. 87 g/cm 3, n. El Stored as IRON * Printing FLUKA materials. . . * Printing FLUKA compounds. . . * G 4 Physical. Volume. Store (0 x 40213 a 60) has 5 * Printing ASSIGNMAT. . . * Printing Magnetic Field. . . No field found. . . Nav. Hist. With. Count: : Nav. Hist. With. Count created Building Proto. . . Nav. Hist. With. Count: : Nav. Hist. With. Count created total_W_plates = 30 Nav. Hist. With. Count: : Nav. Hist. With. Count created Mix. Densite = 2. 15747 g/cm 3 ==> Flugg. Navigator: : Update. Navigator. History(0 xcbb Mix->Get. Radlen()= 75. 0202 mm +f. History=History depth=0 Half. Detector. X= 362. 4, Half. Detector. Y= 83. 775, Half. Detector. Z= 171. 4 Level=[0]: Phys Name=[World. Physical] Type=[N] Proto done. Fe->Get. Radlen() = 17. 5854 mm Al->Get. Radlen() = 88. 9253 mm Ar->Get. Radlen() = 109831 mm Air->Get. Radlen() = 285161 mm Si. O 2 ->Get. Radlen() = 122. 948 mm Epoxy->Get. Radlen() = 315. 979 mm G 10 ->Get. Radlen() = 162. 017 mm Kapton->Get. Radlen() = 285. 597 mm Polystyrene->Get. Radlen() = 425. 434 mm B->Get. Radlen() = 224. 928 mm B 2 O 3 ->Get. Radlen() = 164. 123 mm Na 2 O->Get. Radlen() = 126. 811 mm Pyrex->Get. Radlen() = 126. 518 mm TFE->Get. Radlen() = 77265. 4 mm Isobutane->Get. Radlen() = 169446 mm RPCGAS 1 ->Get. Radlen() = 89867. 3 mm CGAGeometry. Manager starting the detector construction: Asking for the model P 55: Detector construction done. * G 4 Physical. Volume. Store (0 x 40213 a 60) has 561 volumes. Iterating. . . * Storing information. . . Nigel Watson / RAL PPD-Birmingham 8 Imperial, 18/03/2003

Operation * MATERIAL LOW-MAT MATERIAL MATERIAL MATERIAL MATERIAL MATERIAL MATERIAL MATERIAL MATERIAL * * 74. 0 29. 0 14. 0 26. 0 13. 0 183. 840 1. 930 e+01 183. 840 1. 100 e+01 63. 546 28. 090 55. 850 26. 980 39. 950 7. 0 8. 0 14. 010 16. 000 14. 0 28. 090 1. 0 6. 0 1. 010 12. 010 17. 0 35. 450 5. 0 10. 000 11. 0 22. 990 9. 0 18. 998 8. 960 e+00 2. 330 e+00 7. 870 e+00 2. 700 e+00 1. 780 e-03 1. 290 e-03 9. 990 e-01 2. 200 e+00 9. 990 e-01 1. 300 e+00 9. 990 e-01 1. 700 e+00 9. 990 e-01 1. 420 e+00 1. 032 e+00 2. 340 e+00 9. 990 e-01 2. 340 e+00 2. 300 e+00 9. 990 e-01 2. 230 e+00 4. 550 e-03 9. 990 e-01 2. 670 e-03 4. 032 e-03 2. 157 e+00 Nigel Watson / RAL PPD-Birmingham 3. 0 4. 0 TUNGSTEN TUNGST 02 TUNGSTEN COPPER SILICIUM IRON ALUMINUM ARGON AIR NITROGEN OXIGEN QUARTZ SILICON EPOXY HYDROGEN CARBON G 10 CHLORINE KAPTON POLYSTYR BORON 10. 0 B 10 11. 0 B 11 B 2 O 3 SODA SODIUM PYREX TETRAFLU FLUORINE ISOBUTAN RPCGAS 1 MIX 5. 0 6. 0 7. 0 8. 0 9. 0 10. 0 11. 0 12. 0 13. 0 14. 0 15. 0 16. 0 17. 0 18. 0 19. 0 20. 0 21. 0 22. 0 23. 0 24. 0 25. 0 26. 0 27. 0 28. 0 29. 0 30. 0 31. 0 32. 0 33. 0 9 Imperial, 18/03/2003

Operation *********** GEANT 4 VOLUMES *********** * *. . . +. . 1. . +. . 2. . +. . 3. . +. . 4. . +. . 5. . +. . 6. . +. . 7. . . * 546 Layer. Phys 1 World. Physical 547 Layer. Phys 102 Waffer. Physical 548 Layer. Phys 103 Cu. Phys 549 Layer. Phys 104 Cu. Phys 550 Layer. Phys 2 Sens. Waffer. Phys 551 Layer. Phys 3 WSlab. Phys 552 Layer. Phys 4 Waffer. Physical 553 Layer. Phys 5 Waffer. Physical 554 Layer. Phys 6 Waffer. Physical 555 Layer. Phys 7 Waffer. Physical 556 Layer. Phys 8 Waffer. Physical 557 Layer. Phys 9 Waffer. Physical 558 Layer. Phys 10 Waffer. Physical 559 Layer. Phys 11 Waffer. Physical 560 Layer. Phys 12 Waffer. Physical 561 Layer. Phys 13 Waffer. Physical 14 Waffer. Physical 15 Waffer. Physical & etc. . . Nigel Watson / RAL PPD-Birmingham 10 Imperial, 18/03/2003

Ongoing work < Understand runtime warnings, e. g. Nav. Hist. With. Count: : Nav. Hist. With. Count created ==> Flugg. Navigator: : Update. Navigator. History(0 xcbb 7278) +f. History=History depth=0 Level=[0]: Phys Name=[World. Physical] Type=[N] ). . . <== Flugg. Navigator: : Update. Navigator. History(0 xcbb 7278) 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 ! < Check again Fluka user provided inputs, e. g. } low energy neutron cross-section } PEMF, etc. , < Extract information using CGA to Mokka format < Compare with DRW’s work Nigel Watson / RAL PPD-Birmingham 11 Imperial, 18/03/2003

Summary < Identified “easy” way of comparing G 4/Fluka } Cannot use G-Fluka, completely obsolete } Don’t 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 Nigel Watson / RAL PPD-Birmingham 12 Imperial, 18/03/2003