Esempio n. 1
0
def process_one_evaluation(yi,
                           ZA,
                           libraryName='endl2008.2',
                           libraryVersion='',
                           libPath=default_lib_path,
                           outputDir=None,
                           workDir='work',
                           verbose=False,
                           dryrun=False):
    '''process all reactions for an evaluation that we possibley can'''

    try:
        e = endlProject(database=os.path.join(libPath, libraryName,
                                              libraryVersion),
                        projectile=yi,
                        readOnly=True)
        target = e.readZA(ZA)
    except:
        print 'target za%06d yi%02d doesnt exist' % (ZA, yi)
        return []

    if verbose:
        print '---- process_one_evaluation ----'
        print 'yi:', yi
        print '    za:', ZA

    if verbose: target.ll()
    target.read(I=0)
    CList = filter(lambda x: x in default_C_list, target.CList())

    cutOffEnergy = 0.010  # 10 keV cutoff, if threshold > this value, then don't make a TDF file

    #if CList == []: print "No C's to process for yi =", yi, "za =", ZA,"!"
    reactionnames = []
    for C in CList:
        #target.read( I = 0 )
        #print " yi %02d za%06d  C %02d \n"% (yi,ZA,C)
        if target.findDatas(I=0, C=C)[0].data[0][0] > cutOffEnergy:
            #print  "yi %02d za%06d  C %02d  : "%(yi,ZA,C), "first energy point too high \n"
            continue
        if verbose: print "        C:", C
        try:
            reactionname = processTDF_Reaction(target,
                                               C,
                                               workDir=workDir,
                                               outputDir=outputDir,
                                               libraryName=libraryName,
                                               libraryVersion=libraryVersion,
                                               verbose=verbose,
                                               dryrun=dryrun)
            reactionnames.append(reactionname)
            print 'adding to tdf list :', reactionname
        except RuntimeError, err:
            print err
            print "\n.... continuing\n"
Esempio n. 2
0
def process_one_reaction(yi,
                         ZA,
                         C,
                         libraryName='endl2009.0',
                         libraryVersion='',
                         libPath=default_lib_path,
                         outputDir=None,
                         workDir='work',
                         verbose=False,
                         dryrun=False):
    '''
       process one reaction, use this for testing
    
       projectile designator: yi = ( 1,n; 2,p; 3,d; 4,t; 5,3He; 6,He; 7,g )
       target ZA = 1000 * Z + A                                            
       I = 0 is cross section data.                                        
    '''

    if verbose:
        print '---- process_one_reaction ----'
        print 'yi:', yi
        print '    za:', ZA
        print "        C:", C
    print "       next endlproject call...       "
    e = endlProject(os.path.join(libPath, libraryName, libraryVersion),
                    projectile=yi,
                    readOnly=True)
    target = e.readZA(ZA)
    if verbose: target.ll()
    target.read(I=0)
    reactionname = processTDF_Reaction(target,
                                       C,
                                       workDir=workDir,
                                       outputDir=outputDir,
                                       libraryName=libraryName,
                                       libraryVersion=libraryVersion,
                                       verbose=verbose,
                                       dryrun=dryrun)
    if dryrun:
        return reactionname
Esempio n. 3
0
                    '--output',
                    default='endl2gnd',
                    help='prefix for resulting .endf and .xml files')
parser.add_argument(
    '-v',
    '--version',
    default='1.0.0',
    help=
    "Evaluation version number. Should be of form 'major.minor.patchlevel' (i.e. 2011.0.1)"
)

if (__name__ == '__main__'):

    args = parser.parse_args()

    e = endlProject(args.library, projectile=args.projectile, readOnly=True)

    if args.ZA[-1].isalpha():  # isomer, e.g. 95242m
        za = e.readZA(args.ZA[:-1], suffix=args.ZA[-1])
    else:
        za = e.readZA(args.ZA)
    za.read()
    if not args.includeAverageProductData:
        za.removeFile(I=10)
        za.removeFile(I=13)

    r = za.toGND(evaluationLibrary=os.path.basename(args.library),
                 evaluationVersion=args.version)
    r.saveToFile(args.output + ".xml")

    if (args.toENDF6):
Esempio n. 4
0
 def testReadZA(self):
     e = endlProject(os.sep.join(__path__ + ['test', 'testdb']))
     self.assertEqual(e.ZAList(), ['za001001'])
     za = e.readZA(1001)
     za.read()
     self.assertEqual([x.C for x in za.findDatas(I=0)], [1, 10, 46])
Esempio n. 5
0
 def testInit(self):
     e = endlProject(os.sep.join(__path__ + ['test', 'testdb']))
     self.assertEqual(e.ZAList(), ['za001001'])