Electro and gamma nuclear physics in geant4. P ... - Geant4 - CERN

Electro and gamma nuclear physics in 

geant4. 

P. Degtyarenko, M. Kossov, J.P. Wellisch 

J.P. Wellisch, 

CERN/EP/SFT 


CERN/EP/SFT

Outline 



&URVV VHFWLRQ FDOFXODWLRQV 

„ &URVV 

„ 7KH PRGHOLQJ DQG UHVWULFWLRQV RI 

DSSOLFDELOLW\ 

„ )LQDO VWDWH JHQHUDWLRQ 

„ 7KH PRGHOLQJ DQG UHVWULFWLRQV RI 

DSSOLFDELOLW\ 

„ 9HULILFDWLRQ

Gamma nuclear reaction cross--sections 

cross sections 



„ 0RGHOLQJ WKH IROORZLQJ UHJLRQV 

*LDQW'LSROH5HVRQDQFH UHJLPH 

„ *LDQW'LSROH5HVRQDQFH 

„ 2 0H9 

„ 5RSHU UHJLPH 

„ WKH LQWHUPHGLDWH HQHUJ\ µGHVHUW 

„ 'HOWD UHJLPH 

„ 2 *H9 

„ 5HJJHRQ SRPHURQ UHJLPH

GDR GDR regime 



„ *'5 IURP SRZHU ODZ FKLUDO LQYDULDQW SKDVH 

VSDFH GHFD\ SUHGLFWLRQ DQG QXFOHDU EDUULHU 

UHIOHFWLRQ IXQFWLRQ 

GDR ( e, 

p, 

b, 

c, 

s) 

= T ( e, 

b, 

s) 

exp( c ! pe) 

T ( e, 

b, 

s) 

= 

1 

1+ 

exp( ( b ! e) 

„ +HUH H LV ORJ (BJ 

„ 3DUDPHWHUV S E F V DUH WXQHG RQ 

H[SHULPHQWDO GDWD )URP +H WR 8 + DQG 

GHXWHULXP DUH WUHDWHG DV VSHFLDO FDVHV 

s)

Delta Delta isobar isobar region region 



„ 7KH GHOWD LVREDU UHJLRQ LV PRGHOHG DV D %UHLW 

:LHJQHU IXQFWLRQ ZLWK D SURGXFWLRQ WKUHVKROG 

# 

( e, 

d, 

d " T ( e, 

f , g) 

f , g, 

r, 

q) 

= 

2 

1+ 

r " ( e ! q) 

„ +HUH T FDQ EH YLHZHG DV WKH SRVLWLRQ RI WKH 

GHOWD UHVRQDQFH DQG U DV LWV LQYHUVH ZLGWK 

„ 7KH SDUDPHWHUV DUH WXQHG RQ EHIRUH 

PHQWLRQHG H[SHULPHQWDO GDWD DV D IXQFWLRQ RI 

ORJ $

Roper Roper region region 



„ 7KLV UHJLPH ZDV SDUDPHWUL]HG XVLQJ WKH 

VDPH IXQFWLRQDO IRUP DV IRU WKH 'HOWD 

UHJLRQ EXW GURSSLQJ WKH SLRQ WKUHVKROG 

IDFWRU 

v 

Ro( 

e, 

v, 

w, 

u) 

= 

1+ 

w" 

( e ! u) 

2

Reggeon--Pomeron 

Reggeon Pomeron region region 



„ ,Q WKH UHJJHRQ SRPHURQ UHJLRQ ZH XVH 

RP( 

e, 

h) 

= 

h" 

T ( e, 

7. 

0, 

0. 

2) 

" 

h 

= 

A" 

