From 29ee9a42bb1213697f2cc7223a0d058cce291e11 Mon Sep 17 00:00:00 2001
From: Yero1990 <cyero002@fiu.edu>
Date: Sun, 15 Oct 2017 09:32:13 -0400
Subject: [PATCH] add re-structured SHMS DC Calibration in terms of classes and
methods. This is a version of the SHMS DC Calibration code that seems to
work. Along with the parameter files produce by it. (#228)
---
CALIBRATION/shms_dc_calib/scripts/DC_calib.C | 872 ++++++++++++++++++
CALIBRATION/shms_dc_calib/scripts/DC_calib.h | 177 ++++
.../shms_dc_calib/scripts/main_calib.C | 42 +
PARAM/SHMS/DC/CALIB/pdc_calib.param | 338 +++----
PARAM/SHMS/DC/CALIB/pdc_tzero_per_wire.param | 152 +--
5 files changed, 1336 insertions(+), 245 deletions(-)
create mode 100644 CALIBRATION/shms_dc_calib/scripts/DC_calib.C
create mode 100644 CALIBRATION/shms_dc_calib/scripts/DC_calib.h
create mode 100644 CALIBRATION/shms_dc_calib/scripts/main_calib.C
diff --git a/CALIBRATION/shms_dc_calib/scripts/DC_calib.C b/CALIBRATION/shms_dc_calib/scripts/DC_calib.C
new file mode 100644
index 00000000..08f6e9ab
--- /dev/null
+++ b/CALIBRATION/shms_dc_calib/scripts/DC_calib.C
@@ -0,0 +1,872 @@
+//SHMS DC Calibration: Implementation
+#include<iostream>
+#include "DC_calib.h"
+
+using namespace std;
+
+//_____________________________________________________________
+DC_calib::DC_calib(TString a, TString b, const int c, Long64_t d)
+
+ :spec(a), //set spectrometer to 'HMS', or 'SHMS' ex. DC_Calib(HMS, pdc_replay.C, 488, 50000)
+ ifile_name(b), //initialization list
+ run_NUM(c),
+ num_evts(d)
+
+{
+ //Initialize pointers
+ tree = NULL;
+ in_file = NULL;
+ out_file = NULL;
+ main_dir = NULL;
+ plane_dt = NULL;
+ plane_dt_corr = NULL;
+ cell_dt = NULL;
+ cell_dt_corr = NULL;
+ fitted_cell_dt = NULL;
+ dt_vs_wire = NULL;
+ dt_vs_wire_corr = NULL;
+
+ entries = NULL;
+ t_zero = NULL;
+ t_zero_err = NULL;
+ bin_max = NULL;
+ bin_maxContent = NULL;
+ time_max = NULL;
+ twenty_perc_maxContent = NULL;
+ ref_time = NULL;
+ tZero_fit = NULL;
+
+ graph = NULL;
+ gr1_canv = NULL;
+
+
+
+
+}
+
+
+//____________________________________________________________
+DC_calib::~DC_calib()
+{
+ cout << "calling the destructor " << endl;
+
+ delete in_file; in_file = NULL;
+ delete out_file; out_file = NULL;
+ delete graph; graph = NULL;
+ delete gr1_canv; gr1_canv = NULL;
+ //Delete Arrays/pointers to free up 'heap' space
+
+ delete [] dt_vs_wire; dt_vs_wire = NULL;
+ delete [] dt_vs_wire_corr; dt_vs_wire_corr = NULL;
+ delete [] plane_dt; plane_dt = NULL;
+ delete [] plane_dt_corr; plane_dt_corr = NULL;
+
+ //free 2d dynamic array cell_dt[][]
+ for(int ip = 0; ip<NPLANES; ip++)
+ {
+ delete [] entries[ip];
+ delete [] t_zero[ip];
+ delete [] t_zero_err[ip];
+ delete [] cell_dt[ip];
+ delete [] cell_dt_corr[ip];
+ delete [] fitted_cell_dt[ip];
+ delete [] bin_max[ip];
+ delete [] bin_maxContent[ip];
+ delete [] time_max[ip];
+ delete [] twenty_perc_maxContent[ip];
+ delete [] ref_time[ip];
+ }
+
+ delete [] entries; entries = NULL;
+ delete [] t_zero; t_zero = NULL;
+ delete [] t_zero_err; t_zero_err = NULL;
+ delete [] cell_dt; cell_dt = NULL;
+ delete [] cell_dt_corr; cell_dt_corr = NULL;
+ delete [] fitted_cell_dt; fitted_cell_dt = NULL;
+ delete [] bin_max; bin_max = NULL;
+ delete [] bin_maxContent; bin_maxContent = NULL;
+ delete [] time_max; time_max = NULL;
+ delete [] twenty_perc_maxContent; twenty_perc_maxContent = NULL;
+ delete [] ref_time; ref_time = NULL;
+}
+
+//____________________________________________________________
+void DC_calib::printInitVar()
+{
+ cout << "Initialization variables: \n"
+ "Input File: " << ifile_name << "\n"
+ "Run #: " << run_NUM << "\n"
+ "Events: " << num_evts << endl;
+}
+
+
+
+//___________________________________________________________
+void DC_calib::SetPlaneNames()
+{
+
+ DETECTOR = "dc";
+
+ //initialize DC plane names
+ if(spec=="SHMS")
+ {
+ SPECTROMETER = "P";
+
+ plane_names[0]="1u1", nwires[0] = 107;
+ plane_names[1]="1u2", nwires[1] = 107;
+ plane_names[2]="1x1", nwires[2] = 79;
+ plane_names[3]="1x2", nwires[3] = 79;
+ plane_names[4]="1v1", nwires[4] = 107;
+ plane_names[5]="1v2", nwires[5] = 107;
+ plane_names[6]="2v2", nwires[6] = 107;
+ plane_names[7]="2v1", nwires[7] = 107;
+ plane_names[8]="2x2", nwires[8] = 79;
+ plane_names[9]="2x1", nwires[9] = 79;
+ plane_names[10]="2u2", nwires[10] = 107;
+ plane_names[11]="2u1", nwires[11] = 107;
+
+ }
+
+
+ else if(spec=="HMS")
+ {
+ cout << "hms calibration is a work-in-praw gress" << endl;
+ }
+
+}
+
+
+//____________________________________________________________________________________
+
+void DC_calib::GetDCLeafs()
+{
+ cout << "DC_calib::GetDCLeafs " << endl;
+
+ //open input root file
+ in_file = new TFile(ifile_name, "READ" );
+
+ //Get the tree
+ tree = (TTree*)in_file->Get("T");
+
+
+ //Loop over each plane
+ for (int ip = 0; ip < NPLANES; ip++)
+ {
+
+ base_name = SPECTROMETER+"."+DETECTOR+"."+plane_names[ip];
+ ndatatime = "Ndata."+base_name+".time";
+ ndatawirenum = "Ndata."+base_name+".wirenum";
+ drifttime = base_name+".time";
+ wirenum = base_name+".wirenum";
+
+ //Set Branch Address
+ tree->SetBranchAddress(wirenum, wire_num[ip]);
+ tree->SetBranchAddress(drifttime, drift_time[ip]);
+ tree->SetBranchAddress(ndatawirenum, &ndata_wirenum[ip]);
+ tree->SetBranchAddress(ndatatime, &ndata_time[ip]);
+
+ }
+
+ ntrack = SPECTROMETER + "." + DETECTOR + ".ntrack";
+ etracknorm = SPECTROMETER + ".cal.etracknorm";
+
+ tree->SetBranchAddress(ntrack, &dc_ntrack);
+ tree->SetBranchAddress(etracknorm, &psh_etracknorm);
+
+
+}
+
+void DC_calib::AllocateDynamicArrays()
+{
+
+
+ //Allocate 1D dynamic arrays
+ plane_dt = new TH1F[NPLANES]; //create plane drift time histo 1Darray ( get_pdc_time_histo.C )
+ plane_dt_corr = new TH1F[NPLANES]; //create plane drift times to store after applying tzero correction
+ dt_vs_wire = new TH2F[NPLANES]; //create wire drift time and 'drifttime vs wire' arrays
+ dt_vs_wire_corr = new TH2F[NPLANES];
+
+ //Allocate 2D dynamic arrays
+ entries = new Int_t*[NPLANES];
+ t_zero = new Double_t*[NPLANES];
+ t_zero_err = new Double_t*[NPLANES];
+ cell_dt = new TH1F*[NPLANES]; /*create array to store cell drift times*/
+ cell_dt_corr = new TH1F*[NPLANES];
+ fitted_cell_dt = new TH1F*[NPLANES]; /*create array to store cell drift times*/
+ bin_max = new Int_t*[NPLANES]; /*Array to store the bin number corresponding to the drift time distribution peak*/
+ bin_maxContent = new Int_t*[NPLANES]; /*Array to store the content (# events) corresponding to the bin with maximum content*/
+ time_max = new Double_t*[NPLANES]; /*Array to store the x-axis(drift time (ns)) corresponding to bin_max*/
+ twenty_perc_maxContent = new Double_t*[NPLANES]; /*Array to store 20% of maximum bin content (peak)*/
+ ref_time = new Double_t*[NPLANES]; /*Array to store ref_time(time corresp. to 20% of peak) times for each sense wire*/
+
+
+ for(int ip=0; ip<NPLANES; ip++)
+ {
+ entries[ip] = new Int_t[nwires[ip]];
+ t_zero[ip] = new Double_t[nwires[ip]];
+ t_zero_err[ip] = new Double_t[nwires[ip]];
+ cell_dt[ip] = new TH1F[nwires[ip]];
+ cell_dt_corr[ip] = new TH1F[nwires[ip]];
+ fitted_cell_dt[ip] = new TH1F[nwires[ip]];
+ bin_max[ip] = new Int_t[nwires[ip]];
+ bin_maxContent[ip] = new Int_t[nwires[ip]];
+ time_max[ip] = new Double_t[nwires[ip]];
+ twenty_perc_maxContent[ip] = new Double_t[nwires[ip]];
+ ref_time[ip] = new Double_t[nwires[ip]];
+ }
+
+
+}
+
+//_______________________________________________________________
+void DC_calib::CreateHistoNames()
+{
+
+ for(int ip=0; ip<NPLANES; ip++)
+ {
+
+ //Set-Up plane drift time histo labels
+ plane_dt_name = plane_names[ip]+"_time";
+ plane_dt_title = "SHMS DC, Plane "+plane_names[ip]+" Drift Time";
+
+ plane_dt[ip].SetName(plane_dt_name);
+ plane_dt[ip].SetTitle(plane_dt_title);
+ plane_dt[ip].SetBins(NBINS, MINBIN, MAXBIN);
+ plane_dt[ip].SetXTitle("Drift Time (ns)");
+ plane_dt[ip].SetYTitle("Number of Entries / 1 ns");
+
+
+ plane_dt_name_corr = plane_names[ip]+"corrected_time";
+ plane_dt_title_corr = "SHMS DC, Plane "+plane_names[ip]+" Corrected Drift Time";
+
+ plane_dt_corr[ip].SetName(plane_dt_name_corr);
+ plane_dt_corr[ip].SetTitle(plane_dt_title_corr);
+ plane_dt_corr[ip].SetBins(NBINS, MINBIN, MAXBIN);
+ plane_dt_corr[ip].SetXTitle("Drift Time (ns)");
+ plane_dt_corr[ip].SetYTitle("Number of Entries / 1 ns");
+
+
+
+ //Set-Up Drift Time vs. Wire Number Histos labels
+ dt_vs_wire_name = "dt_vs_wire_plane_"+plane_names[ip];
+ dt_vs_wire_title = "SHMS Drift Time vs. Wire: Plane "+plane_names[ip];
+
+ dt_vs_wire[ip].SetName(dt_vs_wire_name);
+ dt_vs_wire[ip].SetTitle(dt_vs_wire_title);
+ dt_vs_wire[ip].SetBins(nwires[ip], 0., nwires[ip], NBINS, MINBIN, MAXBIN);
+ dt_vs_wire[ip].SetXTitle("Wire Number");
+ dt_vs_wire[ip].SetYTitle("Drift Time (ns) / 1 ns");
+
+ dt_vs_wire_corr[ip].SetName(dt_vs_wire_name);
+ dt_vs_wire_corr[ip].SetTitle(dt_vs_wire_title);
+ dt_vs_wire_corr[ip].SetBins(nwires[ip], 0., nwires[ip], NBINS, MINBIN, MAXBIN);
+ dt_vs_wire_corr[ip].SetXTitle("Wire Number");
+ dt_vs_wire_corr[ip].SetYTitle("Drift Time (ns) / 1 ns");
+
+ //Set-Up cell drift time histos names
+ for (wire = 0; wire < nwires[ip]; wire++)
+ {
+ cell_dt_name = Form("Wire_%d", wire+1);
+ cell_dt_title = "SHMS DC Plane " +plane_names[ip] + Form(": Wire_%d", wire+1);
+
+ cell_dt[ip][wire].SetName(cell_dt_name);
+ cell_dt[ip][wire].SetTitle(cell_dt_title);
+ cell_dt[ip][wire].SetBins(NBINS, MINBIN, MAXBIN);
+ cell_dt[ip][wire].SetXTitle("Drift Time (ns)");
+ cell_dt[ip][wire].SetYTitle("Number of Entries / 1 ns");
+
+ fitted_cell_dt_name = Form("Wire_%d", wire+1);
+ fitted_cell_dt_title = "SHMS DC Plane " +plane_names[ip] + Form(": Wire_%d", wire+1);
+
+ fitted_cell_dt[ip][wire].SetName(fitted_cell_dt_name);
+ fitted_cell_dt[ip][wire].SetTitle(fitted_cell_dt_title);
+ fitted_cell_dt[ip][wire].SetBins(NBINS, MINBIN, MAXBIN);
+ fitted_cell_dt[ip][wire].SetXTitle("Drift Time (ns)");
+ fitted_cell_dt[ip][wire].SetYTitle("Number of Entries / 1 ns");
+
+ cell_dt_corr[ip][wire].SetName(cell_dt_name);
+ cell_dt_corr[ip][wire].SetTitle(cell_dt_title);
+ cell_dt_corr[ip][wire].SetBins(NBINS, MINBIN, MAXBIN);
+ cell_dt_corr[ip][wire].SetXTitle("Drift Time (ns)");
+ cell_dt_corr[ip][wire].SetYTitle("Number of Entries / 1 ns");
+
+ }
+
+ }
+
+}
+
+
+//________________________________________________________________
+void DC_calib::EventLoop()
+{
+
+
+ //Loop over all entries
+ for(Long64_t i=0; i<num_evts; i++)
+ {
+ //cout << "entry: " << i << endl;
+ tree->GetEntry(i);
+
+ for(Int_t ip=0; ip<NPLANES; ip++)
+ {
+ // cout << "PLANE: " << ip << endl;
+
+//-----------------------------------------------------------------------------------------
+
+ //Loop over number of hits for each trigger in each DC plane
+ for(Int_t j = 0, k = 0; j < ndata_time[ip], k < ndata_wirenum[ip]; j++, k++)
+ {
+
+
+ //get wire hit for ith event in 'ip' plane
+ wire = int(wire_num[ip][k]);
+
+ // if (dc_ntrack == 1 && psh_etracknorm > 0.6) {
+ //Fill uncorrected plane drift times (from: get_pdc_time_histo.C )
+ plane_dt[ip].Fill(drift_time[ip][j]);
+
+
+
+ //cout << "wire hit: " << wire << endl;
+ //Fill uncorrected wire (cell) drift times and dt-vs-wire Histograms
+ dt_vs_wire[ip].Fill(wire_num[ip][k], drift_time[ip][j]);
+ cell_dt[ip][wire-1].Fill(drift_time[ip][j]);
+ fitted_cell_dt[ip][wire-1].Fill(drift_time[ip][j]);
+ // cout << "drift time: " << drift_time[ip][j] << endl;
+
+ //} // end cut
+
+
+//-----------------------------------------------------------------------------------------
+
+ }
+
+
+
+
+ } //end loop over planes
+
+ } //end loop over events
+
+}
+
+
+//_________________________________________________________________________
+/*
+void DC_calib::CalcT0Historical()
+{
+ //Loop over DC PLANES
+ for (int ip = 0; ip<NPLANES; ip++)
+ {
+ //Loop over DC wires
+ for (wire = 0; wire < nwires[ip]; wire++)
+ {
+ double tempt0 = 0;
+
+ if (cell_dt[ip][wire].GetEntries()>500){
+ bin_max[ip][wire] = cell_dt[ip][wire].GetMaximumBin()+10;
+ double histMax = 0;
+ for (Int_t bin = 1; bin <= bin_max[ip][wire]+10; bin++){
+ double avgMax = 0;
+ double sumMax = 0;
+ int nseq = 0;
+ int nmin = bin+8;
+
+ for (int ii = bin; ii<=nmin; ii++){
+ double cc = cell_dt[ip][wire].GetBinContent(ii);
+ if (cc>10){
+ sumMax += cc;
+ nseq++;
+ }
+ }
+
+ if (nseq>6){
+ avgMax = sumMax/(double)nseq;
+ }
+ if (avgMax>histMax){
+ histMax = avgMax;
+ }
+ }//end loop bins
+ twenty_perc_maxContent[ip][wire] = 0.2*histMax;
+
+ //Use the 20% to find the t0 offset
+ for (Int_t bin = 10; bin <= MAXBIN-8; bin++){
+ double binMax = 0;
+ int nseq = 0;
+ int nmin = bin+8;
+
+ for (int ii = bin; ii<=nmin; ii++){
+ double cc = cell_dt[ip][wire].GetBinContent(ii);
+ if (cc>twenty_perc_maxContent[ip][wire]){
+ nseq++;
+ }
+ }
+
+ if (nseq>7){
+ int binRef = 0;
+ cell_dt[ip][wire].GetBinWithContent(twenty_perc_maxContent[ip][wire], binRef, 1,bin_max[ip][wire],10.0);
+ tempt0 = cell_dt[ip][wire].GetXaxis()->GetBinCenter(binRef);
+ break;
+ }
+ }//end loop bins
+ }//end entries>500
+
+ t_zero[ip][wire] = tempt0;
+
+ }//end loop wires
+ }//end loop planes
+}
+*/
+
+//_________________________________________________________________________
+void DC_calib::GetTwentyPercent_Peak()
+{
+
+ //Loop over DC PLANES
+ for (int ip = 0; ip<NPLANES; ip++)
+ {
+
+ //Loop over DC wires
+ for (wire = 0; wire < nwires[ip]; wire++)
+ {
+
+ bin_max[ip][wire] = cell_dt[ip][wire].GetMaximumBin(); //Get bin with Maximum Content
+ bin_maxContent[ip][wire] = cell_dt[ip][wire].GetBinContent(bin_max[ip][wire]); //Get content of bin_max
+ time_max[ip][wire] = cell_dt[ip][wire].GetXaxis()->GetBinCenter(bin_max[ip][wire]); //Get time (ns) [x-axis] corresponding to bin_max
+ twenty_perc_maxContent[ip][wire] = bin_maxContent[ip][wire] * 0.20;
+ //Calculate 20% of max bin content
+ //ref_time[ip][wire] = cell_dt[ip][wire].GetBinCenter(cell_dt[ip][wire].FindBin(twenty_perc_maxContent[ip][wire]));
+ //Loop over DC drift time bins
+ for (Int_t bin = 1; bin <= bin_max[ip][wire]; bin++)
+ {
+ content_bin = cell_dt[ip][wire].GetBinContent(bin); //get bin content for all bins in a wire
+
+ content.push_back(content_bin); //add bin content to array
+ bin_num.push_back(bin); //add bin number to array
+
+ // check if 2 bin contents have been stored and examine if these contents exceed or not 20% of peak
+ if (content.size() == 5) {
+ //initialize counter to count how many bin contents >= 20%
+ counts = 0;
+
+ // Loop over 2 bin contents stored in array content
+ for (Int_t j=0; j<5; j++)
+ {
+ if(content[j] >= twenty_perc_maxContent[ip][wire])
+ {
+ counts = counts+1;
+ if(counts >= 4) { goto stop;}
+ }
+
+ content.clear();
+ bin_num.clear();
+
+ }
+
+ }
+
+ }
+ stop:
+ ref_time[ip][wire] = cell_dt[ip][wire].GetXaxis()->GetBinCenter(bin_num[0]); //Get time corresponding ~20% Max BIN CONTENT
+
+
+ }
+
+ }
+
+}
+
+
+//____________________________________________________________________________________
+void DC_calib::FitWireDriftTime()
+{
+
+ //Loop over planes
+ for (Int_t ip = 0; ip < NPLANES; ip++)
+ {
+
+ //Loop over DC sense wires
+ for (wire = 0; wire < nwires[ip]; wire++)
+ {
+ //Get Bin no. associated with the reference time
+ binx = cell_dt[ip][wire].GetXaxis()->FindBin(ref_time[ip][wire]);
+
+ //Get time corresponding to bin (fit range)
+ time_init = cell_dt[ip][wire].GetXaxis()->GetBinCenter(binx - 10); //choose bin range over which to fit
+ time_final = cell_dt[ip][wire].GetXaxis()->GetBinCenter(binx + 10);
+
+ //Create Fit Function
+ tZero_fit = new TF1("tZero_fit", "[0]*x + [1]", time_init, time_final);
+
+ //Set Parameter Names and Values
+ tZero_fit->SetParName(0, "slope");
+ tZero_fit->SetParName(1, "y-int");
+ tZero_fit->SetParameter(0, 1.0);
+ tZero_fit->SetParameter(1, 1.0);
+
+ entries[ip][wire] = fitted_cell_dt[ip][wire].GetEntries();
+
+ //Avoid mssg due to zero entries: Warning in <Fit>: Fit data is empty
+ // if (entries[ip][wire] !=0)
+ //{
+ //Fit Function in specified range
+ fitted_cell_dt[ip][wire].Fit("tZero_fit", "QR");
+ //}
+ gStyle->SetOptFit(1);
+ //Get Parameters and their errors
+ m = tZero_fit->GetParameter(0);
+ y_int = tZero_fit->GetParameter(1);
+ m_err = tZero_fit->GetParError(0);
+ y_int_err = tZero_fit->GetParError(1);
+ std_dev = fitted_cell_dt[ip][wire].GetStdDev();
+
+ //Require sufficient events and NOT CRAZY! tzero values, otherwis, set t0 to ZERO
+ if (abs(-y_int/m) < std_dev && m > 0.0 && entries[ip][wire]>100)
+ {
+ t_zero[ip][wire] = - y_int/m ;
+ t_zero_err[ip][wire] = sqrt(y_int_err*y_int_err/(m*m) + y_int*y_int*m_err*m_err/(m*m*m*m) );
+
+ }
+ else // (entries[ip][wire] < 1000 || m < 0.0)
+ {
+
+ t_zero[ip][wire] = 0.0;
+ t_zero_err[ip][wire] = 0.0;
+
+ }
+
+
+ } //END LOOP OVER WIRES
+
+
+ }//END LOOP OVER PLANES
+
+
+}
+
+
+//________________________________________________________________
+void DC_calib::Calculate_tZero()
+{
+
+ //CalcT0Historical();
+
+ GetTwentyPercent_Peak();
+ FitWireDriftTime();
+}
+
+
+//________________________________________________________________
+void DC_calib::WriteToFile(Int_t debug = 0)
+{
+
+
+ //create output ROOT file to write UnCALIB./CALIB. histos
+ ofile_name = "SHMS_DC_driftimes.root";
+ out_file = new TFile(ofile_name, "RECREATE");
+
+
+ //***NOTE*** debug = 0 (OFF), 1 (ON)
+ //-------------------------------------------------------------------
+ if (debug == 1)
+ {
+ //------write uncorrected plane drift time histos to a directory on FILE--------
+ main_dir = out_file->mkdir("uncorr_plane_times");
+ main_dir->cd();
+
+ for(Int_t ip=0; ip<NPLANES; ip++)
+ {
+ //write histos to root file
+ plane_dt[ip].Write();
+ }
+
+ //------write corrected plane drift time histos to a directory on FILE--------
+ main_dir = out_file->mkdir("corr_plane_times");
+ main_dir->cd();
+
+ for(Int_t ip=0; ip<NPLANES; ip++)
+ {
+ //write histos to root file
+ plane_dt_corr[ip].Write();
+ }
+
+
+
+
+
+ //----------------------------------------------------------------------
+
+
+ //--------write uncorrected cell drift times histos to FILE--------
+
+ main_dir = out_file->mkdir("uncorr_wire_times");
+
+ for (int ip=0; ip<NPLANES; ip++)
+ {
+
+ // create planes sub-directories to store wire drift times
+ (main_dir->mkdir("plane "+plane_names[ip], ""))->cd();
+
+ dt_vs_wire[ip].Write(); //write 2d drifttimet_vs_wire to FILE
+
+ for (wire = 0; wire < nwires[ip]; wire++)
+ {
+
+ cell_dt[ip][wire].Write(); //Write wire driftimes to FILE
+
+ }
+
+ } //END LOOP OVER PLANES
+
+ //--------write corrected cell drift times histos to FILE--------
+
+ main_dir = out_file->mkdir("corr_wire_times");
+
+ for (int ip=0; ip<NPLANES; ip++)
+ {
+
+ // create planes sub-directories to store wire drift times
+ (main_dir->mkdir("plane "+plane_names[ip], ""))->cd();
+
+ dt_vs_wire_corr[ip].Write(); //write 2d drifttimet_vs_wire to FILE
+
+ for (wire = 0; wire < nwires[ip]; wire++)
+ {
+
+ cell_dt_corr[ip][wire].Write(); //Write wire driftimes to FILE
+
+ }
+
+ } //END LOOP OVER PLANES
+
+
+
+
+ //------------------------------------------------------------------------
+
+ //-------Write Fitted Wire Drift Time histos to FILE----------------------------
+
+
+ main_dir = out_file->mkdir("fitted_wire_drift_times");
+
+ for (int ip=0; ip<NPLANES; ip++)
+ {
+
+ // create planes sub-directories to store fitted wire drift times
+ (main_dir->mkdir("plane "+plane_names[ip], ""))->cd();
+
+ for (wire = 0; wire < nwires[ip]; wire++)
+ {
+ fitted_cell_dt[ip][wire].Write(); //Write wire driftimes to FILE
+
+ }
+
+ } //END LOOP OVER PLANES
+
+ //-----Write 'tzero' values to a TEXT FILE--------------------
+
+ //open a text FILE to write
+
+ for (int ip = 0; ip < NPLANES; ip++)
+ {
+ otxtfile_name = "t_zero_values_"+plane_names[ip]+".dat";
+ out_txtFILE.open(otxtfile_name);
+ out_txtFILE << "#Plane_" + plane_names[ip] << endl;
+ out_txtFILE << "#Wire " << setw(12) << "tzero " << setw(12) << "t_zero_err " << setw(12) << "entries" << endl;
+
+ for (wire = 0; wire < nwires[ip]; wire++)
+ {
+ out_txtFILE << wire+1 << " " << t_zero[ip][wire] << " " << t_zero_err[ip][wire] << " " << entries[ip][wire] << endl;
+ }
+
+ out_txtFILE.close();
+
+ }
+
+//------------------------------------------------------------------------
+
+ main_dir = out_file->mkdir("t0_vs_wire");
+ main_dir->cd();
+ for (int ip=0; ip<NPLANES; ip++)
+ {
+
+ gr1_canv = new TCanvas("gr1", "", 2000, 500);
+ gr1_canv->SetGrid();
+ //write TGraph: tzero v. wire number to root file
+ itxtfile_name = "t_zero_values_"+plane_names[ip]+".dat";
+ graph = new TGraphErrors(itxtfile_name, "%lg %lg %lg");
+ graph->SetName("graph");
+
+ graph_title = "DC "+plane_names[ip]+": t0 vs. Wire";
+ graph->SetTitle(graph_title);
+
+ graph->SetMarkerStyle(20);
+ graph->SetMarkerColor(1);
+
+ graph->GetXaxis()->SetLimits(0., nwires[ip]);
+ graph->GetXaxis()->SetTitle("Wire Number");
+ graph->GetXaxis()->CenterTitle();
+ graph->GetYaxis()->SetTitle("t-Zero (ns)");
+ graph->GetYaxis()->CenterTitle();
+ graph->GetYaxis()->SetRangeUser(-50.0, 50.0);
+
+ graph->Draw("AP");
+ gr1_canv->Update();
+ gr1_canv->Write(graph_title); //write to a root file
+
+ }
+
+ } //END DEBUG
+
+
+}
+
+
+//__________________________________________________________________________
+void DC_calib::WriteTZeroParam()
+{
+ otxtfile_name = "./pdc_tzero_per_wire_"+std::to_string(run_NUM)+".param";
+ out_txtFILE.open(otxtfile_name);
+
+ for (int ip=0; ip<NPLANES; ip++) {
+
+ //write plane headers
+ out_txtFILE << "ptzero"+plane_names[ip] << "=" << endl;
+
+ for (wire=0; wire<nwires[ip]; wire++)
+ {
+ if (wire <= 10)
+ {
+ out_txtFILE << setprecision(6) << t_zero[ip][wire] << fixed << ",";
+ }
+ else if (wire>10 && wire <(nwires[ip]-1))
+ {
+ out_txtFILE << setprecision(6) << t_zero[ip][wire] << ((wire+1) % 16 ? ", " : "\n") << fixed;
+ }
+ else if (wire==78 || wire == 106)
+ {
+ out_txtFILE << setprecision(6) << t_zero[ip][wire] << fixed << endl;
+ }
+
+ } //END LOOP OVER WIRES
+
+ } //END LOOP OVER PLANES
+
+ out_txtFILE.close();
+
+}
+
+
+
+//_______________________________________________________________________
+void DC_calib::ApplyTZeroCorrection()
+{
+
+
+ cout << "ApplyT0Corr "<< endl;
+
+ //Loop over all entries
+ for(Long64_t i=0; i<num_evts; i++)
+ {
+ tree->GetEntry(i);
+ for(Int_t ip=0; ip<NPLANES; ip++)
+ {
+
+ //Loop over number of hits for each trigger in each DC plane
+
+ for(Int_t j = 0, k = 0; j < ndata_time[ip], k < ndata_wirenum[ip]; j++, k++)
+ {
+ //get wire hit for ith event in 'ip' plane
+ wire = int(wire_num[ip][k]);
+
+ //if (dc_ntrack == 1 && psh_etracknorm > 0.6) {
+
+ //Fill corrected plane drift times
+ plane_dt_corr[ip].Fill(drift_time[ip][j] - t_zero[ip][wire-1]);
+ cell_dt_corr[ip][wire-1].Fill(drift_time[ip][j] - t_zero[ip][wire-1]);
+ dt_vs_wire_corr[ip].Fill(wire_num[ip][k], drift_time[ip][j] - t_zero[ip][wire-1]);
+ //} //end cut
+
+ } //end loop over hits
+
+
+ } //end loop over planes
+
+ } //end loop over events
+}
+
+void DC_calib::WriteLookUpTable()
+{
+ otxtfile_name = "./pdc_calib_"+std::to_string(run_NUM)+".param";
+ out_txtFILE.open(otxtfile_name);
+
+ //Set headers for subsequent columns of data
+ out_txtFILE << Form("; Lookup Table: RUN %d", run_NUM) << "\n";
+ out_txtFILE << "; number of bins in time to distance lookup table" << "\n";
+ out_txtFILE << Form("pdriftbins = %d", TOTAL_BINS+1) << "\n";
+ out_txtFILE << "; number of 1st bin in table in ns" << "\n";
+ out_txtFILE << "pdrift1stbin=0" << "\n";
+ out_txtFILE << "; bin size in ns" << "\n";
+ out_txtFILE << "pdriftbinsz=1" << "\n";
+
+//Loop over each plane of shms Drift Chambers (DC1 & DC2)
+
+ for (int ip=0; ip<NPLANES; ip++){
+
+
+ //Get bin corresponding to t0 = 0 ns
+ bin_t0[ip] = plane_dt_corr[ip].GetXaxis()->FindBin(0.0);
+
+ //Get final bin
+ bin_final[ip] = bin_t0[ip] + TOTAL_BINS;
+
+ //Find total BIN Content over entire integration range
+ binContent_TOTAL[ip] = 0.; //set counter to zero
+
+ for (int bin = bin_t0[ip]; bin <= bin_final[ip]; bin ++ ) {
+
+ bin_Content[ip] = plane_dt_corr[ip].GetBinContent(bin);
+
+ binContent_TOTAL[ip] = bin_Content[ip] + binContent_TOTAL[ip];
+
+ // cout << "Bin: " << bin << endl;
+ // cout << "Content " << bin_Content[ip] << endl;
+ // cout << "Content SUM : " << binContent_TOTAL[ip] << endl;
+ }
+
+ TString headers = "pwc" + plane_names[ip] + "fract=";
+ out_txtFILE << headers;
+
+ //Calculate LookUp Value
+
+ binSUM[ip] = 0.0;
+ bin_count = 0;
+
+ for (int bin = bin_t0[ip]; bin <= bin_final[ip]; bin++) {
+
+ bin_Content[ip] = plane_dt_corr[ip].GetBinContent(bin);
+ binSUM[ip] = binSUM[ip] + bin_Content[ip];
+
+
+ lookup_value[ip] = binSUM[ip] / binContent_TOTAL[ip];
+ bin_count = bin_count + 1;
+
+ if (bin_count <= 16 ) {
+ out_txtFILE << setprecision(5) << lookup_value[ip] << fixed << ",";
+ }
+
+ else if (bin_count >16 && bin_count <= TOTAL_BINS) {
+ out_txtFILE << setprecision(5) << lookup_value[ip] << ((bin_count+1) % 20 ? "," : "\n") << fixed;
+ }
+ else {
+ out_txtFILE << setprecision(5) << lookup_value[ip] << fixed << endl;
+ }
+
+ } //END LOOP OVER plane drift time BINS
+
+ } //END LOOP OVER PLANES
+
+ out_txtFILE.close();
+
+}
diff --git a/CALIBRATION/shms_dc_calib/scripts/DC_calib.h b/CALIBRATION/shms_dc_calib/scripts/DC_calib.h
new file mode 100644
index 00000000..1f380d7d
--- /dev/null
+++ b/CALIBRATION/shms_dc_calib/scripts/DC_calib.h
@@ -0,0 +1,177 @@
+//SHMS DC Calibration Class
+#ifndef DC_CALIB_H
+#define DC_CALIB_H
+
+#define NPLANES 12
+#define NBINS 400
+#define MINBIN -50.5
+#define MAXBIN 349.5
+#define TOTAL_BINS 274
+class DC_calib
+{
+ public:
+
+ //consructor and destructor
+ DC_calib(TString a, TString b, const Int_t c, Long64_t d);
+ ~DC_calib();
+
+
+ //Define Functions
+ void printInitVar();
+ void SetPlaneNames();
+ void GetDCLeafs();
+ void AllocateDynamicArrays();
+ void CreateHistoNames();
+ void EventLoop();
+ void WriteToFile(Int_t debug);
+ void CalcT0Historical();
+ void Calculate_tZero();
+ void GetTwentyPercent_Peak();
+ void FitWireDriftTime();
+ void WriteTZeroParam();
+ void ApplyTZeroCorrection();
+ void WriteLookUpTable();
+
+
+ private:
+
+ Int_t run_NUM;
+ Long64_t num_evts;
+
+
+ TTree *tree;
+ Long64_t nentries;
+
+ TString SPECTROMETER;
+ TString spec;
+ TString DETECTOR;
+ TString plane_names[NPLANES];
+
+ TString base_name;
+ TString ndatatime;
+ TString ndatawirenum;
+
+ TString drifttime;
+ TString wirenum;
+
+ TString ntrack;
+ TString etracknorm;
+
+ Double_t dc_ntrack;
+ Double_t psh_etracknorm;
+
+ Int_t wire;
+
+ Int_t ndata_time[NPLANES];
+ Double_t drift_time[NPLANES][1000];
+
+ Int_t ndata_wirenum[NPLANES];
+ Double_t wire_num[NPLANES][1000];
+
+ Int_t nwires[NPLANES];
+
+ //Declare variables to plot and save histo (dt = drift time)
+ TString plane_dt_name;
+ TString plane_dt_title;
+
+ TString plane_dt_name_corr;
+ TString plane_dt_title_corr;
+
+
+ TString cell_dt_name;
+ TString cell_dt_title;
+
+ TString fitted_cell_dt_name;
+ TString fitted_cell_dt_title;
+
+ TString dt_vs_wire_name;
+ TString dt_vs_wire_title;
+
+ TString ifile_name;
+ TString ofile_name;
+
+ TString itxtfile_name;
+ TString otxtfile_name;
+
+ //Declare variables to store histograms
+ TH1F *plane_dt;
+ TH1F *plane_dt_corr;
+ TH1F **cell_dt;
+ TH1F **cell_dt_corr;
+ TH1F **fitted_cell_dt;
+ TH2F *dt_vs_wire;
+ TH2F *dt_vs_wire_corr;
+
+ //Read/Write ROOT files
+ TFile *in_file;
+ TFile *out_file;
+
+ //Read/Write txt files
+ ifstream in_txtFILE;
+ ofstream out_txtFILE;
+
+ TDirectory *main_dir;
+
+
+ //Declare 2d dynamic arrays to get histogram bin properties
+
+ Int_t **bin_max;
+ Int_t **bin_maxContent;
+ Double_t **time_max;
+ Double_t **twenty_perc_maxContent;
+ Double_t **ref_time;
+
+ //variables to be used in loop over bins for wire drift time
+ Int_t content_bin; //stores content for each bin
+ Int_t counts; //a counter used to count the number of bins that have >20% max bin content for a plane
+
+ vector<Int_t> content; //stores bin content
+ vector <Int_t> bin_num; //stores bin number
+
+
+ //Declare 'FIT' related variables
+ Int_t **entries; //wire drift time histo entries
+ Int_t binx;
+ Double_t time_init; //start fit value
+ Double_t time_final; //end fit value
+ TF1 *tZero_fit; //linear fit function
+
+
+ Double_t m; //slope
+ Double_t y_int; //y-intercept
+ Double_t m_err;
+ Double_t y_int_err;
+ Double_t std_dev;
+ Double_t **t_zero;
+ Double_t **t_zero_err;
+
+ //declare variables to make plot of tzero v. wire number
+ TGraphErrors *graph;
+ TString graph_title;
+ TCanvas *gr1_canv;
+
+ //Declare varibales to write LoookUp Table
+ Int_t bin_count;
+ Int_t bin_t0[NPLANES];
+ Int_t bin_final[NPLANES];
+ Int_t bin_Content[NPLANES];
+ Double_t binContent_TOTAL[NPLANES];
+ Double_t binSUM[NPLANES];
+ Double_t lookup_value[NPLANES];
+ TString lookup_table;
+ TString headers;
+
+
+
+};
+
+
+
+
+
+
+
+
+
+
+#endif //DC_CALIBRATION_H
diff --git a/CALIBRATION/shms_dc_calib/scripts/main_calib.C b/CALIBRATION/shms_dc_calib/scripts/main_calib.C
new file mode 100644
index 00000000..69b038aa
--- /dev/null
+++ b/CALIBRATION/shms_dc_calib/scripts/main_calib.C
@@ -0,0 +1,42 @@
+//Main Calibration Code
+#include "DC_calib.h"
+#include "DC_calib.C"
+#include <iostream>
+#include <ctime>
+using namespace std;
+
+int main_calib()
+{
+
+ //prevent root from displaying graphs while executing
+ gROOT->SetBatch(1);
+
+
+ //measure execution time
+ clock_t cl;
+ cl = clock();
+
+
+ DC_calib obj("SHMS", "../../../ROOTfiles/shms_replay_484_dc_uncal.root ", 484,300000);
+
+
+ obj.printInitVar();
+ obj.SetPlaneNames();
+ obj.GetDCLeafs();
+ obj.AllocateDynamicArrays();
+ obj.CreateHistoNames();
+ obj.EventLoop();
+ obj.Calculate_tZero();
+ obj.WriteTZeroParam();
+ obj.ApplyTZeroCorrection();
+ obj.WriteLookUpTable();
+ obj.WriteToFile(1); //set argument to (1) for debugging
+
+
+
+ //stop clock
+ cl = clock() - cl;
+ cout << "execution time: " << cl/(double)CLOCKS_PER_SEC << " sec" << endl;
+
+ return 0;
+}
diff --git a/PARAM/SHMS/DC/CALIB/pdc_calib.param b/PARAM/SHMS/DC/CALIB/pdc_calib.param
index 037c1375..a591fab9 100644
--- a/PARAM/SHMS/DC/CALIB/pdc_calib.param
+++ b/PARAM/SHMS/DC/CALIB/pdc_calib.param
@@ -1,175 +1,175 @@
-; Lookup Table: RUN 488
+; Lookup Table: RUN 484
; number of bins in time to distance lookup table
pdriftbins = 275
; number of 1st bin in table in ns
pdrift1stbin=0
; bin size in ns
pdriftbinsz=1
-pwc1u1fract=0.0011495,0.00243,0.00366,0.00539,0.00743,0.00975,0.01249,0.01605,0.02001,0.02453,0.02972,0.03511,0.04078,0.04675,0.05374,0.06102,0.06878,0.07687,0.08515
-0.09393,0.10314,0.11256,0.12215,0.13197,0.14184,0.15209,0.16243,0.17325,0.18372,0.19456,0.20521,0.21576,0.22655,0.23767,0.24822,0.25912,0.27040,0.28116,0.29233
-0.30330,0.31426,0.32525,0.33599,0.34687,0.35785,0.36891,0.38022,0.39096,0.40197,0.41297,0.42372,0.43471,0.44537,0.45641,0.46706,0.47768,0.48834,0.49927,0.51002
-0.52094,0.53162,0.54188,0.55271,0.56320,0.57384,0.58379,0.59435,0.60486,0.61506,0.62553,0.63582,0.64587,0.65548,0.66522,0.67533,0.68518,0.69477,0.70445,0.71410
-0.72364,0.73299,0.74193,0.75080,0.75949,0.76812,0.77676,0.78512,0.79353,0.80197,0.81042,0.81850,0.82631,0.83400,0.84162,0.84917,0.85645,0.86350,0.87029,0.87680
-0.88293,0.88865,0.89424,0.89980,0.90490,0.90965,0.91441,0.91866,0.92255,0.92590,0.92947,0.93273,0.93577,0.93842,0.94105,0.94333,0.94559,0.94773,0.94988,0.95189
-0.95370,0.95533,0.95691,0.95840,0.95993,0.96128,0.96258,0.96374,0.96505,0.96615,0.96734,0.96832,0.96950,0.97040,0.97134,0.97214,0.97306,0.97392,0.97475,0.97562
-0.97630,0.97705,0.97780,0.97853,0.97916,0.97981,0.98040,0.98103,0.98154,0.98210,0.98272,0.98323,0.98374,0.98421,0.98470,0.98519,0.98564,0.98606,0.98646,0.98683
-0.98715,0.98749,0.98781,0.98810,0.98840,0.98869,0.98896,0.98921,0.98946,0.98967,0.98994,0.99014,0.99031,0.99051,0.99068,0.99082,0.99099,0.99116,0.99136,0.99150
-0.99163,0.99184,0.99198,0.99211,0.99222,0.99230,0.99245,0.99262,0.99279,0.99290,0.99303,0.99317,0.99331,0.99340,0.99352,0.99361,0.99373,0.99382,0.99390,0.99398
-0.99409,0.99417,0.99429,0.99439,0.99447,0.99455,0.99468,0.99480,0.99491,0.99497,0.99506,0.99520,0.99530,0.99541,0.99550,0.99562,0.99572,0.99583,0.99594,0.99603
-0.99615,0.99625,0.99638,0.99650,0.99659,0.99669,0.99676,0.99687,0.99697,0.99710,0.99721,0.99737,0.99747,0.99756,0.99766,0.99776,0.99781,0.99791,0.99800,0.99808
-0.99815,0.99824,0.99837,0.99844,0.99854,0.99862,0.99871,0.99883,0.99893,0.99905,0.99911,0.99920,0.99926,0.99934,0.99940,0.99946,0.99951,0.99962,0.99969,0.99975
-0.99979,0.99983,0.99987,0.99989,0.99991,0.99993,0.99996,0.99997,0.99998,0.99998,0.99998,0.99999,0.99999,0.99999,1.00000,1.00000
-pwc1u2fract=0.00135,0.00257,0.00417,0.00588,0.00782,0.01023,0.01301,0.01658,0.02034,0.02486,0.02967,0.03505,0.04080,0.04737,0.05406,0.06141,0.06910,0.07711,0.08556
-0.09444,0.10369,0.11351,0.12309,0.13328,0.14365,0.15429,0.16460,0.17538,0.18622,0.19716,0.20808,0.21936,0.23026,0.24108,0.25223,0.26330,0.27404,0.28498,0.29598
-0.30734,0.31848,0.32955,0.34039,0.35140,0.36226,0.37284,0.38416,0.39539,0.40603,0.41663,0.42723,0.43832,0.44879,0.45920,0.46992,0.48049,0.49108,0.50160,0.51229
-0.52275,0.53305,0.54350,0.55381,0.56385,0.57419,0.58418,0.59443,0.60429,0.61436,0.62408,0.63407,0.64365,0.65340,0.66319,0.67276,0.68242,0.69184,0.70126,0.71042
-0.71938,0.72855,0.73786,0.74686,0.75577,0.76456,0.77351,0.78218,0.79041,0.79855,0.80691,0.81506,0.82298,0.83061,0.83774,0.84543,0.85242,0.85973,0.86620,0.87232
-0.87850,0.88448,0.89045,0.89625,0.90149,0.90641,0.91118,0.91552,0.91946,0.92342,0.92687,0.93024,0.93358,0.93651,0.93931,0.94172,0.94409,0.94637,0.94846,0.95042
-0.95248,0.95418,0.95581,0.95735,0.95909,0.96045,0.96179,0.96295,0.96431,0.96552,0.96667,0.96784,0.96888,0.96991,0.97076,0.97166,0.97253,0.97332,0.97412,0.97489
-0.97567,0.97639,0.97701,0.97767,0.97837,0.97910,0.97975,0.98036,0.98103,0.98160,0.98211,0.98263,0.98315,0.98362,0.98411,0.98449,0.98494,0.98544,0.98585,0.98620
-0.98653,0.98688,0.98716,0.98751,0.98781,0.98812,0.98839,0.98862,0.98885,0.98912,0.98936,0.98959,0.98982,0.99002,0.99022,0.99043,0.99064,0.99080,0.99094,0.99111
-0.99128,0.99144,0.99161,0.99175,0.99190,0.99207,0.99224,0.99241,0.99252,0.99261,0.99273,0.99285,0.99297,0.99308,0.99319,0.99328,0.99344,0.99356,0.99367,0.99374
-0.99382,0.99396,0.99405,0.99415,0.99424,0.99435,0.99447,0.99455,0.99461,0.99472,0.99482,0.99493,0.99505,0.99516,0.99526,0.99533,0.99545,0.99560,0.99573,0.99584
-0.99591,0.99603,0.99616,0.99626,0.99640,0.99650,0.99663,0.99677,0.99691,0.99703,0.99710,0.99721,0.99730,0.99740,0.99752,0.99758,0.99768,0.99777,0.99787,0.99797
-0.99809,0.99817,0.99826,0.99834,0.99839,0.99851,0.99864,0.99873,0.99884,0.99891,0.99903,0.99912,0.99921,0.99929,0.99939,0.99948,0.99953,0.99962,0.99969,0.99977
-0.99982,0.99988,0.99990,0.99993,0.99995,0.99996,0.99996,0.99997,0.99997,0.99998,0.99999,0.99999,1.00000,1.00000,1.00000,1.00000
-pwc1x1fract=0.00111,0.00229,0.00364,0.00534,0.00739,0.00990,0.01320,0.01695,0.02106,0.02575,0.03089,0.03672,0.04270,0.04964,0.05674,0.06443,0.07242,0.08111,0.09000
-0.09936,0.10918,0.11894,0.12883,0.13917,0.14968,0.16016,0.17122,0.18151,0.19208,0.20313,0.21429,0.22524,0.23627,0.24719,0.25821,0.26924,0.28004,0.29079,0.30171
-0.31238,0.32303,0.33408,0.34483,0.35585,0.36661,0.37748,0.38804,0.39881,0.40944,0.42020,0.43094,0.44189,0.45228,0.46292,0.47366,0.48450,0.49485,0.50548,0.51575
-0.52606,0.53662,0.54701,0.55722,0.56739,0.57748,0.58777,0.59813,0.60804,0.61831,0.62822,0.63784,0.64737,0.65703,0.66647,0.67604,0.68520,0.69452,0.70386,0.71287
-0.72198,0.73125,0.73984,0.74863,0.75716,0.76560,0.77425,0.78259,0.79077,0.79929,0.80753,0.81535,0.82336,0.83134,0.83909,0.84642,0.85357,0.86012,0.86653,0.87314
-0.87913,0.88470,0.89008,0.89517,0.90039,0.90508,0.90948,0.91363,0.91767,0.92112,0.92461,0.92795,0.93083,0.93366,0.93630,0.93873,0.94135,0.94369,0.94563,0.94757
-0.94933,0.95105,0.95271,0.95441,0.95587,0.95726,0.95865,0.95995,0.96120,0.96245,0.96371,0.96492,0.96601,0.96707,0.96816,0.96911,0.97010,0.97097,0.97203,0.97295
-0.97375,0.97463,0.97545,0.97611,0.97690,0.97763,0.97835,0.97905,0.97977,0.98037,0.98099,0.98155,0.98211,0.98267,0.98320,0.98370,0.98412,0.98458,0.98502,0.98540
-0.98578,0.98619,0.98657,0.98692,0.98723,0.98756,0.98785,0.98815,0.98842,0.98873,0.98897,0.98915,0.98942,0.98964,0.98985,0.99008,0.99030,0.99048,0.99065,0.99082
-0.99104,0.99120,0.99136,0.99153,0.99165,0.99180,0.99191,0.99204,0.99221,0.99235,0.99249,0.99265,0.99277,0.99287,0.99299,0.99311,0.99322,0.99332,0.99347,0.99362
-0.99374,0.99385,0.99395,0.99407,0.99419,0.99429,0.99436,0.99448,0.99460,0.99472,0.99484,0.99494,0.99507,0.99518,0.99530,0.99540,0.99554,0.99567,0.99577,0.99585
-0.99596,0.99607,0.99616,0.99623,0.99633,0.99642,0.99651,0.99662,0.99672,0.99679,0.99686,0.99702,0.99712,0.99722,0.99728,0.99742,0.99753,0.99761,0.99774,0.99785
-0.99793,0.99805,0.99815,0.99828,0.99840,0.99850,0.99858,0.99868,0.99878,0.99893,0.99901,0.99909,0.99917,0.99925,0.99931,0.99937,0.99946,0.99957,0.99962,0.99969
-0.99974,0.99979,0.99981,0.99984,0.99987,0.99990,0.99992,0.99994,0.99996,0.99997,0.99997,0.99998,0.99998,0.99999,1.00000,1.00000
-pwc1x2fract=0.00109,0.00232,0.00378,0.00558,0.00777,0.01056,0.01381,0.01788,0.02227,0.02741,0.03302,0.03935,0.04573,0.05274,0.06074,0.06884,0.07767,0.08666,0.09589
-0.10544,0.11522,0.12522,0.13528,0.14545,0.15591,0.16676,0.17736,0.18800,0.19852,0.20900,0.21992,0.23080,0.24178,0.25252,0.26315,0.27373,0.28468,0.29566,0.30644
-0.31731,0.32815,0.33904,0.34992,0.36067,0.37137,0.38205,0.39290,0.40348,0.41405,0.42453,0.43512,0.44584,0.45607,0.46670,0.47717,0.48769,0.49812,0.50821,0.51841
-0.52857,0.53883,0.54912,0.55959,0.56973,0.57991,0.58992,0.59987,0.60954,0.61916,0.62871,0.63842,0.64826,0.65752,0.66693,0.67628,0.68521,0.69432,0.70363,0.71290
-0.72184,0.73072,0.73937,0.74775,0.75667,0.76534,0.77373,0.78221,0.79063,0.79879,0.80672,0.81458,0.82228,0.82966,0.83721,0.84453,0.85112,0.85813,0.86474,0.87104
-0.87696,0.88273,0.88859,0.89376,0.89886,0.90377,0.90809,0.91259,0.91667,0.92022,0.92389,0.92714,0.93024,0.93313,0.93587,0.93848,0.94090,0.94295,0.94504,0.94689
-0.94857,0.95012,0.95159,0.95313,0.95463,0.95607,0.95732,0.95872,0.95994,0.96114,0.96234,0.96347,0.96465,0.96577,0.96696,0.96799,0.96890,0.96985,0.97075,0.97162
-0.97249,0.97334,0.97417,0.97504,0.97579,0.97654,0.97731,0.97798,0.97861,0.97925,0.97994,0.98061,0.98122,0.98175,0.98229,0.98279,0.98328,0.98378,0.98422,0.98463
-0.98514,0.98559,0.98608,0.98650,0.98689,0.98728,0.98755,0.98789,0.98811,0.98842,0.98870,0.98893,0.98917,0.98941,0.98965,0.98985,0.99005,0.99023,0.99048,0.99066
-0.99088,0.99108,0.99127,0.99141,0.99157,0.99173,0.99188,0.99205,0.99220,0.99232,0.99245,0.99260,0.99271,0.99284,0.99297,0.99309,0.99324,0.99335,0.99347,0.99360
-0.99374,0.99386,0.99399,0.99413,0.99424,0.99434,0.99442,0.99451,0.99464,0.99475,0.99489,0.99500,0.99509,0.99524,0.99533,0.99545,0.99562,0.99573,0.99586,0.99594
-0.99603,0.99617,0.99630,0.99640,0.99650,0.99659,0.99674,0.99686,0.99690,0.99702,0.99711,0.99720,0.99729,0.99741,0.99751,0.99758,0.99770,0.99779,0.99788,0.99797
-0.99805,0.99812,0.99821,0.99833,0.99846,0.99854,0.99862,0.99874,0.99882,0.99892,0.99901,0.99911,0.99918,0.99929,0.99935,0.99942,0.99948,0.99955,0.99963,0.99973
-0.99977,0.99984,0.99987,0.99989,0.99992,0.99993,0.99995,0.99996,0.99998,0.99998,0.99999,1.00000,1.00000,1.00000,1.00000,1.00000
-pwc1v1fract=0.00102,0.00222,0.00354,0.00521,0.00721,0.00963,0.01257,0.01620,0.02021,0.02497,0.02989,0.03540,0.04122,0.04771,0.05442,0.06172,0.06926,0.07741,0.08595
-0.09498,0.10438,0.11392,0.12349,0.13377,0.14422,0.15501,0.16564,0.17648,0.18728,0.19804,0.20931,0.22017,0.23112,0.24210,0.25314,0.26456,0.27552,0.28655,0.29772
-0.30856,0.31991,0.33115,0.34208,0.35303,0.36408,0.37494,0.38606,0.39708,0.40797,0.41878,0.42946,0.44033,0.45111,0.46201,0.47288,0.48349,0.49430,0.50466,0.51483
-0.52505,0.53585,0.54633,0.55670,0.56686,0.57694,0.58703,0.59717,0.60736,0.61708,0.62744,0.63721,0.64739,0.65664,0.66657,0.67619,0.68551,0.69488,0.70414,0.71341
-0.72276,0.73184,0.74090,0.74996,0.75860,0.76712,0.77566,0.78394,0.79236,0.80091,0.80904,0.81677,0.82451,0.83207,0.83928,0.84633,0.85345,0.86044,0.86717,0.87345
-0.87937,0.88520,0.89103,0.89640,0.90134,0.90614,0.91049,0.91454,0.91856,0.92234,0.92576,0.92903,0.93221,0.93498,0.93747,0.94003,0.94245,0.94467,0.94686,0.94889
-0.95082,0.95266,0.95429,0.95587,0.95732,0.95875,0.96006,0.96125,0.96259,0.96373,0.96489,0.96616,0.96716,0.96823,0.96928,0.97024,0.97123,0.97214,0.97299,0.97384
-0.97456,0.97536,0.97609,0.97677,0.97752,0.97816,0.97882,0.97943,0.98004,0.98063,0.98126,0.98191,0.98246,0.98300,0.98351,0.98398,0.98446,0.98494,0.98534,0.98574
-0.98608,0.98645,0.98681,0.98720,0.98751,0.98781,0.98813,0.98844,0.98870,0.98894,0.98921,0.98941,0.98963,0.98981,0.99002,0.99021,0.99044,0.99065,0.99077,0.99092
-0.99103,0.99121,0.99140,0.99153,0.99166,0.99180,0.99194,0.99205,0.99219,0.99231,0.99246,0.99258,0.99268,0.99280,0.99291,0.99302,0.99313,0.99331,0.99341,0.99355
-0.99367,0.99377,0.99390,0.99401,0.99409,0.99421,0.99435,0.99448,0.99457,0.99468,0.99478,0.99487,0.99496,0.99512,0.99519,0.99526,0.99537,0.99551,0.99562,0.99573
-0.99585,0.99595,0.99607,0.99618,0.99627,0.99636,0.99646,0.99651,0.99663,0.99676,0.99685,0.99697,0.99704,0.99716,0.99730,0.99739,0.99745,0.99759,0.99772,0.99781
-0.99799,0.99809,0.99819,0.99829,0.99840,0.99849,0.99863,0.99874,0.99883,0.99892,0.99899,0.99911,0.99922,0.99932,0.99941,0.99947,0.99954,0.99960,0.99968,0.99973
-0.99980,0.99982,0.99986,0.99989,0.99990,0.99992,0.99995,0.99996,0.99998,0.99999,1.00000,1.00000,1.00000,1.00000,1.00000,1.00000
-pwc1v2fract=0.00104,0.00216,0.00353,0.00519,0.00721,0.00957,0.01257,0.01632,0.02040,0.02507,0.03017,0.03585,0.04185,0.04844,0.05567,0.06312,0.07116,0.07955,0.08812
-0.09699,0.10653,0.11605,0.12623,0.13608,0.14623,0.15688,0.16758,0.17847,0.18991,0.20037,0.21117,0.22214,0.23287,0.24406,0.25503,0.26607,0.27701,0.28793,0.29904
-0.30952,0.32053,0.33162,0.34266,0.35395,0.36482,0.37570,0.38622,0.39701,0.40767,0.41862,0.42919,0.43999,0.45054,0.46094,0.47174,0.48234,0.49254,0.50277,0.51337
-0.52387,0.53426,0.54453,0.55481,0.56500,0.57522,0.58545,0.59523,0.60521,0.61511,0.62517,0.63476,0.64466,0.65424,0.66369,0.67361,0.68348,0.69298,0.70245,0.71176
-0.72108,0.73023,0.73942,0.74831,0.75729,0.76595,0.77476,0.78349,0.79183,0.80035,0.80852,0.81668,0.82471,0.83234,0.83974,0.84692,0.85381,0.86095,0.86774,0.87435
-0.88053,0.88662,0.89205,0.89745,0.90250,0.90747,0.91195,0.91588,0.91988,0.92385,0.92726,0.93033,0.93349,0.93639,0.93905,0.94136,0.94369,0.94590,0.94786,0.94975
-0.95169,0.95333,0.95488,0.95638,0.95782,0.95912,0.96049,0.96173,0.96296,0.96404,0.96518,0.96638,0.96739,0.96836,0.96927,0.97016,0.97108,0.97201,0.97285,0.97371
-0.97451,0.97526,0.97603,0.97672,0.97741,0.97800,0.97862,0.97935,0.98004,0.98062,0.98117,0.98171,0.98229,0.98275,0.98320,0.98375,0.98420,0.98466,0.98507,0.98543
-0.98585,0.98620,0.98651,0.98682,0.98716,0.98748,0.98777,0.98804,0.98837,0.98862,0.98886,0.98907,0.98927,0.98951,0.98974,0.98991,0.99009,0.99027,0.99045,0.99060
-0.99084,0.99104,0.99123,0.99140,0.99156,0.99170,0.99186,0.99195,0.99211,0.99224,0.99239,0.99251,0.99266,0.99280,0.99290,0.99301,0.99320,0.99332,0.99340,0.99358
-0.99372,0.99384,0.99396,0.99407,0.99416,0.99425,0.99431,0.99440,0.99449,0.99458,0.99472,0.99485,0.99496,0.99506,0.99515,0.99525,0.99535,0.99543,0.99550,0.99563
-0.99577,0.99591,0.99598,0.99607,0.99617,0.99629,0.99637,0.99650,0.99662,0.99673,0.99682,0.99692,0.99703,0.99711,0.99724,0.99733,0.99744,0.99753,0.99769,0.99785
-0.99795,0.99803,0.99812,0.99818,0.99829,0.99839,0.99846,0.99854,0.99866,0.99876,0.99885,0.99897,0.99903,0.99912,0.99922,0.99931,0.99941,0.99949,0.99956,0.99965
-0.99970,0.99974,0.99981,0.99985,0.99989,0.99992,0.99994,0.99995,0.99996,0.99998,0.99999,0.99999,0.99999,1.00000,1.00000,1.00000
-pwc2v2fract=0.00109,0.00224,0.00372,0.00522,0.00718,0.00949,0.01230,0.01586,0.01952,0.02354,0.02831,0.03381,0.03958,0.04599,0.05261,0.06014,0.06789,0.07617,0.08486
-0.09433,0.10403,0.11409,0.12365,0.13417,0.14476,0.15531,0.16643,0.17723,0.18815,0.19907,0.21005,0.22104,0.23207,0.24285,0.25399,0.26480,0.27597,0.28700,0.29795
-0.30868,0.31962,0.33050,0.34158,0.35235,0.36281,0.37376,0.38454,0.39523,0.40607,0.41654,0.42744,0.43801,0.44883,0.45948,0.47004,0.48066,0.49077,0.50091,0.51140
-0.52166,0.53180,0.54187,0.55183,0.56192,0.57156,0.58116,0.59099,0.60080,0.60994,0.61945,0.62882,0.63811,0.64761,0.65678,0.66565,0.67467,0.68368,0.69236,0.70114
-0.70948,0.71783,0.72629,0.73447,0.74276,0.75079,0.75887,0.76663,0.77472,0.78263,0.79018,0.79753,0.80479,0.81182,0.81880,0.82572,0.83241,0.83903,0.84530,0.85125
-0.85700,0.86252,0.86821,0.87351,0.87864,0.88375,0.88846,0.89294,0.89730,0.90148,0.90509,0.90875,0.91233,0.91566,0.91895,0.92185,0.92473,0.92751,0.93007,0.93238
-0.93473,0.93694,0.93907,0.94122,0.94330,0.94519,0.94698,0.94873,0.95034,0.95188,0.95342,0.95489,0.95623,0.95766,0.95889,0.96015,0.96139,0.96242,0.96348,0.96467
-0.96578,0.96670,0.96762,0.96847,0.96933,0.97018,0.97099,0.97190,0.97279,0.97352,0.97426,0.97499,0.97566,0.97636,0.97692,0.97751,0.97804,0.97865,0.97931,0.97997
-0.98052,0.98106,0.98154,0.98206,0.98256,0.98306,0.98350,0.98392,0.98432,0.98465,0.98501,0.98537,0.98576,0.98615,0.98644,0.98672,0.98703,0.98733,0.98762,0.98792
-0.98818,0.98846,0.98878,0.98906,0.98929,0.98949,0.98964,0.98987,0.99005,0.99027,0.99052,0.99076,0.99101,0.99122,0.99140,0.99159,0.99181,0.99199,0.99218,0.99233
-0.99253,0.99269,0.99284,0.99300,0.99313,0.99325,0.99343,0.99354,0.99368,0.99382,0.99397,0.99412,0.99425,0.99437,0.99449,0.99467,0.99482,0.99497,0.99509,0.99521
-0.99533,0.99547,0.99563,0.99574,0.99585,0.99603,0.99615,0.99627,0.99640,0.99650,0.99661,0.99669,0.99684,0.99695,0.99709,0.99724,0.99738,0.99747,0.99760,0.99774
-0.99782,0.99797,0.99807,0.99821,0.99832,0.99844,0.99857,0.99868,0.99879,0.99890,0.99900,0.99911,0.99923,0.99934,0.99942,0.99950,0.99960,0.99964,0.99971,0.99978
-0.99981,0.99984,0.99987,0.99991,0.99993,0.99994,0.99995,0.99996,0.99997,0.99998,1.00000,1.00000,1.00000,1.00000,1.00000,1.00000
-pwc2v1fract=0.00113,0.00252,0.00395,0.00559,0.00760,0.01001,0.01289,0.01637,0.02044,0.02489,0.02979,0.03519,0.04105,0.04776,0.05450,0.06188,0.06951,0.07765,0.08650
-0.09544,0.10509,0.11483,0.12515,0.13550,0.14591,0.15636,0.16715,0.17819,0.18908,0.19994,0.21079,0.22212,0.23340,0.24479,0.25579,0.26711,0.27849,0.28935,0.29993
-0.31083,0.32220,0.33338,0.34476,0.35563,0.36644,0.37735,0.38778,0.39853,0.40918,0.41978,0.43039,0.44107,0.45183,0.46222,0.47259,0.48281,0.49347,0.50384,0.51407
-0.52424,0.53438,0.54477,0.55496,0.56470,0.57467,0.58483,0.59443,0.60423,0.61415,0.62357,0.63302,0.64205,0.65105,0.66052,0.66949,0.67818,0.68745,0.69638,0.70543
-0.71391,0.72216,0.73035,0.73858,0.74674,0.75503,0.76291,0.77072,0.77861,0.78656,0.79405,0.80147,0.80875,0.81567,0.82266,0.82953,0.83614,0.84257,0.84893,0.85509
-0.86085,0.86664,0.87208,0.87749,0.88255,0.88721,0.89209,0.89679,0.90106,0.90527,0.90911,0.91260,0.91582,0.91894,0.92194,0.92488,0.92762,0.93025,0.93275,0.93506
-0.93741,0.93965,0.94163,0.94373,0.94552,0.94729,0.94907,0.95070,0.95210,0.95368,0.95506,0.95654,0.95791,0.95912,0.96038,0.96156,0.96272,0.96399,0.96515,0.96620
-0.96717,0.96809,0.96903,0.97004,0.97087,0.97175,0.97262,0.97335,0.97406,0.97488,0.97564,0.97637,0.97702,0.97762,0.97822,0.97884,0.97939,0.97998,0.98047,0.98101
-0.98157,0.98204,0.98254,0.98296,0.98344,0.98389,0.98430,0.98474,0.98518,0.98557,0.98597,0.98628,0.98664,0.98692,0.98726,0.98757,0.98790,0.98820,0.98848,0.98872
-0.98899,0.98921,0.98942,0.98966,0.98993,0.99014,0.99036,0.99056,0.99076,0.99093,0.99113,0.99128,0.99144,0.99158,0.99169,0.99185,0.99201,0.99216,0.99240,0.99255
-0.99275,0.99290,0.99301,0.99315,0.99332,0.99345,0.99361,0.99377,0.99390,0.99408,0.99420,0.99436,0.99445,0.99460,0.99473,0.99488,0.99497,0.99510,0.99525,0.99540
-0.99552,0.99567,0.99579,0.99590,0.99603,0.99612,0.99624,0.99635,0.99645,0.99655,0.99666,0.99677,0.99690,0.99701,0.99713,0.99723,0.99735,0.99747,0.99759,0.99769
-0.99786,0.99798,0.99809,0.99819,0.99829,0.99840,0.99857,0.99869,0.99880,0.99894,0.99905,0.99913,0.99921,0.99930,0.99941,0.99950,0.99959,0.99967,0.99974,0.99978
-0.99983,0.99988,0.99990,0.99991,0.99992,0.99994,0.99996,0.99996,0.99998,0.99998,0.99998,0.99999,1.00000,1.00000,1.00000,1.00000
-pwc2x2fract=0.00118,0.00252,0.00399,0.00571,0.00792,0.01034,0.01349,0.01725,0.02146,0.02619,0.03152,0.03763,0.04425,0.05116,0.05840,0.06629,0.07474,0.08351,0.09265
-0.10239,0.11238,0.12243,0.13260,0.14300,0.15325,0.16376,0.17467,0.18538,0.19630,0.20729,0.21823,0.22922,0.24009,0.25104,0.26203,0.27288,0.28386,0.29464,0.30537
-0.31580,0.32616,0.33708,0.34787,0.35878,0.36935,0.37969,0.39023,0.40086,0.41177,0.42230,0.43277,0.44322,0.45379,0.46439,0.47485,0.48532,0.49549,0.50570,0.51602
-0.52646,0.53620,0.54607,0.55600,0.56572,0.57586,0.58543,0.59516,0.60474,0.61427,0.62378,0.63327,0.64274,0.65175,0.66074,0.66972,0.67843,0.68744,0.69655,0.70513
-0.71389,0.72216,0.73070,0.73914,0.74757,0.75585,0.76388,0.77183,0.77950,0.78737,0.79521,0.80263,0.80981,0.81700,0.82405,0.83111,0.83791,0.84446,0.85110,0.85741
-0.86344,0.86929,0.87531,0.88078,0.88603,0.89090,0.89566,0.89996,0.90437,0.90841,0.91194,0.91550,0.91866,0.92166,0.92455,0.92722,0.92969,0.93200,0.93406,0.93619
-0.93815,0.94005,0.94190,0.94359,0.94547,0.94711,0.94862,0.95000,0.95144,0.95278,0.95417,0.95546,0.95672,0.95786,0.95914,0.96027,0.96133,0.96259,0.96366,0.96470
-0.96571,0.96674,0.96772,0.96868,0.96964,0.97054,0.97145,0.97232,0.97319,0.97399,0.97467,0.97540,0.97616,0.97689,0.97761,0.97822,0.97894,0.97962,0.98013,0.98077
-0.98131,0.98186,0.98239,0.98292,0.98344,0.98390,0.98429,0.98471,0.98512,0.98555,0.98588,0.98624,0.98663,0.98695,0.98723,0.98755,0.98793,0.98818,0.98846,0.98878
-0.98899,0.98925,0.98950,0.98974,0.98996,0.99022,0.99046,0.99063,0.99083,0.99104,0.99122,0.99140,0.99161,0.99176,0.99193,0.99210,0.99226,0.99244,0.99259,0.99275
-0.99290,0.99303,0.99321,0.99333,0.99348,0.99362,0.99375,0.99388,0.99400,0.99410,0.99424,0.99436,0.99450,0.99464,0.99475,0.99484,0.99497,0.99507,0.99524,0.99540
-0.99553,0.99567,0.99578,0.99585,0.99599,0.99612,0.99624,0.99635,0.99644,0.99655,0.99666,0.99678,0.99688,0.99702,0.99715,0.99727,0.99738,0.99749,0.99761,0.99775
-0.99784,0.99794,0.99804,0.99818,0.99831,0.99843,0.99853,0.99864,0.99872,0.99882,0.99892,0.99900,0.99911,0.99918,0.99927,0.99941,0.99949,0.99957,0.99962,0.99971
-0.99976,0.99981,0.99984,0.99991,0.99993,0.99994,0.99997,0.99998,0.99998,0.99999,0.99999,1.00000,1.00000,1.00000,1.00000,1.00000
-pwc2x1fract=0.00108,0.00227,0.00356,0.00532,0.00752,0.00994,0.01286,0.01659,0.02081,0.02538,0.03048,0.03620,0.04250,0.04919,0.05642,0.06419,0.07233,0.08098,0.09000
-0.09961,0.10909,0.11877,0.12884,0.13897,0.14942,0.15967,0.17030,0.18060,0.19139,0.20228,0.21314,0.22352,0.23397,0.24493,0.25559,0.26598,0.27663,0.28758,0.29814
-0.30903,0.32010,0.33090,0.34093,0.35145,0.36221,0.37273,0.38306,0.39375,0.40451,0.41477,0.42552,0.43603,0.44645,0.45683,0.46705,0.47717,0.48769,0.49812,0.50857
-0.51849,0.52898,0.53902,0.54910,0.55936,0.56946,0.57907,0.58881,0.59858,0.60840,0.61824,0.62760,0.63700,0.64647,0.65579,0.66515,0.67452,0.68357,0.69294,0.70169
-0.71035,0.71905,0.72786,0.73641,0.74510,0.75359,0.76179,0.76978,0.77764,0.78534,0.79333,0.80100,0.80871,0.81608,0.82351,0.83056,0.83755,0.84441,0.85098,0.85724
-0.86387,0.86980,0.87554,0.88093,0.88617,0.89099,0.89606,0.90053,0.90490,0.90895,0.91259,0.91599,0.91931,0.92230,0.92518,0.92787,0.93045,0.93297,0.93526,0.93747
-0.93946,0.94148,0.94334,0.94512,0.94692,0.94860,0.95033,0.95194,0.95337,0.95465,0.95603,0.95722,0.95848,0.95964,0.96088,0.96210,0.96320,0.96427,0.96530,0.96622
-0.96723,0.96822,0.96920,0.97010,0.97096,0.97186,0.97273,0.97353,0.97437,0.97517,0.97593,0.97660,0.97735,0.97797,0.97863,0.97922,0.97983,0.98049,0.98108,0.98163
-0.98216,0.98265,0.98312,0.98361,0.98401,0.98441,0.98480,0.98520,0.98555,0.98595,0.98634,0.98665,0.98699,0.98733,0.98773,0.98800,0.98828,0.98852,0.98877,0.98904
-0.98929,0.98955,0.98979,0.99004,0.99028,0.99050,0.99073,0.99093,0.99117,0.99139,0.99156,0.99172,0.99188,0.99202,0.99218,0.99235,0.99250,0.99267,0.99283,0.99295
-0.99310,0.99321,0.99336,0.99348,0.99363,0.99378,0.99392,0.99407,0.99424,0.99440,0.99455,0.99467,0.99477,0.99489,0.99502,0.99515,0.99524,0.99537,0.99548,0.99562
-0.99575,0.99585,0.99596,0.99604,0.99618,0.99627,0.99638,0.99649,0.99661,0.99669,0.99682,0.99693,0.99704,0.99715,0.99726,0.99734,0.99746,0.99759,0.99772,0.99781
-0.99790,0.99803,0.99816,0.99834,0.99844,0.99857,0.99866,0.99874,0.99881,0.99895,0.99904,0.99912,0.99921,0.99932,0.99942,0.99948,0.99955,0.99965,0.99971,0.99976
-0.99982,0.99983,0.99986,0.99990,0.99992,0.99995,0.99997,0.99999,0.99999,1.00000,1.00000,1.00000,1.00000,1.00000,1.00000,1.00000
-pwc2u2fract=0.00126,0.00244,0.00386,0.00555,0.00768,0.01014,0.01294,0.01637,0.02041,0.02484,0.02991,0.03537,0.04170,0.04818,0.05489,0.06199,0.06986,0.07778,0.08642
-0.09535,0.10466,0.11410,0.12388,0.13375,0.14404,0.15456,0.16511,0.17545,0.18630,0.19721,0.20818,0.21876,0.22943,0.24038,0.25143,0.26231,0.27312,0.28409,0.29525
-0.30609,0.31695,0.32774,0.33850,0.34958,0.36048,0.37134,0.38209,0.39279,0.40322,0.41353,0.42439,0.43519,0.44569,0.45660,0.46741,0.47772,0.48795,0.49834,0.50876
-0.51911,0.52944,0.53974,0.55018,0.56005,0.57007,0.58027,0.59034,0.60014,0.61001,0.61974,0.62973,0.63901,0.64864,0.65804,0.66767,0.67673,0.68597,0.69511,0.70409
-0.71308,0.72202,0.73055,0.73929,0.74772,0.75612,0.76435,0.77265,0.78063,0.78862,0.79643,0.80391,0.81148,0.81888,0.82598,0.83314,0.84017,0.84687,0.85349,0.85997
-0.86636,0.87211,0.87771,0.88312,0.88811,0.89286,0.89751,0.90212,0.90634,0.91014,0.91396,0.91744,0.92077,0.92376,0.92670,0.92936,0.93192,0.93421,0.93660,0.93867
-0.94075,0.94263,0.94439,0.94632,0.94815,0.94987,0.95151,0.95304,0.95437,0.95568,0.95721,0.95860,0.95979,0.96104,0.96220,0.96337,0.96448,0.96559,0.96662,0.96761
-0.96863,0.96959,0.97058,0.97154,0.97239,0.97318,0.97401,0.97488,0.97557,0.97628,0.97701,0.97766,0.97825,0.97887,0.97955,0.98013,0.98072,0.98117,0.98168,0.98222
-0.98267,0.98314,0.98361,0.98406,0.98452,0.98496,0.98532,0.98575,0.98607,0.98636,0.98677,0.98707,0.98741,0.98774,0.98803,0.98835,0.98859,0.98883,0.98907,0.98934
-0.98959,0.98975,0.98995,0.99014,0.99032,0.99052,0.99068,0.99082,0.99099,0.99116,0.99132,0.99146,0.99162,0.99176,0.99188,0.99202,0.99216,0.99231,0.99246,0.99262
-0.99274,0.99286,0.99299,0.99313,0.99329,0.99341,0.99356,0.99372,0.99384,0.99393,0.99409,0.99427,0.99436,0.99450,0.99462,0.99475,0.99487,0.99500,0.99512,0.99524
-0.99534,0.99551,0.99561,0.99570,0.99579,0.99588,0.99601,0.99616,0.99628,0.99638,0.99650,0.99664,0.99676,0.99689,0.99700,0.99712,0.99725,0.99735,0.99746,0.99762
-0.99772,0.99784,0.99796,0.99808,0.99819,0.99831,0.99840,0.99851,0.99861,0.99872,0.99881,0.99895,0.99906,0.99916,0.99924,0.99935,0.99944,0.99952,0.99957,0.99966
-0.99973,0.99979,0.99983,0.99986,0.99991,0.99993,0.99995,0.99996,0.99998,0.99998,0.99999,1.00000,1.00000,1.00000,1.00000,1.00000
-pwc2u1fract=0.00106,0.00239,0.00373,0.00544,0.00738,0.00954,0.01230,0.01551,0.01934,0.02347,0.02813,0.03335,0.03906,0.04498,0.05132,0.05827,0.06588,0.07394,0.08235
-0.09080,0.09976,0.10927,0.11884,0.12834,0.13806,0.14824,0.15850,0.16913,0.17993,0.19062,0.20142,0.21197,0.22275,0.23360,0.24452,0.25539,0.26671,0.27745,0.28853
-0.29972,0.31068,0.32179,0.33235,0.34329,0.35412,0.36485,0.37574,0.38665,0.39749,0.40797,0.41844,0.42916,0.43944,0.45021,0.46090,0.47144,0.48163,0.49214,0.50248
-0.51306,0.52322,0.53324,0.54338,0.55333,0.56337,0.57327,0.58325,0.59348,0.60350,0.61342,0.62334,0.63284,0.64245,0.65206,0.66150,0.67109,0.68040,0.68956,0.69858
-0.70776,0.71664,0.72529,0.73418,0.74304,0.75147,0.76006,0.76847,0.77655,0.78480,0.79284,0.80082,0.80862,0.81628,0.82369,0.83068,0.83797,0.84498,0.85180,0.85825
-0.86462,0.87059,0.87620,0.88177,0.88706,0.89200,0.89693,0.90127,0.90561,0.90954,0.91330,0.91701,0.92044,0.92360,0.92667,0.92951,0.93214,0.93472,0.93697,0.93927
-0.94151,0.94354,0.94541,0.94724,0.94920,0.95094,0.95260,0.95408,0.95569,0.95723,0.95866,0.96000,0.96138,0.96265,0.96398,0.96527,0.96633,0.96734,0.96854,0.96943
-0.97043,0.97138,0.97222,0.97312,0.97389,0.97479,0.97550,0.97625,0.97704,0.97777,0.97840,0.97911,0.97976,0.98038,0.98099,0.98157,0.98207,0.98257,0.98306,0.98352
-0.98404,0.98445,0.98488,0.98525,0.98568,0.98600,0.98636,0.98667,0.98696,0.98724,0.98761,0.98786,0.98811,0.98835,0.98863,0.98888,0.98916,0.98943,0.98967,0.98985
-0.99011,0.99032,0.99054,0.99072,0.99094,0.99112,0.99129,0.99150,0.99166,0.99180,0.99196,0.99213,0.99227,0.99243,0.99258,0.99272,0.99288,0.99300,0.99313,0.99330
-0.99342,0.99358,0.99371,0.99382,0.99397,0.99409,0.99420,0.99431,0.99446,0.99459,0.99472,0.99483,0.99490,0.99499,0.99508,0.99525,0.99535,0.99543,0.99555,0.99564
-0.99575,0.99586,0.99597,0.99610,0.99618,0.99626,0.99636,0.99649,0.99661,0.99670,0.99685,0.99699,0.99706,0.99718,0.99729,0.99742,0.99749,0.99758,0.99768,0.99779
-0.99788,0.99797,0.99808,0.99819,0.99831,0.99845,0.99857,0.99867,0.99880,0.99894,0.99905,0.99914,0.99922,0.99931,0.99939,0.99947,0.99956,0.99963,0.99971,0.99976
-0.99981,0.99984,0.99989,0.99992,0.99993,0.99995,0.99996,0.99998,0.99999,0.99999,1.00000,1.00000,1.00000,1.00000,1.00000,1.00000
+pwc1u1fract=0.00063,0.00146,0.00245,0.00349,0.00494,0.00663,0.00852,0.01089,0.01359,0.01662,0.02025,0.02444,0.02920,0.03420,0.03994,0.04586,0.05230,0.05964,0.06750
+0.07551,0.08355,0.09250,0.10120,0.11053,0.11998,0.12956,0.13926,0.14934,0.15990,0.17005,0.18039,0.19118,0.20226,0.21317,0.22370,0.23490,0.24581,0.25682,0.26789
+0.27882,0.28980,0.30111,0.31220,0.32327,0.33429,0.34561,0.35717,0.36774,0.37876,0.38971,0.40060,0.41149,0.42282,0.43378,0.44466,0.45543,0.46629,0.47696,0.48768
+0.49824,0.50912,0.51982,0.53032,0.54091,0.55150,0.56171,0.57186,0.58233,0.59217,0.60251,0.61280,0.62295,0.63285,0.64292,0.65290,0.66268,0.67255,0.68224,0.69189
+0.70155,0.71125,0.72044,0.72964,0.73883,0.74788,0.75662,0.76576,0.77469,0.78308,0.79165,0.79981,0.80799,0.81576,0.82383,0.83150,0.83922,0.84668,0.85417,0.86095
+0.86740,0.87381,0.88001,0.88619,0.89170,0.89714,0.90233,0.90728,0.91190,0.91606,0.92006,0.92369,0.92726,0.93070,0.93352,0.93645,0.93916,0.94152,0.94386,0.94607
+0.94815,0.95023,0.95205,0.95359,0.95520,0.95664,0.95827,0.95974,0.96105,0.96222,0.96349,0.96472,0.96592,0.96702,0.96797,0.96905,0.96993,0.97094,0.97190,0.97280
+0.97372,0.97451,0.97540,0.97618,0.97681,0.97753,0.97817,0.97876,0.97948,0.98002,0.98050,0.98103,0.98162,0.98213,0.98261,0.98313,0.98362,0.98403,0.98443,0.98485
+0.98527,0.98567,0.98604,0.98638,0.98672,0.98704,0.98740,0.98769,0.98797,0.98817,0.98844,0.98867,0.98887,0.98909,0.98927,0.98947,0.98970,0.98992,0.99013,0.99032
+0.99051,0.99068,0.99082,0.99098,0.99112,0.99125,0.99139,0.99154,0.99170,0.99183,0.99197,0.99212,0.99225,0.99234,0.99248,0.99260,0.99271,0.99281,0.99291,0.99303
+0.99315,0.99327,0.99339,0.99348,0.99357,0.99368,0.99383,0.99395,0.99405,0.99419,0.99430,0.99440,0.99447,0.99461,0.99472,0.99481,0.99493,0.99505,0.99515,0.99526
+0.99534,0.99548,0.99562,0.99571,0.99587,0.99600,0.99609,0.99623,0.99631,0.99645,0.99658,0.99673,0.99685,0.99699,0.99713,0.99725,0.99734,0.99742,0.99755,0.99764
+0.99777,0.99788,0.99795,0.99809,0.99816,0.99824,0.99831,0.99838,0.99851,0.99863,0.99875,0.99884,0.99892,0.99899,0.99906,0.99916,0.99926,0.99935,0.99946,0.99955
+0.99963,0.99969,0.99977,0.99981,0.99985,0.99989,0.99991,0.99994,0.99996,0.99997,0.99997,0.99999,1.00000,1.00000,1.00000,1.00000
+pwc1u2fract=0.00068,0.00136,0.00226,0.00350,0.00477,0.00642,0.00839,0.01054,0.01310,0.01612,0.01975,0.02384,0.02847,0.03353,0.03910,0.04509,0.05171,0.05884,0.06635
+0.07454,0.08304,0.09180,0.10107,0.11023,0.11972,0.12937,0.13946,0.14970,0.16003,0.17062,0.18118,0.19182,0.20260,0.21321,0.22455,0.23563,0.24657,0.25737,0.26831
+0.27909,0.28991,0.30125,0.31218,0.32295,0.33392,0.34486,0.35601,0.36703,0.37743,0.38837,0.39938,0.41008,0.42067,0.43161,0.44207,0.45297,0.46369,0.47407,0.48469
+0.49494,0.50585,0.51662,0.52726,0.53760,0.54801,0.55788,0.56808,0.57818,0.58823,0.59841,0.60830,0.61809,0.62821,0.63811,0.64746,0.65734,0.66745,0.67686,0.68602
+0.69572,0.70495,0.71402,0.72343,0.73245,0.74159,0.75043,0.75934,0.76822,0.77697,0.78550,0.79378,0.80228,0.81013,0.81794,0.82572,0.83366,0.84100,0.84827,0.85535
+0.86214,0.86869,0.87497,0.88095,0.88692,0.89240,0.89757,0.90274,0.90765,0.91181,0.91622,0.92016,0.92391,0.92746,0.93069,0.93356,0.93639,0.93899,0.94140,0.94361
+0.94580,0.94782,0.94983,0.95156,0.95318,0.95479,0.95637,0.95793,0.95934,0.96063,0.96190,0.96313,0.96440,0.96550,0.96660,0.96772,0.96861,0.96968,0.97064,0.97158
+0.97248,0.97324,0.97397,0.97468,0.97544,0.97618,0.97682,0.97752,0.97823,0.97886,0.97946,0.97998,0.98067,0.98122,0.98170,0.98220,0.98268,0.98315,0.98354,0.98391
+0.98435,0.98477,0.98510,0.98547,0.98583,0.98618,0.98653,0.98692,0.98728,0.98760,0.98788,0.98815,0.98837,0.98861,0.98886,0.98901,0.98925,0.98948,0.98968,0.98989
+0.99007,0.99024,0.99040,0.99055,0.99074,0.99090,0.99106,0.99120,0.99134,0.99147,0.99160,0.99172,0.99188,0.99203,0.99216,0.99228,0.99247,0.99257,0.99269,0.99282
+0.99294,0.99309,0.99321,0.99333,0.99346,0.99357,0.99370,0.99382,0.99394,0.99404,0.99415,0.99427,0.99439,0.99450,0.99461,0.99474,0.99487,0.99497,0.99506,0.99516
+0.99526,0.99539,0.99552,0.99562,0.99572,0.99585,0.99597,0.99610,0.99623,0.99636,0.99649,0.99662,0.99673,0.99689,0.99701,0.99712,0.99725,0.99734,0.99745,0.99756
+0.99768,0.99782,0.99792,0.99802,0.99820,0.99830,0.99839,0.99853,0.99864,0.99874,0.99885,0.99894,0.99905,0.99918,0.99924,0.99932,0.99941,0.99950,0.99958,0.99961
+0.99965,0.99970,0.99977,0.99981,0.99985,0.99988,0.99991,0.99992,0.99993,0.99995,0.99996,0.99997,0.99998,0.99999,1.00000,1.00000
+pwc1x1fract=0.00072,0.00158,0.00255,0.00370,0.00517,0.00682,0.00873,0.01113,0.01404,0.01756,0.02140,0.02590,0.03108,0.03695,0.04320,0.04992,0.05750,0.06531,0.07347
+0.08168,0.09051,0.09950,0.10918,0.11895,0.12915,0.13890,0.14909,0.15956,0.17033,0.18116,0.19184,0.20296,0.21399,0.22504,0.23610,0.24674,0.25752,0.26827,0.27930
+0.29007,0.30127,0.31272,0.32337,0.33449,0.34518,0.35576,0.36669,0.37776,0.38847,0.39912,0.40988,0.42053,0.43123,0.44187,0.45239,0.46313,0.47362,0.48433,0.49498
+0.50531,0.51568,0.52568,0.53596,0.54634,0.55634,0.56654,0.57683,0.58657,0.59670,0.60670,0.61662,0.62671,0.63656,0.64634,0.65594,0.66527,0.67445,0.68398,0.69334
+0.70270,0.71176,0.72090,0.72999,0.73862,0.74728,0.75604,0.76492,0.77318,0.78141,0.78965,0.79800,0.80616,0.81397,0.82190,0.82957,0.83699,0.84420,0.85158,0.85845
+0.86521,0.87161,0.87777,0.88365,0.88909,0.89446,0.89953,0.90444,0.90884,0.91328,0.91743,0.92150,0.92478,0.92792,0.93080,0.93359,0.93626,0.93886,0.94103,0.94324
+0.94535,0.94736,0.94929,0.95107,0.95265,0.95408,0.95574,0.95715,0.95853,0.95989,0.96106,0.96223,0.96337,0.96475,0.96587,0.96702,0.96796,0.96900,0.97009,0.97103
+0.97199,0.97284,0.97366,0.97458,0.97544,0.97626,0.97706,0.97766,0.97836,0.97902,0.97966,0.98020,0.98081,0.98134,0.98188,0.98237,0.98293,0.98339,0.98381,0.98427
+0.98466,0.98516,0.98560,0.98598,0.98636,0.98665,0.98694,0.98727,0.98752,0.98781,0.98807,0.98836,0.98862,0.98890,0.98913,0.98938,0.98955,0.98975,0.98991,0.99013
+0.99035,0.99055,0.99068,0.99087,0.99101,0.99114,0.99128,0.99143,0.99156,0.99169,0.99182,0.99197,0.99211,0.99222,0.99237,0.99255,0.99266,0.99277,0.99291,0.99302
+0.99312,0.99324,0.99336,0.99347,0.99362,0.99375,0.99389,0.99399,0.99410,0.99426,0.99434,0.99446,0.99460,0.99475,0.99491,0.99499,0.99511,0.99523,0.99531,0.99543
+0.99551,0.99560,0.99570,0.99579,0.99589,0.99602,0.99610,0.99623,0.99639,0.99650,0.99660,0.99670,0.99683,0.99695,0.99705,0.99717,0.99727,0.99736,0.99746,0.99756
+0.99769,0.99778,0.99789,0.99801,0.99814,0.99828,0.99840,0.99855,0.99864,0.99872,0.99881,0.99889,0.99898,0.99909,0.99919,0.99932,0.99938,0.99945,0.99952,0.99961
+0.99968,0.99971,0.99975,0.99981,0.99986,0.99990,0.99992,0.99993,0.99995,0.99997,0.99997,0.99999,0.99999,0.99999,1.00000,1.00000
+pwc1x2fract=0.00070,0.00151,0.00248,0.00362,0.00514,0.00692,0.00890,0.01145,0.01445,0.01785,0.02184,0.02645,0.03162,0.03725,0.04345,0.05044,0.05776,0.06563,0.07410
+0.08292,0.09191,0.10119,0.11123,0.12114,0.13116,0.14131,0.15185,0.16227,0.17293,0.18372,0.19462,0.20535,0.21647,0.22724,0.23789,0.24830,0.25937,0.27022,0.28145
+0.29243,0.30284,0.31323,0.32376,0.33431,0.34497,0.35561,0.36641,0.37703,0.38789,0.39858,0.40932,0.41980,0.43052,0.44130,0.45195,0.46251,0.47277,0.48341,0.49406
+0.50394,0.51429,0.52445,0.53491,0.54505,0.55523,0.56515,0.57472,0.58452,0.59434,0.60433,0.61393,0.62353,0.63324,0.64280,0.65182,0.66133,0.67138,0.68086,0.69021
+0.69929,0.70906,0.71802,0.72661,0.73550,0.74446,0.75305,0.76161,0.77025,0.77850,0.78643,0.79442,0.80243,0.81038,0.81853,0.82634,0.83411,0.84153,0.84876,0.85576
+0.86267,0.86929,0.87577,0.88174,0.88755,0.89287,0.89778,0.90277,0.90738,0.91151,0.91574,0.91942,0.92292,0.92643,0.92947,0.93242,0.93525,0.93780,0.94002,0.94223
+0.94434,0.94632,0.94825,0.94994,0.95170,0.95328,0.95480,0.95621,0.95750,0.95885,0.96016,0.96136,0.96250,0.96362,0.96468,0.96581,0.96674,0.96775,0.96877,0.96978
+0.97075,0.97158,0.97239,0.97325,0.97403,0.97475,0.97544,0.97628,0.97691,0.97769,0.97832,0.97887,0.97948,0.98004,0.98060,0.98113,0.98162,0.98207,0.98253,0.98298
+0.98344,0.98392,0.98435,0.98477,0.98509,0.98549,0.98585,0.98619,0.98652,0.98688,0.98710,0.98740,0.98773,0.98799,0.98821,0.98846,0.98872,0.98896,0.98917,0.98939
+0.98962,0.98981,0.98998,0.99012,0.99025,0.99041,0.99063,0.99076,0.99091,0.99106,0.99121,0.99135,0.99145,0.99155,0.99168,0.99179,0.99192,0.99202,0.99217,0.99231
+0.99246,0.99260,0.99276,0.99285,0.99297,0.99310,0.99327,0.99339,0.99354,0.99368,0.99379,0.99392,0.99406,0.99415,0.99429,0.99438,0.99451,0.99464,0.99476,0.99489
+0.99504,0.99518,0.99530,0.99539,0.99549,0.99555,0.99568,0.99582,0.99595,0.99607,0.99621,0.99630,0.99646,0.99661,0.99668,0.99677,0.99689,0.99699,0.99710,0.99721
+0.99732,0.99745,0.99757,0.99771,0.99783,0.99792,0.99803,0.99816,0.99826,0.99836,0.99850,0.99860,0.99872,0.99884,0.99895,0.99905,0.99918,0.99931,0.99940,0.99948
+0.99960,0.99966,0.99972,0.99975,0.99981,0.99985,0.99988,0.99990,0.99993,0.99995,0.99996,0.99997,0.99998,0.99998,0.99999,1.00000
+pwc1v1fract=0.00078,0.00170,0.00279,0.00390,0.00525,0.00694,0.00898,0.01145,0.01418,0.01763,0.02136,0.02588,0.03057,0.03586,0.04160,0.04823,0.05486,0.06207,0.07016
+0.07844,0.08715,0.09596,0.10504,0.11487,0.12474,0.13523,0.14555,0.15632,0.16684,0.17763,0.18841,0.19932,0.21015,0.22112,0.23221,0.24340,0.25460,0.26520,0.27641
+0.28752,0.29857,0.30997,0.32122,0.33187,0.34278,0.35381,0.36435,0.37584,0.38652,0.39777,0.40858,0.41917,0.43019,0.44078,0.45176,0.46231,0.47314,0.48409,0.49439
+0.50481,0.51544,0.52619,0.53648,0.54684,0.55767,0.56804,0.57847,0.58859,0.59909,0.60917,0.61880,0.62842,0.63801,0.64773,0.65730,0.66675,0.67633,0.68559,0.69493
+0.70421,0.71321,0.72263,0.73181,0.74091,0.74986,0.75887,0.76743,0.77643,0.78510,0.79336,0.80139,0.80918,0.81723,0.82517,0.83336,0.84077,0.84803,0.85513,0.86171
+0.86834,0.87486,0.88095,0.88688,0.89250,0.89786,0.90296,0.90773,0.91203,0.91600,0.91970,0.92357,0.92702,0.93044,0.93337,0.93616,0.93874,0.94111,0.94341,0.94547
+0.94738,0.94942,0.95121,0.95292,0.95456,0.95596,0.95742,0.95891,0.96036,0.96160,0.96288,0.96412,0.96532,0.96642,0.96758,0.96864,0.96971,0.97065,0.97166,0.97250
+0.97341,0.97429,0.97519,0.97593,0.97659,0.97731,0.97801,0.97867,0.97929,0.97991,0.98059,0.98115,0.98173,0.98224,0.98277,0.98336,0.98380,0.98423,0.98473,0.98510
+0.98554,0.98591,0.98623,0.98660,0.98693,0.98731,0.98761,0.98785,0.98809,0.98837,0.98863,0.98890,0.98911,0.98931,0.98948,0.98967,0.98985,0.99009,0.99024,0.99045
+0.99068,0.99087,0.99105,0.99122,0.99137,0.99152,0.99166,0.99186,0.99196,0.99211,0.99224,0.99239,0.99253,0.99268,0.99282,0.99295,0.99308,0.99321,0.99328,0.99341
+0.99352,0.99363,0.99374,0.99385,0.99396,0.99410,0.99421,0.99434,0.99451,0.99464,0.99477,0.99487,0.99498,0.99509,0.99520,0.99529,0.99539,0.99550,0.99559,0.99568
+0.99579,0.99590,0.99599,0.99609,0.99620,0.99635,0.99650,0.99659,0.99666,0.99676,0.99688,0.99698,0.99708,0.99716,0.99725,0.99739,0.99750,0.99761,0.99771,0.99782
+0.99790,0.99800,0.99811,0.99824,0.99831,0.99846,0.99857,0.99868,0.99874,0.99882,0.99888,0.99897,0.99908,0.99918,0.99926,0.99938,0.99946,0.99954,0.99962,0.99967
+0.99973,0.99977,0.99980,0.99987,0.99991,0.99995,0.99996,0.99997,0.99997,0.99999,0.99999,0.99999,1.00000,1.00000,1.00000,1.00000
+pwc1v2fract=0.00070,0.00155,0.00256,0.00376,0.00522,0.00699,0.00898,0.01110,0.01362,0.01664,0.02026,0.02421,0.02888,0.03393,0.03963,0.04546,0.05221,0.05933,0.06695
+0.07522,0.08384,0.09260,0.10220,0.11175,0.12201,0.13182,0.14191,0.15249,0.16325,0.17338,0.18411,0.19487,0.20560,0.21669,0.22785,0.23869,0.24953,0.26077,0.27138
+0.28240,0.29370,0.30431,0.31524,0.32628,0.33692,0.34793,0.35888,0.36964,0.38041,0.39118,0.40184,0.41267,0.42342,0.43406,0.44474,0.45526,0.46592,0.47614,0.48672
+0.49734,0.50784,0.51843,0.52869,0.53936,0.54973,0.56006,0.56999,0.57990,0.58992,0.60014,0.61048,0.62052,0.63049,0.64045,0.65037,0.66032,0.67005,0.67990,0.68960
+0.69880,0.70776,0.71697,0.72607,0.73519,0.74418,0.75293,0.76172,0.77033,0.77895,0.78755,0.79574,0.80386,0.81209,0.82016,0.82788,0.83532,0.84237,0.84966,0.85650
+0.86312,0.86967,0.87621,0.88234,0.88798,0.89344,0.89862,0.90371,0.90830,0.91275,0.91687,0.92093,0.92499,0.92844,0.93173,0.93461,0.93724,0.93981,0.94246,0.94476
+0.94674,0.94890,0.95089,0.95266,0.95426,0.95586,0.95723,0.95870,0.95989,0.96131,0.96244,0.96352,0.96453,0.96567,0.96670,0.96775,0.96872,0.96963,0.97050,0.97144
+0.97227,0.97311,0.97386,0.97460,0.97533,0.97603,0.97672,0.97733,0.97799,0.97863,0.97923,0.97979,0.98034,0.98090,0.98142,0.98189,0.98230,0.98287,0.98329,0.98374
+0.98415,0.98454,0.98492,0.98522,0.98563,0.98601,0.98631,0.98656,0.98685,0.98709,0.98741,0.98769,0.98791,0.98821,0.98842,0.98866,0.98884,0.98907,0.98928,0.98948
+0.98963,0.98980,0.98999,0.99017,0.99033,0.99049,0.99068,0.99087,0.99106,0.99120,0.99133,0.99151,0.99160,0.99175,0.99190,0.99204,0.99212,0.99227,0.99238,0.99252
+0.99264,0.99279,0.99290,0.99303,0.99317,0.99329,0.99340,0.99354,0.99365,0.99379,0.99394,0.99406,0.99417,0.99426,0.99435,0.99447,0.99462,0.99475,0.99489,0.99501
+0.99512,0.99524,0.99535,0.99549,0.99561,0.99571,0.99580,0.99591,0.99602,0.99616,0.99627,0.99641,0.99653,0.99665,0.99677,0.99688,0.99699,0.99712,0.99722,0.99737
+0.99747,0.99760,0.99772,0.99784,0.99797,0.99807,0.99817,0.99829,0.99840,0.99854,0.99861,0.99872,0.99881,0.99891,0.99902,0.99912,0.99924,0.99931,0.99939,0.99945
+0.99954,0.99956,0.99963,0.99970,0.99974,0.99979,0.99984,0.99988,0.99993,0.99994,0.99997,0.99998,0.99999,0.99999,1.00000,1.00000
+pwc2v2fract=0.00082,0.00174,0.00269,0.00392,0.00528,0.00701,0.00888,0.01119,0.01378,0.01679,0.02050,0.02445,0.02903,0.03449,0.04038,0.04653,0.05333,0.06071,0.06847
+0.07645,0.08529,0.09421,0.10367,0.11335,0.12307,0.13334,0.14396,0.15448,0.16523,0.17613,0.18698,0.19775,0.20879,0.21995,0.23110,0.24243,0.25345,0.26421,0.27543
+0.28644,0.29730,0.30810,0.31887,0.32965,0.34029,0.35136,0.36271,0.37349,0.38404,0.39503,0.40582,0.41655,0.42705,0.43749,0.44764,0.45807,0.46848,0.47866,0.48896
+0.49922,0.50938,0.51943,0.52933,0.53909,0.54917,0.55906,0.56895,0.57891,0.58828,0.59798,0.60746,0.61683,0.62618,0.63550,0.64478,0.65348,0.66253,0.67133,0.68044
+0.68917,0.69778,0.70647,0.71512,0.72364,0.73205,0.74050,0.74843,0.75642,0.76412,0.77170,0.77944,0.78688,0.79393,0.80097,0.80814,0.81523,0.82222,0.82909,0.83582
+0.84238,0.84849,0.85441,0.86044,0.86588,0.87128,0.87652,0.88142,0.88615,0.89071,0.89519,0.89932,0.90322,0.90699,0.91078,0.91390,0.91704,0.91998,0.92277,0.92550
+0.92794,0.93031,0.93259,0.93484,0.93693,0.93911,0.94102,0.94287,0.94474,0.94656,0.94818,0.94972,0.95115,0.95261,0.95406,0.95553,0.95687,0.95816,0.95935,0.96060
+0.96175,0.96295,0.96414,0.96520,0.96627,0.96732,0.96828,0.96916,0.96998,0.97077,0.97160,0.97237,0.97324,0.97396,0.97460,0.97524,0.97589,0.97653,0.97718,0.97788
+0.97847,0.97906,0.97969,0.98023,0.98072,0.98128,0.98178,0.98234,0.98287,0.98334,0.98374,0.98416,0.98461,0.98498,0.98535,0.98571,0.98605,0.98637,0.98667,0.98701
+0.98734,0.98759,0.98787,0.98816,0.98838,0.98867,0.98892,0.98913,0.98939,0.98962,0.98983,0.99004,0.99025,0.99044,0.99062,0.99080,0.99097,0.99114,0.99132,0.99146
+0.99160,0.99174,0.99190,0.99209,0.99225,0.99240,0.99256,0.99273,0.99290,0.99306,0.99326,0.99342,0.99358,0.99377,0.99393,0.99403,0.99414,0.99431,0.99448,0.99464
+0.99481,0.99492,0.99504,0.99519,0.99534,0.99543,0.99553,0.99566,0.99582,0.99597,0.99610,0.99624,0.99638,0.99654,0.99667,0.99678,0.99692,0.99704,0.99715,0.99728
+0.99746,0.99756,0.99771,0.99784,0.99800,0.99815,0.99826,0.99840,0.99852,0.99866,0.99877,0.99888,0.99897,0.99907,0.99917,0.99929,0.99937,0.99946,0.99955,0.99961
+0.99965,0.99972,0.99976,0.99982,0.99987,0.99990,0.99995,0.99996,0.99997,0.99997,0.99999,0.99999,0.99999,0.99999,1.00000,1.00000
+pwc2v1fract=0.00078,0.00160,0.00261,0.00378,0.00511,0.00680,0.00864,0.01090,0.01336,0.01639,0.02006,0.02410,0.02872,0.03365,0.03944,0.04579,0.05256,0.05948,0.06698
+0.07536,0.08395,0.09273,0.10176,0.11106,0.12076,0.13116,0.14130,0.15150,0.16216,0.17310,0.18413,0.19527,0.20642,0.21720,0.22836,0.23946,0.25092,0.26212,0.27335
+0.28437,0.29514,0.30620,0.31713,0.32825,0.33901,0.35003,0.36100,0.37185,0.38282,0.39352,0.40394,0.41444,0.42512,0.43544,0.44593,0.45635,0.46695,0.47737,0.48812
+0.49865,0.50888,0.51925,0.52954,0.53962,0.54983,0.55981,0.56969,0.57951,0.58905,0.59876,0.60850,0.61814,0.62788,0.63708,0.64659,0.65574,0.66480,0.67371,0.68263
+0.69166,0.70054,0.70930,0.71761,0.72613,0.73399,0.74251,0.75051,0.75858,0.76644,0.77439,0.78217,0.78994,0.79728,0.80469,0.81195,0.81881,0.82568,0.83239,0.83902
+0.84532,0.85163,0.85759,0.86321,0.86872,0.87406,0.87937,0.88446,0.88931,0.89384,0.89809,0.90226,0.90609,0.90976,0.91321,0.91646,0.91961,0.92270,0.92562,0.92836
+0.93088,0.93318,0.93531,0.93741,0.93954,0.94163,0.94350,0.94528,0.94701,0.94875,0.95045,0.95191,0.95349,0.95503,0.95638,0.95769,0.95900,0.96021,0.96142,0.96255
+0.96371,0.96475,0.96572,0.96673,0.96766,0.96866,0.96963,0.97044,0.97134,0.97210,0.97292,0.97362,0.97434,0.97504,0.97563,0.97632,0.97694,0.97768,0.97834,0.97888
+0.97949,0.98015,0.98067,0.98117,0.98166,0.98214,0.98268,0.98318,0.98366,0.98406,0.98447,0.98486,0.98522,0.98554,0.98588,0.98621,0.98657,0.98687,0.98722,0.98753
+0.98784,0.98812,0.98840,0.98866,0.98894,0.98919,0.98938,0.98957,0.98977,0.99000,0.99019,0.99037,0.99056,0.99072,0.99090,0.99108,0.99121,0.99139,0.99157,0.99174
+0.99188,0.99201,0.99217,0.99235,0.99251,0.99268,0.99279,0.99292,0.99309,0.99322,0.99337,0.99352,0.99367,0.99380,0.99391,0.99405,0.99419,0.99436,0.99451,0.99464
+0.99479,0.99494,0.99510,0.99519,0.99533,0.99547,0.99562,0.99574,0.99586,0.99602,0.99618,0.99632,0.99644,0.99653,0.99664,0.99677,0.99687,0.99702,0.99713,0.99724
+0.99738,0.99753,0.99764,0.99777,0.99789,0.99804,0.99816,0.99827,0.99841,0.99850,0.99866,0.99874,0.99889,0.99901,0.99914,0.99921,0.99931,0.99945,0.99956,0.99962
+0.99966,0.99972,0.99975,0.99981,0.99985,0.99989,0.99992,0.99993,0.99995,0.99997,0.99998,0.99999,0.99999,0.99999,1.00000,1.00000
+pwc2x2fract=0.00073,0.00160,0.00276,0.00399,0.00543,0.00731,0.00928,0.01161,0.01438,0.01770,0.02163,0.02590,0.03103,0.03642,0.04256,0.04921,0.05631,0.06391,0.07229
+0.08113,0.08993,0.09957,0.10952,0.11921,0.12912,0.13946,0.14992,0.16057,0.17114,0.18184,0.19248,0.20322,0.21418,0.22502,0.23577,0.24689,0.25758,0.26841,0.27917
+0.28978,0.30064,0.31108,0.32152,0.33284,0.34378,0.35463,0.36503,0.37578,0.38635,0.39665,0.40689,0.41771,0.42810,0.43862,0.44904,0.45973,0.46976,0.48002,0.49016
+0.50058,0.51036,0.52037,0.53068,0.54099,0.55100,0.56110,0.57093,0.58066,0.59034,0.59995,0.60888,0.61856,0.62801,0.63752,0.64675,0.65594,0.66509,0.67446,0.68322
+0.69223,0.70064,0.70937,0.71797,0.72664,0.73527,0.74339,0.75190,0.76004,0.76826,0.77631,0.78431,0.79204,0.79982,0.80700,0.81426,0.82124,0.82826,0.83520,0.84194
+0.84826,0.85461,0.86082,0.86695,0.87263,0.87819,0.88334,0.88840,0.89298,0.89737,0.90160,0.90560,0.90932,0.91291,0.91625,0.91935,0.92223,0.92484,0.92752,0.92999
+0.93244,0.93450,0.93659,0.93862,0.94046,0.94229,0.94407,0.94556,0.94721,0.94890,0.95044,0.95185,0.95332,0.95464,0.95597,0.95720,0.95848,0.95962,0.96082,0.96191
+0.96301,0.96406,0.96506,0.96606,0.96708,0.96806,0.96897,0.96984,0.97074,0.97154,0.97244,0.97319,0.97404,0.97484,0.97550,0.97619,0.97690,0.97764,0.97831,0.97896
+0.97953,0.98010,0.98062,0.98118,0.98171,0.98224,0.98266,0.98315,0.98359,0.98405,0.98447,0.98491,0.98528,0.98563,0.98602,0.98640,0.98674,0.98709,0.98746,0.98786
+0.98813,0.98837,0.98869,0.98892,0.98912,0.98938,0.98960,0.98978,0.98996,0.99015,0.99031,0.99054,0.99078,0.99102,0.99121,0.99137,0.99157,0.99177,0.99194,0.99208
+0.99227,0.99240,0.99251,0.99265,0.99282,0.99301,0.99320,0.99335,0.99349,0.99361,0.99374,0.99386,0.99397,0.99411,0.99423,0.99436,0.99453,0.99468,0.99478,0.99488
+0.99500,0.99516,0.99527,0.99537,0.99547,0.99560,0.99576,0.99587,0.99602,0.99615,0.99623,0.99639,0.99648,0.99661,0.99676,0.99688,0.99699,0.99709,0.99723,0.99733
+0.99743,0.99756,0.99772,0.99785,0.99799,0.99814,0.99822,0.99833,0.99839,0.99850,0.99859,0.99870,0.99881,0.99892,0.99903,0.99911,0.99922,0.99932,0.99943,0.99953
+0.99959,0.99967,0.99971,0.99977,0.99984,0.99987,0.99989,0.99992,0.99993,0.99996,0.99998,0.99998,0.99999,1.00000,1.00000,1.00000
+pwc2x1fract=0.00078,0.00154,0.00253,0.00381,0.00508,0.00680,0.00868,0.01088,0.01359,0.01671,0.02053,0.02481,0.02959,0.03502,0.04090,0.04734,0.05437,0.06171,0.06951
+0.07760,0.08619,0.09528,0.10438,0.11385,0.12373,0.13397,0.14375,0.15435,0.16457,0.17535,0.18615,0.19674,0.20764,0.21809,0.22889,0.23955,0.25021,0.26135,0.27215
+0.28297,0.29361,0.30419,0.31499,0.32592,0.33645,0.34737,0.35828,0.36898,0.37970,0.38999,0.40044,0.41119,0.42210,0.43259,0.44295,0.45330,0.46406,0.47454,0.48481
+0.49493,0.50524,0.51578,0.52626,0.53667,0.54677,0.55689,0.56701,0.57708,0.58722,0.59682,0.60656,0.61613,0.62549,0.63513,0.64441,0.65360,0.66286,0.67203,0.68109
+0.69001,0.69875,0.70769,0.71670,0.72521,0.73369,0.74193,0.75051,0.75833,0.76652,0.77472,0.78281,0.79076,0.79848,0.80599,0.81343,0.82089,0.82799,0.83461,0.84166
+0.84860,0.85490,0.86109,0.86725,0.87296,0.87829,0.88345,0.88876,0.89359,0.89830,0.90259,0.90666,0.91062,0.91426,0.91743,0.92078,0.92371,0.92663,0.92909,0.93164
+0.93393,0.93612,0.93826,0.94014,0.94227,0.94413,0.94584,0.94759,0.94923,0.95068,0.95233,0.95355,0.95495,0.95630,0.95753,0.95890,0.95998,0.96125,0.96231,0.96341
+0.96454,0.96562,0.96665,0.96750,0.96844,0.96933,0.97031,0.97121,0.97202,0.97285,0.97363,0.97441,0.97513,0.97587,0.97659,0.97726,0.97792,0.97858,0.97922,0.97979
+0.98046,0.98108,0.98154,0.98207,0.98256,0.98301,0.98341,0.98389,0.98435,0.98471,0.98507,0.98551,0.98592,0.98627,0.98659,0.98692,0.98723,0.98752,0.98785,0.98813
+0.98838,0.98866,0.98889,0.98915,0.98942,0.98966,0.98986,0.99008,0.99028,0.99051,0.99067,0.99086,0.99106,0.99124,0.99137,0.99154,0.99171,0.99189,0.99206,0.99226
+0.99237,0.99251,0.99265,0.99282,0.99299,0.99311,0.99327,0.99341,0.99355,0.99367,0.99381,0.99397,0.99410,0.99420,0.99434,0.99447,0.99463,0.99475,0.99487,0.99498
+0.99509,0.99521,0.99531,0.99542,0.99551,0.99560,0.99573,0.99586,0.99599,0.99610,0.99620,0.99631,0.99643,0.99657,0.99671,0.99682,0.99691,0.99706,0.99717,0.99729
+0.99741,0.99756,0.99767,0.99785,0.99798,0.99809,0.99817,0.99828,0.99842,0.99851,0.99862,0.99873,0.99882,0.99891,0.99899,0.99909,0.99918,0.99927,0.99941,0.99947
+0.99954,0.99959,0.99967,0.99972,0.99978,0.99984,0.99986,0.99989,0.99992,0.99994,0.99995,0.99997,0.99999,1.00000,1.00000,1.00000
+pwc2u2fract=0.00081,0.00173,0.00285,0.00403,0.00559,0.00719,0.00925,0.01152,0.01402,0.01707,0.02074,0.02480,0.02956,0.03478,0.04054,0.04684,0.05369,0.06084,0.06840
+0.07660,0.08505,0.09368,0.10275,0.11227,0.12202,0.13195,0.14195,0.15210,0.16261,0.17326,0.18368,0.19417,0.20476,0.21521,0.22574,0.23676,0.24746,0.25840,0.26938
+0.28026,0.29132,0.30220,0.31295,0.32355,0.33436,0.34526,0.35608,0.36677,0.37758,0.38837,0.39912,0.40966,0.42022,0.43092,0.44122,0.45153,0.46195,0.47221,0.48228
+0.49219,0.50261,0.51293,0.52316,0.53349,0.54360,0.55349,0.56369,0.57363,0.58363,0.59337,0.60344,0.61309,0.62267,0.63227,0.64179,0.65135,0.66074,0.66999,0.67927
+0.68838,0.69782,0.70692,0.71576,0.72436,0.73271,0.74103,0.74942,0.75775,0.76615,0.77412,0.78221,0.79004,0.79797,0.80566,0.81335,0.82085,0.82821,0.83546,0.84216
+0.84888,0.85495,0.86115,0.86719,0.87305,0.87849,0.88401,0.88925,0.89379,0.89829,0.90272,0.90690,0.91060,0.91408,0.91757,0.92082,0.92395,0.92683,0.92963,0.93214
+0.93457,0.93706,0.93915,0.94112,0.94309,0.94479,0.94672,0.94838,0.95001,0.95156,0.95312,0.95467,0.95604,0.95747,0.95884,0.95998,0.96119,0.96242,0.96357,0.96476
+0.96583,0.96690,0.96786,0.96885,0.96975,0.97061,0.97151,0.97242,0.97327,0.97406,0.97484,0.97555,0.97623,0.97699,0.97770,0.97833,0.97901,0.97959,0.98018,0.98078
+0.98122,0.98173,0.98224,0.98270,0.98313,0.98356,0.98401,0.98442,0.98484,0.98523,0.98560,0.98595,0.98629,0.98661,0.98690,0.98718,0.98748,0.98776,0.98810,0.98838
+0.98861,0.98883,0.98905,0.98923,0.98948,0.98967,0.98986,0.99009,0.99027,0.99045,0.99061,0.99080,0.99099,0.99112,0.99127,0.99146,0.99160,0.99173,0.99191,0.99205
+0.99221,0.99235,0.99250,0.99263,0.99279,0.99292,0.99309,0.99322,0.99334,0.99343,0.99353,0.99364,0.99378,0.99390,0.99402,0.99412,0.99423,0.99433,0.99445,0.99462
+0.99475,0.99487,0.99503,0.99518,0.99532,0.99541,0.99555,0.99565,0.99579,0.99591,0.99604,0.99616,0.99629,0.99638,0.99648,0.99660,0.99675,0.99685,0.99700,0.99710
+0.99721,0.99733,0.99746,0.99759,0.99775,0.99787,0.99797,0.99813,0.99824,0.99836,0.99852,0.99861,0.99873,0.99885,0.99897,0.99907,0.99920,0.99929,0.99938,0.99947
+0.99955,0.99963,0.99968,0.99974,0.99981,0.99986,0.99991,0.99993,0.99996,0.99998,0.99999,1.00000,1.00000,1.00000,1.00000,1.00000
+pwc2u1fract=0.00081,0.00157,0.00263,0.00397,0.00541,0.00715,0.00915,0.01153,0.01420,0.01729,0.02116,0.02540,0.03006,0.03497,0.04060,0.04684,0.05335,0.06035,0.06794
+0.07602,0.08431,0.09310,0.10190,0.11112,0.12071,0.13056,0.14098,0.15123,0.16149,0.17189,0.18251,0.19373,0.20427,0.21514,0.22624,0.23697,0.24789,0.25881,0.26949
+0.28039,0.29113,0.30205,0.31294,0.32399,0.33500,0.34601,0.35688,0.36755,0.37822,0.38934,0.40004,0.41079,0.42145,0.43211,0.44278,0.45310,0.46373,0.47402,0.48449
+0.49478,0.50480,0.51569,0.52596,0.53618,0.54613,0.55622,0.56647,0.57645,0.58641,0.59618,0.60611,0.61604,0.62557,0.63498,0.64441,0.65391,0.66323,0.67272,0.68219
+0.69140,0.70057,0.70993,0.71897,0.72797,0.73669,0.74539,0.75406,0.76221,0.77045,0.77887,0.78692,0.79482,0.80251,0.81007,0.81769,0.82540,0.83264,0.83960,0.84638
+0.85317,0.85966,0.86606,0.87189,0.87735,0.88286,0.88818,0.89307,0.89794,0.90249,0.90672,0.91070,0.91453,0.91816,0.92165,0.92491,0.92773,0.93052,0.93324,0.93585
+0.93824,0.94055,0.94256,0.94455,0.94652,0.94844,0.95013,0.95185,0.95346,0.95502,0.95653,0.95806,0.95941,0.96071,0.96195,0.96324,0.96444,0.96558,0.96684,0.96797
+0.96909,0.97007,0.97100,0.97186,0.97275,0.97361,0.97444,0.97516,0.97596,0.97675,0.97747,0.97827,0.97893,0.97953,0.98012,0.98066,0.98122,0.98177,0.98232,0.98280
+0.98330,0.98369,0.98413,0.98451,0.98492,0.98530,0.98568,0.98600,0.98638,0.98669,0.98704,0.98736,0.98763,0.98786,0.98811,0.98835,0.98862,0.98892,0.98917,0.98937
+0.98962,0.98984,0.99005,0.99024,0.99046,0.99068,0.99086,0.99101,0.99120,0.99138,0.99152,0.99164,0.99179,0.99194,0.99206,0.99218,0.99229,0.99242,0.99255,0.99266
+0.99278,0.99292,0.99304,0.99316,0.99328,0.99341,0.99354,0.99368,0.99379,0.99393,0.99408,0.99422,0.99434,0.99448,0.99460,0.99475,0.99486,0.99496,0.99506,0.99519
+0.99530,0.99542,0.99552,0.99565,0.99574,0.99585,0.99595,0.99607,0.99624,0.99640,0.99654,0.99665,0.99675,0.99687,0.99698,0.99708,0.99721,0.99732,0.99748,0.99762
+0.99775,0.99790,0.99802,0.99814,0.99822,0.99839,0.99851,0.99859,0.99871,0.99883,0.99893,0.99903,0.99914,0.99921,0.99929,0.99939,0.99946,0.99952,0.99961,0.99966
+0.99973,0.99977,0.99982,0.99987,0.99990,0.99993,0.99996,0.99998,0.99998,0.99999,0.99999,1.00000,1.00000,1.00000,1.00000,1.00000
diff --git a/PARAM/SHMS/DC/CALIB/pdc_tzero_per_wire.param b/PARAM/SHMS/DC/CALIB/pdc_tzero_per_wire.param
index 9424471e..374fec84 100644
--- a/PARAM/SHMS/DC/CALIB/pdc_tzero_per_wire.param
+++ b/PARAM/SHMS/DC/CALIB/pdc_tzero_per_wire.param
@@ -1,88 +1,88 @@
ptzero1u1=
-0,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000, 0.000000, 0.000000, 0.000000, 0.000000
-0.000000, 0.000000, 1.795920, 5.815350, -3.864150, -5.927540, 10.719800, -4.813430, 0.404404, 0.755655, 0.000000, 1.565370, 4.950330, -0.345324, -5.542140, 2.861590
-1.534980, -15.681300, 1.979280, 4.854760, 6.446210, 1.620820, 0.177650, 0.680271, 1.611940, 5.395360, 2.935390, 1.323300, 5.197490, 4.067040, 0.772863, 1.704960
-2.751230, 6.300950, 0.858680, 4.487330, 3.288680, 2.332650, 4.070710, 1.355750, 3.087380, 2.726840, 3.151510, 2.523060, 1.251760, 4.044590, 1.242660, 1.512210
-3.171260, 3.743360, 4.626910, 5.624360, 6.038540, 7.064390, 3.798730, 4.800300, 5.590810, 8.440760, 7.197160, 3.502090, 5.330660, 9.031360, 6.141750, 5.212600
-10.965900, 2.301490, 8.010930, 7.743700, 3.590260, 6.029060, 0.816870, -5.213110, 9.752780, 8.307690, 9.460380, 8.567180, 0.000000, 5.231090, 11.354300, 7.508670
-6.368630, 9.673470, 8.767530, 13.277100, -16.695200, 7.085900, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000
+0,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000, 0.000000, 11.737706, 0.000000, -18.219366
+0.000000, -8.596565, -19.709377, -3.888760, -10.093304, 2.922193, -5.095818, -1.462486, 1.550307, 1.868138, -1.254447, 1.315511, -0.776645, -3.174524, -2.942045, -2.571896
+3.252089, -0.994770, 1.251117, 2.212229, 0.360121, 0.050024, 2.714125, 3.315600, 3.263798, 2.579078, 3.798322, 2.740594, 1.057505, 0.138616, 0.061684, 0.837149
+0.167397, -0.169214, -0.198894, 0.184975, -0.105104, -0.311205, 0.489515, -0.169327, 0.030759, -0.573357, -0.472877, 0.586034, 0.151413, 0.273465, -0.269756, -0.982273
+3.582255, 4.598095, 2.651357, 2.616282, 0.732567, 4.875063, 1.279132, 1.570648, 2.029647, 3.208437, 4.650519, 3.427307, 4.787766, 3.674315, 3.400645, -0.707859
+0.447735, -0.321946, 4.087402, 2.526979, 0.398617, 3.959752, 2.701463, 7.761576, 3.491549, 2.344939, 2.775912, 6.148535, 0.595573, 3.061118, 3.264426, 7.903000
+4.873721, 6.153708, -4.776695, -3.199948, -0.856307, 3.757691, -7.031933, 30.533007, -0.000000, 0.000000, 0.000000
ptzero1u2=
-0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000, 0.000000, 0.000000, 0.000000, 0.000000
-0.000000, 10.916700, -4.672130, 6.038300, 4.164760, 8.883000, 6.370700, 9.557240, -0.648253, 4.362870, -10.987600, 3.678900, 7.059010, -1.258060, 5.024760, 5.735800
-4.186680, -0.536697, 0.236342, 2.694360, 4.960730, -3.547860, 2.949800, 2.134110, 4.653480, 5.121060, 5.471440, 2.201530, 3.393430, 7.740120, 4.183370, 4.699110
-6.130200, 2.142410, 2.887710, 4.458210, 1.461180, 0.891988, 3.531650, 3.226970, 3.356920, 1.415850, 3.002660, 1.083940, -2.429840, 2.518900, 4.608790, 2.840310
-5.780790, 1.942860, 4.077160, 4.112970, 2.227040, 2.115580, 3.183260, 4.864210, 3.400000, 0.632734, 1.835060, 4.341910, 5.565070, 4.406240, 7.122810, 8.552850
-5.535050, -0.324697, 3.381260, 7.990940, 3.640930, 23.795600, 6.688140, 8.332550, 7.315810, 10.925900, 0.000000, 1.500810, -0.318182, 4.046240, 9.760080, 12.506000
-11.769200, 6.388520, 4.244090, 1.200000, 2.818180, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000
+0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,-10.186916, 22.542977, 2.081081, -11.136029, -25.537133
+-26.600985, 2.898769, -3.754928, -9.435411, 0.000000, 3.518588, -10.345972, 2.828941, 4.088474, -15.231425, -0.869406, 0.013588, -0.838881, 1.393585, 1.179151, -1.945006
+-0.139143, 0.029517, 1.305504, -0.664838, 0.518358, 2.494754, -0.154310, 0.142750, 0.156720, 0.470993, 0.412774, 0.762821, 3.946781, 4.561689, 5.452406, -0.094996
+0.560123, -0.145073, 0.459576, -0.362110, 0.187945, 0.303776, -0.251532, 0.085632, 0.531553, 0.348465, 0.056768, -0.343636, 0.549847, 0.772046, 4.013889, 0.155580
+0.739518, 4.818560, 3.255604, -0.387000, 1.671022, 0.166904, -0.049546, 3.778922, 1.645198, 3.233069, 0.265056, 4.128437, 5.161766, -0.901327, 4.184098, 3.812798
+5.358346, 4.701613, 5.127571, 1.854201, 4.549845, 5.291580, 1.796003, 5.949658, 3.491952, 3.176604, 3.906187, 0.829896, -0.302273, -4.112992, -11.380168, 0.638508
+-11.116805, -8.329832, -2.379425, -25.373400, 3.799371, -13.616718, -10.727546, 0.000000, 0.000000, 0.000000, 0.000000
ptzero1x1=
--1.862060,-5.060360,0.890689,-2.123470,0.266049,-4.394440,-3.181590,-3.011560,-7.340200,-1.131910,-5.684140,-3.339910, -3.603510, -0.171227, 2.670460, -3.477110
-0.696540, 1.928030, -0.766734, 0.935069, -4.320620, -2.370120, -4.209710, 1.379730, -0.911690, -2.316800, -5.893290, -3.126300, -7.080260, 1.282180, 3.116260, 4.467810
-2.244710, 5.122060, 6.816590, 3.063140, 1.574480, 0.466924, 2.411390, 0.143337, 0.770216, -1.611270, -2.379990, -0.308662, 1.031620, 0.181693, 3.392230, 1.548930
-3.553960, 7.047760, 4.506460, 8.317220, 2.723510, 4.768800, 5.221840, 3.216500, 0.182255, 2.475740, -1.318110, -1.395740, 3.349560, 6.663600, 7.749050, 4.382410
-9.047220, 10.447800, 11.142700, 10.615600, 10.124500, 12.012400, 10.867500, 14.273300, 12.222600, 12.880300, 10.129000, 8.611050, 8.293280, 7.103690, -2.305470
+-2.251630,-6.426284,-1.070542,-7.051324,-3.927424,-4.163845,-5.052772,-5.232631,-6.708763,-9.147509,-5.644984,-6.129807, -3.038017, -5.044662, -4.566734, -3.420908
+-0.707483, -4.275923, -4.442430, -3.533496, -5.881593, -3.349601, -5.607669, -4.633543, -4.342309, -4.007149, -3.374276, -4.883010, -3.499129, -4.376873, -0.006874, 7.573937
+4.709190, 1.990800, 0.340001, 3.584672, -0.430216, 0.159015, -4.231939, 0.793594, -0.910754, -5.616183, -4.063920, -4.499126, -1.315427, 0.936095, 0.478367, 0.208568
+4.386366, 0.882352, 0.993212, 0.752556, 0.871640, 0.566292, 0.402137, 0.169559, 0.545422, 0.233999, -0.008753, 0.181312, -0.449078, 0.551866, 4.550678, 1.625968
+7.531332, 9.483513, 9.026458, 5.717941, 10.668072, 10.131106, 10.080386, 9.597945, 11.280614, 10.218127, 7.679981, 4.473191, 4.413067, 6.908367, 4.685325
ptzero1x2=
--0.388060,10.210500,5.598810,6.123810,7.651880,2.782440,-8.042550,2.420080,0.388554,7.726910,3.171390,7.548380, 6.562330, 4.760400, 4.175150, 5.784480
-4.376930, 3.342330, 4.798840, -1.665290, 1.520700, 0.544997, 3.660080, 1.212270, 2.039610, -0.397125, 2.737430, 4.216970, 3.810030, 6.409610, 4.638020, -1.012270
-0.753259, 2.973150, 1.977040, 0.638232, -1.982390, -2.790790, -3.577970, -1.763710, -4.006910, -2.834720, -0.205417, -1.636630, 0.708044, 0.485227, -4.125460, -0.077492
-3.829980, -0.006536, 5.149800, 0.875090, 3.189890, 2.302570, 0.715657, 0.682134, -3.278750, -2.829270, -2.334770, -4.668020, 0.100722, 0.163470, -0.583767, 5.641670
-3.717870, 6.252000, 5.077890, 4.494680, 8.895850, 8.389630, 5.788920, 4.952720, 6.317280, 4.454020, 2.793250, 1.912540, 0.923993, 0.485133, 0.006984
+3.193140,7.469396,4.377968,7.147898,2.179453,3.918345,2.099492,0.258660,0.534490,-0.441888,-0.040630,0.754379, 1.452394, 3.168537, 1.827288, 1.243380
+-1.436232, 1.365202, 0.418930, 0.225693, -0.955069, 0.833650, -2.176226, 3.418739, 0.379277, -1.092438, -0.214427, 1.336594, 0.992301, 2.073226, -0.921482, -0.884764
+-3.464066, -4.211605, -0.738412, -7.099800, -5.031097, -4.678625, -5.641162, -8.050759, -6.261626, -5.781604, -4.880952, -4.959928, -4.310327, -4.948786, -3.650533, -2.824688
+0.453688, 0.868804, -0.169756, 1.031630, -0.067409, -0.556105, -4.134691, -3.501935, -4.817857, -4.149302, -5.041570, -4.466312, -4.418636, -3.175924, -0.044875, 0.941851
+1.196680, 3.608019, 0.755795, 4.792789, 4.420821, 4.738204, 4.839457, 4.053704, 5.024758, 1.087425, 0.216561, -1.528309, 0.553710, 0.635227, -0.159189
ptzero1v1=
-0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,16.000000,0.000000,0.000000, -3.248620, 10.544300, 5.723470, -0.167832
-5.791570, 0.000000, 0.000000, -16.601400, 7.211620, 6.408960, 3.399260, 1.468270, 5.394960, 0.526718, 5.528730, -1.830630, 4.259580, -1.791670, 7.464050, 2.892630
-1.419100, 3.926050, 4.166810, 5.864460, 0.595406, 1.371100, 5.573380, 3.296000, 7.125270, 6.036240, 5.291410, 6.153040, 5.315270, 8.171190, 7.264820, 8.961300
-9.878830, 8.970030, 9.517580, 7.025000, 8.266900, 9.315160, 7.651030, 8.559870, 11.303900, 7.863830, 11.701000, 9.156380, 9.357870, 7.235850, 10.026400, -2.578380
--0.357369, -2.929150, -1.901050, -2.151240, -2.265530, -7.198300, -2.665640, -0.212098, -3.156100, 0.106613, 1.210510, 1.478900, -5.898540, 4.406570, -4.721420, 4.283780
--0.070747, 0.353721, 3.512260, 4.885820, 4.807120, 5.441330, 3.471640, 3.662200, 4.331610, -7.154520, 2.599710, 1.593700, 3.987850, 5.787910, 19.733300, 14.638900
-0.000000, 3.464290, -19.447200, 5.882980, -17.459500, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000
+0.000000,0.000000,0.000000,0.000000,0.000000,-1.551351,-6.266703,0.000000,0.000000,-29.522777,-3.586469,1.250900, 0.074155, 1.977772, -0.007968, 1.475733
+-0.750025, -2.163524, -1.475238, -10.322831, -8.149647, -4.356046, 0.976002, -0.247428, -0.751697, -0.103731, 0.145272, 1.586233, 4.865697, 0.436067, -0.128442, -0.765906
+4.448741, 4.761456, 0.173060, 2.476602, 2.118164, 0.499495, 3.585117, -0.505931, 0.815364, 3.967533, 2.183604, 4.738548, 4.351805, 5.709207, 4.470957, 5.364703
+9.676337, 5.977037, 7.964222, 5.440858, 5.293707, 4.591992, 5.325566, 5.230815, 9.077105, 8.556767, 6.104932, 4.589074, 5.291665, 9.109364, 8.459031, -2.052650
+-3.574798, -4.585319, -4.303433, -4.722113, -4.753678, -2.899253, -3.527449, -4.776971, -0.664041, -1.056626, -1.696099, -0.989595, -1.222239, -0.202124, -0.368117, 0.224515
+3.110124, 0.311299, -1.000451, -1.643235, -1.098514, 0.531364, 1.414292, 0.446002, 2.653499, 1.034833, 3.085075, 2.304502, 5.752169, 3.682686, 5.274006, 7.015613
+2.893743, 7.744536, 5.026565, -5.671143, -24.353572, -0.247839, 0.000000, 16.622843, 0.000000, 0.000000, 0.000000
ptzero1v2=
-0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,24.857100,0.000000,0.000000, 25.060400, 4.133520, 8.519300, 5.222520
-6.879550, 8.261270, 7.433530, 3.431880, 7.536800, -7.283120, 1.879330, 2.621920, 2.068820, 8.948500, 4.047470, 5.543960, 6.173630, 8.318030, 7.155910, 3.704420
-9.330260, 7.558700, 7.531610, 0.084241, 5.609340, 5.630570, 7.967230, 5.188840, 10.209800, 9.106830, 6.858900, 8.087110, 1.440150, 1.414830, 2.399680, -0.808321
-1.099900, 1.120240, 0.300221, -2.169230, 0.451327, -2.473590, -2.036820, 0.887898, 0.369413, -1.016780, -0.898206, -1.164890, -0.681405, -0.286733, -0.384938, -3.059400
--5.289670, 0.488029, -0.363647, -0.928181, -1.304460, -2.630200, 0.118533, 0.318898, -2.213060, 0.853546, -2.025040, 2.121240, 1.674960, 1.830370, 2.004430, 4.969520
-2.716170, 4.110770, 7.479290, 5.380610, 5.901590, 3.801040, 2.820150, 5.148150, 7.095710, 7.548680, 4.854270, -21.165700, 0.000000, 7.825210, 6.105180, 11.006400
-0.380623, 13.549600, -6.625900, 7.572230, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000
+0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,-1.124471,-24.217226,5.717423, 4.510537, 4.150198, 3.705375, 6.283947
+3.613529, 3.282115, 1.538781, 5.011295, 2.809856, 2.887517, -2.034410, 3.366252, 2.603920, -7.384743, 4.586932, 3.992366, 8.671246, 6.883056, 4.651595, 2.561899
+4.507203, 5.264387, 3.703906, 5.192153, 1.943195, 5.468719, 3.584495, 4.824910, 5.841444, 4.662289, 5.829926, 8.063006, -1.152384, -1.192375, -1.392113, 0.266267
+-1.608814, -2.292075, -1.393527, -4.984799, -4.479954, -5.085875, -4.814115, -4.390844, -4.802972, -1.995489, -4.325520, -4.867185, -4.474102, -4.362590, -4.515954, -5.366059
+-5.041128, -4.804581, -4.011879, -5.147962, -4.740239, -3.361182, -4.346762, -0.514819, -1.307673, -3.742392, -0.439914, 0.182234, 1.754654, -0.755360, 1.100271, 4.249456
+0.547761, 3.353873, 1.670908, 0.222098, 4.065605, 4.060698, 3.365647, -0.242948, 2.634990, -1.429639, 4.854711, 3.363185, -1.589860, 2.760188, 2.164497, 6.613450
+3.434854, 0.430412, 5.977401, -2.414894, -12.915452, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000
ptzero2v2=
-0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000, 0.000000, -11.555600, 4.608350, 0.000000
--10.000000, 4.082760, 0.000000, 0.000000, -2.180610, -4.226240, 4.414830, -1.745280, 8.621950, 5.624680, -2.607140, -0.681955, 6.067340, 13.311900, 9.795680, 11.718100
-8.676360, 5.623300, 9.703400, 9.516520, 11.716700, 10.488800, 8.604710, 10.642600, 5.977520, 11.677700, 8.582270, 7.817510, 12.096200, 5.755470, 10.820500, 11.099900
-9.506290, 9.779980, 7.397800, 8.484320, 7.980850, 8.076540, 7.184750, 6.930630, 8.418590, 7.962490, 7.990790, 6.963160, -0.431749, 3.757360, 4.615140, 2.893580
-3.444480, -3.420450, -3.933450, -2.651720, -0.127254, -1.043080, -1.475670, -0.080271, 0.068876, 2.165830, 4.720000, 0.991301, 1.594030, 5.763560, 8.133160, 10.019800
-6.025830, -1.821310, 7.595050, 3.912210, 3.349220, 4.221970, 5.452880, 3.657390, 3.391830, 1.571480, 2.540660, 3.808400, 7.113780, 2.210830, 0.901170, 5.282630
-5.211540, 3.723870, 1.217340, 6.076170, 3.841830, 7.036590, 12.992400, 0.000000, 3.265860, 9.692310, 0.000000
+0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,-2.874082,0.000000,-39.750000,0.000000, -1.344906, -36.051142, 3.701519, -2.707838
+2.423867, -17.837071, -1.788165, -20.077432, -2.796743, -10.625253, -5.297905, 2.682962, -0.105740, 0.781828, -0.758277, 0.094761, 8.090251, 8.852005, 9.782540, 7.776518
+9.834901, 4.375272, 5.502201, 8.052743, 5.356021, 3.640592, 5.474690, 8.460758, 4.399497, 9.242440, 9.539740, 7.771080, 8.652890, 9.227889, 8.870403, 7.720761
+9.178195, 5.353182, 5.425670, 5.649495, 5.285828, 4.981052, 5.387508, 5.626353, 5.424275, 4.827419, 5.127994, 5.083649, -0.385558, -0.566800, 2.431357, -0.809385
+-0.690341, 0.170127, -0.714967, -0.282690, -1.156368, -4.994609, -0.197328, -0.340535, 0.016642, -3.380818, -0.243427, -5.655406, 4.576470, 0.007464, 4.262034, 4.260324
+0.499093, -0.498844, 0.932504, 4.366748, 3.621820, -1.217561, -0.271748, -0.252796, 0.005325, 3.017543, 0.477952, -1.350115, -0.126434, 4.694625, 3.892524, -1.086741
+2.774948, -1.734968, 3.088991, 1.739675, -1.162498, -1.200825, 9.062851, 3.659432, 1.651065, -1.103480, -4.659509
ptzero2v1=
-0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000, 0.000000, 5.092310, 5.335310, 0.000000
--7.125980, -10.940700, 3.124150, -5.600000, -5.434780, 2.326510, -18.818200, -1.521200, -2.955410, 2.316300, -0.641013, -7.180860, -2.633960, 4.182940, 3.540350, 6.843010
-7.705500, 2.497010, 4.755030, 8.680260, 9.939430, 8.246580, 6.258310, 6.410950, 8.326710, 12.216700, 8.361320, 7.109610, 4.429420, 9.940720, 10.031400, 9.626190
-16.148300, 8.402900, 8.574210, 9.887430, 8.434920, 8.565290, 8.827220, 7.485750, 8.360080, 10.124800, 8.100660, 9.172740, 9.754290, 7.803890, 7.765130, 2.989080
-4.471970, 4.843990, 2.961210, 2.757320, 0.446784, 2.244700, 2.224380, 5.063170, -2.196180, 0.000000, 5.043370, 2.638480, 5.035890, -0.083561, 6.014590, -0.178800
--2.936590, 0.924423, 5.063060, 1.800420, 0.567349, -1.334190, 0.000000, 2.203800, -0.084168, 3.704220, 10.890600, 3.668250, 2.581830, 2.658140, 0.882729, -0.658854
-4.596010, 2.404060, 1.996760, 11.972500, 5.221670, 5.335070, 9.104860, 7.442970, -14.638800, -2.467740, 0.000000
+0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,0.000000,15.797927,0.000000,-26.941333,5.787848,-6.235832, 0.000000, 9.605241, -3.262354, 0.000000
+-3.441877, -16.075351, -6.572906, -5.724473, -9.497821, -10.043196, -2.364265, -9.045754, -4.574509, -3.832558, -2.720129, 0.034401, -4.888662, -4.339271, -0.739835, 9.003159
+10.538456, 2.889849, 4.840125, 7.060004, 5.461183, 5.326597, 4.129304, 4.792289, 5.922200, 4.720092, 6.011253, 9.244853, 5.717675, 5.626520, 4.252377, 9.084188
+13.477456, 9.676237, 5.041439, 4.927180, 8.815112, 5.471457, 4.852762, 5.415544, 5.324521, 5.627313, 5.305027, 5.532285, 4.649912, 5.189615, 4.901391, -0.458141
+-0.562961, 0.109839, -0.110310, 0.089921, 0.982263, -0.384807, -0.527779, -1.591894, -0.956725, 0.000000, 1.763876, 1.880695, 4.074187, 3.168092, 4.148799, -0.244918
+-0.251073, -0.242988, -4.432312, 0.264216, 2.270349, -1.683392, -0.765528, -1.876069, -0.262543, -1.418778, 6.369130, -1.086514, -0.127075, 0.274085, -1.187019, 1.176133
+4.480484, 1.892931, 3.462448, 5.157830, 3.195353, 5.002985, -2.113051, 2.115013, 2.610369, 3.877452, -15.224138
ptzero2x2=
--3.061820,1.242440,-0.920429,0.949141,-3.371390,1.538490,1.537160,2.680810,-0.375996,3.624970,4.023850,3.101700, 3.621450, 6.546410, -3.104770, -2.297460
--4.208040, -0.680161, -0.759493, -3.256850, -3.475440, -1.010530, -1.495540, -5.828560, 0.005329, -1.934660, -0.802933, 2.697560, 2.595770, 2.262840, 3.847550, 6.508350
-7.420170, 7.293470, 4.544480, 6.006310, 1.299380, -1.178930, 0.399639, -2.844170, -0.636414, 2.689540, 0.000000, -0.005122, 3.258850, 3.209370, 2.692630, 6.696770
-8.405310, 4.999500, 6.351480, 5.621890, 4.321230, 7.012510, 5.699240, 2.569630, 6.354140, 4.117900, -2.852930, 0.848214, 3.756270, 5.198940, 1.263700, 7.994620
-7.224540, 6.740530, 8.552430, 7.867670, 8.216770, 9.636760, 8.177560, 5.522230, 9.502410, 13.189900, 12.033100, 12.005700, 12.146600, 7.825520, 12.772000
+-3.393311,-4.612718,-2.167408,-4.088590,-6.346786,-0.544223,-0.861993,-1.376263,-0.931888,0.974320,-1.771465,-1.282101, -0.756078, 0.895922, -5.083300, -5.798590
+-5.238743, -5.962301, -5.209907, -6.647940, -5.376005, -5.030277, -5.105275, -4.989542, -4.869607, -4.577173, -3.996694, -1.455519, -0.109249, 0.135541, 0.122785, 3.231639
+3.955449, 0.820883, 1.328178, 0.856624, -1.419408, -0.203201, -4.858058, 0.486545, -0.204672, -0.628556, 0.000000, 1.080338, -0.207922, 2.318720, 0.069046, 5.719123
+4.660274, 3.852756, 4.781462, 4.377001, 0.654567, 4.633240, -0.198329, 4.472335, 4.162579, -0.243058, 0.381767, 0.014187, -0.496195, 0.210439, -0.658795, 5.156046
+3.561575, 5.296722, 3.962189, 6.506880, 5.076435, 4.957450, 5.241003, 4.843687, 9.614483, 8.657085, 9.717657, 9.287841, 9.359161, 9.076187, 14.776648
ptzero2x1=
--0.994652,0.323887,2.662040,-2.946220,0.944695,4.667010,2.647140,2.461640,4.957990,3.113680,4.142770,1.308310, 5.622420, 5.469640, 7.911730, 4.733910
-2.399680, 2.495830, 1.696250, -0.137557, -0.839917, 0.820595, -0.677687, 1.409470, 6.096500, 1.789550, 3.080510, 1.909060, 5.210970, 8.603190, 7.767040, 4.964880
-1.534830, 2.095010, 0.831891, 2.758970, 0.161976, -2.097360, -3.893230, -2.599900, 0.471164, -3.455090, -2.770680, -0.905358, -2.358650, -1.462630, -0.698122, -0.124938
--0.176906, -0.640411, 0.314966, 3.389660, 0.198043, 1.051300, 0.663665, 0.172697, -5.500470, 1.035880, -1.542160, -6.413390, -2.720510, -1.530320, -3.358420, -1.622080
-5.754260, 7.172000, 7.128620, 7.212770, 9.011560, 12.249800, 7.627180, 8.218560, 10.479100, 14.471400, 13.159500, 14.039200, 12.556300, 14.254100, 10.879100
+-0.663303,0.750004,-1.131372,-0.705505,-4.447790,-1.433070,-0.730344,-2.437121,3.683164,1.349678,3.949167,1.954866, 3.345526, -0.206389, 0.814097, -0.072029
+-4.315611, -1.551488, -5.709898, 0.217207, -4.397529, 0.112074, 0.887342, -0.300364, 0.097019, 0.145158, -2.244817, 0.072637, 5.367818, -0.717107, 5.878672, 0.307109
+0.346803, -0.697651, 0.153453, 0.353230, -4.563616, -5.845637, -4.716185, -5.739089, -5.225427, -4.785183, -5.045789, -3.444209, -3.970617, -1.310207, -0.714994, -4.604468
+-4.783984, -3.799670, 0.552776, -4.988872, -0.135217, 0.381891, -4.820884, -4.682457, -4.309765, -3.916776, -5.392778, -5.921871, -4.572518, -4.067077, -4.478334, -5.374981
+3.921848, 7.897123, 5.113699, 4.979102, 5.734252, 8.295544, 9.085984, 8.228026, 9.357422, 8.897633, 10.024837, 10.418947, 10.199800, 10.012554, 10.822783
ptzero2u2=
-0.000000,0.000000,0.000000,2.927810,1.513160,0.000000,0.082767,11.832700,0.289694,3.096540,10.310600,4.395130, -1.204390, 6.047720, 3.493290, 6.893580
-2.216220, 7.989400, 6.820680, 6.566040, 4.801940, 4.916770, -1.832630, 3.958420, 3.803680, 6.329900, 7.061990, 4.710200, 4.066710, 4.554470, 3.114060, 11.153400
-4.604940, 2.830040, 4.262040, 1.597520, 5.888230, 2.311250, 2.476890, 5.373290, 4.946830, 5.122540, 4.026250, 8.359670, -0.690787, 2.955530, 1.806160, -0.799365
-2.603230, -0.796656, 0.954418, -0.788097, -5.090030, -1.691130, 0.121203, -0.729663, -0.713128, -1.170500, 0.482169, -1.926040, 2.590880, 2.291740, 2.106530, 0.515922
-1.563030, -0.969506, 1.715190, 0.730038, 2.797450, 0.343180, 2.983660, 1.569480, 3.351930, 4.273950, 2.693030, 0.502254, 5.028720, 2.408280, 4.831420, 2.509090
-4.811720, 4.595160, 4.270700, 2.027170, 1.296140, 4.402360, 3.805610, 5.697640, 4.209350, 4.205540, 6.840550, 8.688300, 10.346200, 16.544000, 16.701800, 15.950600
-19.425700, -2.883330, 5.303570, 23.328300, -11.923100, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000
+0.000000,0.000000,-14.919180,-7.406128,-4.792694,-14.401015,3.323219,3.442643,1.992985,7.227504,3.569291,2.724853, 8.593646, 4.886981, 6.808398, -0.520637
+-0.170822, 4.938485, 0.093942, -3.273033, 2.980335, -0.557907, 0.070914, 0.522724, 0.757116, -0.891534, 4.115294, -0.062327, -0.220173, 4.527388, 2.563957, 4.865393
+-0.693500, 1.060291, 0.042111, 0.545020, 2.119336, -0.012555, 0.870881, -0.275448, -0.467212, 0.688319, 3.278168, 4.881899, -2.048968, -0.309803, -0.898869, -3.567386
+-0.501581, -1.824903, -4.975555, -5.133580, -5.335531, -5.165667, -1.120129, -1.015706, -0.996176, -4.611436, -4.249178, -3.225271, -0.085437, -0.728502, 0.002872, -0.187782
+0.382141, 0.312615, 0.162156, -3.135492, -0.068732, -1.076046, -0.800665, 0.966061, 0.611495, 0.724621, -5.611031, 0.049093, -1.721456, 3.378318, 0.009991, 3.580431
+3.406525, 1.134893, 2.525603, 3.095577, 0.732492, 1.831709, 3.164246, 4.780390, 5.077633, 3.689338, 5.732011, -0.091538, 9.814865, -0.005154, 12.473843, 13.401516
+7.629707, 10.990937, -25.484783, 6.148513, 7.199727, 7.699411, -20.800384, 0.000000, 0.000000, 0.000000, 0.000000
ptzero2u1=
-0.000000,0.000000,0.000000,0.000000,1.981130,6.887300,9.154860,-16.603200,-16.215000,-1.705760,8.661480,8.515390, 11.496200, 11.843500, 0.198014, -4.724240
-2.248020, -1.088080, 1.330110, -1.029570, 2.240130, 5.972840, -0.982127, 0.718305, -0.225920, -2.431140, -1.316000, 0.185943, 0.340701, 4.595960, 1.248220, 3.011960
-5.469650, 4.392520, -2.113370, 5.228590, 6.901420, 1.807750, 2.789800, 2.019910, 4.260250, 4.067850, 5.996700, 4.556910, 9.895290, 7.413430, 10.164000, 7.858780
-6.970080, 7.429900, 7.051170, 7.430280, 8.005980, 7.404120, 4.407060, 7.474910, 8.819040, 6.623950, 10.455900, 10.604900, 11.910200, 10.554000, 8.699290, -1.134720
-1.416000, 0.929777, 1.344440, -1.600020, -0.365316, 4.031480, 2.137260, 2.109900, 3.030740, 4.618980, 3.938440, 5.779070, 6.430440, 7.514210, 5.483380, 1.555490
-5.777830, 0.909257, 0.741250, 4.815550, 4.298590, 3.555880, 4.731660, 6.200930, -0.525888, 0.160183, -28.817600, 6.296690, 6.827100, 8.975990, 10.564200, 8.914800
-19.375800, 17.450700, -14.710500, 23.005800, 12.757800, 23.291100, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000
+0.000000,0.000000,-4.236657,11.479925,-23.909753,0.000000,1.265017,4.044591,1.066353,4.931711,4.249248,5.740276, 3.239080, 4.417615, 3.400639, -6.085589
+-1.190744, -3.290192, -0.900773, -4.238490, -0.434513, -1.278303, -1.527290, -1.342075, -1.842262, -3.708529, -0.719099, -0.718612, -1.342815, -1.299639, -0.915499, 2.705918
+0.060151, 0.705120, -0.201201, 0.014094, 4.958492, -0.415360, 0.313299, -0.224189, 5.125605, 4.575965, 3.657751, 4.951192, 4.850901, 4.499458, 8.837723, 9.013613
+4.691552, 7.786943, 8.021501, 4.634629, 5.026911, 4.633873, 5.273969, 8.974319, 5.533759, 5.855616, 8.923941, 5.417079, 9.213529, 9.118045, 8.925730, -0.673704
+-0.265132, 0.428069, -2.340199, -1.962257, -0.301580, -1.143035, -0.282083, -1.373554, 0.422814, -1.002997, -1.330944, 4.213611, 4.115524, 1.043258, 2.710389, -0.296118
+-0.277632, -0.649401, 2.389717, 4.570673, 1.151151, -0.260586, -0.724143, 3.098157, 2.811309, 2.157092, 4.525832, 1.288603, -5.493975, 3.756562, 5.825718, 12.645693
+8.619234, 9.777951, 12.869860, 14.297061, 10.122535, 14.864693, -10.130053, 8.509878, 0.000000, 0.000000, 0.000000
--
GitLab