From 17f470082df24370def74328945f679f4cb16370 Mon Sep 17 00:00:00 2001
From: Yero1990 <cyero002@fiu.edu>
Date: Wed, 14 Jun 2017 07:10:35 -0400
Subject: [PATCH] Shms dc calib (#138)

* -added script that will output a parameter file with individual tzero values (NOT per ribbon cable as before)

* -included the execution of 'get_tzero_per_wire.param' in steering script.

* -added precicion to the t0 calculation. Although the calculation involved the ratios of two 'doubles', the t0 came out to be an integer. This imprecision was corrected.

* -add minor fix

* added minor fix to read in corrected t0 values

* minor fix

* most up-to-date shms DC param files
---
 CALIBRATION/shms_dc_calib/run_Cal.C           |   5 +-
 .../get_pdc_time_histo_tzero_corrected.C      |  10 +-
 .../scripts/get_tzero_per_wire_param.C        | 135 +++++++
 .../shms_dc_calib/scripts/get_wire_tzero.C    |  55 +--
 PARAM/SHMS/DC/pdc_tzero_per_wire.param        | 164 ++++-----
 PARAM/SHMS/DC/pdriftmap.param                 | 336 +++++++++---------
 6 files changed, 420 insertions(+), 285 deletions(-)
 create mode 100644 CALIBRATION/shms_dc_calib/scripts/get_tzero_per_wire_param.C

diff --git a/CALIBRATION/shms_dc_calib/run_Cal.C b/CALIBRATION/shms_dc_calib/run_Cal.C
index eb209b24..de9a9856 100644
--- a/CALIBRATION/shms_dc_calib/run_Cal.C
+++ b/CALIBRATION/shms_dc_calib/run_Cal.C
@@ -62,8 +62,9 @@ void run_Cal()
   // gSystem->Exec("root -l -q update_pdcparam.C");
   
   //execute code to update parameter file with new t0 values per wire
-  gSystem->Exec("root -l -q get_tzero_param.C");
-
+  //gSystem->Exec("root -l -q get_tzero_param.C");
+  gSystem->Exec("root -l -q get_tzero_per_wire_param.C"); 
+  
   //execute code to get t0 corrected drift times
   gSystem->Exec("root -l -q get_pdc_time_histo_tzero_corrected.C");
   
diff --git a/CALIBRATION/shms_dc_calib/scripts/get_pdc_time_histo_tzero_corrected.C b/CALIBRATION/shms_dc_calib/scripts/get_pdc_time_histo_tzero_corrected.C
index f46f9357..90d4af44 100644
--- a/CALIBRATION/shms_dc_calib/scripts/get_pdc_time_histo_tzero_corrected.C
+++ b/CALIBRATION/shms_dc_calib/scripts/get_pdc_time_histo_tzero_corrected.C
@@ -30,7 +30,7 @@ void get_pdc_time_histo_tzero_corrected()
 
   //open and read tzero data file 
   ifstream file;
-  file.open("../data_files/"+run+"/tzero_group.dat");
+  file.open("../data_files/"+run+"/tzero_wire_NEW.dat");
   
   string line;
   int counter;
@@ -46,11 +46,11 @@ void get_pdc_time_histo_tzero_corrected()
       {
 	sscanf(line.c_str(), "%lf", &value);  //pass data in file to variable 'value'
 	tzero_offset[counter] = value;  //write all tzero values for all wires in both DC to an array
-	//	cout << tzero_offset[counter] << endl;
+	//cout << tzero_offset[counter] << endl;
 	counter++;
       }
   }
-
+  
 
   //***************************************************************
   
@@ -154,7 +154,7 @@ void get_pdc_time_histo_tzero_corrected()
 
      shift = tzero[ip][sw-1];  //the shift represents how much the drift time histo needs to be offset
 
-     //cout << "sw: " << sw << " :: " << "offset: " << shift << endl;
+	//cout << "sw: " << sw << " :: " << "offset: " << shift << endl;
     
      //************APPLY TZERO OFFSET ON A WIRE-BY-WIRE BASIS TO ALL WIRES IN ALL PLANES***************
 
@@ -200,6 +200,6 @@ void get_pdc_time_histo_tzero_corrected()
  
  
  
-
+  
 
 }
diff --git a/CALIBRATION/shms_dc_calib/scripts/get_tzero_per_wire_param.C b/CALIBRATION/shms_dc_calib/scripts/get_tzero_per_wire_param.C
new file mode 100644
index 00000000..de38c431
--- /dev/null
+++ b/CALIBRATION/shms_dc_calib/scripts/get_tzero_per_wire_param.C
@@ -0,0 +1,135 @@
+#include <iostream>
+#include <fstream>
+#include <math.h>
+#include <string>
+
+#define NPLANES 12
+
+
+void get_tzero_per_wire_param()
+{
+
+
+ //read run number from input file
+  int run_NUM;
+  TString f0 = "input_RUN.txt";
+  ifstream infile(f0);
+  infile >> run_NUM;   
+
+  TString run = Form("run%d", run_NUM);
+
+
+  TString planes[NPLANES] =   {"1u1","1u2","1x1","1x2","1v1","1v2","2v2","2v1","2x2","2x1","2u2","2u1"};   //planes 0, 1, 2, 3, 4, 5 OK :: 
+
+
+int fNWires[NPLANES] = {107, 107, 79, 79, 107, 107, 107, 107, 79, 79, 107, 107};
+
+int sw;
+int ip;
+
+Int_t **wire = new Int_t*[NPLANES];
+Double_t **t0 = new Double_t*[NPLANES];
+//Double_t **t0_err = new Double_t*[NPLANES];
+//Int_t **entries = new Int_t*[NPLANES];
+
+ for (ip = 0; ip < NPLANES; ip++) {
+   t0[ip] = new Double_t[fNWires[ip]];
+   for (sw=0; sw<fNWires[ip]; sw++) {
+     t0[ip][sw] = 0.0;
+   }
+ }
+
+  
+
+ 
+//Loop over each plane
+ for(ip=0; ip<NPLANES; ip++) {
+   
+         //write plane headers
+
+      Int_t nwire;
+      Double_t nt0;
+      Double_t nt0_err;
+      Int_t nentries;
+
+
+
+    //open and read each wire tzero file
+    string line;
+    ifstream input;
+    input.open(Form("../data_files/run%d/hdc_"+planes[ip]+"tzero_run%d_updated.txt", run_NUM, run_NUM ) );
+    
+    
+     sw = 0;  //se wire counter to 0
+    
+      while (getline(input, line))
+	{
+	 
+	  if (line != '#')
+	    {	     
+	      input >> nwire >> nt0 >> nt0_err >> nentries;
+	     // cout << nwire << " :: " << nt0 <<  endl;		       
+	      //cout << ip << " :: " << sw << endl;
+	   
+	      t0[ip][nwire-1] = nt0;
+	      
+	       sw++;		  
+	    }
+	  
+	} 
+   
+      input.close();
+    
+ } //end loop over planes
+ 
+ 
+//create output file stream to write tzero values to a data file
+ ofstream out;
+ TString tzero_file = "../data_files/"+run+"/tzero_wire_NEW.dat";
+ out.open(tzero_file);
+ 
+ //Create an output parameter file to store tzero values 
+ ofstream ofs;
+ TString wire_tzero = Form("../../../PARAM/SHMS/DC/pdc_tzero_per_wire_run%d_NEW.param", run_NUM);
+ ofs.open (wire_tzero);
+ 
+
+ for (ip=0; ip<NPLANES; ip++) { 
+
+   //write plane headers
+   ofs << "ptzero"+planes[ip] << "=" << endl;
+   out << "#PLANE: " << planes[ip] << endl;
+
+   for (sw=0; sw<fNWires[ip]; sw++) {
+
+     out << setprecision(3) << t0[ip][sw] << endl;
+     
+     if (sw <= 10) 
+       { 
+	 ofs << setprecision(3) << t0[ip][sw] << fixed << ",";
+       }
+     else if (sw>10 && sw <(fNWires[ip]-1))
+       {
+	 ofs << setprecision(3) << t0[ip][sw] << ((sw+1) % 16 ? ", " : "\n") << fixed;
+       }
+     else if (sw==78 || sw == 106) {
+       ofs << setprecision(3) << t0[ip][sw] << fixed << endl;
+      }
+
+    
+
+   }
+ }
+}
+ 
+ 
+
+
+
+ 
+
+
+
+
+
+
diff --git a/CALIBRATION/shms_dc_calib/scripts/get_wire_tzero.C b/CALIBRATION/shms_dc_calib/scripts/get_wire_tzero.C
index 4c31f0d1..66d34a54 100644
--- a/CALIBRATION/shms_dc_calib/scripts/get_wire_tzero.C
+++ b/CALIBRATION/shms_dc_calib/scripts/get_wire_tzero.C
@@ -29,6 +29,7 @@ void get_wire_tzero()
 
  //Declare plane names to loop over
  TString plane_names[NPLANES]={"1u1", "1u2", "1x1", "1x2", "1v1", "1v2", "2v2", "2v1", "2x2", "2x1", "2u2", "2u1"};
+ int fNWires[NPLANES] = {107, 107, 79, 79, 107, 107, 107, 107, 79, 79, 107, 107};
 
  //Declare a root file array to store individual DC cell drift times
  TString root_file;
@@ -233,9 +234,9 @@ void get_wire_tzero()
      t_zero = - y_int/m;
      t_zero_err = sqrt(y_int_err*y_int_err/(m*m) + y_int*y_int*m_err*m_err/(m*m*m*m) );
      entries = cell_dt[sensewire-1]->GetEntries();  //number of entries (triggers) per wire
-     
+     //cout << "y_int: " << y_int << " :: " << "m: " << m << " :: " << "t0: " << setprecision(6) << -y_int/m << endl;
      //Write "t0" values to file
-     ofs << sensewire << "          " << t_zero << "          " << t_zero_err << "          " << entries << endl;
+     ofs << sensewire << "          " << setprecision(6) << -y_int/m << "          " <<  t_zero_err << "          " << entries << endl;
      
      //Change to output root file and write fitted histos to file
      g->cd();
