From 97af94f2f9a26705dbe5a295d1bc4e5f80ff7ba7 Mon Sep 17 00:00:00 2001
From: "Stephen A. Wood" <saw@jlab.org>
Date: Tue, 19 Mar 2013 11:35:54 -0400
Subject: [PATCH] Add histograms for drift time and drift distance.

---
 examples/dcdriftdist.C | 31 +++++++++++++++++++++++++++++++
 examples/output.def    | 28 ++++++++++++++++++++++++++++
 2 files changed, 59 insertions(+)
 create mode 100644 examples/dcdriftdist.C

diff --git a/examples/dcdriftdist.C b/examples/dcdriftdist.C
new file mode 100644
index 0000000..fe67618
--- /dev/null
+++ b/examples/dcdriftdist.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] = hdc1x1_dd;
+  h[1] = hdc1y1_dd;
+  h[2] = hdc1u1_dd;
+  h[3] = hdc1v1_dd;
+  h[4] = hdc1y2_dd;
+  h[5] = hdc1x2_dd;
+  h[6] = hdc2x1_dd;
+  h[7] = hdc2y1_dd;
+  h[8] = hdc2u1_dd;
+  h[9] = hdc2v1_dd;
+  h[10] = hdc2y2_dd;
+  h[11] = hdc2x2_dd;
+  
+  for(Int_t ih=0;ih<12;ih++) {
+    Int_t ipad = 2*(ih%6)+1 + ih/6;
+    cout << ipad << endl;
+    TPad *pt = c1->cd(ipad);
+    h[ih]->SetFillColor(kBlue);
+    h[ih]->SetFillStyle(3354);
+    h[ih]->Draw();
+  }
+
+}
diff --git a/examples/output.def b/examples/output.def
index 7cfb0da..8369d6c 100644
--- a/examples/output.def
+++ b/examples/output.def
@@ -65,3 +65,31 @@ TH1F hdc2u1_wm 'HDC 2U1 Wiremap' H.dc.2u1.tdchits 107 0.5 107.5
 TH1F hdc2v1_wm 'HDC 2V1 Wiremap' H.dc.2v1.tdchits 107 0.5 107.5
 TH1F hdc2y2_wm 'HDC 2Y2 Wiremap' H.dc.2y2.tdchits 52 0.5 52.5
 TH1F hdc2x2_wm 'HDC 2X2 Wiremap' H.dc.2x2.tdchits 113 0.5 113.5
+
+# Drift Time and Drift Distance Spectra
+TH1F hdc1x1_dt 'HDC 1X1 Drift Time' H.dc.1x1.time 200 -100 300
+TH1F hdc1y1_dt 'HDC 1Y1 Drift Time' H.dc.1y1.time 200 -100 300
+TH1F hdc1u1_dt 'HDC 1U1 Drift Time' H.dc.1u1.time 200 -100 300
+TH1F hdc1v1_dt 'HDC 1V1 Drift Time' H.dc.1v1.time 200 -100 300
+TH1F hdc1y2_dt 'HDC 1Y2 Drift Time' H.dc.1y2.time 200 -100 300
+TH1F hdc1x2_dt 'HDC 1X2 Drift Time' H.dc.1x2.time 200 -100 300
+TH1F hdc2x1_dt 'HDC 2X1 Drift Time' H.dc.2x1.time 200 -100 300
+TH1F hdc2y1_dt 'HDC 2Y1 Drift Time' H.dc.2y1.time 200 -100 300
+TH1F hdc2u1_dt 'HDC 2U1 Drift Time' H.dc.2u1.time 200 -100 300
+TH1F hdc2v1_dt 'HDC 2V1 Drift Time' H.dc.2v1.time 200 -100 300
+TH1F hdc2y2_dt 'HDC 2Y2 Drift Time' H.dc.2y2.time 200 -100 300
+TH1F hdc2x2_dt 'HDC 2X2 Drift Time' H.dc.2x2.time 200 -100 300
+
+TH1F hdc1x1_dd 'HDC 1X1 Drift Distance' H.dc.1x1.dist 300 -0.1 0.6
+TH1F hdc1y1_dd 'HDC 1Y1 Drift Distance' H.dc.1y1.dist 300 -0.1 0.6
+TH1F hdc1u1_dd 'HDC 1U1 Drift Distance' H.dc.1u1.dist 300 -0.1 0.6
+TH1F hdc1v1_dd 'HDC 1V1 Drift Distance' H.dc.1v1.dist 300 -0.1 0.6
+TH1F hdc1y2_dd 'HDC 1Y2 Drift Distance' H.dc.1y2.dist 300 -0.1 0.6
+TH1F hdc1x2_dd 'HDC 1X2 Drift Distance' H.dc.1x2.dist 300 -0.1 0.6
+TH1F hdc2x1_dd 'HDC 2X1 Drift Distance' H.dc.2x1.dist 300 -0.1 0.6
+TH1F hdc2y1_dd 'HDC 2Y1 Drift Distance' H.dc.2y1.dist 300 -0.1 0.6
+TH1F hdc2u1_dd 'HDC 2U1 Drift Distance' H.dc.2u1.dist 300 -0.1 0.6
+TH1F hdc2v1_dd 'HDC 2V1 Drift Distance' H.dc.2v1.dist 300 -0.1 0.6
+TH1F hdc2y2_dd 'HDC 2Y2 Drift Distance' H.dc.2y2.dist 300 -0.1 0.6
+TH1F hdc2x2_dd 'HDC 2X2 Drift Distance' H.dc.2x2.dist 300 -0.1 0.6
+
-- 
GitLab