From 01986eb60da48c87e5fd654586aaec199d184b4c Mon Sep 17 00:00:00 2001
From: hallc-online <hallc-online@jlab.org>
Date: Mon, 20 Mar 2017 11:53:38 -0400
Subject: [PATCH] Add beta cut.
---
CALIBRATION/shms_cal_calib/THcPShowerCalib.h | 11 +++++++++++
CALIBRATION/shms_cal_calib/pcal_calib.cpp | 1 -
2 files changed, 11 insertions(+), 1 deletion(-)
diff --git a/CALIBRATION/shms_cal_calib/THcPShowerCalib.h b/CALIBRATION/shms_cal_calib/THcPShowerCalib.h
index 6078a098..60f10bdb 100644
--- a/CALIBRATION/shms_cal_calib/THcPShowerCalib.h
+++ b/CALIBRATION/shms_cal_calib/THcPShowerCalib.h
@@ -34,6 +34,9 @@
#define PR_ADC_THR 0
#define SH_ADC_THR 0
+#define BETA_MIN 0.5
+#define BETA_MAX 1.5
+
//#define MAX_TRACKS 10
using namespace std;
@@ -101,6 +104,7 @@ class THcPShowerCalib {
Double_t P_tr_tg_dp;
Double_t P_hgcer_npe[4];
+ Double_t P_tr_beta;
TBranch* b_P_tr_p;
TBranch* b_P_pr_apos_p;
@@ -113,6 +117,7 @@ class THcPShowerCalib {
TBranch* b_P_tr_yp;
TBranch* b_P_tr_tg_dp;
TBranch* b_P_hgcer_npe;
+ TBranch* b_P_tr_beta;
// Quantities for calculations of the calibration constants.
@@ -221,6 +226,8 @@ void THcPShowerCalib::Init() {
fTree->SetBranchAddress("P.hgcer.npe", P_hgcer_npe,&b_P_hgcer_npe);
+ fTree->SetBranchAddress("P.tr.beta", &P_tr_beta,&b_P_tr_beta);
+
// Histogram declarations.
// hEunc = new TH1F("hEunc", "Edep/P uncalibrated", 500, 0., 5.);
@@ -408,6 +415,10 @@ bool THcPShowerCalib::ReadShRawTrack(THcPShTrack &trk, UInt_t ientry) {
P_hgcer_npe[3] > 2 ;
if(!good_hgcer) return 0;
+ bool good_beta = P_tr_beta > BETA_MIN &&
+ P_tr_beta < BETA_MAX ;
+ if(!good_beta) return 0;
+
// cout << " Track is good." << endl << endl;
// getchar();
diff --git a/CALIBRATION/shms_cal_calib/pcal_calib.cpp b/CALIBRATION/shms_cal_calib/pcal_calib.cpp
index 11440019..78f75a73 100644
--- a/CALIBRATION/shms_cal_calib/pcal_calib.cpp
+++ b/CALIBRATION/shms_cal_calib/pcal_calib.cpp
@@ -8,7 +8,6 @@
// A steering Root script for the SHMS calorimeter calibration.
//
-////void pcal_calib(Int_t RunNumber) {
void pcal_calib(string RunNumber) {
// Initialize the analysis clock
--
GitLab