Guest User

Untitled

a guest
Aug 21st, 2019
108
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. import FWCore.ParameterSet.Config as cms
  2. from Configuration.Generator.Pythia8CommonSettings_cfi import *
  3. from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
  4.  
  5. generator = cms.EDFilter("Pythia8GeneratorFilter",
  6. pythiaPylistVerbosity = cms.untracked.int32(0),
  7. pythiaHepMCVerbosity = cms.untracked.bool(False),
  8. comEnergy = cms.double(13000.0),
  9. ##crossSection = cms.untracked.double(54000000000), # Given by PYTHIA after running
  10. ##filterEfficiency = cms.untracked.double(0.004), # Given by PYTHIA after running
  11. maxEventsToPrint = cms.untracked.int32(0),
  12.  
  13. ExternalDecays = cms.PSet(
  14. EvtGen130 = cms.untracked.PSet(
  15. convertPythiaCodes = cms.untracked.bool(False),
  16. decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2014_NOLONGLIFE.DEC'),
  17. list_forced_decays = cms.vstring('theJpsi'),
  18. operates_on_particles = cms.vint32(443,-443),
  19. particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt_2014.pdl'),
  20. user_decay_embedded = cms.vstring(
  21. """
  22. Alias theJpsi J/psi
  23. ChargeConj theJpsi theJpsi
  24. #
  25. Decay theJpsi
  26. 1.00000 mu+ mu- PHOTOS VLL ;
  27. Enddecay
  28. #
  29. End
  30. """
  31. ),
  32. ),
  33. parameterSets = cms.vstring('EvtGen130')
  34. ),
  35.  
  36. PythiaParameters = cms.PSet(pythia8CommonSettingsBlock,
  37. pythia8CP5SettingsBlock,
  38. ## this is original from 2018
  39. processParameters = cms.vstring( 'SoftQCD:nonDiffractive = on',
  40. 'PTFilter:filter = on', # this turn on the filter
  41. 'PTFilter:quarkToFilter = 5', # PDG id of q quark (can be any other)
  42. 'PTFilter:scaleToFilter = 1.0' ),
  43. # ## this is what was there in 2012
  44. # processParameters = cms.vstring(
  45. # 'Bottomonium:all = on', # Quarkonia, MSEL=61
  46. # ## (sara: i think it was msel = 1), and is the equivalent in pythia6 of softQCD:nonDiffractive = on
  47. #
  48. # 'ParticleDecays:allowPhotonRadiation = on', # Turn on QED FSR
  49. # ### from the manual:
  50. # ### Allow or not photon radiations in decays to a lepton pair.
  51. # ### Note: The current default is to have radiation switched off,
  52. # ### in order to avoid double-counting of emissions if you link to
  53. # ### an external QED-radiation program, as is the norm in many collaborations.
  54. #
  55. # 'ParticleDecays:mixB = off',
  56. # '443:onMode = off', # Turn off J/psi decays
  57. # '443:onIfMatch = 13 -13' # just let J/psi -> mu+ mu-
  58. # ),
  59.  
  60. parameterSets = cms.vstring('pythia8CommonSettings',
  61. 'pythia8CP5Settings',
  62. 'processParameters',
  63. )
  64. )
  65. )
  66.  
  67.  
  68.  
  69.  
  70. ###########
  71. # Filters #
  72. ###########
  73. bfilter = cms.EDFilter("PythiaFilter",
  74. ParticleID = cms.untracked.int32(5)
  75. )
  76.  
  77. jpsifilter = cms.EDFilter("PythiaFilter",
  78. Status = cms.untracked.int32(2),
  79. MaxEta = cms.untracked.double(20.0),
  80. MinEta = cms.untracked.double(-20.0),
  81. MinPt = cms.untracked.double(3.0),
  82. ParticleID = cms.untracked.int32(443)
  83. )
  84.  
  85.  
  86. mumugenfilter = cms.EDFilter("MCParticlePairFilter",
  87. Status = cms.untracked.vint32(1, 1),
  88. MinPt = cms.untracked.vdouble(0.5, 0.5),
  89. MinP = cms.untracked.vdouble(0.,0.),
  90. MaxEta = cms.untracked.vdouble(2.5, 2.5),
  91. MinEta = cms.untracked.vdouble(-2.5, -2.5),
  92. ParticleCharge = cms.untracked.int32(-1),
  93. MaxInvMass = cms.untracked.double(4.0),
  94. MinInvMass = cms.untracked.double(2.0),
  95. ParticleID1 = cms.untracked.vint32(13),
  96. ParticleID2 = cms.untracked.vint32(13)
  97. )
  98.  
  99. ProductionFilterSequence = cms.Sequence(generator*bfilter*jpsifilter*mumugenfilter)
RAW Paste Data