Madgraph plus Pythia, is it feasible to generate SM photon in madgraph and customize it to dark photon in pythia
Hello,
I am trying to generate dark photon which is similar to SM photon, but just want it to be massive and change its kinetic fixing factor. I would like to know if it is possible to generate the process with standard model particles: p p > mu+ mu- a by using madgraph, and later on customize the SM photon to dark photon by using pythia?
If it is feasible, then it will be a pythia/Geant4 question. Or if you know some other way of doing it, please let me know.
I have done the SM process generation and generate the lhe file, and I got stuck at customizing it in pythia. I am using CMSSW, and I find some code in SimG4Core/
So far I changed all the pdgId of SM photon to dark photon (22->1072000) in lhe file, and customize the process as following details, but pythia pops up the errors:
*------- PYTHIA Error and Warning Messages Statistics -------
| |
| times message |
| |
| 100 Abort from Pythia::next: parton+hadronLevel failed; giving up |
| 1937 Error in Pythia::check: unknown particle code |
| 1000 Error in Pythia::next: check of event revealed problems |
| 3 Error in SpaceShower:
| 1 Info from SLHAinterface:
| 1 Info from SLHAinterface:
| 2 Warning in SLHAinterface:
| 2 Warning in SLHAinterface:
| 1 Warning in SLHAinterface:
| 13 Warning in SpaceShower:
| 30 Warning in SpaceShower:
| 2 Warning in StringFragmenta
| |
*------- End PYTHIA Error and Warning Messages Statistics -------
~~~~~~~
The code is the followings:
import FWCore.
from Configuration.
process = cms.Process(
# import of standard configurations
process.
process.
process.
process.
process.
process.
process.
process.
process.
process.
process.
process.
process.
process.
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.
)
# Input source
process.source = cms.Source(
dropDescend
fileNames = cms.untracked.
)
process.options = cms.untracked.PSet(
)
# Output definition
process.
SelectEvents = cms.untracked.PSet(
),
compression
compression
dataset = cms.untracked.PSet(
dataTier = cms.untracked.
filterName = cms.untracked.
),
eventAutoFl
fileName = cms.untracked.
outputCommands = process.
splitLevel = cms.untracked.
)
# Additional output definition
# Other statements
process.
process.
from Configuration.
process.GlobalTag = GlobalTag(
from Configuration.
from Configuration.
process.generator = cms.EDFilter(
),
)
# Path and EndPath definitions
process.
process.
process.
process.endjob_step = cms.EndPath(
process.
process.schedule = cms.Schedule(
from PhysicsTools.
associatePatAlg
#Setup FWK for multithreaded
#process.
process.
# filter all path with the production filter sequence
for path in process.paths:
if path in ['lhe_step']: continue
getattr(
# customisation of the process.
from SimG4Core.
process = customiseDarkPh
# End of customisation functions
# Customisation from command line
# Add early deletion of temporary data products to reduce peak memory need
from Configuration.
process = customiseEarlyD
# End adding early deletion
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Yao Yao for more information if necessary.