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