Playing with SPyroot

• brief setup instructions do:

ssh ui02
(now maybe:  atlas001-cern.unige.ch )
zsh
cd /atlas/users/jboyd/june2008/SPyRoot
source setup_unige_ara.sh
cd examples
spyroot MakePlotsPerProdMechn.py

• The code can be seen from CVS browser at this link
• Results - see plots below

SU3

• METcanv.gif:
  

• Meffplot.gif:
  

• NJcanv.gif:
  
• this looks quite likely as we would expect 2 more jets in Gluino Gluino events than SquarkSquark events and 1 more in GluinoSquark than Squark Squark - which is approx what we see...

• Ratios of events with different production mechanisms

GluinoGluino  had  707  events (=  7%)
SquarkSquark  had  3716  events (=  37%)
GauginoGaugino  had  432  events (= 4%)
SleptonSlepton  had  144  events (= 1%)
GluinoSquark  had  4429  events (= 44%)
GluinoGaugino  had  133  events (= 1%)
GluinoSlepton  had  0  events (= 0  %)
SquarkGaugino  had  436  events (=  4%)
SquarkSlepton  had  0  events (= 0  %)
GauginoSlepton  had  0  events (= 0  %)

SU6

• SU6Meffplot.gif:
  

• Ratios of events with different production mechanisms

GluinoGluino  had  687  events (= 6.87  %)
SquarkSquark  had  3471  events (= 34.71  %)
GauginoGaugino  had  852  events (= 8.52  %)
SleptonSlepton  had  102  events (= 1.02  %)
GluinoSquark  had  4341  events (= 43.41  %)
GluinoGaugino  had  157  events (= 1.57  %)
GluinoSlepton  had  0  events (= 0.0  %)
SquarkGaugino  had  390  events (= 3.9  %)
SquarkSlepton  had  0  events (= 0.0  %)
GauginoSlepton  had  0  events (= 0.0  %)

SU4

• SU4Meffplot.gif:
  

• Ratios of events with different production mechanisms

GluinoGluino  had  736  events (= 15.6595744681  %)
SquarkSquark  had  1513  events (= 32.1914893617  %)
GauginoGaugino  had  213  events (= 4.53191489362  %)
SleptonSlepton  had  1  events (= 0.0212765957447  %)
GluinoSquark  had  2082  events (= 44.2978723404  %)
GluinoGaugino  had  35  events (= 0.744680851064  %)
GluinoSlepton  had  0  events (= 0.0  %)
SquarkGaugino  had  119  events (= 2.53191489362  %)
SquarkSlepton  had  0  events (= 0.0  %)
GauginoSlepton  had  0  events (= 0.0  %)

SU1 (only 500 events available it seems...)

• SU1Meffplot.gif:
  
• Ratios of events with different production mechanisms

GluinoGluino  had  39  events (= 7.8  %)
SquarkSquark  had  172  events (= 34.4  %)
GauginoGaugino  had  25  events (= 5.0  %)
SleptonSlepton  had  13  events (= 2.6  %)
GluinoSquark  had  223  events (= 44.6  %)
GluinoGaugino  had  6  events (= 1.2  %)
GluinoSlepton  had  0  events (= 0.0  %)
SquarkGaugino  had  22  events (= 4.4  %)
SquarkSlepton  had  0  events (= 0.0  %)
GauginoSlepton  had  0  events (= 0.0  %)

Comparison of Meff with the true SUSY particle mass

• for this test i ran on a mixed sampl of different SUSY models (including SU1,3,4,6) and then i plot the true mass of the 2 susy particles produced in the hard scatter versus the experimental Meff (calculated as the scalar sum of missingET + Pt of all Cone4H1TopoParticleJets.
• The plot (below) shows good agreement between Meff and the real mass

Issues & to do...

# to do (jamie 5 june 08)
# 1. make it so to add a new variable is easier (just adding name and tree var to a dictionary??) - then add eg. lepton multiplcities, Pt spetra of particles??
# 3 can it run faster somehow (current rate claimed to be 12Hz on 10k evts which is very slow)???

-- JamieBoyd - 05 Jun 2008