Re: Special Cuts

Dear Florian,

The calls to gMC_>Gstpar() have to be done in MCApplication::InitGeometry; the calls from MCApplication::ConstructGeometry() only generates a warning and the cuts are not set. (See Ex03MCApplication::InitGeometry().) Could you send me your modifications of the example E03 (or a tar file with your modified code) to reproduce the error?

Thank you,

Ivana

On 02/09/2011 09:21 AM, Florian Zenoni wrote:
> Dear Ivana,
>
> I'm using physics list LHEP_EMV or CHIPS. I tried to reproduce the
> same code from the example E03, then I set the cuts in my
> DetectorConstruction.cxx at the end of
> DetectorConstruction::ConstructGeometry(), right after
> gMC->SetRootGeometry().
> Eventually, I could reproduce the error in the example E03, but until
> this time only with the physics list LHEP. However in the table that
> resumes the cuts (at the beginning of the run), they seem to be
> correct.
> On the contrary LHEP_EMV and CHIPS, that give me the error in my
> simulation, run well in the example. Of course, in the same table at
> the beginning of my simulation, the cuts are ignored.
>
> Thank you,
> Florian
>
> 2011/2/8 Ivana Hrivnacova<Ivana.Hrivnacova_at_cern.ch>:
>>
>> Dear Florian,
>>
>> Could you, please, give more information?
>> Which physics list do you use? Can you reproduce the error by modifying the
>> E03 example?
>>
>> Thanks,
>>
>> Ivana
>>
>> On 02/07/2011 04:34 PM, Florian Zenoni wrote:
>>>
>>> Hello everybody,
>>>
>>> I am using Root 5.28, Geant4VMC 2.11 and Geant4.9.4.
>>> I'm working on a silicon detector simulation, and I would like to
>>> apply some cuts in Geant3 way, as explained here:
>>> http://root.cern.ch/drupal/content/special-cuts-and-regions.
>>> So, in my g4tgeoConfig.C, I've added "specialCuts" as the third
>>> argument of TG4RunConfiguration. Then I've set something like this:
>>> gMC->Gstpar(mediumId, "CUTGAM", 10.e-06), as found in the example E03.
>>>
>>> Here is the error when I run my simulation:
>>> "G4ProcessTable::Insert : arguments are 0 pointer".
>>>
>>> Has someone any idea?
>>>
>>> Thank you very much.
>>> Florian
>>>
>>
>> --
>> ==============================================================
>> e-mail: Ivana.Hrivnacova_at_cern.ch
>> address: Institut de Physique Nucléaire, 91406 Orsay, France
>> phone: +33 169156594
>> ==============================================================
>>

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