Geant4 Hadronic Physics Working group progress and status. J.P. - - PowerPoint PPT Presentation

J.P. Wellisch, CERN/EP/CMC

Geant4 Hadronic Physics Working group progress and status.

J.P. Wellisch CERN/EP/CMC

J.P. Wellisch, CERN/EP/CMC

Outline Outline

J.P. Wellisch, CERN/EP/CMC

The dry numbers… The dry numbers…

J.P. Wellisch, CERN/EP/CMC

… …

J.P. Wellisch, CERN/EP/CMC

The milestones. The milestones.

• γ

γ

J.P. Wellisch, CERN/EP/CMC

Particle physics relevant requirements Particle physics relevant requirements collected collected

• Specifically from LHC:

Specifically from LHC:

• Have starting physics list (done)

Have starting physics list (done)

• Improve information flow on V&V (done, to be verified)

Improve information flow on V&V (done, to be verified)

• Provide a cascade code (two focused efforts)

Provide a cascade code (two focused efforts)

• Fix known problems in low energy models (released)

Fix known problems in low energy models (released)

J.P. Wellisch, CERN/EP/CMC

A complete sample requirement. A complete sample requirement.

• Name and E

Name and E-

• mail: Dennis Wright, SLAC

mail: Dennis Wright, SLAC

• Title:

Title: Pion Pion and and kaon kaon nuclear nuclear prodution prodution cross cross-

• sections

sections

• Description: 1

Description: 1-

• 5

5 GeV pions GeV pions and and kaons kaons interacting interacting inelastically

• inelastically. 10%

. 10% precision would be great. precision would be great.

• Rationale: Representation of how the shower develops.

Rationale: Representation of how the shower develops.

• Supporting use

Supporting use-

• cases that require this: Trying to model

cases that require this: Trying to model hadronic hadronic interactions in the interactions in the BaBar intsrumented BaBar intsrumented flux return. flux return.

• Responsible

Responsible cathegory cathegory: : hadronics hadronics

• Fulfillment criterion: Comparison to the data from particles int

Fulfillment criterion: Comparison to the data from particles interacting eracting in the beam in the beam-

• pipe and flux return.

pipe and flux return.

• Relevace

Relevace: very highly relevant for : very highly relevant for BaBar BaBar, relevant also for , relevant also for LHCb LHCb

• References: An E

References: An E-

• mail from Dennis pointing to the data.

mail from Dennis pointing to the data.

J.P. Wellisch, CERN/EP/CMC

Requirements collected Requirements collected – – titles only titles only

• During the last geant4 workshop, the users workshop at SLAC,

During the last geant4 workshop, the users workshop at SLAC, and in private mails: and in private mails:

• Ensure that the physics reference manual match the implementatio

Ensure that the physics reference manual match the implementation, n, and the models are mentioned in the applications developers guid and the models are mentioned in the applications developers guide, e, Referencing papers is just fine. Referencing papers is just fine.

• Memory usage for G4NDL cross

Memory usage for G4NDL cross-

• sections

sections

• Energy and momentum conservation should be checked in regression

Energy and momentum conservation should be checked in regression independently by the working group for all models, and publish t independently by the working group for all models, and publish the he test test-

• suite.

suite.

• Use ’well known’ international benchmarks to validate; publish r

Use ’well known’ international benchmarks to validate; publish results. esults.

• More understandable

More understandable hadronic hadronic physics lists physics lists

• Get documentation on which model is good/

Get documentation on which model is good/usefull usefull/required for which /required for which use use-

• case

case

J.P. Wellisch, CERN/EP/CMC

Requirements Cont. Requirements Cont.

• Provide a set of plots, and put them u; i.e. provide a place to

Provide a set of plots, and put them u; i.e. provide a place to put put these plots for users. these plots for users.

• Cross

Cross-

• sections for n,p inelastic scattering below 150

sections for n,p inelastic scattering below 150 MeV MeV in CMS in CMS tracker materials at 10% level of precisions. tracker materials at 10% level of precisions.

• Pion

Pion and and kaon kaon nuclear nuclear prodution prodution cross cross-

• sections: 1

sections: 1-

• 5

5 GeV pions GeV pions and and kaons kaons interacting interacting inelastically

• inelastically. 10% precision would be great.

. 10% precision would be great.

• 10

10-

• 20% level of description for 10

20% level of description for 10-

• 100GeV incident protons for

100GeV incident protons for example on Beryllium or copper. example on Beryllium or copper.

• Enroll a set of users to validate on complete application; as be

Enroll a set of users to validate on complete application; as beta ta testers so to say. testers so to say.

• Compare inclusive and exclusive cross

Compare inclusive and exclusive cross-

• sectios

sectios to data from the to data from the RAL/Durham database RAL/Durham database

J.P. Wellisch, CERN/EP/CMC

Requirements Cont. Requirements Cont.

• Description on how to set cuts, and its effects

Description on how to set cuts, and its effects

• Possibility to stop low energy neutral particles (like neutrons)

Possibility to stop low energy neutral particles (like neutrons)

• Each model should be specified concerning its application

Each model should be specified concerning its application area/use area/use-

• cases

cases

• Provide a list of models per use

Provide a list of models per use-

• case package

case package

• List of models per use

List of models per use-

• case package

case package

• Include physics list samples, once they exists, into the

Include physics list samples, once they exists, into the phsics phsics editor editor

• Ensure

Ensure tracability tracability of data to the primary source

• f data to the primary source
• Well defined process for updating the databases on request.

Well defined process for updating the databases on request.

• Parametrizations

Parametrizations of

• f hadronic

hadronic showers in showers in CsI CsI and Iron. and Iron.

• Parametrizations

Parametrizations of neutron background in LHC experiments

• f neutron background in LHC experiments
• Models for alpha incident inelastic reactions

Models for alpha incident inelastic reactions

J.P. Wellisch, CERN/EP/CMC

Requirements Cont. Requirements Cont.

• Models for alpha incident

Models for alpha incident difractive difractive dissociation reactions dissociation reactions

• Where the

Where the modelling modelling approach allows to produce the residual, it approach allows to produce the residual, it should be provided. should be provided.

• Neutron production by alphas at energies below 10

Neutron production by alphas at energies below 10 MeV MeV; including ; including reaction cross reaction cross-

• sections at 20%precision and kinematics of neutrons

sections at 20%precision and kinematics of neutrons and gammas produced. and gammas produced.

• Include K0

Include K0 oszillations

• szillations
• Provide

Provide muon muon nuclear reactions nuclear reactions

• Provide internal conversion

Provide internal conversion

• Provide neutron elastic scattering, in particular recoil energy

Provide neutron elastic scattering, in particular recoil energy and and momentum distributions for neutrons below 10MeV. momentum distributions for neutrons below 10MeV.

• Dito

Dito for n inelastic scattering off Xenon and SiO2, CaCO3, H2O for n inelastic scattering off Xenon and SiO2, CaCO3, H2O

J.P. Wellisch, CERN/EP/CMC

Requirements Cont. Requirements Cont.

• Dito

Dito for capture for capture

• Provide gamma nuclear reactions for gamma

Provide gamma nuclear reactions for gamma energies of less than 100 energies of less than 100 MeV MeV, including cross , including cross-

• sections.

sections.

• Provide radioactive decay after transmutation.

Provide radioactive decay after transmutation.

• Provide k

Provide k-

• shell excitation in radioactive decay

shell excitation in radioactive decay

• Activation of detector material and environment by

Activation of detector material and environment by shower particles shower particles

J.P. Wellisch, CERN/EP/CMC

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

J.P. Wellisch, CERN/EP/CMC

Sample plots Sample plots energy deposition energy deposition

• BTEV: All distributions are

in the expected energy range

J.P. Wellisch, CERN/EP/CMC

Active tasks 2002 Active tasks 2002

• Write ‘educated guess’ physics lists for major use

Write ‘educated guess’ physics lists for major use-

• cases

cases

• Include at least one test

Include at least one test-

• beam simulation into regular validation

beam simulation into regular validation

• Include a complete radiation benchmark into WG level validation

Include a complete radiation benchmark into WG level validation

• Improved validation suite for the cascade energy range

Improved validation suite for the cascade energy range

• Possibly further extension of the high energy validation suite

Possibly further extension of the high energy validation suite

• Plan to contribute to SATIF

Plan to contribute to SATIF-

• 6

6

• Release fully leading particle biased mars

Release fully leading particle biased mars-

• 5 re

5 re-

• write

write

• Release of cascade part of HETC re

Release of cascade part of HETC re-

• write

write

• Improve gamma nuclear reactions in QGS model

Improve gamma nuclear reactions in QGS model

• Make a validation/verification WWW page

Make a validation/verification WWW page

J.P. Wellisch, CERN/EP/CMC

Active tasks 2002, cont. Active tasks 2002, cont.

• Possibly first release of high energy heavy ion reactions in

Possibly first release of high energy heavy ion reactions in QGS model, with option to further extension to QMD QGS model, with option to further extension to QMD

• Revision of the reaction cross

Revision of the reaction cross-

• sections

sections

• Improve the charge state treatment for recoils.

Improve the charge state treatment for recoils.

• Bring kinetic model to releasable state

Bring kinetic model to releasable state

• Bring

Bring inucl inucl cascade code to releasable state cascade code to releasable state

• Research the use of CHIPS in string fragmentation for

Research the use of CHIPS in string fragmentation for intrinsically 3D fragmentation intrinsically 3D fragmentation

• Provide a generic scattering term for cascade type models

Provide a generic scattering term for cascade type models

• Alternative coherent elastic model (

Alternative coherent elastic model (reggee reggee theory based) theory based)

J.P. Wellisch, CERN/EP/CMC

Active tasks 2002, cont. Active tasks 2002, cont.

• Improve electro

Improve electro-

• nuclear cross

nuclear cross-

• sections to take hard

sections to take hard scattering into account. scattering into account.

• Investigate JENDL2.2, and LA150 neutron data libraries

Investigate JENDL2.2, and LA150 neutron data libraries

• Collect (even more) requirements

Collect (even more) requirements

• Release work, coordination

Release work, coordination

• Contribute to maintenance and user support

Contribute to maintenance and user support

• Contribute to architecture working group

Contribute to architecture working group

• Contribute to process improvement/establishment

Contribute to process improvement/establishment

• Contribute to training

Contribute to training

• Publish work (11 papers in the plan…)

Publish work (11 papers in the plan…)

J.P. Wellisch, CERN/EP/CMC

Conclusions Conclusions

J.P. Wellisch, CERN/EP/CMC

Some hadronic physics highlights of Some hadronic physics highlights of late 2001 and 2002 late 2001 and 2002

J.P. Wellisch, CERN/EP/CMC

Swapping to show a few Swapping to show a few transparencies on pre transparencies on pre-

• compound

compound neutron yields. neutron yields.

J.P. Wellisch, CERN/EP/CMC

Preview on kinetic Preview on kinetic model model

160 MeV p on Pb, forward neutrons 585MeV p on Al, forward And backward n and π

J.P. Wellisch, CERN/EP/CMC

Low energy neutrons: Low energy neutrons: G4NDL0.2, 3.7 G4NDL0.2, 3.7

• Are granular selections of data from (alphabetic)

Are granular selections of data from (alphabetic)

• Brond

Brond 2.1 2.1

• CENDL 2.2

CENDL 2.2

• EFF

EFF-

• 3

3

• ENDF/B (VI.0, VI.1, VI.5)

ENDF/B (VI.0, VI.1, VI.5)

• ENSDF

ENSDF

• FENDL/E2.0

FENDL/E2.0

• JEF 2.2

JEF 2.2

• JENDL (3.1, 3.2, FF, 3.3 currently under study)

JENDL (3.1, 3.2, FF, 3.3 currently under study)

• MENDL

MENDL-

• 2(P)

2(P)

• Large parts of the selection is guided by the

Large parts of the selection is guided by the FENDL FENDL-

• 2 selection

2 selection

• G4NDL0.2 for non

G4NDL0.2 for non-

• thermal application

thermal application

J.P. Wellisch, CERN/EP/CMC

The neutron_hp transport models The neutron_hp transport models

• Simulate the cross

Simulate the cross-

• sections and interactions of

sections and interactions of neutrons with kinetic energies below 20 neutrons with kinetic energies below 20 MeV MeV down down to thermal energies . to thermal energies .

• The upper limit is set only by the evaluated data

The upper limit is set only by the evaluated data libraries the code is based on. libraries the code is based on.

• We consider elastic scattering, fission, capture and

We consider elastic scattering, fission, capture and inelastic scattering as separate models inelastic scattering as separate models

• Neutron_hp sampling codes for the ENDF/B

Neutron_hp sampling codes for the ENDF/B-

• VI

VI derived data formats are completely generic derived data formats are completely generic (not including

(not including general R general R-

• matrix for the time being)

matrix for the time being)

• Note that for fission there is a quite competitive theory driven

Note that for fission there is a quite competitive theory driven alternative model, alternative model, G4ParaFissionModel. G4ParaFissionModel.

J.P. Wellisch, CERN/EP/CMC

Models for neutron interaction and Models for neutron interaction and thermalization thermalization. .

J.P. Wellisch, CERN/EP/CMC

Doppler broadening Doppler broadening

J.P. Wellisch, CERN/EP/CMC

The The doppler doppler bias illustrated bias illustrated for Carbon for Carbon

J.P. Wellisch, CERN/EP/CMC

qgs qgs model for model for π π and K induced reactions and K induced reactions

J.P. Wellisch, CERN/EP/CMC

K K-

• , scattering off Au

, scattering off Au

(for (for pions pions see V&V section) see V&V section)

J.P. Wellisch, CERN/EP/CMC

Photon Evaporation data base Photon Evaporation data base

• ≤

≤

• ≤

≤ ≤ ≤

J.P. Wellisch, CERN/EP/CMC

J.P. Wellisch, CERN/EP/CMC

Preliminary test results Preliminary test results

J.P. Wellisch, CERN/EP/CMC

57 57Co

Co

J.P. Wellisch, CERN/EP/CMC

Chiral Chiral Invariant Phase Invariant Phase-

• space Decay(CHIPS)

space Decay(CHIPS)

J.P. Wellisch, CERN/EP/CMC

Vacuum CHIPS Vacuum CHIPS

• 4

2 − N

M

T M

e

/ −

dM e M dW

T M N / 4 2 − −

∝

