High Pt Photon ID

Introduction

The High pT photon ID was designed using the FLASHgg framework, with the following code https://github.com/cms-analysis/flashgg.

The High pT ID version 1 (V1) and version 2 (V2) were proposed during EXO diphoton meetings.

• V1: https://indico.cern.ch/event/397481/
• V2: https://indico.cern.ch/event/448912/
• V2: https://indico.cern.ch/event/533798/ (updated kappa value)

The most up-to-date version of the High pT ID is V2.

The V2 working points are documented in detail in AN-15-241 (Section 6), but will be repeated here. Different cuts are applied for EB and EE.

ECAL Geometry

The EB and EE cuts are defined with respect to the detector's coordinates, in particular, using "detector" eta (eta_det) or, equivalently, "supercluster" eta (eta_sc). This is in contrast to using the photon's eta coordinate. The cuts are defined as follows

EB: std::abs(eta_sc) < 1.4442

EE: 1.566 < std::abs(eta_sc) < 2.5

We omit the crack between EB and EE, and we do not extend to the crystal at 1.560.

Implementing the High pT ID

We are trying to reproduce what's being done in the FLASHgg framework here.

Photons

This ID was implemented using pat::Photon objects, rather than reco::Photon. Some differences have been noticed between the two. For miniAOD, we use the tag "slimmedPhotons".

H/E

We use the single HCAL tower behind the photon's ECAL energy deposit, from photon->hadTowOverEm().

Charged hadron isolation

We are calling photon->chargedHadronIso(). Note, this ID variable is not pileup (rho) corrected. Should be using a cone of 0.3.

Neutral hadron isolation

No neutral hadron isolation selections are applied.

Photon isolation

We are calling photon->photonIso(). Should be using a cone of 0.3.

Rho

Using "fixedGridRhoAll".

Note: Egamma uses "fixedGridRhoFastjetAll".

"Corrected" photon isolation

The photon isolation is pileup corrected. Instead of putting photon pT dependence in the cut, it is added to the photon isolation as well, further correcting it. The variable is defined as

corPhoIso = alpha + photon->photonIso() - rho*effectiveArea - kappa*photon->pt()

where alpha, kappa, and effectiveArea are numbers defined in the AN.

Sigma_IetaIeta

The full supercluster's 5x5 contribution is used, given by photon->full5x5_sigmaIetaIeta(). This variable's cut changes depending on whether or not the photon is saturated.

Electron veto

The CSEV is applied using photon->passElectronVeto().

Eta

Use detector coordinates, rather than the coordinates of the photon itself. Call photon->superCluster()->eta().

Photo Saturation

The High pT ID was tuned to take into account the affects of saturation. If the photon's seed crystal, or any cell in the 5x5 surrounding the seed with the seed at the center, is marked saturated, then the photon is considered saturated.

High pT ID Cuts

See the AN, for now.

Our code

See PhotonID.h for our implementation. https://github.com/cms-exotica-diphotons/diphoton-analysis/blob/master/CommonClasses/interface/PhotonID.h

Cross check with FLASHgg

UA

https://indico.cern.ch/event/507650/

https://indico.cern.ch/event/520444/ (latest)

Caltech

https://indico.cern.ch/event/525445/

https://indico.cern.ch/event/532212/

https://indico.cern.ch/event/533798/ (EXO-NH Meeting)

https://indico.cern.ch/event/537100/

https://indico.cern.ch/event/537578/

https://indico.cern.ch/event/537546/ (EXO Meeting)

https://indico.cern.ch/event/539613/

https://indico.cern.ch/event/540477/ (EXO Meeting)

https://indico.cern.ch/event/542288/

https://indico.cern.ch/event/542711/ (latest)

-- AndrewBuccilli - 2016-03-07