Hi Peter,
You need not to redefine your geometry via the
VirtualMC interface.
When you define your geometry directly via TGeo,
you have just to notify VMC that geometry is
built via Root:
YourMCApplication::ConstructGeometry()
{
// build your geometry via TGeo
// eg. via importing it from a file
gMC->SetRootGeometry();
}
TGeoMCGeometry class has been defined for backward
compatibility - it implements the functions for building
geometry defined in TVirtualMC interface.
It allows users who have already defined their MC application
via VMC interface to move to TGeo geometry without
rewriting their code.
In case of more questions/problems, you can better
use the VMC mailing list: vmc@root.cern.ch.
Best regards,
Ivana Hrivnacova
==============================================================
e-mail: Ivana.Hrivnacova@cern.ch
address: Institut Physique Nucleaire, 91406 Orsay, France
phone: (33) 01 69 15 65 94
==============================================================
On Tue, 13 Apr 2004, Peter Steinberg wrote:
>
> Hi everyone (esp. the VMC team) -
>
> I'm trying to use the geant3 module in ROOT 4-00-03 but am
> having trouble incorporating my own geometry into the examples.
>
> I've already defined a geometry using the TGeoManager class
> and saved it to a file. I've then hacked in my geometry into
> the E01 example in geant4_vmc, by reading it in from a file
> (using TGeoManager::Import), but when i run the example i get
> this output:
>
> ----------------------------------------------------------------------
>
> MZLINK. Initialize Link Area /GCSLNK/ for Store 0 NL/NS= 100 100
> Geant3 has been created.
> Geometry will be defined via TGeo
> ===> TGeo, Root geometry manager created
> Warning in <TGeoManager::dtor>: deleting previous geometry: TGeo/Root
> geometry manager
> ===> , created
> ### Geometry loaded from file...
> Top volume is TOP. Master volume is TOP
> ### Voxelization retrieved from file
> --- node ID tracking disabled
> ### 195 nodes/ 18 volume UID's in Simple geometry
> ----------------modeler ready----------------
>
> Calculating cross section tables, see gphysi.dat for more information
>
> Cross section calculation concluded successfully
> **** GTRIGI: IEVENT= 1 IDEVT= 1 Random Seeds = 0
> 0
> root [1]
>
> ----------------------------------------------------------------------
>
> i.e. it doesn't crash, but nothing happens! this example works fine
> in its original form, so it's really just the introduction of my
> geometry.
>
> my main suspicion is that i really should be using the TGeoMCGeometry
> interface from the beginning, but if i check the source, i see that
> the MCGeometry is using the GeoManager anyway.
>
> so my question is if i really have to rewrite my geometry using
> the TGeoMCGeometry interface, or if there is a way to "import" an
> existing TGeoManager into the TGeoMCGeometry -- or if I am completely
> misunderstanding something!
>
> Thanks,
> Peter
>
>
--
This archive was generated by hypermail 2b29 : Sun Jan 02 2005 - 05:50:07 MET