From 49050cf22ab59fbc4ca2968a5410a8c0b9fc65e2 Mon Sep 17 00:00:00 2001
From: "Stephen A. Wood" <saw@jlab.org>
Date: Thu, 12 Dec 2013 16:05:12 -0500
Subject: [PATCH] Add some SOS DC histograms and scripts
---
examples/{dcdriftdist.C => dcdriftdist_hms.C} | 0
examples/dcdriftdist_sos.C | 31 +++++++++++++++++
examples/{dchitmaps.C => dchitmaps_hms.C} | 0
examples/dchitmaps_sos.C | 31 +++++++++++++++++
examples/{dcresiduals.C => dcresiduals_hms.C} | 1 +
examples/dcresiduals_sos.C | 34 +++++++++++++++++++
examples/output.def | 27 +++++++++++++++
7 files changed, 124 insertions(+)
rename examples/{dcdriftdist.C => dcdriftdist_hms.C} (100%)
create mode 100644 examples/dcdriftdist_sos.C
rename examples/{dchitmaps.C => dchitmaps_hms.C} (100%)
create mode 100644 examples/dchitmaps_sos.C
rename examples/{dcresiduals.C => dcresiduals_hms.C} (98%)
create mode 100644 examples/dcresiduals_sos.C
diff --git a/examples/dcdriftdist.C b/examples/dcdriftdist_hms.C
similarity index 100%
rename from examples/dcdriftdist.C
rename to examples/dcdriftdist_hms.C
diff --git a/examples/dcdriftdist_sos.C b/examples/dcdriftdist_sos.C
new file mode 100644
index 0000000..935be5a
--- /dev/null
+++ b/examples/dcdriftdist_sos.C
@@ -0,0 +1,31 @@
+ {
+ TFile* f = new TFile("hodtest.root");
+
+ TCanvas *c1 = new TCanvas("c1", "Drift Chamber Drift Distance", 800, 800);
+ c1->Divide(2, 6);
+
+ TH1F* h[12];
+
+ h[0] = sdc1u1_dd;
+ h[1] = sdc1u2_dd;
+ h[2] = sdc1x1_dd;
+ h[3] = sdc1x2_dd;
+ h[4] = sdc1v1_dd;
+ h[5] = sdc1v2_dd;
+ h[6] = sdc2u1_dd;
+ h[7] = sdc2u2_dd;
+ h[8] = sdc2x1_dd;
+ h[9] = sdc2x2_dd;
+ h[10] = sdc2v1_dd;
+ h[11] = sdc2v2_dd;
+
+ for(Int_t ih=0;ih<12;ih++) {
+ Int_t ipad = 2*(ih%6)+1 + ih/6;
+ cout << ipad << endl;
+ c1->cd(ipad);
+ h[ih]->SetFillColor(kBlue);
+ h[ih]->SetFillStyle(3354);
+ h[ih]->Draw();
+ }
+
+}
diff --git a/examples/dchitmaps.C b/examples/dchitmaps_hms.C
similarity index 100%
rename from examples/dchitmaps.C
rename to examples/dchitmaps_hms.C
diff --git a/examples/dchitmaps_sos.C b/examples/dchitmaps_sos.C
new file mode 100644
index 0000000..ec7a1d0
--- /dev/null
+++ b/examples/dchitmaps_sos.C
@@ -0,0 +1,31 @@
+ {
+ TFile* f = new TFile("hodtest.root");
+
+ TCanvas *c1 = new TCanvas("c1", "Drift Chamber Hit Maps", 800, 800);
+ c1->Divide(2, 6);
+
+ TH1F* h[12];
+
+ h[0] = sdc1u1_wm;
+ h[1] = sdc1u2_wm;
+ h[2] = sdc1x1_wm;
+ h[3] = sdc1x2_wm;
+ h[4] = sdc1v1_wm;
+ h[5] = sdc1v2_wm;
+ h[6] = sdc2u1_wm;
+ h[7] = sdc2u2_wm;
+ h[8] = sdc2x1_wm;
+ h[9] = sdc2x2_wm;
+ h[10] = sdc2v1_wm;
+ h[11] = sdc2v2_wm;
+
+ for(Int_t ih=0;ih<12;ih++) {
+ Int_t ipad = 2*(ih%6)+1 + ih/6;
+ cout << ipad << endl;
+ c1->cd(ipad);
+ h[ih]->SetFillColor(kBlue);
+ h[ih]->SetFillStyle(3354);
+ h[ih]->Draw();
+ }
+
+}
diff --git a/examples/dcresiduals.C b/examples/dcresiduals_hms.C
similarity index 98%
rename from examples/dcresiduals.C
rename to examples/dcresiduals_hms.C
index a5759b9..ae86020 100644
--- a/examples/dcresiduals.C
+++ b/examples/dcresiduals_hms.C
@@ -30,4 +30,5 @@
T->Draw("H.dc.residual[10]","H.dc.residual[10] < 999");
c1->cd(12);
T->Draw("H.dc.residual[11]","H.dc.residual[11] < 999");
+ delete h;
}
diff --git a/examples/dcresiduals_sos.C b/examples/dcresiduals_sos.C
new file mode 100644
index 0000000..cf2dd5a
--- /dev/null
+++ b/examples/dcresiduals_sos.C
@@ -0,0 +1,34 @@
+{
+ TFile* f = new TFile("hodtest.root");
+ TTree *T=(TTree*)f->Get("T");
+ TCanvas *c1 = new TCanvas("c1", "Drift Chamber Residuals", 800, 1000);
+ c1->Divide(2,6);
+
+ TH1F* h=new TH1F("h", "Residual", 100, -1.0, 1.0);
+
+ c1->cd(1);
+ T->Draw("S.dc.residual[0]","S.dc.residual[0] < 999");
+ c1->cd(2);
+ T->Draw("S.dc.residual[1]","S.dc.residual[1] < 999");
+ c1->cd(3);
+ T->Draw("S.dc.residual[2]","S.dc.residual[2] < 999");
+ c1->cd(4);
+ T->Draw("S.dc.residual[3]","S.dc.residual[3] < 999");
+ c1->cd(5);
+ T->Draw("S.dc.residual[4]","S.dc.residual[4] < 999");
+ c1->cd(6);
+ T->Draw("S.dc.residual[5]","S.dc.residual[5] < 999");
+ c1->cd(7);
+ T->Draw("S.dc.residual[6]","S.dc.residual[6] < 999");
+ c1->cd(8);
+ T->Draw("S.dc.residual[7]","S.dc.residual[7] < 999");
+ c1->cd(9);
+ T->Draw("S.dc.residual[8]","S.dc.residual[8] < 999");
+ c1->cd(10);
+ T->Draw("S.dc.residual[9]","S.dc.residual[9] < 999");
+ c1->cd(11);
+ T->Draw("S.dc.residual[10]","S.dc.residual[10] < 999");
+ c1->cd(12);
+ T->Draw("S.dc.residual[11]","S.dc.residual[11] < 999");
+ delete h;
+}
diff --git a/examples/output.def b/examples/output.def
index da415c8..1f5853b 100644
--- a/examples/output.def
+++ b/examples/output.def
@@ -175,3 +175,30 @@ TH1F hs1yfptime 'HODO s1y fptime' H.hod.1y.fptime 80 0 80 H.hod.hgoodstarttime
TH1F hs2xfptime 'HODO s2x fptime' H.hod.2x.fptime 80 0 80 H.hod.hgoodstarttime
TH1F hs2yfptime 'HODO s2y fptime' H.hod.2y.fptime 80 0 80 H.hod.hgoodstarttime
TH1F starttime 'HODO start time' H.hod.starttime 80 0 80 H.hod.hgoodstarttime
+
+#SOS
+TH1F sdc1u1_wm 'SDC 1U1 Wiremap' S.dc.1u1.tdchits 48 0.5 48.5
+TH1F sdc1u2_wm 'SDC 1U2 Wiremap' S.dc.1u2.tdchits 48 0.5 48.5
+TH1F sdc1x1_wm 'SDC 1X1 Wiremap' S.dc.1x1.tdchits 64 0.5 64.5
+TH1F sdc1x2_wm 'SDC 1X2 Wiremap' S.dc.1x2.tdchits 64 0.5 64.5
+TH1F sdc1v1_wm 'SDC 1V1 Wiremap' S.dc.1v1.tdchits 48 0.5 48.5
+TH1F sdc1v2_wm 'SDC 1V2 Wiremap' S.dc.1v2.tdchits 48 0.5 48.5
+TH1F sdc2u1_wm 'SDC 2U1 Wiremap' S.dc.2u1.tdchits 48 0.5 48.5
+TH1F sdc2u2_wm 'SDC 2U2 Wiremap' S.dc.2u2.tdchits 48 0.5 48.5
+TH1F sdc2x1_wm 'SDC 2X1 Wiremap' S.dc.2x1.tdchits 64 0.5 64.5
+TH1F sdc2x2_wm 'SDC 2X2 Wiremap' S.dc.2x2.tdchits 64 0.5 64.5
+TH1F sdc2v1_wm 'SDC 2V1 Wiremap' S.dc.2v1.tdchits 48 0.5 48.5
+TH1F sdc2v2_wm 'SDC 2V2 Wiremap' S.dc.2v2.tdchits 48 0.5 48.5
+
+TH1F sdc1u1_dd 'SDC 1U1 Drift Distance' S.dc.1u1.dist 300 -0.1 0.6
+TH1F sdc1u2_dd 'SDC 1U2 Drift Distance' S.dc.1u2.dist 300 -0.1 0.6
+TH1F sdc1x1_dd 'SDC 1X1 Drift Distance' S.dc.1x1.dist 300 -0.1 0.6
+TH1F sdc1x2_dd 'SDC 1X2 Drift Distance' S.dc.1x2.dist 300 -0.1 0.6
+TH1F sdc1v1_dd 'SDC 1V1 Drift Distance' S.dc.1v1.dist 300 -0.1 0.6
+TH1F sdc1v2_dd 'SDC 1V2 Drift Distance' S.dc.1v2.dist 300 -0.1 0.6
+TH1F sdc2u1_dd 'SDC 2U1 Drift Distance' S.dc.2u1.dist 300 -0.1 0.6
+TH1F sdc2u2_dd 'SDC 2U2 Drift Distance' S.dc.2u2.dist 300 -0.1 0.6
+TH1F sdc2x1_dd 'SDC 2X1 Drift Distance' S.dc.2x1.dist 300 -0.1 0.6
+TH1F sdc2x2_dd 'SDC 2X2 Drift Distance' S.dc.2x2.dist 300 -0.1 0.6
+TH1F sdc2v1_dd 'SDC 2V1 Drift Distance' S.dc.2v1.dist 300 -0.1 0.6
+TH1F sdc2v2_dd 'SDC 2V2 Drift Distance' S.dc.2v2.dist 300 -0.1 0.6
--
GitLab