GSI Forum: EMC » Geant 3, Geant 4 VMC, native geant 4 comparison

Home » PANDA » PandaRoot » EMC » Geant 3, Geant 4 VMC, native geant 4 comparison

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

Geant 3, Geant 4 VMC, native geant 4 comparison [message #5940]

Fri, 22 February 2008 12:20

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

first-grade participant

From: *fuw.edu.pl

Hi again,

I put here a couple of plots that demonstrate difference in deposited energy for Geant3, Geant4 VMC and native Geant4.

So cluster_energy_full.jpg - index.php?t=getfile&id=4193&private=0

is energy deposited from 1GeV photon in emc barrel (theta=50, phi=0-360).

And it is comparable with energy deposited in EmcTest setup (5x5 PbWO4 crystals), implemented by Stefano.
index.php?t=getfile&id=4194&private=0

So around 5% of energy is lost somewhere for Geant4 and it doesn't come from complicated emc geometry.

But with native Geant4 application
index.php?t=getfile&id=4195&private=0

deposited energy is similar to Geant3 VMC and close to the value calculated by Jan in Babar framework for the whole EMC.
http://tau.ep1.rub.de/HyperNews/get/SimGenDisc/2159/1/1/1/1/1.html

So the difference doesn't come from the different Geant4 version
4.7 in Babar framework vs. 4.9 in pandaroot (which I also used in native Geant4 application).

So the source of difference is still under investigation,

Dima

Attachment: cluster_energy_full.jpg
(Size: 35.82KB, Downloaded 897 times)
Attachment: cluster_energy_test.jpg
(Size: 34.97KB, Downloaded 930 times)
Attachment: cluster_energy_g4.jpg
(Size: 32.42KB, Downloaded 940 times)

Report message to a moderator

Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5942 is a reply to message #5940]

Fri, 22 February 2008 13:02

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *physik.uni-giessen.de

Which physics list are you using for native Geant4?

Report message to a moderator

Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5944 is a reply to message #5942]

Fri, 22 February 2008 13:30

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

first-grade participant

From: *fuw.edu.pl

I used standard Electromagnetic Physics List only - PhysListEmStandard.cc which I copied from one of Geant4 examples (/geant4/examples/extended/electromagnetic/TestEm14/src/), but it is identical in most of electromagnetic examples. In Geant4 there are several alternative electromagnetic physics lists, but this is recommended.

And I used default cut value - 1mm, which is actually used in Babar framework.

In pandaroot in g4config.C there is a cut not in range but in energy - 1MeV, but when I run simulation in pandaroot whith Geant4 it prints energy cuts which correspond to 1 mm range cut for PWO (84.8 keV for gamma and 1.13 MeV for electrons) but if I put these values instead of 1 MeV the result is the same.

Dima

Report message to a moderator

Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5945 is a reply to message #5942]

Fri, 22 February 2008 13:38

Johan Messchendorp
Messages: 693
Registered: April 2007
Location: University of Groningen

first-grade participant

From: *KVI.nl

Hi Dima and others,

Thanks for looking into this problem. I think we can learn a lot from it.
I had another question. I was wondering whether it would make sense to run the Geant4 mini-EMC detector within the Geant4-VMC templates provided by Ivana in cbmsoft/transport/geant4_vmc/examples/.... This would be a simulation with the "native" Geant4-VMC, which might differ from - lets say - the CBM/PandaRoot-VMC.

Kind wishes,

Johan.

Johan Messchendorp
University of Groningen/KVI
Zernikelaan 25
NL-9747 AA Groningen
The Netherlands
tel. +31-503633558
fax +31-503634003

Report message to a moderator

Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5946 is a reply to message #5945]

Fri, 22 February 2008 13:42

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

first-grade participant

From: *fuw.edu.pl

Hi Johan,

In principle I planned it ("native" Geant4-VMC) as the next step.

Dima

Report message to a moderator

Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5947 is a reply to message #5946]

Fri, 22 February 2008 13:45

Johan Messchendorp
Messages: 693
Registered: April 2007
Location: University of Groningen

first-grade participant

From: *KVI.nl

Great Dima,

Let me know the outcome....

Johan.

Johan Messchendorp
University of Groningen/KVI
Zernikelaan 25
NL-9747 AA Groningen
The Netherlands
tel. +31-503633558
fax +31-503634003

Report message to a moderator

Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5964 is a reply to message #5947]

Mon, 25 February 2008 21:26

Mohammad Al-Turany
Messages: 518
Registered: April 2004
Location: GSI, Germany

first-grade participant

From: *dip.t-dialin.net

Hi all,

I tried to reproduce this problem, so I switched on only the emStandard list for G4 in the VMC which correspond to the PhysListEmStandard.cc used in native G4, i.e:

Config()
{
....
TG4RunConfiguration* runConfiguration
= new TG4RunConfiguration("geomRoot", "emStandard");
...
}

The result was as you can see in the attached figures some how acceptable!! So it seems that all the problem comes from the combination of the physics lists in VMC. So we need to investigate further what is going on here!

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

and drawing both over each other (G4 in red )

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

Attachment: g3_g4.gif
(Size: 15.63KB, Downloaded 701 times)
Attachment: g3_g4_overlap.gif
(Size: 10.58KB, Downloaded 791 times)

Report message to a moderator

Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5966 is a reply to message #5964]

Mon, 25 February 2008 21:43

Johan Messchendorp
Messages: 693
Registered: April 2007
Location: University of Groningen

first-grade participant

From: *xs4all.nl

Hi Mohammad,

Just for completeness. Did you do a simulation with the 5x5 matrix? The reason why I am asking, is that I noticed that the response for G3 from Dima differ from the G3 from your simulation (albeit slightly, but nevertheless significant by looking at the average energy)....

Johan.

Johan Messchendorp
University of Groningen/KVI
Zernikelaan 25
NL-9747 AA Groningen
The Netherlands
tel. +31-503633558
fax +31-503634003

Report message to a moderator

Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5967 is a reply to message #5966]

Mon, 25 February 2008 22:43

Mohammad Al-Turany
Messages: 518
Registered: April 2004
Location: GSI, Germany

first-grade participant

From: *dip.t-dialin.net

Hallo Johann,
and here is the 5x5 matrix results:

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

Attachment: 5x5g3g4.gif
(Size: 15.61KB, Downloaded 790 times)

Report message to a moderator

Previous Topic:	pwo25 test application in emc directory
Next Topic:	Code for APD

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed May 31 12:10:47 CEST 2023

Total time taken to generate the page: 0.00872 seconds