Re: Special Cuts

Dear Florian,

I was able to be able to reproduce the problem in the Geant4 native example and I filled in a bug report (#1173).

The messages do not seem to be critical: they are invoked when passing an inactivated physics process which is represented by the zero value in the particle process table. When using specialCuts, there is a special process attached to all particles and then it is activated only for particles for which it is relevant according to the cuts set by user. That's why you observe the messages only with SpecialCuts activated.

I will let you know when I have a response from the Geant4 bug system.

Best regards,

Ivana

On 02/11/2011 09:00 AM, Florian Zenoni wrote:
> Dear Ivana,
>
> the calls to gMC->Gstpar() were indeed done in the wrong place. Now
> they are in MCApplication::InitGeometry, but the problem still
> persists.
> I send you in attachment (in a separate e-mail) the tar file of the
> example E03 modified.
> Note: in the tar file I have set the physics list as LHEP to obtain
> the error, but in my simulation it does NOT provoke me any error
> (unlike LHEP_EMV, CHIPS, QGSP, etc.).
>
> Thank you,
> Florian
>
> 2011/2/9 Ivana Hrivnacova<Ivana.Hrivnacova_at_cern.ch>:
>> 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
>> ==============================================================
>>

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