workflows 3, 4, 5

workflows/3/WF3.py

@@ -0,0 +1,33 @@
+#!/bin/env python
+"""script to do X"""
+from DIRAC.Core.Base.Script import parseCommandLine
+parseCommandLine()
+
+from ILCDIRAC.Interfaces.API.NewInterface.Applications import DelphesApp
+from ILCDIRAC.Interfaces.API.NewInterface.UserJob import UserJob
+from ILCDIRAC.Interfaces.API.DiracILC import DiracILC
+
+
+def run():
+  """Run The Job."""
+  job = UserJob()
+  job.setConfigPackage("fccConfig", 'key4hep-devel-2')
+
+  delphes = DelphesApp()
+  delphes.setVersion('key4hep-latest')
+  delphes.setExecutable("DelphesPythia8_EDM4HEP")
+  delphes.setDetectorCard('delphes_card_IDEA.tcl')
+  delphes.setOutputCard('edm4hep_output_config.tcl')
+  delphes.setPythia8CardFile('p8_ee_ggqq_ecm91.cmd')
+  delphes.setRandomSeed(12340)
+  delphes.setEnergy(91.2)
+  delphes.setNumberOfEvents(100)
+  delphes.setOutputFile('output.root')
+
+  job.append(delphes)
+  job.submit(DiracILC(), mode='local')
+
+
+
+if __name__ =="__main__":
+  run()

workflows/3/fccprodwf3

@@ -0,0 +1,53 @@
+[Marlin]
+#ApplicationAttributeName=Value
+
+[DDSim]
+#ApplicationAttributeName=Value
+
+[Overlay]
+#ApplicationAttributeName=Value
+
+[Whizard2]
+#ApplicationAttributeName=Value
+
+[KKMC]
+#ApplicationAttributeName=Value
+
+[Production Parameters]
+
+machine = ee
+prodGroup = several
+
+softwareVersion = key4hep-latest
+generatorApplication = delphes
+generatorVersion = key4hep-latest
+executableName = DelphesPythia8_EDM4HEP
+detectorCard = delphes_card_IDEA.tcl
+outputCard = edm4hep_output_config.tcl
+pythia8CardPath = p8_ee_ggqq_ecm91.cmd
+
+configVersion = key4hep-devel-2
+configPackage = fccConfig
+eventsPerJobs = 1000
+
+numberOfTasks = 1, 2
+
+campaign = winter2023
+energies = 91.2
+processes = qqbar
+detectorModel = IDEA
+datatype = delphes
+
+productionLogLevel = VERBOSE
+outputSE = CERN-DIP-4
+
+finalOutputSE = CERN-SRM
+MoveStatus = Stopped
+MoveGroupSize = 10
+
+#Productions to create: Gen, Split, Sim, Rec, RecOver
+ProdTypes = Gen
+move = False
+
+
+

workflows/3/p8_ee_ggqq_ecm91.cmd

@@ -0,0 +1,27 @@
+Random:setSeed = on
+Main:numberOfEvents = 3000         ! number of events to generate
+Main:timesAllowErrors = 5          ! how many aborts before run stops
+
+
+! 2) Settings related to output in init(), next() and stat().
+Init:showChangedSettings = on      ! list changed settings
+Init:showChangedParticleData = off ! list changed particle data
+Next:numberCount = 100            ! print message every n events
+Next:numberShowInfo = 1            ! print event information n times
+Next:numberShowProcess = 1         ! print process record n times
+Next:numberShowEvent = 0           ! print event record n times
+Stat:showPartonLevel = off
+
+! 3) Beam parameter settings. Values below agree with default ones.
+Beams:idA = 11                   ! first beam, e- = 11
+Beams:idB = -11                  ! second beam, e+ = -11
+
+! 4) Hard process : photon collisions at 91.188 GeV
+Beams:eCM = 91.188               ! CM energy of collision
+PhotonCollision:gmgm2qqbar = on
+PhotonCollision:gmgm2ccbar = on
+PhotonCollision:gmgm2bbbar = on
+
+PartonLevel:ISR = on               ! initial-state radiation
+PartonLevel:FSR = on               ! final-state radiation
+

workflows/4/WF4_p8_ee_ggqq.py

@@ -0,0 +1,30 @@
+#!/bin/env python
+"""script to do X"""
+from DIRAC.Core.Base.Script import parseCommandLine
+parseCommandLine()
+
+from ILCDIRAC.Interfaces.API.NewInterface.Applications import GaudiApp, DelphesApp
+from ILCDIRAC.Interfaces.API.NewInterface.UserJob import UserJob
+from ILCDIRAC.Interfaces.API.DiracILC import DiracILC
+
+def run():
+  """Run The Job."""
+  job = UserJob()
+  job.setConfigPackage("fccConfig", 'key4hep-devel-2')
+
+  ga = GaudiApp()
+  ga.setVersion('key4hep-latest')
+  ga.setExecutable("k4run")
+  ga.setSteeringFile('k4simdelphesalg_pythia.py')
+  ga.setPythia8CardFile('p8_ee_ggqq_ecm91.cmd') #p8_ee_ZH_ecm240.cmd p8_ee_WW_ecm240 p8_ee_ZZ_ecm240.cmd
+  ga.setExtraCLIArguments("--GenAlg.PythiaInterface.pythiacard pythia8card.cmd --k4SimDelphesAlg.DelphesCard delphes_card_IDEA.tcl --k4SimDelphesAlg.DelphesOutputSettings edm4hep_output_config.tcl")
+  ga.setEnergy(91.2)
+  ga.setNumberOfEvents(50)
+  ga.setOutputFile('output.root')
+
+  job.append(ga)
+  job.submit(DiracILC(), mode='local')
+
+
+if __name__ =="__main__":
+  run()

workflows/4/fccprodwf4

@@ -0,0 +1,52 @@
+[Marlin]
+#ApplicationAttributeName=Value
+
+[DDSim]
+#ApplicationAttributeName=Value
+
+[Overlay]
+#ApplicationAttributeName=Value
+
+[Whizard2]
+#ApplicationAttributeName=Value
+
+[KKMC]
+#ApplicationAttributeName=Value
+
+[Production Parameters]
+
+machine = ee
+prodGroup = several
+softwareVersion = key4hep-latest
+
+generatorApplication = gaudi
+generatorVersion = key4hep-latest
+executableName = k4run
+gaudiSteeringFile = k4simdelphesalg_pythia.py
+pythia8CardPath = p8_ee_ggqq_ecm91.cmd, p8_ee_WW_ecm240.cmd, p8_ee_ZH_ecm240.cmd, p8_ee_ZZ_ecm240.cmd
+
+randomSeedFlag = --SimG4Svc.seedValue
+
+configVersion = key4hep-devel-2
+configPackage = fccConfig
+eventsPerJobs = 1000, 2000, 3000, 1500
+
+numberOfTasks = 1, 1, 1, 1
+
+campaign = winter2023
+energies = 91.18, 240, 240, 240
+processes = qqbar, WW, ZZ, ZH
+detectorModel = idea
+datatype = delphes
+
+extraCLIArguments = --Pythia8.PythiaInterface.pythiacard pythia8card.cmd --k4SimDelphesAlg.DelphesCard delphes_card_IDEA.tcl --k4SimDelphesAlg.DelphesOutputSettings edm4hep_output_config.tcl
+
+productionLogLevel = VERBOSE
+outputSE = CERN-DIP-4
+
+finalOutputSE = CERN-SRM
+MoveStatus = Stopped
+MoveGroupSize = 10
+
+prodids =
+ProdTypes = Gen

workflows/4/k4simdelphesalg_pythia.py

