GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » EMC » Geant 3, Geant 4 VMC, native geant 4 comparison
Geant 3, Geant 4 VMC, native geant 4 comparison [message #5940] Fri, 22 February 2008 12:20 Go to next message
Dima Melnychuk is currently offline  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
Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5942 is a reply to message #5940] Fri, 22 February 2008 13:02 Go to previous messageGo to next message
StefanoSpataro is currently offline  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?
Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5944 is a reply to message #5942] Fri, 22 February 2008 13:30 Go to previous messageGo to next message
Dima Melnychuk is currently offline  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
Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5945 is a reply to message #5942] Fri, 22 February 2008 13:38 Go to previous messageGo to next message
Johan Messchendorp is currently offline  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
Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5946 is a reply to message #5945] Fri, 22 February 2008 13:42 Go to previous messageGo to next message
Dima Melnychuk is currently offline  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
Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5947 is a reply to message #5946] Fri, 22 February 2008 13:45 Go to previous messageGo to next message
Johan Messchendorp is currently offline  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
Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5964 is a reply to message #5947] Mon, 25 February 2008 21:26 Go to previous messageGo to next message
Mohammad Al-Turany is currently offline  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

Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5966 is a reply to message #5964] Mon, 25 February 2008 21:43 Go to previous messageGo to next message
Johan Messchendorp is currently offline  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
Re: Geant 3, Geant 4 VMC, native geant 4 comparison [message #5967 is a reply to message #5966] Mon, 25 February 2008 22:43 Go to previous message
Mohammad Al-Turany is currently offline  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 859 times)
Previous Topic: pwo25 test application in emc directory
Next Topic: Code for APD
Goto Forum:
  


Current Time: Thu Mar 28 16:15:08 CET 2024

Total time taken to generate the page: 0.00828 seconds