GSI Forum: Monte Carlo Engines » dipole field in Geant4

Home » PANDA » PandaRoot » Monte Carlo Engines » dipole field in Geant4

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

dipole field in Geant4 [message #12233]

Fri, 08 July 2011 12:22

Mathias Michel
Messages: 21
Registered: May 2009
Location: Mainz

occasional visitor

From: *kph.uni-mainz.de

Hi all,

it seems that somehow, if i switch in a simulation from Geant3 to Geant4, its ignoring the dipole field. Attached you can find two event displays created with the same macro, once with Geant3 and once with Geant4. You can see that with Geant4 the "kink" of the dipole field is missing. I attached the macro as well so you can have a look at the parameters.

I tried both external package release's from this year and pandaroot revision 12630 was the latest i tried. My machine is a FedoraCore15 32bit.

Is there something wrong with the macro, a.e. the way the fieldmaps are read, or does anybody else see see a similar behaviour?

Cheers,
Mathias

Attachment: Geant3.jpg
(Size: 17.00KB, Downloaded 470 times)
Attachment: Geant4.jpg
(Size: 16.76KB, Downloaded 449 times)
Attachment: runMvdSim.C
(Size: 5.38KB, Downloaded 400 times)

Report message to a moderator

Re: dipole field in Geant4 [message #12241 is a reply to message #12233]

Tue, 12 July 2011 10:31

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

Hi,
I think I was able to reproduce your problem.. or at least it seems that I have signal in the forward emc with g3, and no signal with g4.

First I would like to know which external packages and trunk versions you are using.
Second, I would like to understand if this is the first time you are trying with g4, or if maybe before it was working and now it does not (trying to find which is the change that is creating the mess, or of this problem exists since some time).
Third, could you please put the macro you are using as event display?
Fourth... I think we have to wait that Mohammad comes back from vacation... Sad

Report message to a moderator

Re: dipole field in Geant4 [message #12243 is a reply to message #12241]

Tue, 12 July 2011 10:55

Mathias Michel
Messages: 21
Registered: May 2009
Location: Mainz

occasional visitor

From: *kph.uni-mainz.de

Hi,

thanks for having a look into this.

I'm using panda-trunk 12630 with may11 external package release. I also tried the feb11 release with same results, but no other up to now.

I used the Mvd-macro because i thought its easier to check. But for the lmd-macros i can say that they worked with both Geant3 and Geant4 up to a few months ago at least. Now in my setup the lmd-macros ignore the dipole-field with both G3 and G4. With Mvd it works with G3 but not with G4 as shown in my first post.
Maybe this is a problem with my local machine since i upgraded to fedora15 recently, but its so strange that i'd like to find out whats happening and found it worth mentioning.

You can find the eventdisplay attached.

Until Mohammad is back i will run the same macros in the same setup on a different machine to check if this is just a problem of my local machine.

Attachment: eventDisplay.C
(Size: 3.58KB, Downloaded 224 times)

Report message to a moderator

Re: dipole field in Geant4 [message #12244 is a reply to message #12243]

Tue, 12 July 2011 14:07

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

first-grade participant

From: *185.x.209.go.com.jo

Hi,

There was a problem in the initialization of the field when using G4. now it is corrected in r12662, please check it.

regards,

Mohammad

Report message to a moderator

Re: dipole field in Geant4 [message #12245 is a reply to message #12244]

Tue, 12 July 2011 14:40

Mathias Michel
Messages: 21
Registered: May 2009
Location: Mainz

occasional visitor

From: *kph.uni-mainz.de

Hi,

with G4 it is now working again (lmd & mvd), thanks!

Unfortunately i now have some trouble with G3. In case of lmd the antiprotons are still unbended by the dipole. And with the mvd-macro from above the eventdisplay now is empty and i get the following message for each track during simulation:

*** GTRACK *** More than 10000 steps, tracking abandoned!

Report message to a moderator

Re: dipole field in Geant4 [message #12255 is a reply to message #12245]

Wed, 13 July 2011 19:26

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

first-grade participant

From: 46.248.220*

Hi,

I will check this when I get back to the GSI.

Mohammad

Report message to a moderator

Previous Topic:	[Solved]Problem with PndEvtGenDirect
Next Topic:	Random seed in Geant4?

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Apr 20 02:36:08 CEST 2024

Total time taken to generate the page: 0.00963 seconds