Skip to content

N
NeuRad_tests
Commit 5ebbb653 authored by Kostyleva D.A's avatar Kostyleva D.A
Browse files
Options

Logbook info updategit add data/rawDataDSR4/LOGBOOK.md macros/histCorrNew.C

parent bd45314e
Show whitespace changes
Inline Side-by-side
Showing with 73 additions and 27 deletions
data/rawDataDSR4/LOGBOOK.md
View file @ 5ebbb653
...@@ -2,20 +2,22 @@ Raw data. Both in binary and root format. ...@@ -2,20 +2,22 @@ Raw data. Both in binary and root format.
Logbook Logbook
PMT amplification: HV = -900 V, threshold ~ -10 mV
File name | Num of events | Comment 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_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 | 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 | neutron source Cf | 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 | 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 | 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 | 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 | 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_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 ????? | 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 | NeuRad_test_081216_2.dat| 10000 | collimated gamma source; two channels only: ch0 and ch1; changed PMT H9500 -ch0, H10966 - ch1 ; internal PMT pixels|
macros/histCorrNew.C
View file @ 5ebbb653
...@@ -4,7 +4,7 @@ ...@@ -4,7 +4,7 @@
void histCorrNew ( void histCorrNew (
const char *filename = "analysis_08_2.root", const char *filename = "analysis_07_8.root",
// const char *foldername = "report", // const char *foldername = "report",
const char *foldername = "8_2", const char *foldername = "8_2",
// const char *ext = ".eps", // const char *ext = ".eps",
...@@ -64,7 +64,7 @@ void histCorrNew ( ...@@ -64,7 +64,7 @@ void histCorrNew (
c1->Update(); 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) { if(can2) {
...@@ -103,7 +103,7 @@ void histCorrNew ( ...@@ -103,7 +103,7 @@ void histCorrNew (
h6->GetYaxis()->SetTitle("Max_ampl_{0} [V]"); h6->GetYaxis()->SetTitle("Max_ampl_{0} [V]");
h6->GetYaxis()->CenterTitle(); 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) { if(can3) {
...@@ -111,22 +111,66 @@ void histCorrNew ( ...@@ -111,22 +111,66 @@ void histCorrNew (
c3->Divide(2,2); c3->Divide(2,2);
c3->cd(1); c3->cd(1);
TH1F *h7 = new TH1F("h7", "tau TF", 100, -5 , 25); TH1F *h10 = new TH1F("h10", "delta t (front edge fit)", 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 >> h10","Ach0.fTimeFront>100. && Ach0.fTimeFront < 150 && Ach1.fTimeFront>100 && Ach1.fTimeFront < 150","");
h7->Fit("gaus", "","",7,14); h10->Fit("gaus", "","",-5,5);
gStyle->SetOptFit(); gStyle->SetOptFit();
h10->GetXaxis()->SetTitle("delta t [ns]");
h10->GetXaxis()->CenterTitle();
c3->cd(2); c3->cd(2);
TH1F *h8 = new TH1F("h8", "tau TF", 100, -5 , 25); TH1F *h7 = new TH1F("h7", "delta t (front edge fit) with edge slope 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",""); // t->Draw("Ach0.fTimeFront - Ach1.fTimeFront >> h7",
h8->Fit("gaus","","",7,14); // "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(); 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();
}
} }
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment