なかのにっき

sroep [--directory=DIR] --cc=CENTER_CHANNEL --hw=HALF_WIDTH [--debug]

• --directory はデータディレクトリ。省略するとカレントディレクトリ。
• --cc は対象となるチャンネル範囲の中心。
• --hw は対象となるチャンネル範囲の幅。cc-hw〜cc+hw が対象になる。
• --debug を指定すると、チャネル範囲の測定データ・gaussfit のコマンドライン・gaussfit への STDIN と結果の STDOUT, STDERR を標準出力に表示する。

gaussfit へのパラメータを決めているあたりのルーチンは以下の通り。 $imax と $imin は範囲での最大値と最小値、 $CCH はコマンドラインで与えた中心チャンネル、 $sigma は範囲の各チャンネルにおける標準偏差の最大値。

  my $AI = sprintf "%f", $imax;
  my $AU = sprintf "%f", 2 * $imax;
  my $AL = sprintf "%f", 0.3 * $imax + 0.7 * $imin;

  my $CI = sprintf "%f", $CCH;
  my $CL = sprintf "%f", $CCH - 2;
  my $CU = sprintf "%f", $CCH + 2;

  my $BI = sprintf "%f", $imin;
  my $BL = sprintf "%f", $imin - $sigma;
  my $BU = sprintf "%f", 0.7 * $imax + 0.3 * $imin;

  my $WI = sprintf "%f", 1.5;
  my $WL = sprintf "%f", 0.5;
  my $WU = sprintf "%f", 3.0;

W については決め打ち。

適当にデータを拾ってきてためしてみたら、 2σを越えて暗電流より信号が小さいデータがやたらたくさんあるのだが... こいつがシャッターを逆にしてたというやつかなあ。

オプションの処理には Getopt::Long を使っており、 --debug の処理は

$DEBUG = 0;
$result = GetOptions ("dir=s" => \$DIR,
                      "cc=i"  => \$CCH,
                      "hw=i"  => \$HWC,
                      "debug"  => \$DEBUG);
  :
if ($DEBUG) {もげもげ}

という感じ。

そのまま sroep というモジュール名で。

セクションを追加して記述。

「あいあまって」かな？と一瞬迷って辞書を引いた。 「相俟って」らしい。

ワームか。ひさびさでんな。