Skip to content
Snippets Groups Projects
Commit b7cad5f8 authored by Frank Winklmeier's avatar Frank Winklmeier
Browse files

TrigT2BeamSpot: delete obsolete job option fragments and scripts

parent 04b20abc
6 merge requests!58791DataQualityConfigurations: Modify L1Calo config for web display,!46784MuonCondInterface: Enable thread-safety checking.,!46776Updated LArMonitoring config file for WD to match new files produced using MT,!45405updated ART test cron job,!42417Draft: DIRE and VINCIA Base Fragments for Pythia 8.3,!38724TriggerFlags: remove doLVL2/EF and other unused flags
...@@ -15,5 +15,3 @@ atlas_add_component( TrigT2BeamSpot ...@@ -15,5 +15,3 @@ atlas_add_component( TrigT2BeamSpot
# Install files from the package: # Install files from the package:
atlas_install_python_modules( python/*.py POST_BUILD_CMD ${ATLAS_FLAKE8} ) atlas_install_python_modules( python/*.py POST_BUILD_CMD ${ATLAS_FLAKE8} )
atlas_install_joboptions( share/*.py )
#----------------------------------------
# Return a file list given a directory
#----------------------------------------
import os,commands
def getDirectoryList(path):
cmd = 'ls %s/' % (path)
file = ["%s/%s" % (path,i) for i in commands.getoutput(cmd).split('\n')]
return file
def getMixedDirectoryList(path):
cmd = 'ls %s/*.root*' % (path)
file = ["%s" % (i) for i in commands.getoutput(cmd).split('\n')]
return file
def getCastorDirectoryList(path):
cmd = 'nsls %s/ ' % (path)
file = ["rfio:%s/%s" % (path,i) for i in commands.getoutput(cmd).split('\n')]
return files
#-----------------
# Input files
#-----------------
if ('Data' in dir()):
if Data == "shift":
PoolRDOInput = getDirectoryList('/u1/fizisist/ttbar_shift')
elif Data == "gashalo":
PoolRDOInput = getDirectoryList('/u1/fizisist/gashalo')
elif Data == "mixed":
PoolRDOInput = getMixedDirectoryList('/u1/iaracena/data/MixedEventsEnhancedb/RDO')
elif Data == "EnhancedBias":
PoolRDOInput = getMixedDirectoryList('/u1/data/user09.SethCaughron.EnhancedBias_10TeV.digit.RDO.e357_s462_d126')
elif Data == "J2":
PoolRDOInput = getDirectoryList('/u1/fizisist/J2')
elif Data == "J3":
PoolRDOInput = getDirectoryList('/u1/data/J3_RDO_rel12')
elif Data == "J1":
PoolRDOInput = getDirectoryList('/u1/fizisist/J1')
elif Data == "J4":
PoolRDOInput = getDirectoryList('/u1/fizisist/J4')
elif Data == "J0":
PoolRDOInput = getDirectoryList('/u1/fizisist/J0')
else:
Data = "shift"
PoolRDOInput = getDirectoryList('/u1/fizisist/ttbar_shift')
#-----------------
# Output files
#-----------------
PoolESDOutput = 'rfio:/castor/cern.ch/user/f/fizisist/beamspot/myTriggerTest2.mixed.ESD.pool.root'
PoolAODOutput = 'rfio:/castor/cern.ch/user/f/fizisist/beamspot/myTriggerTest2.mixed.AOD.pool.root'
if not ('EvtMax' in dir()):
EvtMax=1000
NEvents = "%d" % EvtMax
#-----------------
# Global Flags
#-----------------
from AthenaCommon.GlobalFlags import GlobalFlags
GlobalFlags.Print()
DetDescrVersion = "ATLAS-GEO-02-01-00"
GlobalFlags.DetDescrVersion = DetDescrVersion
#-----------------
# Global Setup
#-----------------
doTAG = False
doEgamma = False
doMuon = False
doJet = False
doTauRec = False
doObjMissingET = False
doMuGirl = False
doMoore = False
doMuonIDStandAlone = False
doMuonIDCombined = False
doCBNT = False
doWriteAOD = False
doWriteESD = False
doWriteTAG = False
doAOD = False
doESD = False
doTAG = False
doTruth = False
doAtlfast = False
#-----------------
# Trigger Setup
#-----------------
doTrigger = True
#-----------------------------------------------------------
include("RecExCommon/RecExCommon_flags.py")
#-----------------------------------------------------------
#-----------------
# InDet Flags
#-----------------
from InDetRecExample.InDetJobProperties import InDetFlags
from AthenaCommon.JobProperties import jobproperties
from InDetTrigRecExample.InDetTrigFlags import InDetTrigFlags
if Data == "gashalo":
InDetFlags.Enabled = True
InDetFlags.doBeamGas = True
InDetFlags.doBeamHalo = True
InDetFlags.doSiSPSeededTrackFinder = True
InDetFlags.doPrintConfigurables = True
InDetTrigFlags.doBeamGas = True
InDetTrigFlags.doBeamHalo = True
#-----------------
# Trigger Flags
#-----------------
from TriggerJobOpts.TriggerFlags import TriggerFlags
from TriggerMenuPython.BeamSpotSliceFlags import *
BeamSpotSliceFlags.setAll()
#-----------------
# BeamSpot Menu
#-----------------
def MyBeamSpotMenu():
if ('Seed' in dir()):
if Seed == "list":
TriggerFlags.BeamSpotSlice.signatures = [ 'beamspot_vertex_list' ]
elif Seed == "jet10":
TriggerFlags.BeamSpotSlice.signatures = [ 'beamspot_vertex_jet10' ]
elif Seed == "bjet18":
TriggerFlags.BeamSpotSlice.signatures = [ 'beamspot_vertex_bjet18' ]
elif Seed == "FSTracks":
TriggerFlags.BeamSpotSlice.signatures = [ 'beamspot_vertex_FSTracks' ]
elif Seed == "tau16i_loose":
TriggerFlags.BeamSpotSlice.signatures = [ 'beamspot_vertex_tau16i_loose' ]
elif Seed == "e10_medium":
TriggerFlags.BeamSpotSlice.signatures = [ 'beamspot_vertex_e10_medium' ]
elif Seed == "e10_medium_SiTrk":
TriggerFlags.BeamSpotSlice.signatures = [ 'beamspot_vertex_e10_medium_SiTrk' ]
elif Seed == "e10_medium_SiTrk_robust":
TriggerFlags.BeamSpotSlice.signatures = [ 'beamspot_vertex_e10_medium_SiTrk_robust' ]
elif Seed == "trk9i_calib":
TriggerFlags.BeamSpotSlice.signatures = [ 'beamspot_vertex_trk9i_calib' ]
else:
TriggerFlags.BeamSpotSlice.signatures = [ 'beamspot_vertex_FSTracks' ]
from TriggerMenuPython.GenerateMenu import GenerateMenu
GenerateMenu.overwriteSignaturesWith(MyBeamSpotMenu)
#-----------------
# Trigger Options
#-----------------
TriggerFlags.doCalo = True
TriggerFlags.doMuon = True
TriggerFlags.doID = True
TriggerFlags.doLVL2 = True
TriggerFlags.doLVL1 = True
TriggerFlags.doEF = False
TriggerFlags.readHLTconfigFromXML = False
TriggerFlags.readLVL1configFromXML = False
TriggerFlags.Slices_all_setOn()
TriggerFlags.doNtuple = True
TriggerFlags.enableMonitoring = [ 'Validation', 'Time', ]
TriggerFlags.triggerMenuSetup = 'default'
#TriggerFlags.triggerMenuSetup='Test_lumi1E31'
TriggerFlags.abortOnConfigurationError = True
TriggerFlags.doHypo = True
TriggerFlags.doFEX = True
#-----------------
# BeamSpot Config
#-----------------
from TrigT2BeamSpot.TrigT2BeamSpotConf import PESA__T2VertexBeamSpot
PESA__T2VertexBeamSpot.TrackAlgoID = 2
PESA__T2VertexBeamSpot.HistoTrackAlgs = True
PESA__T2VertexBeamSpot.SaveAllVertices = False
PESA__T2VertexBeamSpot.TotalNTrackMin = 3
PESA__T2VertexBeamSpot.TrackMinPt = 0.7
PESA__T2VertexBeamSpot.TrackSeedPt = 1.0
PESA__T2VertexBeamSpot.TrackMaxZ0 = 200.0
PESA__T2VertexBeamSpot.TrackMaxD0 = 10.0
PESA__T2VertexBeamSpot.TrackMaxZ0err = 1.0
PESA__T2VertexBeamSpot.TrackMaxD0err = 1.0
PESA__T2VertexBeamSpot.TrackMinQual = 0.0
PESA__T2VertexBeamSpot.TrackMaxQual = 10.0
PESA__T2VertexBeamSpot.TrackMinSiHits = 4
PESA__T2VertexBeamSpot.TrackMinTRTHits = 1
PESA__T2VertexBeamSpot.VertexMinNTrk = 3
PESA__T2VertexBeamSpot.VertexMaxNTrk = 30
PESA__T2VertexBeamSpot.VertexMaxXerr = 0.2
PESA__T2VertexBeamSpot.VertexMaxYerr = 0.2
PESA__T2VertexBeamSpot.VertexMaxQual = 1000.0
#-----------------------------------------------------------
include("RecExCommon/RecExCommon_topOptions.py")
#-----------------------------------------------------------
from GaudiSvc.GaudiSvcConf import THistSvc
ServiceMgr += THistSvc()
ServiceMgr.THistSvc.Output = ["SHIFT DATAFILE='shift-monitoring.root' OPT='RECREATE'",
"EXPERT DATAFILE='expert-monitoring." + SiTrackOpt + "." + Data + "." + Seed + "." + NEvents + "evt." + TrackType + ".root' OPT='RECREATE'",
"RUNSTAT DATAFILE='runstat-monitoring.root' OPT='RECREATE'",
"DEBUG DATAFILE='debug-monitoring.root' OPT='RECREATE'"]
jobproperties.PerfMonFlags.doMonitoring = False
Service ("StoreGateSvc" ).ActivateHistory = True
print AlgSequence
print ServiceMgr
if not ('OutputLevel' in dir()):
OutputLevel = INFO
AlgSequence().TrigSteer_L2.OutputLevel = OutputLevel
for algo in AlgSequence().TrigSteer_L2.getChildren():
algo.OutputLevel = OutputLevel
MessageSvc.debugLimit = 1000000
MessageSvc.OutputLevel = OutputLevel
MessageSvc.Format = "% F%48W%S%7W%R%T %0W%M"
ServiceMgr.StoreGateSvc.Dump = False
#----------------------------------------
# Return a file list given a directory
#----------------------------------------
import os,commands
def getDirectoryList(path):
cmd = 'ls %s/' % (path)
file = ["%s/%s" % (path,i) for i in commands.getoutput(cmd).split('\n')]
return file
def getMixedDirectoryList(path):
cmd = 'ls %s/*.root*' % (path)
file = ["%s" % (i) for i in commands.getoutput(cmd).split('\n')]
return file
def getCastorDirectoryList(path):
cmd = 'nsls %s/ ' % (path)
file = ["rfio:%s/%s" % (path,i) for i in commands.getoutput(cmd).split('\n')]
return files
#-----------------
# Input/Output
#-----------------
if ('Data' in dir()):
if Data == "shift":
PoolRDOInput = getDirectoryList('/u1/fizisist/ttbar_shift')
elif Data == "gashalo":
PoolRDOInput = getDirectoryList('/u1/fizisist/gashalo')
elif Data == "mixed":
PoolRDOInput = getMixedDirectoryList('/u1/iaracena/data/MixedEventsEnhancedb/RDO')
elif Data == "J2":
PoolRDOInput = getDirectoryList('/u1/fizisist/J2')
elif Data == "J3":
PoolRDOInput = getDirectoryList('/u1/fizisist/J3_2')
elif Data == "J1":
PoolRDOInput = getDirectoryList('/u1/fizisist/J1')
elif Data == "J4":
PoolRDOInput = getDirectoryList('/u1/fizisist/J4')
elif Data == "J0":
PoolRDOInput = getDirectoryList('/u1/fizisist/J0')
else:
Data = "shift"
PoolRDOInput = getDirectoryList('/u1/fizisist/ttbar_shift')
#PoolRDOInput=["/afs/cern.ch/atlas/offline/data/testfile/calib1_csc11.005200.T1_McAtNlo_Jimmy.digit.RDO.v12000301_tid003138._00016_extract_10evt.pool.root"]
#PoolRDOInput = ["rfio:/castor/cern.ch/user/f/fizisist/beamspot/misal1_csc11.005013.J4_pythia_jetjet.digit.RDO.v12003103_tid004016._02277.pool.root.1"]
#PoolRDOInput = getCastorDirectoryList('/castor/cern.ch/user/f/fizisist/beamspot/ttbarShift')
#PoolRDOInput = getCastorDirectoryList('/castor/cern.ch/grid/atlas/caf/mc12/mc12.005012.J3_pythia_jetjet.digit.RDO.v12000003')
#PoolRDOInput = getCastorDirectoryList('/castor/cern.ch/grid/atlas/dq2/misal1_csc11/misal1_csc11.005013.J4_pythia_jetjet.digit.RDO.v12003103_tid004016')
#readRDOBS=["/u1/iaracena/data/14.0.X/Mixed_lumi1E31/daq.csc13.0000000.Single.Stream.LB0000.Athena._0001.data"]
#PoolRDOInput = getDirectoryList('/u1/fizisist/ttbar_shift')
#PoolRDOInput = getMixedDirectoryList('/u1/iaracena/data/MixedEventsEnhancedb')
#PoolRDOInput = getDirectoryList('/u1/fizisist/J2')
#PoolRDOInput = getDirectoryList('/u1/fizisist/J3')
PoolESDOutput = 'rfio:/castor/cern.ch/user/f/fizisist/beamspot/myTriggerTest2.mixed.ESD.pool.root'
PoolAODOutput = 'rfio:/castor/cern.ch/user/f/fizisist/beamspot/myTriggerTest2.mixed.AOD.pool.root'
if not ('EvtMax' in dir()):
EvtMax=1000
NEvents = "%d" % EvtMax
#-----------------
# Global Flags
#-----------------
from AthenaCommon.GlobalFlags import GlobalFlags
GlobalFlags.Print()
#GlobalFlags.DetDescrVersion ="ATLAS-CSC-01-02-00"
#DetDescrVersion = "ATLAS-CSC-01-02-00"
#-----------------
# Global Setup
#-----------------
doTAG = False
doEgamma = False
doMuon = False
doJet = False
doTauRec = False
doObjMissingET = False
doMuGirl = False
doMoore = False
doMuonIDStandAlone = False
doMuonIDCombined = False
doCBNT = False
doWriteAOD = False
doWriteESD = False
doWriteTAG = False
doAOD = False
doESD = False
doTAG = False
doTruth = False
doAtlfast = False
#-----------------
# Trigger Setup
#-----------------
doTrigger = True
#-----------------------------------------------------------
include("RecExCommon/RecExCommon_flags.py")
#-----------------------------------------------------------
#-----------------
# InDet Flags
#-----------------
from InDetRecExample.InDetJobProperties import InDetFlags
from AthenaCommon.JobProperties import jobproperties
from InDetTrigRecExample.InDetTrigFlags import InDetTrigFlags
if Data == "gashalo":
InDetFlags.Enabled = True
InDetFlags.doBeamGas = True
InDetFlags.doBeamHalo = True
InDetFlags.doSiSPSeededTrackFinder = True
InDetFlags.doPrintConfigurables = True
InDetTrigFlags.doBeamGas = True
InDetTrigFlags.doBeamHalo = True
#-----------------
# Trigger Flags
#-----------------
TriggerFlags.doCalo = True
TriggerFlags.doMuon = True
TriggerFlags.doID = True
TriggerFlags.doLVL2 = True
TriggerFlags.doLVL1 = True
TriggerFlags.doEF = False
TriggerFlags.readHLTconfigFromXML = False
TriggerFlags.readLVL1configFromXML = False
TriggerFlags.Slices_all_setOn()
TriggerFlags.doNtuple = True
TriggerFlags.enableMonitoring = [ 'Validation', 'Time', ]
TriggerFlags.triggerMenuSetup = 'default'
TriggerFlags.abortOnConfigurationError = True
TriggerFlags.doHypo = True
TriggerFlags.doFEX = True
#-----------------
# SiTrack Setup
#-----------------
#from TrigSiTrack.TrigSiTrack_Config import *
#if ('SiTrackOpt' in dir()):
# if SiTrackOpt == "loose.v2":
# TrigSiTrackFlags.LoosenCuts = True
# else:
# TrigSiTrackFlags.LoosenCuts = False
#else:
SiTrackOpt = "defaults"
# TrigSiTrackFlags.LoosenCuts = False
#-----------------
# Beamspot Setup
#-----------------
from TriggerMenuPython.BeamSpotSliceFlags import *
BeamSpotSliceFlags.unsetAll()
TriggerFlags.BeamSpotSlice.signatures = [ 'vertex_bjet18' ]
from TrigT2BeamSpot.TrigT2BeamSpotConf import PESA__T2VertexBeamSpot
# Track Algorithm: SiTrack = 1, IdScan = 2
if ('TrackType' in dir()):
if TrackType == "SiTrack":
PESA__T2VertexBeamSpot.TrackAlgoID = 1
else:
PESA__T2VertexBeamSpot.TrackAlgoID = 2
else:
TrackType = "IdScan"
PESA__T2VertexBeamSpot.TrackAlgoID = 2
PESA__T2VertexBeamSpot.UseBeamConditionService = True
PESA__T2VertexBeamSpot.DoPerLBHistograms = True
PESA__T2VertexBeamSpot.LBHistoryLength = 10
PESA__T2VertexBeamSpot.MonitoringHistoPath = "/EXPERT/T2VertexBeamSpot_Fex/"
PESA__T2VertexBeamSpot.ExpressHistoPath = "/EXPERT/T2VertexBeamSpot_Fex/"
PESA__T2VertexBeamSpot.HistoTrackAlgs = False
PESA__T2VertexBeamSpot.SaveAllVertices = False
PESA__T2VertexBeamSpot.TotalNTrackMin = 6
PESA__T2VertexBeamSpot.TrackMinPt = 1.0
PESA__T2VertexBeamSpot.TrackSeedPt = 4.0
PESA__T2VertexBeamSpot.TrackMaxZ0 = 200.0
PESA__T2VertexBeamSpot.TrackMaxD0 = 10.0
PESA__T2VertexBeamSpot.TrackMaxZ0err = 1.0
PESA__T2VertexBeamSpot.TrackMaxD0err = 1.0
PESA__T2VertexBeamSpot.TrackMinQual = 0.0
PESA__T2VertexBeamSpot.TrackMaxQual = 10.0
PESA__T2VertexBeamSpot.TrackMinSiHits = 5
PESA__T2VertexBeamSpot.TrackMinTRTHits = -1
PESA__T2VertexBeamSpot.VertexMinNTrk = 6
PESA__T2VertexBeamSpot.VertexMaxNTrk = 30
PESA__T2VertexBeamSpot.VertexMaxXerr = 0.1
PESA__T2VertexBeamSpot.VertexMaxYerr = 0.1
PESA__T2VertexBeamSpot.VertexMaxQual = 30.0
#-----------------------------------------------------------
include("RecExCommon/RecExCommon_topOptions.py")
#-----------------------------------------------------------
from GaudiSvc.GaudiSvcConf import THistSvc
ServiceMgr += THistSvc()
ServiceMgr.THistSvc.Output = ["SHIFT DATAFILE='shift-monitoring.root' OPT='RECREATE'",
"EXPERT DATAFILE='/u1/fizisist/beamspot/rootHistos/expert-monitoring.b18." + SiTrackOpt + "." + Data + "." + NEvents + "evt." + TrackType + ".root' OPT='RECREATE'",
"RUNSTAT DATAFILE='runstat-monitoring.root' OPT='RECREATE'",
"DEBUG DATAFILE='debug-monitoring.root' OPT='RECREATE'"]
jobproperties.PerfMonFlags.doMonitoring = False
Service ("StoreGateSvc" ).ActivateHistory = True
print AlgSequence
print ServiceMgr
if not ('OutputLevel' in dir()):
OutputLevel = INFO
AlgSequence().TrigSteer_L2.OutputLevel = OutputLevel
for algo in AlgSequence().TrigSteer_L2.getChildren():
algo.OutputLevel = INFO
MessageSvc.debugLimit = 1000000
MessageSvc.OutputLevel = OutputLevel
MessageSvc.Format = "% F%48W%S%7W%R%T %0W%M"
ServiceMgr.StoreGateSvc.Dump = True
#!/bin/bash
echo
echo "We will now be running several jobs for testing the beam spot code"
echo
homeDir=/afs/cern.ch/user/f/fizisist
baseDir=/u1/fizisist/beamspot
logDir=${baseDir}
jobOptions=${baseDir}/myTrigger_jobOptions.py
dataList="gashalo J1 J3 mixed shift" # Data
trackList="SiTrack IdScan" # TrackType
optsList="defaults.v8" # SiTrackOpt
#source ${homeDir}/work/14.0.0/cmt/setup.sh -tag=14.0.0,32
cd ${baseDir}
for datatype in ${dataList}; do
for trktype in ${trackList}; do
for cutstype in ${optsList}; do
outDir=${baseDir}/${datatype}/${trktype}/${cutstype}
if [ ! -d ${outDir} ]; then
mkdir ${outDir}
fi
cd ${outDir}
echo "*** ${trktype} (SiTrack=${cutstype}) on ${datatype} events"
#echo "athena.py -c EvtMax=$1;Data='${datatype}';TrackType='${trktype}';SiTrackOpt='${cutstype}'; ${jobOptions} >& ${logDir}/beamspot.${datatype}.$1evt.${trktype}.${cutstype}.log & "
athena.py -c "EvtMax=$1;Data='${datatype}';TrackType='${trktype}';SiTrackOpt='${cutstype}';" ${jobOptions} >& ${logDir}/beamspot.b18.${datatype}.$1evt.${trktype}.${cutstype}.log &
done
done
done
echo
cd ${baseDir}
#!/bin/bash
echo
echo "Running single job"
echo
##--------------
## Seed choices
##--------------
# beamspot_vertex_list
# beamspot_vertex_jet10
# beamspot_vertex_bjet18
# beamspot_vertex_FSTracks
# beamspot_vertex_tau16i_loose
# beamspot_vertex_e10_medium
# beamspot_vertex_e10_medium_SiTrk
# beamspot_vertex_e10_medium_SiTrk_robust
# beamspot_vertex_trk9i_calib
mydate=`date '+%Y-%m-%d-%HH-%MM-%SS'`
evtMax=1000
datatype=EnhancedBias
L2Seed=beamspot_vertex_FSTracks
trktype=IdScan
cutstype=${mydate}
jobOptions=myNewTrigger_jobOptions.py
logDir=/u1/fizisist/beamspot
Output=DEBUG
athena.py -c "EvtMax=${evtMax};OutputLevel=${Output};Data='${datatype}';Seed='${L2Seed}';TrackType='${trktype}';SiTrackOpt='${cutstype}';" ${jobOptions} >& ${logDir}/beamspot.${datatype}.${evtMax}evt.${trktype}.${L2Seed}.${cutstype}.${mydate}.log &
#mydate=`date '+%Y-%m-%d-%HH-%MM-%SS'`
#athenaMT.py -M -W -c isOnline=True -l INFO -f Raw.data -n 10000 testHLT_standalone.py >& /u1/fizisist/beamspot/testMT.raw.withOH.${mydate}.log &
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment