| 1 |
algomez |
1.1 |
#!/bin/bash
|
| 2 |
|
|
|
| 3 |
|
|
sample=(200 250 300 350 400 450 500 750 1000 1250 1500 )
|
| 4 |
|
|
crossSection=(502 139 50.2 19.9 8.26 3.70 1.81 0.0816 0.00632 0.000768 0.000185 )
|
| 5 |
|
|
|
| 6 |
|
|
option=(_less _more)
|
| 7 |
|
|
parp67=(2.0 3.0)
|
| 8 |
|
|
parp71=(2.5 8.0)
|
| 9 |
|
|
|
| 10 |
|
|
sample_elements=${#sample[@]}
|
| 11 |
|
|
option_elements=${#option[@]}
|
| 12 |
|
|
for ((j=0;j<$option_elements;j++));
|
| 13 |
|
|
do
|
| 14 |
|
|
for ((i=0;i<$sample_elements;i++));
|
| 15 |
|
|
do
|
| 16 |
|
|
OUTFILE='RPVgluino_M_'${sample[${i}]}'_8TeV_pythia6D6THC_112_cff'${option[${j}]}'_ISRFSR.py'
|
| 17 |
|
|
/bin/rm -f $OUTFILE
|
| 18 |
|
|
echo "import FWCore.ParameterSet.Config as cms
|
| 19 |
|
|
|
| 20 |
|
|
# from Configuration.Generator.PythiaUEZ2Settings_cfi import *
|
| 21 |
|
|
generator = cms.EDFilter(\"Pythia6GeneratorFilter\",
|
| 22 |
|
|
pythiaHepMCVerbosity = cms.untracked.bool(False),
|
| 23 |
|
|
maxEventsToPrint = cms.untracked.int32(0),
|
| 24 |
|
|
pythiaPylistVerbosity = cms.untracked.int32(0),
|
| 25 |
|
|
filterEfficiency = cms.untracked.double(1.0),
|
| 26 |
|
|
comEnergy = cms.double(8000.0),
|
| 27 |
|
|
crossSection = cms.untracked.double(${crossSection[${i}]}),
|
| 28 |
|
|
PythiaParameters = cms.PSet(
|
| 29 |
|
|
pythiaUESettings = cms.vstring('MSTJ(11)=3 ! Choice of the fragmentation function',
|
| 30 |
|
|
'MSTJ(22)=2 ! Decay those unstable particles',
|
| 31 |
|
|
'PARJ(71)=10 . ! for which ctau 10 mm',
|
| 32 |
|
|
'MSTP(2)=1 ! which order running alphaS',
|
| 33 |
|
|
'MSTP(33)=0 ! no K factors in hard cross sections',
|
| 34 |
|
|
'MSTP(51)=10042 ! CTEQ6L1 structure function chosen',
|
| 35 |
|
|
'MSTP(52)=2 ! work with LHAPDF',
|
| 36 |
|
|
'MSTP(81)=1 ! multiple parton interactions 1 is Pythia default',
|
| 37 |
|
|
'MSTP(82)=4 ! Defines the multi-parton model',
|
| 38 |
|
|
'MSTU(21)=1 ! Check on possible errors during program execution',
|
| 39 |
|
|
'PARP(82)=1.8387 ! pt cutoff for multiparton interactions',
|
| 40 |
|
|
'PARP(89)=1960. ! sqrts for which PARP82 is set',
|
| 41 |
|
|
'PARP(83)=0.5 ! Multiple interactions: matter distrbn parameter',
|
| 42 |
|
|
'PARP(84)=0.4 ! Multiple interactions: matter distribution parameter',
|
| 43 |
|
|
'PARP(90)=0.16 ! Multiple interactions: rescaling power',
|
| 44 |
|
|
#'PARP(67)=2.5 ! amount of initial-state radiation',
|
| 45 |
|
|
'PARP(85)=1.0 ! gluon prod. mechanism in MI',
|
| 46 |
|
|
'PARP(86)=1.0 ! gluon prod. mechanism in MI',
|
| 47 |
|
|
'PARP(62)=1.25 ! ',
|
| 48 |
|
|
'PARP(64)=0.2 ! ',
|
| 49 |
|
|
'MSTP(91)=1 !',
|
| 50 |
|
|
'PARP(91)=2.1 ! kt distribution',
|
| 51 |
|
|
'PARP(93)=15.0 ! ',
|
| 52 |
|
|
#ISR/FSR paramters
|
| 53 |
|
|
'PARP(67)=${parp67[${j}]} ! ',
|
| 54 |
|
|
'PARP(71)=${parp71[${j}]} ! '),
|
| 55 |
|
|
processParameters = cms.vstring('MSEL=39',
|
| 56 |
|
|
'IMSS(1)=1',
|
| 57 |
|
|
'IMSS(3)=1',
|
| 58 |
|
|
'IMSS(51)=0',
|
| 59 |
|
|
'IMSS(52)=0',
|
| 60 |
|
|
'IMSS(53)=3',
|
| 61 |
|
|
'RVLAMB(1,1,2) = 0.005',
|
| 62 |
|
|
'RMSS(1)=10000.0',
|
| 63 |
|
|
'RMSS(2)=10000.0',
|
| 64 |
|
|
'RMSS(4)=10000.0',
|
| 65 |
|
|
'RMSS(3)=${sample[${i}]}',
|
| 66 |
|
|
'RMSS(5)=5.00',
|
| 67 |
|
|
'RMSS(8)=2000.0',
|
| 68 |
|
|
'RMSS(9)=2000.0',
|
| 69 |
|
|
'RMSS(10)=2000.0',
|
| 70 |
|
|
'RMSS(11)=2000.0',
|
| 71 |
|
|
'RMSS(12)=2000.0',
|
| 72 |
|
|
'RMSS(6)=10000.0',
|
| 73 |
|
|
'RMSS(7)=10000.0',
|
| 74 |
|
|
'RMSS(13)=10000.0',
|
| 75 |
|
|
'RMSS(14)=10000.0',
|
| 76 |
|
|
'RMSS(15)=2.4',
|
| 77 |
|
|
'RMSS(16)=-3.3',
|
| 78 |
|
|
'RMSS(17)=1.0',
|
| 79 |
|
|
'RMSS(18)=-0.1020',
|
| 80 |
|
|
'RMSS(19)=10000.0'),
|
| 81 |
|
|
parameterSets = cms.vstring('pythiaUESettings',
|
| 82 |
|
|
'processParameters')
|
| 83 |
|
|
)
|
| 84 |
|
|
)" >> $OUTFILE
|
| 85 |
|
|
done
|
| 86 |
|
|
done
|
