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 GenericHistSFApplier< DAS::RecPhoton >
	GenericHistSFApplier (const std::filesystem::path &filePath, bool correction, bool uncertainties)

virtual	~GenericHistSFApplier () 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 GenericHistSFApplier< 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 GenericHistSFApplier< DAS::RecPhoton >
enum	Interpretation

Constructor & Destructor Documentation

◆ ConversionVetoApplier()

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

inline

         : GenericHistSFApplier(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

Description

Public Member Functions

Protected Member Functions

Private Attributes