*Running the analysis*

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

Topic attachments
I Attachment History Action Size Date Who Comment
PNGpng SelectedTtilde_M_1000.png r1 manage 10.5 K 2013-06-10 - 15:27 TimothyAndeen M = 1000 GeV
PNGpng SelectedTtilde_M_1500.png r1 manage 9.7 K 2013-06-10 - 15:30 TimothyAndeen TTildeM1500
C source code fileC Ttilde.C r1 manage 19.5 K 2013-06-10 - 15:38 TimothyAndeen  
C source code fileC runAnalysis.C r1 manage 0.6 K 2013-06-10 - 15:38 TimothyAndeen  
Edit | Attach | Watch | Print version | History: r2 < r1 | Backlinks | Raw View | Raw edit | More topic actions
Topic revision: r2 - 2013-06-19 - TimothyAndeen
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2023 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback