IVGammaORTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2018 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef GAMMAORTOOLS_IVGAMMAORTOOL_H
#define GAMMAORTOOLS_IVGAMMAORTOOL_H
 
#include <vector>
#include <map>
 
#include <xAODTruth/TruthParticleContainer.h>
#include <xAODTruth/TruthParticle.h>
#include "AsgTools/IAsgTool.h"
 
/****************************************************************************
*   This tool can be used to identify the overlap of a V+jets and V+gamma sample
*   or in fact between most V^n+gamma^m and V^n+gamma^(m-1) samples
*   (e.g. V+gammagamma vs V+gamma or VV+gamma vs VV).
*   It can either decide whether an event is in the overlap region for fixed cuts on photons and lepton
*   or be used to extract relevant information from the event so that the overlap removal (OR) can be performed
*   at a later stage with a simple 'if' statement.
*
*   The OR functions (inOverlap(result), photonPtsOutsideDr(result), or photonPtsOutsideDrs(result)),
*   set the result via reference and return a status code.
*   Usually the tool will take leptons and photons from the current event but this behavior can be
*   overwritten by explicitly setting leptons and photons as additional arguments to the functions.
*
*   Further documented at
*   https://twiki.cern.ch/twiki/bin/viewauth/AtlasProtected/VGammaORTool
*
*********************************************************************************/
 
class IVGammaORTool: virtual public asg::IAsgTool {
public:
  ASG_TOOL_INTERFACE(IVGammaORTool)
 
  
  //    four vectors (can be used if the full truth record is not available or broken)
  virtual StatusCode inOverlap(bool& result,
                   const std::vector<TLorentzVector>* leptons=0,
                   const std::vector<TLorentzVector>* photons=0,
                   const std::vector<int>* lepton_origins=0,
                   const std::vector<int>* photon_origins=0) const =0;

  //    four vectors (can be used if the full truth record is not available or broken)
  virtual StatusCode photonPtsOutsideDr(std::vector<float>& result,
                    const std::vector<TLorentzVector>* leptons=0,
                    const std::vector<TLorentzVector>* photons=0,
                    const std::vector<int>* lepton_origins=0,
                    const std::vector<int>* photon_origins=0) const=0;

  //    four vectors (can be used if the full truth record is not available or broken)
  virtual StatusCode photonPtsOutsideDrs(std::map<float, std::vector<float> >& result,
                     const std::vector<TLorentzVector>* leptons=0,
                     const std::vector<TLorentzVector>* photons=0,
                     const std::vector<int>* lepton_origins=0,
                     const std::vector<int>* photon_origins=0) const=0;

  virtual bool frixioneIsolated(const xAOD::TruthParticle& photon,
                const xAOD::TruthParticleContainer& truthParticles,
                float dR0, float exponent, float epsilon) const=0;
};
 
 
 
#endif //> !GAMMAORTOOLS_IVGAMMAORTOOL_H