1 / 12

Status of Mokka / Fluka Integration

Status of Mokka / Fluka Integration. Aim Method Current Status. 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

jmaher
Download Presentation

Status of Mokka / Fluka Integration

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Status of Mokka / Fluka Integration • Aim • Method • Current Status Nigel Watson / RAL PPD-Birmingham

  2. 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

  3. Method • Geomety and physics decoupled in G4 and Fluka • Wrappers for f77/C++ • Fluka authors’ tests • Simple detectors • Identical results Fluka,Fluka+G4 • T36 calorimeter: 81 layers Pb (10mm)-scint.(2.5mm) • Consistent results • My first test • Use T36 calorimeter as above [From ATL-SOFT-98-039] Nigel Watson / RAL PPD-Birmingham

  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 - Nigel Watson / RAL PPD-Birmingham

  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… Nigel Watson / RAL PPD-Birmingham

  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 Nigel Watson / RAL PPD-Birmingham

  7. Operation /home/nkw/lc/calice/flugg05/examples/fluggEx/emptyEx/src: wildcard *.cc -rw-r--r-- 1 nkw opal 23998 Mar 13 12:11 Control.cc -rw-r--r-- 1 nkw opal 4275 Mar 6 16:13 CGAInit.cc -rw-r--r-- 1 nkw opal 10139 Mar 6 09:46 TPC00.cc -rw-r--r-- 1 nkw opal 13878 Mar 5 11:46 SD.cc -rw-r--r-- 1 nkw opal 27482 Mar 5 11:44 Proto01.cc -rw-r--r-- 1 nkw opal 20384 Mar 5 11:40 Proto00.cc -rw-r--r-- 1 nkw opal 17489 Mar 5 10:46 CGAGeometryManager.cc -rw-r--r-- 1 nkw opal 24553 Mar 5 10:46 Ecal02.cc -rw-r--r-- 1 nkw opal 20648 Mar 5 10:46 Hcal02.cc -rw-r--r-- 1 nkw opal 30122 Mar 5 10:45 Hcal03.cc -rw-r--r-- 1 nkw opal 20115 Mar 5 10:45 HcalWMod.cc -rw-r--r-- 1 nkw opal 5052 Mar 5 10:36 Hodoscope00.cc -rw-r--r-- 1 nkw opal 10789 Mar 5 09:40 HECSD.cc -rw-r--r-- 1 nkw opal 6664 Mar 4 15:19 ECSD.cc -rw-r--r-- 1 nkw opal 5277 Mar 4 12:44 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 Nigel Watson / RAL PPD-Birmingham

  8. + 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... NavHistWithCount::NavHistWithCount createdNavHistWithCount::NavHistWithCount createdNavHistWithCount::NavHistWithCount created==> Flugg FluggNavigator::UpdateNavigatorHistory(0xcbb6f18) +fHistory=History depth=0Level=[0]: Phys Name=[WorldPhysical] Type=[N] First pass, G4  Fluka conversion Connecting to the database models00Building sub_detector P55, geometry db P55, driver proto00: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 waferConnecting to the database P55Sub_detector P55 DONE!Fe->GetRadlen() = 17.5854 mmAl->GetRadlen() = 88.9253 mmAr->GetRadlen() = 109831 mmAir->GetRadlen() = 285161 mmSiO2->GetRadlen() = 122.948 mmEpoxy->GetRadlen() = 315.979 mmG10->GetRadlen() = 162.017 mmKapton->GetRadlen() = 285.597 mmPolystyrene->GetRadlen() = 425.434 mmB->GetRadlen() = 224.928 mmB2O3->GetRadlen() = 164.123 mmNa2O->GetRadlen() = 126.811 mmPyrex->GetRadlen() = 126.518 mmTFE->GetRadlen() = 77265.4 mmIsobutane->GetRadlen() = 169446 mmRPCGAS1->GetRadlen() = 89867.3 mmCGAGeometryManager starting the detector construction: Asking for the model P55:Building Proto...total_W_plates = 30MixDensite = 2.15747 g/cm3Mix->GetRadlen()= 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... NavHistWithCount::NavHistWithCount createdNavHistWithCount::NavHistWithCount createdNavHistWithCount::NavHistWithCount created==> Flugg FluggNavigator::UpdateNavigatorHistory(0xcbb6f18) +fHistory=History depth=0Level=[0]: Phys Name=[WorldPhysical] Type=[N] Nigel Watson / RAL PPD-Birmingham

  9. Operation *MATERIAL 74.0 183.840 1.930e+01 3.0 TUNGSTENMATERIAL 74.0 183.840 1.100e+01 4.0 TUNGST02LOW-MAT 4.0 TUNGSTENMATERIAL 29.0 63.546 8.960e+00 5.0 COPPERMATERIAL 14.0 28.090 2.330e+00 6.0 SILICIUMMATERIAL 26.0 55.850 7.870e+00 7.0 IRONMATERIAL 13.0 26.980 2.700e+00 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.200e+00 13.0 QUARTZMATERIAL 14.0 28.090 9.990e-01 14.0 SILICONMATERIAL 1.300e+00 15.0 EPOXYMATERIAL 1.0 1.010 9.990e-01 16.0 HYDROGENMATERIAL 6.0 12.010 9.990e-01 17.0 CARBONMATERIAL 1.700e+00 18.0 G10MATERIAL 17.0 35.450 9.990e-01 19.0 CHLORINEMATERIAL 1.420e+00 20.0 KAPTONMATERIAL 1.032e+00 21.0 POLYSTYRMATERIAL 2.340e+00 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.340e+00 25.0 B2O3MATERIAL 2.300e+00 26.0 SODAMATERIAL 11.0 22.990 9.990e-01 27.0 SODIUMMATERIAL 2.230e+00 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.157e+00 33.0 MIX** Nigel Watson / RAL PPD-Birmingham

  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 Nigel Watson / RAL PPD-Birmingham

  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 DRW’s work NavHistWithCount::NavHistWithCount createdNavHistWithCount::NavHistWithCount created==> Flugg FluggNavigator::UpdateNavigatorHistory(0xcbb7278) +fHistory=History depth=0Level=[0]: Phys Name=[WorldPhysical] Type=[N]) ...<== Flugg FluggNavigator::UpdateNavigatorHistory(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 ! Nigel Watson / RAL PPD-Birmingham

  12. Summary • Identified “easy” way of comparing G4/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

More Related