makePolycone Namespace Reference

Functions
def	usage ()
def	ReadPoints (fileName)
def	FindRadius (z, points)
def	FindRange (inner, outer, step)
def	SegmentIsOK (p1, p2, middle, thres)

Variables
string	vectorName = "polyconeShape"
float	boundaryStep = 0.3
float	precision = 0.3
float	minimumSegment = 0.5
	innerPoints = ReadPoints(args[0])
	outerPoints = ReadPoints(args[1])
list	boundary = []
	z = boundaryStep*iZ
	inner = FindRadius(z,innerPoints)
	outer = FindRadius(z,outerPoints)
list	savedPoints = [boundary[0]]
list	runningPoints = []

Function Documentation

def makePolycone.FindRadius	(		z,
			points
	)

Definition at line 84 of file makePolycone.py.

def FindRadius(z,points):
    radii=[]
    there=points[-1]
    for here in points:
        if ((there[0] <= z and z < here[0])
            or (here[0] <= z and z < there[0])):
            r = (here[1]-there[1])/(here[0]-there[0])*(z-there[0]) + there[1]
            radii.append(r)
        there = here
    if len(radii) == 0: return 0.0
    average = reduce(lambda x, y: x+y, radii)/len(radii)
    r = 0.0
    for p in radii:
        r = r + abs(p-average)
    return r/len(radii)

# Find the total range along the vertical axis of the cryostat that
# needs to be used to define the polycone.

def makePolycone.FindRange	(		inner,
			outer,
			step
	)

Definition at line 102 of file makePolycone.py.

def FindRange(inner,outer,step):
    minZ = min(innerPoints[0][0],outerPoints[0][0])
    maxZ = max(innerPoints[0][0],outerPoints[0][0])
    for p in inner:
        minZ = min(p[0], minZ)
        maxZ = max(p[0], maxZ)
    for p in outer:
        minZ = min(p[0], minZ)
        maxZ = max(p[0], maxZ)
    return xrange(int(minZ/step)-1, int(maxZ/step)+2)

# Check to see if a new polycone segment needs to be started.

def makePolycone.ReadPoints

(

fileName

)

Definition at line 68 of file makePolycone.py.

def ReadPoints(fileName):
    print "// Filename: ", fileName
    points = []
    for line in open(fileName):
        if line[0] == "/":
            print line.rstrip()
            continue
        if len(line) < 5: continue
        spl = line.split()
        points.append([float(spl[0]),float(spl[1])])
    return points

# Find the radius of the cryostat at a particular vertical position.
# This is called "z" since it's going to be used to buill a G4Polycone
# along the Z axis.  In the Z axis here is equivalent to the X axis in
# the input file.

def makePolycone.SegmentIsOK	(		p1,
			p2,
			middle,
			thres
	)

Definition at line 114 of file makePolycone.py.

def SegmentIsOK(p1,p2,middle,thres):
    if abs(p1[0]-p2[0]) < minimumSegment: return True;
    maxThres = 0.0
    for p in middle:
        r1 = (p[0]-p2[0])*(p1[1]-p2[1])/(p1[0]-p2[0]) + p2[1]
        r2 = (p[0]-p2[0])*(p1[2]-p2[2])/(p1[0]-p2[0]) + p2[2]
        maxThres = max(maxThres, abs(r1-p[1]))
        maxThres = max(maxThres, abs(r2-p[2]))
        if maxThres > thres: return False;
    return True

def makePolycone.usage

(

)

Definition at line 7 of file makePolycone.py.

def usage():
    print """
    makePolycone <inner-surface.txt> <outer-surface.txt>

    Outputs text to the standard output that is suitable to be included
    into a C++ source file defining the inner and outer surfaces of a
    G4Polycone.

    --vector-name <name>
           The name of a vector to be used in the include file.

    --step
           The step size along the Z axis to use while building the
           polycone.  This can be small (but it will make this run slowly),
           since only the steps where the precision is bad will be saved.

    --precision
           Set the maximum distance the wall of the cryostat can be from the
           interpolated line.
        """
    
try:
    options, args = getopt.gnu_getopt(sys.argv[1:], 
                                      "h",
                                      ["vector-name=",
                                       "precision=",
                                       "step=",
                                       "minimum-segment=",
                                       "help"])
except:
    usage()
    print "ERROR: Invalid option."
    sys.exit(1)

Variable Documentation

list makePolycone.boundary = []

Definition at line 151 of file makePolycone.py.

makePolycone.boundaryStep = 0.3

Definition at line 47 of file makePolycone.py.

makePolycone.inner = FindRadius(z,innerPoints)

Definition at line 154 of file makePolycone.py.

makePolycone.innerPoints = ReadPoints(args[0])

Definition at line 144 of file makePolycone.py.

makePolycone.minimumSegment = 0.5

Definition at line 49 of file makePolycone.py.

makePolycone.outer = FindRadius(z,outerPoints)

Definition at line 155 of file makePolycone.py.

makePolycone.outerPoints = ReadPoints(args[1])

Definition at line 146 of file makePolycone.py.

makePolycone.precision = 0.3

Definition at line 48 of file makePolycone.py.

list makePolycone.runningPoints = []

Definition at line 163 of file makePolycone.py.

list makePolycone.savedPoints = [boundary[0]]

Definition at line 162 of file makePolycone.py.

makePolycone.vectorName = "polyconeShape"

Definition at line 46 of file makePolycone.py.

makePolycone.z = boundaryStep*iZ

Definition at line 153 of file makePolycone.py.