@@ -243,7 +244,7 @@ void get_wire_tzero()
      
    }
    
-   // Make Plot of t0 versus Wire Number 
+   /*    // Make Plot of t0 versus Wire Number 
    
    TCanvas *t = new TCanvas("t", "", 2000,500);
    t->SetGrid();
@@ -266,12 +267,12 @@ void get_wire_tzero()
    t->Write(title);   //write to a root file
    
    //close dat file
-   ofs.close();
+   
    //save plots
    //TString tzero_plots = "plots/"+run_NUM +"/hdc"+plane_names[ip]+Form("TESTING_tzero_v_wire_%d.eps", run);
    //t->SaveAs(tzero_plots);
-   
-   
+   */
+   ofs.close();
    //*****************************************************************************************//
    //        CALCULATE THE "t0s" WEIGHTED AVERAGE FOR WIRE DRIFT TIMES WITH ENTRIES > = 300   //
    //*****************************************************************************************//
@@ -293,6 +294,9 @@ void get_wire_tzero()
    double weighted_AVG;
    double weighted_AVG_err; 
   
+    int counter;
+    double t0_corr;
+    double t0_corr_err;
    //set them to zero to start sum inside while loop 
    sum_NUM = 0.0;
    sum_DEN = 0.0;
@@ -300,21 +304,27 @@ void get_wire_tzero()
    weighted_AVG =0.0 ;
    weighted_AVG_err= 0.0; 
    
