Functions
def	local2Global (n)

None	print_info (str path, node, int level, int currentlevel, print_sub_det_info=False)

Variables
ArgumentParser	parser = ArgumentParser()

	help

	required

	default

	action

ArgumentParser	options = parser.parse_args()

ArgumentParser	fname = options.geometry

ROOT	fgeom = ROOT.TFile.Open(fname)

ROOT	fGeo = fgeom.Get("FAIRGeom")

ROOT	top = fGeo.GetTopVolume()

	else :

int	currentlevel = 0

Function Documentation

◆ local2Global()

def getGeoInformation.local2Global ( n )

Definition at line 16 of file getGeoInformation.py.

def local2Global(n):
    Info = {}
    nav = ROOT.gGeoManager.GetCurrentNavigator()
    nav.cd(n)
    Info["node"] = nav.GetCurrentNode()
    Info["path"] = n
    tmp = Info["node"].GetVolume().GetShape()
    Info["material"] = Info["node"].GetVolume().GetMaterial().GetName()
    if options.moreInfo:
        x = ROOT.gGeoManager.GetVerboseLevel()
        ROOT.gGeoManager.SetVerboseLevel(0)
        Info["weight"] = Info["node"].GetVolume().Weight()  # kg
        Info["cubicmeter"] = Info["node"].GetVolume().Capacity() / 1000000.0  #
        ROOT.gGeoManager.SetVerboseLevel(x)
    o = [tmp.GetOrigin()[0], tmp.GetOrigin()[1], tmp.GetOrigin()[2]]
    Info["locorign"] = o
    local = array("d", o)
    globOrigin = array("d", [0, 0, 0])
    nav.LocalToMaster(local, globOrigin)
    Info["origin"] = globOrigin
    shifts = [
        [-tmp.GetDX() + o[0], o[1], o[2]],
        [tmp.GetDX() + o[0], o[1], o[2]],
        [o[0], -tmp.GetDY() + o[1], o[2]],
        [o[0], tmp.GetDY() + o[1], o[2]],
        [o[0], o[1], -tmp.GetDZ() + o[2]],
        [o[0], o[1], tmp.GetDZ() + o[2]],
    ]
    shifted = []
    for s in shifts:
        local = array("d", s)
        glob = array("d", [0, 0, 0])
        nav.LocalToMaster(local, glob)
        shifted.append([glob[0], glob[1], glob[2]])
    Info["boundingbox"] = {}
    for j in range(3):
        jmin = 1e30
        jmax = -1e30
        for s in shifted:
            if s[j] < jmin:
                jmin = s[j]
            if s[j] > jmax:
                jmax = s[j]
        Info["boundingbox"][j] = [jmin, jmax]
    return Info
 
 

◆ print_info()

None getGeoInformation.print_info	(	str	path,
			node,
		int	level,
		int	currentlevel,
			print_sub_det_info = `False`
	)

Definition at line 63 of file getGeoInformation.py.

def print_info(path: str, node, level: int, currentlevel: int, print_sub_det_info=False) -> None:
    sub_nodes = {}
    fullInfo = {}
    for subnode in node.GetNodes():
        name = subnode.GetName()
        fullInfo[name] = local2Global(path + "/" + name)
        sub_nodes[name] = fullInfo[name]["origin"][2]
 
    for name, _ in sorted(sub_nodes.items(), key=operator.itemgetter(1)):
        boundingbox = fullInfo[name]["boundingbox"]
 
        format_string = (
            "{:<28s}: z={:10.4F}cm  dZ={:10.4F}cm  [{:10.4F}   {:10.4F}]"
            + " dx={:10.4F}cm [{:10.4F}   {:10.4F}] dy={:10.4F}cm [{:10.4F}   {:10.4F}] {:>20s}"
        )
 
        format_variable = [
            "   " * int(currentlevel) + name,
            fullInfo[name]["origin"][2],
            abs(boundingbox[2][0] - boundingbox[2][1]) / 2.0,
            boundingbox[2][0],
            boundingbox[2][1],
            abs(boundingbox[0][0] - boundingbox[0][1]) / 2.0,
            boundingbox[0][0],
            boundingbox[0][1],
            abs(boundingbox[1][0] - boundingbox[1][1]) / 2.0,
            boundingbox[1][0],
            boundingbox[1][1],
            fullInfo[name]["material"],
        ]
 
        if options.moreInfo:
            cubicmeter = fullInfo[name]["cubicmeter"]
            weight = fullInfo[name]["weight"]
            format_string += " {:10.1F}kg {:10.1F}m3"
            format_variable.extend([weight, cubicmeter])
 
        print(format_string.format(*format_variable))
 
        if options.volume in ["", name]:
            print_sub_det_info = True
 
        if print_sub_det_info and currentlevel < level and fullInfo[name]["node"].GetNodes():
            print_info(fullInfo[name]["path"], fullInfo[name]["node"], level, currentlevel + 1, print_sub_det_info)
 
        if currentlevel == 0:
            print_sub_det_info = False
 
 

Variable Documentation

◆ action

getGeoInformation.action

Definition at line 116 of file getGeoInformation.py.

◆ currentlevel

int getGeoInformation.currentlevel = 0

Definition at line 136 of file getGeoInformation.py.

◆ default

getGeoInformation.default

Definition at line 114 of file getGeoInformation.py.

◆ else

getGeoInformation.else :

Definition at line 130 of file getGeoInformation.py.

◆ fGeo

ROOT getGeoInformation.fGeo = fgeom.Get("FAIRGeom")

Definition at line 121 of file getGeoInformation.py.

◆ fgeom

ROOT getGeoInformation.fgeom = ROOT.TFile.Open(fname)

Definition at line 120 of file getGeoInformation.py.

◆ fname

ArgumentParser getGeoInformation.fname = options.geometry

Definition at line 119 of file getGeoInformation.py.

◆ help

getGeoInformation.help

Definition at line 113 of file getGeoInformation.py.

◆ options

ArgumentParser getGeoInformation.options = parser.parse_args()

Definition at line 118 of file getGeoInformation.py.

◆ parser

ArgumentParser getGeoInformation.parser = ArgumentParser()

Definition at line 112 of file getGeoInformation.py.

◆ required

getGeoInformation.required