fisx 1.1.7

许可证

此代码在麻省理工学院的许可证下转发，详见许可证文件。

安装

要安装python库，只需使用pip install fisx。如果要从代码源存储库中构建python使用的库，只需使用pip install .或python setup.py install方法之一。安装cython>；=0.17很方便（但不是强制性的）。

测试

要在安装后运行测试，请运行：

python -m fisx.tests.testAll

示例

有一个web application正在使用此库计算预期的X射线计数率。

这段python代码展示了如何通过其python绑定使用库。

fromfisximportElementsfromfisximportMaterialfromfisximportDetectorfromfisximportXRFelementsInstance=Elements()elementsInstance.initializeAsPyMca()# After the slow initialization (to be made once), the rest is fairly fast.xrf=XRF()xrf.setBeam(16.0)# set incident beam as a single photon energy of 16 keVxrf.setBeamFilters([["Al1",2.72,0.11,1.0]])# Incident beam filters# Steel composition of Schoonjans et al, 2012 used to generate table Isteel={"C":0.0445,"N":0.04,"Si":0.5093,"P":0.02,"S":0.0175,"V":0.05,"Cr":18.37,"Mn":1.619,"Fe":64.314,# calculated by subtracting the sum of all other elements"Co":0.109,"Ni":12.35,"Cu":0.175,"As":0.010670,"Mo":2.26,"W":0.11,"Pb":0.001}SRM_1155=Material("SRM_1155",1.0,1.0)SRM_1155.setComposition(steel)elementsInstance.addMaterial(SRM_1155)xrf.setSample([["SRM_1155",1.0,1.0]])# Sample, density and thicknessxrf.setGeometry(45.,45.)# Incident and fluorescent beam anglesdetector=Detector("Si1",2.33,0.035)# Detector Material, density, thicknessdetector.setActiveArea(0.50)# Area and distance in consistent unitsdetector.setDistance(2.1)# expected cm2 and cm.xrf.setDetector(detector)Air=Material("Air",0.0012048,1.0)Air.setCompositionFromLists(["C1","N1","O1","Ar1","Kr1"],[0.0012048,0.75527,0.23178,0.012827,3.2e-06])elementsInstance.addMaterial(Air)xrf.setAttenuators([["Air",0.0012048,5.0,1.0],["Be1",1.848,0.002,1.0]])# Attenuatorsfluo=xrf.getMultilayerFluorescence(["Cr K","Fe K","Ni K"],elementsInstance,secondary=2,useMassFractions=1)print("Element   Peak          Energy       Rate      Secondary  Tertiary")forkeyinfluo:forlayerinfluo[key]:peakList=list(fluo[key][layer].keys())peakList.sort()forpeakinpeakList:# energy of the peakenergy=fluo[key][layer][peak]["energy"]# expected measured raterate=fluo[key][layer][peak]["rate"]# primary photons (no attenuation and no detector considered)primary=fluo[key][layer][peak]["primary"]# secondary photons (no attenuation and no detector considered)secondary=fluo[key][layer][peak]["secondary"]# tertiary photons (no attenuation and no detector considered)tertiary=fluo[key][layer][peak].get("tertiary",0.0)# correction due to secondary excitationenhancement2=(primary+secondary)/primaryenhancement3=(primary+secondary+tertiary)/primaryprint("%s%s%.4f%.3g%.5g%.5g"% \
                               (key,peak+(13-len(peak))*" ",energy,rate,enhancement2,enhancement3))