GWのガイガーカウンタ自作チャレンジ（２台目） ネクスト（１２） 1台目の動作テスト

MacBook Proで、１分間、32Kbpsの音声録音して、非圧縮のAIFFファイルに保存。この音声ファイルから、プログラムでパルスを計数する。この前、使ったプログラムを更に修正して、CK1026でのCPMとμsV/hの換算は、このサイトのデータを使用して、3100CMP=1μsV/hに修正。パルス電圧の閾値も、修正。

【パルス計数プログラム】

#include 
#include 
#include 

struct HEADER_CNK {
	u_int32_t ckID;
	int32_t ckSize;
	u_int32_t formType;
};

struct HEADER_COM {
	u_int32_t ckID;
	int32_t ckSize;
	short numChs;
	u_int32_t numSam;
	short samSize;
	char samRate[10];
};

struct HEADER_SSND {
	u_int32_t ckID;
	int32_t ckSize;
	u_int32_t oft;
	u_int32_t blkSize;
};

main( int argc, char **argv )
{
	int fp;
	int ret;
	struct HEADER_CNK hCnk;
	struct HEADER_COM hCom;
	struct HEADER_SSND hSsnd;
	char sd;
	char sd_prev;
	char sd_prev_prev;
	char sd_max;
	char sd_min;

	double svh = 0.0;

	int cnt = 0;
	sd_max = sd_min = 0;

	fp = open( argv[1], O_RDONLY );
	if (fp == -1) {
		fprintf( stderr, "ERROR(%d), %s\n", __LINE__, argv[1] );
		return -1;
	}

	ret = read( fp, &hCnk, sizeof(hCnk) );
	ret = read( fp, &hCom, sizeof(hCom) );
	ret = read( fp, &hSsnd, sizeof(hSsnd) );
	sd = sd_prev = sd_prev_prev = 0;

	for ( cnt=0;; ) {
		ret = read( fp, &sd, sizeof(sd)  );
		if ( ret < 1 ) break;

		if ( sd > sd_max ) sd_max = sd;
		if ( sd < sd_min ) sd_min = sd;

		if ( sd < 0 ) sd =0;
		if ( (sd > 20 ) && (sd_prev > 0) && (sd_prev_prev == 0) ) 
		{
			cnt++;
		} else {
		}
		sd_prev_prev = sd_prev;
		sd_prev = sd;
		
	}
	close( fp );

	svh = (double)cnt / (double)3100.0 ; // CK1026 by http://einstlab.web.fc2.com/geiger/geiger3.html
	fprintf( stderr, "%d, %d, %d, cnt = %d, usv/h = %e\n", 
		hSsnd.ckSize, hSsnd.oft , hCom.samSize, cnt, svh );

	fprintf( stderr, "Max: %d, Min: %d\n", sd_max, sd_min );
}

【実行結果例】

luna:GM_LOG digi_ponta$ ./a.out 11060501.aiff 
2164, 0, 3392, cnt = 163, usv/h = 5.258065e-02
Max: 58, Min: -12
luna:GM_LOG digi_ponta$ 

0.053μsV/hと出てる。