From 4777da884a55ef1cfd67d3908c00d0214ee143c1 Mon Sep 17 00:00:00 2001
From: Yero1990 <cyero002@fiu.edu>
Date: Tue, 6 Jun 2017 12:02:32 -0400
Subject: [PATCH] Files: Changes Made: (#131)
run_Cal.C commented out some boolean statements producing
that were producing a compilaton error. Commented out
execution of script update_pdcparam.C, as it is obsolete,
since we will be using a new paramerer file containing
individual wire tzero.
get_pdc_time_histo.C commented out obsolete lines that created some directories
wire_drift_times.C declared 'wirenum' variable outside of loop to solve the
"use of undeclared identifier" compilation error
get_wire_tzero.C moved several variable declarations outside the loop
over sensewires to solve "use of undeclared identifier"
compilation error
get_LookUp_Values.C minor syntax error fix.
---
CALIBRATION/shms_dc_calib/run_Cal.C | 28 ++++---
.../shms_dc_calib/scripts/get_LookUp_Values.C | 2 +-
.../scripts/get_pdc_time_histo.C | 4 +-
.../shms_dc_calib/scripts/get_wire_tzero.C | 75 +++++++++++--------
.../shms_dc_calib/scripts/wire_drift_times.C | 5 +-
5 files changed, 67 insertions(+), 47 deletions(-)
diff --git a/CALIBRATION/shms_dc_calib/run_Cal.C b/CALIBRATION/shms_dc_calib/run_Cal.C
index 8d08d80b..f16ff2dd 100644
--- a/CALIBRATION/shms_dc_calib/run_Cal.C
+++ b/CALIBRATION/shms_dc_calib/run_Cal.C
@@ -1,8 +1,11 @@
//SCRIPT TO RUN OVER ALL HMS DC CALIBRATION SCRIPTS AT ONCE, AND UPDATE THE
//NECESSARY PARAMTER FILES hdriftmap.param and hdc.param
+
+
+using namespace std;
void run_Cal()
{
-
+
//User Input Run
int run_NUM;
cout << "Enter Run Number: " << endl;
@@ -13,25 +16,26 @@ void run_Cal()
fout.open("scripts/input_RUN.txt");
fout << run_NUM << endl;
fout.close();
-
+
//Create root and data files Directories if they dont exist
- char *dir_root = "mkdir ./root_files/";
- char *dir_data = "mkdir ./data_files/";
+ const char *dir_root = "mkdir ./root_files/";
+ const char *dir_data = "mkdir ./data_files/";
- if (system(dir_root || dir_data) != 0) {
+ // if (std::system(dir_root || dir_data) != 0) {
system(dir_root);
system(dir_data);
- }
+ // }
+
//Create run Directories if they dont exist
char *dir0 = Form("mkdir ./root_files/run%d", run_NUM);
char *dir1 = Form("mkdir ./data_files/run%d", run_NUM);
- if (system(dir0 || dir1) != 0) {
+ //if (system(dir0 || dir1) != 0) {
system(dir0);
system(dir1);
- }
+ // }
@@ -48,21 +52,21 @@ void run_Cal()
-gSystem->cd("./scripts");
+ gSystem->cd("./scripts");
//execute code to get t0 from each wire in each plane
gSystem->Exec("root -l -q -b get_wire_tzero.C");
//execute code to update pdc parameter file
- gSystem->Exec("root -l -q update_pdcparam.C");
+ // gSystem->Exec("root -l -q update_pdcparam.C");
//execute code to get t0 corrected drift times
gSystem->Exec("root -l -q get_pdc_time_histo_tzero_corrected.C");
//execute code to update LookUp Table
gSystem->Exec("root -l -q get_LookUp_Values.C");
-
-
}
+
+
diff --git a/CALIBRATION/shms_dc_calib/scripts/get_LookUp_Values.C b/CALIBRATION/shms_dc_calib/scripts/get_LookUp_Values.C
index 52f7457e..8e0515e8 100644
--- a/CALIBRATION/shms_dc_calib/scripts/get_LookUp_Values.C
+++ b/CALIBRATION/shms_dc_calib/scripts/get_LookUp_Values.C
@@ -98,7 +98,7 @@ void get_LookUp_Values() {
lookup_value[ip] = binSUM[ip] / binContent_TOTAL[ip];
bin_count = bin_count + 1;
- if (bin_count < = 8 ) {
+ if (bin_count <= 8 ) {
ofs << setprecision(5) << lookup_value[ip] << fixed << ",";
}
diff --git a/CALIBRATION/shms_dc_calib/scripts/get_pdc_time_histo.C b/CALIBRATION/shms_dc_calib/scripts/get_pdc_time_histo.C
index 87669b2b..2c5e4984 100644
--- a/CALIBRATION/shms_dc_calib/scripts/get_pdc_time_histo.C
+++ b/CALIBRATION/shms_dc_calib/scripts/get_pdc_time_histo.C
@@ -14,14 +14,16 @@ void get_pdc_time_histo()
TString file_name = "../../../ROOTfiles/shms_replay_%d.root";
+ /*
//Create RUN Directories if they dont exist
- char *dir0 = Form("mkdir -p ../root_files/run%d", run_NUM);
+ char *dir0 = Form("mkdir -p ../root_files/run%d", run_NUM);
char *dir1 = Form("mkdir -p ../data_files/run%d", run_NUM);
if (system(dir0 || dir1) != 0) {
system(dir0);
system(dir1);
}
+ */
//open file
TFile *f = new TFile(Form(file_name.Data(), run_NUM), "READ");
diff --git a/CALIBRATION/shms_dc_calib/scripts/get_wire_tzero.C b/CALIBRATION/shms_dc_calib/scripts/get_wire_tzero.C
index 1c6974e2..592d9e06 100644
--- a/CALIBRATION/shms_dc_calib/scripts/get_wire_tzero.C
+++ b/CALIBRATION/shms_dc_calib/scripts/get_wire_tzero.C
@@ -5,7 +5,7 @@
a certain number of bins and this fit is extrapolated to y=0(x-axis). The extrapolated value is take to be t0*/
#include <vector>
-#include <TMath>
+#include "TMath.h"
#define NPLANES 12
@@ -34,7 +34,7 @@ void get_wire_tzero()
TString root_file;
TFile *f[NPLANES];
- int total_wires; //integer to store total sense wires for a plane chosen by the user
+ // int total_wires; //integer to store total sense wires for a plane chosen by the user
//Loop over all planes
for (int ip = 0; ip < NPLANES; ip++){
@@ -57,40 +57,61 @@ void get_wire_tzero()
f[ip]->cd(); //change to file containing the wire drift times histos
- int total_wires; //integer to store total sense wires for a plane chosen by the user
-
+ // int total_wires; //integer to store total sense wires for a plane chosen by the user
+
+ //INITIALIZE VARIABLES
+ int total_wires;
+ int sensewire;
+ TH1F *cell_dt[107];
+ Int_t *bin_max;
+ Int_t *bin_maxContent;
+ Double_t *time_max;
+ Double_t *twenty_perc_maxContent;
+ Double_t *ref_time;
+
+
+
+ //Declarations
+ int time_init; //start fit value
+ int time_final; //end fit value
+ int t_zero;
+ int entries; //entries for each wire
+
+ double m; //slope
+ double y_int; //y-intercept
+ double m_err;
+ double y_int_err;
+ double t_zero_err;
+
//Set variables depending on which plane is being studied
if(ip == 0 || ip == 1 || ip == 4 || ip == 5 || ip == 6 || ip == 7 || ip == 10 || ip == 11) {
- TH1F *cell_dt[107]; //declare array of histos to store drift times
+ //declare array of histos to store drift times
total_wires=107;
//Declare bin properties for given sense wires in a plane
- int bin_max[107]; /*Array to store the bin number corresponding to the drift time distribution peak*/
- int bin_maxContent[107]; /*Array to store the content (# events) corresponding to the bin with maximum content*/
- double time_max[107]; /*Array to store the x-axis(drift time (ns)) corresponding to bin_max*/
- double twenty_perc_maxContent[107]; /*Array to store 20% of maximum bin content (peak)*/
- double ref_time[107]; /*Array to store reference times for each sense wire*/
+ bin_max = new Int_t[total_wires]; /*Array to store the bin number corresponding to the drift time distribution peak*/
+ bin_maxContent= new Int_t[total_wires]; /*Array to store the content (# events) corresponding to the bin with maximum content*/
+ time_max= new Double_t[total_wires]; /*Array to store the x-axis(drift time (ns)) corresponding to bin_max*/
+ twenty_perc_maxContent= new Double_t[total_wires]; /*Array to store 20% of maximum bin content (peak)*/
+ ref_time= new Double_t[total_wires]; /*Array to store reference times for each sense wire*/
}
-
+
else if(ip == 2 || ip == 3 || ip == 8 || ip == 9) {
- TH1F *cell_dt[79];
total_wires=79;
-
- int bin_max[79];
- int bin_maxContent[79];
- double time_max[79];
- double twenty_perc_maxContent[79];
- double ref_time[79];
-
+ bin_max = new Int_t[total_wires];
+ bin_maxContent= new Int_t[total_wires];
+ time_max= new Double_t[total_wires];
+ twenty_perc_maxContent= new Double_t[total_wires];
+ ref_time= new Double_t[total_wires];
}
/*Get wire histos from root file and loop over each
sense wire of a plane in shms Drift Chambers (DC1 or DC2)*/
- for (int sensewire=1; sensewire<=total_wires; sensewire++){
+ for (sensewire=1; sensewire<=total_wires; sensewire++){
//Get title of histos in root file
TString drift_time_histo = Form("wire_%d", sensewire);
@@ -153,7 +174,7 @@ void get_wire_tzero()
// Loop over 2 bin contents stored in array content
for (j=0; j<2; j++){
- if(content[j] > = twenty_perc_maxContent[sensewire-1]){
+ if(content[j] >= twenty_perc_maxContent[sensewire-1]){
counts = counts+1;
if(counts >= 2) { goto stop;}
@@ -184,17 +205,7 @@ void get_wire_tzero()
//*******Extract the "t0" Using a Fitting Procedure********//
//*********************************************************//
- //Declarations
- int time_init; //start fit value
- int time_final; //end fit value
- int t_zero;
- int entries; //entries for each wire
-
- double m; //slope
- double y_int; //y-intercept
- double m_err;
- double y_int_err;
- double t_zero_err;
+
//Get time corresponding to bin (fit range)
time_init = cell_dt[sensewire-1] -> GetXaxis() -> GetBinCenter(bin_num[0]-5); //choose bin range over which to fit
diff --git a/CALIBRATION/shms_dc_calib/scripts/wire_drift_times.C b/CALIBRATION/shms_dc_calib/scripts/wire_drift_times.C
index 7d8cb945..eaf1ade2 100644
--- a/CALIBRATION/shms_dc_calib/scripts/wire_drift_times.C
+++ b/CALIBRATION/shms_dc_calib/scripts/wire_drift_times.C
@@ -66,6 +66,9 @@ void wire_drift_times::Loop()
//If specific planes are encountered, treat them as follows:
+ //variable declarations
+ int wirenum;
+
/*PLANE 1U1, 1V1, 2U1, 2V1*/
//If specific planes are encountered, treat them as follows:
if(ip == 0 || ip == 1 || ip == 4 || ip == 5 || ip == 6 || ip == 7 || ip == 10 || ip == 11) {
@@ -74,7 +77,7 @@ void wire_drift_times::Loop()
TH2F *wire_vs_dt = new TH2F("wire_vs_dt", "", 200., -50., 350., 107., 0.,107.);
//Initialize wire drift time histograms
- for (int wirenum=1; wirenum<=total_wires_uv; wirenum++){
+ for (wirenum=1; wirenum<=total_wires_uv; wirenum++){
cell_dt[wirenum-1] = new TH1F(Form("wire_%d", wirenum), "", 200., -50., 350.);
}
--
GitLab