Below is a work in progress, but I wanted to make it available.

First, download and setup Delphes (Delphes 3.0.7), 3.0.9 (the most recent) works fine also. Install it following the instructions on the delphes website.

Then copy the files attached here: runAnalysis.C and Ttilde.C, to the Delphes-3.0.7/ directory. runAnalysis is short and controls the input, so figure out where you have copied the root files (ie files like Ts_ht_unweighted_events_33TEV_L3_M1500_*_NoPileUp.root, which you can get here: Input files) and edit the input string in runAnalysis.C. It will loop over twenty files for each mass point automatically, but of course you can change this.

Ttilde.C is based on the example analysis included in the delphes examples/ directory. It takes the reconstructed objets and attempts to do a single lepton (muon, electron) analysis as in arXiv:1207.0830. The functions fitAlpha, candidatesFromWMass_Rotation and delta2_fcn are used to solve for the best neutrino pz momentum, using the W mass. This requires TMinuit, and I have ROOT version 5.34.07. Nearly all the cuts described in the paper are applied and buried in the code is a switch to select the cuts tuned for 1000 GeV Ttilde mass and 1500 GeV Ttilde mass, which should be extended. It should be clear how histograms are created and filled (definied in BookHistograms, filled with plots->fYourHistName->Fill(value)), and much should be added here.

You can run it as it is with:

root runAnalysis.C

There will be a couple of warnings, but no errors and on a fast machine I can run over ~200k events in <1 minute. The output will be several png files (such as below) and a root file with the histograms.

Sample outout:

Reconstructed TTilde mass after selection ( 1000 GeV, no pileup, T~->ht):

Reconstructed TTilde mass after selection ( 1500 GeV, no pileup, T~->ht)

SVN area : svn.cern.ch/reps/atlas-tandeen/SnowmassAnalysis/

-- TimothyAndeen - 2013-06-10