Electro and gamma nuclear physics in geant4. P ... - Geant4 - CERN
Electro and gamma nuclear physics in geant4. P ... - Geant4 - CERN
Electro and gamma nuclear physics in geant4. P ... - Geant4 - CERN
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
<strong>Electro</strong> <strong>and</strong> <strong>gamma</strong> <strong>nuclear</strong> <strong>physics</strong> <strong>in</strong><br />
<strong>geant4.</strong><br />
P. Degtyarenko, M. Kossov, J.P. Wellisch<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT
Outl<strong>in</strong>e<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
&URVV VHFWLRQ FDOFXODWLRQV<br />
„ &URVV<br />
„ 7KH PRGHOLQJ DQG UHVWULFWLRQV RI<br />
DSSOLFDELOLW\<br />
„ )LQDO VWDWH JHQHUDWLRQ<br />
„ 7KH PRGHOLQJ DQG UHVWULFWLRQV RI<br />
DSSOLFDELOLW\<br />
„ 9HULILFDWLRQ
Gamma <strong>nuclear</strong> reaction cross--sections<br />
cross sections<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ 0RGHOLQJ WKH IROORZLQJ UHJLRQV<br />
*LDQW'LSROH5HVRQDQFH UHJLPH<br />
„ *LDQW'LSROH5HVRQDQFH<br />
„ 2 0H9<br />
„ 5RSHU UHJLPH<br />
„ WKH LQWHUPHGLDWH HQHUJ\ µGHVHUW<br />
„ 'HOWD UHJLPH<br />
„ 2 *H9<br />
„ 5HJJHRQ SRPHURQ UHJLPH
GDR GDR regime<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ *'5 IURP SRZHU ODZ FKLUDO LQYDULDQW SKDVH<br />
VSDFH GHFD\ SUHGLFWLRQ DQG QXFOHDU EDUULHU<br />
UHIOHFWLRQ IXQFWLRQ<br />
GDR ( e,<br />
p,<br />
b,<br />
c,<br />
s)<br />
= T ( e,<br />
b,<br />
s)<br />
exp( c ! pe)<br />
T ( e,<br />
b,<br />
s)<br />
=<br />
1<br />
1+<br />
exp( ( b ! e)<br />
„ +HUH H LV ORJ (BJ<br />
„ 3DUDPHWHUV S E F V DUH WXQHG RQ<br />
H[SHULPHQWDO GDWD )URP +H WR 8 + DQG<br />
GHXWHULXP DUH WUHDWHG DV VSHFLDO FDVHV<br />
s)
Delta Delta isobar isobar region region<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ 7KH GHOWD LVREDU UHJLRQ LV PRGHOHG DV D %UHLW<br />
:LHJQHU IXQFWLRQ ZLWK D SURGXFWLRQ WKUHVKROG<br />
#<br />
( e,<br />
d,<br />
d " T ( e,<br />
f , g)<br />
f , g,<br />
r,<br />
q)<br />
=<br />
2<br />
1+<br />
r " ( e ! q)<br />
„ +HUH T FDQ EH YLHZHG DV WKH SRVLWLRQ RI WKH<br />
GHOWD UHVRQDQFH DQG U DV LWV LQYHUVH ZLGWK<br />
„ 7KH SDUDPHWHUV DUH WXQHG RQ EHIRUH<br />
PHQWLRQHG H[SHULPHQWDO GDWD DV D IXQFWLRQ RI<br />
ORJ $
Roper Roper region region<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ 7KLV UHJLPH ZDV SDUDPHWUL]HG XVLQJ WKH<br />
VDPH IXQFWLRQDO IRUP DV IRU WKH 'HOWD<br />
UHJLRQ EXW GURSSLQJ WKH SLRQ WKUHVKROG<br />
IDFWRU<br />
v<br />
Ro(<br />
e,<br />
v,<br />
w,<br />
u)<br />
=<br />
1+<br />
w"<br />
( e ! u)<br />
2
Reggeon--Pomeron<br />
Reggeon Pomeron region region<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ ,Q WKH UHJJHRQ SRPHURQ UHJLRQ ZH XVH<br />
RP(<br />
e,<br />
h)<br />
=<br />
h"<br />
T ( e,<br />
7.<br />
0,<br />
0.<br />
2)<br />
"<br />
h<br />
=<br />
A"<br />
exp( ! log( A)<br />
"<br />
( 0.<br />
885<br />
( 0.<br />
0116<br />
+<br />
"<br />
0.<br />
0048<br />
exp( 0.<br />
16<br />
" log( A)))<br />
" e)<br />
+<br />
0.<br />
4<br />
" exp( ! 0.<br />
2"<br />
e)),<br />
with
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT
Implementation restrictions:<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ 1RQH<br />
„ 'HWDLOV DUH FXUUHQWO\ EHLQJ SUHSDUHG IRU<br />
SXEOLFDWLRQ LQ D UHIHUHHG MRXUQDO
F<strong>in</strong>al F<strong>in</strong>al state state generation:<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ (QHUJLHV DERYH *H9 .DLGDORYV TXDUN<br />
JOXRQ VWULQJ PRGHO VHH WDON SUHVHQWHG<br />
E\ *XQWHU )ROJHU<br />
„ (QHUJLHV EHORZ *H9 &KLUDO LQYDULDQW<br />
SKDVH VSDFH GHD\
Chiral Invariant Phase--space Phase space Decay.<br />
„ $ TXDUN OHYHO GLPHQVLRQDO HYHQW JHQHUDWRU IRU<br />
„<br />
IUDJPHQWDWLRQ RI H[FLWHG KDGURQLF V\VWHPV<br />
TXDVPRQV LQWR KDGURQV<br />
%DVHG RQ WKH 4&' LGHD RI DV\PSWRWLF IUHHGRP<br />
„ /RFDO FKLUDO LQYDULDQFH UHVWRUDWLRQ OHWV XV FRQVLGHU<br />
TXDUN SDUWRQV PDVVOHVV DQG ZH FDQ LQWHJUDWH WKH<br />
LQYDULDQW SKDVH VSDFH GLVWULEXWLRQ RI TXDUN<br />
SDUWRQV DQG TXDUN H[FKDQJH IXVLRQ PHFKDQLVP<br />
RI KDGURQL]DWLRQ<br />
„ 7KH RQO\ QRQ NLQHPDWLFDO FRQFHSW XVHG LV WKDW RI D<br />
WHPSHUDWXUH RI WKH TXDVPRQ<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT
Vacuum CHIPS<br />
„ 2U KRZ WR FDOFXODWH WKH GHFD\ RI IUHH H[FLWHG<br />
KDGURQLF V\VWHPV<br />
„ ,Q DQ ILQLWH WKHUPDOL]HG V\VWHP RI 1 SDUWRQV ZLWK<br />
WRWDO PDVV 0 WKH LQYDULDQW SKDVH VSDFH LQWHJUDO LV<br />
2 N ! 4<br />
SURSRUWLRQDO WR M DQG WKH VWDWLVWLFDO<br />
! M / T<br />
GHQVLW\ RI VWDWHV LV SURSRUWLRQDO WR e +HQFH<br />
ZH FDQ ZULWH WKH SUREDELOLW\ WR ILQG 1 SDUWRQV ZLWK<br />
WHPSHUDWXUH 7 LQ D VWDWH ZLWK PDVV 0 DV<br />
dW<br />
„ 1RWH WKDW IRU WKLV GLVWULEXWLRQ WKH PHDQ PDVV<br />
VTXDUH LV<br />
2<br />
M<br />
2<br />
=<br />
2N ( 2N<br />
! 2)<br />
T<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
$<br />
M<br />
e<br />
2 N ! 4 ! M / T<br />
dM
Vacuum CHIPS<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ :H XVH WKLV IRUPXOD WR FDOFXODWH WKH QXPEHU<br />
RI SDUWRQV LQ WKH TXDVPRQ DQG REWDLQ WKH<br />
SDUWRQ VSHFWUXP<br />
dW<br />
kdk<br />
* 2k<br />
'<br />
$ ( 1!<br />
%<br />
) M &<br />
„ 7R REWDLQ WKH SUREDELOLW\ IRU TXDUN IXVLRQ LQWR<br />
KDGURQV ZH FDQ FRPSXWH WKH SUREDELOLW\ WR<br />
ILQG WZR SDUWRQV ZLWK PRPHQWD T DQG N ZLWK<br />
WKH LQYDULDQW PDVV µ<br />
N<br />
! 4<br />
+ % *<br />
' *<br />
'<br />
(<br />
2q<br />
% ( 2 2kq(<br />
1!<br />
cos,<br />
)<br />
P(<br />
k,<br />
M , µ ) = 1!<br />
× - µ !<br />
qdqd cos,<br />
(<br />
% (<br />
) M 1!<br />
2k<br />
M & ) 1!<br />
2k<br />
M &<br />
N ! 3
Vacuum CHIPS<br />
„ 8VLQJ WKH GHOWD IXQFWLRQ WR SHUIRUP WKH LQWHJUDWLRQ<br />
DQG WKH PDVV FRQVWUDLQW ZH ILQG WKH WRWDO<br />
NLQHPDWLFDO SUREDELOLW\ RI KDGURQL]DWLRQ RI D SDUWRQ<br />
ZLWK PRPHQWXP N LQWR D KDGURQ ZLWK PDVV µ:<br />
„ $FFRXQWLQJ IRU VSLQ DQG TXDUN FRQWHQW RI WKH ILQDO<br />
VWDWH KDGURQ DGGV IDFWRU<br />
V DQG D FRPELQDWRULDO<br />
„ $W WKLV OHYHO RI WKH ODQJXDJH &+,36 FDQ EH DSSOLHG<br />
WR S SEDU DQQLKLODWLRQ<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
M<br />
4k(<br />
N<br />
! N<br />
2 k<br />
!<br />
! 3)<br />
( ) 3<br />
2<br />
1!<br />
µ<br />
2kM
For For more more <strong>in</strong>formation see: see:<br />
Eur.Phys.Journal Eur.Phys.Journal<br />
A8,217(2000)<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT
Nuclear CHIPS<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ ,Q RUGHU WR DSSO\ &+,36 IRU DQ H[FLWHG KDGURQLF<br />
V\VWHP ZLWKLQ QXFOHL ZH KDYH WR DGG SDUWRQ<br />
H[FKDQJH ZLWK QXFOHDU FOXVWHUV WR WKH PRGHO<br />
„ 7KH NLQHPDWLFDO SLFWXUH LV WKDW D FRORU QHXWUDO<br />
TXDVPRQ HPLWV D SDUWRQ ZKLFK LV DEVRUEHG E\ D<br />
QXFOHRQ RU D QXFOHDU FOXVWHU 7KLV UHVXOWV LQ D<br />
FRORUHG UHVLGXDO TXDVPRQ DQG D FRORUHG<br />
FRPSRXQG<br />
„ 7KH FRORUHG FRPSRXQG WKHQ GHFD\V LQWR DQ<br />
RXWJRLQJ QXFOHDU IUDJPHQW DQG D µUHFRLO TXDUN<br />
WKDW LV LQFRUSRUDWHG E\ WKH FRORUHG TXDVPRQ
The The parton exchange diagram<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT
Nuclear CHIPS<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ $SSO\LQJ PHFKDQLVPV DQDORJXH WR YDFXXP &+,36<br />
ZH FDQ ZULWH WKH SUREDELOLW\ RI HPLVVLRQ RI D<br />
QXFOHDU IUDJPHQW ZLWK PDVV µ DV D UHVXOW RI WKH<br />
WUDQVLWLRQ RI D SDUWRQ ZLWK PRPHQWXP N IURP WKH<br />
TXDVPRQ WR D IUDJPHQW ZLWK PDVV µ DV<br />
*<br />
'<br />
(<br />
2(<br />
k ! # )<br />
%<br />
µ ’ ( k ! # )<br />
P ( k,<br />
µ ’ , µ ) = + 1!<br />
d cos,<br />
(<br />
% [ + ! ] 2<br />
) µ ’ + k(<br />
1!<br />
cos,<br />
kq ) & 2 µ ’ k(<br />
1 cos,<br />
kq )<br />
„ +HUH Q LV WKH QXPEHU RI TXDUN SDUWRQV LQ WKH<br />
QXFOHDU FOXVWHU DQG # LV WKH FRYDULDQW ELQGLQJ<br />
HQHUJ\ RI WKH FOXVWHU DQG WKH LQWHJUDO LV RYHU WKH<br />
DQJOH EHWZHHQ SDUWRQ DQG UHFRLO SDUWRQ<br />
n!<br />
3<br />
kq
Nuclear CHIPS<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ 7R FDOFXODWH WKH IUDJPHQW \LHOGV LW LV<br />
QHFHVVDU\ WR FDOFXODWH WKH SUREDELOLW\ WR ILQG D<br />
FOXVWHU RI . QXFOHRQV ZLWKLQ D QXFOHXV :H GR<br />
WKLV XVLQJ WKH IROORZLQJ DVVXPSWLRQV<br />
„ $ IUDFWLRQ /1 RI DOO QXFOHRQV LV QRW FOXVWHULVLQJ<br />
„ $ IUDFWLRQ /2 RI WKH QXFOHRQV LQ WKH SHULSKHU\ RI<br />
WKH QXFOHXV LV FOXVWHULQJ LQWR WZR QXFOHRQ FOXVWHUV<br />
„ 7KHUH LV D VLQJOH FOXVWHUL]DWLRQ SUREDELOLW\ 0<br />
„ DQG ILQG ZLWK D EHLQJ WKH QXPEHU RI<br />
QXFOHRQV LQYROYHG LQ FOXVWHUL]DWLRQ<br />
a ! 1<br />
C 0<br />
P = a!<br />
1<br />
( 1+<br />
0)<br />
.<br />
.<br />
.
Nuclear CHIPS<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ $W WKLV OHYHO RI WKH ODQJXDJH &+,36 FDQ<br />
EH DSSOLHG WR UHDFWLRQV RI SLRQV DQG<br />
SKRWR QXFOHDU UHDFWLRQV
A A few few results: For For more more <strong>in</strong>formation see: see:<br />
Eur.Phys.Journal Eur.Phys.Journal<br />
A9,411(2000) <strong>and</strong> <strong>and</strong><br />
Eur.Phys.Journal Eur.Phys.Journal<br />
A9,421(2000)<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT
<strong>Electro</strong>--<strong>nuclear</strong> <strong>Electro</strong> <strong>nuclear</strong> scatter<strong>in</strong>g:<br />
reaction cross--sections<br />
cross sections<br />
„ %DVHG RQ )HUPLV PHWKRG RI HTXLYDOHQW<br />
SKRWRQV DV GHYHORSHG E\ :HLW]VDHFNHU DQG<br />
:LOOLDPV<br />
dN<br />
1<br />
23<br />
= ! log( y)<br />
d log( y),<br />
2<br />
„ )ROGLQJ WKLV IOX[ ZLWK WKH JDPPD UHDFWLRQ<br />
FURVV VHFWLRQ DV GHVFULEHG DERYH DQG<br />
LQWHJUDWLQJ WKH JDPPD VSHFWUXP ZH REWDLQ<br />
23<br />
23<br />
( ( . ) d log( . )<br />
4 eA 5<br />
X ) = log( E ) + e 4 1A(<br />
. ) d log( . ) ! + log( . ) 4 1A<br />
J.P. Wellisch,<br />
2<br />
2<br />
<strong>CERN</strong>/EP/SFT<br />
y<br />
=<br />
.<br />
E<br />
e
Restrictions<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ 7KH PRGHOLQJ DVVXPSWLRQV LQ WKH<br />
HTXLYDOHQW SKRWRQ VSHFWUXP<br />
„ 7KH ',6 SDUW ZDV UHFHQWO\ DGGHG<br />
„ 7R VPRRWKHQ WKH HOHFWUR QXFOHDU FURVV<br />
VHFWLRQ DW *H9 WKH IROORZLQJ 35 WHUP<br />
ZDV XVHG<br />
RP ( e,<br />
h)<br />
=<br />
h " ( 0.<br />
0116"<br />
exp( 0.<br />
16"<br />
e)<br />
+ 1.<br />
0"<br />
exp( ! 0.<br />
26"<br />
e))
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
<strong>Electro</strong>--<strong>nuclear</strong> <strong>Electro</strong> <strong>nuclear</strong> scatter<strong>in</strong>g: f<strong>in</strong>al f<strong>in</strong>al<br />
state state model<strong>in</strong>g<br />
„ 7KH IRUPXODV SUHVHQWHG LQ WKH FURVV VHFWLRQ<br />
VHFWLRQ FDQ EH XVHG WR FDOFXODWH WKH<br />
SUREDELOLW\ GLVWULEXWLRQ RI WKH HTXLYDOHQW<br />
SKRWRQV<br />
„ 7KLV GLVWULEXWLRQ LV VDPSOHG WR FDOFXODWH WKH<br />
HQHUJ\ WUDQVIHU RI WKH QXFOHDU UHDFWLRQ<br />
„ 7KH HQHUJ\ LV DVVXPHG WR EH WUDQVIHUUHG E\<br />
JDPPD H[FKDQJH DQG WKH PRGHOV IRU JDPPD<br />
QXFOHDU UHDFWLRQ &+,36 4*6WULQJ DUH XVHG
Implementation restrictions<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ $VVXPSWLRQ RI HTXLYDOHQW SKRWRQ<br />
„ 2WKHUZLVH QRQH<br />
„ )RU HQHUJ\ WUDQVIHUV EHORZ *H9 SHU<br />
FROOLVLRQ &+,36 LV XVHG<br />
„ DERYH TXDUN JOXRQ VWULQJ PRGHO LV XVHG<br />
„ ',6 LV LQFOXGHG
Hard Hard scatter<strong>in</strong>g <strong>in</strong> <strong>in</strong> electro--<strong>nuclear</strong><br />
electro <strong>nuclear</strong><br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT
Hard Hard scatter<strong>in</strong>g <strong>in</strong> <strong>in</strong> electro--<strong>nuclear</strong><br />
electro <strong>nuclear</strong><br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT
Conclusions<br />
J.P. Wellisch,<br />
<strong>CERN</strong>/EP/SFT<br />
„ *HDQW SK\VLFV PRGHOLQJ LV QRZ<br />
HTXLSSHG WR GHDO ZLWK OLQHDU FROOLGHU RU<br />
-HIIHUVRQ /DE UDGLDWLRQ HQYLURQPHQW<br />
DQG /+& SUHFLVLRQ JDPPD LGHQWLILFDWLRQ<br />
XVH FDVHV