ConversionVetoApplier

Description

Applies the Conversion Safe Electron Veto and corrections.

Modifies the weights of the reconstructed photons to apply the efficiency correction. Scale factors are taken from a ROOT file in the format provided by Egamma for UL scale factors.

Inheritance diagram for ConversionVetoApplier:

Collaboration diagram for ConversionVetoApplier:

Public Member Functions
	ConversionVetoApplier (const fs::path &filePath, const string &ID, bool correction, bool uncertainties, const string &type)
Public Member Functions inherited from GenericSFApplier< DAS::RecPhoton >
	GenericSFApplier (const std::filesystem::path &filePath, bool correction, bool uncertainties)
virtual	~GenericSFApplier () noexcept=default
void	operator() (DAS::RecPhoton &object, const Context &... ctx) const
void	operator() (std::vector< DAS::RecPhoton > &objects, const Context &... ctx) const
std::vector< std::string >	weightNames () const

Protected Member Functions
bool	passes (const RecPhoton &photon) const override
int	binIndex (const RecPhoton &photon, const std::unique_ptr< TH1 > &hist) const override
Protected Member Functions inherited from GenericSFApplier< DAS::RecPhoton >
void	loadNominal (const std::string &histPath)
void	loadBinWiseUnc (const std::string &name, const std::string &histPath, Interpretation intp=UseBinContent)
void	loadGlobalUnc (const std::string &name, const std::string &histPath, Interpretation intp=UseBinContent)
virtual bool	passes (const DAS::RecPhoton &obj, const Context &... ctx) const=0
virtual int	binIndex (const DAS::RecPhoton &obj, const Context &... ctx, const std::unique_ptr< TH1 > &hist) const=0

Private Attributes
RecPhoton::Identification	selector

Additional Inherited Members
Protected Types inherited from GenericSFApplier< DAS::RecPhoton >
enum	Interpretation

Constructor & Destructor Documentation

ConversionVetoApplier()

ConversionVetoApplier ( const fs::path &  filePath, const string &  ID, bool  correction, bool  uncertainties, const string &  type  )

inline

        : GenericSFApplier(filePath, correction, uncertainties)
    {
        /*
         * Meta info: Loose | Medium | Tight | MVA80 | MVA90
         * Egamma conveners believe that the 'MVA' working point in the files
         * corresponds to MVA90. Reference:
         *   https://mattermost.web.cern.ch/cmseg/pl/yewftned9b8bxg4m3y7gdck4ze
         */
        string IDName = ID + "ID";
        if (ID == "MVA80")
            BOOST_THROW_EXCEPTION(std::invalid_argument(
                "There are no CSEV scale factors for the MVA80 ID."));
        else if (ID == "MVA90")
            IDName = "MVAID";
 
        /*
         * root [1] .ls
         * TFile**         /nfs/dust/cms/user/mourelou/darwin/Darwin/tables/Egamma/Photon_CSEV_UL18.root
         *  TFile*         /nfs/dust/cms/user/mourelou/darwin/Darwin/tables/Egamma/Photon_CSEV_UL18.root
         *   KEY: TDirectoryFile   LooseID;1       LooseID
         *   KEY: TDirectoryFile   MediumID;1      MediumID
         *   KEY: TDirectoryFile   TightID;1       TightID
         *   KEY: TDirectoryFile   MVAID;1 MVAID
         * root [2] LooseID->ls()
         * TDirectoryFile*         LooseID LooseID
         *  KEY: TH1F      SF_CSEV_LooseID;1
         *  KEY: TH1F      Staunc_CSEV_LooseID;1
         *  KEY: TH1F      PUunc_CSEV_LooseID;1
         *  KEY: TH1F      Modelunc_CSEV_LooseID;1
         */
 
        loadNominal(IDName + "/SF_" + type + "_" + IDName);
        loadBinWiseUnc(type + "_stat", IDName + "/Staunc_" + type + "_" + IDName);
        loadGlobalUnc(type + "_PU",    IDName + "/PUunc_" + type + "_" + IDName);
        loadGlobalUnc(type + "_model", IDName + "/Modelunc_" + type + "_" + IDName);
 
        if (type == "CSEV")
            selector = RecPhoton::ConversionSafeElectronVeto;
        else if (type == "HasPix")
            selector = RecPhoton::PixelSeedVeto;
        else
            BOOST_THROW_EXCEPTION(invalid_argument("Type `"
                                  + type + "`is unknown"));
        cout << "Initialisation of " << __func__ << " done" << endl;
    }

Member Function Documentation

binIndex()

int binIndex ( const RecPhoton &  photon, const std::unique_ptr< TH1 > &  hist  ) const

inlineoverrideprotected

    {
        /*
         * The axis is categorical:
         *   0 = EB inclusive
         *   1 = EB high R9
         *   2 = EB low R9
         *   3 = EE inclusive
         *   4 = EE high R9
         *   5 = EE low R9
         * The threshold for high R9 is 0.96 (using the "shifted" value).
         * Since we don't have R9 in the ntuples we use the inclusive.
         */
        return hist->FindBin(abs(photon.scEta) < 1.5 ? 0 : 3);
    }

passes()

bool passes ( const RecPhoton &  photon ) const

inlineoverrideprotected

    {
        return selector & photon.selectors;
    }

Member Data Documentation

selector

RecPhoton::Identification selector

private

expecting PSV or CSEV

---

The documentation for this class was generated from the following file:

• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Photons/bin/applyPhotonConversionVeto.cc