Skip to content

Instantly share code, notes, and snippets.

@paulromano

paulromano/ReconstructResonances.ipynb

Created March 26, 2014 23:57
Show Gist options
  • Save paulromano/9796630 to your computer and use it in GitHub Desktop.
Save paulromano/9796630 to your computer and use it in GitHub Desktop.
Download ZIP
Display the source blob
Display the rendered blob
Raw
ReconstructResonances.ipynb
{
"metadata": {
"name": ""
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"from fudge.legacy.converting.endfFileToGND import endfFileToGND"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"eval, cov = endfFileToGND('n-094_Pu_239.endf')"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
" total fission neutron data: LNU = 2\n",
" Delayed fission neutron data (MT=455) LDG=0 LNU=2 : MF=5, LF=5 : MF=5, LF=5 : MF=5, LF=5 : MF=5, LF=5"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" : MF=5, LF=5 : MF=5, LF=5\n",
" prompt fission neutron data: LNU = 2\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" WARNING: distributions for MT=3 (nonelastic) are not supported and have been ignored\n",
" 2 [3, 4, 33] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 51 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 52 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 53 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 54 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 55 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 56 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 57 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 58 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 59 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 60 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 61 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 62 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 63 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 64 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 65 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 66 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 67 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 68 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 69 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 70 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 71 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 72 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 73 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 74 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 75 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 76 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 77 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 78 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 79 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 80 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 81 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 82 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 83 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 84 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 85 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 86 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 87 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 88 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 89 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 90 [3, 4] : MF=4, LTT = 1\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 91 [3, 6] : MF=6 : ZAP=1, LAW=1, LANG=2\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 4 [3, 12, 14, 15, 33] : MF=12 LO=1 : ZAP=0 : MF=14 : MF=15 \n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 16 [3, 6, 33] : MF=6 : ZAP=1, LAW=1, LANG=2\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 17 [3, 6, 33] : MF=6 : ZAP=1, LAW=1, LANG=2\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 37 [3, 6, 33] : MF=6 : ZAP=1, LAW=1, LANG=2\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 18 [3, 4, 5, 12, 14, 15, 33, 35] : MF=4, LTT = 0 : MF=5, LF=1 : MF=12 LO=1 : ZAP=0 "
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
": MF=14 : MF=15 \n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
" WARNING: distributions for MT=3 (nonelastic) are not supported and have been ignored\n",
" WARNING: have prompt fission nu_bar so not including total\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" WARNING: have prompt fission nu_bar so not including total\n",
" 19 [3]\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 20 [3]\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 21 [3]\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 38 [3]\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 102 [3, 12, 14, 15, 33] : MF=12 LO=1 : ZAP=0 : MF=14 : MF=15 \n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 1 [3, 33]\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" 3 [3]\n"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
" WARNING: multiple resolved energy intervals are deprecated!\n"
]
},
{
"output_type": "stream",
"stream": "stderr",
"text": [
" WARNING: multiple resolved energy intervals are deprecated!\n"
]
}
],
"prompt_number": 2
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"eval.reconstructResonances()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"ename": "AttributeError",
"evalue": "'module' object has no attribute 'pointwiseXY_C'",
"output_type": "pyerr",
"traceback": [
"\u001b[1;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[1;31mAttributeError\u001b[0m Traceback (most recent call last)",
"\u001b[1;32m<ipython-input-3-9133e57a5f27>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[0meval\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mreconstructResonances\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m",
"\u001b[1;32m/home/romano/.local/lib/python2.7/site-packages/fudge/gnd/reactionSuite.pyc\u001b[0m in \u001b[0;36mreconstructResonances\u001b[1;34m(self, accuracy, verbose)\u001b[0m\n\u001b[0;32m 434\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 435\u001b[0m \u001b[0mbackground\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mchannel\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcrossSection\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mforms\u001b[0m\u001b[1;33m[\u001b[0m \u001b[0mchannel\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mcrossSection\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mnativeData\u001b[0m \u001b[1;33m]\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtabulatedData\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 436\u001b[1;33m \u001b[0mbackground\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mbackground\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtoPointwiseLinear\u001b[0m\u001b[1;33m(\u001b[0m \u001b[0mepsilon\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mepsilon\u001b[0m \u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 437\u001b[0m \u001b[0mRRxsec\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mRRxsec\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtoPointwiseLinear\u001b[0m\u001b[1;33m(\u001b[0m \u001b[0mepsilon\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mepsilon\u001b[0m \u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 438\u001b[0m \u001b[1;31m# before adding to another pointwise region, y-value of upper point must be 0:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/romano/.local/lib/python2.7/site-packages/fudge/gnd/reactionData/crossSection.pyc\u001b[0m in \u001b[0;36mtoPointwiseLinear\u001b[1;34m(self, lowerEps, upperEps)\u001b[0m\n\u001b[0;32m 563\u001b[0m \u001b[0maccuracy\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;36m1e-3\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 564\u001b[0m \u001b[1;32mif\u001b[0m\u001b[1;33m(\u001b[0m \u001b[0mlen\u001b[0m\u001b[1;33m(\u001b[0m \u001b[0mself\u001b[0m \u001b[1;33m)\u001b[0m \u001b[1;33m>\u001b[0m \u001b[1;36m0\u001b[0m \u001b[1;33m)\u001b[0m \u001b[1;33m:\u001b[0m \u001b[0maccuracy\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mgetAccuracy\u001b[0m\u001b[1;33m(\u001b[0m \u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 565\u001b[1;33m \u001b[0mxys\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mregions\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mregionsXYs\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mtoPointwiseLinear\u001b[0m\u001b[1;33m(\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0maccuracy\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mlowerEps\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mupperEps\u001b[0m \u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 566\u001b[0m \u001b[1;32mreturn\u001b[0m\u001b[1;33m(\u001b[0m \u001b[0mlinear\u001b[0m\u001b[1;33m(\u001b[0m \u001b[0mxys\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0maxes\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mxys\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0maccuracy\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mxys\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mgetAccuracy\u001b[0m\u001b[1;33m(\u001b[0m \u001b[1;33m)\u001b[0m \u001b[1;33m)\u001b[0m \u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 567\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;32m/home/romano/.local/lib/python2.7/site-packages/fudge/core/math/xData/regions.pyc\u001b[0m in \u001b[0;36mtoPointwiseLinear\u001b[1;34m(self, accuracy, lowerEps, upperEps, removeOverAdjustedPoints, axes_)\u001b[0m\n\u001b[0;32m 158\u001b[0m \u001b[1;32mraise\u001b[0m \u001b[0mException\u001b[0m\u001b[1;33m(\u001b[0m \u001b[1;34m'Adjustment at %s makes new x = %s >= prior x = %s; eps = %s'\u001b[0m \u001b[1;33m%\u001b[0m \u001b[1;33m(\u001b[0m \u001b[0mx12\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mx\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mx11\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mlowerEps\u001b[0m \u001b[1;33m)\u001b[0m \u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 159\u001b[0m \u001b[1;32melse\u001b[0m \u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 160\u001b[1;33m \u001b[0mxys\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;33m-\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m]\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m \u001b[0mx\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mXYs\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mpointwiseXY_C\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0minterpolatePoint\u001b[0m\u001b[1;33m(\u001b[0m \u001b[0minterpolationStr\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mx\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mx11\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0my11\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mx12\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0my12\u001b[0m \u001b[1;33m)\u001b[0m \u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 161\u001b[0m \u001b[1;32mif\u001b[0m\u001b[1;33m(\u001b[0m \u001b[0mupperEps\u001b[0m \u001b[1;33m!=\u001b[0m \u001b[1;36m0.\u001b[0m \u001b[1;33m)\u001b[0m \u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 162\u001b[0m \u001b[0mx22\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0my22\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mregion_\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m1\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
"\u001b[1;31mAttributeError\u001b[0m: 'module' object has no attribute 'pointwiseXY_C'"
]
}
],
"prompt_number": 3
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment