Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2022-08-12 01:59:51

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 from Configuration.Generator.PythiaUESettings_cfi import *
0004 
0005 collisionParametersRHIC200GeV = cms.PSet(
0006         fAw             = cms.double(197.0),    # beam/target atomic number
0007         fSqrtS          = cms.double(200.0),    # 
0008     fMuB        = cms.double(0.0285),   # Chemical baryon potential per unit charge, fMuB [GeV]
0009     fMuS        = cms.double(0.007),    # Chemical strangeness potential per unit charge, fMuS [GeV]  
0010     fMuI3       = cms.double(-0.001),   # Chemical isospin potential per unit charge, fMuI3 [GeV]   
0011 
0012     ### Thermodinamic parameters at thermal freez-out ###
0013     fThFO       = cms.double(0.1),  # Temperature at thermal freeze-out, fTthFO [GeV]
0014     fMu_th_pip  = cms.double(0.06),     # Chemical potential of pi+ at thermal freeze-out, fMu_th_pip [GeV] 
0015 
0016         ### Maximal longitudinal flow rapidity at thermal freeze-out ###
0017         fYlmax          = cms.double(3.3),      # Maximal longitudinal flow rapidity at thermal freeze-out, fYlmax
0018         fUmax           = cms.double(1.1),      # Maximal transverse flow rapidity at thermal freeze-out for central collisions, fUmax
0019         fPtmin          = cms.double(3.55),     # Minimal pt of parton-parton scattering in PYTHIA event, fPtmin [GeV/c] 
0020         fT0             = cms.double(0.3),      # Initial QGP temperature for central Pb+Pb collisions in mid-rapidity, fT0 [GeV]; allowed range [0.2,2.0]GeV;
0021 
0022         ### Volume parameters at thermal freeze-out ###
0023         fTau            = cms.double(8.),       # Proper time proper at thermal freeze-out for central collisions, fTau [fm/c]
0024         fR              = cms.double(10.),      # Maximal transverse radius at thermal freeze-out for central collisions, fR [fm]
0025         fSigmaTau       = cms.double(2.),       # Duration of emission at thermal freeze-out for central collisions, fSigmaTau [fm/c]
0026 
0027     fWeakDecay  = cms.double(0.),   # Low decay width threshold fWeakDecay[GeV]: width<fWeakDecay decay off, width>=fDecayWidth decay on; can be used to switch off weak decays
0028 )
0029 
0030 collisionParameters2760GeV = cms.PSet(
0031     fAw         = cms.double(208.0),    # beam/target atomic number
0032     fSqrtS      = cms.double(2760.0),   #
0033     fMuB        = cms.double(0.),   # Chemical baryon potential per unit charge, fMuB [GeV]
0034     fMuS        = cms.double(0.),   # Chemical strangeness potential per unit charge, fMuS [GeV]  
0035     fMuI3       = cms.double(0.),   # Chemical isospin potential per unit charge, fMuI3 [GeV]   
0036 
0037     ### Thermodinamic parameters at thermal freez-out ###
0038     fThFO       = cms.double(0.105),    # Temperature at thermal freeze-out, fTthFO [GeV]
0039     fMu_th_pip  = cms.double(0.),   # Chemical potential of pi+ at thermal freeze-out, fMu_th_pip [GeV] 
0040 
0041         ### Maximal longitudinal flow rapidity at thermal freeze-out ###
0042         fYlmax          = cms.double(4.5),      # Maximal longitudinal flow rapidity at thermal freeze-out, fYlmax
0043 
0044     fUmax       = cms.double(1.265),    # Maximal transverse flow rapidity at thermal freeze-out for central collisions, fUmax
0045     fPtmin      = cms.double(8.2),  # Minimal pt of parton-parton scattering in PYTHIA event, fPtmin [GeV/c] 
0046     fT0         = cms.double(1.),   # Initial QGP temperature for central Pb+Pb collisions in mid-rapidity, fT0 [GeV]; allowed range [0.2,2.0]GeV;
0047 
0048     ### Volume parameters at thermal freeze-out ###
0049     fTau        = cms.double(12.2),     # Proper time proper at thermal freeze-out for central collisions, fTau [fm/c]
0050     fR      = cms.double(13.45),    # Maximal transverse radius at thermal freeze-out for central collisions, fR [fm]
0051     fSigmaTau   = cms.double(3.5),  # Duration of emission at thermal freeze-out for central collisions, fSigmaTau [fm/c]
0052 
0053     fWeakDecay  = cms.double(0.00000000000001), # Low decay width threshold fWeakDecay[GeV]: width<fWeakDecay decay off, width>=fDecayWidth decay on; can be used to switch off weak decays
0054 
0055 )
0056 
0057 collisionParameters5020GeV = cms.PSet(
0058     fAw         = cms.double(208.0),    # beam/target atomic number
0059     fSqrtS      = cms.double(5020.0),
0060     fMuB        = cms.double(0.),   # Chemical baryon potential per unit charge, fMuB [GeV]
0061     fMuS        = cms.double(0.),   # Chemical strangeness potential per unit charge, fMuS [GeV]  
0062     fMuI3       = cms.double(0.),   # Chemical isospin potential per unit charge, fMuI3 [GeV]   
0063 
0064     ### Thermodinamic parameters at thermal freez-out ###
0065     fThFO       = cms.double(0.105),    # Temperature at thermal freeze-out, fTthFO [GeV]
0066     fMu_th_pip  = cms.double(0.),   # Chemical potential of pi+ at thermal freeze-out, fMu_th_pip [GeV] 
0067 
0068         ### Maximal longitudinal flow rapidity at thermal freeze-out ###
0069         fYlmax          = cms.double(3.99),      # Maximal longitudinal flow rapidity at thermal freeze-out, fYlmax
0070 
0071 
0072     fUmax       = cms.double(1.280),    # Maximal transverse flow rapidity at thermal freeze-out for central collisions, fUmax
0073     fPtmin      = cms.double(9.06),     # Minimal pt of parton-parton scattering in PYTHIA event, fPtmin [GeV/c] 
0074     fT0         = cms.double(1.1),  # Initial QGP temperature for central Pb+Pb collisions in mid-rapidity, fT0 [GeV]; allowed range [0.2,2.0]GeV;
0075 
0076     ### Volume parameters at thermal freeze-out ###
0077     fTau        = cms.double(11.5),     # Proper time proper at thermal freeze-out for central collisions, fTau [fm/c]
0078     fR      = cms.double(16.),  # Maximal transverse radius at thermal freeze-out for central collisions, fR [fm]
0079     fSigmaTau   = cms.double(2.),   # Duration of emission at thermal freeze-out for central collisions, fSigmaTau [fm/c]
0080 
0081     fWeakDecay  = cms.double(0.00000000000001), # Low decay width threshold fWeakDecay[GeV]: width<fWeakDecay decay off, width>=fDecayWidth decay on; can be used to switch off weak decays
0082 )
0083 
0084 qgpParametersLHC = cms.PSet(
0085     fTau0       = cms.double(0.1),  # Proper QGP formation time in fm/c, fTau0 (0.01<fTau0<10)
0086     fNf         = cms.int32(0),     # Number of active quark flavours in QGP, fNf (0, 1, 2 or 3)
0087 )
0088 
0089 qgpParametersRHIC = cms.PSet(
0090         fTau0           = cms.double(0.4),      # Proper QGP formation time in fm/c, fTau0 (0.01<fTau0<10)
0091         fNf             = cms.int32(2),         # Number of active quark flavours in QGP, fNf (0, 1, 2 or 3)
0092 )
0093 
0094 
0095 hydjet2Parameters = cms.PSet(
0096     ### Thermodinamic parameters at chemical freez-out ###
0097     fTMuType    = cms.double(0.),   # Flag to use calculated T_ch, mu_B and mu_S as a function of fSqrtS, fTMuType (=0 user's ones, >0 calculated) 
0098     fT      = cms.double(0.165),    # Temperature at chemical freeze-out, fT [GeV]
0099     fMuC        = cms.double(0.),   # Chemical charm potential per unit charge, fMuC [GeV] (used if charm production is turned on)  
0100     ### Strangeness suppression factor ###
0101     fCorrS      = cms.double(1.),   # Strangeness supression factor gamma_s with fCorrS value (0<fCorrS <=1, if fCorrS <= 0 then it is calculated)
0102 
0103 
0104         fRhou2      = cms.double(0.),   # Parameter to specify/switch ON(>0)/OFF(0) elliptic modulations of collective velocity profile for the alternative treatment of anisotropic flow 
0105     
0106     fRhou3          = cms.double(0.),        # Parameter to specify/switch ON(>0)/OFF(0) triangular modulations of collective velocity profile for the alternative treatment of anisotropic flow
0107 
0108     fRhou4          = cms.double(0.),        # Parameter to specify/switch ON(>0)/OFF(0) quadrangular modulations of collective velocity profile for the alternative treatment of anisotropic flow
0109 
0110     ### Anizotropy parameter at thermal freeze-out ###
0111     fIfDeltaEpsilon = cms.double(1.),   # Flag to specify fDelta and fEpsilon values, fIfDeltaEpsilon (=0 user's ones, >=1 calculated) 
0112     fDelta      = cms.double(0.1),  # Momentum azimuthal anizotropy parameter at thermal freeze-out, fDelta
0113     fEpsilon    = cms.double(0.05),     # Spatial azimuthal anisotropy parameter at thermal freeze-out, fEpsilon
0114 
0115     fKeps2      = cms.double(0.34), # Parameter to specify/switch ON(>0)/OFF(0) elliptic flow fluctuations
0116 
0117         fKeps3          = cms.double(0.52),     # Parameter to specify/switch ON(>0)/OFF(0) triangular flow fluctuations
0118 
0119     ### Decays ###
0120     fDecay      = cms.int32(1),     # Flag to switch on/off hadron decays, fDecay (=0 decays off, >=1 decays on)
0121     
0122     ### Charm ###
0123     fCharmProd  = cms.int32(1),     # Flag to include thermal charm production, fCharmProd (=0 no charm production, >=1 charm production) 
0124     fCorrC      = cms.double(-1.),  # Charmness enhancement factor gamma_c with fCorrC value (fCorrC >0, if fCorrC<0 then it is calculated)
0125     fEtaType    = cms.double(1.),   # Flag to choose longitudinal flow rapidity distribution, fEtaType (=0 uniform, >0 Gaussian with the dispersion Ylmax)
0126     fIshad      = cms.int32(1),     # Flag to switch on/off nuclear shadowing, fIshad (0 shadowing off, 1 shadowing on) 
0127     fPyhist     = cms.int32(0),     # Flag to suppress the output of particle history from PYTHIA, fPyhist (=1 only final state particles; =0 full particle history from PYTHIA)
0128     fIenglu     = cms.int32(0),     # Flag to fix type of partonic energy loss, fIenglu (0 radiative and collisional loss, 1 radiative loss only, 2 collisional loss only)
0129     fIanglu     = cms.int32(0),     # Flag to fix type of angular distribution of in-medium emitted gluons, fIanglu (0 small-angular, 1 wide-angular, 2 collinear).
0130     embeddingMode   = cms.int32(0),
0131     rotateEventPlane= cms.bool(True)
0132 )
0133 
0134 PythiaDefaultBlock = cms.PSet(
0135     pythiaUESettingsBlock,
0136     TDB = cms.vstring(
0137         'PARJ(14)=0.'   # ! replacing the same parametr from pythiaQuarkoniaSettings block to avoid producing h'_1 (pdg - 10333)
0138     ),
0139     hydjet2PythiaDefault = cms.vstring(
0140         'MSEL=1',   # ! type of hard QCD production process
0141         'MSTU(21) = 1', # ! controle parameter to avoid stopping run
0142         'PARU(14)=1.',  # ! tolerance parameter to adjust fragmentation'
0143         'MSTP(81)=1',   # ! pp multiple scattering on (UE model)
0144         'MSTJ(21) = 1', # ! hadron decays on (if off - decays by FASTMC decayer) 
0145         'MSTP(2) = 1',  # ! which order running alphaS 
0146         'MSTP(33) = 0', # ! inclusion of k factor in cross section (on/off)
0147 
0148         #Mod for proQ20
0149                 'parp(67)=1.',  # ! ISR Q2max factor (amount of initial-state radiation)                            
0150                 'parp(82)=2.',  # ! UE IR cutoff at reference ecm    
0151                 'mstj(11)=3',   # ! HAD choice of fragmentation function(s)      
0152 
0153         #'MSTJ(22)=2',  # ! particle decays if lifetime < parj(71)
0154         #'PARJ(71)=10.',# ! ctau=10 mm 
0155         #'MSTP(52) = 1',# ! NO LAPDF        
0156         #'mstp(122)=0'  # ! no printout of Pythia initialization information hereinafter 
0157     ),
0158     ProQ2Otune = cms.vstring(  
0159         'mstp(51)=7',   # ! PDF set: structure function chosen - CTEQ5M pdf                                      
0160             'mstp(3)=2',    # ! QCD switch for choice of LambdaQCD           
0161             'parp(62)=2.9', # ! ISR IR cutoff                                
0162             'parp(64)=0.14',# ! ISR renormalization scale prefactor          
0163             #'parp(67)=2.65',# ! ISR Q2max factor                             
0164             'mstp(68)=3',   # ! ISR phase space choice & ME corrections      
0165             'parp(71)=4.',  # ! FSR Q2max factor for non-s-channel procs     
0166             'parj(81)=0.29',# ! FSR Lambda_QCD scale                         
0167             'parj(82)=1.65',# ! FSR IR cutoff                                
0168             'mstp(33)=0',   # ! "K" switch for K-factor on/off & type        
0169             'mstp(81)=1',   # ! UE model                                     
0170             #'parp(82)=1.9',# ! UE IR cutoff at reference ecm                
0171             'parp(89)=1800.',# ! UE IR cutoff reference ecm                   
0172             'parp(90)=0.22',# ! UE IR cutoff ecm scaling power               
0173             'mstp(82)=4',   # ! UE hadron transverse mass distribution       
0174             'parp(83)=0.83',# ! UE mass distribution parameter               
0175             'parp(84)=0.6', # ! UE mass distribution parameter              
0176             'parp(85)=0.86',# ! UE gg colour correlated fraction             
0177             'parp(86)=0.93',# ! UE total gg fraction                         
0178             'mstp(91)=1',   # ! BR primordial kT distribution                
0179             'parp(91)=2.1', # ! BR primordial kT width <|kT|>                
0180             'parp(93)=5.',  # ! BR primordial kT UV cutoff               
0181             #'mstj(11)=5',  # ! HAD choice of fragmentation function(s)      
0182             'parj(1)=0.073',# ! HAD diquark suppression                      
0183             'parj(2)=0.2',  # ! HAD strangeness suppression                  
0184             'parj(3)=0.94', # ! HAD strange diquark suppression              
0185             'parj(4)=0.032',# ! HAD vector diquark suppression               
0186             'parj(11)=0.31',# ! HAD P(vector meson), u and d only            
0187             'parj(12)=0.4', # ! HAD P(vector meson), contains s              
0188             'parj(13)=0.54',# ! HAD P(vector meson), heavy quarks            
0189             'parj(21)=0.325',# ! HAD fragmentation pT                         
0190             'parj(25)=0.63',# ! HAD eta0 suppression                        
0191             'parj(26)=0.12',# ! HAD eta0' suppression                       
0192             'parj(41)=0.5', # ! HAD string parameter a                       
0193             'parj(42)=0.6', # ! HAD string parameter b                       
0194             'parj(46)=1.',  # ! HAD Lund(=0)-Bowler(=1) rQ (rc)              
0195             'parj(47)=0.67' # ! HAD Lund(=0)-Bowler(=1) rb            
0196 
0197     ),
0198     ppJets = cms.vstring('MSEL=1'),     # ! QCD hight pT processes
0199     customProcesses = cms.vstring('MSEL=0'),# ! User processes
0200     pythiaJets = cms.vstring(
0201         'MSUB(11)=1',   # ! q+q->q+q
0202         'MSUB(12)=1',   # ! q+qbar->q+qbar
0203         'MSUB(13)=1',   # ! q+qbar->g+g
0204         'MSUB(28)=1',   # ! q+g->q+g
0205         'MSUB(53)=1',   # ! g+g->q+qbar
0206         'MSUB(68)=1'    # ! g+g->g+g
0207     ),
0208     pythiaPromptPhotons = cms.vstring(
0209         'MSUB(14)=1',   # ! q+qbar->g+gamma
0210         'MSUB(18)=1',   # ! q+qbar->gamma+gamma
0211         'MSUB(29)=1',   # ! q+g->q+gamma
0212         'MSUB(114)=1',  # ! g+g->gamma+gamma
0213         'MSUB(115)=1'   # ! g+g->g+gamma
0214     ),
0215     pythiaWeakBosons = cms.vstring(
0216         'MSUB(1)=1',
0217         'MSUB(2)=1'
0218     ),
0219     pythiaZjets = cms.vstring(
0220         'MSUB(15)=1',
0221         'MSUB(30)=1'
0222     ),
0223     pythiaCharmoniumNRQCD = cms.vstring(
0224         'MSUB(421) = 1',
0225         'MSUB(422) = 1',
0226         'MSUB(423) = 1',
0227         'MSUB(424) = 1',
0228         'MSUB(425) = 1',
0229         'MSUB(426) = 1',
0230         'MSUB(427) = 1',
0231         'MSUB(428) = 1',
0232         'MSUB(429) = 1',
0233         'MSUB(430) = 1',
0234         'MSUB(431) = 1',
0235         'MSUB(432) = 1',
0236         'MSUB(433) = 1',
0237         'MSUB(434) = 1',
0238         'MSUB(435) = 1',
0239         'MSUB(436) = 1',
0240         'MSUB(437) = 1',
0241         'MSUB(438) = 1',
0242         'MSUB(439) = 1'
0243     ),
0244     pythiaBottomoniumNRQCD = cms.vstring(
0245         'MSUB(461) = 1',
0246         'MSUB(462) = 1',
0247         'MSUB(463) = 1',
0248         'MSUB(464) = 1',
0249         'MSUB(465) = 1',
0250         'MSUB(466) = 1',
0251         'MSUB(467) = 1',
0252         'MSUB(468) = 1',
0253         'MSUB(469) = 1',
0254         'MSUB(470) = 1',
0255         'MSUB(471) = 1',
0256         'MSUB(472) = 1',
0257         'MSUB(473) = 1',
0258         'MSUB(474) = 1',
0259         'MSUB(475) = 1',
0260         'MSUB(476) = 1',
0261         'MSUB(477) = 1',
0262         'MSUB(478) = 1',
0263         'MSUB(479) = 1',
0264     ),
0265     pythiaQuarkoniaSettings = cms.vstring(
0266         'PARP(141)=1.16', # Matrix Elements
0267         'PARP(142)=0.0119',
0268         'PARP(143)=0.01',
0269         'PARP(144)=0.01',
0270         'PARP(145)=0.05',
0271         'PARP(146)=9.28',
0272         'PARP(147)=0.15',
0273         'PARP(148)=0.02',
0274         'PARP(149)=0.02',
0275         'PARP(150)=0.085',
0276         # Meson spin
0277         'PARJ(13)=0.60',
0278         'PARJ(14)=0.162',
0279         'PARJ(15)=0.018',
0280         'PARJ(16)=0.054',
0281         # Polarization
0282         'MSTP(145)=0',
0283         'MSTP(146)=0',
0284         'MSTP(147)=0',
0285         'MSTP(148)=1',
0286         'MSTP(149)=1',
0287         # Chi_c branching ratios
0288         'BRAT(861)=0.202',
0289         'BRAT(862)=0.798',
0290         'BRAT(1501)=0.013',
0291         'BRAT(1502)=0.987',
0292         'BRAT(1555)=0.356',
0293         'BRAT(1556)=0.644'
0294     ),
0295     pythiaZtoMuons = cms.vstring(
0296         "MDME(174,1)=0", # !Z decay into d dbar,
0297         "MDME(175,1)=0", # !Z decay into u ubar,
0298         "MDME(176,1)=0", # !Z decay into s sbar,
0299         "MDME(177,1)=0", # !Z decay into c cbar,
0300         "MDME(178,1)=0", # !Z decay into b bbar,
0301         "MDME(179,1)=0", # !Z decay into t tbar,
0302         "MDME(182,1)=0", # !Z decay into e- e+,
0303         "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
0304         "MDME(184,1)=1", # !Z decay into mu- mu+,
0305         "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
0306         "MDME(186,1)=0", # !Z decay into tau- tau+,
0307         "MDME(187,1)=0"  # !Z decay into nu_tau nu_taubar
0308     ),
0309     pythiaZtoElectrons = cms.vstring(
0310         "MDME(174,1)=0", # !Z decay into d dbar,
0311         "MDME(175,1)=0", # !Z decay into u ubar,
0312         "MDME(176,1)=0", # !Z decay into s sbar,
0313         "MDME(177,1)=0", # !Z decay into c cbar,
0314         "MDME(178,1)=0", # !Z decay into b bbar,
0315         "MDME(179,1)=0", # !Z decay into t tbar,
0316         "MDME(182,1)=1", # !Z decay into e- e+,
0317         "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
0318         "MDME(184,1)=0", # !Z decay into mu- mu+,
0319         "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
0320         "MDME(186,1)=0", # !Z decay into tau- tau+,
0321         "MDME(187,1)=0"  # !Z decay into nu_tau nu_taubar
0322     ),
0323     pythiaZtoMuonsAndElectrons = cms.vstring(
0324         "MDME(174,1)=0", # !Z decay into d dbar,
0325         "MDME(175,1)=0", # !Z decay into u ubar,
0326         "MDME(176,1)=0", # !Z decay into s sbar,
0327         "MDME(177,1)=0", # !Z decay into c cbar,
0328         "MDME(178,1)=0", # !Z decay into b bbar,
0329         "MDME(179,1)=0", # !Z decay into t tbar,
0330         "MDME(182,1)=1", # !Z decay into e- e+,
0331         "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
0332         "MDME(184,1)=1", # !Z decay into mu- mu+,
0333         "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
0334         "MDME(186,1)=0", # !Z decay into tau- tau+,
0335         "MDME(187,1)=0"  # !Z decay into nu_tau nu_taubar
0336     ),
0337     pythiaUpsilonToMuons = cms.vstring(
0338         'BRAT(1034) = 0 ', # switch off',
0339         'BRAT(1035) = 1 ', # switch on',
0340         'BRAT(1036) = 0 ', # switch off',
0341         'BRAT(1037) = 0 ', # switch off',
0342         'BRAT(1038) = 0 ', # switch off',
0343         'BRAT(1039) = 0 ', # switch off',
0344         'BRAT(1040) = 0 ', # switch off',
0345         'BRAT(1041) = 0 ', # switch off',
0346         'BRAT(1042) = 0 ', # switch off',
0347         'MDME(1034,1) = 0 ', # switch off',
0348         'MDME(1035,1) = 1 ', # switch on',
0349         'MDME(1036,1) = 0 ', # switch off',
0350         'MDME(1037,1) = 0 ', # switch off',
0351         'MDME(1038,1) = 0 ', # switch off',
0352         'MDME(1039,1) = 0 ', # switch off',
0353         'MDME(1040,1) = 0 ', # switch off',
0354         'MDME(1041,1) = 0 ', # switch off',
0355         'MDME(1042,1) = 0 ', # switch off'
0356     ),
0357     pythiaJpsiToMuons = cms.vstring(
0358         'BRAT(858) = 0 ', # switch off',
0359         'BRAT(859) = 1 ', # switch on',
0360         'BRAT(860) = 0 ', # switch off',
0361         'MDME(858,1) = 0 ', # switch off',
0362         'MDME(859,1) = 1 ', # switch on',
0363         'MDME(860,1) = 0 ', # switch off'
0364     ),
0365     pythiaBToJpsi = cms.vstring(
0366         'BRAT(889) = 0 ', # switch off B0->J/Psi K0',
0367         'BRAT(890) = 0 ', # switch off B0->J/Psi K*0',
0368         'BRAT(934) = 0 ', # switch off B+->J/Psi K+',
0369         'BRAT(935) = 0 ', # switch off B+->J/Psi K*+',
0370         'BRAT(980) = 0 ', # switch off B_s0->J/Psi eta',
0371         'BRAT(981) = 0 ', # switch off B_s0->J/Psi eta"',
0372         'BRAT(982) = 0 ', # switch off B_s0->J/Psi phi',
0373         'BRAT(1001) = 0 ', # switch off B_c+>J/Psi nu_e e+',
0374         'BRAT(1003) = 0 ', # switch off B_c+>J/Psi nu_mu mu+',
0375         'BRAT(1005) = 0 ', # switch off B_c+>J/Psi nu_tau tau+',
0376 
0377         'MDME(889,1) = 0 ', # switch off',
0378         'MDME(890,1) = 0 ', # switch off',
0379         'MDME(934,1) = 0 ', # switch off'
0380         'MDME(935,1) = 0 ', # switch off'
0381         'MDME(980,1) = 0 ', # switch off'
0382         'MDME(981,1) = 0 ', # switch off'
0383         'MDME(982,1) = 0 ', # switch off'
0384         'MDME(1001,1) = 0 ', # switch off'
0385         'MDME(1003,1) = 0 ', # switch off'
0386         'MDME(1005,1) = 0 ', # switch off'
0387     ),
0388     pythiaXToJpsi = cms.vstring(
0389         'BRAT(1228) = 0 ', # switch off Lambda_b0->J/Psi Lambda0',
0390         'BRAT(1501) = 0 ', # switch off chi_0c->J/Psi gamma',
0391         'BRAT(1555) = 0 ', # switch off chi_1c->J/Psi gamma',
0392         'BRAT(1570) = 0 ', # switch off psi"->J/Psi pi+ pi-',
0393         'BRAT(1571) = 0 ', # switch off psi"->J/Psi pi0 pi0',
0394         'BRAT(1572) = 0 ', # switch off psi"->J/Psi eta',
0395         'BRAT(1573) = 0 ', # switch off psi"->J/Psi pi0',
0396 
0397         'MDME(1228,1) = 0 ', # switch off'
0398         'MDME(1501,1) = 0 ', # switch off'
0399         'MDME(1555,1) = 0 ', # switch off'
0400         'MDME(1570,1) = 0 ', # switch off'
0401         'MDME(1571,1) = 0 ', # switch off'
0402         'MDME(1572,1) = 0 ', # switch off'
0403         'MDME(1573,1) = 0 ', # switch off'
0404     ),
0405     pythiaPromptJpsi = cms.vstring(
0406         'BRAT(4285) = 1 ', # switch off cc~[3S18]->J/Psi g',
0407         'BRAT(4286) = 1 ', # switch off cc~[1S08]->J/Psi g',
0408         'BRAT(4287) = 1 ', # switch off cc~[3P08]->J/Psi g',
0409 
0410         'MDME(4285,1) = 1 ', # switch off'
0411         'MDME(4286,1) = 1 ', # switch off'
0412         'MDME(4287,1) = 1 ', # switch off'
0413     ),
0414     pythiaMuonCandidates = cms.vstring(
0415         'CKIN(3)=20',
0416         'MSTJ(22)=2',
0417         'PARJ(71)=40.'
0418     ),
0419     myParameters = cms.vstring('MDCY(310,1)=0')
0420 )