@@ -0,0 +1,54 @@
+from Gaudi.Configuration import *
+from GaudiKernel import SystemOfUnits as units
+
+from Configurables import ApplicationMgr
+app = ApplicationMgr()
+app.EvtMax = 100
+app.EvtSel = "NONE"
+
+from Configurables import k4DataSvc
+podioevent = k4DataSvc("EventDataSvc")
+app.ExtSvc += [podioevent]
+
+
+from Configurables import PythiaInterface
+pythia8gentool = PythiaInterface()
+### Example of pythia configuration file to generate events
+# take from $K4GEN if defined, locally if not
+path_to_pythiafile = os.environ.get("K4GEN", "")
+pythiafilename = "Pythia_standard.cmd"
+pythiafile = os.path.join(path_to_pythiafile, pythiafilename)
+# Example of pythia configuration file to read LH event file
+#pythiafile="options/Pythia_LHEinput.cmd"
+pythia8gentool.pythiacard = pythiafile
+pythia8gentool.doEvtGenDecays = False
+pythia8gentool.printPythiaStatistics = True
+
+from Configurables import GenAlg
+gun = GenAlg()
+gun.SignalProvider = pythia8gentool
+gun.hepmc.Path = "hepmc"
+ApplicationMgr().TopAlg += [gun]
+
+
+from Configurables import HepMCToEDMConverter
+hepmc_converter = HepMCToEDMConverter()
+hepmc_converter.hepmc.Path="hepmc"
+hepmc_converter.GenParticles.Path = "GenParticles"
+ApplicationMgr().TopAlg += [hepmc_converter]
+
+from Configurables import k4SimDelphesAlg
+delphesalg = k4SimDelphesAlg()
+delphesalg.DelphesCard = "delphes_card_IDEA.tcl"
+delphesalg.DelphesOutputSettings = "edm4hep_output_config.tcl"
+delphesalg.GenParticles.Path = "GenParticles"
+delphesalg.OutputLevel = VERBOSE
+ApplicationMgr().TopAlg += [delphesalg]
+
+
+from Configurables import PodioOutput
+out = PodioOutput("out", filename = "output_k4SimDelphes_pythia.root")
+out.outputCommands = ["keep *"]
+ApplicationMgr().TopAlg += [out]
+
+

workflows/4/p8_ee_WW_ecm240.cmd

@@ -0,0 +1,22 @@
+Random:setSeed = on
+Main:numberOfEvents = 1000         ! number of events to generate
+Main:timesAllowErrors = 5          ! how many aborts before run stops
+
+! 2) Settings related to output in init(), next() and stat().
+Init:showChangedSettings = on      ! list changed settings
+Init:showChangedParticleData = off ! list changed particle data
+Next:numberCount = 100             ! print message every n events
+Next:numberShowInfo = 1            ! print event information n times
+Next:numberShowProcess = 1         ! print process record n times
+Next:numberShowEvent = 0           ! print event record n times
+
+Beams:idA = 11                   ! first beam, e+ = 11
+Beams:idB = -11                   ! second beam, e- = -11
+
+! 3) Hard process : WW at 240 GeV
+Beams:eCM = 240  ! CM energy of collision
+WeakDoubleBoson:ffbar2WW = on
+
+! 4) Settings for the event generation process in the Pythia8 library.
+PartonLevel:ISR = on               ! no initial-state radiation
+PartonLevel:FSR = on               ! no final-state radiation

workflows/4/p8_ee_ZH_ecm240.cmd

@@ -0,0 +1,25 @@
+Random:setSeed = on
+Main:numberOfEvents = 1000         ! number of events to generate
+Main:timesAllowErrors = 5          ! how many aborts before run stops
+
+! 2) Settings related to output in init(), next() and stat().
+Init:showChangedSettings = on      ! list changed settings
+Init:showChangedParticleData = off ! list changed particle data
+Next:numberCount = 100             ! print message every n events
+Next:numberShowInfo = 1            ! print event information n times
+Next:numberShowProcess = 1         ! print process record n times
+Next:numberShowEvent = 0           ! print event record n times
+
+Beams:idA = 11                   ! first beam, e+ = 11
+Beams:idB = -11                   ! second beam, e- = -11
+
+! 3) Hard process : ZH at 240 GeV
+Beams:eCM = 240.  ! CM energy of collision
+HiggsSM:ffbar2HZ = on
+
+! 4) Settings for the event generation process in the Pythia8 library.
+PartonLevel:ISR = on               ! initial-state radiation
+PartonLevel:FSR = on               ! final-state radiation
+
+! 5) Non-standard settings; exemplifies tuning possibilities.
+25:m0        = 125.0               ! Higgs mass

