Script to display SOS hodoscope hit maps

examples/hitmaps_sos.C

+{
+  TFile* f = new TFile("hodtest.root");
+ 
+  TCanvas *c1 = new TCanvas("c1", "SOS Scintillator Hit Maps", 800, 800); 
+  c1->Divide(2, 4);
+
+  TH1F* h[16];
+
+  h[0] = spostdc1;
+  h[1] = snegtdc1;
+  h[4] = spostdc2;
+  h[5] = snegtdc2;
+  h[8] = spostdc3;
+  h[9] = snegtdc3;
+  h[12] = spostdc4;
+  h[13] = snegtdc4;
+
+  h[2] = sposadc1;
+  h[3] = snegadc1;
+  h[6] = sposadc2;
+  h[7] = snegadc2;
+  h[10] = sposadc3;
+  h[11] = snegadc3;
+  h[14] = sposadc4;
+  h[15] = snegadc4;
+  
+  for(Int_t ipad=1;ipad<=8;ipad++) {
+    TPad *pt = c1->cd(ipad);
+    Int_t ihp = 2*(ipad-1);
+    Int_t ihn = ihp+1;
+    h[ihp]->SetFillColor(kGreen);
+    h[ihp]->SetFillStyle(3345);
+    h[ihp]->Draw();
+    h[ihn]->SetFillColor(kBlue);
+    h[ihn]->SetFillStyle(3354);
+    h[ihn]->Draw("same");
+  }
+
+}

examples/output.def

@@ -177,6 +177,26 @@ 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
+# TDC hits per paddle
+TH1F spostdc1 'SOS s1x+ TDC hits' S.hod.1x.postdchits 9 0.5 9.5
+TH1F snegtdc1 'SOS s1x- TDC hits' S.hod.1x.negtdchits 9 0.5 9.5
+TH1F spostdc2 'SOS s1y+ TDC hits' S.hod.1y.postdchits 9 0.5 9.5
+TH1F snegtdc2 'SOS s1y- TDC hits' S.hod.1y.negtdchits 9 0.5 9.5
+TH1F spostdc3 'SOS s2x+ TDC hits' S.hod.2x.postdchits 16 0.5 16.5
+TH1F snegtdc3 'SOS s2x- TDC hits' S.hod.2x.negtdchits 16 0.5 16.5
+TH1F spostdc4 'SOS s2y+ TDC hits' S.hod.2y.postdchits 9 0.5 9.5
+TH1F snegtdc4 'SOS s2y- TDC hits' S.hod.2y.negtdchits 9 0.5 9.5
+
+# ADC hits per paddle
+TH1F sposadc1 'SOS s1x+ ADC hits' S.hod.1x.posadchits 9 0.5 9.5
+TH1F snegadc1 'SOS s1x- ADC hits' S.hod.1x.negadchits 9 0.5 9.5
+TH1F sposadc2 'SOS s1y+ ADC hits' S.hod.1y.posadchits 9 0.5 9.5
+TH1F snegadc2 'SOS s1y- ADC hits' S.hod.1y.negadchits 9 0.5 9.5
+TH1F sposadc3 'SOS s2x+ ADC hits' S.hod.2x.posadchits 16 0.5 16.5
+TH1F snegadc3 'SOS s2x- ADC hits' S.hod.2x.negadchits 16 0.5 16.5
+TH1F sposadc4 'SOS s2y+ ADC hits' S.hod.2y.posadchits 9 0.5 9.5
+TH1F snegadc4 'SOS s2y- ADC hits' S.hod.2y.negadchits 9 0.5 9.5
+
 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