GSI Forum: Monte Carlo Engines » Difference of neutral candidate between Geant4 and Geant3

Home » PANDA » PandaRoot » Monte Carlo Engines » Difference of neutral candidate between Geant4 and Geant3

Show: Today's Messages :: Polls :: Message Navigator

Difference of neutral candidate between Geant4 and Geant3 [message #15624]

Thu, 31 October 2013 10:12

donghee
Messages: 385
Registered: January 2009
Location: Germnay

first-grade participant

From: *kph.uni-mainz.de

Dear all,

I am trying to use Geant4 in the simulation, and have seen quite different behaviour for Geant4 and Geant3 at the multiplicity of neutral candidate.

Quote:

PANDAroot apr13 released version.
Beam momentum = 15 GeV
MC : DPM
Cut on the neutral candidate : E > 10 MeV

I compared two cases, one is the charged candidate and another one is the neutral candidate. please find the plot in attachment.

For charged candidate, both Geant3 and Geant4 look pretty much same, but for neutral case, the size of neutral candidates significantly different.
I didn't touch any option in the g4config.C

Which engine describe better for photon simulation?
Do I need some fine tuning from option control to use Geant4?

Best wishes,
Donghee

Attachment: multiplicity_dpm_dc_geant.eps
(Size: 14.38KB, Downloaded 462 times)

Report message to a moderator

Re: Difference of neutral candidate between Geant4 and Geant3 [message #15633 is a reply to message #15624]

Fri, 01 November 2013 12:41

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: 2.235.190*

Can you please switch off the specialcuts and specialcontrols in your g4config.C, and check again the g4 outcome?

Report message to a moderator

Re: Difference of neutral candidate between Geant4 and Geant3 [message #15706 is a reply to message #15633]

Mon, 11 November 2013 13:46

donghee
Messages: 385
Registered: January 2009
Location: Germnay

first-grade participant

From: *kph.uni-mainz.de

Hi Stefano,

Hier is the plot for the multiplicity from single photon of box generator using geant3 and geant4.
There are no difference between them.
Energy distributions are perfectly identical for both.

And then in second page I have tried with 6 different particle sources
with gamma, electron, muon, pion, kaon and proton.
They are produced with box generator, each set have been generated with multiplicity 6.
Every hadron sources make the difference between Geant3 and Geant4, while lepton or photon source show same distribution.

Best wishes,
Donghee

Attachment: Pandaroot_meeting_2.pdf
(Size: 97.69KB, Downloaded 344 times)

Report message to a moderator

Re: Difference of neutral candidate between Geant4 and Geant3 [message #15708 is a reply to message #15706]

Mon, 11 November 2013 18:04

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

Could you please try to remove the special cuts/controls in your gconfig/g4Config.C and try again?

TG4RunConfiguration* runConfiguration
        = new TG4RunConfiguration("geomRoot", "QGSP_BERT_EMV", "");

Report message to a moderator

Re: Difference of neutral candidate between Geant4 and Geant3 [message #15710 is a reply to message #15708]

Tue, 12 November 2013 10:53

donghee
Messages: 385
Registered: January 2009
Location: Germnay

first-grade participant

From: *scc.kit.edu

Hi Stefano,

I tried your suggestion.

Quote:

TG4RunConfiguration* runConfiguration = new TG4RunConfiguration("geomRoot", "QGSP_BERT_EMV", "");

Single pion had been simulated (see attached file).
Same behaviour showed up in the neutral multiplicity(candidate).
6 pion simulation is now running. but assume that will be same.

I think that Geant4 and Geant3 define the cluster differently, specially for the geometrical size of cluster.
Or Geant3 produce more photon from the interaction in the front of EMC at DIRC and DISC.

I hope that this discrepancy is not so big trouble in the analysis level.
I need more time to investigate further.

Ciao,
Donghee

Attachment: multiplicity_for_single_pion_without_option.eps
(Size: 9.37KB, Downloaded 315 times)

Report message to a moderator

Re: Difference of neutral candidate between Geant4 and Geant3 [message #15724 is a reply to message #15710]

Fri, 15 November 2013 15:34

donghee
Messages: 385
Registered: January 2009
Location: Germnay

first-grade participant

From: *dip0.t-ipconnect.de

Hi Stefano,

I have tried to look more detail the dependency of multiplicity with EMC cluster for Geant3 and Geant4 including physics option or excluding options.
Please have a look all distributions in attached file.

Main difference at Multiplicity between Geant3 and Geant4 start already from the EMC cluster size.

With option at Geant4 many of photon production are restricted. but no option produce huge amount of photon mainly due to secondary interaction.
That make main difference at the size of EMC cluster.

I think that this issue is surely related with EMC simulation.
Is someone already reported about Geant4 simulation for EMC part? or will be tested in the future?

Best wishes,
Donghee

Attachment: Pandaroot_meeting_2.pdf
(Size: 321.57KB, Downloaded 324 times)

Report message to a moderator

Re: Difference of neutral candidate between Geant4 and Geant3 [message #15727 is a reply to message #15724]

Fri, 15 November 2013 16:15

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

Hi,
the last try would be to use different physics lists in g4, and see if there are changes. I am not able to understand if the problems are in the vmc side or in our side. As far as I know g4 response was already tested for emc, but is seems there is something wrong somewhere.

Report message to a moderator

Re: Difference of neutral candidate between Geant4 and Geant3 [message #15835 is a reply to message #15624]

Thu, 06 February 2014 11:56

Dima Melnychuk
Messages: 213
Registered: April 2004
Location: National Centre for Nucle...

first-grade participant

From: *fuw.edu.pl

Hi,

As Stefano suggested I tried to compare the number of neutral candidate between Geant 3 and Geant 4 with different "hadronic process" setting for Geant 3. The default value for "hadronic process" in /gconfig/SetCuts.C is 1.
gMC->SetProcess("HADR",1);

I tried "3" and "5" and with "5" simulation stuck. But actually between "hadronic process" "1" and "3" with respect to number of neutral candidates there is no much difference.

So the results for 15 GeV DPM events.
index.php?t=getfile&id=7679&private=0

Multiplicity for Geant 4 is much lower, and almost no difference between different "hadronic process" in Geant3.

Distribution of cluster size (number of crystals per cluster)
index.php?t=getfile&id=7680&private=0

indicate 2 times difference in single cluster crystals between Geant3 and Geant4.

And if we consider the multiplicity of neutrals with more than 2 crystals:
index.php?t=getfile&id=7681&private=0

there is already a reasonable agreement between Geant3 and Geant4.

And still there is no conclusion which engine is more correct Geant3 or Geant4.

Dima

Attachment: number_neutral.png
(Size: 16.77KB, Downloaded 738 times)
Attachment: number_crystals.png
(Size: 14.48KB, Downloaded 676 times)
Attachment: number_neutral_gt2.png
(Size: 14.34KB, Downloaded 765 times)

Report message to a moderator

Re: Difference of neutral candidate between Geant4 and Geant3 [message #15844 is a reply to message #15835]

Fri, 07 February 2014 19:46

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

Thanks Dima for your check.
I also tried with gcalor (5) and with DPM it stucks at the second event, but with evtgen it goes up to the end.
Maybe there are some problems in gcalr with some kind of particle, maybe.

Report message to a moderator

Re: Difference of neutral candidate between Geant4 and Geant3 [message #17721 is a reply to message #15844]

Fri, 09 January 2015 22:21

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: 93.49.202*

Hi,
I did some recent changed triggered by the fact that currently geant4 simulation is faster than geant3.

I simulated the standard psi(2S)->J/psipi+pi- events and compared the emc digi multiplicity with geant3, geant4, and geant4 setting off the "special cuts" and "special controls", i.e. the emulation of geant3 cuts in energy.
This is what I get:

index.php?t=getfile&id=8230&private=0

I understand that the lower multiplicity of geant4 comes not from geant itself but from the special cuts in energy that we use by default. If we remove the flag from gconfig/g4Config.C the computing time increases a lot, become a factor 2 slower than g3, and the multiplicity is higher than geant3 as expected.

At the end, I understand this is just a fact of cuts, and only a comparison with experimental data can help to find the better thresholds.

Attachment: emcdigi_mult.gif
(Size: 11.74KB, Downloaded 802 times)

Report message to a moderator

Previous Topic:	MC true match with geant4
Next Topic:	Pythia as external decayer for Geant3

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Mon Dec 02 05:33:07 CET 2024

Total time taken to generate the page: 0.00829 seconds