workflows/4/p8_ee_ZZ_ecm240.cmd

@@ -0,0 +1,23 @@
+Random:setSeed = on
+Main:numberOfEvents = 1000         ! number of events to generate
+Main:timesAllowErrors = 5          ! how many aborts before run stops
+
+! 2) Settings related to output in init(), next() and stat().
+Init:showChangedSettings = on      ! list changed settings
+Init:showChangedParticleData = off ! list changed particle data
+Next:numberCount = 100             ! print message every n events
+Next:numberShowInfo = 1            ! print event information n times
+Next:numberShowProcess = 1         ! print process record n times
+Next:numberShowEvent = 0           ! print event record n times
+
+Beams:idA = 11                   ! first beam, e+ = 11
+Beams:idB = -11                   ! second beam, e- = -11
+
+! 3) Hard process : ZZ at 240 GeV
+Beams:eCM = 240  ! CM energy of collision
+WeakDoubleBoson:ffbar2gmZgmZ = on
+
+! 4) Settings for the event generation process in the Pythia8 library.
+PartonLevel:ISR = on               ! no initial-state radiation
+PartonLevel:FSR = on               ! no final-state radiation
+

workflows/4/p8_ee_ggqq_ecm91.cmd

@@ -0,0 +1,27 @@
+Random:setSeed = on
+Main:numberOfEvents = 3000         ! number of events to generate
+Main:timesAllowErrors = 5          ! how many aborts before run stops
+
+
+! 2) Settings related to output in init(), next() and stat().
+Init:showChangedSettings = on      ! list changed settings
+Init:showChangedParticleData = off ! list changed particle data
+Next:numberCount = 100            ! print message every n events
+Next:numberShowInfo = 1            ! print event information n times
+Next:numberShowProcess = 1         ! print process record n times
+Next:numberShowEvent = 0           ! print event record n times
+Stat:showPartonLevel = off
+
+! 3) Beam parameter settings. Values below agree with default ones.
+Beams:idA = 11                   ! first beam, e- = 11
+Beams:idB = -11                  ! second beam, e+ = -11
+
+! 4) Hard process : photon collisions at 91.188 GeV
+Beams:eCM = 91.188               ! CM energy of collision
+PhotonCollision:gmgm2qqbar = on
+PhotonCollision:gmgm2ccbar = on
+PhotonCollision:gmgm2bbbar = on
+
+PartonLevel:ISR = on               ! initial-state radiation
+PartonLevel:FSR = on               ! final-state radiation
+

workflows/5/WF5.py

@@ -0,0 +1,30 @@
+#!/bin/env python
+"""script to do X"""
+from DIRAC.Core.Base.Script import parseCommandLine
+parseCommandLine()
+
+from ILCDIRAC.Interfaces.API.NewInterface.Applications import GaudiApp
+from ILCDIRAC.Interfaces.API.NewInterface.UserJob import UserJob
+from ILCDIRAC.Interfaces.API.DiracILC import DiracILC
+
+def run():
+  """Run The Job."""
+  job = UserJob()
+  job.setConfigPackage("fccConfig", 'key4hep-devel-2')
+
+  ga = GaudiApp()
+  ga.setVersion('key4hep-latest')
+  ga.setExecutable("k4run")
+  ga.setSteeringFile('geant_fullsim_fccee_lar_pgun.py')
+  ga.setRandomSeedFlag('--SimG4Svc.seedValue')
+  ga.setRandomSeed(100)
+  ga.setEnergy(10)
+  ga.setNumberOfEvents(25)
+  ga.setOutputFile('output.root')
+
+  job.append(ga)
+  job.submit(DiracILC(), mode='local')
+
+
+if __name__ =="__main__":
+  run()