Remove ALPACA generator

CHANGELOG.md

@@ -24,6 +24,7 @@ it in future.
 
 * Remove Millepede
 * Remove outdated example shipEvent_ex.py
+* Remove ALPACA generator
 
 ## 24.11
 

macro/run_simScript.py

@@ -3,7 +3,6 @@
 import sys
 import ROOT
 ROOT.gSystem.Load('libEGPythia8')
-import makeALPACAEvents
 
 import shipunit as u
 import shipRoot_conf
@@ -86,7 +85,6 @@
 parser.add_argument("--Genie",   dest="genie",   help="Genie for reading and processing neutrino interactions", required=False, action="store_true")
 parser.add_argument("--NuRadio", dest="nuradio", help="misuse GenieGenerator for neutrino radiography and geometry timing test", required=False, action="store_true")
 parser.add_argument("--Ntuple",  dest="ntuple",  help="Use ntuple as input", required=False, action="store_true")
-group.add_argument("--ALPACA",  dest="ALPACA",  help="Use ALPACA as input", required=False, action="store_true")
 parser.add_argument("--MuonBack",dest="muonback",  help="Generate events from muon background file, --Cosmics=0 for cosmic generator data", required=False, action="store_true")
 parser.add_argument("--FollowMuon",dest="followMuon", help="Make muonshield active to follow muons", required=False, action="store_true")
 parser.add_argument("--FastMuon",  dest="fastMuon",  help="Only transport muons for a fast muon only background estimate", required=False, action="store_true")
@@ -155,7 +153,6 @@
 if options.genie:    simEngine = "Genie"
 if options.nuradio:  simEngine = "nuRadiography"
 if options.ntuple:   simEngine = "Ntuple"
-if options.ALPACA:   simEngine = "ALPACA"
 if options.muonback: simEngine = "MuonBack"
 if options.nuage:    simEngine = "Nuage"
 if options.mudis:    simEngine = "muonDIS"
@@ -327,24 +324,6 @@
 # P8gen.SetId(-211)
  primGen.AddGenerator(P8gen)
 
-if simEngine == "ALPACA":
-  print('Generating ALP events of mass {} GeV with the photon coupling coefficient {} GeV^-1'.format(options.theMass, options.theDPepsilon))
-  target     = ship_geo.target
-  startZ     = target.z0
-  lengthZ    = target.length
-  endZ       = startZ + lengthZ
-  SmearBeam  = 1*u.cm # finite beam size
-  Lmin       = ((ship_geo.Chamber1.z - ship_geo.chambers.Tub1length) - ship_geo.target.z0)/100.
-  Lmax       = (ship_geo.TrackStation1.z - ship_geo.target.z0)/100.
-  print('ALPACA is initializing.')
-  inputFile  = makeALPACAEvents.runEvents(options.theMass,options.theDPepsilon,options.nEvents,Lmin,Lmax,startZ,endZ,SmearBeam)
-  if inputFile: print('ALPACA is done.')
-  ut.checkFileExists(inputFile)
-  ALPACAgen = ROOT.ALPACAGenerator()
-  ALPACAgen.Init(inputFile)
-  print('ALPACAGenerator is reading the ALPACA events')
-  primGen.AddGenerator(ALPACAgen)
-
 if simEngine == "FixedTarget":
  P8gen = ROOT.FixedTargetGenerator()
  P8gen.SetTarget("volTarget_1",0.,0.)

shipgen/CMakeLists.txt

@@ -46,7 +46,6 @@ MuonBackGenerator.cxx
 CosmicsGenerator.cxx
 MuDISGenerator.cxx
 FixedTargetGenerator.cxx
-ALPACAGenerator.cxx
 EvtCalcGenerator.cxx
 )
 

shipgen/GenLinkDef.h

@@ -15,6 +15,5 @@
 #pragma link C++ class  MuDISGenerator+;
 #pragma link C++ class  NuageGenerator+;
 #pragma link C++ class  FixedTargetGenerator+;
-#pragma link C++ class  ALPACAGenerator+;
 #pragma link C++ class  EvtCalcGenerator+;
 #endif