GSI Forum: General » Barrel DIRC in Geant 4 - Cherenkov photon propagation

Home » PANDA » PandaRoot » General » Barrel DIRC in Geant 4 - Cherenkov photon propagation - [FIXED]

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

Barrel DIRC in Geant 4 - Cherenkov photon propagation - [FIXED] [message #10988]

Tue, 14 September 2010 18:38

Jochen Schwiening
Messages: 84
Registered: June 2009
Location: GSI, Darmstadt

continuous participant

From: *gsi.de

Hi folks,
in an effort to study our frustrating problems with the photon gun usage in G3 (see Maria's post titled "Cherenkov photons generation") I tried to run our usual macro /macro/drc/sim_dirc.C using Geant 4. I enabled the optical process ("+optical" in the g4 config file) and ran a couple of events, looked at them with the event display.
I see that Cherenkov photons are produced but I also see that they are all lost/killed/absorbed as soon as they hit the first bar surface. I am attaching a screen snapshot showing a side view, a muon track producing many photons that all die as soon as the first surface is hit. It looks as if G4 doesn't know what to do with the photons, how to propagate them inside the material.
I should mention that I see an enormous number of warning messages
"Warning: G4MaterialPropertyVector::GetProperty==> attempt to Retrieve Property above range"
which may very well indicate a fatal problem - sadly, I don't know where to look to fix the issue.
It's the first time I tried G4 in PandaROOT and so far I only heard that the barrel DIRC "doesn't work" in G4. It'd be nice to change that.
Any ideas?
Thanks, Jochen

Attachment: dirc_geant4_prop.png
(Size: 9.67KB, Downloaded 308 times)

[Updated on: Wed, 15 September 2010 15:09]

Report message to a moderator

Re: Barrel DIRC in Geant 4 - Cherenkov photon propagation [message #10990 is a reply to message #10988]

Wed, 15 September 2010 13:49

Jochen Schwiening
Messages: 84
Registered: June 2009
Location: GSI, Darmstadt

continuous participant

From: *gsi.de

Just a quick update: the photon propagation issues seems to be solved.
The photons were actually stopped by our own DIRC code PndDrc.cxx where we stop photons that leave the bar volume in a location other than the bar end. This was done to avoid tracking Cherenkov photons that left the bar because they were not internally reflected but then keep getting scattered around PANDA and use up tons of CPU time.
This cut worked well in Geant 3 but in Geant 4 there seems to be a volume change step even during internal reflections (seems wrong but at least that what it looks like), which caused all photons to be stopped by our PndDrc code.
After I removed this explicit stop command I get photons propagating down the bar and registering on the detector pixels, as can be seen in the attached snapshot. Looks fine to me.
(The photons seen leaving the bar close to where the track hits the DIRC - those are the ones we wanted to stop. Keeping them increases the time per event by a factor of 5...)

By adding additional entries for short and long wavelengths into the media_pnd.geo for DIRCairNoSens I was able to make most of the annoying warning message go away.

I'll follow up on the photon gun issue in G4 in Maria's thread.

Cheers, Jochen

Attachment: dirc_geant4_prop_ok.png
(Size: 24.37KB, Downloaded 211 times)

Report message to a moderator

Re: Barrel DIRC in Geant 4 - Cherenkov photon propagation [message #10991 is a reply to message #10990]

Wed, 15 September 2010 14:25

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *to.infn.it

Hi,
yes you are right. Ages ago with Annalisa we have introduced that cut to speed up the simulation, but I had completely forgotten about its existence.

Nice to know that there is one problem less Smile

Could you please add "[FIXED]" to the topic name of your first message? Thanks.

Report message to a moderator

Re: Barrel DIRC in Geant 4 - Cherenkov photon propagation [message #10992 is a reply to message #10990]

Wed, 15 September 2010 14:27

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

first-grade participant

From: *dip.t-dialin.net

Hallo Jochen,

I got a mail from Peter Koch in Giessen (I forward it at that time to Dipanwita and Carsten) about a bug in the reflection of Cherenkov photons in G4.9.3 (http://bugzilla-geant4.kek.jp/show_bug.cgi?id=1103) this is corrected in the beta of G4.9.4 but it still not released. may be this is related to what you see with G4.

The other topic I am not sure about, I will answer Maria's mail and post there a macro because with the box generator I can create Cherenkov photon without problems!.

regards

Mohammad

Report message to a moderator

Re: Barrel DIRC in Geant 4 - Cherenkov photon propagation [message #11128 is a reply to message #10992]

Wed, 27 October 2010 13:33

Maria Patsyuk
Messages: 58
Registered: April 2010

continuous participant

From: *gsi.de

As Jochen wrote in the topic "Barrel DIRC in GEANT4 - Cherenkov photons propagation [FIXED]" we are now able to generate cherenkov photons with box generator in GEANT4, but there are a lot of the following messages while the simulation:

*** G4Exception : ZeroPolarization
issued by : G4PrimaryTransformer::GenerateSingleTrack
Polarization of the optical photon is null. Random polarization is assumed.
*** This is just a warning message.
This warning message is issued up to 10 times.

Is it a harmless warning or should I do something about it?

Report message to a moderator

Re: Barrel DIRC in Geant 4 - Cherenkov photon propagation [message #11135 is a reply to message #11128]

Wed, 27 October 2010 18:06

StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino

first-grade participant

From: *vpn.unito.it

No idea,
but I fear that you have to wait for the next release of external packages (once new g4 v 9.4 will be ready) to have proper photon propagation in geant4...

Report message to a moderator

Previous Topic:	PndMicroWriter???
Next Topic:	Problems with PandaRootCharmoniumTutorial

Goto Forum:

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Apr 24 07:22:39 CEST 2024

Total time taken to generate the page: 0.00568 seconds