From 5ebbb6532708e94a87f0e9366a4cc507dfea2441 Mon Sep 17 00:00:00 2001 From: "Kostyleva D.A" Date: Thu, 2 Feb 2017 20:05:11 +0300 Subject: [PATCH] Logbook info updategit add data/rawDataDSR4/LOGBOOK.md macros/histCorrNew.C --- data/rawDataDSR4/LOGBOOK.md | 30 ++++++++-------- macros/histCorrNew.C | 70 ++++++++++++++++++++++++++++++------- 2 files changed, 73 insertions(+), 27 deletions(-) diff --git a/data/rawDataDSR4/LOGBOOK.md b/data/rawDataDSR4/LOGBOOK.md index 0bfe215..4e38a9f 100644 --- a/data/rawDataDSR4/LOGBOOK.md +++ b/data/rawDataDSR4/LOGBOOK.md @@ -2,20 +2,22 @@ Raw data. Both in binary and root format. Logbook +PMT amplification: HV = -900 V, threshold ~ -10 mV + File name | Num of events | Comment -------------------------|---------------|---------------------------------------------------------------------| -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| 1039 | neutron source Cf | -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 | -------------------------|---------------|---------------------------------------------------------------------| +------------------------|---------------|--------------------------------------------------------------------------------------------------------------------| +NeuRad_test_071216_1.dat| 500 | collimated gamma source; ch0 and ch1 have coincidence; trigger ch0; 2 PMT H9500; periphery PMT pixels | +NeuRad_test_071216_2.dat| 5000 | collimated gamma source; ch0 and ch1 have coincidence; trigger ch0; 2 PMT H9500; periphery PMT pixels | +NeuRad_test_071216_3.dat| 1039 | non-collimated neutron source Cf; 2 PMT H9500; periphery PMT pixels | +NeuRad_test_071216_4.dat| 2000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back; internal PMT pixels | +NeuRad_test_071216_5.dat| 2000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back; internal PMT pixels | +NeuRad_test_071216_6.dat| 2000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back; internal PMT pixels | +NeuRad_test_071216_7.dat| 2000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back; internal PMT pixels | +------------------------|---------------|--------------------------------------------------------------------------------------------------------------------| -Change of PMT amplification: increasing voltage up to 1100 mV, threshold 25 mV +Change of PMT amplification: increasing HV up to -1100 V, threshold -25 mV -------------------------|---------------|---------------------------------------------------------------------| -NeuRad_test_071216_8.dat| 10000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back | -NeuRad_test_081216_1.dat| 2204 | trigger ch2 ????? | -NeuRad_test_081216_2.dat| 10000 | collimated gamma source; two channels only: ch0 and ch1 | +------------------------|---------------|--------------------------------------------------------------------------------------------------------------------| +NeuRad_test_071216_8.dat| 10000 | source 60Co; trigger ch0; ch0-front, ch1,ch2,ch3 - back; internal PMT pixels !!!! | +NeuRad_test_081216_1.dat| 2204 | trigger ch2; internal PMT pixels ????? | +NeuRad_test_081216_2.dat| 10000 | collimated gamma source; two channels only: ch0 and ch1; changed PMT H9500 -ch0, H10966 - ch1 ; internal PMT pixels| diff --git a/macros/histCorrNew.C b/macros/histCorrNew.C index a1ead1c..eea6209 100644 --- a/macros/histCorrNew.C +++ b/macros/histCorrNew.C @@ -4,7 +4,7 @@ void histCorrNew ( - const char *filename = "analysis_08_2.root", + const char *filename = "analysis_07_8.root", // const char *foldername = "report", const char *foldername = "8_2", // const char *ext = ".eps", @@ -64,7 +64,7 @@ void histCorrNew ( c1->Update(); - c1->Print(Form("../macros/picsDRS4/file%s/time_front_corr1%s", foldername, ext)); + // c1->Print(Form("../macros/picsDRS4/file%s/time_front_corr1%s", foldername, ext)); } if(can2) { @@ -103,7 +103,7 @@ void histCorrNew ( h6->GetYaxis()->SetTitle("Max_ampl_{0} [V]"); h6->GetYaxis()->CenterTitle(); - c2->Print(Form("../macros/picsDRS4/file%s/time_front_corr2%s", foldername, ext)); + // c2->Print(Form("../macros/picsDRS4/file%s/time_front_corr2%s", foldername, ext)); } if(can3) { @@ -111,22 +111,66 @@ void histCorrNew ( c3->Divide(2,2); c3->cd(1); - TH1F *h7 = new TH1F("h7", "tau TF", 100, -5 , 25); - t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h7","Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150 && Ach0.fEdgeSlope > 0.01 && Ach1.fEdgeSlope > 0.01",""); - h7->Fit("gaus", "","",7,14); + TH1F *h10 = new TH1F("h10", "delta t (front edge fit)", 100, -5 , 25); + t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h10","Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150",""); + h10->Fit("gaus", "","",-5,5); gStyle->SetOptFit(); + h10->GetXaxis()->SetTitle("delta t [ns]"); + h10->GetXaxis()->CenterTitle(); c3->cd(2); - TH1F *h8 = new TH1F("h8", "tau TF", 100, -5 , 25); - t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h8","Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150 && Ach0.fChargeTF > 0.05e-9 && Ach1.fChargeTF > 0.025e-9",""); - h8->Fit("gaus","","",7,14); + TH1F *h7 = new TH1F("h7", "delta t (front edge fit) with edge slope correction", 100, -5 , 25); + // t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h7", + // "Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150 && Ach0.fEdgeSlope > 0.01 && Ach1.fEdgeSlope > 0.01",""); + t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h7","Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150 && Ach0.fEdgeSlope > 0.005 && Ach1.fEdgeSlope > 0.001",""); + h7->Fit("gaus", "","",-5,5); gStyle->SetOptFit(); + h7->GetXaxis()->SetTitle("delta t [ns]"); + h7->GetXaxis()->CenterTitle(); + + c3->cd(3); + TH1F *h8 = new TH1F("h8", "delta t (front edge fit) with charge correction", 100, -5 , 25); + t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h8", + "Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150 && Ach0.fChargeTF > 0.05e-9 && Ach1.fChargeTF > 0.025e-9",""); + h8->Fit("gaus","","",-5,5); + gStyle->SetOptFit(); + h8->GetXaxis()->SetTitle("delta t [ns]"); + h8->GetXaxis()->CenterTitle(); + + c3->cd(4); + TH1F *h9 = new TH1F("h9", "delta t (front edge fit) with rise time correction", 100, -5 , 25); + t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h9", + "Ach0.fTimeFront> 100 && Ach0.fTimeFront<150 && Ach1.fTimeFront> 100 && Ach1.fTimeFront<150 && (Ach0.fTime90 - Ach0.fTime10) > 1.5 && (Ach0.fTime90 - Ach0.fTime10)<3 && (Ach1.fTime90 - Ach1.fTime10) > 1.5 && (Ach1.fTime90 - Ach1.fTime10)<3",""); + h9->Fit("gaus","","",-5,5); + gStyle->SetOptFit(); + h9->GetXaxis()->SetTitle("delta t [ns]"); + h9->GetXaxis()->CenterTitle(); - c3->Print(Form("../macros/picsDRS4/file%s/time_front_corr3%s", foldername, ext)); - - + // c3->Print(Form("../macros/picsDRS4/file%s/time_front_corr3%s", foldername, ext)); } + if(can4) { + TCanvas *c4 = new TCanvas ("c4", "Correlation pictures for collimated data: time difference"); + c4->Divide(2,2); + + c4->cd(1); + TH1F *h11 = new TH1F("h11", "rise time", 100, 0 , 5); + t->Draw("Ach0.fTime90 - Ach0.fTime10 >> h11","",""); + //gStyle->SetLineColor(kRed); + + c4->cd(2); + TH1F *h12 = new TH1F("h12", "rise time", 100, 0 , 5); + t->Draw("Ach1.fTime90 - Ach1.fTime10 >> h12","",""); + // h10->Fit("gaus", "","",7,14); + // gStyle->SetOptFit(); + // h10->GetXaxis()->SetTitle("delta t [ns]"); + // h10->GetXaxis()->CenterTitle(); + + c4->cd(3); + TH1F *h13 = new TH1F("h13", "rise time", 100, -10, 10); + t->Draw("( Ach0.fTime10 + (Ach0.fTime90 - Ach0.fTime10)*0.5 ) - ( Ach1.fTime10 + (Ach1.fTime90 - Ach1.fTime10)*0.5 )>> h13","Ach0.fEdgeSlope > 0.03 && Ach1.fEdgeSlope > 0.006",""); + h13->Fit("gaus","","",-5,5); + gStyle->SetOptFit(); - + } } -- 2.18.1