From d66cbcb7602b58f1c3e42618367fb74476dfda7e Mon Sep 17 00:00:00 2001 From: Joosep Pata Date: Fri, 10 May 2024 16:18:18 +0300 Subject: [PATCH 1/8] wip for MLPF sim --- CLDConfig/CLDReconstruction.py | 8 ++-- CLDConfig/p8_ee_tt_ecm365.cmd | 33 +++++++++++++++ CLDConfig/pythia.py | 73 ++++++++++++++++++++++++++++++++++ CLDConfig/run_sim.sh | 37 +++++++++++++++++ 4 files changed, 147 insertions(+), 4 deletions(-) create mode 100644 CLDConfig/p8_ee_tt_ecm365.cmd create mode 100644 CLDConfig/pythia.py create mode 100644 CLDConfig/run_sim.sh diff --git a/CLDConfig/CLDReconstruction.py b/CLDConfig/CLDReconstruction.py index b82f84a..87631a5 100644 --- a/CLDConfig/CLDReconstruction.py +++ b/CLDConfig/CLDReconstruction.py @@ -17,7 +17,7 @@ # limitations under the License. # import os -from Gaudi.Configuration import INFO, WARNING, DEBUG +from Gaudi.Configuration import INFO, WARNING, INFO from Configurables import k4DataSvc, MarlinProcessorWrapper from k4MarlinWrapper.inputReader import create_reader, attach_edm4hep2lcio_conversion @@ -468,7 +468,7 @@ JetClusteringAndRefiner.OutputLevel = WARNING JetClusteringAndRefiner.ProcessorType = "LcfiplusProcessor" JetClusteringAndRefiner.Parameters = { - "Algorithms": ["JetClustering", "JetVertexRefiner"], + "Algorithms": ["JetClustering"], "JetClustering.AlphaParameter": ["1.0"], "JetClustering.BetaParameter": ["1.0"], "JetClustering.GammaParameter": ["1.0"], @@ -1097,7 +1097,7 @@ "DropCollectionNames": [], "DropCollectionTypes": ["MCParticle", "LCRelation", "SimCalorimeterHit", "CalorimeterHit", "SimTrackerHit", "TrackerHit", "TrackerHitPlane", "Track", "ReconstructedParticle", "LCFloatVec"], "FullSubsetCollections": ["EfficientMCParticles", "InefficientMCParticles", "MCPhysicsParticles"], - "KeepCollectionNames": ["MCParticlesSkimmed", "MCPhysicsParticles", "RecoMCTruthLink", "SiTracks", "SiTracks_Refitted", "PandoraClusters", "PandoraPFOs", "SelectedPandoraPFOs", "LooseSelectedPandoraPFOs", "TightSelectedPandoraPFOs", "RefinedVertexJets", "RefinedVertexJets_rel", "RefinedVertexJets_vtx", "RefinedVertexJets_vtx_RP", "BuildUpVertices", "BuildUpVertices_res", "BuildUpVertices_RP", "BuildUpVertices_res_RP", "BuildUpVertices_V0", "BuildUpVertices_V0_res", "BuildUpVertices_V0_RP", "BuildUpVertices_V0_res_RP", "PrimaryVertices", "PrimaryVertices_res", "PrimaryVertices_RP", "PrimaryVertices_res_RP", "RefinedVertices", "RefinedVertices_RP"], + "KeepCollectionNames": ["MCParticlesSkimmed", "MCPhysicsParticles", "RecoMCTruthLink", "SiTracks", "SiTracks_Refitted", "PandoraClusters", "PandoraPFOs", "SelectedPandoraPFOs", "LooseSelectedPandoraPFOs", "TightSelectedPandoraPFOs", "BuildUpVertices", "BuildUpVertices_res", "BuildUpVertices_RP", "BuildUpVertices_res_RP", "BuildUpVertices_V0", "BuildUpVertices_V0_res", "BuildUpVertices_V0_RP", "BuildUpVertices_V0_res_RP", "PrimaryVertices", "PrimaryVertices_res", "PrimaryVertices_RP", "PrimaryVertices_res_RP", "RefinedVertices", "RefinedVertices_RP"], "LCIOOutputFile": [f"{reco_args.outputBasename}_DST.slcio"], "LCIOWriteMode": ["WRITE_NEW"] } @@ -1111,7 +1111,7 @@ lcioConvTool.collNameMapping = { "MCParticle": "MCParticles" } - lcioConvTool.OutputLevel = DEBUG + lcioConvTool.OutputLevel = INFO # attach to the last non output processor EventNumber.Lcio2EDM4hepTool = lcioConvTool diff --git a/CLDConfig/p8_ee_tt_ecm365.cmd b/CLDConfig/p8_ee_tt_ecm365.cmd new file mode 100644 index 0000000..448739c --- /dev/null +++ b/CLDConfig/p8_ee_tt_ecm365.cmd @@ -0,0 +1,33 @@ +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 + +Beams:allowMomentumSpread = off + +! Vertex smearing : +Beams:allowVertexSpread = on +Beams:sigmaVertexX = 27.3e-3 +Beams:sigmaVertexY = 48.8E-6 +Beams:sigmaVertexZ = 1.33 +Beams:sigmaTime = 1.94 ! 6.46 ps + + +! 3) Hard process : tt at 365 GeV +Beams:eCM = 365 ! CM energy of collision +Top:ffbar2ttbar(s:gmZ) = 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 + diff --git a/CLDConfig/pythia.py b/CLDConfig/pythia.py new file mode 100644 index 0000000..b54e4d7 --- /dev/null +++ b/CLDConfig/pythia.py @@ -0,0 +1,73 @@ +""" +Pythia8, integrated in the FCCSW framework. + +Generates according to a pythia .cmd file and saves them in fcc edm format. + +""" + +import os +from GaudiKernel import SystemOfUnits as units +from Gaudi.Configuration import * + +from Configurables import ApplicationMgr +ApplicationMgr().EvtSel = 'NONE' +ApplicationMgr().EvtMax = 1000 +ApplicationMgr().OutputLevel = INFO +ApplicationMgr().ExtSvc +=["RndmGenSvc"] + +#### Data service +from Configurables import k4DataSvc +podioevent = k4DataSvc("EventDataSvc") +ApplicationMgr().ExtSvc += [podioevent] + +from Configurables import GaussSmearVertex +smeartool = GaussSmearVertex() +smeartool.xVertexSigma = 0.5*units.mm +smeartool.yVertexSigma = 0.5*units.mm +smeartool.zVertexSigma = 40.0*units.mm +smeartool.tVertexSigma = 180.0*units.picosecond + +from Configurables import PythiaInterface +pythia8gentool = PythiaInterface() +### Example of pythia configuration file to generate events +# take from $K4GEN if defined, locally if not +pythia8gentool.pythiacard = "p8_ee_Z_Ztautau_ecm125.cmd" +pythia8gentool.doEvtGenDecays = False +pythia8gentool.printPythiaStatistics = True + +from Configurables import GenAlg +pythia8gen = GenAlg("Pythia8") +pythia8gen.SignalProvider = pythia8gentool +pythia8gen.VertexSmearingTool = smeartool +pythia8gen.hepmc.Path = "hepmc" +ApplicationMgr().TopAlg += [pythia8gen] + +### Reads an HepMC::GenEvent from the data service and writes a collection of EDM Particles +from Configurables import HepMCToEDMConverter +hepmc_converter = HepMCToEDMConverter() +ApplicationMgr().TopAlg += [hepmc_converter] + +### Filters generated particles +# accept is a list of particle statuses that should be accepted +from Configurables import GenParticleFilter +genfilter = GenParticleFilter("StableParticles") +genfilter.accept = [1] +genfilter.GenParticles.Path = "GenParticles" +genfilter.GenParticlesFiltered.Path = "GenParticlesStable" +ApplicationMgr().TopAlg += [genfilter] + +from Configurables import EDMToHepMCConverter +edm_converter = EDMToHepMCConverter("BackConverter") +edm_converter.GenParticles.Path="GenParticleStable" + +from Configurables import HepMCFileWriter +dumper = HepMCFileWriter("Dumper") +dumper.hepmc.Path = "hepmc" +ApplicationMgr().TopAlg += [dumper] + +#from Configurables import PodioOutput +#out = PodioOutput("out") +#out.filename = "out_pythia.root" +#out.outputCommands = ["keep *"] + + diff --git a/CLDConfig/run_sim.sh b/CLDConfig/run_sim.sh new file mode 100644 index 0000000..baebb66 --- /dev/null +++ b/CLDConfig/run_sim.sh @@ -0,0 +1,37 @@ +#!/bin/bash +#SBATCH -p main +#SBATCH --mem-per-cpu=6G +#SBATCH --cpus-per-task=1 +#SBATCH -o logs/slurm-%x-%j-%N.out +set -e +set -x + +env +df -h + +NEV=100 +NUM=$1 #random seed +SAMPLE=$2 #main card + +#Change these as needed +OUTDIR=/local/joosep/cld_edm4hep/2024_05/ +SIMDIR=/home/joosep/key4hep-sim/cld/CLDConfig/CLDConfig/ +WORKDIR=/scratch/local/$USER/${SAMPLE}_${SLURM_JOB_ID} +FULLOUTDIR=${OUTDIR}/${SAMPLE} + +mkdir -p $FULLOUTDIR + +mkdir -p $WORKDIR +cd $WORKDIR + +cp $SIMDIR/${SAMPLE}.cmd card.cmd +cp $SIMDIR/pythia.py ./ +cp $SIMDIR/cld_steer.py ./ +cp -R $SIMDIR/PandoraSettingsCLD ./ +cp -R $SIMDIR/CLDReconstruction.py ./ + +source /cvmfs/sw-nightlies.hsf.org/key4hep/setup.sh + +k4run pythia.py -n 100 --Dumper.Filename out.hepmc --Pythia8.PythiaInterface.pythiacard card.cmd +ddsim -I out.hepmc -N -1 -O out_SIM.root --compactFile $K4GEO/FCCee/CLD/compact/CLD_o2_v05/CLD_o2_v05.xml --steeringFile cld_steer.py +k4run CLDReconstruction.py --inputFiles out_SIM.root --outputBasename out_RECO --num-events -1 From 5de045cd1af7a89834d9c2976d040b773a85d27a Mon Sep 17 00:00:00 2001 From: Joosep Pata Date: Fri, 10 May 2024 16:23:32 +0300 Subject: [PATCH 2/8] up --- CLDConfig/run_sim.sh | 13 +++++++++++-- 1 file changed, 11 insertions(+), 2 deletions(-) mode change 100644 => 100755 CLDConfig/run_sim.sh diff --git a/CLDConfig/run_sim.sh b/CLDConfig/run_sim.sh old mode 100644 new mode 100755 index baebb66..a66dcf2 --- a/CLDConfig/run_sim.sh +++ b/CLDConfig/run_sim.sh @@ -15,7 +15,7 @@ SAMPLE=$2 #main card #Change these as needed OUTDIR=/local/joosep/cld_edm4hep/2024_05/ -SIMDIR=/home/joosep/key4hep-sim/cld/CLDConfig/CLDConfig/ +SIMDIR=/home/joosep/key4hep-sim/cld/CLDConfig/CLDConfig WORKDIR=/scratch/local/$USER/${SAMPLE}_${SLURM_JOB_ID} FULLOUTDIR=${OUTDIR}/${SAMPLE} @@ -30,8 +30,17 @@ cp $SIMDIR/cld_steer.py ./ cp -R $SIMDIR/PandoraSettingsCLD ./ cp -R $SIMDIR/CLDReconstruction.py ./ +echo "Random:seed=${NUM}" >> card.cmd +cat card.cmd + source /cvmfs/sw-nightlies.hsf.org/key4hep/setup.sh -k4run pythia.py -n 100 --Dumper.Filename out.hepmc --Pythia8.PythiaInterface.pythiacard card.cmd +k4run pythia.py -n $NEV --Dumper.Filename out.hepmc --Pythia8.PythiaInterface.pythiacard card.cmd ddsim -I out.hepmc -N -1 -O out_SIM.root --compactFile $K4GEO/FCCee/CLD/compact/CLD_o2_v05/CLD_o2_v05.xml --steeringFile cld_steer.py k4run CLDReconstruction.py --inputFiles out_SIM.root --outputBasename out_RECO --num-events -1 + +#Copy the outputs +cp out_RECO_edm4hep.root $FULLOUTDIR/reco_${SAMPLE}_${NUM}.root +cp out.hepmc $FULLOUTDIR/sim_${SAMPLE}_${NUM}.hepmc + +rm -Rf $WORKDIR From 17a274a3c3c4673056926e5b0168179c6c928c62 Mon Sep 17 00:00:00 2001 From: Joosep Pata Date: Fri, 10 May 2024 16:28:25 +0300 Subject: [PATCH 3/8] up --- CLDConfig/run_sim.sh | 46 -------------------------------------------- 1 file changed, 46 deletions(-) delete mode 100755 CLDConfig/run_sim.sh diff --git a/CLDConfig/run_sim.sh b/CLDConfig/run_sim.sh deleted file mode 100755 index a66dcf2..0000000 --- a/CLDConfig/run_sim.sh +++ /dev/null @@ -1,46 +0,0 @@ -#!/bin/bash -#SBATCH -p main -#SBATCH --mem-per-cpu=6G -#SBATCH --cpus-per-task=1 -#SBATCH -o logs/slurm-%x-%j-%N.out -set -e -set -x - -env -df -h - -NEV=100 -NUM=$1 #random seed -SAMPLE=$2 #main card - -#Change these as needed -OUTDIR=/local/joosep/cld_edm4hep/2024_05/ -SIMDIR=/home/joosep/key4hep-sim/cld/CLDConfig/CLDConfig -WORKDIR=/scratch/local/$USER/${SAMPLE}_${SLURM_JOB_ID} -FULLOUTDIR=${OUTDIR}/${SAMPLE} - -mkdir -p $FULLOUTDIR - -mkdir -p $WORKDIR -cd $WORKDIR - -cp $SIMDIR/${SAMPLE}.cmd card.cmd -cp $SIMDIR/pythia.py ./ -cp $SIMDIR/cld_steer.py ./ -cp -R $SIMDIR/PandoraSettingsCLD ./ -cp -R $SIMDIR/CLDReconstruction.py ./ - -echo "Random:seed=${NUM}" >> card.cmd -cat card.cmd - -source /cvmfs/sw-nightlies.hsf.org/key4hep/setup.sh - -k4run pythia.py -n $NEV --Dumper.Filename out.hepmc --Pythia8.PythiaInterface.pythiacard card.cmd -ddsim -I out.hepmc -N -1 -O out_SIM.root --compactFile $K4GEO/FCCee/CLD/compact/CLD_o2_v05/CLD_o2_v05.xml --steeringFile cld_steer.py -k4run CLDReconstruction.py --inputFiles out_SIM.root --outputBasename out_RECO --num-events -1 - -#Copy the outputs -cp out_RECO_edm4hep.root $FULLOUTDIR/reco_${SAMPLE}_${NUM}.root -cp out.hepmc $FULLOUTDIR/sim_${SAMPLE}_${NUM}.hepmc - -rm -Rf $WORKDIR From 3a682b73088ccb6d54925dd815138417ec1ac4cf Mon Sep 17 00:00:00 2001 From: Joosep Pata Date: Wed, 3 Jul 2024 14:55:08 +0300 Subject: [PATCH 4/8] keep/drop --- CLDConfig/CLDReconstruction.py | 18 +++++++++++++++++- 1 file changed, 17 insertions(+), 1 deletion(-) diff --git a/CLDConfig/CLDReconstruction.py b/CLDConfig/CLDReconstruction.py index 87631a5..89d0c2f 100644 --- a/CLDConfig/CLDReconstruction.py +++ b/CLDConfig/CLDReconstruction.py @@ -1117,7 +1117,23 @@ from Configurables import PodioOutput out = PodioOutput("PodioOutput", filename = f"{reco_args.outputBasename}_edm4hep.root") - out.outputCommands = ["keep *"] + out.outputCommands = [ + "drop *", + "keep *MCParticles*", + "keep *Pandora*", + "keep *SiTracks*", + "keep *MCTruthLink*", + "keep *Clusters*", + "keep *RecoParticles*", + "keep *MUON*", + "keep *ECALBarrel*", + "keep *ECALEndcap*", + "keep *HCALBarrel*", + "keep *HCALEndcap*", + "keep *HCALOther*", + "keep *TrackerHit*", + "keep *Vertices*", + ] algList.append(out) # We need to convert the inputs in case we have EDM4hep input From 982a1601e111feca4ccf4c4fcc6571d9a8f19d87 Mon Sep 17 00:00:00 2001 From: Joosep Pata Date: Mon, 7 Oct 2024 10:19:39 +0300 Subject: [PATCH 5/8] up --- CLDConfig/p8_ee_tt_ecm240.cmd | 33 +++++++++++++++++++++++++++++++++ 1 file changed, 33 insertions(+) create mode 100644 CLDConfig/p8_ee_tt_ecm240.cmd diff --git a/CLDConfig/p8_ee_tt_ecm240.cmd b/CLDConfig/p8_ee_tt_ecm240.cmd new file mode 100644 index 0000000..f146992 --- /dev/null +++ b/CLDConfig/p8_ee_tt_ecm240.cmd @@ -0,0 +1,33 @@ +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 + +Beams:allowMomentumSpread = off + +! Vertex smearing : +Beams:allowVertexSpread = on +Beams:sigmaVertexX = 27.3e-3 +Beams:sigmaVertexY = 48.8E-6 +Beams:sigmaVertexZ = 1.33 +Beams:sigmaTime = 1.94 ! 6.46 ps + + +! 3) Hard process : tt at 365 GeV +Beams:eCM = 240 ! CM energy of collision +Top:ffbar2ttbar(s:gmZ) = 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 + From de66e8725759d708829cf845251a69e11fa0d6b8 Mon Sep 17 00:00:00 2001 From: hart Date: Tue, 16 Sep 2025 22:26:52 +0200 Subject: [PATCH 6/8] update config to retain everything --- CLDConfig/CLDReconstruction.py | 16 +--------------- CLDConfig/CLDReconstruction.xml | 10 ---------- 2 files changed, 1 insertion(+), 25 deletions(-) diff --git a/CLDConfig/CLDReconstruction.py b/CLDConfig/CLDReconstruction.py index 89d0c2f..297b780 100644 --- a/CLDConfig/CLDReconstruction.py +++ b/CLDConfig/CLDReconstruction.py @@ -1118,21 +1118,7 @@ from Configurables import PodioOutput out = PodioOutput("PodioOutput", filename = f"{reco_args.outputBasename}_edm4hep.root") out.outputCommands = [ - "drop *", - "keep *MCParticles*", - "keep *Pandora*", - "keep *SiTracks*", - "keep *MCTruthLink*", - "keep *Clusters*", - "keep *RecoParticles*", - "keep *MUON*", - "keep *ECALBarrel*", - "keep *ECALEndcap*", - "keep *HCALBarrel*", - "keep *HCALEndcap*", - "keep *HCALOther*", - "keep *TrackerHit*", - "keep *Vertices*", + "keep *", ] algList.append(out) diff --git a/CLDConfig/CLDReconstruction.xml b/CLDConfig/CLDReconstruction.xml index 65272cc..ccb43c6 100644 --- a/CLDConfig/CLDReconstruction.xml +++ b/CLDConfig/CLDReconstruction.xml @@ -1474,16 +1474,6 @@ MCPhysicsParticles - MCParticle - LCRelation - SimCalorimeterHit - CalorimeterHit - SimTrackerHit - TrackerHit - TrackerHitPlane - Track - ReconstructedParticle - LCFloatVec MCParticlesSkimmed From b0cbcdf5e94f5f0ca84d3188efe1d7b613fc8fe2 Mon Sep 17 00:00:00 2001 From: Gabriel Facini Date: Thu, 30 Oct 2025 15:20:06 +0100 Subject: [PATCH 7/8] more samples --- CLDConfig/p8_ee_tt_dilep_ecm365.cmd | 44 +++++++++++++++++++++++++++ CLDConfig/p8_ee_ww_lvqq_ecm365.cmd | 47 +++++++++++++++++++++++++++++ CLDConfig/p8_ee_zh_vvbb_ecm365.cmd | 44 +++++++++++++++++++++++++++ CLDConfig/p8_ee_zz_vvqq_ecm365.cmd | 43 ++++++++++++++++++++++++++ 4 files changed, 178 insertions(+) create mode 100644 CLDConfig/p8_ee_tt_dilep_ecm365.cmd create mode 100644 CLDConfig/p8_ee_ww_lvqq_ecm365.cmd create mode 100644 CLDConfig/p8_ee_zh_vvbb_ecm365.cmd create mode 100644 CLDConfig/p8_ee_zz_vvqq_ecm365.cmd diff --git a/CLDConfig/p8_ee_tt_dilep_ecm365.cmd b/CLDConfig/p8_ee_tt_dilep_ecm365.cmd new file mode 100644 index 0000000..a45ca30 --- /dev/null +++ b/CLDConfig/p8_ee_tt_dilep_ecm365.cmd @@ -0,0 +1,44 @@ +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 + +Beams:allowMomentumSpread = off + +! Vertex smearing : +Beams:allowVertexSpread = on +Beams:sigmaVertexX = 27.3e-3 +Beams:sigmaVertexY = 48.8E-6 +Beams:sigmaVertexZ = 1.33 +Beams:sigmaTime = 1.94 ! 6.46 ps + + +! 3) Hard process : tt at 365 GeV +Beams:eCM = 365 ! CM energy of collision +Top:ffbar2ttbar(s:gmZ) = on + +! Control W decay modes +24:onMode = off +24:onIfMatch = 11 -12 ! W+ → e⁺ νₑ +24:onIfMatch = 13 -14 ! W+ → μ⁺ ν_μ +! 24:onIfAny = 1 2 3 4 ! W→qq + +-24:onMode = off +-24:onIfMatch = -11 12 ! W− → e⁻ anti-νₑ +-24:onIfMatch = -13 14 ! W− → μ⁻ anti-ν_μ +! -24:onIfAny = 1 2 3 4 ! W→qq + +! 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 + diff --git a/CLDConfig/p8_ee_ww_lvqq_ecm365.cmd b/CLDConfig/p8_ee_ww_lvqq_ecm365.cmd new file mode 100644 index 0000000..bdd0e84 --- /dev/null +++ b/CLDConfig/p8_ee_ww_lvqq_ecm365.cmd @@ -0,0 +1,47 @@ +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 + +Beams:allowMomentumSpread = off + +! Vertex smearing : +Beams:allowVertexSpread = on +Beams:sigmaVertexX = 27.3e-3 +Beams:sigmaVertexY = 48.8E-6 +Beams:sigmaVertexZ = 1.33 +Beams:sigmaTime = 1.94 ! 6.46 ps + + +! 3) Hard process : WW at 365 GeV +Beams:eCM = 365 ! CM energy of collision +WeakDoubleBoson:ffbar2gmZgmZ = on +WeakDoubleBoson:ffbar2WW = on +WeakDoubleBoson:ffbar2ZW = off + +! Control W decay modes +! => W → (mu,nu) or (udsc) +| W-> muv = 0.1, W->qq = 0.67 +! 0.2% muvmuv, 23% muvqq, 76% qqqq +24:onMode = off +24:onIfMatch = 13 -14 ! W -> mu+ nu_mu (and charge-conjugate) +24:onIfAny = 1 2 3 4 ! W -> qq′ (u,d,s,c) + +-24:onMode = off +-24:onIfMatch = -13 14 ! W− → μ⁻ anti-ν_μ +-24:onIfAny = 1 2 3 4 ! W→qq + +! 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 + diff --git a/CLDConfig/p8_ee_zh_vvbb_ecm365.cmd b/CLDConfig/p8_ee_zh_vvbb_ecm365.cmd new file mode 100644 index 0000000..bc3528a --- /dev/null +++ b/CLDConfig/p8_ee_zh_vvbb_ecm365.cmd @@ -0,0 +1,44 @@ +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 + +Beams:allowMomentumSpread = off + +! Vertex smearing : +Beams:allowVertexSpread = on +Beams:sigmaVertexX = 27.3e-3 +Beams:sigmaVertexY = 48.8E-6 +Beams:sigmaVertexZ = 1.33 +Beams:sigmaTime = 1.94 ! 6.46 ps + + +! 3) Hard process : ZH at 365 GeV (not the Higgs factory of 240 GeV) +Beams:eCM = 365 ! CM energy of collision +HiggsSM:ffbar2HZ = on +HiggsSM:gg2H = off +HiggsSM:ffbar2H = off + +! Control Z decay modes +23:onMode = off +23:onIfAny = 12 14 16 ! Z->vv + +! Control H decay modes +25:onMode = off +25:onIfMatch = 5 -5 ! H->bb +25:onIfMatch = 4 -4 ! H->cc + +! 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 + diff --git a/CLDConfig/p8_ee_zz_vvqq_ecm365.cmd b/CLDConfig/p8_ee_zz_vvqq_ecm365.cmd new file mode 100644 index 0000000..c26dc07 --- /dev/null +++ b/CLDConfig/p8_ee_zz_vvqq_ecm365.cmd @@ -0,0 +1,43 @@ +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 + +Beams:allowMomentumSpread = off + +! Vertex smearing : +Beams:allowVertexSpread = on +Beams:sigmaVertexX = 27.3e-3 +Beams:sigmaVertexY = 48.8E-6 +Beams:sigmaVertexZ = 1.33 +Beams:sigmaTime = 1.94 ! 6.46 ps + + +! 3) Hard process : ZZ at 365 GeV +Beams:eCM = 365 ! CM energy of collision +WeakDoubleBoson:ffbar2gmZgmZ = on +WeakDoubleBoson:ffbar2WW = off +WeakDoubleBoson:ffbar2ZW = off + +! Control Z decay modes +! => Z → (u,d,s,c,b) or (v_e,v_mu,v_tau) will get qq+qq, qq+vv, vv+vv +! Z->qq = 0.7, Z->vv = 0.2 +! 35% vvqq, 5% vvvv, 60% qqqq +23:onMode = off +23:onIfAny = 1 2 3 4 5 12 14 16 + + +! 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 + From 945d60c27b10f721829db26057ef237940259e8e Mon Sep 17 00:00:00 2001 From: gfacini Date: Thu, 11 Dec 2025 12:31:32 +0100 Subject: [PATCH 8/8] Update p8_ee_ww_lvqq_ecm365.cmd turn of zz process --- CLDConfig/p8_ee_ww_lvqq_ecm365.cmd | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/CLDConfig/p8_ee_ww_lvqq_ecm365.cmd b/CLDConfig/p8_ee_ww_lvqq_ecm365.cmd index bdd0e84..bfb9727 100644 --- a/CLDConfig/p8_ee_ww_lvqq_ecm365.cmd +++ b/CLDConfig/p8_ee_ww_lvqq_ecm365.cmd @@ -25,7 +25,7 @@ Beams:sigmaTime = 1.94 ! 6.46 ps ! 3) Hard process : WW at 365 GeV Beams:eCM = 365 ! CM energy of collision -WeakDoubleBoson:ffbar2gmZgmZ = on +WeakDoubleBoson:ffbar2gmZgmZ = off WeakDoubleBoson:ffbar2WW = on WeakDoubleBoson:ffbar2ZW = off