Re: Special Cuts

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
> ==============================================================
>
Received on Fri Feb 11 2011 - 09:00:51 CET