Update production replay scripts

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);

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");

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);

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);
 
    
 

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);