Dc calib develop (#135)

* -added the lines to read run number to the top of the script.

* -Modified SHMS DC DEF FILE to display drift times with binning scheme same as used in calib.
-Added the most up-to-date shms dc parameters obtained from the modified shms dc calib. script

* -set option to use per_wire tzero ON.

CALIBRATION/shms_dc_calib/scripts/get_pdc_time_histo_tzero_corrected.C

@@ -4,6 +4,13 @@ using namespace std;
 void get_pdc_time_histo_tzero_corrected()
 {
   
+  
+  int run_NUM;
+  TString f0 = "input_RUN.txt";
+  ifstream infile(f0);
+  infile >> run_NUM;   
+
+  TString run = Form("run%d", run_NUM);
 
   //this script reads all tzero values inside tzero_group.dat and assigns them to their corresponding planes. 
   //these values will then be used to shift the drift time histos, wire by wire (t0-correction wire-by-wire)
@@ -23,7 +30,7 @@ void get_pdc_time_histo_tzero_corrected()
 
   //open and read tzero data file 
   ifstream file;
-  file.open("../data_files/run484/tzero_group.dat");
+  file.open("../data_files/"+run+"/tzero_group.dat");
   
   string line;
   int counter;
@@ -86,13 +93,6 @@ void get_pdc_time_histo_tzero_corrected()
   //THIS SCRIPT WILL READ the drift times array tzero[ip][sw] on a wire basis and apply the tzero correction.
 
 
-  
-  int run_NUM;
-  TString f0 = "input_RUN.txt";
-  ifstream infile(f0);
-  infile >> run_NUM;   
-
-  TString run = Form("run%d", run_NUM);
 
   int i;
  

DEF-files/SHMS/GEN/pstackana_cuts.def

@@ -9,5 +9,19 @@ Decode_master 1
 Block: CoarseTracking
 CoarseTracking_master 1
 
+#add cuts to stricter times and apply them to drift distances 
+time_cut1                            P.dc.1u1.time > 0 && P.dc.1u1.time < 250
+time_cut2                            P.dc.1u2.time > 0 && P.dc.1u2.time < 250
+time_cut3                            P.dc.1x1.time > 0 && P.dc.1x1.time < 250
+time_cut4                            P.dc.1x2.time > 0 && P.dc.1x2.time < 250
+time_cut5                            P.dc.1v1.time > 0 && P.dc.1v1.time < 250
+time_cut6                            P.dc.1v2.time > 0 && P.dc.1v2.time < 250
+time_cut7                            P.dc.2v2.time > 0 && P.dc.2v2.time < 250
+time_cut8                            P.dc.2v1.time > 0 && P.dc.2v1.time < 250
+time_cut9                            P.dc.2x2.time > 0 && P.dc.2x2.time < 250
+time_cut10                            P.dc.2x1.time > 0 && P.dc.2x1.time < 250
+time_cut11                            P.dc.2u2.time > 0 && P.dc.2u2.time < 250
+time_cut12                            P.dc.2u1.time > 0 && P.dc.2u1.time < 250
+
 Block: CoarseReconstruct
 CourseReconstruct 0

PARAM/SHMS/DC/pdc.param

@@ -18,4 +18,4 @@ pdc_chamber_planes = 1, 1, 1, 1, 1, 1
 pdc_nrwire = 107, 107, 79, 79, 107, 107
              107, 107, 79, 79, 107, 107
 ; 1 means true
-p_using_tzero_per_wire =0
+p_using_tzero_per_wire = 1

PARAM/SHMS/DC/pdc_tracking.param

@@ -27,8 +27,14 @@ pdc_tdc_time_per_channel = -0.10
 ;                      1283.94, 1283.34, 1286.95, 1286.46, 1286.37, 1284.47
 
 ;Zero time correction for each plane in ns that is added to TDC time.
-pdc_plane_time_zero = 1293.66, 1293.5, 1286.86, 1294.36, 1291.36, 1293.05
-                      1289.26, 1290.03, 1287.02, 1286.41, 1290.52, 1288.88
+;pdc_plane_time_zero = 1293.66, 1293.5, 1286.86, 1294.36, 1291.36, 1293.05
+;                      1289.26, 1290.03, 1287.02, 1286.41, 1290.52, 1288.88
+
+;Zero time correction for each plane in ns that is added to TDC time.
+pdc_plane_time_zero = 1290.00, 1290.00, 1290.00, 1290.00, 1290.00, 1290.00
+                      1290.00, 1290.00, 1290.00, 1290.00, 1290.00, 1290.00
+
+
 
 
 ; For wire velocity corrections.