From 9e813dd7ad7b7b456b11e67170025effeced5593 Mon Sep 17 00:00:00 2001
From: Mark Jones <jones@jlab.org>
Date: Mon, 20 Aug 2018 20:13:18 -0400
Subject: [PATCH] Add BPM epics X and Y histograms
1) Add BPM EPICS X and Y histograms
to HMS/PRODUCTION/RASTER/hraster_histos.def
and SHMS/PRODUCTION/RASTER/praster_histos.def
2) Add script UTIL/BEAMLINE/plot_beam.C
to plot the mean of the BPM histograms
as a function of Z_beam to show the
beam track on target.
3) Add call plot_beam.C to the production
config files so that it is included in
the online analysis.
---
.../HMS/PRODUCTION/RASTER/hraster_histos.def | 13 +++++++-
.../SHMS/PRODUCTION/RASTER/praster_histos.def | 12 ++++++-
.../HMS/PRODUCTION/hms_coin_production.cfg | 5 +++
.../CONFIG/HMS/PRODUCTION/hms_production.cfg | 5 +++
.../SHMS/PRODUCTION/shms_coin_production.cfg | 5 +++
.../SHMS/PRODUCTION/shms_production.cfg | 5 +++
onlineGUI/UTIL/BEAMLINE/plot_beam.C | 33 +++++++++++++++++++
7 files changed, 76 insertions(+), 2 deletions(-)
create mode 100644 onlineGUI/UTIL/BEAMLINE/plot_beam.C
diff --git a/DEF-files/HMS/PRODUCTION/RASTER/hraster_histos.def b/DEF-files/HMS/PRODUCTION/RASTER/hraster_histos.def
index 43c5a239..9a06630f 100644
--- a/DEF-files/HMS/PRODUCTION/RASTER/hraster_histos.def
+++ b/DEF-files/HMS/PRODUCTION/RASTER/hraster_histos.def
@@ -34,4 +34,15 @@ TH2F hFRAraw_XvsY 'HMS FRA Raw X vs Y; FRXA RAW ADC (Volts); FRYA RAW ADC (Volts
TH2F hFRBraw_XvsY 'HMS FRB Raw X vs Y; FRXB RAW ADC (Volts); FRYB RAW ADC (Volts)' FRXBraw_V FRYBraw_V 100 0 1.0 100 0 1.0 H.dc.ntrack>0
TH2F hFRApos_XvsY 'HMS FRA X vs Y; FRA X Position (cm); FRA Y Position (cm)' H.rb.raster.fr_xa H.rb.raster.fr_ya 40 -0.2 0.2 40 -0.2 0.2 H.dc.ntrack>0
-TH2F hFRBpos_XvsY 'HMS FRB X vs Y; FRB X Position (cm); FRB Y Position (cm)' H.rb.raster.fr_xb H.rb.raster.fr_yb 40 -0.2 0.2 40 -0.2 0.2 H.dc.ntrack>0
\ No newline at end of file
+TH2F hFRBpos_XvsY 'HMS FRB X vs Y; FRB X Position (cm); FRB Y Position (cm)' H.rb.raster.fr_xb H.rb.raster.fr_yb 40 -0.2 0.2 40 -0.2 0.2 H.dc.ntrack>0
+
+TH1F hBPMA_X 'Beam BPMA_X ; BPMA_X (CM); Counts' H.rb.raster.fr_xbpmA 300 -.3 .3
+TH1F hBPMB_X 'Beam BPMB_X ; BPMB_X (CM); Counts' H.rb.raster.fr_xbpmB 300 -.3 .3
+TH1F hBPMC_X 'Beam BPMC_X ; BPMC_X (CM); Counts' H.rb.raster.fr_xbpmC 300 -.3 .3
+TH1F hBPMT_X 'Beam BPMT_X ; BPM_X projected to target (CM); Counts' H.rb.raster.fr_xbpm_tar 300 -.3 .3
+TH1F hBPMA_Y 'Beam BPMA_Y ; BPMA_Y (CM); Counts' H.rb.raster.fr_ybpmA 300 -.3 .3
+TH1F hBPMB_Y 'Beam BPMB_Y ; BPMB_Y (CM); Counts' H.rb.raster.fr_ybpmB 300 -.3 .3
+TH1F hBPMC_Y 'Beam BPMC_Y ; BPMC_Y (CM); Counts' H.rb.raster.fr_ybpmC 300 -.3 .3
+TH1F hBPMT_Y 'Beam BPMT_Y ; BPM_Y projected to target (CM); Counts' H.rb.raster.fr_ybpm_tar 300 -.3 .3
+
+
diff --git a/DEF-files/SHMS/PRODUCTION/RASTER/praster_histos.def b/DEF-files/SHMS/PRODUCTION/RASTER/praster_histos.def
index 20bf847d..37f53cd8 100644
--- a/DEF-files/SHMS/PRODUCTION/RASTER/praster_histos.def
+++ b/DEF-files/SHMS/PRODUCTION/RASTER/praster_histos.def
@@ -34,4 +34,14 @@ TH2F pFRAraw_XvsY 'SHMS FRA Raw X vs Y; FRXA RAW ADC (Volts); FRYA RAW ADC (Volt
TH2F pFRBraw_XvsY 'SHMS FRB Raw X vs Y; FRXB RAW ADC (Volts); FRYB RAW ADC (Volts)' FRXBraw_V FRYBraw_V 100 0 1.0 100 0 1.0 P.dc.ntrack>0
TH2F pFRApos_XvsY 'SHMS FRA X vs Y; FRA X Position (cm); FRA Y Position (cm)' P.rb.raster.fr_xa P.rb.raster.fr_ya 40 -0.2 0.2 40 -0.2 0.2 P.dc.ntrack>0
-TH2F pFRBpos_XvsY 'SHMS FRB X vs Y; FRB X Position (cm); FRB Y Position (cm)' P.rb.raster.fr_xb P.rb.raster.fr_yb 40 -0.2 0.2 40 -0.2 0.2 P.dc.ntrack>0
\ No newline at end of file
+TH2F pFRBpos_XvsY 'SHMS FRB X vs Y; FRB X Position (cm); FRB Y Position (cm)' P.rb.raster.fr_xb P.rb.raster.fr_yb 40 -0.2 0.2 40 -0.2 0.2 P.dc.ntrack>0
+
+TH1F hBPMA_X 'Beam BPMA_X ; BPMA_X (CM); Counts' P.rb.raster.fr_xbpmA 300 -.3 .3
+TH1F hBPMB_X 'Beam BPMB_X ; BPMB_X (CM); Counts' P.rb.raster.fr_xbpmB 300 -.3 .3
+TH1F hBPMC_X 'Beam BPMC_X ; BPMC_X (CM); Counts' P.rb.raster.fr_xbpmC 300 -.3 .3
+TH1F hBPMT_X 'Beam BPMT_X ; BPM_X projected to target (CM); Counts' P.rb.raster.fr_xbpm_tar 300 -.3 .3
+TH1F hBPMA_Y 'Beam BPMA_Y ; BPMA_Y (CM); Counts' P.rb.raster.fr_ybpmA 300 -.3 .3
+TH1F hBPMB_Y 'Beam BPMB_Y ; BPMB_Y (CM); Counts' P.rb.raster.fr_ybpmB 300 -.3 .3
+TH1F hBPMC_Y 'Beam BPMC_Y ; BPMC_Y (CM); Counts' P.rb.raster.fr_ybpmC 300 -.3 .3
+TH1F hBPMT_Y 'Beam BPMT_Y ; BPM_Y projected to target (CM); Counts' P.rb.raster.fr_ybpm_tar 300 -.3 .3
+
diff --git a/onlineGUI/CONFIG/HMS/PRODUCTION/hms_coin_production.cfg b/onlineGUI/CONFIG/HMS/PRODUCTION/hms_coin_production.cfg
index 54de8f81..a39ae506 100644
--- a/onlineGUI/CONFIG/HMS/PRODUCTION/hms_coin_production.cfg
+++ b/onlineGUI/CONFIG/HMS/PRODUCTION/hms_coin_production.cfg
@@ -142,6 +142,11 @@ title HMS Fast Raster
hFRApos_XvsY -nostat
hFRBpos_XvsY -nostat
+newpage 2 1
+title HMS EPICS BPM
+macro UTIL/BEAMLINE/plot_beam.C("X")
+macro UTIL/BEAMLINE/plot_beam.C("Y")
+
newpage 2 2
title HMS Kinematics
hdc_trk_mom -nostat
diff --git a/onlineGUI/CONFIG/HMS/PRODUCTION/hms_production.cfg b/onlineGUI/CONFIG/HMS/PRODUCTION/hms_production.cfg
index fc11477d..94efb217 100644
--- a/onlineGUI/CONFIG/HMS/PRODUCTION/hms_production.cfg
+++ b/onlineGUI/CONFIG/HMS/PRODUCTION/hms_production.cfg
@@ -142,6 +142,11 @@ title HMS Fast Raster
hFRApos_XvsY -nostat
hFRBpos_XvsY -nostat
+newpage 2 1
+title HMS EPICS BPM
+macro UTIL/BEAMLINE/plot_beam.C("X")
+macro UTIL/BEAMLINE/plot_beam.C("Y")
+
newpage 2 2
title HMS Kinematics
hdc_trk_mom -nostat
diff --git a/onlineGUI/CONFIG/SHMS/PRODUCTION/shms_coin_production.cfg b/onlineGUI/CONFIG/SHMS/PRODUCTION/shms_coin_production.cfg
index fffc6907..b0283f02 100644
--- a/onlineGUI/CONFIG/SHMS/PRODUCTION/shms_coin_production.cfg
+++ b/onlineGUI/CONFIG/SHMS/PRODUCTION/shms_coin_production.cfg
@@ -172,6 +172,11 @@ title SHMS Fast Raster
pFRApos_XvsY -nostat
pFRBpos_XvsY -nostat
+newpage 2 1
+title SHMS EPICS BPM
+macro UTIL/BEAMLINE/plot_beam.C("X")
+macro UTIL/BEAMLINE/plot_beam.C("Y")
+
newpage 2 2
title SHMS Kinematics
pdc_trk_mom -nostat
diff --git a/onlineGUI/CONFIG/SHMS/PRODUCTION/shms_production.cfg b/onlineGUI/CONFIG/SHMS/PRODUCTION/shms_production.cfg
index a53c4b2a..51f1bad3 100644
--- a/onlineGUI/CONFIG/SHMS/PRODUCTION/shms_production.cfg
+++ b/onlineGUI/CONFIG/SHMS/PRODUCTION/shms_production.cfg
@@ -172,6 +172,11 @@ title SHMS Fast Raster
pFRApos_XvsY -nostat
pFRBpos_XvsY -nostat
+newpage 2 1
+title SHMS EPICS BPM
+macro UTIL/BEAMLINE/plot_beam.C("X")
+macro UTIL/BEAMLINE/plot_beam.C("Y")
+
newpage 2 2
title SHMS Kinematics
pdc_trk_mom -nostat
diff --git a/onlineGUI/UTIL/BEAMLINE/plot_beam.C b/onlineGUI/UTIL/BEAMLINE/plot_beam.C
new file mode 100644
index 00000000..0c366514
--- /dev/null
+++ b/onlineGUI/UTIL/BEAMLINE/plot_beam.C
@@ -0,0 +1,33 @@
+void plot_beam(TString dir) {
+ TH1F* h1d[4];
+ if (dir=="Y") {
+ h1d[0] = (TH1F*) gDirectory->Get("hBPMA_Y");
+ h1d[1] = (TH1F*) gDirectory->Get("hBPMB_Y");
+ h1d[2] = (TH1F*) gDirectory->Get("hBPMC_Y");
+ h1d[3] = (TH1F*) gDirectory->Get("hBPMT_Y");
+ } else {
+ h1d[0] = (TH1F*) gDirectory->Get("hBPMA_X");
+ h1d[1] = (TH1F*) gDirectory->Get("hBPMB_X");
+ h1d[2] = (TH1F*) gDirectory->Get("hBPMC_X");
+ h1d[3] = (TH1F*) gDirectory->Get("hBPMT_X");
+ }
+ Double_t mean_x[4];
+ Double_t zbpm[4]={-320.82,-224.86,-129.44,0.};
+ for (Int_t nb=0;nb<4;nb++) {
+ mean_x[nb]=h1d[nb]->GetMean();
+ }
+ TGraph *grbeamx = new TGraph(4,zbpm,mean_x);
+ if (dir=="Y") grbeamx->SetTitle("; Z BPM (cm); Y BPM (cm) (+Y up) ");
+ if (dir=="X") grbeamx->SetTitle("; Z BPM (cm); X BPM (cm) (+Y up) ");
+ gPad->SetGridy();
+ gPad->SetGridx();
+ grbeamx->SetMarkerStyle(20);
+ grbeamx->SetMarkerSize(1);
+ grbeamx->GetXaxis()->SetLimits(-400,50);
+ grbeamx->SetMinimum(-.5);
+ grbeamx->SetMaximum(.5);
+ grbeamx->Draw("AP");
+
+
+ //
+}
--
GitLab