Skip to content
Snippets Groups Projects
Commit 0db4241b authored by Thomas Strebler's avatar Thomas Strebler Committed by Ahmed Tarek
Browse files

Clean up legacy ALFA_BeamTransport configs 

Clean up legacy ALFA_BeamTransport configs
parent 8f35817c
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 1 addition and 206 deletions
ForwardDetectors/ALFA/ALFA_BeamTransport/CMakeLists.txt
+ 1
− 4
View file @ 0db4241b
# Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
# Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
# Declare the package name:
atlas_subdir( ALFA_BeamTransport )
......@@ -12,6 +12,3 @@ atlas_add_component( ALFA_BeamTransport
src/components/*.cxx
INCLUDE_DIRS ${ROOT_INCLUDE_DIRS}
LINK_LIBRARIES ${ROOT_LIBRARIES} AtlasHepMCLib TruthUtils AthenaBaseComps FPTracker GaudiKernel StoreGateLib xAODEventInfo GeneratorObjects FPTracker )
# Install files from the package:
atlas_install_joboptions( share/*.py )
ForwardDetectors/ALFA/ALFA_BeamTransport/share/ALFA_BeamTransportConfig.py deleted 100644 → 0
+ 0
− 25
View file @ 8f35817c
import os
base_path = os.getenv('TwissFilesPATH') + '/0090.00m'
from ALFA_BeamTransport.ALFA_BeamTransportConf import ALFA_BeamTransport
myBeamTransport = ALFA_BeamTransport("ALFA_BeamTransport")
myBeamTransport.ConfDir = base_path
myBeamTransport.UseALFA = True # use the ALFA oftics hast to be true in oure case
myBeamTransport.Debug = True # write debug output
myBeamTransport.IP = 1 # which interaction point is used in our case always 1
myBeamTransport.UseAper = False
myBeamTransport.apermb = 0.
myBeamTransport.xcol1 = 999.
myBeamTransport.xcol2 = 999.
myBeamTransport.BeamEnergy = 3500. #set beam energy in GeV
myBeamTransport.RPDistance = 237.398 #RP position
myBeamTransport.absZMagMax = 500. # untill this point the magnets are read has to be bigger then the RP pos
myBeamTransport.EtaCut = 8.
myBeamTransport.XiCut = 0.2
myBeamTransport.FPOutputBeam1 = "FPTrackerBeam1.txt"
myBeamTransport.FPOutputBeam2 = "FPTrackerBeam2.txt"
myBeamTransport.McEventCollectionName = "TruthEvent"
topSeq += myBeamTransport
ForwardDetectors/ALFA/ALFA_BeamTransport/share/jobOptions.py deleted 100644 → 0
+ 0
− 143
View file @ 8f35817c
# -*- coding: utf-8 -*-
# -*- coding: utf-8 -*-
###############################################################
#
# JOB OPTIONS FILE. THIS FILE IS TO DEMONSTRATE HOW THE USER CAN
# CHANGE THE PYTHIA PARAMETERS. IT DOES NOT NECESARILY COMPLY TO
# THE ATLAS DEFAULTS FOR PHYSICS EVENT GENERATION.
#==============================================================
#--------------------------------------------------------------
# Private Application Configuration option
from AthenaCommon.AppMgr import ServiceMgr
ServiceMgr.MessageSvc.OutputLevel = DEBUG
#--------------------------------------------------------------
# Event related parameters
#--------------------------------------------------------------
# Number of events to be processed (default is 10)
theApp.EvtMax = 10000
#--------------------------------------------------------------
# Algorithms Private Options
#--------------------------------------------------------------
from RngComps.RngCompsConf import AtRndmGenSvc
ServiceMgr += AtRndmGenSvc()
from AthenaCommon.AlgSequence import AlgSequence
job=AlgSequence()
from Pythia8_i.Pythia8_iConf import Pythia8_i
job +=Pythia8_i()
include("pythia8.common.py")
include("pythia8.elastic.py")
#include("pythia8.single_diffraction.py")
#include("MC9.105001.pythia_minbias.py")
#sj#from TruthExamples.TruthExamplesConf import DumpMC
#sj#job += DumpMC()
# start my own
#from HepMCAnalysis_i.HepMCAnalysis_iConf import HepMCAnalysis
from ALFA_BeamTransport.ALFA_BeamTransportConf import ALFA_BeamTransport
myBeamTransport = ALFA_BeamTransport("ALFA_BeamTransport")
#myBeamTransport.somthing = what
myBeamTransport.ConfDir = "../../../FPTracker/share" #configuration folder with twiss files
myBeamTransport.UseALFA = True # use the ALFA oftics hast to be true in oure case
myBeamTransport.Debug = True # write debug output
myBeamTransport.IP = 1 # which interaction point is used in our case always 1
#some FPTracker suff don't change
myBeamTransport.UseAper = False
myBeamTransport.apermb = 0.
myBeamTransport.xcol1 = 999.
myBeamTransport.xcol2 = 999.
myBeamTransport.BeamEnergy = 7000. #set beam energy in GeV
myBeamTransport.RPDistance = 237.398 #RP position
myBeamTransport.absZMagMax = 500. # untill this point the magnets are read has to be bigger then the RP pos
#some cut for Production
myBeamTransport.EtaCut = 8.
myBeamTransport.XiCut = 0.2
#the two output files if debug i on
myBeamTransport.FPOutputBeam1 = "Beam1.dat1_new"
myBeamTransport.FPOutputBeam2 = "Beam2.dat1_new"
myBeamTransport.McEventCollectionName = "GEN_EVENT"
job += myBeamTransport
from ALFA_BeamAnalysis.ALFA_BeamAnalysisConf import ALFA_BeamAnalysisOtherProcess
myALFA_BeamAnalysisOtherProcess = ALFA_BeamAnalysisOtherProcess("ALFA_BeamAnalysisOtherProcess")
myALFA_BeamAnalysisOtherProcess.OutputFile="ALFA_BeamAnalysisOtherProcess.root"
myALFA_BeamAnalysisOtherProcess.FileOutput=True
myALFA_BeamAnalysisOtherProcess.BeamEnergy=7000.
myALFA_BeamAnalysisOtherProcess.Debug=True
# these cuts need to be the same as in ALFA_BeamTransport
myALFA_BeamAnalysisOtherProcess.EtaCut = 8.
myALFA_BeamAnalysisOtherProcess.XiCut = 0.2
job += myALFA_BeamAnalysisOtherProcess
#from ALFA_BeamAnalysis.ALFA_BeamAnalysisConf import ALFA_BeamAnalysis
#myALFA_BeamAnalysis = ALFA_BeamAnalysis("ALFA_BeamAnalysis")
#myALFA_BeamAnalysis.OutputFile="ALFA_BeamAnalysis.root"
#myALFA_BeamAnalysis.FileOutput=True
#myALFA_BeamAnalysis.BeamEnergy=7000.
#myALFA_BeamAnalysis.Debug=True
#job += myALFA_BeamAnalysis
#from HepMCAnalysis_i.HepMCAnalysis_iConfig import HepMCAnalysis_i
#myHepMCAnalysis = HepMCAnalysis_i("HepMCAnalysis_i", file = 'Pythia8_Process_elastic.root')
#myHepMCAnalysis.JetFinder=False
#myHepMCAnalysis.UserAnalysis=True
#job += myHepMCAnalysis
#---------------------------------------------------------------
# Pool Persistency
#---------------------------------------------------------------
#from AthenaPoolCnvSvc.WriteAthenaPool import AthenaPoolOutputStream
#Stream1 = AthenaPoolOutputStream( "Stream1" )
#Stream1.OutputFile = "elastic.pool.root"
#Stream1.ItemList += [ 'EventInfo#*', 'McEventCollection#*' ]
# end my own
#
# End of job options file
#
###############################################################
ForwardDetectors/ALFA/ALFA_BeamTransport/share/pythia8.common.py deleted 100755 → 0
+ 0
− 5
View file @ 8f35817c
# minimal changed settings
job.Pythia8_i.Commands += ["Main:writeAllSettings = on",
"Main:showAllStatistics = on",
"Random:setSeed = on",
"Random:seed = 0"]
ForwardDetectors/ALFA/ALFA_BeamTransport/share/pythia8.elastic.py deleted 100644 → 0
+ 0
− 29
View file @ 8f35817c
# -*- coding: utf-8 -*-
# elastic scattering with Pythia8
# Z decays into leptons and tau decays disabled
job.Pythia8_i.Commands += ["SoftQCD:elastic = on", # elastic scattering
"SigmaElastic:setOwn = on", # t-generation accoring to elastic formula
"SigmaTotal:setOwn = on", # use own total cross section
"SigmaTotal:sigmaTot = 100.", # actual value of sig_tot
"SigmaElastic:bSlope = 18.", # nuclear slope
"SigmaElastic:rho = 0.12.", # rho-parameter
"SigmaElastic:lambda = 0.72", # form factor parameter
"SigmaElastic:phaseConst = 0.577", # Coulomb phase constant (gamma_E)
"SigmaElastic:tAbsMin = 4e-5", # t_min for event generation
"Beams:allowMomentumSpread = on", # set up beam class (incoming protons)
"Beams:sigmaPxA = 0.010045", # set divergence in Px beam 1 (0.226 murad)
"Beams:sigmaPyA = 0.009485", # set divergence in Py beam 1
"Beams:sigmaPzA = 0.3955", # set energy smearing beam 1 in GeV (1.13E-4%)
"Beams:sigmaPxB = 0.010045", # set divergence in Px beam 2 (0.226 murad)
"Beams:sigmaPyB = 0.009485", # set divergence in Py beam 2
"Beams:sigmaPzB = 0.3955", # set energy smearing beam 2 in GeV (1.13E-4%)
"Beams:maxDevA = 100.0", # cut off smearing at 100 sigma
"Beams:maxDevB = 100.0", # cut off smearing at 100 sigma
"Beams:allowVertexSpread = on", # beam class vertex
"Beams:sigmaVertexX = 0.250 ", # spread in x
"Beams:sigmaVertexY = 0.237", # spread in y
"Beams:sigmaVertexZ = 0.0", # don't spread in z
"Beams:maxDevVertex = 10.0", # cut off spread at 10 sigma
"Random:setSeed=on", #from here added#########################################################
"Random:seed=4211",
]
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment