Skip to Content.
Sympa Menu

sphenix-software-l - [Sphenix-software-l] GEANT4 and hadronic showers

sphenix-software-l AT lists.bnl.gov

Subject: sPHENIX discussion of software

List archive

Chronological Thread  
  • From: John Lajoie <lajoie AT iastate.edu>
  • To: "sphenix-hcal-l AT lists.bnl.gov" <sphenix-hcal-l AT lists.bnl.gov>, <sphenix-software-l AT lists.bnl.gov>
  • Subject: [Sphenix-software-l] GEANT4 and hadronic showers
  • Date: Wed, 21 Oct 2015 11:33:05 -0500

Dear sPHENIX Folks:

    After the spirited discussions in both the simulations and HCAL meeting yesterday regarding the energy resolution we are seeing in the MC, John Haggerty pointed me to a very nice paper from the CALICE collaboration:

http://arxiv.org/abs/1510.04063

I have been pretty vocal about voicing skepticism about the low-energy neutron response in GEANT4, but they show time structure data that seem to show the late time (neutron) hits are well-reproduced with the QGSP_BERT_HP physics list in a tungsten calorimeter.  They also claim (but don't show) that "for steel the high precision modelling has much less influence", which is consistent with what we have seen when comparing HCAL simulations with QGSP_BERT and QGSP_BERT_HP physics lists (QGSP_BERT is our default).  OK, I'll shut up about the neutrons.

Something else that was discussed in the HCAL meeting in particular was whether or not what I am showing is inconsistent with what has been shown in the proposal. It certainly is (see the second of the attached slides).  The low energy response is significantly better in the new simulations, which is reflected in the improved stochastic term when fitting the full distribution. However, the fit over all momentum does have a very poor Chi2.  If I fit only to points >16GeV (last slide) the Chi2 is much better and the projection to lower energies overshoots the data.

Edward suggested that some of this might be due to the fact that the EMCal is quite different in the new simulations (current simulations use the 2D spacal) and low-energy hadronic showers deposit a fair fraction of their energy in the EMCal.  The next-to-last slide shows the EMCal energy fraction distributions as a function of energy for pi- at eta=0 and a selection of energies.  There is certainly a transition below 10GeV.

Over the weekend I'll try to get resolution results with towering and clustering for comparison.

Regards,
John

--

John Lajoie

PHENIX Deputy Spokesperson

Professor of Physics

Iowa State University

 

(515) 294-6952

lajoie AT iastate.edu


Contact me: john.lajoie

Attachment: compare2.pdf
Description: Adobe PDF document




Archive powered by MHonArc 2.6.24.

Top of Page