SHARE
TWEET

Untitled

a guest Aug 21st, 2019 94 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
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top