geometry for shms (#217)

* added files used for shms dc alignment

* moved some dc alignment files around

* removed script in production

* updated shms geo param

* updated plots for shms dc monitoring

PARAM/SHMS/DC/GEOM/pdc_geom.param

@@ -26,22 +26,24 @@ pdc_sigma = 0.02, 0.02, 0.02, 0.02, 0.02, 0.02  ; TODO - SOS values for now
             0.02, 0.02, 0.02, 0.02, 0.02, 0.02
 
 ; Wire velocity corrections.
-pdc_wire_velocity  = 12.0;
-pdc_central_time   = 0, 0, 0, 0, 0, 0  ; TODO - no correction for now
+pdc_wire_velocity  = 12.0; -this should be checked (old HMS value for now)
+pdc_central_time   = 0, 0, 0, 0, 0, 0  ; -not used in new SHMS chambers
                      0, 0, 0, 0, 0, 0
-pdc_drifttime_sign = 0, 0, 0, 0, 0, 0  ; TODO - no correction for now
+pdc_drifttime_sign = 0, 0, 0, 0, 0, 0  ; -not used in new SHMS chambers
                      0, 0, 0, 0, 0, 0
-;pdc_readoutLR = -1.0, 1.0, -1.0, -1.0, 1.0, -1.0
-;                 1.0, -1.0, 1.0, 1.0, -1.0, 1.0
-;pdc_readoutTB = 1.0, -1.0, 1.0, -1.0, 1.0, -1.0
-;                -1.0, 1.0, -1.0, 1.0, -1.0, 1.0             
+
+pdc_readoutLR = -1.0, 1.0, -1.0, -1.0, 1.0, -1.0
+                 1.0, -1.0, 1.0, 1.0, -1.0, 1.0
+pdc_readoutTB = 1.0, -1.0, 1.0, -1.0, 1.0, -1.0
+                -1.0, 1.0, -1.0, 1.0, -1.0, 1.0             
 
 ; CHAMBER POSITION PARAMETERS
 
 ; x, y and z position of center for each wire chamber in cm.
 ; From survey C1771.
 pdc_xcenter = -0.036, 0.008 ;
-pdc_ycenter =  0.019, 0.013 ;
+pdc_ycenter = 0.019, 0.013 ;
+
 pdc_1_zpos = -40 - 0.656 ;
 pdc_2_zpos =  40 - 0.668 ;
 
@@ -99,7 +101,6 @@ pdc_ypos = pdc_ycenter[0]-0.007
            pdc_ycenter[1]+0.0035
            pdc_ycenter[1]-0.006
 
-
 ; z position in cm of each plane measured from focal plane.
 ; From technical drawings.
 pdc_zpos = pdc_1_zpos - 0.68701*cminch

onlineGUI/Config/SHMS/DC/pdc_stand.cfg

-protorootfile ../ROOTfiles/pdc_replay_XXXXX.root
+protorootfile ../ROOTfiles/shms_replay_production_XXXXX_1000000.root
 guicolor orange
-canvassize 1600 800
-#canvassize 800 800
-newpage 6 2
-title SHMS DC Wire Map
-pdc1u1_wirenum -nostat
-pdc1u2_wirenum -nostat
-pdc1x1_wirenum -nostat
-pdc1x2_wirenum -nostat
-pdc1v1_wirenum -nostat
-pdc1v2_wirenum -nostat
-pdc2v2_wirenum -nostat
-pdc2v1_wirenum -nostat
-pdc2x2_wirenum -nostat
-pdc2x1_wirenum -nostat
-pdc2u2_wirenum -nostat
-pdc2u1_wirenum -nostat
-newpage 6 2 
-title SHMS DC Plane Raw TDC 
-pdc1u1_raw_tdc -nostat
-pdc1u2_raw_tdc -nostat
-pdc1x1_raw_tdc -nostat
-pdc1x2_raw_tdc -nostat
-pdc1v1_raw_tdc -nostat
-pdc1v2_raw_tdc -nostat
-pdc2v2_raw_tdc -nostat
-pdc2v1_raw_tdc -nostat
-pdc2x2_raw_tdc -nostat
-pdc2x1_raw_tdc -nostat
-pdc2u2_raw_tdc -nostat
-pdc2u1_raw_tdc -nostat
+canvassize 1200 1000
 
-newpage 6 2 logz
-title SHMS DC Raw TDC vs. Wire Number
-pdc1u1_raw_tdc_vs_wirenum -nostat
-pdc1u2_raw_tdc_vs_wirenum -nostat
-pdc1x1_raw_tdc_vs_wirenum -nostat
-pdc1x2_raw_tdc_vs_wirenum -nostat
-pdc1v1_raw_tdc_vs_wirenum -nostat
-pdc1v2_raw_tdc_vs_wirenum -nostat
-pdc2v2_raw_tdc_vs_wirenum -nostat
-pdc2v1_raw_tdc_vs_wirenum -nostat
-pdc2x2_raw_tdc_vs_wirenum -nostat
-pdc2x1_raw_tdc_vs_wirenum -nostat
-pdc2u2_raw_tdc_vs_wirenum -nostat
-pdc2u1_raw_tdc_vs_wirenum -nostat
-newpage 6 2 
-title SHMS DC Plane Drift Time 
-pdc1u1_drifttime -nostat
-pdc1u2_drifttime -nostat
-pdc1x1_drifttime -nostat
-pdc1x2_drifttime -nostat
-pdc1v1_drifttime -nostat
-pdc1v2_drifttime -nostat
-pdc2v2_drifttime -nostat
-pdc2v1_drifttime -nostat
-pdc2x2_drifttime -nostat
-pdc2x1_drifttime -nostat
-pdc2u2_drifttime -nostat
-pdc2u1_drifttime -nostat
-newpage 6 2 logz
-title SHMS DC Drift Time vs. Wire Number
-pdc1u1_drifttime_vs_wirenum -nostat
-pdc1u2_drifttime_vs_wirenum -nostat
-pdc1x1_drifttime_vs_wirenum -nostat
-pdc1x2_drifttime_vs_wirenum -nostat
-pdc1v1_drifttime_vs_wirenum -nostat
-pdc1v2_drifttime_vs_wirenum -nostat
-pdc2v2_drifttime_vs_wirenum -nostat
-pdc2v1_drifttime_vs_wirenum -nostat
-pdc2x2_drifttime_vs_wirenum -nostat
-pdc2x1_drifttime_vs_wirenum -nostat
-pdc2u2_drifttime_vs_wirenum -nostat
-pdc2u1_drifttime_vs_wirenum -nostat
-newpage 6 2 
-title SHMS DC Plane Drift Distance
-pdc1u1_ddist -nostat
-pdc1u2_ddist -nostat
-pdc1x1_ddist -nostat
-pdc1x2_ddist -nostat
-pdc1v1_ddist -nostat
-pdc1v2_ddist -nostat
-pdc2v2_ddist -nostat
-pdc2v1_ddist -nostat
-pdc2x2_ddist -nostat
-pdc2x1_ddist -nostat
-pdc2u2_ddist -nostat
-pdc2u1_ddist -nostat
-newpage 6 2 logz
-title SHMS DC Drift Distance vs. Wire Number
-pdc1u1_wirenum_vs_ddist -nostat
-pdc1u2_wirenum_vs_ddist -nostat
-pdc1x1_wirenum_vs_ddist -nostat
-pdc1x2_wirenum_vs_ddist -nostat
-pdc1v1_wirenum_vs_ddist -nostat
-pdc1v2_wirenum_vs_ddist -nostat
-pdc2v2_wirenum_vs_ddist -nostat
-pdc2v1_wirenum_vs_ddist -nostat
-pdc2x2_wirenum_vs_ddist -nostat
-pdc2x1_wirenum_vs_ddist -nostat
-pdc2u2_wirenum_vs_ddist -nostat
-pdc2u1_wirenum_vs_ddist -nostat
+newpage 2 3
+title SHMS DC raw TDC
+macro    Macros/EXAMPLES/overlay2.C("pdc1u1_raw_tdc","pdc2v1_raw_tdc","1U1","2V1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1u2_raw_tdc","pdc2v2_raw_tdc","1U2","2V2")
+macro    Macros/EXAMPLES/overlay2.C("pdc1x1_raw_tdc","pdc2x1_raw_tdc","1X1","2X1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1x2_raw_tdc","pdc2x2_raw_tdc","1X2","2X2")
+macro    Macros/EXAMPLES/overlay2.C("pdc1v1_raw_tdc","pdc2u1_raw_tdc","1V1","2U1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1v2_raw_tdc","pdc2u2_raw_tdc","1V2","2U2")
 
-newpage 6 2 
-title SHMS DC Plane Residuals 
-pdc1u1_residuals -nostat
-pdc1u2_residuals -nostat
-pdc1x1_residuals -nostat
-pdc1x2_residuals -nostat
-pdc1v1_residuals -nostat
-pdc1v2_residuals -nostat
-pdc2v2_residuals -nostat
-pdc2v1_residuals -nostat
-pdc2x2_residuals -nostat
-pdc2x1_residuals -nostat
-pdc2u2_residuals -nostat
-pdc2u1_residuals -nostat
-newpage 6 2 logz
-title SHMS DC Residuals vs. Wire Number
-pdc1u1_residuals_vs_wirenum -nostat
-pdc1u2_residuals_vs_wirenum -nostat
-pdc1x1_residuals_vs_wirenum -nostat
-pdc1x2_residuals_vs_wirenum -nostat
-pdc1v1_residuals_vs_wirenum -nostat
-pdc1v2_residuals_vs_wirenum -nostat
-pdc2v2_residuals_vs_wirenum -nostat
-pdc2v1_residuals_vs_wirenum -nostat
-pdc2x2_residuals_vs_wirenum -nostat
-pdc2x1_residuals_vs_wirenum -nostat
-pdc2u2_residuals_vs_wirenum -nostat
-pdc2u1_residuals_vs_wirenum -nostat
-newpage 1 1 logz
-title X vs. Y at F.P.
+newpage 3 4 logz
+title SHMS DC raw TDC vs wire
+pdc1u1_raw_tdc_vs_wirenum -colz
+pdc1u2_raw_tdc_vs_wirenum -colz
+pdc1x1_raw_tdc_vs_wirenum -colz
+pdc1x2_raw_tdc_vs_wirenum -colz
+pdc1v1_raw_tdc_vs_wirenum -colz
+pdc1v2_raw_tdc_vs_wirenum -colz
+pdc2v2_raw_tdc_vs_wirenum -colz
+pdc2v1_raw_tdc_vs_wirenum -colz
+pdc2x2_raw_tdc_vs_wirenum -colz
+pdc2x1_raw_tdc_vs_wirenum -colz
+pdc2u2_raw_tdc_vs_wirenum -colz
+pdc2u1_raw_tdc_vs_wirenum -colz
+
+newpage 2 3
+title SHMS DC wire hits
+macro    Macros/EXAMPLES/overlay2.C("pdc1u1_wirenum","pdc2v1_wirenum","1U1","2V1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1u2_wirenum","pdc2v2_wirenum","1U2","2V2")
+macro    Macros/EXAMPLES/overlay2.C("pdc1x1_wirenum","pdc2x1_wirenum","1X1","2X1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1x2_wirenum","pdc2x2_wirenum","1X2","2X2")
+macro    Macros/EXAMPLES/overlay2.C("pdc1v1_wirenum","pdc2u1_wirenum","1V1","2U1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1v2_wirenum","pdc2u2_wirenum","1V2","2U2")
+
+newpage 2 3
+title SHMS DC drift distance
+macro    Macros/EXAMPLES/overlay2.C("pdc1u1_ddist","pdc2v1_ddist","1U1","2V1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1u2_ddist","pdc2v2_ddist","1U2","2V2")
+macro    Macros/EXAMPLES/overlay2.C("pdc1x1_ddist","pdc2x1_ddist","1X1","2X1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1x2_ddist","pdc2x2_ddist","1X2","2X2")
+macro    Macros/EXAMPLES/overlay2.C("pdc1v1_ddist","pdc2u1_ddist","1V1","2U1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1v2_ddist","pdc2u2_ddist","1V2","2U2")
+
+newpage 3 4 logz
+title SHMS DC drift dist vs wire
+pdc1u1_wirenum_vs_ddist -colz
+pdc1u2_wirenum_vs_ddist -colz
+pdc1x1_wirenum_vs_ddist -colz
+pdc1x2_wirenum_vs_ddist -colz
+pdc1v1_wirenum_vs_ddist -colz
+pdc1v2_wirenum_vs_ddist -colz
+pdc2v2_wirenum_vs_ddist -colz
+pdc2v1_wirenum_vs_ddist -colz
+pdc2x2_wirenum_vs_ddist -colz
+pdc2x1_wirenum_vs_ddist -colz
+pdc2u2_wirenum_vs_ddist -colz
+pdc2u1_wirenum_vs_ddist -colz
+
+newpage 2 3
+title SHMS DC drift time
+macro    Macros/EXAMPLES/overlay2.C("pdc1u1_drifttime","pdc2v1_drifttime","1U1","2V1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1u2_drifttime","pdc2v1_drifttime","1U2","2V2")
+macro    Macros/EXAMPLES/overlay2.C("pdc1x1_drifttime","pdc2x1_drifttime","1X1","2X1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1x2_drifttime","pdc2x2_drifttime","1X2","2X2")
+macro    Macros/EXAMPLES/overlay2.C("pdc1v1_drifttime","pdc2u1_drifttime","1V1","2U1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1v2_drifttime","pdc2u2_drifttime","1V2","2U2")
+
+newpage 3 4 logz
+title SHMS DC drift time vs wire
+pdc1u1_drifttime_vs_wirenum -colz
+pdc1u2_drifttime_vs_wirenum -colz
+pdc1x1_drifttime_vs_wirenum -colz
+pdc1x2_drifttime_vs_wirenum -colz
+pdc1v1_drifttime_vs_wirenum -colz
+pdc1v2_drifttime_vs_wirenum -colz
+pdc2v2_drifttime_vs_wirenum -colz
+pdc2v1_drifttime_vs_wirenum -colz
+pdc2x2_drifttime_vs_wirenum -colz
+pdc2x1_drifttime_vs_wirenum -colz
+pdc2u2_drifttime_vs_wirenum -colz
+pdc2u1_drifttime_vs_wirenum -colz
+
+newpage 2 3
+title SHMS DC residuals
+macro    Macros/EXAMPLES/overlay2.C("pdc1u1_residuals","pdc2v1_residuals","1U1","2V1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1u2_residuals","pdc2v2_residuals","1U2","2V2")
+macro    Macros/EXAMPLES/overlay2.C("pdc1x1_residuals","pdc2x1_residuals","1X1","2X1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1x2_residuals","pdc2x2_residuals","1X2","2X2")
+macro    Macros/EXAMPLES/overlay2.C("pdc1v1_residuals","pdc2u1_residuals","1V1","2U1")
+macro    Macros/EXAMPLES/overlay2.C("pdc1v2_residuals","pdc2u2_residuals","1V2","2U2")
+
+newpage 3 4 logz
+title SHMS DC residuals vs wire 
+pdc1u1_residuals_vs_wirenum -colz
+pdc1u2_residuals_vs_wirenum -colz
+pdc1x1_residuals_vs_wirenum -colz
+pdc1x2_residuals_vs_wirenum -colz
+pdc1v1_residuals_vs_wirenum -colz
+pdc1v2_residuals_vs_wirenum -colz
+pdc2v2_residuals_vs_wirenum -colz
+pdc2v1_residuals_vs_wirenum -colz
+pdc2x2_residuals_vs_wirenum -colz
+pdc2x1_residuals_vs_wirenum -colz
+pdc2u2_residuals_vs_wirenum -colz
+pdc2u1_residuals_vs_wirenum -colz
+
+newpage 2 2 
+title SHMS DC track at focal plane
 pdc_xfp_vs_yxp -nostat
-newpage 1 1 logz
-title Xp vs. Yp at F.P
-pdc_xfpp_vs_yfpp -nostat
+pdc_xpfp_vs_ypfp -nostat
+pdc_xfp_vs_ypfp -nostat
+pdc_xfp_vs_xpfp -nostat
\ No newline at end of file

onlineGUI/Macros/EXAMPLES/overlay2.C

+/*******************************
+ * This macro overlays two histograms.
+ * The label1 and label2 are optional
+ * and can be used to label the
+ * different histograms as they are 
+ * drawn.
+ *
+ * hszumila@jlab.org
+ *******************************/
+
+void overlay2(TString hist1name, TString hist2name, TString label1="", TString label2="", Bool_t golden=false,TString samestr=""){
+  TH1F* H1;
+  TH1F* H2;
+  
+  H1 = (TH1F*) gDirectory->Get(hist1name);
+  H2 = (TH1F*) gDirectory->Get(hist2name);
+  if (H1 && H2) {
+  if(golden){
+    H1->SetLineColor(30);
+    H1->SetFillColor(30);
+    H1->SetFillStyle(3027);
+    H2->SetLineColor(46);
+    H2->SetFillColor(46);
+    H2->SetFillStyle(3027);
+  } else {
+    H1->SetLineColor(4);
+    H1->SetLineWidth(2);
+    H2->SetLineWidth(2);
+    H2->SetLineColor(2);
+    H1->SetStats(0);
+    H2->SetStats(0);
+    H1->SetMinimum(0);
+    H2->SetMinimum(0);
+    H1->SetTitle("");
+    H2->SetTitle("");
+    H1->GetXaxis()->SetTitleSize(0.05);
+    H1->GetYaxis()->SetTitleSize(0.045);
+    H2->GetXaxis()->SetTitleSize(0.05);
+    H2->GetYaxis()->SetTitleSize(0.045);
+    H1->GetXaxis()->SetLabelSize(0.04);
+    H1->GetYaxis()->SetLabelSize(0.04);
+    H2->GetXaxis()->SetLabelSize(0.04);
+    H2->GetYaxis()->SetLabelSize(0.04);
+
+  }
+ 
+  double xpos = H1->GetXaxis()->GetBinCenter(H1->GetXaxis()->GetFirst()+4);
+  double ypos = 0.95*H1->GetMaximum();
+
+  TText* t1 = new TText(xpos,ypos,label1);
+  t1->SetTextColor(4);
+  TText* t2 = new TText(xpos,0.9*ypos,label2);
+  t2->SetTextColor(2);
+  
+  if(golden){
+    if(!samestr.Contains("sames")) samestr += "sames";
+    H1->Draw(samestr);
+    H2->Draw(samestr);
+  } else {
+    if(!golden && !samestr.Contains("sames")){
+      // pick which hist to draw first:
+      if(H1->GetMaximum() > H2->GetMaximum()){
+	H1->Draw(samestr);
+	if(!samestr.Contains("sames")) samestr += "sames";
+	H2->Draw(samestr);
+	t1->Draw();
+	t2->Draw();
+      } else {
+	H2->Draw(samestr);
+	if(!samestr.Contains("sames")) samestr += "sames";
+	H1->Draw(samestr);
+	t1->Draw();
+	t2->Draw();
+      }
+    } else {
+      H1->Draw(samestr);
+      if(!samestr.Contains("sames")) samestr += "sames";
+      H2->Draw(samestr);
+      t1->Draw();
+      t2->Draw();
+
+    H2->GetXaxis()->SetTitleOffset(.6);
+    H2->GetXaxis()->SetTitleSize(0.08);
+    H2->GetYaxis()->SetTitleOffset(.6);
+    H2->GetYaxis()->SetTitleSize(0.08);
+    }
+  }
+  } else {
+    if (!H1) cout << " Histogram " << hist1name << " does not exist" << endl;
+    if (!H2) cout << " Histogram " << hist2name << " does not exist" << endl;
+  }
+}