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],