Heavy Ion Photon reconstruction

Complete: 5

Goal of this page

The aim of this page is to document the photon reconstruction used in heavy ion events, where the code lives, how to run it, the input and output, and details of the implementation.

Code and tags

The code lives in the package RecoHI/HiEgammaAlgos.

The most recent version of the code will be tagged in CMSSW_3_4_0_pre4.

Before CMSSW_3_4_x release, if you want to check out the latest working version yourself, you can check out the release by the following command:

scramv1 project CMSSW CMSSW_3_3_0
cd CMSSW_3_3_0/src
cvs co -r V00-00-06 RecoHI/HiEgammaAlgos

and compile:

scramv1 b

This will prepare the working environment for you.

Configuration files

Compared to the p+p photon reconstruction algorithm, the photon reconstruction in heavy ion collision use Island SuperCluster algorithms for both barrel and endcap regions while in p+p one uses the hybrid algorithm for the barrel and multi5x5 algorithm for the endcap. This choice is found to cache less background contribution compared to hybrid and multi5x5 algorithm.

The reconstruction sequence is defined in RecoHI/HiEgammaAlgos/python/HiEgamma_cff.py.

Also the primary vertex for photon reconstruction is taken from heavy ion algorithm.

How to run the photon reconstruction

Do you need to run the photon reconstruction?

By default the sequences

are called in RecoHI/Configuration/python/Reconstruction_HI_cff.py during the standard heavy ion reconstruction process. If you want to include the egamma isolation variables, one can include the
in the scquence to produce heavy ion isolation maps.

Running the photon reconstruction only

In your configuration file, include the HiEgammaAlgos_cff by

# Egamma

and put hiEcalClusters and hiEgammaSequence into your cms.Sequence. This will reconstruct photons with heavy ion specific algorithm up to AOD photon level.

Running the full heavy ion reconstruction

In your configuration file, include the Reconstruction_HI_cff.py by


and put globalRecoPbPb into your cms.Sequence.

Heavy ion specific isolation variables:

The heavyion specific isolation variables are categorized as four types:

  • isoCCx: background subtracted isolation cone energy in ECAL with cone size x*0.1
  • isoCRx: background subtracted isolation cone energy in HCAL with cone size x*0.1
  • isoTxy: number of reconstructed tracks in a cone of x*0.1 with pt > y * 0.4 GeV/c
  • isoRxy: in a cone size of isoRxy, the number of reconstructed tracks with pt > y * 0.4 GeV/c is smaller than x.

They are accessible in the PAT level as user floats. Before the PAT level, one can access them directly as value maps to AOD Photons.

PAT Photon reconstruction

After running the full heavy ion reconstruction, one can include UserCode/CmsHi/HiPatAlgos/python/HiPatPhotonSequence_cfi.py and run the heavy ion PAT photon reconstruction by using the

in the configuration file.

This sequence will process the reconstructed photon and embed the heavy ion photon isolation information into PAT::Photon as user floats. At the same time, p+p default photon ID and isolation are also stored for comparisons / studies.

Accessing the heavy ion photons isolation variable

A test module is prepared in RecoHI/HiEgammaAlgos/test/patTest/PATHIPhotonTestModule.cc, which works on reconstructed heavyion PAT Photon.

One can access the information by the follow example code:

    edm::Handle > photons;
   std::auto_ptr > output(new std::vector());

to get the PAT photon collection and loop over the photon candidates by:

  for (edm::View::const_iterator photon = photons->begin(), end = photons->end(); photon != end; ++photon) {

Inside the loop, one can access the embedded isolation variables by

   Float_t isoCC1 = photon->userFloat("isoCC1");
   Float_t isoCR1 = photon->userFloat("isoCR1");
   Float_t isoT12 = photon->userFloat("isoT12");
   Float_t isoDR23 = photon->userFloat("isoDR23");

for your own analysis.


Responsible: YenJieLee

Edit | Attach | Watch | Print version | History: r8 < r7 < r6 < r5 < r4 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r6 - 2009-11-17 - YenJieLee
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    CMSPublic All webs login

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