+   counter = 0;
    //read line bt line the t_zero_file
    while(getline(ifs, line)) {
-     if(!line.length()|| line[0] == '#')
-       continue;
-     //	sensewire = 0, t_zero = 0.0, t_zero_err = 0.0, entries = 0 ; //set values to zero
-     
-     sscanf(line.c_str(), "%d %d %lf %d", &sensewire, &t_zero, &t_zero_err, &entries); //assign each of the variables above a data in the t_zero_file
-     
+
+     if(line!='#') {
+       
+       sscanf(line.c_str(), "%d %lf %lf %d", &sensewire, &t0_corr, &t0_corr_err, &entries); //assign each of the variables above a data in the t_zero_file
+ 
+       if(sensewire<=fNWires[ip]){
+	
+
      //Check if entries for each sensewire exceeds a certain number of events
      
-     if (entries>300 && t_zero < 30) {
+     if (entries>300 && t_zero < 30) 
+{
+
 	
+     	
        //Calculate the weighted average of t0s
-       sum_NUM = sum_NUM + t_zero/(t_zero_err*t_zero_err);
-       sum_DEN = sum_DEN + 1.0/(t_zero_err*t_zero_err);      
+       sum_NUM = sum_NUM + t0_corr/(t0_corr_err*t0_corr_err);
+       sum_DEN = sum_DEN + 1.0/(t0_corr_err*t0_corr_err);      
        
        //cout << "sum_NUM : " << sum_NUM << endl;  
        //cout << "sum_DEN : " << sum_DEN << endl;  
@@ -322,14 +332,15 @@ void get_wire_tzero()
     
 
 
-       ofs << sensewire << "        " << t_zero << "        " << t_zero_err << "        " << entries << endl;
-
+       ofs << sensewire << "        " << t0_corr << "        " << t0_corr_err << "        " << entries << endl;
+       //cout << "TZERO: " << t_zero << endl;
        
        
      }
      
      else { ofs << sensewire << "        " << 0.0 << "       " << 0.0 << "        " << entries << endl;}
-     
+       }
+     } //end if statement
    }
    
    
@@ -355,8 +366,8 @@ void get_wire_tzero()
    
    ifs.close();
 
-   // Make Plot of t0 versus Wire Number for entries > 300 events
-
+    // Make Plot of t0 versus Wire Number for entries > 300 events
+   /*
    TCanvas *t1 = new TCanvas("t1", "", 2000,500);
    t1->SetGrid();
 
@@ -388,9 +399,9 @@ void get_wire_tzero()
    ltx1->DrawLatex(t1->GetUxmax()*0.75,40, Form("Weighted Average = %lf #pm %lf ns", weighted_AVG, weighted_AVG_err) );
    
    t1->Write(title1);   //write canvas to a root file
-   
+   */
    ofs.close();  //close data file
-
+   
    
 
 
diff --git a/PARAM/SHMS/DC/pdc_tzero_per_wire.param b/PARAM/SHMS/DC/pdc_tzero_per_wire.param
index d8ec3cce..3e8bd248 100644
--- a/PARAM/SHMS/DC/pdc_tzero_per_wire.param
+++ b/PARAM/SHMS/DC/pdc_tzero_per_wire.param
@@ -1,100 +1,88 @@
 ptzero1u1=
-0, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000
-1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526, 1.7847526
-2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375, 2.2139375
-2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357, 2.3320357
-3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510, 3.2052510
-3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186, 3.7202186
-7.3742278, 7.3742278, 7.3742278, 7.3742278, 7.3742278, 7.3742278, 7.3742278, 7.3742278, 7.3742278, 7.3742278, 7.3742278, 7.3742278
-
+0,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000, 0.000, 0.000, 0.000, 0.000
+0.000, 0.000, 1.298, -1.640, 1.136, -1.287, 1.637, -1.226, -6.179, -4.345, -5.184, 2.044, -2.931, -1.992, 0.391, 0.745
+-1.456, 1.647, 1.724, -2.288, -1.218, -0.510, -1.422, -1.248, -0.866, -0.958, -1.121, -0.975, -0.644, 1.606, -1.298, -0.868
+1.148, 1.575, -1.189, -1.024, -1.393, -1.589, 1.701, -1.379, -1.276, -1.139, -1.097, -0.825, -1.861, -1.312, -1.546, -1.257
+1.098, 1.541, 2.157, 1.863, -0.965, 4.191, -0.232, 2.516, 3.822, 1.190, -0.338, 0.309, 2.887, 5.467, 1.324, 3.491
+3.631, -1.121, 1.147, 2.870, 2.803, 1.791, 1.561, 2.720, 5.830, 3.990, -2.890, -2.679, -0.333, 4.363, 2.898, 7.216
+1.607, 4.099, -18.641, -5.702, -7.432, -3.305, 0.000, 0.000, 0.000, 0.000, 0.000
 ptzero1u2=
-0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 2.6770076, 2.6770076, 2.6770076
-2.6770076, 2.6770076, 2.6770076, 2.6770076, 2.6770076, 2.6770076, 2.6770076, 2.6770076, 2.6770076, 2.6770076, 2.6770076, 2.6770076, 2.6770076, 2.2077937, 2.2077937, 2.2077937
-2.2077937, 2.2077937, 2.2077937, 2.2077937, 2.2077937, 2.2077937, 2.2077937, 2.2077937, 2.2077937, 2.2077937, 2.2077937, 2.2077937, 2.2077937, 1.8137457, 1.8137457, 1.8137457
-1.8137457, 1.8137457, 1.8137457, 1.8137457, 1.8137457, 1.8137457, 1.8137457, 1.8137457, 1.8137457, 1.8137457, 1.8137457, 1.8137457, 1.8137457, 2.8732654, 2.8732654, 2.8732654
-2.8732654, 2.8732654, 2.8732654, 2.8732654, 2.8732654, 2.8732654, 2.8732654, 2.8732654, 2.8732654, 2.8732654, 2.8732654, 2.8732654, 2.8732654, 3.4385692, 3.4385692, 3.4385692
-3.4385692, 3.4385692, 3.4385692, 3.4385692, 3.4385692, 3.4385692, 3.4385692, 3.4385692, 3.4385692, 3.4385692, 3.4385692, 3.4385692, 3.4385692, 4.6403486, 4.6403486, 4.6403486
-4.6403486, 4.6403486, 4.6403486, 4.6403486, 4.6403486, 4.6403486, 4.6403486, 4.6403486, 4.6403486, 4.6403486, 4.6403486, 4.6403486
-
+0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000, 0.000, 0.000, 0.000, 0.000
+0.000, -2.259, -6.423, 0.866, -6.373, 3.260, 0.228, 3.475, -0.986, -0.091, -2.053, -2.161, 0.986, -0.667, -1.335, -1.144
+-0.988, 2.450, -0.783, 0.930, -0.494, -0.744, -1.135, -0.185, -0.864, 2.173, -0.715, -1.646, 1.488, 2.565, 2.685, 2.411
+2.057, -1.281, -0.742, -0.953, -1.529, -1.142, -1.328, -1.374, -1.311, -1.011, -1.281, -2.159, -0.906, 1.915, 2.023, 2.316
+2.860, -0.980, 1.076, -0.988, -0.983, -2.039, -1.289, -0.999, 0.713, -1.011, 2.025, 0.366, 1.548, 1.253, 1.612, 0.951
+2.049, 1.652, -2.061, 3.358, 3.554, -0.096, 0.003, 5.227, 1.397, 8.266, 0.395, -0.846, 3.159, 5.009, 4.247, -2.114
+3.470, 4.427, 1.356, -1.029, 1.766, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000
 ptzero1x1=
--9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514, -9.1481514
--9.1481514, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401, -6.2416401
--6.2416401, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815, -3.8353815
--3.8353815, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967, -3.3290967
--3.3290967, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141, 3.9013141
-
+-4.100,-3.863,-4.900,-6.566,-5.828,-9.011,-5.034,-7.786,-8.327,-5.643,-8.956,-6.258, -5.092, -5.224, -2.088, -5.075
+-2.166, -2.446, -7.240, -5.218, -6.554, -5.572, -6.060, -4.970, -5.041, -5.193, -7.812, -5.146, -4.917, -5.068, -4.150, 0.241
+-1.206, 1.889, 1.826, -1.596, -1.040, -3.413, -1.433, -4.678, -3.043, -4.870, -4.915, -2.659, -1.679, -1.041, -0.771, -1.072
+-0.514, 2.935, 2.008, 2.193, 0.712, 2.606, 2.235, -1.133, -1.551, -1.184, -2.828, -1.533, -1.522, 2.317, 3.103, 1.915
+6.329, 6.605, 7.082, 9.057, 7.130, 9.753, 9.218, 10.309, 10.729, 10.147, 5.405, 5.596, 5.891, 2.166, 2.515
 ptzero1x2=
-3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923, 3.1489923
-3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498, 3.0414498
--2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186, -2.1861186
-0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274, 0.7342274
-3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818, 3.9505818
-
+-0.054,6.495,3.634,3.292,2.924,-1.446,-3.039,-0.984,-6.139,0.263,1.217,2.159, -0.675, -1.141, -0.846, 3.090
+-0.284, 1.032, -1.310, -3.983, -2.561, -2.450, 0.766, -2.104, -3.275, -1.230, -0.012, -1.139, -0.020, 3.428, 3.339, -4.780
+-1.179, -2.102, -4.915, -5.032, -5.952, -6.334, -8.801, -8.209, -9.309, -8.443, -6.099, -2.252, -4.817, -4.801, -4.513, -2.579
+-0.571, -4.221, -1.441, -5.170, -0.800, -0.934, -1.825, -4.839, -5.031, -5.267, -4.813, -6.362, -4.978, -4.699, -4.462, -0.623
+2.132, 2.098, -0.626, 2.738, 3.790, 2.873, 2.680, 2.984, 2.754, 2.277, -1.814, -1.704, -1.497, -0.875, -2.401
 ptzero1v1=
-2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196, 2.7230196
-0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673, 0.4787673
-1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379, 1.2140379
-5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516, 5.2155516
--3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404, -3.3917404
-0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137, 0.7047137
-3.6539902, 3.6539902, 3.6539902, 3.6539902, 3.6539902, 3.6539902, 3.6539902, 3.6539902, 3.6539902, 3.6539902, 3.6539902, 3.6539902
-
+0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,-10.341,-12.552,-0.041, 5.189, 6.299, -3.498, 4.101
+0.505, -1.908, 1.121, 2.678, 2.466, -5.564, 1.003, 1.987, 3.703, 0.045, 2.090, -1.412, -2.056, 0.135, 1.815, -1.742
+1.419, -1.178, -1.702, 1.961, -1.917, -1.774, 0.281, -1.629, 2.937, 1.794, 2.601, 3.408, 2.533, 3.120, 2.624, 5.775
+7.080, 6.698, 6.144, 6.784, 6.305, 5.765, 5.959, 6.806, 7.314, 6.504, 6.804, 6.410, 6.020, 6.503, 5.746, -5.496
+-5.735, -8.333, -5.148, -3.950, -5.182, -7.476, -2.617, -5.788, -5.493, -2.629, -5.362, -4.541, -4.430, -1.447, -5.489, -1.323
+0.520, -2.334, -1.785, 0.915, -0.521, 1.826, -1.457, 0.657, 1.843, 2.513, 2.258, 3.729, 4.797, 1.715, 2.997, 1.746
+0.643, 7.899, 2.017, 5.716, 11.513, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000
 ptzero1v2=
-1.9417105, 1.9417105, 1.9417105, 1.9417105, 1.9417105, 1.9417105, 1.9417105, 1.9417105, 1.9417105, 1.9417105, 1.9417105, 1.9417105, 2.8306079, 2.8306079, 2.8306079
-2.8306079, 2.8306079, 2.8306079, 2.8306079, 2.8306079, 2.8306079, 2.8306079, 2.8306079, 2.8306079, 2.8306079, 2.8306079, 2.8306079, 2.8306079, 5.0124527, 5.0124527, 5.0124527
-5.0124527, 5.0124527, 5.0124527, 5.0124527, 5.0124527, 5.0124527, 5.0124527, 5.0124527, 5.0124527, 5.0124527, 5.0124527, 5.0124527, 5.0124527, -1.6862142, -1.6862142, -1.6862142
--1.6862142, -1.6862142, -1.6862142, -1.6862142, -1.6862142, -1.6862142, -1.6862142, -1.6862142, -1.6862142, -1.6862142, -1.6862142, -1.6862142, -1.6862142, -1.5509324, -1.5509324, -1.5509324
--1.5509324, -1.5509324, -1.5509324, -1.5509324, -1.5509324, -1.5509324, -1.5509324, -1.5509324, -1.5509324, -1.5509324, -1.5509324, -1.5509324, -1.5509324, 2.7159725, 2.7159725, 2.7159725
-2.7159725, 2.7159725, 2.7159725, 2.7159725, 2.7159725, 2.7159725, 2.7159725, 2.7159725, 2.7159725, 2.7159725, 2.7159725, 2.7159725, 2.7159725, 4.8231869, 4.8231869, 4.8231869
-4.8231869, 4.8231869, 4.8231869, 4.8231869, 4.8231869, 4.8231869, 4.8231869, 4.8231869, 4.8231869, 4.8231869, 4.8231869, 4.8231869
-
+0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,-11.494,-7.211,-0.034, -0.627, 2.679, 2.464, 2.600
+3.408, 3.271, 2.324, -1.483, 1.714, -0.356, 2.042, 3.849, -3.658, -0.373, 2.503, -0.316, 3.062, 5.337, 2.624, 1.666
+6.537, 3.899, 5.882, 4.294, 3.847, 1.913, 3.168, 2.690, 3.548, 7.212, 6.422, 2.220, -4.340, -5.046, -1.268, -3.064
+-1.663, -2.421, -2.072, -5.132, -4.890, -2.886, -5.145, -2.192, -4.648, -4.866, -1.935, -3.456, -4.907, -4.300, -4.372, -4.803
+-4.935, -1.951, -5.922, -4.586, -3.678, -4.533, -5.104, -4.386, -4.944, -4.571, -4.442, -1.702, -1.503, -1.794, 2.135, -0.690
+1.184, 2.213, 2.331, 4.006, 2.482, 1.780, 3.119, -0.157, 4.467, 2.560, 1.708, 2.754, -0.304, 4.914, 3.849, 8.307
+1.586, 8.185, -1.600, 7.045, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000
 ptzero2v2=
-0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, -4.0823132, -4.0823132, -4.0823132
--4.0823132, -4.0823132, -4.0823132, -4.0823132, -4.0823132, -4.0823132, -4.0823132, -4.0823132, -4.0823132, -4.0823132, -4.0823132, -4.0823132, -4.0823132, 2.9209175, 2.9209175, 2.9209175
-2.9209175, 2.9209175, 2.9209175, 2.9209175, 2.9209175, 2.9209175, 2.9209175, 2.9209175, 2.9209175, 2.9209175, 2.9209175, 2.9209175, 2.9209175, 3.1686963, 3.1686963, 3.1686963
-3.1686963, 3.1686963, 3.1686963, 3.1686963, 3.1686963, 3.1686963, 3.1686963, 3.1686963, 3.1686963, 3.1686963, 3.1686963, 3.1686963, 3.1686963, -3.3746987, -3.3746987, -3.3746987
--3.3746987, -3.3746987, -3.3746987, -3.3746987, -3.3746987, -3.3746987, -3.3746987, -3.3746987, -3.3746987, -3.3746987, -3.3746987, -3.3746987, -3.3746987, -1.7125143, -1.7125143, -1.7125143
--1.7125143, -1.7125143, -1.7125143, -1.7125143, -1.7125143, -1.7125143, -1.7125143, -1.7125143, -1.7125143, -1.7125143, -1.7125143, -1.7125143, -1.7125143, -1.5667796, -1.5667796, -1.5667796
--1.5667796, -1.5667796, -1.5667796, -1.5667796, -1.5667796, -1.5667796, -1.5667796, -1.5667796, -1.5667796, -1.5667796, -1.5667796, -1.5667796
-
+0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000, 0.000, -5.411, -1.731, -3.606
+1.443, 1.881, 0.041, -4.095, 1.264, -2.293, 2.071, -3.527, 0.572, 1.145, -6.359, -1.203, 5.445, 10.642, 5.838, 5.287
+6.655, 3.184, 8.361, 5.883, 6.061, 4.524, 2.691, 6.552, 3.021, 5.879, 2.375, 3.164, 6.280, 6.164, 5.666, 6.755
+6.645, 6.457, 6.699, 6.170, 5.950, 6.264, 2.720, 6.361, 6.641, 6.320, 2.764, 6.278, -1.478, -1.164, -1.837, -1.754
+-1.052, -1.878, -2.219, -1.587, -3.191, -2.260, -5.202, -4.525, -2.050, -1.333, 0.886, -2.618, -1.560, -0.690, 1.936, -1.629
+-2.073, -0.183, 1.784, 2.129, 2.055, 0.814, 1.494, -1.156, -0.388, -1.481, 2.086, -2.859, 4.009, -0.725, -0.763, -1.119
+2.252, -0.919, -1.008, -0.193, 1.821, 1.667, 2.035, -3.199, 3.255, 9.044, 0.000
 ptzero2v1=
-3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896, 3.4548896
--3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468, -3.9597468
-3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243, 3.4164243
-3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274, 3.7355274
--2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547, -2.0650547
--2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959, -2.9003959
--0.0782946, -0.0782946, -0.0782946, -0.0782946, -0.0782946, -0.0782946, -0.0782946, -0.0782946, -0.0782946, -0.0782946, -0.0782946, -0.0782946
-
+0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000,0.000, 0.000, -0.904, 8.184, -7.572
+-5.314, -7.689, 0.846, -3.625, -1.955, -8.984, -1.496, -5.006, -3.879, -2.740, -1.798, -2.678, -2.002, -2.171, -3.622, 2.854
+7.242, 3.488, 2.421, 3.132, 2.582, 2.150, 3.531, 3.035, 2.429, 5.682, 6.277, 4.023, 3.475, 5.642, 5.130, 6.548
+15.273, 6.371, 3.119, 6.425, 6.286, 6.229, 6.090, 2.728, 6.531, 6.805, 3.372, 6.168, 6.342, 3.058, 6.619, 0.969
+1.982, -1.019, -1.608, -0.664, -1.892, -0.929, -1.313, -1.413, -1.819, 0.000, 0.103, -1.860, 1.485, 0.250, 1.576, -1.899
+0.385, -2.457, -1.562, -1.442, -1.845, -2.866, -3.172, -1.167, -2.214, -2.223, 5.055, -0.695, -4.876, -6.590, 2.167, 0.130
+2.710, 2.687, -0.306, 2.288, 0.013, 1.922, -0.973, 2.878, 5.169, 4.275, 0.000
 ptzero2x2=
--5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139, -5.1079139
--6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936, -6.3685936
--3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037, -3.4858037
--2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431, -2.8346431
-0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300, 0.8723300
-
+-4.559,-3.178,-4.888,-4.677,-5.489,-5.792,-1.467,-1.747,-1.469,0.163,2.311,-1.721, 1.002, 1.895, -5.946, -9.063
+-8.602, -5.900, -6.015, -5.249, -5.762, -5.250, -5.439, -8.507, -5.577, -3.931, -4.176, -1.581, 0.552, -0.905, -1.879, -1.867
+2.846, 1.787, 2.456, 3.017, -1.747, -4.832, -0.761, -2.882, -4.517, -1.409, 0.000, -0.497, -1.885, -0.782, -1.349, 2.647
+3.111, 2.211, 5.915, 3.830, 2.812, 2.738, 2.139, 2.035, 2.642, -1.135, -2.111, -1.772, -1.293, -1.361, -5.522, 2.853
+2.935, 2.492, 1.651, 3.103, 3.489, 5.792, 5.717, 3.037, 5.862, 6.673, 9.245, 9.401, 7.316, 7.497, 9.499
 ptzero2x1=
--3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482, -3.6266482
--3.6266482, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425, -4.9649425
--4.9649425, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056, -6.5772056
--6.5772056, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077, -7.5255077
--7.5255077, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270, 1.7846270
-
+-2.216,-1.267,-2.240,-4.060,-1.156,-2.173,1.824,-1.343,1.661,-0.110,-1.303,0.401, 2.048, 1.581, 2.781, -1.031
+-4.712, -5.003, -4.612, -2.864, -5.192, -4.253, -5.211, -1.950, -0.743, -0.983, -0.777, -1.464, 1.795, 1.697, 5.364, 2.630
+-2.834, 0.332, 2.143, -1.714, -2.364, -5.359, -5.379, -5.389, -5.522, -5.700, -5.753, -4.270, -5.461, -5.184, -5.052, -5.263
+-5.572, -1.997, -0.972, -1.853, -4.781, -1.685, -5.385, -5.325, -4.684, -1.547, -3.355, -9.483, -6.114, -5.379, -9.109, -5.059
+2.542, 6.458, 1.813, 3.042, 6.244, 6.822, 6.729, 6.696, 7.682, 11.285, 6.999, 10.294, 10.066, 10.425, 10.184
 ptzero2u2=
-1.6466604, 1.6466604, 1.6466604, 1.6466604, 1.6466604, 1.6466604, 1.6466604, 1.6466604, 1.6466604, 1.6466604, 1.6466604, 1.6466604, -0.3276615, -0.3276615, -0.3276615
--0.3276615, -0.3276615, -0.3276615, -0.3276615, -0.3276615, -0.3276615, -0.3276615, -0.3276615, -0.3276615, -0.3276615, -0.3276615, -0.3276615, -0.3276615, 2.1571865, 2.1571865, 2.1571865
-2.1571865, 2.1571865, 2.1571865, 2.1571865, 2.1571865, 2.1571865, 2.1571865, 2.1571865, 2.1571865, 2.1571865, 2.1571865, 2.1571865, 2.1571865, -3.6805938, -3.6805938, -3.6805938
--3.6805938, -3.6805938, -3.6805938, -3.6805938, -3.6805938, -3.6805938, -3.6805938, -3.6805938, -3.6805938, -3.6805938, -3.6805938, -3.6805938, -3.6805938, -1.5710936, -1.5710936, -1.5710936
--1.5710936, -1.5710936, -1.5710936, -1.5710936, -1.5710936, -1.5710936, -1.5710936, -1.5710936, -1.5710936, -1.5710936, -1.5710936, -1.5710936, -1.5710936, -0.7569959, -0.7569959, -0.7569959
--0.7569959, -0.7569959, -0.7569959, -0.7569959, -0.7569959, -0.7569959, -0.7569959, -0.7569959, -0.7569959, -0.7569959, -0.7569959, -0.7569959, -0.7569959, 9.1694019, 9.1694019, 9.1694019
-9.1694019, 9.1694019, 9.1694019, 9.1694019, 9.1694019, 9.1694019, 9.1694019, 9.1694019, 9.1694019, 9.1694019, 9.1694019, 9.1694019
-
+0.000,0.000,0.000,0.405,-1.487,-6.910,2.197,1.827,-1.380,-1.979,2.264,-0.928, 0.935, 1.441, 2.840, 1.403
+-2.756, 1.790, 4.070, 1.049, -1.328, -2.054, -0.090, -1.564, -1.539, -0.924, 1.952, 0.746, 0.399, 2.088, -0.792, 7.827
+-1.372, -1.014, 0.569, -0.177, 2.453, -1.166, 2.156, -1.408, 1.645, 1.381, 2.878, 5.599, -2.666, -1.567, -2.774, -2.003
+-5.376, -1.820, -4.650, -5.297, -5.144, -5.963, -5.051, -5.532, -4.563, -4.362, -5.216, -2.565, -1.218, -1.364, -1.859, -1.300
+-0.992, -1.824, -0.538, -1.168, -2.146, -2.601, -1.508, -1.481, 0.541, -0.428, -1.309, -4.917, 0.433, 1.950, 1.556, -1.300
+-1.593, -1.856, 1.892, 1.684, -0.575, -0.236, 3.766, -1.318, -0.228, 0.032, 2.795, 2.527, 10.708, 11.407, 11.900, 12.005
+4.221, 5.191, 12.781, 17.758, -1.475, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000
 ptzero2u1=
-2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888, 2.5827888
--4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750, -4.3624750
-0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103, 0.0707103
-3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520, 3.0557520
--3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296, -3.1667296
--1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877, -1.7278877
-7.3810489, 7.3810489, 7.3810489, 7.3810489, 7.3810489, 7.3810489, 7.3810489, 7.3810489, 7.3810489, 7.3810489, 7.3810489, 7.3810489
-
+0.000,0.000,0.000,0.000,3.246,2.838,3.743,-0.838,1.993,1.961,2.057,2.645, 8.569, 3.941, 4.707, -5.476
+-4.864, -4.584, -2.936, -2.783, -1.279, 0.203, -4.515, -6.547, -3.904, -3.416, -2.928, -3.791, -3.295, 1.113, -3.542, 0.590
+2.172, -0.502, 1.490, 1.708, -1.400, 1.829, -1.859, 2.870, 2.363, 2.226, 2.136, 2.076, 5.586, 5.368, 5.663, 5.562
+3.044, 5.675, 4.096, 3.263, 2.668, 3.230, 2.614, 5.988, 6.117, 6.291, 6.631, 5.792, 6.271, 6.099, 7.125, -5.071
+-2.270, -1.127, -0.956, -4.889, -5.398, -2.239, -1.186, -0.990, -2.110, -1.802, 0.020, -1.039, 3.813, 6.148, 1.914, -0.971
+-0.340, -0.983, -1.516, 2.058, 2.030, -1.306, -0.926, 2.455, 1.172, -0.657, -0.300, 6.149, 2.051, 4.677, 4.324, 10.896
+15.255, 13.527, 11.018, -135.385, 13.268, 19.795, 0.000, 0.000, 0.000, 0.000, 0.000
diff --git a/PARAM/SHMS/DC/pdriftmap.param b/PARAM/SHMS/DC/pdriftmap.param
index 7799e910..3236157a 100644
--- a/PARAM/SHMS/DC/pdriftmap.param
+++ b/PARAM/SHMS/DC/pdriftmap.param
@@ -5,171 +5,171 @@ pdriftbins = 138
 pdrift1stbin=0
 ; bin size in ns
 pdriftbinsz=2
-pwc1u1fract=0.0012918,0.00393,0.00820,0.01502,0.02425,0.03662,0.05029,0.06599,0.08300
-0.10141,0.12121,0.14156,0.16225,0.18309,0.20438,0.22640,0.24785,0.26943,0.29094
-0.31248,0.33413,0.35621,0.37791,0.39938,0.42126,0.44260,0.46433,0.48604,0.50716
-0.52843,0.54937,0.57034,0.59094,0.61141,0.63218,0.65184,0.67234,0.69158,0.71067
-0.72970,0.74780,0.76594,0.78341,0.80072,0.81751,0.83289,0.84881,0.86334,0.87660
-0.88864,0.89945,0.90916,0.91760,0.92493,0.93150,0.93681,0.94170,0.94612,0.94989
-0.95348,0.95672,0.95983,0.96253,0.96502,0.96715,0.96934,0.97131,0.97314,0.97484
-0.97664,0.97817,0.97953,0.98087,0.98209,0.98319,0.98425,0.98510,0.98594,0.98672
-0.98745,0.98816,0.98877,0.98923,0.98972,0.99012,0.99053,0.99090,0.99118,0.99152
-0.99175,0.99202,0.99227,0.99254,0.99287,0.99313,0.99338,0.99362,0.99382,0.99400
-0.99416,0.99435,0.99458,0.99475,0.99495,0.99517,0.99541,0.99560,0.99581,0.99598
-0.99620,0.99641,0.99664,0.99689,0.99708,0.99726,0.99744,0.99758,0.99779,0.99793
-0.99813,0.99833,0.99855,0.99871,0.99890,0.99913,0.99933,0.99953,0.99974,0.99992
-0.99995,0.99996,0.99997,0.99997,0.99998,0.99998,0.99999,1.00000,1.00000
-pwc1u2fract=0.00084,0.00261,0.00608,0.01174,0.01994,0.03047,0.04371,0.05919,0.07632
-0.09403,0.11382,0.13511,0.15574,0.17712,0.19842,0.22085,0.24207,0.26366,0.28503
-0.30704,0.32860,0.35011,0.37177,0.39324,0.41472,0.43579,0.45711,0.47832,0.49939
-0.52040,0.54102,0.56149,0.58188,0.60190,0.62203,0.64158,0.66105,0.68067,0.69954
-0.71840,0.73671,0.75494,0.77257,0.78944,0.80656,0.82260,0.83811,0.85270,0.86690
-0.88010,0.89219,0.90275,0.91204,0.91989,0.92699,0.93297,0.93856,0.94325,0.94756
-0.95127,0.95465,0.95768,0.96047,0.96318,0.96575,0.96794,0.96997,0.97195,0.97370
-0.97553,0.97699,0.97841,0.97983,0.98107,0.98213,0.98313,0.98408,0.98497,0.98583
-0.98656,0.98720,0.98779,0.98842,0.98894,0.98946,0.98986,0.99028,0.99067,0.99108
-0.99136,0.99164,0.99196,0.99221,0.99252,0.99274,0.99298,0.99319,0.99341,0.99363
-0.99385,0.99404,0.99420,0.99439,0.99459,0.99484,0.99502,0.99528,0.99553,0.99577
-0.99602,0.99621,0.99645,0.99670,0.99693,0.99715,0.99736,0.99758,0.99772,0.99791
-0.99811,0.99834,0.99852,0.99870,0.99891,0.99913,0.99936,0.99953,0.99978,0.99992
-0.99997,0.99998,0.99998,0.99999,1.00000,1.00000,1.00000,1.00000,1.00000
-pwc1x1fract=0.00199,0.00555,0.01149,0.02037,0.03148,0.04552,0.06196,0.07957,0.09820
-0.11791,0.13888,0.15965,0.18097,0.20207,0.22416,0.24500,0.26649,0.28752,0.30895
-0.33010,0.35174,0.37289,0.39431,0.41543,0.43666,0.45765,0.47858,0.49959,0.52043
-0.54111,0.56187,0.58255,0.60258,0.62249,0.64258,0.66215,0.68081,0.69942,0.71745
-0.73556,0.75320,0.77041,0.78783,0.80499,0.82153,0.83660,0.85099,0.86446,0.87735
-0.88915,0.89926,0.90829,0.91603,0.92285,0.92891,0.93453,0.93921,0.94374,0.94759
-0.95092,0.95406,0.95697,0.95973,0.96233,0.96488,0.96711,0.96945,0.97132,0.97311
-0.97487,0.97652,0.97795,0.97951,0.98083,0.98180,0.98295,0.98405,0.98483,0.98569
-0.98647,0.98717,0.98776,0.98832,0.98883,0.98935,0.98988,0.99022,0.99063,0.99104
-0.99139,0.99168,0.99202,0.99228,0.99254,0.99280,0.99302,0.99332,0.99354,0.99374
-0.99394,0.99416,0.99439,0.99455,0.99479,0.99501,0.99523,0.99550,0.99573,0.99595
-0.99619,0.99642,0.99660,0.99674,0.99697,0.99712,0.99733,0.99754,0.99775,0.99800
-0.99818,0.99841,0.99865,0.99887,0.99907,0.99926,0.99940,0.99957,0.99974,0.99983
-0.99990,0.99996,0.99998,0.99999,0.99999,0.99999,0.99999,1.00000,1.00000
-pwc1x2fract=0.00104,0.00304,0.00682,0.01294,0.02230,0.03465,0.04998,0.06633,0.08490
-0.10426,0.12405,0.14473,0.16554,0.18667,0.20787,0.22900,0.24992,0.27112,0.29198
-0.31371,0.33521,0.35640,0.37773,0.39895,0.42012,0.44129,0.46246,0.48265,0.50281
-0.52309,0.54299,0.56371,0.58384,0.60377,0.62314,0.64259,0.66198,0.68112,0.69986
-0.71785,0.73552,0.75291,0.77037,0.78734,0.80354,0.81980,0.83515,0.84997,0.86352
-0.87618,0.88764,0.89807,0.90743,0.91548,0.92242,0.92847,0.93374,0.93838,0.94251
-0.94621,0.94948,0.95253,0.95545,0.95840,0.96102,0.96352,0.96584,0.96805,0.97004
-0.97189,0.97352,0.97526,0.97678,0.97830,0.97968,0.98094,0.98199,0.98319,0.98417
-0.98511,0.98601,0.98668,0.98737,0.98799,0.98855,0.98911,0.98956,0.99001,0.99042
-0.99081,0.99115,0.99145,0.99175,0.99206,0.99230,0.99261,0.99285,0.99313,0.99338
-0.99367,0.99392,0.99418,0.99443,0.99463,0.99488,0.99511,0.99533,0.99556,0.99577
-0.99602,0.99628,0.99647,0.99668,0.99689,0.99712,0.99728,0.99747,0.99774,0.99791
-0.99811,0.99827,0.99849,0.99868,0.99884,0.99902,0.99920,0.99942,0.99963,0.99978
-0.99989,0.99993,0.99996,0.99998,0.99998,0.99998,0.99999,0.99999,1.00000
-pwc1v1fract=0.00124,0.00390,0.00803,0.01445,0.02371,0.03499,0.04899,0.06461,0.08196
-0.10025,0.11984,0.14093,0.16231,0.18420,0.20581,0.22720,0.24977,0.27193,0.29343
-0.31529,0.33767,0.35914,0.38069,0.40209,0.42362,0.44576,0.46682,0.48812,0.50932
-0.53014,0.55069,0.57083,0.59087,0.61156,0.63185,0.65123,0.67066,0.68930,0.70817
-0.72701,0.74476,0.76252,0.77952,0.79657,0.81329,0.82890,0.84484,0.85921,0.87279
-0.88484,0.89560,0.90519,0.91376,0.92125,0.92775,0.93332,0.93841,0.94289,0.94685
-0.95067,0.95400,0.95696,0.95975,0.96253,0.96492,0.96731,0.96945,0.97143,0.97323
-0.97496,0.97646,0.97794,0.97928,0.98054,0.98181,0.98279,0.98391,0.98489,0.98575
-0.98644,0.98716,0.98778,0.98845,0.98898,0.98945,0.98995,0.99032,0.99062,0.99092
-0.99121,0.99150,0.99179,0.99206,0.99236,0.99261,0.99279,0.99305,0.99329,0.99355
-0.99380,0.99399,0.99422,0.99443,0.99468,0.99486,0.99508,0.99530,0.99550,0.99574
-0.99596,0.99618,0.99638,0.99659,0.99676,0.99695,0.99718,0.99740,0.99757,0.99782
-0.99805,0.99829,0.99856,0.99878,0.99899,0.99920,0.99941,0.99962,0.99977,0.99987
-0.99993,0.99996,0.99997,0.99998,0.99998,0.99998,0.99999,1.00000,1.00000
-pwc1v2fract=0.00120,0.00356,0.00791,0.01501,0.02503,0.03780,0.05265,0.06963,0.08795
-0.10693,0.12710,0.14864,0.16955,0.19124,0.21313,0.23459,0.25639,0.27811,0.29960
-0.32117,0.34208,0.36349,0.38478,0.40636,0.42740,0.44819,0.46989,0.49108,0.51235
-0.53328,0.55372,0.57372,0.59371,0.61376,0.63356,0.65358,0.67295,0.69235,0.71095
-0.72951,0.74771,0.76565,0.78345,0.80015,0.81689,0.83271,0.84776,0.86191,0.87554
-0.88784,0.89884,0.90833,0.91679,0.92419,0.93048,0.93578,0.94064,0.94494,0.94863
-0.95195,0.95525,0.95794,0.96073,0.96320,0.96563,0.96789,0.96985,0.97169,0.97331
-0.97497,0.97652,0.97787,0.97923,0.98036,0.98162,0.98279,0.98379,0.98473,0.98550
-0.98626,0.98690,0.98746,0.98808,0.98860,0.98896,0.98946,0.98987,0.99023,0.99063
-0.99096,0.99129,0.99161,0.99191,0.99218,0.99248,0.99277,0.99302,0.99323,0.99348
-0.99371,0.99390,0.99412,0.99438,0.99461,0.99479,0.99502,0.99521,0.99538,0.99559
-0.99579,0.99600,0.99623,0.99646,0.99666,0.99690,0.99711,0.99729,0.99753,0.99773
-0.99791,0.99813,0.99840,0.99864,0.99888,0.99909,0.99924,0.99948,0.99965,0.99977
-0.99988,0.99999,0.99999,0.99999,0.99999,1.00000,1.00000,1.00000,1.00000
-pwc2v2fract=0.00116,0.00336,0.00699,0.01286,0.02132,0.03216,0.04616,0.06217,0.07977
-0.09901,0.11949,0.14089,0.16228,0.18440,0.20633,0.22805,0.25032,0.27173,0.29323
-0.31432,0.33572,0.35738,0.37866,0.39962,0.42101,0.44223,0.46310,0.48401,0.50486
-0.52500,0.54535,0.56530,0.58499,0.60417,0.62324,0.64198,0.66008,0.67800,0.69562
-0.71280,0.72981,0.74627,0.76238,0.77790,0.79326,0.80792,0.82222,0.83574,0.84859
-0.86046,0.87176,0.88209,0.89143,0.89989,0.90732,0.91420,0.92004,0.92537,0.93038
-0.93495,0.93933,0.94327,0.94693,0.95011,0.95319,0.95601,0.95872,0.96117,0.96340
-0.96559,0.96757,0.96943,0.97114,0.97280,0.97425,0.97559,0.97695,0.97820,0.97941
-0.98054,0.98155,0.98253,0.98346,0.98414,0.98497,0.98584,0.98650,0.98707,0.98760
-0.98819,0.98877,0.98921,0.98970,0.99016,0.99056,0.99099,0.99142,0.99182,0.99217
-0.99252,0.99283,0.99317,0.99351,0.99380,0.99406,0.99434,0.99464,0.99489,0.99513
-0.99537,0.99563,0.99590,0.99618,0.99641,0.99665,0.99686,0.99711,0.99732,0.99755
-0.99782,0.99811,0.99838,0.99863,0.99890,0.99913,0.99939,0.99960,0.99976,0.99987
-0.99995,0.99997,0.99997,0.99998,0.99999,0.99999,0.99999,0.99999,1.00000
-pwc2v1fract=0.00106,0.00267,0.00550,0.01041,0.01777,0.02803,0.04043,0.05556,0.07182
-0.09063,0.11136,0.13254,0.15353,0.17503,0.19729,0.21878,0.24043,0.26245,0.28389
-0.30556,0.32716,0.34865,0.36989,0.39129,0.41230,0.43366,0.45482,0.47541,0.49591
-0.51616,0.53686,0.55710,0.57690,0.59662,0.61616,0.63480,0.65363,0.67195,0.69012
-0.70795,0.72464,0.74095,0.75724,0.77374,0.78943,0.80408,0.81842,0.83256,0.84582
-0.85808,0.86985,0.88080,0.89065,0.89953,0.90748,0.91436,0.92049,0.92606,0.93111
-0.93550,0.93951,0.94356,0.94726,0.95046,0.95363,0.95664,0.95937,0.96176,0.96411
-0.96636,0.96829,0.97015,0.97185,0.97350,0.97492,0.97638,0.97774,0.97892,0.98009
-0.98117,0.98227,0.98316,0.98410,0.98496,0.98573,0.98648,0.98718,0.98785,0.98835
-0.98890,0.98935,0.98984,0.99022,0.99070,0.99103,0.99135,0.99173,0.99203,0.99239
-0.99273,0.99307,0.99334,0.99363,0.99392,0.99419,0.99445,0.99471,0.99494,0.99521
-0.99547,0.99575,0.99601,0.99622,0.99648,0.99676,0.99698,0.99721,0.99747,0.99766
-0.99786,0.99812,0.99836,0.99858,0.99884,0.99909,0.99932,0.99957,0.99970,0.99984
-0.99995,0.99997,0.99997,0.99997,0.99997,0.99998,0.99999,1.00000,1.00000
-pwc2x2fract=0.00170,0.00463,0.00947,0.01681,0.02696,0.03988,0.05573,0.07310,0.09202
-0.11181,0.13244,0.15300,0.17465,0.19592,0.21768,0.23898,0.26022,0.28119,0.30216
-0.32338,0.34460,0.36564,0.38665,0.40722,0.42847,0.44916,0.46986,0.49049,0.51145
-0.53172,0.55170,0.57151,0.59052,0.60980,0.62823,0.64707,0.66570,0.68388,0.70135
-0.71885,0.73548,0.75232,0.76833,0.78451,0.80004,0.81480,0.82904,0.84284,0.85647
-0.86841,0.87929,0.88996,0.89896,0.90691,0.91380,0.91993,0.92528,0.93005,0.93439
-0.93824,0.94203,0.94541,0.94851,0.95147,0.95439,0.95701,0.95941,0.96197,0.96408
-0.96619,0.96829,0.97024,0.97201,0.97366,0.97537,0.97677,0.97820,0.97941,0.98069
-0.98173,0.98281,0.98385,0.98468,0.98548,0.98625,0.98699,0.98762,0.98819,0.98878
-0.98927,0.98977,0.99022,0.99062,0.99107,0.99145,0.99186,0.99218,0.99254,0.99281
-0.99311,0.99336,0.99365,0.99393,0.99420,0.99448,0.99471,0.99498,0.99517,0.99543
-0.99564,0.99584,0.99610,0.99635,0.99656,0.99679,0.99704,0.99732,0.99755,0.99779
-0.99801,0.99824,0.99846,0.99867,0.99888,0.99909,0.99933,0.99955,0.99976,0.99992
-0.99995,0.99998,0.99998,0.99999,0.99999,0.99999,0.99999,1.00000,1.00000
-pwc2x1fract=0.00184,0.00503,0.01008,0.01807,0.02845,0.04114,0.05585,0.07310,0.09138
-0.11098,0.13140,0.15219,0.17311,0.19406,0.21464,0.23548,0.25666,0.27740,0.29860
-0.31939,0.34055,0.36146,0.38242,0.40345,0.42427,0.44451,0.46480,0.48584,0.50670
-0.52676,0.54726,0.56755,0.58735,0.60709,0.62627,0.64521,0.66372,0.68223,0.70029
-0.71820,0.73547,0.75264,0.76901,0.78511,0.80101,0.81633,0.83086,0.84462,0.85762
-0.86981,0.88128,0.89147,0.90046,0.90807,0.91518,0.92136,0.92660,0.93143,0.93593
-0.94017,0.94390,0.94736,0.95058,0.95351,0.95631,0.95884,0.96135,0.96375,0.96597
-0.96794,0.96999,0.97189,0.97355,0.97517,0.97672,0.97804,0.97947,0.98057,0.98176
-0.98269,0.98357,0.98451,0.98524,0.98595,0.98666,0.98729,0.98795,0.98846,0.98907
-0.98956,0.99005,0.99053,0.99095,0.99133,0.99167,0.99199,0.99236,0.99270,0.99296
-0.99319,0.99348,0.99374,0.99405,0.99429,0.99457,0.99483,0.99503,0.99530,0.99556
-0.99577,0.99600,0.99626,0.99644,0.99666,0.99692,0.99717,0.99739,0.99759,0.99778
-0.99803,0.99828,0.99853,0.99880,0.99897,0.99924,0.99949,0.99966,0.99975,0.99986
-0.99995,0.99997,0.99997,0.99998,0.99998,0.99999,0.99999,0.99999,1.00000
-pwc2u2fract=0.00168,0.00467,0.00944,0.01694,0.02662,0.03880,0.05370,0.07025,0.08822
-0.10728,0.12710,0.14783,0.16849,0.18970,0.21092,0.23261,0.25397,0.27587,0.29758
-0.31901,0.34025,0.36146,0.38317,0.40436,0.42527,0.44598,0.46709,0.48780,0.50889
-0.52946,0.54962,0.57009,0.59034,0.61033,0.62965,0.64944,0.66862,0.68709,0.70499
-0.72295,0.74034,0.75733,0.77377,0.79026,0.80563,0.82104,0.83583,0.84956,0.86274
-0.87481,0.88577,0.89588,0.90421,0.91169,0.91822,0.92403,0.92915,0.93380,0.93801
-0.94188,0.94561,0.94912,0.95232,0.95522,0.95788,0.96051,0.96301,0.96531,0.96754
-0.96956,0.97132,0.97321,0.97478,0.97631,0.97764,0.97889,0.98017,0.98122,0.98228
-0.98333,0.98416,0.98495,0.98581,0.98648,0.98711,0.98777,0.98834,0.98892,0.98938
-0.98975,0.99014,0.99049,0.99081,0.99116,0.99142,0.99174,0.99205,0.99238,0.99269
-0.99299,0.99324,0.99349,0.99381,0.99411,0.99433,0.99459,0.99479,0.99502,0.99522
-0.99544,0.99566,0.99595,0.99618,0.99640,0.99662,0.99688,0.99713,0.99734,0.99756
-0.99782,0.99814,0.99833,0.99851,0.99877,0.99899,0.99923,0.99946,0.99964,0.99981
-0.99993,0.99997,0.99997,0.99998,0.99998,0.99998,0.99998,0.99999,1.00000
-pwc2u1fract=0.00138,0.00386,0.00848,0.01523,0.02445,0.03599,0.04963,0.06559,0.08327
-0.10229,0.12169,0.14196,0.16239,0.18355,0.20505,0.22665,0.24828,0.27077,0.29198
-0.31368,0.33562,0.35694,0.37839,0.39925,0.42009,0.44114,0.46206,0.48350,0.50421
-0.52463,0.54510,0.56520,0.58545,0.60526,0.62487,0.64472,0.66373,0.68222,0.70071
-0.71896,0.73686,0.75407,0.77085,0.78745,0.80368,0.81939,0.83448,0.84887,0.86212
-0.87423,0.88547,0.89526,0.90402,0.91192,0.91873,0.92510,0.93040,0.93531,0.93988
-0.94404,0.94779,0.95112,0.95426,0.95731,0.96017,0.96273,0.96526,0.96752,0.96950
-0.97161,0.97338,0.97522,0.97671,0.97806,0.97941,0.98055,0.98180,0.98277,0.98377
-0.98458,0.98534,0.98614,0.98682,0.98746,0.98802,0.98854,0.98903,0.98950,0.98991
-0.99042,0.99082,0.99117,0.99148,0.99184,0.99215,0.99246,0.99276,0.99308,0.99335
-0.99364,0.99389,0.99410,0.99436,0.99460,0.99484,0.99507,0.99525,0.99548,0.99573
-0.99589,0.99613,0.99639,0.99659,0.99680,0.99699,0.99724,0.99746,0.99771,0.99790
-0.99812,0.99834,0.99858,0.99880,0.99903,0.99926,0.99946,0.99963,0.99979,0.99986
-0.99995,0.99997,0.99997,0.99997,0.99998,0.99998,0.99999,0.99999,1.00000
+pwc1u1fract=0.0017378,0.00422,0.00765,0.01260,0.01932,0.02812,0.03902,0.05177,0.06669
+0.08263,0.10049,0.11930,0.13923,0.15950,0.18074,0.20249,0.22407,0.24555,0.26803
+0.28965,0.31145,0.33338,0.35490,0.37716,0.39879,0.42088,0.44228,0.46433,0.48559
+0.50713,0.52862,0.54989,0.57105,0.59115,0.61217,0.63224,0.65245,0.67225,0.69196
+0.71088,0.73006,0.74806,0.76552,0.78279,0.79960,0.81632,0.83219,0.84753,0.86146
+0.87474,0.88720,0.89830,0.90840,0.91747,0.92527,0.93193,0.93782,0.94288,0.94753
+0.95143,0.95514,0.95809,0.96102,0.96363,0.96593,0.96813,0.97003,0.97178,0.97357
+0.97528,0.97684,0.97836,0.97961,0.98100,0.98202,0.98307,0.98407,0.98506,0.98602
+0.98679,0.98746,0.98812,0.98868,0.98920,0.98969,0.99010,0.99054,0.99085,0.99117
+0.99150,0.99179,0.99212,0.99232,0.99258,0.99288,0.99310,0.99331,0.99354,0.99379
+0.99399,0.99417,0.99435,0.99457,0.99479,0.99498,0.99515,0.99539,0.99559,0.99581
+0.99602,0.99627,0.99647,0.99667,0.99685,0.99712,0.99730,0.99751,0.99768,0.99784
+0.99800,0.99819,0.99840,0.99861,0.99879,0.99896,0.99917,0.99937,0.99949,0.99961
+0.99976,0.99990,0.99994,0.99999,0.99999,1.00000,1.00000,1.00000,1.00000
+pwc1u2fract=0.00177,0.00440,0.00781,0.01292,0.01970,0.02875,0.03940,0.05267,0.06756
+0.08431,0.10256,0.12179,0.14228,0.16294,0.18442,0.20643,0.22811,0.25017,0.27252
+0.29440,0.31647,0.33873,0.36051,0.38232,0.40408,0.42577,0.44671,0.46822,0.48972
+0.51039,0.53098,0.55188,0.57206,0.59293,0.61329,0.63285,0.65182,0.67127,0.69003
+0.70859,0.72710,0.74516,0.76299,0.78018,0.79735,0.81373,0.82935,0.84419,0.85830
+0.87189,0.88432,0.89597,0.90625,0.91539,0.92306,0.92992,0.93609,0.94149,0.94609
+0.95023,0.95391,0.95734,0.96026,0.96306,0.96540,0.96762,0.96974,0.97152,0.97305
+0.97479,0.97635,0.97763,0.97905,0.98034,0.98157,0.98269,0.98366,0.98458,0.98536
+0.98616,0.98698,0.98756,0.98805,0.98862,0.98914,0.98958,0.99003,0.99043,0.99078
+0.99109,0.99143,0.99175,0.99205,0.99234,0.99255,0.99280,0.99307,0.99328,0.99349
+0.99371,0.99388,0.99409,0.99430,0.99454,0.99474,0.99493,0.99510,0.99533,0.99553
+0.99579,0.99599,0.99625,0.99648,0.99674,0.99701,0.99718,0.99740,0.99759,0.99777
+0.99796,0.99815,0.99834,0.99849,0.99865,0.99891,0.99913,0.99931,0.99946,0.99963
+0.99977,0.99987,0.99994,0.99998,0.99999,1.00000,1.00000,1.00000,1.00000
+pwc1x1fract=0.00180,0.00464,0.00820,0.01325,0.02098,0.03045,0.04231,0.05646,0.07196
+0.09003,0.10897,0.12895,0.14981,0.17129,0.19283,0.21453,0.23628,0.25858,0.28020
+0.30182,0.32320,0.34495,0.36664,0.38861,0.40969,0.43119,0.45271,0.47374,0.49492
+0.51589,0.53706,0.55749,0.57793,0.59822,0.61816,0.63811,0.65723,0.67624,0.69443
+0.71301,0.73116,0.74860,0.76584,0.78273,0.79942,0.81560,0.83119,0.84638,0.86032
+0.87305,0.88494,0.89546,0.90490,0.91357,0.92121,0.92820,0.93404,0.93922,0.94369
+0.94777,0.95126,0.95456,0.95756,0.96027,0.96276,0.96493,0.96720,0.96943,0.97124
+0.97300,0.97469,0.97623,0.97772,0.97913,0.98052,0.98172,0.98273,0.98370,0.98462
+0.98546,0.98625,0.98697,0.98760,0.98824,0.98876,0.98930,0.98973,0.99016,0.99051
+0.99090,0.99121,0.99155,0.99183,0.99209,0.99240,0.99264,0.99286,0.99313,0.99336
+0.99362,0.99389,0.99408,0.99427,0.99447,0.99468,0.99492,0.99516,0.99538,0.99562
+0.99582,0.99601,0.99622,0.99642,0.99659,0.99680,0.99700,0.99722,0.99744,0.99765
+0.99786,0.99808,0.99828,0.99849,0.99867,0.99887,0.99906,0.99922,0.99940,0.99953
+0.99967,0.99979,0.99986,0.99990,0.99996,0.99998,0.99999,0.99999,1.00000
+pwc1x2fract=0.00169,0.00405,0.00738,0.01245,0.01958,0.02873,0.04068,0.05480,0.07137
+0.08913,0.10809,0.12786,0.14864,0.16965,0.19098,0.21228,0.23377,0.25576,0.27678
+0.29868,0.32045,0.34223,0.36362,0.38504,0.40647,0.42789,0.44835,0.46969,0.49054
+0.51111,0.53162,0.55186,0.57194,0.59204,0.61164,0.63127,0.65005,0.66885,0.68779
+0.70576,0.72413,0.74215,0.75942,0.77621,0.79249,0.80839,0.82422,0.83890,0.85304
+0.86663,0.87886,0.89011,0.90029,0.90958,0.91760,0.92478,0.93124,0.93691,0.94144
+0.94551,0.94919,0.95236,0.95519,0.95772,0.96042,0.96272,0.96483,0.96710,0.96937
+0.97106,0.97269,0.97432,0.97614,0.97755,0.97890,0.98022,0.98143,0.98247,0.98338
+0.98437,0.98533,0.98611,0.98693,0.98768,0.98827,0.98873,0.98926,0.98969,0.99014
+0.99053,0.99095,0.99124,0.99163,0.99196,0.99224,0.99250,0.99275,0.99303,0.99327
+0.99349,0.99375,0.99399,0.99426,0.99442,0.99470,0.99494,0.99515,0.99541,0.99561
+0.99584,0.99609,0.99630,0.99657,0.99678,0.99696,0.99711,0.99733,0.99754,0.99770
+0.99788,0.99804,0.99828,0.99849,0.99868,0.99886,0.99904,0.99919,0.99935,0.99950
+0.99964,0.99978,0.99988,0.99994,0.99998,0.99999,1.00000,1.00000,1.00000
+pwc1v1fract=0.00187,0.00456,0.00826,0.01369,0.02145,0.03105,0.04242,0.05584,0.07176
+0.08884,0.10732,0.12725,0.14792,0.16946,0.19162,0.21311,0.23527,0.25760,0.28018
+0.30273,0.32482,0.34658,0.36911,0.39091,0.41314,0.43433,0.45563,0.47722,0.49846
+0.51967,0.54043,0.56097,0.58134,0.60122,0.62132,0.64103,0.66093,0.67963,0.69873
+0.71739,0.73583,0.75340,0.77063,0.78751,0.80415,0.82016,0.83569,0.85019,0.86350
+0.87643,0.88850,0.89908,0.90883,0.91743,0.92469,0.93117,0.93666,0.94172,0.94625
+0.95015,0.95372,0.95685,0.95957,0.96211,0.96446,0.96675,0.96902,0.97085,0.97261
+0.97435,0.97579,0.97727,0.97847,0.97978,0.98104,0.98218,0.98331,0.98428,0.98527
+0.98609,0.98692,0.98757,0.98820,0.98870,0.98918,0.98967,0.99007,0.99047,0.99079
+0.99106,0.99134,0.99162,0.99187,0.99214,0.99240,0.99266,0.99292,0.99312,0.99333
+0.99360,0.99385,0.99404,0.99424,0.99450,0.99469,0.99490,0.99504,0.99532,0.99558
+0.99580,0.99600,0.99622,0.99642,0.99663,0.99684,0.99707,0.99722,0.99741,0.99767
+0.99791,0.99813,0.99838,0.99859,0.99884,0.99905,0.99920,0.99940,0.99956,0.99969
+0.99978,0.99988,0.99991,0.99995,0.99999,1.00000,1.00000,1.00000,1.00000
+pwc1v2fract=0.00177,0.00449,0.00797,0.01283,0.02026,0.02979,0.04130,0.05509,0.07113
+0.08830,0.10623,0.12579,0.14646,0.16806,0.18953,0.21155,0.23327,0.25517,0.27696
+0.29907,0.32061,0.34278,0.36481,0.38671,0.40796,0.42884,0.45050,0.47215,0.49299
+0.51367,0.53488,0.55563,0.57555,0.59562,0.61538,0.63551,0.65489,0.67446,0.69346
+0.71243,0.73059,0.74849,0.76585,0.78324,0.80015,0.81663,0.83213,0.84704,0.86121
+0.87453,0.88736,0.89806,0.90764,0.91635,0.92399,0.93079,0.93680,0.94180,0.94620
+0.95019,0.95370,0.95675,0.95955,0.96211,0.96428,0.96643,0.96837,0.97036,0.97215
+0.97387,0.97537,0.97684,0.97813,0.97948,0.98068,0.98180,0.98284,0.98375,0.98469
+0.98556,0.98619,0.98695,0.98762,0.98812,0.98863,0.98910,0.98955,0.98993,0.99028
+0.99071,0.99102,0.99138,0.99171,0.99193,0.99220,0.99248,0.99276,0.99303,0.99327
+0.99356,0.99381,0.99405,0.99424,0.99444,0.99463,0.99484,0.99500,0.99521,0.99538
+0.99564,0.99584,0.99608,0.99626,0.99647,0.99669,0.99685,0.99704,0.99728,0.99752
+0.99779,0.99800,0.99821,0.99839,0.99855,0.99874,0.99896,0.99912,0.99935,0.99948
+0.99964,0.99972,0.99984,0.99991,0.99996,0.99998,0.99999,1.00000,1.00000
+pwc2v2fract=0.00199,0.00487,0.00855,0.01361,0.02082,0.02974,0.04077,0.05386,0.06965
+0.08703,0.10601,0.12604,0.14744,0.16901,0.19107,0.21310,0.23492,0.25730,0.27900
+0.30033,0.32194,0.34358,0.36555,0.38695,0.40840,0.42988,0.45133,0.47226,0.49317
+0.51393,0.53393,0.55421,0.57407,0.59307,0.61246,0.63148,0.65014,0.66834,0.68610
+0.70328,0.72031,0.73655,0.75263,0.76819,0.78394,0.79911,0.81346,0.82721,0.84019
+0.85279,0.86411,0.87467,0.88487,0.89409,0.90240,0.90982,0.91645,0.92279,0.92809
+0.93288,0.93737,0.94160,0.94555,0.94902,0.95228,0.95510,0.95780,0.96034,0.96258
+0.96484,0.96681,0.96867,0.97039,0.97210,0.97368,0.97509,0.97635,0.97765,0.97882
+0.97999,0.98108,0.98216,0.98319,0.98408,0.98480,0.98550,0.98622,0.98681,0.98740
+0.98801,0.98856,0.98910,0.98951,0.98994,0.99041,0.99083,0.99123,0.99164,0.99200
+0.99233,0.99266,0.99294,0.99322,0.99355,0.99386,0.99416,0.99443,0.99470,0.99500
+0.99521,0.99548,0.99574,0.99598,0.99623,0.99648,0.99672,0.99697,0.99723,0.99747
+0.99770,0.99794,0.99818,0.99840,0.99870,0.99896,0.99915,0.99934,0.99948,0.99964
+0.99978,0.99986,0.99993,0.99996,0.99998,0.99999,1.00000,1.00000,1.00000
+pwc2v1fract=0.00200,0.00479,0.00875,0.01389,0.02117,0.03014,0.04155,0.05532,0.07088
+0.08843,0.10721,0.12697,0.14792,0.16955,0.19149,0.21366,0.23618,0.25856,0.28091
+0.30271,0.32477,0.34710,0.36895,0.39050,0.41176,0.43291,0.45428,0.47499,0.49590
+0.51652,0.53655,0.55660,0.57690,0.59677,0.61594,0.63488,0.65354,0.67179,0.68943
+0.70692,0.72373,0.74039,0.75669,0.77268,0.78826,0.80302,0.81741,0.83072,0.84383
+0.85628,0.86779,0.87875,0.88888,0.89819,0.90628,0.91354,0.92010,0.92581,0.93122
+0.93590,0.94022,0.94424,0.94789,0.95109,0.95411,0.95696,0.95968,0.96212,0.96440
+0.96655,0.96856,0.97037,0.97203,0.97355,0.97505,0.97655,0.97783,0.97907,0.98020
+0.98119,0.98225,0.98317,0.98402,0.98486,0.98567,0.98637,0.98699,0.98766,0.98828
+0.98879,0.98934,0.98976,0.99022,0.99057,0.99097,0.99128,0.99160,0.99197,0.99230
+0.99255,0.99286,0.99315,0.99353,0.99384,0.99411,0.99436,0.99462,0.99488,0.99516
+0.99543,0.99565,0.99592,0.99615,0.99635,0.99659,0.99683,0.99706,0.99729,0.99752
+0.99775,0.99801,0.99821,0.99842,0.99869,0.99894,0.99916,0.99936,0.99954,0.99965
+0.99981,0.99990,0.99994,0.99998,0.99999,0.99999,1.00000,1.00000,1.00000
+pwc2x2fract=0.00199,0.00459,0.00825,0.01336,0.02047,0.03001,0.04141,0.05551,0.07156
+0.08938,0.10884,0.12828,0.14914,0.17047,0.19202,0.21365,0.23563,0.25709,0.27892
+0.30046,0.32153,0.34326,0.36458,0.38581,0.40698,0.42799,0.44932,0.47030,0.49106
+0.51150,0.53138,0.55158,0.57098,0.59056,0.60954,0.62830,0.64689,0.66533,0.68352
+0.70087,0.71829,0.73517,0.75181,0.76792,0.78347,0.79867,0.81347,0.82784,0.84137
+0.85453,0.86673,0.87786,0.88858,0.89792,0.90628,0.91387,0.92019,0.92580,0.93077
+0.93514,0.93913,0.94276,0.94618,0.94945,0.95230,0.95493,0.95743,0.95989,0.96202
+0.96404,0.96618,0.96816,0.97004,0.97190,0.97359,0.97504,0.97648,0.97791,0.97926
+0.98051,0.98165,0.98271,0.98359,0.98456,0.98535,0.98609,0.98686,0.98751,0.98811
+0.98861,0.98912,0.98963,0.99010,0.99056,0.99099,0.99133,0.99171,0.99204,0.99237
+0.99263,0.99295,0.99328,0.99358,0.99384,0.99402,0.99432,0.99454,0.99480,0.99508
+0.99535,0.99556,0.99585,0.99609,0.99629,0.99652,0.99674,0.99695,0.99726,0.99746
+0.99768,0.99789,0.99813,0.99837,0.99859,0.99879,0.99897,0.99915,0.99932,0.99952
+0.99965,0.99980,0.99987,0.99993,0.99996,0.99998,0.99999,1.00000,1.00000
+pwc2x1fract=0.00170,0.00438,0.00799,0.01297,0.01998,0.02965,0.04099,0.05502,0.07039
+0.08804,0.10681,0.12637,0.14712,0.16766,0.18906,0.21084,0.23247,0.25366,0.27466
+0.29600,0.31723,0.33868,0.35980,0.38097,0.40165,0.42255,0.44389,0.46472,0.48535
+0.50587,0.52639,0.54694,0.56706,0.58653,0.60627,0.62558,0.64466,0.66316,0.68165
+0.69964,0.71699,0.73426,0.75143,0.76767,0.78368,0.79903,0.81422,0.82876,0.84266
+0.85582,0.86820,0.87986,0.89032,0.89970,0.90810,0.91524,0.92147,0.92739,0.93243
+0.93697,0.94111,0.94490,0.94827,0.95145,0.95434,0.95701,0.95935,0.96177,0.96404
+0.96607,0.96797,0.96978,0.97170,0.97339,0.97502,0.97646,0.97780,0.97910,0.98028
+0.98143,0.98252,0.98347,0.98433,0.98514,0.98586,0.98658,0.98729,0.98788,0.98846
+0.98894,0.98947,0.98995,0.99042,0.99092,0.99125,0.99162,0.99197,0.99231,0.99263
+0.99292,0.99318,0.99348,0.99379,0.99405,0.99435,0.99460,0.99489,0.99513,0.99534
+0.99558,0.99577,0.99599,0.99622,0.99643,0.99668,0.99688,0.99712,0.99735,0.99757
+0.99778,0.99803,0.99826,0.99853,0.99870,0.99890,0.99911,0.99930,0.99946,0.99960
+0.99972,0.99983,0.99990,0.99996,0.99998,1.00000,1.00000,1.00000,1.00000
+pwc2u2fract=0.00190,0.00455,0.00790,0.01242,0.01937,0.02810,0.03897,0.05192,0.06650
+0.08302,0.10075,0.11980,0.13975,0.16109,0.18206,0.20342,0.22533,0.24730,0.26868
+0.29072,0.31234,0.33443,0.35580,0.37740,0.39883,0.42047,0.44175,0.46293,0.48377
+0.50468,0.52532,0.54612,0.56583,0.58612,0.60570,0.62539,0.64457,0.66335,0.68173
+0.70029,0.71815,0.73531,0.75268,0.76919,0.78504,0.80067,0.81593,0.83047,0.84446
+0.85760,0.86984,0.88124,0.89130,0.90059,0.90904,0.91611,0.92246,0.92829,0.93336
+0.93771,0.94172,0.94555,0.94919,0.95243,0.95541,0.95810,0.96059,0.96286,0.96515
+0.96727,0.96925,0.97114,0.97283,0.97443,0.97603,0.97743,0.97867,0.97984,0.98102
+0.98199,0.98291,0.98390,0.98476,0.98564,0.98635,0.98699,0.98757,0.98822,0.98876
+0.98921,0.98963,0.99006,0.99044,0.99077,0.99109,0.99137,0.99168,0.99199,0.99225
+0.99254,0.99285,0.99307,0.99338,0.99365,0.99391,0.99419,0.99444,0.99469,0.99490
+0.99516,0.99541,0.99559,0.99582,0.99609,0.99629,0.99656,0.99680,0.99707,0.99728
+0.99755,0.99781,0.99803,0.99828,0.99846,0.99866,0.99886,0.99905,0.99930,0.99945
+0.99964,0.99976,0.99987,0.99993,0.99998,0.99999,0.99999,1.00000,1.00000
+pwc2u1fract=0.00193,0.00468,0.00817,0.01287,0.01956,0.02850,0.03902,0.05167,0.06618
+0.08269,0.10040,0.11982,0.13946,0.16005,0.18114,0.20293,0.22453,0.24651,0.26857
+0.29092,0.31319,0.33479,0.35672,0.37829,0.39949,0.42087,0.44214,0.46275,0.48382
+0.50467,0.52496,0.54552,0.56559,0.58529,0.60511,0.62477,0.64422,0.66308,0.68211
+0.70019,0.71795,0.73570,0.75295,0.77001,0.78631,0.80236,0.81788,0.83214,0.84598
+0.85942,0.87170,0.88297,0.89315,0.90241,0.91065,0.91757,0.92400,0.92988,0.93482
+0.93934,0.94352,0.94726,0.95063,0.95381,0.95685,0.95950,0.96206,0.96440,0.96662
+0.96858,0.97051,0.97213,0.97373,0.97521,0.97674,0.97804,0.97930,0.98039,0.98141
+0.98229,0.98321,0.98408,0.98479,0.98547,0.98613,0.98672,0.98724,0.98777,0.98829
+0.98883,0.98922,0.98963,0.99012,0.99054,0.99092,0.99123,0.99153,0.99187,0.99223
+0.99258,0.99286,0.99317,0.99347,0.99377,0.99405,0.99431,0.99454,0.99472,0.99502
+0.99528,0.99549,0.99576,0.99603,0.99630,0.99649,0.99669,0.99698,0.99724,0.99743
+0.99762,0.99782,0.99806,0.99831,0.99858,0.99879,0.99905,0.99927,0.99947,0.99963
+0.99974,0.99986,0.99991,0.99994,0.99997,0.99999,1.00000,1.00000,1.00000
-- 
GitLab