From 0abb4fa6861c30ad91f06179398c62a304e4f0b2 Mon Sep 17 00:00:00 2001
From: Yero1990 <cyero002@fiu.edu>
Date: Fri, 9 Jun 2017 07:04:34 -0400
Subject: [PATCH] -added a few lines of code to produce a weighted tzero data
file (#133)
per wire basis, where the weight is by groups (ribbon cables)
---
.../shms_dc_calib/scripts/get_tzero_param.C | 20 ++++++++++++++-----
1 file changed, 15 insertions(+), 5 deletions(-)
diff --git a/CALIBRATION/shms_dc_calib/scripts/get_tzero_param.C b/CALIBRATION/shms_dc_calib/scripts/get_tzero_param.C
index 6feff232..6f8c06aa 100644
--- a/CALIBRATION/shms_dc_calib/scripts/get_tzero_param.C
+++ b/CALIBRATION/shms_dc_calib/scripts/get_tzero_param.C
@@ -92,7 +92,7 @@ for (ip=0; ip<NPLANES; ip++){
dc_group_max[ip] = new Int_t[group_size[ip]];
}
-//Fill the 2d array with 1st sensewire corresponding to each group in a given plane
+//Fill the 2d array with last sensewire corresponding to each group in a given plane
for (ip=0; ip<NPLANES; ip++) {
@@ -234,8 +234,12 @@ for (ip=0; ip<NPLANES; ip++){
} //end loop over planes
+ //create output file stream to write tzero values to a data file
+ ofstream out;
+ TString tzero_file = "../data_files/"+run+"/tzero_group.dat";
+ out.open(tzero_file);
- //Create an output file to store tzero values values
+ //Create an output parameter file to store tzero values
ofstream ofs;
TString wire_tzero = Form("./pdc_tzero_per_wire_run%d.param", run_NUM);
ofs.open (wire_tzero);
@@ -245,17 +249,18 @@ for (ip=0; ip<NPLANES; ip++){
{
//write plane headers
ofs << "ptzero"+planes[ip] << "=" << endl;
-
+ out << "#PLANE: " << planes[ip] << endl;
//loop over group
for (group = 0; group < group_size[ip]; group++)
{
-
+ out << "#GROUP: " << group+1 << endl;
+
//loop over sensewires
for (sw = 1; sw<=fNWires[ip]; sw++) {
if (sw >= dc_group_min[ip][group] && sw <= dc_group_max[ip][group])
{
-
+ out << setprecision(7) << weighted_avg[ip][group] << endl;
if (sw<fNWires[ip])
{
ofs << setprecision(7) << weighted_avg[ip][group] << ((sw+1) % 16 ? ", " : "\n") << fixed;
@@ -271,6 +276,11 @@ for (ip=0; ip<NPLANES; ip++){
}// end loop over groups
}//end loop over planes
+
+ //close files
+ ofs.close();
+ out.close();
+
}
--
GitLab