後期、4個目のSACLAビームタイムが終了

イメージ
厳くんのビームタイム（100テスラXFEL回折）があり、先週の月から金でSACLAに滞在していた。実験は順調に進み、3種類ものサンプルのデータがとれた。大きな変化という意味では一つのサンプルが大成功で、あとはちょっと磁場が足りないことが判明した。データクオリティはめちゃ高かったので雰囲気も良かった。11 keVシードビームの安定性も高い。 今回は今まで自分がやっていた作業も含めて、学生さんたちにかなり任せて活躍してもらった。学生さんたちはとても良くやってくれました。（自分でやるよりもずっと早くて上手だと判明） フルショットを打ちまくって、検出器にデブリがとびまくった。久保田さんが心配そうにしているので、新プロテクターを作らなくてはだめだ。 2.5日のビームタイムで破壊ショット20発。2発はミスでデータなし。前回・前々回に引き続き非常に良いビームタイムでした。PINK-02の実験が共用実験っぽくなってきてよかったよかった。後期で4つめのPINK-SACLAのビームタイムが終了したのでやっと安心。低温・強磁場・X線の準備が整っていたのでうまくいったと思う。年間ではPINK-SACLAが5件。その他、野尻・プラズマも合わせると、7件のSACLA実験があった。それ以外に関西研にも行った。ナノテラスも3月に行くかも。今年度は今までで一番多かった。 来年度は前期に3つ。後期はどうなるか。SPring-8や関西研もありそう。少し時間があったので、SPring-8のビームラインで打ち合わせがてら装置を見たり、兵庫県立大にも短い時間滞在した。効率的に過ごせた。 週末はアートの展覧会にお呼ばれして行き、ちょっとおしゃべりもしてきた。動画を後日見せてもらったが、さすがに、伸びシロしかない・・・。 連休中、調布ではトトロが再上映。家族で観覧した。 ここ1年ほどは映画館で新作・リバイバルを見るのがブーム 国宝 チェンソーマン　レゼ編（2~3回） PERFECT BLUE もののけ姫　IMAX（2回） エヴァンゲリオン　Air/まごころを君に となりのトトロ 時計じかけのオレンジ その他にこどもと おまえうまそうだな すみっコぐらし クリーミーマミ ズートピア2 時計じかけのオレンジに影響されて、赤ワインを飲んだり、聖書を買って読んだりした。 聖書は意外に読みやすく、最後まで目を通した。めっち...
コメントを投稿

Python: Load data function for general files

ほとんどどんなオシロスコープのデータでも読み込めるfunction

import numpy as np
import matplotlib.pyplot as plt
import itertools


def load_scope_data(flist, j):
    # This loads data from oscilloscope
    # It takes file name and returns data list
    # NaN data is ...
    
    fname1 = flist[j]
    for i, k, m in itertools.product(range(100)
                                      , [None, '\t', ',']
                                      , ['utf-8', 'shift-jis']):
        try:
            dat = np.loadtxt(fname1
                             , encoding = m
                             , skiprows = i
                             , delimiter= k).T
        except ValueError:
            # print(f'fail')
            # print(i, k, m)
            pass
        else: 
            print(f'load successful: {fname1}')
            print(f'skiprows = {i}, encoding = {m}, delimiter = {k}')
            print(dat.shape)
            break

    num = len(dat)
    fig, ax = plt.subplots(1, num, tight_layout = True, figsize = (num * 2, 2))
    for i in range(num):
        ax[i].plot(dat[i])
    
    # plt.suptitle(fname1)
    plt.show()
    return dat

コメント

コメントを投稿

このブログの人気の投稿