Direct Photon Analysis

This page documents the work on the prompt photon production in the Standard Model group. The current goal is to have a CSC-like note available as soon as possible.

Meetings in indico: Direct Photon Meetings
Mailing list: Direct Photon HyperNews

-- AndrewHamilton - 19 Sep 2007

Interested People

The list of people who have stated an interested in direct photons at ATLAS:

  • Tere Dova (UNLP)
  • Martin Tripiana (UNLP)
  • Andrew Hamilton (Geneva)
  • Francesca Bucci (Geneva)
  • John Wilson (Birmingham)
  • Mark Stockton (Birmingham)
  • Paul Newman (Birmingham)
  • Bob Blair (ANL)
  • Larry Price (ANL)
  • Sergei Chekanov (ANL)
  • Valeria Perez-Reale (CERN)

Photon Note

The note is accessible via the DirectPhoton package in CVS. To edit it do the following (email AndrewHamilton if you want commit rights):

  • set CVSROOT environment variable
    export CVSROOT=:gserver:isscvs.cern.ch:/local/reps/atlas
    (or just set up ATHENA)
  • cvs co -d DIRNAME groups/StandardModel/DirectPhoton/PhotonNote (the DIRNAME directory will be created with content of the CVS PhotonNote directory)
  • commit changes and add files with the usual CVS commands:
    • cvs add -m "blah, blah, message" filename (if file doesn't already exist)
    • cvs update
    • cvs commit -m "blah, blah, message" filename
  • Please turn "Word Wrapping in Text Modes" off - this helps significantly with the readability of the tex file! If you leave it on, then save the file, it will add newlines with whatever width your emacs window happened to be.

For consistency of plots throughout the note, there is a rootlogon.C and a rootalias.C file on CVS. Please put them in your working directory, ROOT will load them by default.

Relevant Datasets

MartinAANs
MarkAODs

Efficiency Plots

Please post any efficiency plots that you feel are relevant for the note. This way we can have a collection of plots to compare and be sure everybody is getting similar results.

kt kick and pt range

Where should we put the lower threshold of the photon pt to avoid (assuming we want to) theoretical uncertainties in the kt of the photon. The link below is a preliminary summary of some work by Fred Olness. Here is a quote from Fred's email: "The short answer is for gamma/Z production the curves merge at about PT=20 for Tevatron, and PT=30 for LHC. What does this mean for direct photons? At Tevatron, the curves merged in the range of PT=30 to 40 GeV. If I scale by the same factor, I'm probably around 60 GeV at LHC. So, if you are interested in the gluon, then you'll need to definitely look above 40 and as high as 60 or more. Working down to lower PT, this will be an excellent test our our ability to really understand (theoretically) the soft physics." olness_kt.pdf: preliminary study by Fred Olness.

Proposed analysis stategy: Use the region PT > 60 for PDF measurements, use region PT < 60 to look at kT problem. Note that we will also need the diphoton measurement to study the kT problem. Joey Huston has written summary of PDF measurements at LHC which looks at this kT problem, although somewhat old, here it is: huston_pdfsAtLHC_1999.pdf.

Another consideration are trigger thresholds, see TriggerPhysicsMenu. The 10^31 menu contains a prescaled g10 and unprescaled g15, g15i, g20, g20i, g25, g105 thresholds. Maybe Martin has some wisdom regarding trigger efficiencies...

Pi0 rejection

Plotted pi0 variables against photon variables (using 20GeV CSC samples): gamVsPi0.ps

The variables plotted are from the PhotonParamDefs.h:

  etcone       =  ET in a cone of R=0.45 in EM calo
  etcone20   =  ET in a cone of R=0.20 in EM calo
  etcone30   =  ET in a cone of R=0.30 in EM calo
  etcone40   =  ET in a cone of R=0.40 in EM calo
  etaBE2       = eta from the second sample
  et37          =  ET in 3x7
  e237         = uncor. energy in 3x7 cells in em sampling 2
  e277         = uncor. energy in 7x7 cells in em sampling 2
  ethad1      =  ET leakage into 1st sampling of had calo
  weta1        =  corrected width in 3 strips in the 1st samp.
  weta2        = corrected width in 3x5 cells in the 2nd samp
  f1              = fraction of energy found in 1st em sampling
  e2tsts1      = 2nd maximum in strips
  emins1      = energy of strip with minimum between max 1 & 2
  wtots1       = total width in em sampling 1 in 20 strips
  fracs1        = energy outside core (E(+-7)-E(+-3))/E(+-7)

Useful Papers

Software to Facilitate Getting High Statistics Background Samples

The package PhotonGenFilter has been checked into CVS (groups/StandardModel/DirectPhoton/PhotonGenFilter). This is a generator filter that finds a cluster of photons and checks to see if they satisfy an isolation cut in both and in the fraction of energy around the cluster. A fragment of python for setting up the filter is here:
#--------------------------------------------------------------
# Filter Options (multi photon filter with isolation for photon bkg. sample)
#--------------------------------------------------------------
GamBkgFilter = Algorithm( "GamBkgFilter" )
GamBkgFilter.RClus = 0.18 # the region of the gamma like cluster
GamBkgFilter.RIso = 0.4  # the size of the isolation cone
GamBkgFilter.IsoCut = 2000. # iso Et cut if Et is less then this in MeV it passes
GamBkgFilter.FracIsoCut = 0.1 # fractional iso cut if fraction is less it passes
GamBkgFilter.FilterMin = 1 # minimum number of candidates
GamBkgFilter.FilterMax = 9999 # maximum number of candidates
GamBkgFilter.CandPtMin = 20000. # minimum pt for a candidate in MeV

#--------------------------------------------------------------
# Pool Output Options
#--------------------------------------------------------------
Stream1 = Algorithm ("Stream1")
Stream1.RequireAlgs += ["GamBkgFilter"]
Note that it passes events that have from "FilterMin" to "FilterMax" candidate photon clusters. In this way it can be used to select events with one or more candidate clusters or as an anti-selection to prevent double counting events when adding a set of filtered events to an unfiltered set (allowing them to be exclusive and permitting weighting of each set by 1/ L).

Topic attachments
I Attachment History Action Size Date Who Comment
Postscriptps Vcdf5721_photon_xsec_thesis.ps r1 manage 4063.9 K 2007-09-27 - 18:10 RobertBlair  
Postscriptps cdf8708_dipho_mass_prl_v2.ps r1 manage 541.6 K 2007-09-18 - 17:30 AndrewHamilton CDF high mass diphotons (with RS limit)
PDFpdf diphoton_cdf.pdf r1 manage 385.9 K 2006-12-09 - 21:32 AndrewHamilton CDF diphoton measurement
Postscriptps gamVsPi0.ps r1 manage 361.0 K 2006-12-12 - 16:35 AndrewHamilton gamma vs. pi0
PDFpdf hepph9808467_kTeffects.pdf r1 manage 408.4 K 2007-09-20 - 09:44 AndrewHamilton Huston et. al. on kT effects on Direct photon production
PDFpdf huston_pdfsAtLHC_1999.pdf r1 manage 1721.7 K 2007-10-31 - 10:09 AndrewHamilton Huston PDFs at LHC
PDFpdf olness_kt.pdf r1 manage 394.6 K 2007-10-03 - 09:47 AndrewHamilton  
Compressed Zip archivetgz photon_note.tgz r2 r1 manage 145.1 K 2007-09-20 - 10:01 AndrewHamilton  
PDFpdf yanwen_thesis.pdf r1 manage 3715.6 K 2007-04-19 - 17:36 AndrewHamilton Yanwen Liu's CDF diphoton thesis
Edit | Attach | Watch | Print version | History: r32 < r31 < r30 < r29 < r28 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r32 - 2009-02-25 - AndrewHamilton
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback