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.

Jack2.png Jack2c.png

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:

Process sigma(LO) sigma(NLO)
st_1 st_1 2.0 pb 3.1 pb
sb_1 sb_1 2.1 pb 3.4 pb
X1+- X20 32+15 fb 40+18 fb
X1+ X1- 26 fb 33 fb
X10 X20 25 fb 31 fb
X1+- X30 0.11+0.04 fb 0.13 + 0.05 fb
X20 X30 0.07 fb 0.09 fb
gluino gluino 2.1 fb 6.3 fb

Compressed:

Process sigma(LO) sigma(NLO)
st_1 st_1 22 fb 37 fb
sb_1 sb_1 22 fb 38 pb
X1+- X20 26+14 fb 32+18 fb
X1+ X1- 26 fb 32 fb
X10 X20 24 fb 30 fb
X1+- X30 0.13+0.05 fb 0.16 + 0.04 fb
X20 X30 0.09 fb 0.11 fb
gluino gluino 1.6 fb 5.1 fb

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

Topic attachments
I Attachment Action Size Date Who Comment
pngpng Jack2.png manage 47.7 K 2013/06/06 - 14:10 FkW  
elseslha Jack2.slha manage 36.4 K 2013/06/06 - 14:07 FkW  
pngpng Jack2c.png manage 48.3 K 2013/06/06 - 14:10 FkW  
elseslha Jack2c.slha manage 36.6 K 2013/06/06 - 14:06 FkW  
Topic revision: r2 - 2013/07/26 - 01:19:57 - FkW
 
This site is powered by the TWiki collaboration platformCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback