doPlan = readerPlan.DataPlan(pathPlan,
nStart=nStart,
nTot=nStop,
nPar=512 * 512 * 512)
# read splitter ids and data
splitterIds = np.load(pathCT + "splitterIds.npy", allow_pickle=True)
splitterData = np.load(pathCT + "splitterData.npy", allow_pickle=True)
clumpObjList = []
ctFileNameList = os.listdir(pathCT)
for fileName in ctFileNameList:
if fileName[0:8] != "splitter":
inArr = np.load(pathCT + fileName, allow_pickle=True)
#print(fileName)
#print(inArr)
clumpObjList.append(
readerCt.Clump(inArr, doPlan, splitterIds, splitterData))
################################################################################
# get initial mass spectrum
nSplitters = 0
nStats = 0
mp = []
for clump in clumpObjList:
m0 = clump.massDict[clump.n0]
addToList = True
# apply conditions to clumps to make the list that get their stats taken
#if len(clump.massDict.keys())<10: addToList = False
#if m0 < 1.e-4: addToList = False
#if m0 > 1.0: addToList = False
if clump.splitter:
addToList = False
nSplitters += 1
################################################################################
# read peak files
doPlan = readerPlan.DataPlan(pathPlan, nStart=nStart, nTot=nStop, nPar=512*512*512)
# read splitter ids and data
splitterIds = np.load(pathCT + "splitterIds.npy", allow_pickle=True)
splitterData = np.load(pathCT + "splitterData.npy", allow_pickle=True)
clumpObjList = []
ctFileNameList = os.listdir(pathCT)
for fileName in ctFileNameList:
if fileName[0:8] != "splitter":
inArr = np.load(pathCT + fileName, allow_pickle=True)
#print(fileName)
#print(inArr)
clumpObjList.append(readerCt.Clump(inArr, doPlan, splitterIds, splitterData))
################################################################################
# get initial mass spectrum
nSplitters=0; nStats=0;
mp = []
for clump in clumpObjList:
m0 = clump.massDict[clump.n0]
addToList = True
# apply conditions to clumps to make the list that get their stats taken
#if len(clump.massDict.keys())<10: addToList = False
#if m0 < 1.e-4: addToList = False
#if m0 > 1.0: addToList = False
if clump.splitter: addToList = False; nSplitters+=1;
if addToList: mp.append(m0); nStats+=1;
print("total clumps ever: " + str(len(clumpObjList)))