Skip to content

Instantly share code, notes, and snippets.

@pn11
pn11 / print_tree.cc
Last active Aug 29, 2015
指定したROOTファイルに含まれるTTreeの中身を表示するだけのプログラム。
View print_tree.cc
#include<iostream>
#include<TList.h>
#include<TFile.h>
int main(int argc, char* argv[]){
using namespace std;
TFile *f1 = new TFile(argv[1]);
TList *l1 = f1->GetListOfKeys();
@pn11
pn11 / chi2_test_normalize_check.cc
Created Sep 10, 2014
ROOTのTH1:Chi2Test() (http://root.cern.ch/root/html/TH1.html#TH1:Chi2Test) がnormalizeしてるか確認した。
View chi2_test_normalize_check.cc
{
TH1D *h1 = new TH1D("h1", "h1", 3, 0, 3);
h1->Fill(1);
h1->Fill(1);
h1->Fill(2);
h1->Fill(2);
h1->Fill(2);
h1->Fill(2);
h1->Fill(3);
h1->Fill(3);
@pn11
pn11 / check_rms.cc
Created Sep 11, 2014
ROOTのTMath::RMS() (http://root.cern.ch/root/html/TMath.html#TMath:RMS) はRMSじゃなくて標準偏差を計算してるけど、それが標本標準偏差(?)なのか不偏標準偏差(?)なのか(標本数nで割ってるのかn-1で割ってるのか)調べた。結果、n-1で割ってた。
View check_rms.cc
{
const int ndata = 10000;
double data[ndata];
TRandom *gaus = new TRandom();
for (int idata = 0; idata < ndata; idata++){
data[idata] = gaus->Gaus(0., 1.);
}
@pn11
pn11 / draw_Kolmogorov.cc
Created Sep 16, 2014
TMath::KolmogorovProb() (http://root.cern.ch/root/html/src/TMath.cxx.html#RDBIQ) のソースが良くわからない(c1, c2, c3という謎の係数が出てくる)のでちゃんと正しい値なのか描いてみた件。
View draw_Kolmogorov.cc
{
Int_t npoint = 10000;
Double_t x, xstart = 0., xend = 2.0;
TGraph *g1 = new TGraph();
for (Int_t ipoint=0; ipoint<npoint; ipoint++){
x = xstart + (xend - xstart) / double(npoint) * double(ipoint);
g1->SetPoint(ipoint, x, TMath::KolmogorovProb(x));
}
g1->SetTitle("Kolmogorov distribution");
@pn11
pn11 / central_limit_theorem.cc
Last active Aug 29, 2015
中心極限定理を確かめて遊んでみた。(script for ROOT)
View central_limit_theorem.cc
// 中心極限定理を確かめて遊んでみた。(script for ROOT)
{
int n = 10000;
int m = 2;
int nbin = 100;
double limit = 8.0;
double binwidth = 2.0*limit/double(nbin);
TH1D *h = new TH1D("h", "hist", nbin, - limit, limit);
@pn11
pn11 / hadd_and_rename.sh
Created Sep 30, 2014
haddしてから足しあわせた数をファイル名に追加するだけのスクリプト。while read 以下のloopの中でexprを使ってnadd+1をやろうとしたのになぜかできなかった。しかたないのでリダイレクトしてごにょごにょ。
View hadd_and_rename.sh
#!/bin/bash
#hadd_and_rename.sh target.root sources.root
if [ $# -lt 2 ]
then
echo "Number of argument is wrong. Exit hadd_and_rename.sh..."
exit
fi
@pn11
pn11 / plot_Lindhard.cc
Created Nov 3, 2014
ROOT script for just calculating quenching factor of Si-28 by Lindhard eq. Lindhardの式でシリコンのQFを計算するだけのスクリプト。
View plot_Lindhard.cc
{
double z = 14.; // atomic number
double a = 28.; // mass number
int npoint = 1000;
// double range[2] = {0., 1000.}; // 0 - 1 MeV
// double range[2] = {0., 100.}; // 0 - 100 keV
double range[2] = {0., 10.}; // 0 - 10 keV
@pn11
pn11 / ConvertMollweide.cc
Last active Aug 29, 2015
緯度、経度からモルワイデ図法での座標を出すプログラムをクラスの勉強をしながらつくった。定義はWikipediaに従う。テータはニュートン法で求める。 http://ja.wikipedia.org/wiki/%E3%83%A2%E3%83%AB%E3%83%AF%E3%82%A4%E3%83%87%E5%9B%B3%E6%B3%95
View ConvertMollweide.cc
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <iostream>
#include <fstream>
using namespace std;
class ConvertMollweide
@pn11
pn11 / myAm.mac
Last active Aug 29, 2015
Am-241 source defined with GPS.
View myAm.mac
#/tracking/verbose 0
/tracking/verbose 1
#/control/verbose 1
# /vis/viewer/set/targetPoint 0. 0. 0. cm
# /vis/viewer/set/targetPoint 12.5 0. 0. cm
#/vis/viewer/set/viewpointVector 1 1 1
/vis/viewer/set/viewpointVector 0 0 1
# Those method are from
@pn11
pn11 / backup_android.md
Created Feb 24, 2015
Mac (Homebrew) でAndroidをバックアップする
View backup_android.md

#Mac (Homebrew) でAndroidをバックアップする

##adbを入れる brew install android-platform-tools

##接続確認 AndroidのUSBデバッグモードをONにしてつないで

adb devices