exp( ! log( A) 

" 

( 0. 

885 

( 0. 

0116 

+ 

" 

0. 

0048 

exp( 0. 

16 

" log( A))) 

" e) 

+ 

0. 

4 

" exp( ! 0. 

2" 

e)), 

with



Implementation restrictions: 



„ 1RQH 

„ 'HWDLOV DUH FXUUHQWO\ EHLQJ SUHSDUHG IRU 

SXEOLFDWLRQ LQ D UHIHUHHG MRXUQDO

Final Final state state generation: 



„ (QHUJLHV DERYH *H9 .DLGDORYV TXDUN 

JOXRQ VWULQJ PRGHO VHH WDON SUHVHQWHG 

E\ *XQWHU )ROJHU 

„ (QHUJLHV EHORZ *H9 &KLUDO LQYDULDQW 

SKDVH VSDFH GHD\

Chiral Invariant Phase--space Phase space Decay. 

„ $ TXDUN OHYHO GLPHQVLRQDO HYHQW JHQHUDWRU IRU 

„ 

IUDJPHQWDWLRQ RI H[FLWHG KDGURQLF V\VWHPV 

TXDVPRQV LQWR KDGURQV 

%DVHG RQ WKH 4&' LGHD RI DV\PSWRWLF IUHHGRP 

„ /RFDO FKLUDO LQYDULDQFH UHVWRUDWLRQ OHWV XV FRQVLGHU 

TXDUN SDUWRQV PDVVOHVV DQG ZH FDQ LQWHJUDWH WKH 

LQYDULDQW SKDVH VSDFH GLVWULEXWLRQ RI TXDUN 

SDUWRQV DQG TXDUN H[FKDQJH IXVLRQ PHFKDQLVP 

RI KDGURQL]DWLRQ 

„ 7KH RQO\ QRQ NLQHPDWLFDO FRQFHSW XVHG LV WKDW RI D 

WHPSHUDWXUH RI WKH TXDVPRQ 



Vacuum CHIPS 

„ 2U KRZ WR FDOFXODWH WKH GHFD\ RI IUHH H[FLWHG 

KDGURQLF V\VWHPV 

„ ,Q DQ ILQLWH WKHUPDOL]HG V\VWHP RI 1 SDUWRQV ZLWK 

WRWDO PDVV 0 WKH LQYDULDQW SKDVH VSDFH LQWHJUDO LV 

2 N ! 4 

SURSRUWLRQDO WR M DQG WKH VWDWLVWLFDO 

! M / T 

GHQVLW\ RI VWDWHV LV SURSRUWLRQDO WR e +HQFH 

ZH FDQ ZULWH WKH SUREDELOLW\ WR ILQG 1 SDUWRQV ZLWK 

WHPSHUDWXUH 7 LQ D VWDWH ZLWK PDVV 0 DV 

dW 

„ 1RWH WKDW IRU WKLV GLVWULEXWLRQ WKH PHDQ PDVV 

VTXDUH LV 

2 

M 

2 

= 

2N ( 2N 

! 2) 

T 



$ 

M 

e 

2 N ! 4 ! M / T 

dM

Vacuum CHIPS 



„ :H XVH WKLV IRUPXOD WR FDOFXODWH WKH QXPEHU 

RI SDUWRQV LQ WKH TXDVPRQ DQG REWDLQ WKH 

SDUWRQ VSHFWUXP 

dW 

kdk 

* 2k 

' 

$ ( 1! 

% 

) M & 

„ 7R REWDLQ WKH SUREDELOLW\ IRU TXDUN IXVLRQ LQWR 

KDGURQV ZH FDQ FRPSXWH WKH SUREDELOLW\ WR 

ILQG WZR SDUWRQV ZLWK PRPHQWD T DQG N ZLWK 

WKH LQYDULDQW PDVV µ 

N 

! 4 

+ % * 

' * 

' 

( 

2q 

% ( 2 2kq( 

1! 

cos, 

) 

P( 

k, 

M , µ ) = 1! 

× - µ ! 

qdqd cos, 

( 

% ( 

) M 1! 

2k 

M & ) 1! 

2k 

M & 

N ! 3

Vacuum CHIPS 

„ 8VLQJ WKH GHOWD IXQFWLRQ WR SHUIRUP WKH LQWHJUDWLRQ 

DQG WKH PDVV FRQVWUDLQW ZH ILQG WKH WRWDO 

NLQHPDWLFDO SUREDELOLW\ RI KDGURQL]DWLRQ RI D SDUWRQ 

ZLWK PRPHQWXP N LQWR D KDGURQ ZLWK PDVV µ: 

„ $FFRXQWLQJ IRU VSLQ DQG TXDUN FRQWHQW RI WKH ILQDO 

VWDWH KDGURQ DGGV IDFWRU 

V DQG D FRPELQDWRULDO 

„ $W WKLV OHYHO RI WKH ODQJXDJH &+,36 FDQ EH DSSOLHG 

WR S SEDU DQQLKLODWLRQ 



M 

4k( 

N 

! N 

2 k 

! 

! 3) 

( ) 3 

2 

1! 

µ 

2kM

For For more more information see: see: 

Eur.Phys.Journal Eur.Phys.Journal 

A8,217(2000) 



Nuclear CHIPS 



„ ,Q RUGHU WR DSSO\ &+,36 IRU DQ H[FLWHG KDGURQLF 

V\VWHP ZLWKLQ QXFOHL ZH KDYH WR DGG SDUWRQ 

H[FKDQJH ZLWK QXFOHDU FOXVWHUV WR WKH PRGHO 

„ 7KH NLQHPDWLFDO SLFWXUH LV WKDW D FRORU QHXWUDO 

TXDVPRQ HPLWV D SDUWRQ ZKLFK LV DEVRUEHG E\ D 

QXFOHRQ RU D QXFOHDU FOXVWHU 7KLV UHVXOWV LQ D 

FRORUHG UHVLGXDO TXDVPRQ DQG D FRORUHG 

FRPSRXQG 

„ 7KH FRORUHG FRPSRXQG WKHQ GHFD\V LQWR DQ 

RXWJRLQJ QXFOHDU IUDJPHQW DQG D µUHFRLO TXDUN 

WKDW LV LQFRUSRUDWHG E\ WKH FRORUHG TXDVPRQ

The The parton exchange diagram 



Nuclear CHIPS 



„ $SSO\LQJ PHFKDQLVPV DQDORJXH WR YDFXXP &+,36 

ZH FDQ ZULWH WKH SUREDELOLW\ RI HPLVVLRQ RI D 

QXFOHDU IUDJPHQW ZLWK PDVV µ DV D UHVXOW RI WKH 

WUDQVLWLRQ RI D SDUWRQ ZLWK PRPHQWXP N IURP WKH 

TXDVPRQ WR D IUDJPHQW ZLWK PDVV µ DV 

* 

' 

( 

2( 

k ! # ) 

% 

µ ’ ( k ! # ) 

P ( k, 

µ ’ , µ ) = + 1! 

d cos, 

( 

% [ + ! ] 2 

) µ ’ + k( 

1! 

cos, 

kq ) & 2 µ ’ k( 

1 cos, 

kq ) 

„ +HUH Q LV WKH QXPEHU RI TXDUN SDUWRQV LQ WKH 

QXFOHDU FOXVWHU DQG # LV WKH FRYDULDQW ELQGLQJ 

HQHUJ\ RI WKH FOXVWHU DQG WKH LQWHJUDO LV RYHU WKH 

DQJOH EHWZHHQ SDUWRQ DQG UHFRLO SDUWRQ 

n! 

3 

kq

Nuclear CHIPS 



„ 7R FDOFXODWH WKH IUDJPHQW \LHOGV LW LV 

QHFHVVDU\ WR FDOFXODWH WKH SUREDELOLW\ WR ILQG D 

FOXVWHU RI . QXFOHRQV ZLWKLQ D QXFOHXV :H GR 

WKLV XVLQJ WKH IROORZLQJ DVVXPSWLRQV 

„ $ IUDFWLRQ /1 RI DOO QXFOHRQV LV QRW FOXVWHULVLQJ 

„ $ IUDFWLRQ /2 RI WKH QXFOHRQV LQ WKH SHULSKHU\ RI 

WKH QXFOHXV LV FOXVWHULQJ LQWR WZR QXFOHRQ FOXVWHUV 

„ 7KHUH LV D VLQJOH FOXVWHUL]DWLRQ SUREDELOLW\ 0 

„ DQG ILQG ZLWK D EHLQJ WKH QXPEHU RI 

QXFOHRQV LQYROYHG LQ FOXVWHUL]DWLRQ 

a ! 1 

C 0 

P = a! 

1 

( 1+ 

0) 

. 

. 

.

Nuclear CHIPS 



„ $W WKLV OHYHO RI WKH ODQJXDJH &+,36 FDQ 

EH DSSOLHG WR UHDFWLRQV RI SLRQV DQG 

SKRWR QXFOHDU UHDFWLRQV

A A few few results: For For more more information see: see: 


A9,411(2000) and and 


A9,421(2000) 





Electro--nuclear Electro nuclear scattering: 

reaction cross--sections 

cross sections 

„ %DVHG RQ )HUPLV PHWKRG RI HTXLYDOHQW 

SKRWRQV DV GHYHORSHG E\ :HLW]VDHFNHU DQG 

:LOOLDPV 

dN 

1 

23 

= ! log( y) 

d log( y), 

2 

„ )ROGLQJ WKLV IOX[ ZLWK WKH JDPPD UHDFWLRQ 

FURVV VHFWLRQ DV GHVFULEHG DERYH DQG 

LQWHJUDWLQJ WKH JDPPD VSHFWUXP ZH REWDLQ 

23 

23 

( ( . ) d log( . ) 

4 eA 5 

X ) = log( E ) + e 4 1A( 

. ) d log( . ) ! + log( . ) 4 1A 


2 

2 


y 

= 

. 

E 

e

Restrictions 



„ 7KH PRGHOLQJ DVVXPSWLRQV LQ WKH 

HTXLYDOHQW SKRWRQ VSHFWUXP 

„ 7KH ',6 SDUW ZDV UHFHQWO\ DGGHG 

„ 7R VPRRWKHQ WKH HOHFWUR QXFOHDU FURVV 

VHFWLRQ DW *H9 WKH IROORZLQJ 35 WHUP 

ZDV XVHG 

RP ( e, 

h) 

= 

h " ( 0. 

0116" 

exp( 0. 

16" 

e) 

+ 1. 

0" 

exp( ! 0. 

26" 

e))



Electro--nuclear Electro nuclear scattering: final final 

state state modeling 

„ 7KH IRUPXODV SUHVHQWHG LQ WKH FURVV VHFWLRQ 

VHFWLRQ FDQ EH XVHG WR FDOFXODWH WKH 

SUREDELOLW\ GLVWULEXWLRQ RI WKH HTXLYDOHQW 

SKRWRQV 

„ 7KLV GLVWULEXWLRQ LV VDPSOHG WR FDOFXODWH WKH 

HQHUJ\ WUDQVIHU RI WKH QXFOHDU UHDFWLRQ 

„ 7KH HQHUJ\ LV DVVXPHG WR EH WUDQVIHUUHG E\ 

JDPPD H[FKDQJH DQG WKH PRGHOV IRU JDPPD 

QXFOHDU UHDFWLRQ &+,36 4*6WULQJ DUH XVHG

Implementation restrictions 



„ $VVXPSWLRQ RI HTXLYDOHQW SKRWRQ 

„ 2WKHUZLVH QRQH 

„ )RU HQHUJ\ WUDQVIHUV EHORZ *H9 SHU 

FROOLVLRQ &+,36 LV XVHG 

„ DERYH TXDUN JOXRQ VWULQJ PRGHO LV XVHG 

„ ',6 LV LQFOXGHG

Hard Hard scattering in in electro--nuclear 

electro nuclear 



Hard Hard scattering in in electro--nuclear 

electro nuclear 



Conclusions 



„ *HDQW SK\VLFV PRGHOLQJ LV QRZ 

HTXLSSHG WR GHDO ZLWK OLQHDU FROOLGHU RU 

-HIIHUVRQ /DE UDGLDWLRQ HQYLURQPHQW 

DQG /+& SUHFLVLRQ JDPPD LGHQWLILFDWLRQ 

XVH FDVHV

Electro and gamma nuclear physics in geant4. P ... - Geant4 - CERN

Create successful ePaper yourself

Delete template?

Save as template?