def test_transmission(self): GPSANS() SetSampleDetectorDistance(6000) DataPath(TEST_DIR) DirectBeamCenter("BioSANS_empty_cell.xml") TimeNormalization() DirectBeamTransmission(sample_file="BioSANS_sample_trans.xml", empty_file="BioSANS_empty_trans.xml") AzimuthalAverage(binning="0.01,0.001,0.11", error_weighting=True) AppendDataFile("BioSANS_test_data.xml") Reduce1D() ws = AnalysisDataService.retrieve("BioSANS_test_data_Iq") data = ws.dataY(0) check = [ 0.514758, 0.520759, 0.516451, 0.51932, 0.524206, 0.519275, 0.520125, 0.518997, 0.518729, 0.516198, 0.518718, 0.51072, 0.512816, 0.513449, 0.509453, 0.502968, 0.503003, 0.505098, 0.503835, 0.503088, 0.500716, 0.499304, 0.495777, 0.49332, 0.489926, 0.486497, 0.491656, 0.491858, 0.48546, 0.490808, 0.483111, 0.482176, 0.482359, 0.483098, 0.477528, 0.474279, 0.474485, 0.468472, 0.470305, 0.467228, 0.467934, 0.467971, 0.465358, 0.463885, 0.461762, 0.462352, 0.461285, 0.458322, 0.458118, 0.453064, 0.44927, 0.448151, 0.446129, 0.444207, 0.442629, 0.439792, 0.437958, 0.434826, 0.434443, 0.432655, 0.430731, 0.428771, 0.425893, 0.424477, 0.422421, 0.419886, 0.416942, 0.415876, 0.414037, 0.41339, 0.414353, 0.409062, 0.408431, 0.408712, 0.419282, 0.412833, 0.41062, 0.414427, 0.400056, 0.400141, 0.394724, 0.393821, 0.390721, 0.38932, 0.387497, 0.383062, 0.381603, 0.380016, 0.374635, 0.374214, 0.369733, 0.370353, 0.366464, 0.364109, 0.362184, 0.361299, 0.355246, 0.354339, 0.352412, 0.349748 ] deltas = map(_diff_iq, data, check) delta = reduce(_add, deltas) / len(deltas) self.assertTrue(math.fabs(delta) < 0.001)
def test_reduction_2(self): GPSANS() SetSampleDetectorDistance(6000) DataPath(TEST_DIR) DirectBeamCenter("BioSANS_empty_cell.xml") AppendDataFile("BioSANS_test_data.xml") DarkCurrent("BioSANS_dark_current.xml") AzimuthalAverage(binning="0.01,0.001,0.11", error_weighting=True) Reduce1D() ws = AnalysisDataService.retrieve("BioSANS_test_data_Iq") data = ws.dataY(0) check = [0.268942, 0.272052, 0.269806, 0.27129, 0.273852, 0.271301, 0.271732, 0.271103, 0.270996, 0.269677, 0.27098, 0.266802, 0.26789, 0.268222, 0.266125, 0.262736, 0.262752, 0.263827, 0.26315, 0.262775, 0.261541, 0.260818, 0.258955, 0.257675, 0.255908, 0.254088, 0.256778, 0.256883, 0.253568, 0.25636, 0.252323, 0.251833, 0.251914, 0.252298, 0.249375, 0.247718, 0.247768, 0.244636, 0.245604, 0.243996, 0.244332, 0.244363, 0.242985, 0.242234, 0.241118, 0.241411, 0.24084, 0.239293, 0.2392, 0.236565, 0.234557, 0.233974, 0.232905, 0.231898, 0.231085, 0.229586, 0.22862, 0.227001, 0.226783, 0.225837, 0.224835, 0.223807, 0.222296, 0.221557, 0.220464, 0.219139, 0.217611, 0.217049, 0.21606, 0.215739, 0.216233, 0.213467, 0.213141, 0.213275, 0.219695, 0.216121, 0.215502, 0.21792, 0.209364, 0.209368, 0.2064, 0.205844, 0.20431, 0.203443, 0.202442, 0.200195, 0.199408, 0.19853, 0.195654, 0.195514, 0.193086, 0.193388, 0.19137, 0.190122, 0.189119, 0.18864, 0.185473, 0.184958, 0.183981, 0.182581] deltas = list(map(_diff_iq, data, check)) delta = reduce(_add, deltas) / len(deltas) self.assertLess(math.fabs(delta), 0.00001)
def test_transmission_by_hand_w_sensitivity(self): GPSANS() SetSampleDetectorDistance(6000) DataPath(self._work_dir) DirectBeamCenter("BioSANS_empty_cell.xml") AppendDataFile("BioSANS_test_data.xml") SetTransmission(0.51944, 0.011078) SensitivityCorrection("BioSANS_flood_data.xml") AzimuthalAverage(binning="0.01,0.001,0.11", error_weighting=True) Reduce1D() ws = AnalysisDataService.retrieve("BioSANS_test_data_Iq") data = ws.dataY(0) check = [ 0.374914, 0.393394, 0.414756, 0.443152, 0.460175, 0.456802, 0.477264, 0.478456, 0.488523, 0.489758, 0.489871, 0.487127, 0.497585, 0.499346, 0.49526, 0.489273, 0.486082, 0.486923, 0.494208, 0.494531, 0.492264, 0.494608, 0.478766, 0.487872, 0.48357, 0.474654, 0.48052, 0.483367, 0.485269, 0.480079, 0.482254, 0.47413, 0.48245, 0.471207, 0.476589, 0.474701, 0.472014, 0.465479, 0.468236, 0.462524, 0.46773, 0.458851, 0.457653, 0.461929, 0.465216, 0.451887, 0.45733, 0.450281, 0.45045, 0.447508, 0.446209, 0.445063, 0.446328, 0.445735, 0.444096, 0.438758, 0.43707, 0.432302, 0.437903, 0.430176, 0.426317, 0.427858, 0.433131, 0.423087, 0.418146, 0.421584, 0.417606, 0.420891, 0.414255, 0.412448, 0.413393, 0.405706, 0.404541, 0.403016, 0.402806, 0.400023, 0.394248, 0.396725, 0.389808, 0.386475, 0.386525, 0.386674, 0.386575, 0.382081, 0.373986, 0.370391, 0.37367, 0.373686, 0.367479, 0.36732, 0.36031, 0.366588, 0.362994, 0.361712, 0.351433, 0.349867, 0.3479, 0.352355, 0.344987, 0.340605 ] # Check that I(q) is the same for both data sets deltas = list(map(_diff_iq, data, check)) delta = reduce(_add, deltas) / len(deltas) self.assertLess(math.fabs(delta), 0.00001)