eic-opticks/ana/gfind.py

0001 #!/usr/bin/env python
0002
0003 import logging, textwrap
0004 log = logging.getLogger(__name__)
0005
0006 from opticks.analytic.gdml import GDML
0007 specs_ = lambda s:filter(lambda s:s[0] != "#", filter(None,textwrap.dedent(s).split("\n")))
0008
0009
0010 if __name__ == '__main__':
0011
0012
0013     lvx = specs_(r"""
0014     PMT_3inch_log
0015     NNVTMCPPMTlMaskVirtual
0016     HamamatsuR12860lMaskVirtual
0017     mask_PMT_20inch_vetolMaskVirtual
0018     """
0019     )
0020     ilv = 3
0021
0022     labels = specs_(r"""
0023     tds_ngt
0024     tds_ngt_pcnk
0025     """
0026     )
0027
0028
0029     g0 = GDML.parse("$OPTICKS_PREFIX/%s.gdml" % labels[0])
0030     g0.smry()
0031     #g0.volume_summary()
0032
0033     lv0 = g0.find_one_volume(lvx[ilv])
0034     print("lv0:%s [%s]" % (lv0, lvx[ilv]) )
0035
0036
0037     if len(labels)>1:
0038         g1 = GDML.parse("$OPTICKS_PREFIX/%s.gdml" % labels[1])
0039         g1.smry()
0040         #g1.volume_summary()
0041
0042         lv1 = g1.find_one_volume(lvx[ilv])
0043         print("lv1:%s" % lv1 )
0044     pass
0045
0046
0047
0048