diff --git a/CALIBRATION/shms_dc_calib/README.md b/CALIBRATION/shms_dc_calib/README.md index e16f7ed1345a6cc624e66c24770da5a82888b4bd..725ac087546fdf2a78efb9465c4e5e1f2f4e2216 100644 --- a/CALIBRATION/shms_dc_calib/README.md +++ b/CALIBRATION/shms_dc_calib/README.md @@ -26,7 +26,7 @@ Running code * From the calibration results, two parameter files will be produced in: * hallc_replay/PARAM/SHMS/DC/pdc_tzero_per_wire_run%d_NEW.param. %d=run_number - * hallc_replay/PARAM/SHMS/DC/pdriftmap_new.param + * hallc_replay/PARAM/SHMS/DC/pdriftmap_run%d_NEW.param * Rename the new parameter files as follows: * copy: pdc_tzero_per_wire_run%d_NEW.param to pdc_tzero_per_wire.param diff --git a/CALIBRATION/shms_dc_calib/scripts/get_LookUp_Values.C b/CALIBRATION/shms_dc_calib/scripts/get_LookUp_Values.C index a0791bdd2126e19fbfbde909b1f97c428088bb5e..ab747052bf6b818d01e3bbd679784808ac57fa47 100644 --- a/CALIBRATION/shms_dc_calib/scripts/get_LookUp_Values.C +++ b/CALIBRATION/shms_dc_calib/scripts/get_LookUp_Values.C @@ -38,7 +38,7 @@ void get_LookUp_Values() { //Create an output file to store lookup values ofstream ofs; - TString lookup_table = "../../../PARAM/SHMS/DC/pdriftmap_new.param"; + TString lookup_table = Form("../../../PARAM/SHMS/DC/pdriftmap_run%d_NEW.param", run_NUM); ofs.open (lookup_table); diff --git a/CALIBRATION/shms_dc_calib/scripts/get_tzero_per_wire_param.C b/CALIBRATION/shms_dc_calib/scripts/get_tzero_per_wire_param.C index b6bb687d096c67bcb7ef76cc945d5e4cdf9f09f6..1633167400588da3f2e08d3d37a0aec7e2433ecb 100644 --- a/CALIBRATION/shms_dc_calib/scripts/get_tzero_per_wire_param.C +++ b/CALIBRATION/shms_dc_calib/scripts/get_tzero_per_wire_param.C @@ -69,15 +69,16 @@ Double_t **t0 = new Double_t*[NPLANES]; // cout << line << endl; if (line[0]!='#') { - if(!input.good()) {break;} + // if(!input.good()) {break;} sscanf(line.c_str(), "%d %lf", &nwire, &nt0); - if(sw>0) { - - t0[ip][nwire-1] = nt0; - - } - sw++; + + // if(sw>0) { + + t0[ip][nwire-1] = nt0; + + // } + // sw++; } @@ -111,18 +112,18 @@ Double_t **t0 = new Double_t*[NPLANES]; for (sw=0; sw<fNWires[ip]; sw++) { - out << setprecision(3) << t0[ip][sw] << endl; + out << setprecision(6) << t0[ip][sw] << endl; if (sw <= 10) { - ofs << setprecision(3) << t0[ip][sw] << fixed << ","; + ofs << setprecision(6) << t0[ip][sw] << fixed << ","; } else if (sw>10 && sw <(fNWires[ip]-1)) { - ofs << setprecision(3) << t0[ip][sw] << ((sw+1) % 16 ? ", " : "\n") << fixed; + ofs << setprecision(6) << t0[ip][sw] << ((sw+1) % 16 ? ", " : "\n") << fixed; } else if (sw==78 || sw == 106) { - ofs << setprecision(3) << t0[ip][sw] << fixed << endl; + ofs << setprecision(6) << t0[ip][sw] << fixed << endl; }