&gamma&gamma+MET

&gamma&gamma+ME_T Analysis Codes
Published PRL: Phys. Rev. Lett. 104, 011801 (2010)

The following instruction describes how to run &gamma&gamma+ME_T analysis codes. All analysis techniques and physics are well described in my thesis and cdfnote 9575. The first thing you need to do is to install Stntuple package on a local machine. See Stntuple Documentation tiki page for more information on this. Once you have successfully installed Stntuple package then do the step-by-step instruction below:

How to set up and run main analysis codes

The tarball of the analysis code package is located in here. You can download it over your local area.

1. Untar and copy all files in the tarball into Your_Stntuple_release/Stntuple/ana/

2. If you don't have a jet corrections package do this:
>addpakg JetUser jetCorr12
>gmake JetUser.lib LINK_SHARED_ROOT=yes
>cvsroot zoom
>addpkg ZMutility
>addpkg Exceptions
>gmake ZMutility.all USESHLIBS=1
>gmake Exceptions.all USESHLIBS=1

3. Modify rootlogon.C according to instructions on Jet Correction webpage

4. The main script to run all analysis codes is myDiPho_ggXsignal_GMSB_fnl.C. Copy this code over your local analysis area, your_directory

5. In myDiPho_ggXsignal_GMSB_fnl.C replace /home/elee/Stn_614_new/ with an actual path to your release

6. In myDiPho_ggXsignal_GMSB_fnl.C replace all input data/MC stntuple root file locations with an actual location of your files (For 2.6 fb^-1 all data and MC root files are located at txpc15.fnal.gov:/raid2/elee/GMSB/ and GMSB signal sample for (140,0) at txpc15.fnal.gov:/raid2/elee/emtestfinal/0-11/.)

7. For a test run with 10 events, for example, our optimal cuts (H_T>200 GeV, ME_T-significance>3, &Delta&phi(&gamma₁,&gamma₂) < &pi -0.35) do the following:
>root.exe
root>.x myDiPho_ggXsignal_GMSB_fnl.C(1,10)

8. For the entire sample with the preselection cuts (no kinematic cuts):
>root.exe
root>.x myDiPho_ggXsignal_GMSB_fnl.C(0,-1)

Note that the first parameter controls the kinematic cuts and the second does the number of events you want to run over.

9. After running the script you will get small output root files that contains lots of histograms for plots: Exmaple histograms of some kinematic variables for the preselection diphoton data sample are shown below

ME_T ME_T-significance H_T &Delta&phi

10. Here is the example code to produce the ME_T-significance plot for the diphoton presample.

For comments or questions, please e-mail to Eunsin Lee