Merge pull request #15 from pooser/pooser-fix-replay-for-new-data

Pooser fix replay for new data

DBASE/standard.database

 0-99999
 g_ctp_parm_filename ="DBASE/general.param"
-g_decode_map_filename ="MAPS/hmschambers.map"
+g_decode_map_filename ="MAPS/HMS/hdc.map"
 g_ctp_kinematics_filename ="DBASE/standard.kinematics"

DEF-files/hdcana.def

@@ -17,16 +17,44 @@ TH1F hdc2v1_wm 'HDC 2V1 Wiremap' H.dc.2v1.wirenum 107 0.5 107.5
 TH1F hdc2y2_wm 'HDC 2Y2 Wiremap' H.dc.2y2.wirenum  52 0.5  52.5
 TH1F hdc2x2_wm 'HDC 2X2 Wiremap' H.dc.2x2.wirenum 113 0.5 113.5
 
-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
+TH2F hdc1x1_wire_num_vs_raw_tdc 'HDC 1X1 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.1x1.wirenum H.dc.1x1.rawtdc 113 0.5 113.5 5000 0 50000
+TH2F hdc1y1_wire_num_vs_raw_tdc 'HDC 1Y1 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.1y1.wirenum H.dc.1y1.rawtdc 52  0.5  52.5 5000 0 50000
+TH2F hdc1u1_wire_num_vs_raw_tdc 'HDC 1U1 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.1u1.wirenum H.dc.1u1.rawtdc 107 0.5 107.5 5000 0 50000
+TH2F hdc1v1_wire_num_vs_raw_tdc 'HDC 1V1 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.1v1.wirenum H.dc.1v1.rawtdc 107 0.5 107.5 5000 0 50000
+TH2F hdc1y2_wire_num_vs_raw_tdc 'HDC 1Y2 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.1y2.wirenum H.dc.1y2.rawtdc 52  0.5  52.5 5000 0 50000
+TH2F hdc1x2_wire_num_vs_raw_tdc 'HDC 1X2 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.1x2.wirenum H.dc.1x2.rawtdc 113 0.5 113.5 5000 0 50000
+
+TH2F hdc2x1_wire_num_vs_raw_tdc 'HDC 2X1 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.2x1.wirenum H.dc.2x1.rawtdc 113 0.5 113.5 5000 0 50000
+TH2F hdc2y1_wire_num_vs_raw_tdc 'HDC 2Y1 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.2y1.wirenum H.dc.2y1.rawtdc 52  0.5  52.5 5000 0 50000
+TH2F hdc2u1_wire_num_vs_raw_tdc 'HDC 2U1 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.2u1.wirenum H.dc.2u1.rawtdc 107 0.5 107.5 5000 0 50000
+TH2F hdc2v1_wire_num_vs_raw_tdc 'HDC 2V1 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.2v1.wirenum H.dc.2v1.rawtdc 107 0.5 107.5 5000 0 50000
+TH2F hdc2y2_wire_num_vs_raw_tdc 'HDC 2Y2 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.2y2.wirenum H.dc.2y2.rawtdc 52  0.5  52.5 5000 0 50000
+TH2F hdc2x2_wire_num_vs_raw_tdc 'HDC 2X2 Wiremap; Wire Number; Raw TDC Value (clicks)' H.dc.2x2.wirenum H.dc.2x2.rawtdc 113 0.5 113.5 5000 0 50000
+
+TH2F hdc1x1_wire_num_vs_tdc_time 'HDC 1X1 Wiremap; Wire Number; TDC Time (ns)' H.dc.1x1.wirenum H.dc.1x1.time 113 0.5 113.5 6000 -1000 5000
+TH2F hdc1y1_wire_num_vs_tdc_time 'HDC 1Y1 Wiremap; Wire Number; TDC Time (ns)' H.dc.1y1.wirenum H.dc.1y1.time 52  0.5  52.5 6000 -1000 5000
+TH2F hdc1u1_wire_num_vs_tdc_time 'HDC 1U1 Wiremap; Wire Number; TDC Time (ns)' H.dc.1u1.wirenum H.dc.1u1.time 107 0.5 107.5 6000 -1000 5000
+TH2F hdc1v1_wire_num_vs_tdc_time 'HDC 1V1 Wiremap; Wire Number; TDC Time (ns)' H.dc.1v1.wirenum H.dc.1v1.time 107 0.5 107.5 6000 -1000 5000
+TH2F hdc1y2_wire_num_vs_tdc_time 'HDC 1Y2 Wiremap; Wire Number; TDC Time (ns)' H.dc.1y2.wirenum H.dc.1y2.time 52  0.5  52.5 6000 -1000 5000
+TH2F hdc1x2_wire_num_vs_tdc_time 'HDC 1X2 Wiremap; Wire Number; TDC Time (ns)' H.dc.1x2.wirenum H.dc.1x2.time 113 0.5 113.5 6000 -1000 5000
+
+TH2F hdc2x1_wire_num_vs_tdc_time 'HDC 2X1 Wiremap; Wire Number; TDC Time (ns)' H.dc.2x1.wirenum H.dc.2x1.time 113 0.5 113.5 6000 -1000 5000
+TH2F hdc2y1_wire_num_vs_tdc_time 'HDC 2Y1 Wiremap; Wire Number; TDC Time (ns)' H.dc.2y1.wirenum H.dc.2y1.time 52  0.5  52.5 6000 -1000 5000
+TH2F hdc2u1_wire_num_vs_tdc_time 'HDC 2U1 Wiremap; Wire Number; TDC Time (ns)' H.dc.2u1.wirenum H.dc.2u1.time 107 0.5 107.5 6000 -1000 5000
+TH2F hdc2v1_wire_num_vs_tdc_time 'HDC 2V1 Wiremap; Wire Number; TDC Time (ns)' H.dc.2v1.wirenum H.dc.2v1.time 107 0.5 107.5 6000 -1000 5000
+TH2F hdc2y2_wire_num_vs_tdc_time 'HDC 2Y2 Wiremap; Wire Number; TDC Time (ns)' H.dc.2y2.wirenum H.dc.2y2.time 52  0.5  52.5 6000 -1000 5000
+TH2F hdc2x2_wire_num_vs_tdc_time 'HDC 2X2 Wiremap; Wire Number; TDC Time (ns)' H.dc.2x2.wirenum H.dc.2x2.time 113 0.5 113.5 6000 -1000 5000
+
+TH2F hdc1x1_wire_num_vs_ddist 'HDC 1X1 Wiremap; Wire Number; Drift Distance (cm)' H.dc.1x1.wirenum H.dc.1x1.dist 113 0.5 113.5 300 -0.1 0.6
+TH2F hdc1y1_wire_num_vs_ddist 'HDC 1Y1 Wiremap; Wire Number; Drift Distance (cm)' H.dc.1y1.wirenum H.dc.1y1.dist 52  0.5  52.5 300 -0.1 0.6
+TH2F hdc1u1_wire_num_vs_ddist 'HDC 1U1 Wiremap; Wire Number; Drift Distance (cm)' H.dc.1u1.wirenum H.dc.1u1.dist 107 0.5 107.5 300 -0.1 0.6
+TH2F hdc1v1_wire_num_vs_ddist 'HDC 1V1 Wiremap; Wire Number; Drift Distance (cm)' H.dc.1v1.wirenum H.dc.1v1.dist 107 0.5 107.5 300 -0.1 0.6
+TH2F hdc1y2_wire_num_vs_ddist 'HDC 1Y2 Wiremap; Wire Number; Drift Distance (cm)' H.dc.1y2.wirenum H.dc.1y2.dist 52  0.5  52.5 300 -0.1 0.6
+TH2F hdc1x2_wire_num_vs_ddist 'HDC 1X2 Wiremap; Wire Number; Drift Distance (cm)' H.dc.1x2.wirenum H.dc.1x2.dist 113 0.5 113.5 300 -0.1 0.6
+
+TH2F hdc2x1_wire_num_vs_ddist 'HDC 2X1 Wiremap; Wire Number; Drift Distance (cm)' H.dc.2x1.wirenum H.dc.2x1.dist 113 0.5 113.5 300 -0.1 0.6
+TH2F hdc2y1_wire_num_vs_ddist 'HDC 2Y1 Wiremap; Wire Number; Drift Distance (cm)' H.dc.2y1.wirenum H.dc.2y1.dist 52  0.5  52.5 300 -0.1 0.6
+TH2F hdc2u1_wire_num_vs_ddist 'HDC 2U1 Wiremap; Wire Number; Drift Distance (cm)' H.dc.2u1.wirenum H.dc.2u1.dist 107 0.5 107.5 300 -0.1 0.6
+TH2F hdc2v1_wire_num_vs_ddist 'HDC 2V1 Wiremap; Wire Number; Drift Distance (cm)' H.dc.2v1.wirenum H.dc.2v1.dist 107 0.5 107.5 300 -0.1 0.6
+TH2F hdc2y2_wire_num_vs_ddist 'HDC 2Y2 Wiremap; Wire Number; Drift Distance (cm)' H.dc.2y2.wirenum H.dc.2y2.dist 52  0.5  52.5 300 -0.1 0.6
+TH2F hdc2x2_wire_num_vs_ddist 'HDC 2X2 Wiremap; Wire Number; Drift Distance (cm)' H.dc.2x2.wirenum H.dc.2x2.dist 113 0.5 113.5 300 -0.1 0.6

