[Gate-users] Gate 9.0 does not fill Hits Tree
Jan Gajewski
jan.gajewski at ifj.edu.pl
Fri Sep 4 09:07:37 CEST 2020
Hi David
>> was it the same with previous Gate version? The way the trees are filled
changed from 9.0
no. in gate 8.2 Hits were filled and also the Phase Space Actor.
Jan: is it similar with other output type (non ascii, not hits but singles
etc) ?
The Hits Tree is not filled in any format (ASCII/root/numpy) in Gate 9.0.
Also Singles is empty. The output is also empty if I define it as:
/gate/output/root/enable
/gate/output/root/setFileName tmp/out
/gate/output/root/setRootHitFlag 1
/gate/output/root/setRootSinglesFlag 1
or
/gate/output/tree/enable
/gate/output/tree/addFileName Results/avTPXHits.npy
/gate/output/tree/addFileName Results/avTPXHits.root
/gate/output/tree/hits/enable
best
---
Jan Gajewski
*jan.gajewski at ifj.edu.pl <jan.gajewski at ifj.edu.pl>**Institute of Nuclear
Physics* (IFJ PAN)
Radzikowskiego 152, 31-342 Kraków, Poland
phone: (+48) 12 662 8331
On Fri, 4 Sep 2020 at 07:48, David Sarrut <David.Sarrut at creatis.insa-lyon.fr>
wrote:
> Hello,
>
> anyone to have a look at this issue ?
> https://github.com/OpenGATE/Gate/issues/335
>
> Maybe Mateusz or Ane ?
>
> It may be related to SensitiveDetector. It seems that when using the
> Actor, the SD of the HitsTree is not enabled.
>
> Jan : was it the same with previous Gate version? The way the trees are
> filled changed from 9.0
> Jan: is it similar with other output type (non ascii, not hits but singles
> etc) ?
>
> thanks,
> David
>
> ---------- Forwarded message ---------
> From: Jan Gajewski <jan.gajewski at ifj.edu.pl>
> Date: Fri, Sep 4, 2020 at 7:17 AM
> Subject: Re: [Gate-users] Gate 9.0 does not fill Hits Tree
> To: <gate-users at lists.opengatecollaboration.org>
>
>
> Thanks David for the response. I saw the issue on GitHub. It is actually
> the same as https://github.com/OpenGATE/Gate/issues/335. I hope the issue
> will be solved soon.
>
> best
>
>
> ---
> Jan Gajewski
>
> *jan.gajewski at ifj.edu.pl <jan.gajewski at ifj.edu.pl>**Institute of Nuclear
> Physics* (IFJ PAN)
> Radzikowskiego 152, 31-342 Kraków, Poland
> phone: (+48) 12 662 8331
>
>
>
> On Thu, 3 Sep 2020 at 21:39, David Boersma <david.boersma at acmit.at> wrote:
>
>> Hi Jan,
>>
>>
>> It's possible that the problem reported by you is related to issue #350
>> on GitHub: https://github.com/OpenGATE/Gate/issues/350 and hopefully
>> this will be fixed in the not too distant future. :)
>>
>>
>> HTH, David
>> ------------------------------
>> *Von:* Gate-users <gate-users-bounces at lists.opengatecollaboration.org>
>> im Auftrag von Jan Gajewski <jan.gajewski at ifj.edu.pl>
>> *Gesendet:* Montag, 31. August 2020 10:00:22
>> *An:* gate-users at lists.opengatecollaboration.org
>> *Betreff:* [Gate-users] Gate 9.0 does not fill Hits Tree
>>
>> Hi gate users,
>>
>> I have a problem with the newest gate 9 version. The Hits Tree is not
>> filled when any other actor is attached to the sensitive volume.
>>
>> A simple example below. I am defining a scanner with a single crystal and
>> I wanted to attach a phase space actor to it. The phase space output is
>> correct but the Hits Tree is empty. The Hits Tree is correct only if no
>> actors are attached to the crystal or any ancestor.
>> I found out that attaching even SimulationStatisticActor to the crystal
>> or to World causes that the Hits Tree is empty. I suppose that any other
>> actor attached to the crystal leads to the same. It is the same regardless
>> of the form of output definition ( /gate/output/root/enable vs
>> /gate/output/tree/enable).
>>
>> #=====================================================
>> # VISUALISATION
>> #=====================================================
>> /vis/disable
>>
>> #=====================================================
>> # VERBOSITY
>> #=====================================================
>> /gate/verbose Physic 0
>> /gate/verbose Step 0
>> /gate/verbose Cuts 0
>> /gate/verbose SD 0
>> /gate/verbose Actions 0
>> /gate/verbose Actor 0
>> /gate/verbose Error 0
>> /gate/verbose Warning 0
>> /gate/verbose Output 0
>> /gate/verbose Beam 0
>> /gate/verbose Volume 0
>> /gate/verbose Image 0
>> /gate/verbose Geometry 0
>>
>> #=====================================================
>> # MATERIALS
>> #=====================================================
>> /gate/geometry/setMaterialDatabase GateMaterials.db
>>
>> #=====================================================
>> # GEOMETRY
>> #=====================================================
>> # WORLD INITIALIZATION
>> /gate/world/geometry/setXLength 10 cm
>> /gate/world/geometry/setYLength 10 cm
>> /gate/world/geometry/setZLength 10 cm
>> /gate/world/setMaterial Vacuum
>>
>> # scanner
>> /gate/world/daughters/name scanner
>> /gate/world/daughters/insert box
>> /gate/scanner/setMaterial Vacuum
>> /gate/scanner/geometry/setXLength 20 mm
>> /gate/scanner/geometry/setYLength 20 mm
>> /gate/scanner/geometry/setZLength 10 mm
>> /gate/scanner/vis/setColor blue
>> /gate/scanner/vis/forceWireframe
>>
>> # detector
>> /gate/scanner/daughters/name TPX
>> /gate/scanner/daughters/insert box
>> /gate/TPX/setMaterial Silicon
>> /gate/TPX/placement/setRotationAxis 1 0 0
>> /gate/TPX/geometry/setXLength 15 mm
>> /gate/TPX/geometry/setYLength 15 mm
>> /gate/TPX/geometry/setZLength 5 mm
>> /gate/TPX/vis/setColor yellow
>> /gate/TPX/vis/forceSolid 0
>> /gate/systems/scanner/level1/attach TPX
>> /gate/TPX/attachCrystalSD
>>
>>
>> #=====================================================
>> # PHYSICS
>> #=====================================================
>> # physics list with electromagnetic opt. 4
>> /gate/physics/addPhysicsList QGSP_BIC_HP_EMZ
>>
>> #=====================================================
>> # ACTORS
>> #=====================================================
>> # phase space actor in TPX
>> /gate/actor/addActor PhaseSpaceActor psActorTPX
>> /gate/actor/psActorTPX/save
>> Results/psActorTPX.root
>> /gate/actor/psActorTPX/attachTo TPX
>>
>> # statistics actor
>> /gate/actor/addActor SimulationStatisticActor simStat
>> # /gate/actor/simStat/attachTo world
>> /gate/actor/simStat/saveEveryNSeconds 60
>> /gate/actor/simStat/save Results/simStat.out
>>
>> #=====================================================
>> # OUTPUT
>> #=====================================================
>> # /gate/output/root/enable
>> # /gate/output/root/setFileName Results/avTPXHits
>> # /gate/output/root/setRootHitFlag 1
>> # /gate/output/root/setRootSinglesFlag 0
>> # /gate/output/root/setRootNtupleFlag 0
>>
>> /gate/output/tree/enable
>> # /gate/output/tree/addFileName Results/avTPXHits.npy
>> /gate/output/tree/addFileName Results/avTPXHits.root
>> # /gate/output/tree/hits/branches/trackLocalTime/disable
>> # /gate/output/tree/hits/branches/volumeIDs/disable
>> /gate/output/tree/hits/enable
>>
>> #=====================================================
>> # INITIALISATION
>> #=====================================================
>> /gate/run/initialize
>> #/gate/physics/print Results/physics.txt
>>
>> #=====================================================
>> # SOURCE
>> #=====================================================
>> /gate/source/addSource beam PencilBeam
>> /gate/source/beam/setParticleType proton
>> /gate/source/beam/setEnergy 50 MeV
>> /gate/source/beam/setSigmaEnergy 1 MeV
>> /gate/source/beam/setPosition 0. 0. -200 cm
>> /gate/source/beam/setSigmaX 13 mm
>> /gate/source/beam/setSigmaY 7.1 mm
>> /gate/source/beam/setSigmaTheta 6 mrad
>> /gate/source/beam/setSigmaPhi 4 mrad
>> /gate/source/beam/setEllipseXThetaEmittance 77 mm*mrad
>> /gate/source/beam/setEllipseXThetaRotationNorm positive
>> /gate/source/beam/setEllipseYPhiEmittance 50 mm*mrad
>> /gate/source/beam/setEllipseYPhiRotationNorm positive
>>
>>
>> #=====================================================
>> # RANDOM
>> #=====================================================
>> /gate/random/setEngineName MersenneTwister
>> /gate/random/setEngineSeed 591460577
>> /gate/random/verbose 0
>>
>> #=====================================================
>> # ACQUISITION SETTINGS
>> #=====================================================
>> /gate/application/setTotalNumberOfPrimaries 1E2
>> /gate/application/setTimeStart 0. s
>> /gate/application/setTimeStop 100. s
>> /gate/application/startDAQ
>>
>>
>> Thanks in advance for any help.
>>
>>
>> ---
>> Jan Gajewski
>>
>> *jan.gajewski at ifj.edu.pl <jan.gajewski at ifj.edu.pl> **Institute of
>> Nuclear Physics* (IFJ PAN)
>> Radzikowskiego 152, 31-342 Kraków, Poland
>> phone: (+48) 12 662 8331
>>
>> _______________________________________________
>> Gate-users mailing list
>> Gate-users at lists.opengatecollaboration.org
>> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
>
>
> --
> David Sarrut, Phd
> Directeur de recherche CNRS
> CREATIS, UMR CNRS 5220, Inserm U1206
> Centre de lutte contre le cancer Léon Bérard
> 28 rue Laënnec, 69373 Lyon cedex 08
> Tel : 04 78 78 51 51 / 06 74 72 05 42
> http://www.creatis.insa-lyon.fr/~dsarrut
> _________________________________
> "2 + 2 = 5, for extremely large values of 2"
> _________________________________
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200904/e8d85558/attachment-0001.html>
More information about the Gate-users
mailing list