PhotonPlots.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#include "PhotonPlots.h"
#include "xAODEgamma/EgammaxAODHelpers.h"
#include "xAODEgamma/PhotonxAODHelpers.h"
#include "MCTruthClassifier/MCTruthClassifierDefs.h"
using CLHEP::GeV;
 
namespace Egamma{
 
 
PhotonPlots::PhotonPlots(PlotBase* pParent, const std::string& sDir, const std::string& sParticleType):PlotBase(pParent, sDir),
                                             // inclusive plots, no isolation , ID or truth-matching requirements 
                                             m_oKinAllPlots(this, "All/KinPlots/", "Reco " + sParticleType +" Photon"), 
                                             m_oShowerShapesAllPlots(this, "All/ShowerShapesPlots/", sParticleType  +" Photon"),
                                             m_oShowerShapesAll10GeVPlots(this, "All/ShowerShapes10GeVPlots/", sParticleType  +" Photon"),
                                             m_oPhotonConversionAllPlots(this, "All/PhotonConversionPlots/", sParticleType  +" Photon"),
                                             // truth-matchint "Iso" requirement
                                             m_oKinIsoPlots(this, "Iso/KinPlots/", "Reco " + sParticleType +" Photon"),
                                             m_oShowerShapesIsoPlots(this, "Iso/ShowerShapesPlots/", sParticleType  +" Photon"),
                                             m_oShowerShapesIso10GeVPlots(this, "Iso/ShowerShapes10GeVPlots/", sParticleType  +" Photon"),      
                                             m_oPhotonConversionIsoPlots(this, "Iso/PhotonConversionPlots/", sParticleType  +" Photon"),
                                             m_oKinIsoUncPlots(this, "Iso/KinPlotsUnc/", "Reco " + sParticleType +" Photon"),
                                             m_oKinIsoConvPlots(this, "Iso/KinPlotsConv/", "Reco " + sParticleType +" Photon"),
                                             //
                                             m_oKinIsoLoosePlots(this, "IsoLoose/KinPlots/", "Reco " + sParticleType +" Photon"),
                                             m_oKinIsoLooseUncPlots(this, "IsoLoose/KinPlotsUnc/", "Reco " + sParticleType +" Photon"),
                                             m_oKinIsoLooseConvPlots(this, "IsoLoose/KinPlotsConv/", "Reco " + sParticleType +" Photon"),
                                             //
                                             m_oKinIsoTightPlots(this, "IsoTight/KinPlots/", "Reco " + sParticleType +" Photon"),
                                             m_oKinIsoTightUncPlots(this, "IsoTight/KinPlotsUnc/", "Reco " + sParticleType +" Photon"),
                                             m_oKinIsoTightConvPlots(this, "IsoTight/KinPlotsConv/", "Reco " + sParticleType +" Photon"),
                                             m_sParticleType(sParticleType),
                                             m_nParticles(nullptr),
                                             m_nParticles_weighted(nullptr)
{}
 
void PhotonPlots::initializePlots(){
  m_nParticles = Book1D("n", "Number of "+ m_sParticleType + "s;#" + m_sParticleType + "s;Events", 15, 0., 15.);
  m_nParticles_weighted = Book1D("n_weighted", "Number of "+ m_sParticleType + "s;#" + m_sParticleType + "s;Events", 15, 0., 15.);
}
 
  void PhotonPlots::fill(const xAOD::Photon& photon, const xAOD::EventInfo& eventInfo, bool isPrompt) {
  m_oKinAllPlots.fill(photon,eventInfo);
  m_oShowerShapesAllPlots.fill(photon,eventInfo);
  if (photon.pt() > 10 * Gaudi::Units::GeV) {
    m_oShowerShapesAll10GeVPlots.fill(photon,eventInfo);
  }
  m_oPhotonConversionAllPlots.fill(photon,eventInfo);
  
  if(!isPrompt) return;
 
  m_oKinIsoPlots.fill(photon,eventInfo);
  m_oShowerShapesIsoPlots.fill(photon,eventInfo);
  if (photon.pt() > 10 * Gaudi::Units::GeV) {
    m_oShowerShapesIso10GeVPlots.fill(photon,eventInfo);
  }
  
  m_oPhotonConversionIsoPlots.fill(photon,eventInfo);
  if (!xAOD::EgammaHelpers::isConvertedPhoton(&photon)) m_oKinIsoUncPlots.fill(photon,eventInfo);
  else m_oKinIsoConvPlots.fill(photon,eventInfo);
  
  bool val_loose=false;    
 
  photon.passSelection(val_loose, "Loose");
  if(val_loose) {
    m_oKinIsoLoosePlots.fill(photon,eventInfo);
    if (!xAOD::EgammaHelpers::isConvertedPhoton(&photon)) m_oKinIsoLooseUncPlots.fill(photon,eventInfo);
    else m_oKinIsoLooseConvPlots.fill(photon,eventInfo);
  }
  
  bool val_tight=false;    
 
  photon.passSelection(val_tight, "Tight");
  if(val_tight) {
    m_oKinIsoTightPlots.fill(photon,eventInfo);
    if (!xAOD::EgammaHelpers::isConvertedPhoton(&photon)) m_oKinIsoTightUncPlots.fill(photon,eventInfo);
    else m_oKinIsoTightConvPlots.fill(photon,eventInfo);
  }  
}
 
}