From 1f51bb6577ee7b2cad43cfd00c1aeb0993ec5b73 Mon Sep 17 00:00:00 2001
From: Yero1990 <cyero002@fiu.edu>
Date: Wed, 9 May 2018 07:04:01 -0400
Subject: [PATCH] Hms hodo calib work (#463)
* Added method to write hms hodo params into a .param file
* added runNUM argument to timeWalkCalib.C, so that it would be included in the output param filename :: Added a few lines to the hms output param file
* added the new hms hodo param file
---
CALIBRATION/hms_hodo_calib/timeWalkCalib.C | 123 ++++++++++++++++++++-
PARAM/HMS/HODO/hhodo_calib_1267.param | 81 ++++++++++++++
2 files changed, 200 insertions(+), 4 deletions(-)
create mode 100644 PARAM/HMS/HODO/hhodo_calib_1267.param
diff --git a/CALIBRATION/hms_hodo_calib/timeWalkCalib.C b/CALIBRATION/hms_hodo_calib/timeWalkCalib.C
index 27ccf0dc..f582a0ff 100644
--- a/CALIBRATION/hms_hodo_calib/timeWalkCalib.C
+++ b/CALIBRATION/hms_hodo_calib/timeWalkCalib.C
@@ -1,13 +1,18 @@
// Macro to perform time-walk fits and extract the calibration parameters
// Author: Eric Pooser, pooser@jlab.org
+
+
// Declare ROOT files
TFile *histoFile, *outFile;
+// Declare Output Parameter File
+ofstream outParam;
+
// Declare constants
static const UInt_t nPlanes = 4;
static const UInt_t nSides = 2;
-static const UInt_t nBarsMax = 21;
+static const UInt_t nBarsMax = 16;
static const UInt_t nTwFitPars = 3;
static const Double_t tdcThresh = 120.0; // 30 mV in units of FADC channels
@@ -17,6 +22,12 @@ static const Double_t c0twParInit = -15.0;
static const Double_t c1twParInit = 5.0;
static const Double_t c2twParInit = 0.25;
+//Parameter values to be written to param file
+Double_t c0[nPlanes][nSides][nBarsMax] = {0.};
+Double_t c1[nPlanes][nSides][nBarsMax] = {0.};
+Double_t c2[nPlanes][nSides][nBarsMax] = {0.};
+
+
static const Double_t fontSize = 0.05;
static const Double_t yTitleOffset = 0.75;
static const Double_t markerSize = 2.0;
@@ -222,11 +233,113 @@ void drawParams(UInt_t iplane) {
return;
} // drawParams
+//Add a method to Get Fit Parameters
+void WriteFitParam(int runNUM)
+{
+
+ TString outPar_Name = Form("../../PARAM/HMS/HODO/hhodo_calib_%d.param", runNUM);
+ outParam.open(outPar_Name);
+ outParam << ";HMS Hodoscopes Output Parameter File" << endl;
+ outParam << " " << endl;
+ outParam << " ; use htofusinginvadc=1 if want invadc_offset " << endl;
+ outParam << " ; invadc_linear, and invadc_adc to be used " << endl;
+ outParam << "htofusinginvadc=1 " << endl;
+ outParam << " " << endl;
+
+ //Fill 3D Par array
+ for (UInt_t iplane=0; iplane < nPlanes; iplane++)
+ {
+
+ for (UInt_t iside=0; iside < nSides; iside++) {
+
+
+ for(UInt_t ipaddle = 0; ipaddle < nbars[iplane]; ipaddle++) {
+
+ c0[iplane][iside][ipaddle] = twFit[iplane][iside][ipaddle]->GetParameter("c_{0}");
+ c1[iplane][iside][ipaddle] = twFit[iplane][iside][ipaddle]->GetParameter("c_{1}");
+ c2[iplane][iside][ipaddle] = twFit[iplane][iside][ipaddle]->GetParameter("c_{2}");
+
+ } //end paddle loop
+
+ } //end side loop
+
+ } //end plane loop
+
+ //Wrtie to Param FIle
+
+ outParam << ";Param c1-Pos" << endl;
+ outParam << "; " << setw(12) << "1x " << setw(15) << "1y " << setw(15) << "2x " << setw(15) << "2y " << endl;
+ outParam << "c1_Pos = ";
+ //Loop over all paddles
+ for(UInt_t ipaddle = 0; ipaddle < nBarsMax; ipaddle++) {
+ //Write c1-Pos values
+ if(ipaddle==0){
+ outParam << c1[0][0][ipaddle] << setw(15) << c1[1][0][ipaddle] << setw(15) << c1[2][0][ipaddle] << setw(15) << c1[3][0][ipaddle] << fixed << endl;
+ }
+ else {
+ outParam << setw(17) << c1[0][0][ipaddle] << setw(15) << c1[1][0][ipaddle] << setw(15) << c1[2][0][ipaddle] << setw(15) << c1[3][0][ipaddle] << fixed << endl;
+ }
+ } //end loop over paddles
+
+ outParam << " " << endl;
+ outParam << ";Param c1-Neg" << endl;
+ outParam << "; " << setw(12) << "1x " << setw(15) << "1y " << setw(15) << "2x " << setw(15) << "2y " << endl;
+ outParam << "c1_Neg = ";
+ //Loop over all paddles
+ for(UInt_t ipaddle = 0; ipaddle < nBarsMax; ipaddle++) {
+ //Write c1-Neg values
+ if(ipaddle==0){
+ outParam << c1[0][1][ipaddle] << setw(15) << c1[1][1][ipaddle] << setw(15) << c1[2][1][ipaddle] << setw(15) << c1[3][1][ipaddle] << fixed << endl;
+ }
+ else {
+ outParam << setw(17) << c1[0][1][ipaddle] << setw(15) << c1[1][1][ipaddle] << setw(15) << c1[2][1][ipaddle] << setw(15) << c1[3][1][ipaddle] << fixed << endl;
+ }
+} //end loop over paddles
+
+ outParam << " " << endl;
+ outParam << ";Param c2-Pos" << endl;
+ outParam << "; " << setw(12) << "1x " << setw(15) << "1y " << setw(15) << "2x " << setw(15) << "2y " << endl;
+ outParam << "c2_Pos = ";
+ //Loop over all paddles
+ for(UInt_t ipaddle = 0; ipaddle < nBarsMax; ipaddle++) {
+ //Write c2-Pos values
+ if(ipaddle==0)
+ {
+ outParam << c2[0][0][ipaddle] << setw(15) << c2[1][0][ipaddle] << setw(15) << c2[2][0][ipaddle] << setw(15) << c2[3][0][ipaddle] << fixed << endl;
+ }
+ else {
+ outParam << setw(17) << c2[0][0][ipaddle] << setw(15) << c2[1][0][ipaddle] << setw(15) << c2[2][0][ipaddle] << setw(15) << c2[3][0][ipaddle] << fixed << endl;
+ }
+ } //end loop over paddles
+
+ outParam << " " << endl;
+ outParam << ";Param c2-Neg" << endl;
+ outParam << "; " << setw(12) << "1x " << setw(15) << "1y " << setw(15) << "2x " << setw(15) << "2y " << endl;
+ outParam << "c2_Neg = ";
+ //Loop over all paddles
+ for(UInt_t ipaddle = 0; ipaddle < nBarsMax; ipaddle++) {
+ //Write c2-Neg values
+ if (ipaddle==0){
+ outParam << c2[0][1][ipaddle] << setw(15) << c2[1][1][ipaddle] << setw(15) << c2[2][1][ipaddle] << setw(15) << c2[3][1][ipaddle] << fixed << endl;
+ }
+ else{
+ outParam << setw(17) << c2[0][1][ipaddle] << setw(15) << c2[1][1][ipaddle] << setw(15) << c2[2][1][ipaddle] << setw(15) << c2[3][1][ipaddle] << fixed << endl;
+ }
+ } //end loop over paddles
+
+ outParam.close();
+} //end method
+
//=:=:=:=:=
//=: Main
//=:=:=:=:=
-void timeWalkCalib() {
+void timeWalkCalib(int run) {
+
+using namespace std;
+
+//prevent root from displaying graphs while executing
+ gROOT->SetBatch(1);
// ROOT settings
gStyle->SetTitleFontSize(fontSize);
@@ -285,9 +398,11 @@ void timeWalkCalib() {
// Draw the time-walk parameter graphs
drawParams(iplane);
} // Plane loop
- return;
-} // timeWalkFits()
+
+ //Write to a param file
+ WriteFitParam(run);
+}
diff --git a/PARAM/HMS/HODO/hhodo_calib_1267.param b/PARAM/HMS/HODO/hhodo_calib_1267.param
new file mode 100644
index 00000000..24bd6f72
--- /dev/null
+++ b/PARAM/HMS/HODO/hhodo_calib_1267.param
@@ -0,0 +1,81 @@
+;HMS Hodoscopes Output Parameter File
+
+ ; use htofusinginvadc=1 if want invadc_offset
+ ; invadc_linear, and invadc_adc to be used
+htofusinginvadc=1
+
+;Param c1-Pos
+; 1x 1y 2x 2y
+c1_Pos = 0.819027 1.41186 2.32872 1.46428
+ 2.207128 3.187273 2.074593 1.750816
+ 7.983453 2.126108 1.906823 2.117386
+ 3.073156 1.582985 2.086830 2.449768
+ 3.768713 4.136275 1.676022 1.394880
+ 3.509375 2.185954 2.514874 2.000962
+ 3.491312 2.763838 1.759023 3.858148
+ 2.440381 2.145233 2.139788 2.088941
+ 5.276700 0.404955 1.699047 2.313945
+ 1.771944 1.581550 1.682843 1.542734
+ 14.931448 0.000000 2.084104 0.000000
+ 9.956889 0.000000 3.557773 0.000000
+ 11.323196 0.000000 2.622266 0.000000
+ 2.716784 0.000000 2.258604 0.000000
+ 2.012462 0.000000 1.659209 0.000000
+ 1.380740 0.000000 1.363198 0.000000
+
+;Param c1-Neg
+; 1x 1y 2x 2y
+c1_Neg = 2.491850 1.219653 1.514545 1.359626
+ 1.417311 2.595391 2.340332 1.662741
+ 8.952472 2.093015 2.531239 2.902317
+ 4.904775 2.003444 3.784816 2.047449
+ 1.451639 2.527254 4.614227 1.227058
+ 16.113918 2.874847 5.681276 1.824491
+ 8.758572 2.298264 3.905578 1.153105
+ 3.328502 2.699598 5.064324 1.794633
+ 15.977368 2.265645 6.660858 1.658787
+ 5.869666 1.341041 2.633762 0.954160
+ 13.911405 0.000000 6.624087 0.000000
+ 7.222109 0.000000 3.399204 0.000000
+ 19.733347 0.000000 3.646857 0.000000
+ 3.554593 0.000000 3.719042 0.000000
+ 2.501336 0.000000 2.009574 0.000000
+ 1.514542 0.000000 3.041217 0.000000
+
+;Param c2-Pos
+; 1x 1y 2x 2y
+c2_Pos = 0.800818 0.688159 0.361689 0.547470
+ 0.424856 0.274696 0.414385 0.473044
+ 0.111285 0.409492 0.436868 0.378090
+ 0.313996 0.512042 0.407730 0.352738
+ 0.237399 0.212177 0.509222 0.604016
+ 0.265618 0.388077 0.353205 0.408538
+ 0.260990 0.313623 0.461489 0.219490
+ 0.400690 0.415163 0.408679 0.400332
+ 0.171046 1.400490 0.510330 0.392314
+ 0.511252 0.584952 0.498493 0.653343
+ 0.063082 0.000000 0.408709 0.000000
+ 0.087474 0.000000 0.255524 0.000000
+ 0.086283 0.000000 0.346336 0.000000
+ 0.324846 0.000000 0.393374 0.000000
+ 0.429733 0.000000 0.545845 0.000000
+ 0.596471 0.000000 0.573301 0.000000
+
+;Param c2-Neg
+; 1x 1y 2x 2y
+c2_Neg = 0.392238 0.639491 0.502620 0.587681
+ 0.594349 0.362502 0.372293 0.497877
+ 0.098630 0.401988 0.349340 0.305566
+ 0.197926 0.524730 0.245058 0.443191
+ 0.657146 0.370732 0.209873 0.652018
+ 0.061999 0.318930 0.158678 0.501406
+ 0.111300 0.390369 0.252101 0.716545
+ 0.283347 0.343220 0.190034 0.505662
+ 0.059642 0.420855 0.144001 0.523623
+ 0.164370 0.628903 0.357733 0.773707
+ 0.073184 0.000000 0.148840 0.000000
+ 0.140914 0.000000 0.271973 0.000000
+ 0.053692 0.000000 0.275449 0.000000
+ 0.271366 0.000000 0.256495 0.000000
+ 0.393127 0.000000 0.490331 0.000000
+ 0.584891 0.000000 0.290152 0.000000
--
GitLab