TopPartnerAnalysis < Main

<strong>*Running the analysis*</strong>

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

First, download and setup Delphes ([[http://cp3.irmp.ucl.ac.be/downloads/Delphes-3.0.7.tar.gz][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: [[%ATTACHURL%/runAnalysis.C][runAnalysis.C]] and [[%ATTACHURL%/Ttilde.C][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: [[http://brux2.hep.brown.edu/~narain/VQDelphes/][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-&gt;f<em>YourHistName</em>-&gt;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 &lt;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~-&gt;ht):

<img alt="SelectedTtilde_M_1000.png" src="%ATTACHURL%/SelectedTtilde_M_1000.png" />

Reconstructed TTilde mass after selection ( *1500 GeV*, no pileup, T~-&gt;ht)

<img alt="SelectedTtilde_M_1500.png" src="%ATTACHURL%/SelectedTtilde_M_1500.png" />

SVN area : <a target="_blank" href="http://svn.cern.ch/reps/atlas-tandeen">svn.cern.ch/reps/atlas-tandeen</a>/SnowmassAnalysis/

-- Nevis.TimothyAndeen - 2013-06-10

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
C	runAnalysis.C	r1	manage	0.6 K	2013-06-10 - 15:38	TimothyAndeen
png	SelectedTtilde_M_1000.png	r1	manage	10.5 K	2013-06-10 - 15:27	TimothyAndeen	M = 1000 GeV
png	SelectedTtilde_M_1500.png	r1	manage	9.7 K	2013-06-10 - 15:30	TimothyAndeen	TTildeM1500
C	Ttilde.C	r1	manage	19.5 K	2013-06-10 - 15:38	TimothyAndeen

This topic: Main > TWikiUsers > TimothyAndeen > TopPartnerLheFiles > TopPartnerAnalysis
Topic revision: r2 - 2013-06-19 - TimothyAndeen