diff --git a/Simulation/FastSimulation/FastChainPileup/test/test_TrackingOverlay_JZ5W_wPU.sh b/Simulation/FastSimulation/FastChainPileup/test/test_TrackingOverlay_JZ5W_wPU.sh
new file mode 100755
index 0000000000000000000000000000000000000000..bb44ba6e05c558bbf32ac20eb1890b529cd58a29
--- /dev/null
+++ b/Simulation/FastSimulation/FastChainPileup/test/test_TrackingOverlay_JZ5W_wPU.sh
@@ -0,0 +1,41 @@
+#!/bin/sh
+#
+# art-description: Run the full reconstruction on HITS with pile-up for JZ5W samples (leading jet pt: 800-1300 GeV) for validation of tracking overlay
+# art-include: 21.0/Athena
+# art-include: master/Athena
+# art-type: grid
+
+export inputHitsFile=/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/TrackingOverlay/JZ5W_HITS
+export inputHighPtMinbiasHitsFile=/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/FastChainPileup/mc16_13TeV.361239.Pythia8EvtGen_A3NNPDF23LO_minbias_inelastic_high.simul.HITS.e4981_s3087_s3161
+export inputLowPtMinbiasHitsFile=/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/FastChainPileup/mc16_13TeV.361238.Pythia8EvtGen_A3NNPDF23LO_minbias_inelastic_low.simul.HITS.e4981_s3087_s3161
+export outputRDOFile=RDO.JZ5W_wPU.root
+export outputESDFile=ESD.JZ5W_wPU.root
+export outputAODFile=AOD.JZ5W_wPU.root
+export outputNTUP_PHYSVALFile=NTUP_PHYSVAL.JZ5W_wPU.root
+
+# Copied from r11627 (asetup Athena,21.0.97)
+# Running also with asetup 21.0,Athena,2020-02-28T2145
+Reco_tf.py --inputHITSFile $inputHitsFile/*root.? \
+           --inputHighPtMinbiasHitsFile $inputHighPtMinbiasHitsFile/*root.? \
+           --inputLowPtMinbiasHitsFile $inputLowPtMinbiasHitsFile/*.root.? \
+           --triggerConfig "RDOtoRDOTrigger=MCRECO:DBF:TRIGGERDBMC:2232,86,278" \
+           --steering "doRDO_TRIG" \
+           --autoConfiguration everything \
+           --conditionsTag "default:OFLCOND-MC16-SDR-25" \
+           --geometryVersion "default:ATLAS-R2-2016-01-00-01" \
+           --postExec "all:CfgMgr.MessageSvc().setError+=['HepMcParticleLink']" 'svcMgr.AthenaPoolCnvSvc.MaxFileSizes=["15000000000"]' \
+           --postInclude "default:PyJobTransforms/UseFrontier.py" \
+           --preExec "all:rec.Commissioning.set_Value_and_Lock(True);from AthenaCommon.BeamFlags import jobproperties;jobproperties.Beam.numberOfCollisions.set_Value_and_Lock(20.0);from LArROD.LArRODFlags import larRODFlags;larRODFlags.NumberOfCollisions.set_Value_and_Lock(20);larRODFlags.nSamples.set_Value_and_Lock(4);larRODFlags.doOFCPileupOptimization.set_Value_and_Lock(True);larRODFlags.firstSample.set_Value_and_Lock(0);larRODFlags.useHighestGainAutoCorr.set_Value_and_Lock(True); from LArDigitization.LArDigitizationFlags import jobproperties;jobproperties.LArDigitizationFlags.useEmecIwHighGain.set_Value_and_Lock(False)" 'HITtoRDO:userRunLumiOverride={"run":310000,"lb":61,"starttstamp":1550003600,"mu":60.500};' "from InDetPhysValMonitoring.InDetPhysValJobProperties import InDetPhysValFlags; InDetPhysValFlags.doValidateTightPrimaryTracks.set_Value_and_Lock(True);" \
+           --preInclude "HITtoRDO:Digitization/ForceUseOfPileUpTools.py,SimulationJobOptions/preInlcude.PileUpBunchTrainsMC16c_2017_Config1.py,RunDependentSimData/configLumi_user.py" \
+           --valid=True \
+           --validationFlags doEgamma,doInDet,doZee,doJet \
+           --outputRDOFile $outputRDOFile \
+           --outputESDFile $outputESDFile \
+           --outputAODFile $outputAODFile \
+           --outputNTUP_PHYSVALFile $outputNTUP_PHYSVALFile \
+           --jobNumber=2 \
+           --maxEvents=2000 \
+           --numberOfCavernBkg 0 \
+           --numberOfHighPtMinBias 0.2595392 \
+           --numberOfLowPtMinBias 99.2404608 \
+           --pileupFinalBunch 6
diff --git a/Simulation/FastSimulation/FastChainPileup/test/test_TrackingOverlay_JZ5W_woPU.sh b/Simulation/FastSimulation/FastChainPileup/test/test_TrackingOverlay_JZ5W_woPU.sh
new file mode 100755
index 0000000000000000000000000000000000000000..844acca7c7545f9b1739989c6e4fc1d078e1d513
--- /dev/null
+++ b/Simulation/FastSimulation/FastChainPileup/test/test_TrackingOverlay_JZ5W_woPU.sh
@@ -0,0 +1,32 @@
+#!/bin/sh
+#
+# art-description: Run the full reconstruction on HITS without pile-up JZ5W samples (leading jet pt: 800-1300 GeV) for validation of tracking overlay
+# art-include: 21.0/Athena
+# art-include: master/Athena
+# art-type: grid
+
+export inputHitsFile=/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/TrackingOverlay/JZ5W_HITS
+export outputRDOFile=RDO.JZ5W_woPU.root
+export outputESDFile=ESD.JZ5W_woPU.root
+export outputAODFile=AOD.JZ5W_woPU.root
+export outputNTUP_PHYSVALFile=NTUP_PHYSVAL.JZ5W_woPU.root
+
+# Copied from r11627 (asetup Athena,21.0.97)
+# Running also with asetup 21.0,Athena,2020-02-28T2145
+Reco_tf.py --inputHITSFile $inputHitsFile/*root.? \
+           --triggerConfig "RDOtoRDOTrigger=MCRECO:DBF:TRIGGERDBMC:2232,86,278" \
+           --steering "doRDO_TRIG" \
+           --autoConfiguration everything \
+           --conditionsTag "default:OFLCOND-MC16-SDR-25" \
+           --geometryVersion "default:ATLAS-R2-2016-01-00-01" \
+           --postExec "all:CfgMgr.MessageSvc().setError+=['HepMcParticleLink']" \
+           --postInclude "default:PyJobTransforms/UseFrontier.py" \
+           --preExec "all:rec.Commissioning.set_Value_and_Lock(True);from AthenaCommon.BeamFlags import jobproperties;jobproperties.Beam.numberOfCollisions.set_Value_and_Lock(20.0);from LArROD.LArRODFlags import larRODFlags;larRODFlags.NumberOfCollisions.set_Value_and_Lock(20);larRODFlags.nSamples.set_Value_and_Lock(4);larRODFlags.doOFCPileupOptimization.set_Value_and_Lock(True);larRODFlags.firstSample.set_Value_and_Lock(0);larRODFlags.useHighestGainAutoCorr.set_Value_and_Lock(True); from LArDigitization.LArDigitizationFlags import jobproperties;jobproperties.LArDigitizationFlags.useEmecIwHighGain.set_Value_and_Lock(False)" 'HITtoRDO:userRunLumiOverride={"run":310000,"lb":61,"starttstamp":1550003600,"mu":60.500};' "from InDetPhysValMonitoring.InDetPhysValJobProperties import InDetPhysValFlags; InDetPhysValFlags.doValidateTightPrimaryTracks.set_Value_and_Lock(True);" \
+           --valid=True \
+           --validationFlags doEgamma,doInDet,doZee,doJet \
+           --outputRDOFile $outputRDOFile \
+           --outputESDFile $outputESDFile \
+           --outputAODFile $outputAODFile \
+           --outputNTUP_PHYSVALFile $outputNTUP_PHYSVALFile \
+           --jobNumber=2 \
+           --maxEvents=2000