JP Wellisch - LCG Project - Applications Area - Cern
JP Wellisch - LCG Project - Applications Area - Cern
JP Wellisch - LCG Project - Applications Area - Cern
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
Geant4 Hadronic Physics Status and<br />
Validation for large HEP Detectors<br />
J.P. <strong>Wellisch</strong><br />
CERN/EP/SFT
Outline:<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ %ULHI RYHUYLHZ PHWKRGRORJ\ DQG<br />
XQGHUO\LQJ SK\VLFV<br />
„ &XUUHQW SLFWXUH RI YHULILFDWLRQ IRU +(3<br />
„ UREXVWQHVV &38 SK\VLFV<br />
„ &RQIURQWDWLRQ ZLWK UHTXLUHPHQWV LQ +(3<br />
„ WHVW EHDP VLPXODWLRQ H[DPSOHV
Principal software considerations:<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ Use an open and extendible architecture.<br />
„ Use implementation frameworks to achieve maximum<br />
extendibility and flexibility in the physics modeling<br />
„ We invite collaboration with all experts in the field, to render<br />
the physics base of the geant4 toolkit as rich as possible<br />
„ Use abstract interfaces to de-couple de couple and enable independent<br />
development by experts<br />
„ Ensure that the physics implementations are both accessible<br />
to and transparent for the user, so results can more easily be<br />
verified, and if deviations are found, they can be traced back,<br />
and the modeling can be adapted to the usage.
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Modeling overview<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ 7KUHH FDWHJRULHV RI PRGHOLQJ<br />
DSSURDFKHV<br />
3DUDPHWUL]DWLRQ GULYHQ PRGHOLQJ<br />
„ 3DUDPHWUL]DWLRQ<br />
„ 'DWD GULYHQ PRGHOLQJ<br />
„ 7KHRU\ GULYHQ PRGHOLQJ
Parameterisation driven models<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ 7ZR GRPDLQV<br />
„ KLJK HQHUJ\ LQHODVWLF $DFKHQ &06<br />
„ ORZ HQHUJ\ LQHODVWLF HODVWLF ILVVLRQ FDSWXUH 75,80)<br />
8%& &(51<br />
„ 6WRSSLQJ SDUWLFOHV<br />
„ EDVH OLQH 75,80) &+$26<br />
„ PX 75,80) ),'81$<br />
„ SL ,1)1 &(51 75,80)<br />
„ . &U\VWDO %DUUHO 75,80)<br />
„ DQWL SURWRQV -/$% &(51<br />
„ (OHFWURPDJQHWLF WUDQVLWLRQV RI WKH H[RWLF DWRP SULRU WR<br />
FDSWXUH HIIHFWV RI DWRPLF ELQGLQJ 1RYRVLELUVN (6$
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
Data driven models:<br />
Low energy neutron (J.P. <strong>Wellisch</strong>)<br />
„ %DVHG RQ HYDOXDWHG GDWD LQWHUQDWLRQDO VWDQGDUGV<br />
„ (1') -HI -(1'/ &(1'/ (16') %URQG ,5')<br />
)(1'/ 0(1'/<br />
„ 6DPSOLQJ FRGHV IRU (1') % 9, GHULYHG GDWD IRUPDWV<br />
„ 8VH WKH ILOH V\VWHP WR HQVXUH JUDQXODU DQG WUDQVSDUHQW<br />
DFFHVV XVDJH RI GDWD VHWV<br />
„ 'RSSOHU EURDGHQLQJ QRW VWDWLF RQ LQSXW GDWD EXW RQ WKH IO\<br />
IURP . GDWD<br />
„ %\ QRZ D YHU\ ZLGHO\ XVHG FRGH<br />
„ 5DGLRDFWLYH GHFD\ '(5$ (6$ SKRWRQ HYDSRUDWLRQ<br />
LQWHUQDO FRQYHUVLRQ HODVWLF VFDWWHULQJ 6$,' LQWHUQDO<br />
FRQYHUVLRQ HWF
Theory driven models<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
8OWUD KLJK HQHUJ\ PRGHOV<br />
„ 8OWUD<br />
„ 3DUWRQ WUDQVSRUW PRGHO LQ GLVFXVVLRQ<br />
„ +LJK HQHUJ\ PRGHOV<br />
„ µ)ULWMRI )ULWMRI W\SH VWULQJ PRGHO &(51<br />
„ 4XDUN JOXRQ 6WULQJ &(51<br />
„ 3\WKLD LQWHUIDFH /XQG &(51<br />
„ ,QWUD QXFOHDU WUDQVSRUW PRGHOV RU UHSODFHPHQWV<br />
„ +DGURQLF FDVFDGH SUH HTXLOLEULXP +,3 &(51<br />
„ %LQDU\ FDVFDGHV &(51 ,7 1RYRVLELUVN 6/$&<br />
„ 40' W\SH PRGHOV &(51 ,QVW 7K 3K\V )UDQNIXUW<br />
„ &KLUDO LQYDULDQW SKDVH VSDFH GHFD\ -/$% &(51 ,7(3<br />
„ 3DUWLDO 0DUV UHZULWH .\RWR LQ FROODERUDWLRQ ZLWK 8YLF $QG )1$/<br />
„ 'H H[FLWDWLRQ<br />
„ (YDSRUDWLRQ ILVVLRQ PXOWL IUD<strong>JP</strong>HQWDWLRQ IHUPL EUHDN XS &06
:KDW LV QHZ LQ JHDQW IRU KDGURQLF<br />
SK\VLFV<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ 7KHRUHWLFDO PRGHOV<br />
„ )LUVW UHOHDVH UHOHDVH<br />
RI %LQDU\ FDVFDGH %&<br />
„ )LUVW UHOHDVH UHOHDVH<br />
RI FODVVLFDO %HUWLQL FDVFDGH &&<br />
„ )LUVW UHOHDVH RI D JHQHULF VFDWWHULQJ WHUP IRU<br />
FDVFDGH PRGHOV<br />
„ )LUVW UHOHDVH RI LQWHUQDO FRQYHUVLRQ<br />
„ &KLUDO LQYDULDQW SKDVH VSDFH GHFD\ &+,36 IRU<br />
HOHFWUR QXFOHDU ZLWK 4 !<br />
„ 4XDUN JOXRQ VWULQJ PRGHO IRU UHDO DQG YLUWXDO<br />
JDPPD UHDFWLRQV<br />
„ $XJHU HOHFWURQV HOHFWURQV<br />
DGGHG
:KDW LV QHZ LQ JHDQW IRU KDGURQLF<br />
SK\VLFV<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
'DWD<br />
„ 1HZ SKRWRQ HYDSRUDWLRQ GDWD DYDLODEOH<br />
„ 1HZ UDGLRDFWLYH UDGLRDFWLYH<br />
GHFD\ GDWD DYDLODEOH<br />
„ 1HXWURQ WUDQVSRUW FRQVLVWHQW ZLWK * 1'/<br />
WKURXJK<br />
3DUDPHWUL]HG PHWUL]HG PRGHO<br />
„ /RZ HQHUJ\ SDUDPHWUL]HG SDUDPHWUL]H PRGHOV %HWWHU HQHUJ\<br />
FRQVHUYDWLRQ IRU DQWL SDUWLFOHV<br />
„ 0DUV UHZULWH 1HXWURQ 1HXWURQ<br />
VSHFWUD H[WHQG EHORZ<br />
0H9<br />
„ 'DWD<br />
„ 3DUD
:KDW LV QHZ LQ JHDQW IRU KDGURQLF<br />
SK\VLFV<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ &URVV VHFWLRQV VHFWLRQV<br />
„ SLRQ FURVV VHFWLRQV $ QHZ GDWDVHW FODVV<br />
„ SURWHFWLRQ IRU ORZ QXFOHDU PDVVHV<br />
„ SURWHFWLRQ IRU HOHFWUR QXFOHDU IL[HV &06 SUREOHP<br />
„ PDQDJHPHQW<br />
„ PRGHO FODVVHV FODVVHV<br />
JHW GHOHWHG QR PDWWHU ZKDW WKH<br />
SK\VLFV OLVW ORRNV OLNH<br />
„ WUDFLQJ RI RULJLQDWRU PRGHO QRZ SRVVLEOH<br />
„ 6WRSSLQJ<br />
„ QHXWULQR IODYRU IL[HG IRU VWRSSLQJ PX<br />
„ 2WKHU EXJIL[HV
Example: the binary cascade<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ $ GLPHQVLRQDO QXFOHDU PRGHO LV EXLOW XS DV IRU<br />
TXDUN JOXRQ VWULQJ PRGHO<br />
„ 7KH LPSDFW SDUDPHWHU LV VDPSOHG<br />
„ 7KH WLPH WR FORVHVW DSSURDFK WR WKH QXFOHRQV LQ<br />
WKH QXFOHXV LV FDOFXODWHG<br />
„ 5HDFWLRQV DUH VHOHFWHG DFFRUGLQJ WR KDGURQ<br />
QXFOHRQ FURVV VHFWLRQ<br />
„ 7KH V\VWHP LV WUDQVSRUWHG LQ WKH QXFOHDU SRWHQWLDO<br />
WR WKH WLPH RI WKH HDUOLHVW UHDFWLRQ RU GHFD\<br />
„ 3URGXFWV DUH FUHDWHG DQG VXEMHFW WR SDXOLV<br />
H[FOXVLRQ SULQFLSOH DFFHSWHG<br />
„ 7KH SURFHGXUH LV UHSHDWHG IRU DOO SURGXFWV<br />
DOORZLQJ IRU UH VFDWWHULQJ
Why binary cascade?<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ 7KH QDPH ELQDU\ FDVFDGH FRPHV IURP WKH<br />
IDFW WKDW RQO\ ELQDU\ FROOLVLRQV DQG GHFD\<br />
DUH FRQVLGHUHG OLNH<br />
pp<br />
→<br />
„ 1R GHWDLOV RQ WKH PDWKHPDWLFV EXW WKLV WKH<br />
QXFOHRQ DQG GHOWD UHVRQDQFHV WDNHQ LQWR<br />
FRQVLGHUDWLRQ DUH WKHVH<br />
1232, , ∆1600 ∆1600,<br />
, ∆1620 ∆1620,<br />
, ∆1700 ∆1700,<br />
, ∆1900 ∆1900,<br />
, ∆1905 ∆1905,<br />
, ∆1910 ∆1910,<br />
, ∆1920 ∆1920,<br />
, ∆1930 ∆1930,<br />
, ∆1950<br />
1400, , Ν1520 Ν1520,<br />
, Ν1535 Ν1535,<br />
, Ν1650 Ν1650,<br />
, Ν1675 Ν1675,<br />
, Ν1680 Ν1680,<br />
, Ν1700 Ν1700,<br />
, Ν1710 Ν1710,<br />
, Ν1720 Ν1720,<br />
, Ν1900 Ν1900,<br />
, Ν1990 Ν1990,<br />
,<br />
Ν2090 2090, , Ν2190 Ν2190,<br />
, Ν2220 Ν2220,<br />
, Ν2250<br />
„ ∆1232<br />
„ Ν1400<br />
∆<br />
1232<br />
N<br />
*<br />
2190<br />
or<br />
K<br />
− p<br />
→<br />
Σ<br />
0<br />
1775
The potentials, transport equation, cross-<br />
section calculation, pauli-principle, pauli principle, etc..<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ $OO QHHGHG RI FRXUVH EXW IRU WKH VDNH RI<br />
EULHIQHVV MXVW D IHZ FRPSDULVRQV ZLWK<br />
H[LVWLQJ GDWD
Proton induced reactions cross-<br />
sections<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Binary cascade<br />
prediction<br />
Sample the<br />
Impact parameter over<br />
A large area.<br />
Make the ratio of ‘hits’<br />
To trials, times the area<br />
sampled<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Forward peaks in proton<br />
scattering (256 MeV) MeV<br />
Beryllium<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
Aluminum<br />
Iron<br />
Lead
One verification example:<br />
597 MeV p on Pb (PRC 22, p1184)<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
30 degrees<br />
60 degrees<br />
90 degrees<br />
120 degrees<br />
150 degrees<br />
Neutron production<br />
At 30, 60, 90, 120<br />
And 150 degrees
One verification example:<br />
597 MeV p on Pb (PRC 22, p1184)<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
Pi+ production at<br />
22.5, 45, 60, 90,<br />
And 130 degrees
One verification example:<br />
597 MeV p on Pb (PRC 22, p1184)<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
Pi- production at<br />
22.5, 45, 60, 90,<br />
And 130 degrees
Model validation<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ )RXU WLHU VWUDWHJ\<br />
$XWKRU YDOLGDWLRQ SORWV IRU WKH LQGLYLGXDO PRGHOV<br />
„ 3UHFRQGLWLRQ IRU PRGHO WR EH D FDQGLGDWH IRU LQFOXVLRQ<br />
„ $XWKRU YDOLGDWLRQ SORWV IRU WKH LQGLYLGXDO PRGHOV<br />
3UHFRQGLWLRQ IRU PRGHO WR EH D FDQGLGDWH IRU LQFOXVLRQ<br />
„ ,QGHSHQGHQW YDOLGDWLRQ RQ WKLQ WDUJHW GDWD ZLWK<br />
UHJUHVVLRQ VXLWHV E\ WKH ZRUNLQJ JURXSV<br />
„ 9HULILHG EHIRUH HYHU\ UHOHDVH<br />
„ ,QGHSHQGHQW YDOLGDWLRQ RQ EHQFKPDUNV ZKHUH<br />
WKHVH DUH DYDLODEOH<br />
„ 9HULILHG EHIRUH HYHU\ UHOHDVH<br />
„ 9DOLGDWLRQ RQ IXOO VLPXODWLRQ SURJUDPV<br />
„ JHDQW WDNHV PRGHO YDOLGDWLRQ PXFK PRUH<br />
VHULRXVO\ WKDQ LW ZDV LQ WKH WLPHV RI JHDQW
Author validation<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ $XWKRU YDOLGDWLRQ<br />
„ &RPSDULVRQV W\SLFDOO\ ZLWK PHDVXUHPHQWV IURP<br />
WKLQ WDUJHW GDWD<br />
„ /RRNLQJ DW SDUWLFOH \LHOGV DQG GLVWULEXWLRQV HWD DQG<br />
SW GLVWULEXWLRQV LQYDULDQW FURVV VHFWLRQV [I<br />
GLVWULEXWLRQV SDUWLFOH UDWLRQV HWF<br />
„ 5HTXHVWHG E\ WKH ZRUNLQJ JURXS ZKHQ PD\RU<br />
FKDQJHV WR D PRGHO RFFXU<br />
„ 2ZQHG E\ WKH DXWKRU OLNH WKH WHVW EHDP UHVXOW RI<br />
DQ H[SHULPHQWDO JURXS
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
Working group validation<br />
„ :RUNLQJ JURXS YDOLGDWLRQ VXLWHV<br />
„ )RU HWD SW [I PXOW G σ G S Gσ G7 QBSURQJ<br />
FKDUJH UDWLRV Gσ GΩG( G( HWF LQ SODFH IRU WKH<br />
YDULRXV HQHUJ\ UHJLPHV ,V DOUHDG\ TXLWH<br />
VDWLVIDFWRU\<br />
„ 7ULYLDO TXDQWLWLHV DOVR DUH FKHFNHG<br />
„ 1RWH WKDW WKLV FDQ EH GRQH RQO\ ZLWK WKH<br />
FRQVHQW RI WKH DXWKRU<br />
„ 7KLV OHYHO RI YDOLGDWLRQ ZDV QHYHU SHUIRUPHG LQ<br />
DQ\ GHSWK IRU JHDQW<br />
„ $W SUHVHQW DERXW SORWV
Anti proton annihilation<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Anti proton annihilation<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Stopping pion minus<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Low energy neutron capture<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Doppler broadening (low energy neutrons)<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Neutron induced isotope production<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Isotope production<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Proton induced reactions<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Example WG test<br />
results<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Test-beams Test beams<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
+DGURQLF WHVW EHDP FRPSDULVRQV FRPH IURP FROODERUDWLRQ RI<br />
H[SHULPHQWV GHWHFWRU JURXSV ZLWK µFRUH JHDQW SHUVRQQHO<br />
„ +DGURQLF<br />
„ $7/$6 7LOH WHVW EHDP<br />
„ &06 7LOH WHVW EHDP<br />
„ $7/$6 +(& WHVW EHDP<br />
„ $7/$6 )&$/ WHVW EHDP<br />
„ %7(9 FU\VWDO WHVW EHDP<br />
„ &06 FRPELQHG WHVW EHDP<br />
„ &V, WHVW EHDP EHQFKPDUN<br />
„ */$67 VWDUWLQJ WHVW EHDP<br />
„ + IRUZDUG EDUUHO VWDUWLQJ<br />
VWDUWLQJ<br />
„ $7/$6 FRPELQHG HQGFDS VWDUWLQJ
7HVW EHDP VDPSOH UHVXOW<br />
D FRXUWHV\ RI WKH GHWHFWRU JURXSV<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
ATLAS TILE prelim. CMS combined prelim.
A few recent test-beam test beam results (courtesy of<br />
ATLAS)<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
What can we profit from 5.0 for<br />
physics list updates<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ 0RVW RI WKH XSJUDGH ZLOO KDYH LWV HIIHFW<br />
EXW LQ SDUWLFXODU<br />
„ 7KH FDVFDGH FRGHV<br />
„ 7KH SLRQ FURVV VHFWLRQ UHYLHZ<br />
„ 'HSHQGLQJ RQ WKH DQDO\VLV WKH QHZ<br />
HOHFWUR QXFOHDU SK\VLFV<br />
„ ([DPSOH ,PSDFW RI WKH FURVV VHFWLRQ VHFWLRQ<br />
UHYLHZ RQ ORQJLWXGLQDO VKRZHU VKDSHV
Studying physics-lists physics lists updates for G4 5.0<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ 1HZ PRGHOLQJ SRVVLELOLWLHV DUH FXUUHQWO\<br />
XQGHU LQYHVWLJDWLRQ IRU LQFOXVLRQ LQWR<br />
HGXFDWHG JXHVV SK\VLFV OLVWV IRU JHDQW<br />
„ $ SUHYLHZ RQ FXUUHQW ZRUN ,PSDFW RI<br />
WKH FURVV VHFWLRQ VHFWLRQ<br />
UHYLHZ RQ ORQJLWXGLQDO<br />
VKRZHU VKDSHV
HEC G4 5.0 (true geometry, my toy analysis),<br />
what for the rest?<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
HEC G4 5.0 (true geometry, my toy analysis),<br />
what for the rest?<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
HEC shower shapes G4 5.0 (true geometry, my toy analysis)<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT
Summary<br />
J.P. <strong>Wellisch</strong>,<br />
CERN/EP/SFT<br />
„ :H KDYH FUHDWHG ZKDW LV WKH PRVW FRPSOHWH<br />
KDGURQLF VLPXODWLRQ HQJLQH LQ WKH ZRUOG<br />
„ ,W LV XVHG LQ HDUQHVW LQ PDQ\ +(3 DQG QRQ<br />
+(3 DSSOLFDWLRQV<br />
„ *($17 YHULILFDWLRQ KDV EHFRPH WHHWK<br />
„ ,Q WKH TXDQWLWLHV VWXGLHG *($17 KDGURQLF<br />
SK\VLFV LV HTXLYDOHQW RU EHWWHU WKDQ<br />
*($17<br />
„ 'HYHORSPHQW ZLOO FRQWLQXH DV PRUH DQG<br />
PRUH XVH FDVHV XQIROG