Full Model Studies at 14TeV in preparation for phase 2 Upgrade
If we were to see an excess in the jets + MET + btags (+ leptons) final state
during Run2 of CMS, then we would almost certainly need all the data we can
get with the phase 2 upgrade, i.e. the nominal 3000/fb, in order to study the excess, and
start understanding its origin.
To get a crude sense of what this may entail, we picked two "full model" spectra
motivated by "natural SUSY" and consistent with a 125GeV Higgs that can be
discovered in Run2.
This page describes these two model points, and provides pointers to
studies on these two points.
The Spectra, SLHA, and LHE files
There are two spectra (
slha compressed,
slha uncompressed), one with
"compressed" stop/sbottom, and one with
"uncompressed" stop/sbottom.
Other than the stop/sbottom mass, the two spectra are essentially identical. The Figures below show the spectra with all transitions that have more than 20% branching fraction.
The lines are darker the larger the branching fraction.
We produced LHE files from these SLHA files using pythia. For each, we produced 100k events to start with. They can be found
here.
The same location also hosts the Pythia logfile, and thus a record of the LO cross section as calculated by Pythia.
The LO and NLO cross sections for all production processes
Prospino2, CTEQ6.6M NLO PDF.
Compressed:
Compressed:
Simple Object multiplicity study
Simple Study of Kinematics at generator level
Getting Started for Milu
Log into uaf-9
ssh mcw@uaf-9.t2.ucsd.edu
copy the basic code you'll need from my directory:
cd Project
cp ~fkw/public_html/LHERead.tgz .
tar -xvzf LHERead.tgz
Now you have the software you need.
Up to hear you have to do only once.
From here on down, you need to do this every time to start a new terminal and log into uaf-9.
setup your environment after ssh into uaf-9
cd ~fkw/CMS/cmsswWork
source source_me
cd CMSSW_4_2_7/src
eval `scramv1 run -sh`
At this point you are all set to start the project for real.
cd LHERead
follow the README, except not quite. Do as follows instead:
python python/readT2bwHE.py 8TeV_T2bw_2j_500_100_30000.lhe test.root
This fails with the following error:
Traceback (most recent call last):
File "python/readT2bwHE.py", line 26, in
myLHEevent.fillEvent(oneEvent)
File "/home/users/fkw/junk/LHERead/python/LHEevent.py", line 22, in fillEvent
self.Particles.append(self.readParticle(lheLines[i]))
File "/home/users/fkw/junk/LHERead/python/LHEevent.py", line 30, in readParticle
return {'ID': int(dataINgood[0]),
ValueError: invalid literal for int() with base 10: '#'
I think the problem is that the event block in the *.lhe file contains two lines per event that start with "#".
E.g.:
18 4 1.000000E+00 5.196250E+02 7.957747E-02 1.018780E-01
21 -1 0 0 101 102 0.0000000000E+00 0.0000000000E+00 5.6189293180E+02 5.6189293180E+02 0.0000000000E+00 0. 9.
21 -1 0 0 102 103 0.0000000000E+00 0.0000000000E+00 -7.0757767881E+02 7.0757767881E+02 0.0000000000E+00 0. 9.
1000006 2 1 2 104 0 6.3010697192E+01 -1.0485284718E+02 -2.3525552390E+02 5.6595903478E+02 5.0000000000E+02 0. 9.
-1000006 2 1 2 0 105 -3.9751046406E+01 1.4021546852E+02 7.9394609393E+01 5.2682447485E+02 5.0000000000E+02 0. 9.
21 1 1 2 105 103 6.3886919639E+00 -4.3762453436E+01 -7.4466465781E+01 8.6609596707E+01 0.0000000000E+00 0. 9.
21 1 1 2 101 104 -2.9648342750E+01 8.3998320982E+00 8.4642633281E+01 9.0077504271E+01 0.0000000000E+00 0. 9.
1000024 2 3 0 0 0 4.1439285836E+01 -3.1080255145E+00 6.2658420434E+00 3.0295153540E+02 3.0002249547E+02 0. 9.
5 1 3 0 104 0 2.1571411356E+01 -1.0174482167E+02 -2.4152136594E+02 2.6300749938E+02 4.8000000000E+00 0. 9.
1000022 1 7 0 0 0 -6.5133332677E+01 -6.0714044295E+01 -5.6021391778E+01 1.4514455738E+02 1.0000000000E+02 0. 9.
24 2 7 0 0 0 1.0657261851E+02 5.7606018781E+01 6.2287233821E+01 1.5780697802E+02 7.9669105666E+01 0. 9.
-15 1 10 0 0 0 4.6253232643E+00 5.2641183257E+00 -1.6986715196E+01 1.8461061161E+01 1.7770000000E+00 0. 9.
16 1 10 0 0 0 1.0194729525E+02 5.2341900455E+01 7.9273949017E+01 1.3934591686E+02 0.0000000000E+00 0. 9.
-1000024 2 4 0 0 0 1.1093583828E+02 4.9945899498E+01 1.2206960048E+02 3.4599492203E+02 3.0001690888E+02 0. 9.
-5 1 4 0 0 105 -1.5068688469E+02 9.0269569022E+01 -4.2674991090E+01 1.8082955282E+02 4.8000000000E+00 0. 9.
1000022 1 13 0 0 0 4.7813983164E+01 -7.8402732720E+01 -8.0278779972E+00 1.3600592748E+02 1.0000000000E+02 0. 9.
-24 2 13 0 0 0 6.3121855120E+01 1.2834863222E+02 1.3009747848E+02 2.0998899455E+02 8.1928529467E+01 0. 9.
1 1 16 0 106 0 1.1353657568E+01 3.1326985510E+01 -7.1886900295E+00 3.4089173144E+01 3.3000000000E-01 0. 9.
-2 1 16 0 0 106 5.1768197552E+01 9.7021646707E+01 1.3728616851E+02 1.7589982141E+02 3.3000000000E-01 0. 9.
# 0.5196250E+03 0.5196250E+03 0.4406665E+02 0.3081528E+02
# model T2bw_500_100_30000 0.0855847
See the last two lines.
The code that reads this and fails is in LHEevent.py here:
def readParticle(self, lheLine):
dataIN = lheLine[:-1].split(" ")
dataINgood = []
for entry in dataIN:
if entry != "": dataINgood.append(entry)
return {'ID': int(dataINgood[0]),
'mIdx': int(dataINgood[2])-1,
'Px' : float(dataINgood[6]),
'Py' : float(dataINgood[7]),
'Pz' : float(dataINgood[8]),
'E' : float(dataINgood[9]),
'M' : float(dataINgood[10])}
It looks to me like it fails when trying to convert the hashtag into an integer.
So you need to put an if in here to avoid the hastag being read at all.
Basically, lines with hastag should be ignored completely.
Good luck.
--
FkW - 2013/06/06