From ff5871f746cb95b991bc94fb6714612ff3eb7602 Mon Sep 17 00:00:00 2001
From: Mark Jones <jones@jlab.org>
Date: Fri, 8 Dec 2017 15:00:32 -0500
Subject: [PATCH] Update production replay scripts
---
SCRIPTS/COIN/PRODUCTION/replay_production_coin.C | 8 ++++++++
SCRIPTS/HMS/PRODUCTION/replay_production_hms.C | 4 ++++
SCRIPTS/HMS/PRODUCTION/replay_production_hms_kpp.C | 10 ++++++++--
SCRIPTS/SHMS/PRODUCTION/replay_production_shms.C | 4 ++++
SCRIPTS/SHMS/PRODUCTION/replay_production_shms_kpp.C | 4 ++++
5 files changed, 28 insertions(+), 2 deletions(-)
diff --git a/SCRIPTS/COIN/PRODUCTION/replay_production_coin.C b/SCRIPTS/COIN/PRODUCTION/replay_production_coin.C
index 4f1087d3..ed3e8ba1 100644
--- a/SCRIPTS/COIN/PRODUCTION/replay_production_coin.C
+++ b/SCRIPTS/COIN/PRODUCTION/replay_production_coin.C
@@ -78,9 +78,13 @@ void replay_production_coin (Int_t RunNumber = 0, Int_t MaxEvent = 0) {
gHaApps->Add(pbeam);
THaReactionPoint* prp= new THaReactionPoint("P.react"," SHMS reaction point","P","P.rb");
gHaPhysics->Add(prp);
+ THcExtTarCor* pext = new THcExtTarCor("P.extcor"," HMS extended target corrections","P","P.react");
+ gHaPhysics->Add(pext);
// Add Physics Module to calculate primary (scattered) beam kinematics
THcPrimaryKine* pkin_primary = new THcPrimaryKine("P.kin.primary", "SHMS Single Arm Kinematics", "P", "P.rb");
gHaPhysics->Add(pkin_primary);
+ THcHodoEff* peff = new THcHodoEff("phodeff"," SHMS hodo efficiency","P.hod");
+ gHaPhysics->Add(peff);
// Add Physics Module to calculate secondary (scattered) beam kinematics
// THcSecondaryKine* pkin_secondary = new THcSecondaryKine("P.kin.secondary", "SHMS Single Arm Kinematics", "P", "H.kin.primary");
// gHaPhysics->Add(pkin_secondary);
@@ -124,12 +128,16 @@ void replay_production_coin (Int_t RunNumber = 0, Int_t MaxEvent = 0) {
gHaApps->Add(hbeam);
THaReactionPoint* hrp= new THaReactionPoint("H.react"," HMS reaction point","H","H.rb");
gHaPhysics->Add(hrp);
+ THcExtTarCor* hext = new THcExtTarCor("H.extcor"," HMS extended target corrections","H","H.react");
+ gHaPhysics->Add(hext);
// Add Physics Module to calculate primary (scattered) beam kinematics
// THcPrimaryKine* hkin_primary = new THcPrimaryKine("H.kin.primary", "HMS Single Arm Kinematics", "H", "H.rb");
// gHaPhysics->Add(hkin_primary);
// Add Physics Module to calculate secondary (scattered) beam kinematics
THcSecondaryKine* hkin_secondary = new THcSecondaryKine("H.kin.secondary", "HMS Single Arm Kinematics", "H", "P.kin.primary");
gHaPhysics->Add(hkin_secondary);
+ THcHodoEff* heff = new THcHodoEff("hhodeff"," HMS hodo efficiency","H.hod");
+ gHaPhysics->Add(heff);
// Add event handler for scaler events
THcScalerEvtHandler *hscaler = new THcScalerEvtHandler("H", "Hall C scaler event type 4");
hscaler->AddEvtType(2);
diff --git a/SCRIPTS/HMS/PRODUCTION/replay_production_hms.C b/SCRIPTS/HMS/PRODUCTION/replay_production_hms.C
index ed8f98d9..bff5c9c1 100644
--- a/SCRIPTS/HMS/PRODUCTION/replay_production_hms.C
+++ b/SCRIPTS/HMS/PRODUCTION/replay_production_hms.C
@@ -77,12 +77,16 @@ void replay_production_hms(Int_t RunNumber=0, Int_t MaxEvent=0) {
gHaApps->Add(beam);
THaReactionPoint* hrp= new THaReactionPoint("H.react"," HMS reaction point","H","H.rb");
gHaPhysics->Add(hrp);
+ THcExtTarCor* hext = new THcExtTarCor("H.extcor"," HMS extended target corrections","H","H.react");
+ gHaPhysics->Add(hext);
// Add Ideal Beam Apparatus
// THaApparatus* beam = new THaIdealBeam("IB", "Ideal Beamline");
// gHaApps->Add(beam);
// Add Physics Module to calculate primary (scattered) beam kinematics
THcPrimaryKine* hkin = new THcPrimaryKine("H.kin", "HMS Single Arm Kinematics", "H", "H.rb");
gHaPhysics->Add(hkin);
+ THcHodoEff* heff = new THcHodoEff("hhodeff"," HMS hodo efficiency","H.hod");
+ gHaPhysics->Add(heff);
// Add handler for prestart event 125.
THcConfigEvtHandler* ev125 = new THcConfigEvtHandler("HC", "Config Event type 125");
diff --git a/SCRIPTS/HMS/PRODUCTION/replay_production_hms_kpp.C b/SCRIPTS/HMS/PRODUCTION/replay_production_hms_kpp.C
index a4d8f0a4..cdca9822 100644
--- a/SCRIPTS/HMS/PRODUCTION/replay_production_hms_kpp.C
+++ b/SCRIPTS/HMS/PRODUCTION/replay_production_hms_kpp.C
@@ -72,10 +72,16 @@ void replay_production_hms_kpp(Int_t RunNumber=0, Int_t MaxEvent=0) {
THaGoldenTrack* gtr = new THaGoldenTrack("H.gtr", "HMS Golden Track", "H");
gHaPhysics->Add(gtr);
// Add Ideal Beam Apparatus
- THaApparatus* beam = new THaIdealBeam("IB", "Ideal Beamline");
+ // THaApparatus* beam = new THaIdealBeam("IB", "Ideal Beamline");
+// Add Rastered Beam Apparatus
+ THaApparatus* beam = new THcRasteredBeam("H.rb", "Rastered Beamline");
gHaApps->Add(beam);
+ THaReactionPoint* hrp= new THaReactionPoint("H.react"," HMS reaction point","H","H.rb");
+ gHaPhysics->Add(hrp);
+ THcExtTarCor* hext = new THcExtTarCor("H.extcor"," HMS extended target corrections","H","H.react");
+ gHaPhysics->Add(hext);
// Add Physics Module to calculate primary (scattered) beam kinematics
- THcPrimaryKine* hkin = new THcPrimaryKine("H.kin", "HMS Single Arm Kinematics", "H", "IB");
+ THcPrimaryKine* hkin = new THcPrimaryKine("H.kin", "HMS Single Arm Kinematics", "H", "H.rb");
gHaPhysics->Add(hkin);
THcHodoEff* heff = new THcHodoEff("hodeff"," HMS hodo efficiency","H.hod");
gHaPhysics->Add(heff);
diff --git a/SCRIPTS/SHMS/PRODUCTION/replay_production_shms.C b/SCRIPTS/SHMS/PRODUCTION/replay_production_shms.C
index f8f6900e..a1896596 100644
--- a/SCRIPTS/SHMS/PRODUCTION/replay_production_shms.C
+++ b/SCRIPTS/SHMS/PRODUCTION/replay_production_shms.C
@@ -83,9 +83,13 @@ void replay_production_shms (Int_t RunNumber = 0, Int_t MaxEvent = 0) {
gHaApps->Add(beam);
THaReactionPoint* prp= new THaReactionPoint("P.react"," SHMS reaction point","P","P.rb");
gHaPhysics->Add(prp);
+ THcExtTarCor* pext = new THcExtTarCor("P.extcor"," HMS extended target corrections","P","P.react");
+ gHaPhysics->Add(pext);
// Add Physics Module to calculate primary (scattered beam - usually electrons) kinematics
THcPrimaryKine* kin = new THcPrimaryKine("P.kin", "SHMS Single Arm Kinematics", "P", "P.rb");
gHaPhysics->Add(kin);
+ THcHodoEff* peff = new THcHodoEff("phodeff"," SHMS hodo efficiency","P.hod");
+ gHaPhysics->Add(peff);
diff --git a/SCRIPTS/SHMS/PRODUCTION/replay_production_shms_kpp.C b/SCRIPTS/SHMS/PRODUCTION/replay_production_shms_kpp.C
index b1065cd6..53629967 100644
--- a/SCRIPTS/SHMS/PRODUCTION/replay_production_shms_kpp.C
+++ b/SCRIPTS/SHMS/PRODUCTION/replay_production_shms_kpp.C
@@ -82,6 +82,10 @@ void replay_production_shms_kpp (Int_t RunNumber = 0, Int_t MaxEvent = 0) {
// Add Rastered Beam Apparatus
THaApparatus* beam = new THcRasteredBeam("P.rb", "Rastered Beamline");
gHaApps->Add(beam);
+ THaReactionPoint* prp= new THaReactionPoint("P.react"," SHMS reaction point","P","P.rb");
+ gHaPhysics->Add(prp);
+ THcExtTarCor* pext = new THcExtTarCor("P.extcor"," HMS extended target corrections","P","P.react");
+ gHaPhysics->Add(pext);
// Add Physics Module to calculate primary (scattered beam - usually electrons) kinematics
THcPrimaryKine* kin = new THcPrimaryKine("P.kin", "SHMS Single Arm Kinematics", "P", "P.rb");
gHaPhysics->Add(kin);
--
GitLab