2 2

) 2 2 ( 2 T N N M − =

J.P. Wellisch, CERN/EP/CMC

Vacuum CHIPS Vacuum CHIPS

• µ

µ

• 3

2 1

−

      − ∝

N

M k kdk dW

∫

        − − − ×         − − =

−

θ θ µ δ µ cos 2 1 ) cos 1 ( 2 2 1 2 1 ) , , (

2 4

qdqd M k kq M k M q M k P

N

J.P. Wellisch, CERN/EP/CMC

Vacuum CHIPS Vacuum CHIPS

• µ:

µ:

• (

)

3 2 2

1 ) 3 ( 4 2

−

− − −

N

kM N k k M µ

J.P. Wellisch, CERN/EP/CMC

Anti proton annihilation Anti proton annihilation

J.P. Wellisch, CERN/EP/CMC

Anti proton annihilation Anti proton annihilation

J.P. Wellisch, CERN/EP/CMC

Nuclear CHIPS Nuclear CHIPS

J.P. Wellisch, CERN/EP/CMC

Nuclear CHIPS Nuclear CHIPS

• µ

µ

• µ

µ

• ∆

∆

• [

]

∫

− + ∆ −         − + ∆ − − =

− kq kq n kq

d k k k k k P θ θ µ µ θ µ µ µ cos ) cos 1 ( ’ 2 ) ( ’ ) cos 1 ( ’ ) ( 2 1 ) , ’ , (

2 3

J.P. Wellisch, CERN/EP/CMC

Nuclear CHIPS Nuclear CHIPS

• ν

ν

• ε1

ε1

• ε2

ε2

• ω

ω

• 1

1

) 1 (

− −

+ =

a a

C P ω ων

ν ν

J.P. Wellisch, CERN/EP/CMC

Nuclear CHIPS Nuclear CHIPS

J.P. Wellisch, CERN/EP/CMC

Intra Intra-

• nuclear CHIPS

nuclear CHIPS

J.P. Wellisch, CERN/EP/CMC

Hard scattering in electro Hard scattering in electro-

• nuclear

nuclear

J.P. Wellisch, CERN/EP/CMC

Hard scattering in electro Hard scattering in electro-

• nuclear

nuclear

J.P. Wellisch, CERN/EP/CMC

A propagation test for QMD development A propagation test for QMD development

J.P. Wellisch, CERN/EP/CMC

The support The support process process – – static view static view

J.P. Wellisch, CERN/EP/CMC

The support process The support process – – dynamic view dynamic view

J.P. Wellisch, CERN/EP/CMC

The hopefully no longer dry numbers… The hopefully no longer dry numbers…

J.P. Wellisch, CERN/EP/CMC

Conclusions Conclusions

J.P. Wellisch, CERN/EP/CMC

GHAD Validation&Verification GHAD Validation&Verification

• Our validation strategy is deployed since spring 1999. It was

Our validation strategy is deployed since spring 1999. It was also submitted as paper to CHEP2001 at the time of the last also submitted as paper to CHEP2001 at the time of the last

• review. My apologies for not having presented it last time.
• review. My apologies for not having presented it last time.
• It was subsequently presented again in CMS and ATLAS, at

It was subsequently presented again in CMS and ATLAS, at the LHC the LHC-

• geant4 validation meeting, and the recent ACAT

geant4 validation meeting, and the recent ACAT conference in Moscow. conference in Moscow.

• I would have been pleased to also present it in invited talks at

I would have been pleased to also present it in invited talks at the SATIF workshop, and the IDM2002 workshop, but I had to the SATIF workshop, and the IDM2002 workshop, but I had to turn these down due to lack of travel turn these down due to lack of travel-

• money.

money.

J.P. Wellisch, CERN/EP/CMC

Model validation Model validation

J.P. Wellisch, CERN/EP/CMC

Author validation Author validation

J.P. Wellisch, CERN/EP/CMC

Working group validation Working group validation

• σ

σ

• σ

σ

• σ

σ

• Ω

Ω

J.P. Wellisch, CERN/EP/CMC

Anti proton annihilation Anti proton annihilation

J.P. Wellisch, CERN/EP/CMC

Anti proton annihilation Anti proton annihilation

J.P. Wellisch, CERN/EP/CMC

Stopping Stopping pion pion minus minus

J.P. Wellisch, CERN/EP/CMC

Low energy neutron capture Low energy neutron capture

J.P. Wellisch, CERN/EP/CMC

Doppler broadening Doppler broadening

J.P. Wellisch, CERN/EP/CMC

Neutron induced isotope production Neutron induced isotope production

J.P. Wellisch, CERN/EP/CMC

Isotope production Isotope production

J.P. Wellisch, CERN/EP/CMC

Proton induced reactions Proton induced reactions

J.P. Wellisch, CERN/EP/CMC

Example WG test Example WG test results results

J.P. Wellisch, CERN/EP/CMC

‘phi’ plots (in ‘phi’ plots (in Pb Pb) )

Nucleon phi distributions For incident π+, π-,p-bar,p At energies 50MeV-40GeV

J.P. Wellisch, CERN/EP/CMC

More ‘phi’ distributions (in lead) More ‘phi’ distributions (in lead)

J.P. Wellisch, CERN/EP/CMC

‘Trivial’ plots ‘Trivial’ plots energy deposition energy deposition

• BTEV: All distributions are

in the expected energy range

J.P. Wellisch, CERN/EP/CMC

Validation in complete applications Validation in complete applications

J.P. Wellisch, CERN/EP/CMC

Benchmark comparisons Benchmark comparisons

J.P. Wellisch, CERN/EP/CMC

Radiation benchmarks Radiation benchmarks – – example example

J.P. Wellisch, CERN/EP/CMC

A sample source spectra A sample source spectra

J.P. Wellisch, CERN/EP/CMC

25 cm 25 cm schielding schielding 43MeV 43MeV

J.P. Wellisch, CERN/EP/CMC

50 cm 50 cm schielding schielding 68MeV 68MeV

J.P. Wellisch, CERN/EP/CMC

50 cm 50 cm schielding schielding 68MeV 68MeV

J.P. Wellisch, CERN/EP/CMC

Test Test-

• beams

beams

J.P. Wellisch, CERN/EP/CMC

Test Test-

• beam sample result

beam sample result

Courtesy of ATLAS TILE prelim. Courtesy of CMS prelim.

J.P. Wellisch, CERN/EP/CMC

A test beams study in regression A test beams study in regression

• ATLAS HEC as a calorimeter benchmark set

ATLAS HEC as a calorimeter benchmark set-

• up

up

• Detailed description of the detector

Detailed description of the detector

• Very constructive help from the ATLAS calorimeter

Very constructive help from the ATLAS calorimeter community community

• A toy analysis: E=E_front + 2E_back

A toy analysis: E=E_front + 2E_back

• Results from the ATLAS test

Results from the ATLAS test-

• beam analysis are

beam analysis are

• verlaid , and labeled as ‘org.’.
• verlaid , and labeled as ‘org.’.
• Data are taken from CALOR 2002 presentation

Data are taken from CALOR 2002 presentation

J.P. Wellisch, CERN/EP/CMC

The physics lists studied in test The physics lists studied in test-

• beam

beam

J.P. Wellisch, CERN/EP/CMC

The overall parameters The overall parameters

J.P. Wellisch, CERN/EP/CMC

ATLAS HEC, CALOR 2002 Relative timing of geant3 and geant4 for pion test-beam

J.P. Wellisch, CERN/EP/CMC

J.P. Wellisch, CERN/EP/CMC

J.P. Wellisch, CERN/EP/CMC

J.P. Wellisch, CERN/EP/CMC

J.P. Wellisch, CERN/EP/CMC

Conclusions Conclusions

J.P. Wellisch, CERN/EP/CMC

Other areas of known usage (likely Other areas of known usage (likely incomplete) incomplete)

J.P. Wellisch, CERN/EP/CMC

Collaboration with 3 Collaboration with 3rd

rd parties

parties Some of the reasoning: Some of the reasoning:

• Geant3 had used two strategies. There were shower

Geant3 had used two strategies. There were shower packages released with geant3, and there were interfaces packages released with geant3, and there were interfaces released with geant3; the latter were interfacing to external released with geant3; the latter were interfacing to external

• packages. The first was a working model, for the latter,
• packages. The first was a working model, for the latter,

geant3 always was claimed to be obsolete. geant3 always was claimed to be obsolete.

• GISMO: the no physics situation, but only interfacing to

GISMO: the no physics situation, but only interfacing to external packages. They never really got support for the use external packages. They never really got support for the use

• f these codes with GISMO.
• f these codes with GISMO.
• MCNPX: Gets it right. They encourage and help 3

MCNPX: Gets it right. They encourage and help 3rd

rd parties to

parties to release MCNP interfaces with their 3 release MCNP interfaces with their 3rd

rd party code. It solves

party code. It solves the support question. the support question.

J.P. Wellisch, CERN/EP/CMC

Collaboration with 3 Collaboration with 3rd

rd parties

parties

• Basis:

Basis: We provide a set of well defined, published, and highly stable

We provide a set of well defined, published, and highly stable interfaces that interfaces that allows interested 3 allows interested 3rd

rd parties to release adapters to use their code, or to use geant4

parties to release adapters to use their code, or to use geant4 physics implementations within their infrastructure. physics implementations within their infrastructure.

• EGS:

EGS: geant4 chips code for

geant4 chips code for γ γ-

• nuclear reactions also in EGS

nuclear reactions also in EGS

• HETC:

HETC: Being re

Being re-

• written to become natively available in G4

written to become natively available in G4

• INUCL:

INUCL: Being integrated to become natively available in G4

Being integrated to become natively available in G4

• UrQMD

UrQMD: : In the process of being re

In the process of being re-

• engineered to become natively

engineered to become natively available in geant4 available in geant4

• MCNP:

MCNP: Discussion on using the geant4 interfaces in MCNP

Discussion on using the geant4 interfaces in MCNP

• G

G-

• FLUKA:

FLUKA: Interfaced by ‘air shower’ users for their own use.

Interfaced by ‘air shower’ users for their own use.

• Liege

Liege Cascade code: Discussion in progress. We hope that they will

Cascade code: Discussion in progress. We hope that they will release a G4 interface soon, and are of course happy to help. release a G4 interface soon, and are of course happy to help.