def __init__(self, gParms=[], dParms=[]):
"""Constructor for detectorInfo"""
#
if gParms: # maintain this keyword for compatability with Don's usage
self.detector = detector.newDetector('ge', gParms=gParms, dParms=dParms)
else:
self.detector = detector.newDetector('ge')
self.mrbImages = []
self.fitParams = []
#
# Refinement variables
#
self.refineFlags = list(DetectorInfo.DFLT_REFINE_FLAGS)
#
self.calInp = self.rdrInp = None
#
return
def __init__(self, gParms=[], dParms=[]):
"""Constructor for detectorInfo"""
#
if gParms: # maintain this keyword for compatability with Don's usage
self.detector = detector.newDetector('ge', gParms=gParms, dParms=dParms)
else:
self.detector = detector.newDetector('ge')
self.mrbImages = []
self.fitParams = []
#
# Refinement variables
#
self.refineFlags = list(DetectorInfo.DFLT_REFINE_FLAGS)
#
self.calInp = self.rdrInp = None
#
return
def __init__(self):
"""Constructor for detectorInfo"""
#
self.detector = detector.newDetector("ge")
self.mrbImages = []
self.fitParams = []
#
# Refinement variables
#
self.refineFlags = list(DetectorInfo.DFLT_REFINE_FLAGS)
#
self.calInp = self.rdrInp = None
#
return
deltaOme = (omeMax - omeMin) / nOme # may be negative
omegas = num.linspace(omeMin + 0.5 * deltaOme, omeMax - 0.5 * deltaOme,
nOme)
num.random.seed(seed=opts.randomSeed)
######################################################################
# ... add generic detector capability, so that can have smaller number of pixels -- detector.newGenericDetector
# ... add detector geom parameters
#detectorGeom = detector.newDetector('ge')
dgType, argsStr, kwargsStr = opts.detectorGeomStr.split(';')
args = eval(argsStr)
kwargs = eval(kwargsStr)
kwargs.setdefault('readerKWArgs', {"doFlip": False})
detectorGeom = detector.newDetector(dgType, *args, **kwargs)
detectorGeom.workDist = 2000.
tThMax = detectorGeom.getTThMax(func=num.max)
print tThMax / degToRad
planeData.tThMax = tThMax
######################################################################
if opts.tests.count("makeSynthStack") > 0:
tic = time.time()
gList = makeRandGrains(opts.numGrains, omegaMM, detectorGeom)
'''
spotParamsList = [
[(0.5*tThMax,0.5*num.pi,omegas[1]), (0.002, 0.05, 0.05), 1e2], # center, fwhm, A
[(0.5*tThMax,0.0*num.pi,omegas[0]), (0.002, 0.05, 0.05), 1e2],
omegaMM = ( min(omeMin, omeMax), max(omeMin, omeMax) )
deltaOme = (omeMax-omeMin)/nOme # may be negative
omegas = num.linspace(omeMin + 0.5*deltaOme, omeMax - 0.5*deltaOme, nOme)
num.random.seed(seed=opts.randomSeed)
######################################################################
# ... add generic detector capability, so that can have smaller number of pixels -- detector.newGenericDetector
# ... add detector geom parameters
#detectorGeom = detector.newDetector('ge')
dgType, argsStr, kwargsStr = opts.detectorGeomStr.split(';')
args = eval(argsStr); kwargs = eval(kwargsStr);
kwargs.setdefault('readerKWArgs',{"doFlip":False})
detectorGeom = detector.newDetector(dgType, *args, **kwargs)
detectorGeom.workDist = 2000.
tThMax = detectorGeom.getTThMax(func=num.max)
print tThMax/degToRad
planeData.tThMax = tThMax
######################################################################
if opts.tests.count("makeSynthStack") > 0:
tic = time.time()
gList = makeRandGrains(opts.numGrains, omegaMM, detectorGeom)
'''
spotParamsList = [
[(0.5*tThMax,0.5*num.pi,omegas[1]), (0.002, 0.05, 0.05), 1e2], # center, fwhm, A
[(0.5*tThMax,0.0*num.pi,omegas[0]), (0.002, 0.05, 0.05), 1e2],
