Re: gMC->Gstpar crash

Hi Filimon,

I cannot reproduce the crash when I put your three lines in the E03 example, in Ex03DetectorConstruction::SetCuts(). It runs ok with Geant3 and there are issued warnings when running with Geant4, saying that this parameter is not implemented.
Did you put these calls in your MCApplication::InitGeometry()?

Best regards,

Ivana

On 04/20/2011 08:50 PM, Filimon Roukoutakis wrote:
> Hi, I am trying to enable Birk's attenuation under G3 through VMC. The
> simplest way as far as I know would be to do
>
> gMC->Gstpar(3, "BIRK1", 1);
> gMC->Gstpar(3, "BIRK2", 0.013);
> gMC->Gstpar(3, "BIRK3", 9.6e-6);
>
> I suppose then that GBIRK is going to be called appropriately in the
> step manager and the value gMC->Edep() is going to be the attenuated one
> in principle right? (please verify) I would happily check this if the
> crash below would not exist. The simulation seems to work ok if these
> statements are not involved. Any ideas? Medium 3 is the expected one as
> far as I can tell (but I am open to some suggestion to really verify
> this also if needed). Thanks,
> filimon
>
> ===========================================================
> #10 0xb5f03098 in g3stpar_ (itmed=0xbff12824, chpar=0xa43917c,
> parval=0xbff1280c, __g77_length_chpar=5) at gcons/gstpar.F:63
> #11 0xb60db5a2 in TGeant3::Gstpar (this=0xa36ee88, itmed=3,
> param=0xa43917c "BIRK2", parval=0.012999999999999999) at
> TGeant3/TGeant3.cxx:3641
> ===========================================================
>

-- 
==============================================================
e-mail:  Ivana.Hrivnacova_at_cern.ch
address: Institut de Physique Nucléaire, 91406 Orsay, France
phone:   +33 169156594
==============================================================