Commit 4a205a29 authored by Carlos Vazquez Sierra's avatar Carlos Vazquez Sierra 👁‍🗨
Browse files

Add HiddenValley decay files for various modes

parent 56bec29b
Pipeline #1724127 failed with stages
in 27 minutes and 53 seconds
# EventType: 44102000
#
# Descriptor: pivDiag -> K+K-
#
# NickName: HiddenValley_HiggsPortal_KK_1.02GeV_5ps
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/HiddenValleyDiKaon
#
# InsertPythonCode:
# from Gaudi.Configuration import *
# Generation().PileUpTool = "FixedLuminosityForRareProcess"
# importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
# from Configurables import Pythia8Production
# Generation().Special.ProductionTool = "Pythia8Production"
# Generation().Special.addTool( Pythia8Production )
#
# ## Pythia8 production commands
# GlobalParameters = ["PartonLevel:MPI = on",
# "PartonLevel:ISR = on",
# "PartonLevel:FSR = on",
# "HadronLevel:Hadronize = on"]
# HiggsProduction = ["HiggsSM:all = on",
# "25:onMode = off",
# "25:addChannel = 1 1 103 4900101 -4900101"]
# HidValleyParameters = ["HiddenValley:FSR = on",
# "HiddenValley:fragment = on",
# "HiddenValley:alphaOrder = 1",
# "HiddenValley:alphaFSR = 0.7",
# "HiddenValley:Ngauge = 3",
# "HiddenValley:Lambda = 4",
# "HiddenValley:pTminFSR = 1",
# "HiddenValley:nFlav = 1",
# "HiddenValley:probVector = 0"]
# HidPartParameters = ["4900101:m0 = 1", # qv mass
# "4900113:m0 = 1.02", # Rhov mass
# "4900111:m0 = 1.02", # Piv mass (GeV)
# "4900111:tau0 = 1.49896", # Piv lifetime (5 ps) in mm (c = 0.299792 mm/ps)
# "4900111:onMode = on",
# "4900111:mayDecay = on",
# "4900111:addChannel 1 1 1 321 -321", # K+ K-
# "321:mayDecay = on"]
# Generation().Special.Pythia8Production.Commands += GlobalParameters
# Generation().Special.Pythia8Production.Commands += HiggsProduction
# Generation().Special.Pythia8Production.Commands += HidValleyParameters
# Generation().Special.Pythia8Production.Commands += HidPartParameters
#
# ## Property Particle table update:
# from Configurables import LHCb__ParticlePropertySvc as PPS
# PPS().Particles = ["gv 4900021 4900021 0.0 0.0 0.0 gv 4900021 0",
# "Zv 4900023 4900023 0.0 125.0 6.58e-23 Zv 4900023 1e-2",
# "qv 4900101 4900101 0.0 1.0 0.0 qv 4900101 0",
# "qvbar -4900101 -4900101 0.0 1.0 0.0 qvbar -4900101 0",
# "pivDiag 4900111 4900111 0.0 1.02 5e-12 pivDiag 4900111 2.0e-11",
# "pivUp 4900211 4900211 0.0 1.02 0.0 pivUp 4900211 0",
# "pivDn -4900211 -4900211 0.0 1.02 0.0 pivDn -4900211 0",
# "rhovDiag 4900113 4900113 0.0 1.02 0.0 rhovDiag 4900113 0",
# "rhovUp 4900213 4900213 0.0 1.02 0.0 rhovUp 4900213 0",
# "rhovDn -4900213 -4900213 0.0 1.02 0.0 rhovDn -4900213 0"]
#
# ## Geant4 classes:
# from Configurables import GenerationToSimulation
# from Configurables import GiGaHiddenValleyParticles, GiGaPhysListModular
# GenerationToSimulation("GenToSim").KeepCode = "in_list( GABSID, [ 'pivDiag', 'gv', 'qv', 'qvbar' ] )"
# gigaHVPart = GiGaHiddenValleyParticles()
# gigaHVPart.HiddenValley = ["pivDiag", "gv", "qv", "qvbar"]
# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHVPart ]
#
# ## Generator level cuts:
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "HiddenValleyDiKaon" )
# tracksInAcc = Generation().HiddenValleyDiKaon
# tracksInAcc.Code = " count ( isGoodHV ) > 1 "
# ### - HepMC::IteratorRange::descendants 4
# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodKaon = ( ( GP > 2.*GeV ) & ( GTHETA < 400.0*mrad ) & ( 'K+' == GABSID ) )"
# , "isGoodHV = ( ( 'pivDiag' == GID ) & ( GNINTREE( isGoodKaon, 4 ) > 1 ) )" ]
# EndInsertPythonCode
#
# Documentation: HV pion decaying into KK, at least two HV pions in the event, 1.02 GeV and 5 ps.
# EndDocumentation
#
# PhysicsWG: Exotica
# Tested: Yes
# Responsible: Carlos Vazquez Sierra
# Email: carlos.vazquez.sierra@cern.ch
# Date: 20200617
# CPUTime: < 1 min
#
End
#
# EventType: 44102001
#
# Descriptor: pivDiag -> K+K-
#
# NickName: HiddenValley_HiggsPortal_KK_1.08GeV_5ps
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/HiddenValleyDiKaon
#
# InsertPythonCode:
# from Gaudi.Configuration import *
# Generation().PileUpTool = "FixedLuminosityForRareProcess"
# importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
# from Configurables import Pythia8Production
# Generation().Special.ProductionTool = "Pythia8Production"
# Generation().Special.addTool( Pythia8Production )
#
# ## Pythia8 production commands
# GlobalParameters = ["PartonLevel:MPI = on",
# "PartonLevel:ISR = on",
# "PartonLevel:FSR = on",
# "HadronLevel:Hadronize = on"]
# HiggsProduction = ["HiggsSM:all = on",
# "25:onMode = off",
# "25:addChannel = 1 1 103 4900101 -4900101"]
# HidValleyParameters = ["HiddenValley:FSR = on",
# "HiddenValley:fragment = on",
# "HiddenValley:alphaOrder = 1",
# "HiddenValley:alphaFSR = 0.7",
# "HiddenValley:Ngauge = 3",
# "HiddenValley:Lambda = 4",
# "HiddenValley:pTminFSR = 1",
# "HiddenValley:nFlav = 1",
# "HiddenValley:probVector = 0"]
# HidPartParameters = ["4900101:m0 = 1", # qv mass
# "4900113:m0 = 1.08", # Rhov mass
# "4900111:m0 = 1.08", # Piv mass (GeV)
# "4900111:tau0 = 1.49896", # Piv lifetime (5 ps) in mm (c = 0.299792 mm/ps)
# "4900111:onMode = on",
# "4900111:mayDecay = on",
# "4900111:addChannel 1 1 1 321 -321", # K+ K-
# "321:mayDecay = on"]
# Generation().Special.Pythia8Production.Commands += GlobalParameters
# Generation().Special.Pythia8Production.Commands += HiggsProduction
# Generation().Special.Pythia8Production.Commands += HidValleyParameters
# Generation().Special.Pythia8Production.Commands += HidPartParameters
#
# ## Property Particle table update:
# from Configurables import LHCb__ParticlePropertySvc as PPS
# PPS().Particles = ["gv 4900021 4900021 0.0 0.0 0.0 gv 4900021 0",
# "Zv 4900023 4900023 0.0 125.0 6.58e-23 Zv 4900023 1e-2",
# "qv 4900101 4900101 0.0 1.0 0.0 qv 4900101 0",
# "qvbar -4900101 -4900101 0.0 1.0 0.0 qvbar -4900101 0",
# "pivDiag 4900111 4900111 0.0 1.08 5e-12 pivDiag 4900111 2.0e-11",
# "pivUp 4900211 4900211 0.0 1.08 0.0 pivUp 4900211 0",
# "pivDn -4900211 -4900211 0.0 1.08 0.0 pivDn -4900211 0",
# "rhovDiag 4900113 4900113 0.0 1.08 0.0 rhovDiag 4900113 0",
# "rhovUp 4900213 4900213 0.0 1.08 0.0 rhovUp 4900213 0",
# "rhovDn -4900213 -4900213 0.0 1.08 0.0 rhovDn -4900213 0"]
#
# ## Geant4 classes:
# from Configurables import GenerationToSimulation
# from Configurables import GiGaHiddenValleyParticles, GiGaPhysListModular
# GenerationToSimulation("GenToSim").KeepCode = "in_list( GABSID, [ 'pivDiag', 'gv', 'qv', 'qvbar' ] )"
# gigaHVPart = GiGaHiddenValleyParticles()
# gigaHVPart.HiddenValley = ["pivDiag", "gv", "qv", "qvbar"]
# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHVPart ]
#
# ## Generator level cuts:
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "HiddenValleyDiKaon" )
# tracksInAcc = Generation().HiddenValleyDiKaon
# tracksInAcc.Code = " count ( isGoodHV ) > 1 "
# ### - HepMC::IteratorRange::descendants 4
# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodKaon = ( ( GP > 2.*GeV ) & ( GTHETA < 400.0*mrad ) & ( 'K+' == GABSID ) )"
# , "isGoodHV = ( ( 'pivDiag' == GID ) & ( GNINTREE( isGoodKaon, 4 ) > 1 ) )" ]
# EndInsertPythonCode
#
# Documentation: HV pion decaying into KK, at least two HV pions in the event, 1.08 GeV and 5 ps.
# EndDocumentation
#
# PhysicsWG: Exotica
# Tested: Yes
# Responsible: Carlos Vazquez Sierra
# Email: carlos.vazquez.sierra@cern.ch
# Date: 20200617
# CPUTime: < 1 min
#
End
#
# EventType: 44102002
#
# Descriptor: pivDiag -> K+K-
#
# NickName: HiddenValley_HiggsPortal_KK_1.3GeV_5ps
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/HiddenValleyDiKaon
#
# InsertPythonCode:
# from Gaudi.Configuration import *
# Generation().PileUpTool = "FixedLuminosityForRareProcess"
# importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
# from Configurables import Pythia8Production
# Generation().Special.ProductionTool = "Pythia8Production"
# Generation().Special.addTool( Pythia8Production )
#
# ## Pythia8 production commands
# GlobalParameters = ["PartonLevel:MPI = on",
# "PartonLevel:ISR = on",
# "PartonLevel:FSR = on",
# "HadronLevel:Hadronize = on"]
# HiggsProduction = ["HiggsSM:all = on",
# "25:onMode = off",
# "25:addChannel = 1 1 103 4900101 -4900101"]
# HidValleyParameters = ["HiddenValley:FSR = on",
# "HiddenValley:fragment = on",
# "HiddenValley:alphaOrder = 1",
# "HiddenValley:alphaFSR = 0.7",
# "HiddenValley:Ngauge = 3",
# "HiddenValley:Lambda = 4",
# "HiddenValley:pTminFSR = 1",
# "HiddenValley:nFlav = 1",
# "HiddenValley:probVector = 0"]
# HidPartParameters = ["4900101:m0 = 1", # qv mass
# "4900113:m0 = 1.3", # Rhov mass
# "4900111:m0 = 1.3", # Piv mass (GeV)
# "4900111:tau0 = 1.49896", # Piv lifetime (5 ps) in mm (c = 0.299792 mm/ps)
# "4900111:onMode = on",
# "4900111:mayDecay = on",
# "4900111:addChannel 1 1 1 321 -321", # K+ K-
# "321:mayDecay = on"]
# Generation().Special.Pythia8Production.Commands += GlobalParameters
# Generation().Special.Pythia8Production.Commands += HiggsProduction
# Generation().Special.Pythia8Production.Commands += HidValleyParameters
# Generation().Special.Pythia8Production.Commands += HidPartParameters
#
# ## Property Particle table update:
# from Configurables import LHCb__ParticlePropertySvc as PPS
# PPS().Particles = ["gv 4900021 4900021 0.0 0.0 0.0 gv 4900021 0",
# "Zv 4900023 4900023 0.0 125.0 6.58e-23 Zv 4900023 1e-2",
# "qv 4900101 4900101 0.0 1.0 0.0 qv 4900101 0",
# "qvbar -4900101 -4900101 0.0 1.0 0.0 qvbar -4900101 0",
# "pivDiag 4900111 4900111 0.0 1.3 5e-12 pivDiag 4900111 2.0e-11",
# "pivUp 4900211 4900211 0.0 1.3 0.0 pivUp 4900211 0",
# "pivDn -4900211 -4900211 0.0 1.3 0.0 pivDn -4900211 0",
# "rhovDiag 4900113 4900113 0.0 1.3 0.0 rhovDiag 4900113 0",
# "rhovUp 4900213 4900213 0.0 1.3 0.0 rhovUp 4900213 0",
# "rhovDn -4900213 -4900213 0.0 1.3 0.0 rhovDn -4900213 0"]
#
# ## Geant4 classes:
# from Configurables import GenerationToSimulation
# from Configurables import GiGaHiddenValleyParticles, GiGaPhysListModular
# GenerationToSimulation("GenToSim").KeepCode = "in_list( GABSID, [ 'pivDiag', 'gv', 'qv', 'qvbar' ] )"
# gigaHVPart = GiGaHiddenValleyParticles()
# gigaHVPart.HiddenValley = ["pivDiag", "gv", "qv", "qvbar"]
# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHVPart ]
#
# ## Generator level cuts:
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "HiddenValleyDiKaon" )
# tracksInAcc = Generation().HiddenValleyDiKaon
# tracksInAcc.Code = " count ( isGoodHV ) > 1 "
# ### - HepMC::IteratorRange::descendants 4
# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodKaon = ( ( GP > 2.*GeV ) & ( GTHETA < 400.0*mrad ) & ( 'K+' == GABSID ) )"
# , "isGoodHV = ( ( 'pivDiag' == GID ) & ( GNINTREE( isGoodKaon, 4 ) > 1 ) )" ]
# EndInsertPythonCode
#
# Documentation: HV pion decaying into KK, at least two HV pions in the event, 1.3 GeV and 5 ps.
# EndDocumentation
#
# PhysicsWG: Exotica
# Tested: Yes
# Responsible: Carlos Vazquez Sierra
# Email: carlos.vazquez.sierra@cern.ch
# Date: 20200617
# CPUTime: < 1 min
#
End
#
# EventType: 44102003
#
# Descriptor: pivDiag -> K+K-
#
# NickName: HiddenValley_HiggsPortal_KK_1.8GeV_5ps
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/HiddenValleyDiKaon
#
# InsertPythonCode:
# from Gaudi.Configuration import *
# Generation().PileUpTool = "FixedLuminosityForRareProcess"
# importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
# from Configurables import Pythia8Production
# Generation().Special.ProductionTool = "Pythia8Production"
# Generation().Special.addTool( Pythia8Production )
#
# ## Pythia8 production commands
# GlobalParameters = ["PartonLevel:MPI = on",
# "PartonLevel:ISR = on",
# "PartonLevel:FSR = on",
# "HadronLevel:Hadronize = on"]
# HiggsProduction = ["HiggsSM:all = on",
# "25:onMode = off",
# "25:addChannel = 1 1 103 4900101 -4900101"]
# HidValleyParameters = ["HiddenValley:FSR = on",
# "HiddenValley:fragment = on",
# "HiddenValley:alphaOrder = 1",
# "HiddenValley:alphaFSR = 0.7",
# "HiddenValley:Ngauge = 3",
# "HiddenValley:Lambda = 4",
# "HiddenValley:pTminFSR = 1",
# "HiddenValley:nFlav = 1",
# "HiddenValley:probVector = 0"]
# HidPartParameters = ["4900101:m0 = 1", # qv mass
# "4900113:m0 = 1.8", # Rhov mass
# "4900111:m0 = 1.8", # Piv mass (GeV)
# "4900111:tau0 = 1.49896", # Piv lifetime (5 ps) in mm (c = 0.299792 mm/ps)
# "4900111:onMode = on",
# "4900111:mayDecay = on",
# "4900111:addChannel 1 1 1 321 -321", # K+ K-
# "321:mayDecay = on"]
# Generation().Special.Pythia8Production.Commands += GlobalParameters
# Generation().Special.Pythia8Production.Commands += HiggsProduction
# Generation().Special.Pythia8Production.Commands += HidValleyParameters
# Generation().Special.Pythia8Production.Commands += HidPartParameters
#
# ## Property Particle table update:
# from Configurables import LHCb__ParticlePropertySvc as PPS
# PPS().Particles = ["gv 4900021 4900021 0.0 0.0 0.0 gv 4900021 0",
# "Zv 4900023 4900023 0.0 125.0 6.58e-23 Zv 4900023 1e-2",
# "qv 4900101 4900101 0.0 1.0 0.0 qv 4900101 0",
# "qvbar -4900101 -4900101 0.0 1.0 0.0 qvbar -4900101 0",
# "pivDiag 4900111 4900111 0.0 1.8 5e-12 pivDiag 4900111 2.0e-11",
# "pivUp 4900211 4900211 0.0 1.8 0.0 pivUp 4900211 0",
# "pivDn -4900211 -4900211 0.0 1.8 0.0 pivDn -4900211 0",
# "rhovDiag 4900113 4900113 0.0 1.8 0.0 rhovDiag 4900113 0",
# "rhovUp 4900213 4900213 0.0 1.8 0.0 rhovUp 4900213 0",
# "rhovDn -4900213 -4900213 0.0 1.8 0.0 rhovDn -4900213 0"]
#
# ## Geant4 classes:
# from Configurables import GenerationToSimulation
# from Configurables import GiGaHiddenValleyParticles, GiGaPhysListModular
# GenerationToSimulation("GenToSim").KeepCode = "in_list( GABSID, [ 'pivDiag', 'gv', 'qv', 'qvbar' ] )"
# gigaHVPart = GiGaHiddenValleyParticles()
# gigaHVPart.HiddenValley = ["pivDiag", "gv", "qv", "qvbar"]
# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHVPart ]
#
# ## Generator level cuts:
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "HiddenValleyDiKaon" )
# tracksInAcc = Generation().HiddenValleyDiKaon
# tracksInAcc.Code = " count ( isGoodHV ) > 1 "
# ### - HepMC::IteratorRange::descendants 4
# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodKaon = ( ( GP > 2.*GeV ) & ( GTHETA < 400.0*mrad ) & ( 'K+' == GABSID ) )"
# , "isGoodHV = ( ( 'pivDiag' == GID ) & ( GNINTREE( isGoodKaon, 4 ) > 1 ) )" ]
# EndInsertPythonCode
#
# Documentation: HV pion decaying into KK, at least two HV pions in the event, 1.8 GeV and 5 ps.
# EndDocumentation
#
# PhysicsWG: Exotica
# Tested: Yes
# Responsible: Carlos Vazquez Sierra
# Email: carlos.vazquez.sierra@cern.ch
# Date: 20200617
# CPUTime: < 1 min
#
End
#
# EventType: 44102004
#
# Descriptor: pivDiag -> K+K-
#
# NickName: HiddenValley_HiggsPortal_KK_2.5GeV_5ps
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/HiddenValleyDiKaon
#
# InsertPythonCode:
# from Gaudi.Configuration import *
# Generation().PileUpTool = "FixedLuminosityForRareProcess"
# importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
# from Configurables import Pythia8Production
# Generation().Special.ProductionTool = "Pythia8Production"
# Generation().Special.addTool( Pythia8Production )
#
# ## Pythia8 production commands
# GlobalParameters = ["PartonLevel:MPI = on",
# "PartonLevel:ISR = on",
# "PartonLevel:FSR = on",
# "HadronLevel:Hadronize = on"]
# HiggsProduction = ["HiggsSM:all = on",
# "25:onMode = off",
# "25:addChannel = 1 1 103 4900101 -4900101"]
# HidValleyParameters = ["HiddenValley:FSR = on",
# "HiddenValley:fragment = on",
# "HiddenValley:alphaOrder = 1",
# "HiddenValley:alphaFSR = 0.7",
# "HiddenValley:Ngauge = 3",
# "HiddenValley:Lambda = 4",
# "HiddenValley:pTminFSR = 1",
# "HiddenValley:nFlav = 1",
# "HiddenValley:probVector = 0"]
# HidPartParameters = ["4900101:m0 = 1", # qv mass
# "4900113:m0 = 2.5", # Rhov mass
# "4900111:m0 = 2.5", # Piv mass (GeV)
# "4900111:tau0 = 1.49896", # Piv lifetime (5 ps) in mm (c = 0.299792 mm/ps)
# "4900111:onMode = on",
# "4900111:mayDecay = on",
# "4900111:addChannel 1 1 1 321 -321", # K+ K-
# "321:mayDecay = on"]
# Generation().Special.Pythia8Production.Commands += GlobalParameters
# Generation().Special.Pythia8Production.Commands += HiggsProduction
# Generation().Special.Pythia8Production.Commands += HidValleyParameters
# Generation().Special.Pythia8Production.Commands += HidPartParameters
#
# ## Property Particle table update:
# from Configurables import LHCb__ParticlePropertySvc as PPS
# PPS().Particles = ["gv 4900021 4900021 0.0 0.0 0.0 gv 4900021 0",
# "Zv 4900023 4900023 0.0 125.0 6.58e-23 Zv 4900023 1e-2",
# "qv 4900101 4900101 0.0 1.0 0.0 qv 4900101 0",
# "qvbar -4900101 -4900101 0.0 1.0 0.0 qvbar -4900101 0",
# "pivDiag 4900111 4900111 0.0 2.5 5e-12 pivDiag 4900111 2.0e-11",
# "pivUp 4900211 4900211 0.0 2.5 0.0 pivUp 4900211 0",
# "pivDn -4900211 -4900211 0.0 2.5 0.0 pivDn -4900211 0",
# "rhovDiag 4900113 4900113 0.0 2.5 0.0 rhovDiag 4900113 0",
# "rhovUp 4900213 4900213 0.0 2.5 0.0 rhovUp 4900213 0",
# "rhovDn -4900213 -4900213 0.0 2.5 0.0 rhovDn -4900213 0"]
#
# ## Geant4 classes:
# from Configurables import GenerationToSimulation
# from Configurables import GiGaHiddenValleyParticles, GiGaPhysListModular
# GenerationToSimulation("GenToSim").KeepCode = "in_list( GABSID, [ 'pivDiag', 'gv', 'qv', 'qvbar' ] )"
# gigaHVPart = GiGaHiddenValleyParticles()
# gigaHVPart.HiddenValley = ["pivDiag", "gv", "qv", "qvbar"]
# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHVPart ]
#
# ## Generator level cuts:
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "HiddenValleyDiKaon" )
# tracksInAcc = Generation().HiddenValleyDiKaon
# tracksInAcc.Code = " count ( isGoodHV ) > 1 "
# ### - HepMC::IteratorRange::descendants 4
# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodKaon = ( ( GP > 2.*GeV ) & ( GTHETA < 400.0*mrad ) & ( 'K+' == GABSID ) )"
# , "isGoodHV = ( ( 'pivDiag' == GID ) & ( GNINTREE( isGoodKaon, 4 ) > 1 ) )" ]
# EndInsertPythonCode
#
# Documentation: HV pion decaying into KK, at least two HV pions in the event, 2.5 GeV and 5 ps.
# EndDocumentation
#
# PhysicsWG: Exotica
# Tested: Yes
# Responsible: Carlos Vazquez Sierra
# Email: carlos.vazquez.sierra@cern.ch
# Date: 20200617
# CPUTime: < 1 min
#
End
#
# EventType: 44102009
#
# Descriptor: pivDiag -> mu+mu-
#
# NickName: HiddenValley_HiggsPortal_mumu_0.3GeV_5ps
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/HiddenValleyDiMuon
#
# InsertPythonCode:
# from Gaudi.Configuration import *
# Generation().PileUpTool = "FixedLuminosityForRareProcess"
# importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
# from Configurables import Pythia8Production
# Generation().Special.ProductionTool = "Pythia8Production"
# Generation().Special.addTool( Pythia8Production )
#
# ## Pythia8 production commands
# GlobalParameters = ["PartonLevel:MPI = on",
# "PartonLevel:ISR = on",
# "PartonLevel:FSR = on",
# "HadronLevel:Hadronize = on"]
# HiggsProduction = ["HiggsSM:all = on",
# "25:onMode = off",
# "25:addChannel = 1 1 103 4900101 -4900101"]
# HidValleyParameters = ["HiddenValley:FSR = on",
# "HiddenValley:fragment = on",
# "HiddenValley:alphaOrder = 1",
# "HiddenValley:alphaFSR = 0.7",
# "HiddenValley:Ngauge = 3",
# "HiddenValley:Lambda = 4",
# "HiddenValley:pTminFSR = 1",
# "HiddenValley:nFlav = 1",
# "HiddenValley:probVector = 0"]
# HidPartParameters = ["4900101:m0 = 1", # qv mass
# "4900113:m0 = 0.3", # Rhov mass
# "4900111:m0 = 0.3", # Piv mass (GeV)
# "4900111:tau0 = 1.49896", # Piv lifetime (5 ps) in mm (c = 0.299792 mm/ps)
# "4900111:onMode = on",
# "4900111:mayDecay = on",
# "4900111:addChannel 1 1 1 13 -13", # mu+ mu-
# "13:mayDecay = on"]
# Generation().Special.Pythia8Production.Commands += GlobalParameters
# Generation().Special.Pythia8Production.Commands += HiggsProduction
# Generation().Special.Pythia8Production.Commands += HidValleyParameters
# Generation().Special.Pythia8Production.Commands += HidPartParameters
#
# ## Property Particle table update:
# from Configurables import LHCb__ParticlePropertySvc as PPS
# PPS().Particles = ["gv 4900021 4900021 0.0 0.0 0.0 gv 4900021 0",
# "Zv 4900023 4900023 0.0 125.0 6.58e-23 Zv 4900023 1e-2",
# "qv 4900101 4900101 0.0 1.0 0.0 qv 4900101 0",
# "qvbar -4900101 -4900101 0.0 1.0 0.0 qvbar -4900101 0",
# "pivDiag 4900111 4900111 0.0 0.3 5e-12 pivDiag 4900111 2.0e-11",
# "pivUp 4900211 4900211 0.0 0.3 0.0 pivUp 4900211 0",
# "pivDn -4900211 -4900211 0.0 0.3 0.0 pivDn -4900211 0",
# "rhovDiag 4900113 4900113 0.0 0.3 0.0 rhovDiag 4900113 0",
# "rhovUp 4900213 4900213 0.0 0.3 0.0 rhovUp 4900213 0",
# "rhovDn -4900213 -4900213 0.0 0.3 0.0 rhovDn -4900213 0"]
#
# ## Geant4 classes:
# from Configurables import GenerationToSimulation
# from Configurables import GiGaHiddenValleyParticles, GiGaPhysListModular
# GenerationToSimulation("GenToSim").KeepCode = "in_list( GABSID, [ 'pivDiag', 'gv', 'qv', 'qvbar' ] )"
# gigaHVPart = GiGaHiddenValleyParticles()
# gigaHVPart.HiddenValley = ["pivDiag", "gv", "qv", "qvbar"]
# GiGaPhysListModular("ModularPL").PhysicsConstructors += [ gigaHVPart ]
#
# ## Generator level cuts:
# from Configurables import LoKi__FullGenEventCut
# Generation().addTool( LoKi__FullGenEventCut, "HiddenValleyDiMuon" )
# tracksInAcc = Generation().HiddenValleyDiMuon
# tracksInAcc.Code = " count ( isGoodHV ) > 1 "
# ### - HepMC::IteratorRange::descendants 4
# tracksInAcc.Preambulo += [ "from GaudiKernel.SystemOfUnits import GeV, mrad"
# , "isGoodPion = ( ( GP > 2.*GeV ) & ( GTHETA < 400.0*mrad ) & ( 'mu+' == GABSID ) )"
# , "isGoodHV = ( ( 'pivDiag' == GID ) & ( GNINTREE( isGoodPion, 4 ) > 1 ) )" ]
# EndInsertPythonCode
#
# Documentation: HV pion decaying into mumu, at least two HV pions in the event, 0.3 GeV and 5 ps.
# EndDocumentation
#
# PhysicsWG: Exotica
# Tested: Yes
# Responsible: Carlos Vazquez Sierra
# Email: carlos.vazquez.sierra@cern.ch
# Date: 20200617
# CPUTime: < 1 min
#
End
#
# EventType: 44102010
#
# Descriptor: pivDiag -> mu+mu-
#
# NickName: HiddenValley_HiggsPortal_mumu_0.9GeV_5ps
#
# Cuts: None
# FullEventCuts: LoKi::FullGenEventCut/HiddenValleyDiMuon
#
# InsertPythonCode:
# from Gaudi.Configuration import *
# Generation().PileUpTool = "FixedLuminosityForRareProcess"
# importOptions( "$DECFILESROOT/options/SwitchOffAllPythiaProcesses.py" )
# from Configurables import Pythia8Production
# Generation().Special.ProductionTool = "Pythia8Production"
# Generation().Special.addTool( Pythia8Production )
#
# ## Pythia8 production commands
# GlobalParameters = ["PartonLevel:MPI = on",
# "PartonLevel:ISR = on",
# "PartonLevel:FSR = on",
# "HadronLevel:Hadronize = on"]
# HiggsProduction = ["HiggsSM:all = on",
# "25:onMode = off",
# "25:addChannel = 1 1 103 4900101 -4900101"]
# HidValleyParameters = ["HiddenValley:FSR = on",
# "HiddenValley:fragment = on",
# "HiddenValley:alphaOrder = 1",
# "HiddenValley:alphaFSR = 0.7",
# "HiddenValley:Ngauge = 3",
# "HiddenValley:Lambda = 4",
# "HiddenValley:pTminFSR = 1",
# "HiddenValley:nFlav = 1",
# "HiddenValley:probVector = 0"]
# HidPartParameters = ["4900101:m0 = 1", # qv mass
# "4900113:m0 = 0.9", # Rhov mass
# "4900111:m0 = 0.9", # Piv mass (GeV)
# "4900111:tau0 = 1.49896", # Piv lifetime (5 ps) in mm (c = 0.299792 mm/ps)
# "4900111:onMode = on",
# "4900111:mayDecay = on",
# "4900111:addChannel 1 1 1 13 -13", # mu+ mu-
# "13:mayDecay = on"]
# Generation().Special.Pythia8Production.Commands += GlobalParameters
# Generation().Special.Pythia8Production.Commands += HiggsProduction
# Generation().Special.Pythia8Production.Commands += HidValleyParameters
# Generation().Special.Pythia8Production.Commands += HidPartParameters
#
# ## Property Particle table update:
# from Configurables import LHCb__ParticlePropertySvc as PPS
# PPS().Particles = ["gv 4900021 4900021 0.0 0.0 0.0 gv 4900021 0",
# "Zv 4900023 4900023 0.0 125.0 6.58e-23 Zv 4900023 1e-2",
# "qv 4900101 4900101 0.0 1.0 0.0 qv 4900101 0",
# "qvbar -4900101 -4900101 0.0 1.0 0.0 qvbar -4900101 0",
# "pivDiag 4900111 4900111 0.0 0.9 5e-12 pivDiag 4900111 2.0e-11",