From bae4d8e64c8f3447ff03b0b96ce3b4ecc132ca86 Mon Sep 17 00:00:00 2001 From: "Kostyleva D.A" Date: Fri, 13 Jan 2017 15:30:31 +0300 Subject: [PATCH] Scripts for drawing CFD histos added --- data/rawDataDSR4/README.md | 18 +++++++++--------- macros/analyse.C | 4 ++-- macros/histCFDf7_8.C | 16 ++++++++++++++++ macros/histCFDf8_2.C | 16 ++++++++++++++++ 4 files changed, 43 insertions(+), 11 deletions(-) create mode 100644 macros/histCFDf7_8.C create mode 100644 macros/histCFDf8_2.C diff --git a/data/rawDataDSR4/README.md b/data/rawDataDSR4/README.md index a0eea58..d4ac709 100644 --- a/data/rawDataDSR4/README.md +++ b/data/rawDataDSR4/README.md @@ -4,18 +4,18 @@ Logbook File name | Num of events | Comment ------------------------|---------------|---------------------------------------------------------------------| -NeuRad_test_071216_1.dat| 500 | -NeuRad_test_071216_2.dat| 500 | collimated gamma source; ch0 and ch1 have coincidence; ?trigger ch0?| +NeuRad_test_071216_1.dat| 500 | collimated gamma source; ch0 and ch1 have coincidence; trigger ch0 | +NeuRad_test_071216_2.dat| 5000 | collimated gamma source; ch0 and ch1 have coincidence; trigger ch0 | NeuRad_test_071216_3.dat| | neutron source Cf | -NeuRad_test_071216_4.dat| 2000 | source 60Co; trigger ch1; ch1-front, ch2,ch3,ch4 - back | -NeuRad_test_071216_5.dat| 2000 | source 60Co; trigger ch1; ch1-front, ch2,ch3,ch4 - back | -NeuRad_test_071216_6.dat| 2000 | source 60Co; trigger ch1; ch1-front, ch2,ch3,ch4 - back | -NeuRad_test_071216_7.dat| 2000 | source 60Co; trigger ch1; ch1-front, ch2,ch3,ch4 - back | +NeuRad_test_071216_4.dat| 2000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back | +NeuRad_test_071216_5.dat| 2000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back | +NeuRad_test_071216_6.dat| 2000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back | +NeuRad_test_071216_7.dat| 2000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back | ------------------------|---------------|---------------------------------------------------------------------| Change of PMT amplification: increasing voltage up to 1100 mV, threshold 25 mV ------------------------|---------------|---------------------------------------------------------------------| -NeuRad_test_071216_8.dat| 10000 | source 60Co; trigger ch1; ch1-front, ch2,ch3,ch4 - back ????? | -NeuRad_test_081216_1.dat| | trigger ch3 ????? | -NeuRad_test_081216_2.dat| | collimated gamma source; two channels only: ch1 and ch2 | +NeuRad_test_071216_8.dat| 10000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back | +NeuRad_test_081216_1.dat| | trigger ch2 ????? | +NeuRad_test_081216_2.dat| 10000 | collimated gamma source; two channels only: ch0 and ch1 | diff --git a/macros/analyse.C b/macros/analyse.C index 4921699..02d13c6 100644 --- a/macros/analyse.C +++ b/macros/analyse.C @@ -2,7 +2,7 @@ void analyse() { gSystem->Load("../libData.so"); - TFile *f = new TFile("../data/rawDataDSR4/NeuRad_test_07_8.root"); + TFile *f = new TFile("../data/rawDataDSR4/NeuRad_test_08_2.root"); TTree *tr = (TTree*)f->Get("rtree"); const Int_t noBranches = 4; @@ -19,7 +19,7 @@ void analyse() // tr->SetMakeClass(1); - TFile *fw = new TFile("../data/dataDSR4/analysis_07_8.root", "RECREATE"); //create .root file with somehow analyzed data + TFile *fw = new TFile("../data/dataDSR4/analysis_08_2.root", "RECREATE"); //create .root file with somehow analyzed data TTree *tw = new TTree("atree", "title of drs4 analysis tree"); //create analysis tree atree in it AEvent *wevent[noBranches]; // pointer to the array (of AEvent class) in which analyzed data for each channel will be put diff --git a/macros/histCFDf7_8.C b/macros/histCFDf7_8.C new file mode 100644 index 0000000..f4bf1e3 --- /dev/null +++ b/macros/histCFDf7_8.C @@ -0,0 +1,16 @@ +void histCFDf7_8() +{ + gSystem->Load("../libData.so"); + TFile *f = new TFile("../data/dataDSR4/analysis_07_8.root"); + TTree *t = (TTree*)f->Get("atree"); + TCanvas *c = new TCanvas("CFD canvas"); + c->Divide(2,2); + c->cd(1); + t->Draw("Ach0.fTimeCFD","",""); + c->cd(2); + t->Draw("Ach2.fTimeCFD","",""); + c->cd(3); + t->Draw("Ach0.fTimeCFD - Ach2.fTimeCFD","",""); + c->cd(4); + +} diff --git a/macros/histCFDf8_2.C b/macros/histCFDf8_2.C new file mode 100644 index 0000000..f10cfdb --- /dev/null +++ b/macros/histCFDf8_2.C @@ -0,0 +1,16 @@ +void histCFDf8_2() +{ + gSystem->Load("../libData.so"); + TFile *f = new TFile("../data/dataDSR4/analysis_08_2.root"); + TTree *t = (TTree*)f->Get("atree"); + TCanvas *c = new TCanvas("CFD canvas"); + c->Divide(2,2); + c->cd(1); + t->Draw("Ach0.fTimeCFD","",""); + c->cd(2); + t->Draw("Ach1.fTimeCFD","",""); + c->cd(3); + t->Draw("Ach0.fTimeCFD - Ach1.fTimeCFD","",""); + c->cd(4); + +} -- 2.18.1