def compute(self): try: sys.stdout = EmittingStream(textWritten=self.writeStdOut) self.checkFields() if not self.DESCRIPTOR == 18: # GRAPHITE tmp = bragg(interactive=False, DESCRIPTOR=self.crystals[self.DESCRIPTOR], H_MILLER_INDEX=self.H_MILLER_INDEX, K_MILLER_INDEX=self.K_MILLER_INDEX, L_MILLER_INDEX=self.L_MILLER_INDEX, TEMPERATURE_FACTOR=self.TEMPERATURE_FACTOR, E_MIN=self.E_MIN, E_MAX=self.E_MAX, E_STEP=self.E_STEP, SHADOW_FILE=congruence.checkFileName( self.SHADOW_FILE)) else: OWxsh_bragg.new_bragg( H_MILLER_INDEX=self.H_MILLER_INDEX, K_MILLER_INDEX=self.K_MILLER_INDEX, L_MILLER_INDEX=self.L_MILLER_INDEX, TEMPERATURE_FACTOR=self.TEMPERATURE_FACTOR, E_MIN=self.E_MIN, E_MAX=self.E_MAX, E_STEP=self.E_STEP, SHADOW_FILE=congruence.checkFileName(self.SHADOW_FILE)) self.send("PreProcessor_Data", ShadowPreProcessorData(bragg_data_file=self.SHADOW_FILE)) except Exception as exception: QMessageBox.critical(self, "Error", str(exception), QMessageBox.Ok) if self.IS_DEVELOP: raise exception
def compute(self): try: sys.stdout = EmittingStream(textWritten=self.writeStdOut) self.checkFields() tmp = bragg(interactive=False, DESCRIPTOR=self.crystals[self.DESCRIPTOR], H_MILLER_INDEX=self.H_MILLER_INDEX, K_MILLER_INDEX=self.K_MILLER_INDEX, L_MILLER_INDEX=self.L_MILLER_INDEX, TEMPERATURE_FACTOR=self.TEMPERATURE_FACTOR, E_MIN=self.E_MIN, E_MAX=self.E_MAX, E_STEP=self.E_STEP, SHADOW_FILE=congruence.checkFileName(self.SHADOW_FILE)) self.send("PreProcessor_Data", ShadowPreProcessorData(bragg_data_file=self.SHADOW_FILE)) except Exception as exception: QMessageBox.critical(self, "Error", str(exception), QMessageBox.Ok)
def test_kb_with_external_errors(): # # initialize shadow3 source (oe0) and beam # beam = Shadow.Beam() oe0 = Shadow.Source() # # preprocessors # prerefl(interactive=0, SYMBOL="Pt", DENSITY=21.45, FILE="Pt5_55.dat", E_MIN=5000.0, E_MAX=55000.0, E_STEP=100.0) bragg(interactive=0, DESCRIPTOR="Si", H_MILLER_INDEX=1, K_MILLER_INDEX=1, L_MILLER_INDEX=1, TEMPERATURE_FACTOR=1.0, E_MIN=5000.0, E_MAX=53000.0, E_STEP=50, SHADOW_FILE="si5_55.111") dm = dabam.dabam() dm.set_input_shadowCalc(1) dm.set_input_shadowNx(10) dm.set_input_shadowWidth(10.0) dm.load(21) oe0.BENER = 6.03999996 oe0.EPSI_X = 5.5200001e-08 oe0.EPSI_Z = 1.07999998e-09 oe0.FDISTR = 4 oe0.FSOURCE_DEPTH = 4 oe0.F_COLOR = 3 oe0.F_PHOT = 0 oe0.HDIV1 = 0.00079999998 oe0.HDIV2 = 0.00079999998 oe0.ISTAR1 = 567656 oe0.NCOL = 0 oe0.NPOINT = 500000 oe0.N_COLOR = 0 oe0.PH1 = 10000.0 oe0.PH2 = 10010.0 oe0.POL_DEG = 0.0 oe0.R_ALADDIN = 2685.56747 oe0.R_MAGNET = 26.8556747 oe0.SIGDIX = 0.0 oe0.SIGDIZ = 0.0 oe0.SIGMAX = 0.00230000005 oe0.SIGMAY = 0.0 oe0.SIGMAZ = 0.000360000005 oe0.VDIV1 = 2.20000002e-05 oe0.VDIV2 = 2.20000002e-05 oe0.WXSOU = 0.0 oe0.WYSOU = 0.0 oe0.WZSOU = 0.0 bl = Shadow.CompoundOE(name="BL") bl.append_monochromator_double_crystal(p0=3153, q0=503, photon_energy_ev=10005, separation=40.0, reflectivity_file="si5_55.111") bl.append_kb(p0=498.0, q0=95.0, separation=30.0, focal_positions=[4154.0, 95.0], dimensions1=[6.0, 30.0], dimensions2=[6.0, 30.0], reflectivity_kind=[1, 1], reflectivity_files=["Pt5_55.dat", "Pt5_55.dat"], surface_error_files=["Shadow.dat", "Shadow.dat"]) # surface_error_files=["/mnt/scisoft/users/srio/Working/rt/CRG-LISA/MICROFUOCO_v01/mirror.dat", # "/mnt/scisoft/users/srio/Working/rt/CRG-LISA/MICROFUOCO_v01/mirror.dat"]) beam.genSource(oe0) beam.traceCompoundOE(bl) return beam
def test_kb_with_external_errors(): # # initialize shadow3 source (oe0) and beam # beam = Shadow.Beam() oe0 = Shadow.Source() # # preprocessors # prerefl(interactive=0,SYMBOL="Pt",DENSITY=21.45,FILE="Pt5_55.dat",E_MIN=5000.0,E_MAX=55000.0,E_STEP=100.0) bragg(interactive=0,DESCRIPTOR="Si",H_MILLER_INDEX=1,K_MILLER_INDEX=1,L_MILLER_INDEX=1,TEMPERATURE_FACTOR=1.0, E_MIN=5000.0,E_MAX=53000.0,E_STEP=50,SHADOW_FILE="si5_55.111") dm = dabam.dabam() dm.set_input_shadowCalc(1) dm.set_input_shadowNx(10) dm.set_input_shadowWidth(10.0) dm.load(21) oe0.BENER = 6.03999996 oe0.EPSI_X = 5.5200001e-08 oe0.EPSI_Z = 1.07999998e-09 oe0.FDISTR = 4 oe0.FSOURCE_DEPTH = 4 oe0.F_COLOR = 3 oe0.F_PHOT = 0 oe0.HDIV1 = 0.00079999998 oe0.HDIV2 = 0.00079999998 oe0.ISTAR1 = 567656 oe0.NCOL = 0 oe0.NPOINT = 500000 oe0.N_COLOR = 0 oe0.PH1 = 10000.0 oe0.PH2 = 10010.0 oe0.POL_DEG = 0.0 oe0.R_ALADDIN = 2685.56747 oe0.R_MAGNET = 26.8556747 oe0.SIGDIX = 0.0 oe0.SIGDIZ = 0.0 oe0.SIGMAX = 0.00230000005 oe0.SIGMAY = 0.0 oe0.SIGMAZ = 0.000360000005 oe0.VDIV1 = 2.20000002e-05 oe0.VDIV2 = 2.20000002e-05 oe0.WXSOU = 0.0 oe0.WYSOU = 0.0 oe0.WZSOU = 0.0 bl = Shadow.CompoundOE(name="BL") bl.append_monochromator_double_crystal(p0=3153,q0=503,photon_energy_ev=10005,separation=40.0, reflectivity_file="si5_55.111") bl.append_kb(p0=498.0,q0=95.0,separation=30.0,focal_positions=[4154.0,95.0], dimensions1=[6.0,30.0],dimensions2=[6.0,30.0],reflectivity_kind=[1,1], reflectivity_files=["Pt5_55.dat", "Pt5_55.dat"], surface_error_files=["Shadow.dat", "Shadow.dat"]) # surface_error_files=["/mnt/scisoft/users/srio/Working/rt/CRG-LISA/MICROFUOCO_v01/mirror.dat", # "/mnt/scisoft/users/srio/Working/rt/CRG-LISA/MICROFUOCO_v01/mirror.dat"]) beam.genSource(oe0) beam.traceCompoundOE(bl) return beam
def write_bragg_preprocessor_file(): bragg(interactive=False, DESCRIPTOR="Si",H_MILLER_INDEX=1,K_MILLER_INDEX=1,L_MILLER_INDEX=1, TEMPERATURE_FACTOR=1.0,E_MIN=5000.0,E_MAX=55000.0,E_STEP=100.0,SHADOW_FILE="si5_55.111")