Geant4 Hadronic Physics Working group progress and status.

Transcription

1 Geant4 Hadronic Physics Working group progress and status. J.P. Wellisch

2 Outline 6WDWXVRQPLOHVWRQHVDQGUHFHQW GHYHORSPHQWV 9DOLGDWLRQYHULILFDWLRQ 1HZVRQRXWVLGHFRQWDFWV UG SDUWLHV

3 The dry numbers 1XPEHURISDFNDJHV 5HOHDVHG 7RWDO 1XPEHURIFODVVHV 5HOHDVHG 7RWDO /LQHVRIVRXUFHFRGH 5HOHDVHGa 7RWDOa

4 $WSUHVHQWDERXWSHRSOHDUH FRQWULEXWLQJWRWKLVHIIRUWZLWKVRPHRI WKHLUWLPHFUHDWLYLW\RUH[SHUWLVH 7KHQXPEHURIXVHFDVHSDFNDJHV FRQVLGHUHGLVFXUUHQWO\DQGZH SURYLGHDWRWDORISK\VLFVOLVWVIRUWKH YDULRXVDUHDVRIDSSOLFDELOLW\

5 The milestones. 'LVWULEXWHµHGXFDWHGJXHVV SK\VLFVOLVWVIRUPDMRUXVHFDVHV FDVHV,QFOXGHDWOHDVWRQHWHVWEHDPVLPXODWLRQLQSUH EHDPVLPXODWLRQLQSUHUHOHDVH:*OHYHOUHOHDVH:*OHYHO YDOLGDWLRQ,PSURYHGYHULILFDWLRQVXLWHIRUWKHFDVFDGHHQHUJ\UDQJH 5HOHDVHELDVHG0$56UHZULWHIRUHQHUJLHVEHORZ ZULWHIRUHQHUJLHVEHORZ*H9,QFOXGHγQXFOHDUUHDFWLRQVLQTXDUN QXFOHDUUHDFWLRQVLQTXDUNJOXRQVWULQJPRGHO,PSURYHWKHFKDUJHVWDWHWUHDWPHQWIRUUHFRLOVUHVLGXDOV %ULQJNLQHWLFPRGHOWRDUHOHDVDEOHVWDWHRQJRLQJ 5HOHDVHRIDFDVFDGHFRGHIURP+(7&PLOHVWRQHRQJRLQJ 3URYLGHDJHQHULFVFDWWHULQJWHUPIRUFDVFDGHW\SHPRGHOV,PSURYHHOHFWURQXFOHDUFURVV QXFOHDUFURVVVHFWLRQWRLQFOXGHKDUGVFDWWHULQJVHFWLRQWRLQFOXGHKDUGVFDWWHULQJ 3XEOLVKZRUN DWOHDVWSDSHUVVXEPLWWHGWRUHIHUHHGMRXUQDOV GRQH UHDFKHG

6 Particle physics relevant requirements collected Specifically from LHC: Have starting physics list (done) Improve information flow on V&V (done, to be verified) Provide a cascade code (two focused efforts) Fix known problems in low energy models (released)

7 A complete sample requirement. Name and E Dennis Wright, SLAC Title: Pion and kaon nuclear prodution cross-sections sections Description: 1-51 GeV pions and kaons interacting inelastically.. 10% precision would be great. Rationale: Representation of how the shower develops. Supporting use-cases that require this: Trying to model hadronic interactions in the BaBar intsrumented flux return. Responsible cathegory: hadronics Fulfillment criterion: Comparison to the data from particles interacting in the beam-pipe and flux return. Relevace: : very highly relevant for BaBar,, relevant also for LHCb References: An E from Dennis pointing to the data.

8 Requirements collected titles only During the last geant4 workshop, the users workshop at SLAC, and in private mails: Ensure that the physics reference manual match the implementation, n, and the models are mentioned in the applications developers guide, Referencing papers is just fine. Memory usage for G4NDL cross-sections sections Energy and momentum conservation should be checked in regression independently by the working group for all models, and publish the t test-suite. suite. Use well known international benchmarks to validate; publish results. r More understandable hadronic physics lists Get documentation on which model is good/usefull usefull/required for which use-case

9 Requirements Cont. Provide a set of plots, and put them u; i.e. provide a place to put these plots for users. Cross-sections sections for n,p inelastic scattering below 150 MeV in CMS tracker materials at 10% level of precisions. Pion and kaon nuclear prodution cross-sections: sections: 1-51 GeV pions and kaons interacting inelastically.. 10% precision would be great % level of description for GeV 100GeV incident protons for example on Beryllium or copper. Enroll a set of users to validate on complete application; as beta testers so to say. Compare inclusive and exclusive cross-sectios sectios to data from the RAL/Durham database

10 Requirements Cont. Description on how to set cuts, and its effects Possibility to stop low energy neutral particles (like neutrons) Each model should be specified concerning its application area/use-cases Provide a list of models per use-case package List of models per use-case package Include physics list samples, once they exists, into the phsics editor Ensure tracability of data to the primary source Well defined process for updating the databases on request. Parametrizations of hadronic showers in CsI and Iron. Parametrizations of neutron background in LHC experiments Models for alpha incident inelastic reactions

11 Requirements Cont. Models for alpha incident difractive dissociation reactions Where the modelling approach allows to produce the residual, it should be provided. Neutron production by alphas at energies below 10 MeV; ; including reaction cross-sections sections at 20%precision and kinematics of neutrons and gammas produced. Include K0 oszillations Provide muon nuclear reactions Provide internal conversion Provide neutron elastic scattering, in particular recoil energy and momentum distributions for neutrons below 10MeV. Dito for n inelastic scattering off Xenon and SiO2, CaCO3, H2O

12 Requirements Cont. Dito for capture Provide gamma nuclear reactions for gamma energies of less than 100 MeV,, including cross- sections. Provide radioactive decay after transmutation. Provide k-shell k excitation in radioactive decay Activation of detector material and environment by shower particles

13 Note: Almost all of these requirements are by now fulfilled. Many were fulfilled when the issue was entered as a requirement, so only information was to be provided.the rest are to be scheduled for being addressed, according to priorities.

14 Requirements cont. IUHVKUHTXLUHPHQWVKDUYHVWHGGXULQJ WKHJHDQWZRUNVKRSODVWZHHN 'LUHFWLQWHUDFWLRQZLWKWKHH[SHUWVLV YHU\SURGXFWLYHLQWKLV

15 Sample plots energy deposition BTEV: All distributions are in the expected energy range

16 Active tasks 2002 Write educated guess physics lists for major use-cases Include at least one test-beam simulation into regular validation Include a complete radiation benchmark into WG level validation Improved validation suite for the cascade energy range Possibly further extension of the high energy validation suite Plan to contribute to SATIF-6 Release fully leading particle biased mars-5 5 re-write Release of cascade part of HETC re-write Improve gamma nuclear reactions in QGS model Make a validation/verification WWW page

17 Active tasks 2002, cont. Possibly first release of high energy heavy ion reactions in QGS model, with option to further extension to QMD Revision of the reaction cross-sections sections Improve the charge state treatment for recoils. Bring kinetic model to releasable state Bring inucl cascade code to releasable state Research the use of CHIPS in string fragmentation for intrinsically 3D fragmentation Provide a generic scattering term for cascade type models Alternative coherent elastic model (reggee( theory based)

18 Active tasks 2002, cont. Improve electro-nuclear cross-sections sections to take hard scattering into account. Investigate JENDL2.2, and LA150 neutron data libraries Collect (even more) requirements Release work, coordination Contribute to maintenance and user support Contribute to architecture working group Contribute to process improvement/establishment Contribute to training Publish work (11 papers in the plan )

19 Conclusions :HKDYHDW7,0(12:PHWPLOHVWRQHVIRU :HJHWDFRQVWDQWIORZRIUHTXLUHPHQWVIHHEDFNDQG WHVWEHDPUHVXOWVIURPWKHGHWHFWRUJURXSVWKLVLV H[FHOOHQWQHZV :HKDYHDJRRGWHDPDQGZHDUHQRUPDOO\DEOHWR DWWUDFWWKHH[SHUWLVHQHHGHGIRUPRGHOLQJ :HIXOO\GHSHQGRQYLVLWRUDQGWUDYHOPRQH\IURP&(51 )URQWOLQHVXSSRUWPDQ OLQHVXSSRUWPDQSRZHUQRQWULYLDOWRILQG

20 Some hadronic physics highlights of late 2001 and HXWURQVSHFWUDIURPSUHHTXLOLEULXPGHFD\DQG HTXLOLEULXPGHFD\DQG DVQHDNSUHYLHZRQNLQHWLFPRGHOSHUIRUPDQFH TJV PRGHOIRU SLRQ DQGNDRQ NDRQDQGJDPPD LQGXFHGUHDFWLRQV 'RSSOHUEURDGHQLQJRQWKHIO\,QWHUQDOFRQYHUVLRQDQGDQHZSKRWRQ HYDSRUDWLRQGDWDEDVH &KLUDO LQYDULDQWSKDVHVSDFHGHFD\ VSDFHGHFD\ $SURSDJDWLRQWHVWIRUTXDQWXPPROHFXODU G\QDPLFV

21 Swapping to show a few transparencies on pre-compound neutron yields.

22 Preview on kinetic model 160 MeV p on Pb, forward neutrons 585MeV p on Al, forward And backward n and π

23 Low energy neutrons: G4NDL0.2, 3.7 Are granular selections of data from (alphabetic) Brond 2.1 CENDL 2.2 EFF-3 ENDF/B (VI.0, VI.1, VI.5) ENSDF FENDL/E2.0 JEF 2.2 JENDL (3.1, 3.2, FF, 3.3 currently under study) MENDL-2(P) Large parts of the selection is guided by the FENDL-2 2 selection G4NDL0.2 for non-thermal application

24 The neutron_hp transport models Simulate the cross-sections sections and interactions of neutrons with kinetic energies below 20 MeV down to thermal energies. The upper limit is set only by the evaluated data libraries the code is based on. We consider elastic scattering, fission, capture and inelastic scattering as separate models Neutron_hp sampling codes for the ENDF/B-VI derived data formats are completely generic (not including (not including general R-matrix R for the time being) Note that for fission there is a quite competitive theory driven alternative model, G4ParaFissionModel.

25 Models for neutron interaction and thermalization. QHXWURQBKSPRGHOVDQGFURVVVHFWLRQV VHFWLRQV 8VHVWKHXQL[ ILOHV\VWHPWRHQVXUH JUDQXODUDQGWUDQVSDUHQWDFFHVVXVDJHRI GDWDVHWV 0RUHWKDQAHYHQWVUXQ 8VHVSRLQWZLVHFURVV ZLVHFURVVVHFWLRQVÎ QR DUWLIDFWVGXHWRPXOWLJURXSVWUXFWXUH

26 Doppler broadening 'RHVH[DFW GRSSOHU EURDGHQLQJRQWKHIO\ EDVHGRQ.GDWDÎ QRSUHIRUPDWWLQJRI GDWDWRIL[HGWHPSHUDWXUHVDQGHDV\ VLPXODWLRQRIVHWXSVZLWKPL[HG XSVZLWKPL[HG WHPSHUDWXUHV $GGVWKH GRSSOHU ELDVWRWKHQXFOHDU PRPHQWXPGLVWULEXWLRQ 3RLQWRQHLVWRWKHEHVWRIRXUNQRZOHGJHQRW DYDLODEOHIURPDQ\RWKHUWUDQVSRUWFRGHWKH VHFRQGLVDOVRLQ0&13

27 The doppler bias illustrated for Carbon

28 qgs model for π and K induced reactions 3RPHURQWUDMHFWRU\DQGYHUWH[ SDUDPHWHUVWXQHGWRGHVFULEHHODVWLF WRWDODQGGLIIUDFWLYHDVVXPHG FURVVVHFWLRQVIRU VHFWLRQVIRUNDRQDQGSLRQ VFDWWHULQJRIIQXFOHRQV 1RWXQLQJRQILQDOVWDWHGLVWULEXWLRQV $IHZSORWVWRLOOXVWUDWHWKHTXDOLW\RI SUHGLFWLRQ

29 K-,, scattering off Au (for pions see V&V section)

30 Photon Evaporation data base 2ULJLQDOO\FRQWDLQLQJDGRSWHGOHYHODQGJDPPDUD\ UD\ WUDQVLWLRQHQHUJLHVSKRWRQLQWHQVLW\PXOWLSRODULW\ PXOWLSRODULW\ KDOIOLIHDQGVSLQSDULW\IRULVRWRSHVXSWR= OLIHDQGVSLQSDULW\IRULVRWRSHVXSWR= $ ([SDQGHGWRLQFOXGHSUREDELOLW\RILQWHUQDOFRQYHUVLRQ DQGLQWHUQDOFRQYHUVLRQFRHIILFLHQWV,&&IURPVKHOOV.///00000DQG1 %DVHGRQ(16')GDWDIURP/%1/DQGWDEXODWHG WKHRUHWLFDO,&&GDWDIURP%DQGHWDO HWDO XVHGIRU= DQG DQG5 VHOHWDOHWDO XVHG =

31 ,&&V DUHFDOFXODWHGE\ FXELFVSOLQH LQWHUSRODWLRQ XVLQJDERYHWDEOHVDWWKH UHTXLUHGJDPPDUD\ UD\ HQHUJ\,&&FDOFXODWHGIRU0L[HG PXOWLSRODULW\0(LI PL[LQJUDWLRDYDLODEOH 6RPHFKDQJHVZHUH LQWURGXFHGLQWKHIRUPDW RIWKHGDWDEDVHHQWULHV WRNHHSWKHVL]HRIWKH ILOHVGRZQGDWDEDVHLV QRZWLPHVODUJHU

32 Preliminary test results (16')GHFD\GDWDSURFHVVHGZLWK5$'/,67 %1/FRGHDQG*HDQWIRUGHFD\V &V

33 57 Co

34 A sample development: Chiral Invariant Phase-space space Decay. $TXDUNOHYHOGLPHQVLRQDOHYHQWJHQHUDWRUIRU IUDJPHQWDWLRQRIH[FLWHGKDGURQLF KDGURQLFV\VWHPVLQWR KDGURQV %DVHGRQWKH4&'LGHDRIDV\PSWRWLFIUHHGRP /RFDOFKLUDO FKLUDOLQYDULDQFHUHVWRUDWLRQOHWVXVFRQVLGHULQYDULDQFHUHVWRUDWLRQOHWVXVFRQVLGHU TXDUNSDUWRQVPDVVOHVV SDUWRQVPDVVOHVVDQGZHFDQLQWHJUDWHWKH LQYDULDQWSKDVHVSDFHGLVWULEXWLRQRITXDUN VSDFHGLVWULEXWLRQRITXDUN SDUWRQVDQGTXDUNH[FKDQJHIXVLRQPHFKDQLVP DQGTXDUNH[FKDQJHIXVLRQPHFKDQLVP RIKDGURQL]DWLRQ 7KHRQO\QRQNLQHPDWLFDOFRQFHSWXVHGLVWKDWRID WHPSHUDWXUHRIWKHKDGURQLF KDGURQLFV\VWHPTXDVPRQ

35 Vacuum CHIPS 7KLVDOORZVWRFDOFXODWHWKHGHFD\RIIUHHH[FLWHG KDGURQLF V\VWHPV,QDQILQLWHWKHUPDOL]HG WKHUPDOL]HG V\VWHPRI1SDUWRQV ZLWK WRWDOPDVV0WKHLQYDULDQWSKDVHVSDFHLQWHJUDOLV VSDFHLQWHJUDOLV 2N 4 SURSRUWLRQDOWR M DQGWKHVWDWLVWLFDO M / T GHQVLW\RIVWDWHVLVSURSRUWLRQDOWR+HQFH e ZHFDQZULWHWKHSUREDELOLW\WRILQG1SDUWRQV SDUWRQVZLWK WHPSHUDWXUH7LQDVWDWHZLWKPDVV0DV dw 1RWHWKDWIRUWKLVGLVWULEXWLRQWKHPHDQPDVV VTXDUHLV 2 2 M = 2N(2N 2) T M e 2N 4 M / T dm

36 Vacuum CHIPS :HXVHWKLVIRUPXODWRFDOFXODWHWKHQXPEHU RI SDUWRQV LQDQH[FLWHGWKHUPDOL]HG WKHUPDOL]HG KDGURQLF V\VWHPDQGREWDLQWKHSDUWRQ SDUWRQVSHFWUXP dw kdk 1 2k M 7RREWDLQWKHSUREDELOLW\IRUTXDUNIXVLRQLQWR KDGURQVZHFDQQRZFRPSXWHWKHSUREDELOLW\ WRILQGWZRSDUWRQV SDUWRQVZLWKPRPHQWD TDQGN ZLWKWKHLQYDULDQWPDVVµ N 4 2q 2 2kq(1 cosθ ) P( k, M, µ ) = 1 δ µ qdqd cosθ M 1 2k M 1 2k M N 3

37 Vacuum CHIPS 8VLQJWKHGHOWDIXQFWLRQWRSHUIRUPWKHLQWHJUDWLRQ DQGWKHPDVVFRQVWUDLQWZHILQGWKHWRWDO NLQHPDWLFDOSUREDELOLW\RIKDGURQL]DWLRQ KDGURQL]DWLRQRIDSDUWRQ ZLWKPRPHQWXPNLQWRDKDGURQ KDGURQ ZLWKPDVVµ: $FFRXQWLQJIRUVSLQDQGTXDUNFRQWHQWRIWKHILQDO VWDWHKDGURQ KDGURQDGGVVDQGDFRPELQDWRULDODGGVVDQGDFRPELQDWRULDO IDFWRU $WWKLVOHYHORIWKHODQJXDJH&+,36FDQEHDSSOLHG WRSSEDU SEDUDQQLKLODWLRQDQQLKLODWLRQ M 4k( N k 3) 2 ( 2 ) N 1 µ 2kM 3

38 Anti proton annihilation

39 Anti proton annihilation

40 Nuclear CHIPS,QRUGHUWRDSSO\&+,36IRUDQH[FLWHGKDGURQLF V\VWHPZLWKLQQXFOHLZHKDYHWRDGGSDUWRQ H[FKDQJHZLWKQXFOHDUFOXVWHUVWRWKHPRGHO 7KHNLQHPDWLFDOSLFWXUHLVWKDWDFRORUQHXWUDO TXDVPRQHPLWVD HPLWVDSDUWRQZKLFKLVDEVRUEHGE\D QXFOHRQRUDQXFOHDUFOXVWHU7KLVUHVXOWVLQD FRORUHGUHVLGXDOTXDVPRQ TXDVPRQDQGDFRORUHG FRPSRXQG 7KHFRORUHGFRPSRXQGWKHQGHFD\VLQWRDQ RXWJRLQJQXFOHDUIUDJPHQWDQGDµUHFRLO TXDUN WKDWLVLQFRUSRUDWHGE\WKHFRORUHGTXDVPRQ TXDVPRQ

41 Nuclear CHIPS $SSO\LQJPHFKDQLVPVDQDORJXHWRYDFXXP&+,36 ZHFDQZULWHWKHSUREDELOLW\RIHPLVVLRQRID QXFOHDUIUDJPHQWZLWKPDVVµ DVDUHVXOWRIWKH WUDQVLWLRQRIDSDUWRQ SDUWRQZLWKPRPHQWXPNIURPWKH TXDVPRQWRDIUDJPHQWZLWKPDVV µ DV P 2( k ) µ ( k ) k, µ, µ ) = 1 d cosθ [ + ] µ + k(1 cosθ kq ) 2 µ k(1 cosθ kq ) +HUHQLVWKHQXPEHURITXDUNSDUWRQV SDUWRQVLQWKH QXFOHDUFOXVWHUDQG LVWKHFRYDULDQWELQGLQJ HQHUJ\RIWKHFOXVWHUDQGWKHLQWHJUDOLVRYHUWKH DQJOHEHWZHHQSDUWRQ SDUWRQDQGUHFRLOSDUWRQ n 3 ( 2 kq

42 Nuclear CHIPS 7RFDOFXODWHWKHIUDJPHQW\LHOGVLWLV QHFHVVDU\WRFDOFXODWHWKHSUREDELOLW\WRILQGD FOXVWHURIν QXFOHRQVZLWKLQDQXFOHXV:HGR WKLVXVLQJWKHIROORZLQJDVVXPSWLRQV $IUDFWLRQε1 RIDOOQXFOHRQVLVQRWFOXVWHULVLQJ FOXVWHULVLQJ $IUDFWLRQε2 RIWKHQXFOHRQVLQWKHSHULSKHU\RI WKHQXFOHXVLVFOXVWHULQJLQWRWZRQXFOHRQFOXVWHUV 7KHUHLVDVLQJOH FOXVWHUL]DWLRQ SUREDELOLW\ω DQGILQGZLWKDEHLQJWKHQXPEHURI QXFOHRQVLQYROYHGLQFOXVWHUL]DWLRQ FOXVWHUL]DWLRQ a 1 C ω ν ν Pν = a 1 (1 + ω)

43 Nuclear CHIPS $WWKLVOHYHORIWKHODQJXDJH&+,36FDQ EHDSSOLHGWRFDSWXUHRISLRQV SLRQVDQG SKRWRQXFOHDUUHDFWLRQV

44 Intra-nuclear nuclear CHIPS ([WHQVLRQVWRLQFOXGHWKHEHKDYLRURI PXOWLSOH TXDVPRQV ZLWKLQRQHQXFOHXV KDYHEHHQDGGHG

45 Hard scattering in electro-nuclear

46 Hard scattering in electro-nuclear

47 A propagation test for QMD development 6RPHFKDUDFWHULVWLFVRI4'0 $NLQHPDWLFDOFDVFDGHZLWKGHWDLOHGPRGHOLQJRI WKHQXFOHXV 1XFOHDU+DPLOWRQLDQFDOFXODWHGIURPDQGERG\ SRWHQWLDOVRIDOOKDGURQVSUHVHQWLQWKHV\VWHP 6ROYLQJWKHHTXDWLRQRIPRWLRQE\LQWHJUDWLQJWKLV WLPHGHSHQGHQW+DPLOWRQLDQ 6FDWWHULQJWHUPLQWHUPVRIORFDOL]HGLQWHUDFWLRQV DQGGHFD\V (WF

48 The support process static view

49 The support process dynamic view

50 Test-beams +DGURQLF WHVWEHDPFRPSDULVRQVFRPHIURPFROODERUDWLRQRI H[SHULPHQWV GHWHFWRUJURXSVZLWKµFRUH JHDQWSHUVRQQHO $7/$67LOHWHVWEHDP &067LOHWHVWEHDP $7/$6+(&WHVWEHDP $7/$6)&$/WHVWEHDP %7(9FU\VWDOWHVWEHDP &06FRPELQHGWHVWEHDP EHDP &V, WHVWEHDPEHQFKPDUN */$67VWDUWLQJWHVWEHDP 3ORWVEHLQJVROLFLWHGDVFRXUWHV\RIWKHH[SHULPHQWDOJURXSV

51 Other areas of known usage (likely incomplete) 7UDFNHUSHUIRUPDQFH $7/$6&06%D%DU 0HGLFDO 8SSVDOD7(5$ 1HXWURQGRVLPHWU\ GRVLPHWU\PHDVXUHPHQWEHDPOLQHV 612/RV$ODPRV&(5136'R' 'R'&DQHWF 5DGLDWLRQVFKLHOGLQJ VFKLHOGLQJDFWLYDWLRQWKHUPDOL]DWLRQWKHUPDOL]DWLRQ '<1$0,;0(&2$/,&("&06(6$HWF 2LOVHDUFKDQGVLPLODU 0LWVXELVKL*HQHUDOHOHFWULFV(;;21$/&$7(/«

52 The hopefully no longer dry numbers 1XPEHURISDFNDJHV 5HOHDVHG 7RWDO 1XPEHURIFODVVHV 5HOHDVHG 7RWDO /LQHVRIVRXUFHFRGH 5HOHDVHGa 7RWDOa

53 Conclusions :HKDYHDJRRGWHDP :HDUHQRUPDOO\DEOHWRDWWUDFWWKH H[SHUWLVHZHQHHG :HIXOO\GHSHQGRQYLVLWRUDQGWUDYHO PRQH\IURP&(51 )URQWOLQHVXSSRUWPDQ OLQHVXSSRUWPDQSRZHUQRQ WULYLDOWRILQG

54 Conclusions 7KHPDLQIRFXVRIDOOWKHVHGHYHORSPHQWVLVRI FRXUVHRQ/+&DQG%D%DU VKRZHUSK\VLFVDQG GRVLPHWU\ $//HIIRUWVWKDWZDQWWRFRQWULEXWHWRSK\VLFVLQ WKHJHDQWFRQWH[WDUHZHOFRPH 3K\VLFVPRGHOLQJSK\VLFV9 9DQGSK\VLFV UHVHDUFKLVERWKWKHVFRSHDQGFRQFHUQRIWKH JHDQWKDGURQLF KDGURQLFZRUNLQJJURXS 1RWHWKDWZHVWULYHWRPDNHVXUHWKDWLQGLYLGXDO DFWLYLWLHVDUHLQWHJUDWHGWRDYRLGGXSOLFDWLRQRIZRUN EXWDOVRWULYLDOPLVWDNHV

55 Collaboration with 3 rd parties Some of the reasoning: Geant3 had used two strategies. There were shower packages released with geant3, and there were interfaces released with geant3; the latter were interfacing to external packages. The first was a working model, for the latter, geant3 always was claimed to be obsolete. GISMO: the no physics situation, but only interfacing to external packages. They never really got support for the use of these codes with GISMO. MCNPX: Gets it right. They encourage and help 3 rd parties to release MCNP interfaces with their 3 rd party code. It solves the support question.

56 Collaboration with 3 rd parties Basis: Basis: We provide a set of well defined, published, and highly stable interfaces that allows interested 3 rd parties to release adapters to use their code, or to use geant4 physics implementations within their infrastructure. EGS: geant4 chips code for γ-nuclear reactions also in EGS HETC: Being re-written to become natively available in G4 INUCL: Being integrated to become natively available in G4 UrQMD: In the process of being re-engineered engineered to become natively available in geant4 MCNP: Discussion on using the geant4 interfaces in MCNP FLUKA: Interfaced by air shower users for their own use. Liege Cascade code: Discussion in progress. We hope that they will release a G4 interface soon, and are of course happy to help. EGS: HETC: INUCL: UrQMD MCNP: G-FLUKA: Liege

57 Conclusions,WLVYHU\LPSRUWDQWWKDWLQGLYLGXDOFRQWULEXWRUV DUHHQIUDQFKLVHGWRMRLQWKHFROODERUDWLRQLQ SDUWLFXODULQWKHDUHDRISK\VLFVPRGHOLQJ 7KH\PXVWIHHODVVXUHGWKDWWKH\DUHZHOO SURWHFWHGIURPDQ\DWWHPSWWRGHSULYHWKHPRII RUFRS\VWHDOVXEOWLOLVHU VXEOWLOLVHUWKHZRUNWKDWEXLOWWKHLU FDUHHUV,HWKHLUFRGH 7KH\RWKHUZLVHZRXOGEHDVNHGWRFRQWULEXWHDW WKHLURZQSHULO :HVKRXOGH[SOLFLWO\VWDWHDSROLF\HQVXULQJWKLVLQ WKH0R8 0R8UHYLVLRQ