PARAM/hdc_test_stand.param

+; Custom parameter file which should loaded when aiming to
+; analyze HMS drift chamber data with no tracking.
+
+
+; Increase HMS TDC window size
+  hdc_tdc_min_win = -50000,-50000,-50000,-50000,-50000,-50000
+  		    -50000,-50000,-50000,-50000,-50000,-50000
+  hdc_tdc_max_win = 50000,50000,50000,50000,50000,50000
+                    50000,50000,50000,50000,50000,50000
+; Remove hodoscope dependence on tracking
+  hsel_using_scin  = 0
+;  hsel_using_prune = 0
+
+;  Adjust HMS DC plane t = 0
+;   hdc_plane_time_zero = 1500.0,1500.0,1500.0,1500.0,1500.0,1500.0
+;  	                 1500.0,1500.0,1500.0,1500.0,1500.0,1500.0
+

SCRIPTS/replay_hms.C

@@ -17,7 +17,7 @@ void replay_hms(Int_t RunNumber=0, Int_t MaxEvent=0) {
 
   // Create file name patterns.
   const char* RunFileNamePattern = "raw/test_%d.dat";
-  const char* ROOTFileNamePattern = "ROOTfiles/hms1190_%d.root";
+  const char* ROOTFileNamePattern = "ROOTfiles/test_%d.root";
   // Add variables to global list.
   gHcParms->Define("gen_run_number", "Run Number", RunNumber);
   gHcParms->AddString("g_ctp_database_filename", "DBASE/standard.database");
@@ -29,6 +29,9 @@ void replay_hms(Int_t RunNumber=0, Int_t MaxEvent=0) {
   gHcParms->Load(gHcParms->GetString("g_ctp_kinematics_filename"), RunNumber);
   gHcParms->Load(gHcParms->GetString("g_ctp_parm_filename"));
 
+  // Load params for HMS DC test stand configuration
+  gHcParms->Load("PARAM/hdc_test_stand.param");
+
   // Generate db_cratemap to correspond to map file contents via Pearl script.
   //char command[100];
   //sprintf(