EcalIsolationCorrector Class Reference

#include <EcalIsolationCorrector.h>

Public Types
enum	RunRange { RunAB, RunC, RunD }

Public Member Functions
float	correctForHLTDefinition (reco::GsfElectron e, bool isData=false, float intL_AB=5.5, float intL_C=6.7, float intL_D=7.3)
float	correctForHLTDefinition (reco::GsfElectron e, int runNumber, bool isData=false)
float	correctForHLTDefinition (reco::GsfElectron e, std::string runName, bool isData=false)
float	correctForHLTDefinition (reco::Photon p, bool isData=false, float intL_AB=5.5, float intL_C=6.7, float intL_D=7.3)
float	correctForHLTDefinition (reco::Photon p, int runNumber, bool isData=false)
float	correctForHLTDefinition (reco::Photon p, std::string runName, bool isData=false)
float	correctForHLTDefinition (pat::Electron e, bool isData=false, float intL_AB=5.5, float intL_C=6.7, float intL_D=7.3)
float	correctForHLTDefinition (pat::Electron e, int runNumber, bool isData=false)
float	correctForHLTDefinition (pat::Electron e, std::string runName, bool isData=false)
float	correctForHLTDefinition (pat::Photon p, bool isData=false, float intL_AB=5.5, float intL_C=6.7, float intL_D=7.3)
float	correctForHLTDefinition (pat::Photon p, int runNumber, bool isData=false)
float	correctForHLTDefinition (pat::Photon p, std::string runName, bool isData=false)
float	correctForNoise (reco::GsfElectron e, bool isData=false, float intL_AB=5.5, float intL_C=6.7, float intL_D=7.3)
float	correctForNoise (reco::GsfElectron e, int runNumber, bool isData=false)
float	correctForNoise (reco::GsfElectron e, std::string runName, bool isData=false)
float	correctForNoise (reco::Photon p, bool isData=false, float intL_AB=5.5, float intL_C=6.7, float intL_D=7.3)
float	correctForNoise (reco::Photon p, int runNumber, bool isData=false)
float	correctForNoise (reco::Photon p, std::string runName, bool isData=false)
float	correctForNoise (pat::Electron e, bool isData=false, float intL_AB=5.5, float intL_C=6.7, float intL_D=7.3)
float	correctForNoise (pat::Electron e, int runNumber, bool isData=false)
float	correctForNoise (pat::Electron e, std::string runName, bool isData=false)
float	correctForNoise (pat::Photon p, bool isData=false, float intL_AB=5.5, float intL_C=6.7, float intL_D=7.3)
float	correctForNoise (pat::Photon p, int runNumber, bool isData=false)
float	correctForNoise (pat::Photon p, std::string runName, bool isData=false)
	EcalIsolationCorrector (bool forElectrons)
	~EcalIsolationCorrector ()

Protected Member Functions
RunRange	checkRunRange (int runNumber)
float	correctForHLTDefinition (float iso, bool isBarrel, RunRange runRange)
float	correctForNoise (float iso, bool isBarrel, RunRange runRange, bool isData)

Private Attributes
bool	isElectron_

Detailed Description

Definition at line 17 of file EcalIsolationCorrector.h.

Member Enumeration Documentation

enum EcalIsolationCorrector::RunRange

Enumerator
RunAB
RunC
RunD

Definition at line 20 of file EcalIsolationCorrector.h.

{RunAB, RunC, RunD};

Constructor & Destructor Documentation

EcalIsolationCorrector::EcalIsolationCorrector

(

bool

forElectrons

)

Definition at line 7 of file EcalIsolationCorrector.cc.

                                                               :isElectron_(forElectrons)
{}

EcalIsolationCorrector::~EcalIsolationCorrector ( )

inline

Definition at line 23 of file EcalIsolationCorrector.h.

{};

Member Function Documentation

EcalIsolationCorrector::RunRange EcalIsolationCorrector::checkRunRange ( int  runNumber )

protected

Definition at line 10 of file EcalIsolationCorrector.cc.

References RunAB, RunC, and RunD.

Referenced by correctForHLTDefinition(), and correctForNoise().

                                                                                  {
  
  EcalIsolationCorrector::RunRange runRange = RunAB;
  
  if (runNumber <= 203755 && runNumber > 197770)
    runRange = RunC;
  else if (runNumber > 203755)
    runRange = RunD;
  
  return runRange;
}

float EcalIsolationCorrector::correctForHLTDefinition	(	reco::GsfElectron	e,
		bool	isData = false,
		float	intL_AB = 5.5,
		float	intL_C = 6.7,
		float	intL_D = 7.3
	)

Definition at line 380 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), isElectron_, RunAB, RunC, and RunD.

Referenced by correctForHLTDefinition().

                                                                                                                               {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
  
  float iso = e.dr03EcalRecHitSumEt();
  if (!isData)
    iso = correctForNoise(e, isData, intL_AB, intL_C, intL_D);
  
  float combination = (intL_AB * correctForHLTDefinition(iso, e.isEB(), RunAB) +
               intL_C  * correctForHLTDefinition(iso, e.isEB(), RunC ) +
               intL_D  * correctForHLTDefinition(iso, e.isEB(), RunD ))/(intL_AB + intL_C + intL_D);
  
  return combination;
}

float EcalIsolationCorrector::correctForHLTDefinition	(	reco::GsfElectron	e,
		int	runNumber,
		bool	isData = false
	)

Definition at line 342 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForHLTDefinition(), correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), and isElectron_.

                                                                                                   {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
  
  float iso = e.dr03EcalRecHitSumEt();
  EcalIsolationCorrector::RunRange runRange = checkRunRange(runNumber);
  
  if (!isData)
    iso = correctForNoise(iso, e.isEB(), runRange, false);
  
  return correctForHLTDefinition(iso, e.isEB(), runRange);
}

float EcalIsolationCorrector::correctForHLTDefinition	(	reco::GsfElectron	e,
		std::string	runName,
		bool	isData = false
	)

Definition at line 356 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForHLTDefinition(), correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                                       {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
  
  float iso = e.dr03EcalRecHitSumEt();
  EcalIsolationCorrector::RunRange runRange = RunAB;
  if (runName == "RunAB") 
    runRange = RunAB;
  else if (runName == "RunC")
    runRange = RunC;
  else if (runName == "RunD")
    runRange = RunD;
  else {
    std::cerr << "Error: Unknown run range " << runName << std::endl;
    abort();
  }
  
  if (!isData)
    iso = correctForNoise(iso, e.isEB(), runRange, false);
  
  return correctForHLTDefinition(iso, e.isEB(), runRange);
}

float EcalIsolationCorrector::correctForHLTDefinition	(	reco::Photon	p,
		bool	isData = false,
		float	intL_AB = 5.5,
		float	intL_C = 6.7,
		float	intL_D = 7.3
	)

Definition at line 491 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForHLTDefinition(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                                                          {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  
  if (!isData)
    iso = correctForNoise(p, isData, intL_AB, intL_C, intL_D);
  
  float combination = (intL_AB * correctForHLTDefinition(iso, p.isEB(), RunAB) +
               intL_C  * correctForHLTDefinition(iso, p.isEB(), RunC ) +
               intL_D  * correctForHLTDefinition(iso, p.isEB(), RunD ))/(intL_AB + intL_C + intL_D);
  
  return combination;
}

float EcalIsolationCorrector::correctForHLTDefinition	(	reco::Photon	p,
		int	runNumber,
		bool	isData = false
	)

Definition at line 452 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForHLTDefinition(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), and isElectron_.

                                                                                              {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  EcalIsolationCorrector::RunRange runRange = checkRunRange(runNumber);
  
  if (!isData)
    iso = correctForNoise(iso, p.isEB(), runRange, false);
  
  return correctForHLTDefinition(iso, p.isEB(), runRange);
}

float EcalIsolationCorrector::correctForHLTDefinition	(	reco::Photon	p,
		std::string	runName,
		bool	isData = false
	)

Definition at line 466 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForHLTDefinition(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                                  {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  
  EcalIsolationCorrector::RunRange runRange = RunAB;
  if (runName == "RunAB") 
    runRange = RunAB;
  else if (runName == "RunC")
    runRange = RunC;
  else if (runName == "RunD")
    runRange = RunD;
  else {
    std::cerr << "Error: Unknown run range " << runName << std::endl;
    abort();
  }
  
  if (!isData)
    iso = correctForNoise(iso, p.isEB(), runRange, false);
  
  return correctForHLTDefinition(iso, p.isEB(), runRange);
}

float EcalIsolationCorrector::correctForHLTDefinition	(	pat::Electron	e,
		bool	isData = false,
		float	intL_AB = 5.5,
		float	intL_C = 6.7,
		float	intL_D = 7.3
	)

Definition at line 435 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForHLTDefinition(), correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                                                           {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
  
  float iso = e.dr03EcalRecHitSumEt();
  if (!isData)
    iso = correctForNoise(e, isData, intL_AB, intL_C, intL_D);
  
  float combination = (intL_AB * correctForHLTDefinition(iso, e.isEB(), RunAB) +
               intL_C  * correctForHLTDefinition(iso, e.isEB(), RunC ) +
               intL_D  * correctForHLTDefinition(iso, e.isEB(), RunD ))/(intL_AB + intL_C + intL_D);
  
  return combination;
}

float EcalIsolationCorrector::correctForHLTDefinition	(	pat::Electron	e,
		int	runNumber,
		bool	isData = false
	)

Definition at line 397 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForHLTDefinition(), correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), and isElectron_.

                                                                                               {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
  
  float iso = e.dr03EcalRecHitSumEt();
  EcalIsolationCorrector::RunRange runRange = checkRunRange(runNumber);
  
  if (!isData)
    iso = correctForNoise(iso, e.isEB(), runRange, false);
  
  return correctForHLTDefinition(iso, e.isEB(), runRange);
}

float EcalIsolationCorrector::correctForHLTDefinition	(	pat::Electron	e,
		std::string	runName,
		bool	isData = false
	)

Definition at line 411 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForHLTDefinition(), correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                                   {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
  
  float iso = e.dr03EcalRecHitSumEt();
  EcalIsolationCorrector::RunRange runRange = RunAB;
  if (runName == "RunAB") 
    runRange = RunAB;
  else if (runName == "RunC")
    runRange = RunC;
  else if (runName == "RunD")
    runRange = RunD;
  else {
    std::cerr << "Error: Unknown run range " << runName << std::endl;
    abort();
  }
  
  if (!isData)
    iso = correctForNoise(iso, e.isEB(), runRange, false);
  
  return correctForHLTDefinition(iso, e.isEB(), runRange);
}

float EcalIsolationCorrector::correctForHLTDefinition	(	pat::Photon	p,
		bool	isData = false,
		float	intL_AB = 5.5,
		float	intL_C = 6.7,
		float	intL_D = 7.3
	)

Definition at line 548 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForHLTDefinition(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                                                         {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  
  if (!isData)
    iso = correctForNoise(p, isData, intL_AB, intL_C, intL_D);
  
  float combination = (intL_AB * correctForHLTDefinition(iso, p.isEB(), RunAB) +
               intL_C  * correctForHLTDefinition(iso, p.isEB(), RunC ) +
               intL_D  * correctForHLTDefinition(iso, p.isEB(), RunD ))/(intL_AB + intL_C + intL_D);
  
  return combination;
}

float EcalIsolationCorrector::correctForHLTDefinition	(	pat::Photon	p,
		int	runNumber,
		bool	isData = false
	)

Definition at line 509 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForHLTDefinition(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), and isElectron_.

                                                                                             {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  EcalIsolationCorrector::RunRange runRange = checkRunRange(runNumber);
  
  if (!isData)
    iso = correctForNoise(iso, p.isEB(), runRange, false);
  
  return correctForHLTDefinition(iso, p.isEB(), runRange);
}

float EcalIsolationCorrector::correctForHLTDefinition	(	pat::Photon	p,
		std::string	runName,
		bool	isData = false
	)

Definition at line 523 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForHLTDefinition(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                                 {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  
  EcalIsolationCorrector::RunRange runRange = RunAB;
  if (runName == "RunAB") 
    runRange = RunAB;
  else if (runName == "RunC")
    runRange = RunC;
  else if (runName == "RunD")
    runRange = RunD;
  else {
    std::cerr << "Error: Unknown run range " << runName << std::endl;
    abort();
  }
  
  if (!isData)
    iso = correctForNoise(iso, p.isEB(), runRange, false);
  
  return correctForHLTDefinition(iso, p.isEB(), runRange);
}

float EcalIsolationCorrector::correctForHLTDefinition ( float  iso, bool  isBarrel, EcalIsolationCorrector::RunRange  runRange  )

protected

Definition at line 297 of file EcalIsolationCorrector.cc.

References isElectron_, mps_fire::result, RunAB, RunC, and RunD.

                                                                                                                       {
  
  float result = iso;
  
  if (!isElectron_) {
    if (runRange == RunAB) {
      if (!isBarrel)
    result = iso*0.8499-0.6510;
      else
    result = iso*0.8504-0.5658;
    } else if (runRange == RunC) {
      if (!isBarrel)
    result = iso*0.9346-0.9987;
      else
    result = iso*0.8529-0.6816;
    } else if (runRange == RunD) {
      if (!isBarrel) 
    result = iso*0.8318-0.9999;
      else
    result = iso*0.8853-0.8783;
    }
  } else {
    if (runRange == RunAB) {
      if (!isBarrel)
    result = iso*0.9849-0.6871;
      else
    result = iso*0.8542-0.3558;
    } else if (runRange == RunC) {
      if (!isBarrel)
    result = iso*0.9996-0.8485;
      else
    result = iso*0.9994-0.5085;
    } else if (runRange == RunD) {
      if (!isBarrel) 
    result = iso*0.9467-0.9998;
      else
    result = iso*0.8574-0.4862;
    }
  }
  
  return result;
}

float EcalIsolationCorrector::correctForNoise	(	reco::GsfElectron	e,
		bool	isData = false,
		float	intL_AB = 5.5,
		float	intL_C = 6.7,
		float	intL_D = 7.3
	)

Definition at line 108 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), isElectron_, RunAB, RunC, and RunD.

Referenced by correctForHLTDefinition(), and correctForNoise().

                                                                                                                       {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
  
  float iso = e.dr03EcalRecHitSumEt();
  float combination = (intL_AB * correctForNoise(iso, e.isEB(), RunAB, isData) +
               intL_C  * correctForNoise(iso, e.isEB(), RunC,  isData) +
               intL_D  * correctForNoise(iso, e.isEB(), RunD,  isData))/(intL_AB + intL_C + intL_D);
  
  return combination;
}

float EcalIsolationCorrector::correctForNoise	(	reco::GsfElectron	e,
		int	runNumber,
		bool	isData = false
	)

Definition at line 75 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), and isElectron_.

                                                                                           {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
  
  float iso = e.dr03EcalRecHitSumEt();
  EcalIsolationCorrector::RunRange runRange = checkRunRange(runNumber);
  
  return correctForNoise(iso, e.isEB(), runRange, isData);
}

float EcalIsolationCorrector::correctForNoise	(	reco::GsfElectron	e,
		std::string	runName,
		bool	isData = false
	)

Definition at line 86 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                               {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;  
  
  float iso = e.dr03EcalRecHitSumEt();
  EcalIsolationCorrector::RunRange runRange = RunAB;
  
  if (runName == "RunAB") 
    runRange = RunAB;
  else if (runName == "RunC")
    runRange = RunC;
  else if (runName == "RunD")
    runRange = RunD;
  else {
    std::cerr << "Error: Unknown run range " << runName << std::endl;
    abort();
  }
  
  return correctForNoise(iso, e.isEB(), runRange, isData);
}

float EcalIsolationCorrector::correctForNoise	(	reco::Photon	p,
		bool	isData = false,
		float	intL_AB = 5.5,
		float	intL_C = 6.7,
		float	intL_D = 7.3
	)

Definition at line 202 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                                                  {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  float combination = (intL_AB * correctForNoise(iso, p.isEB(), RunAB, isData) +
               intL_C  * correctForNoise(iso, p.isEB(), RunC,  isData) +
               intL_D  * correctForNoise(iso, p.isEB(), RunD,  isData))/(intL_AB + intL_C + intL_D);
  
  return combination;
}

float EcalIsolationCorrector::correctForNoise	(	reco::Photon	p,
		int	runNumber,
		bool	isData = false
	)

Definition at line 169 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), and isElectron_.

                                                                                      {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  EcalIsolationCorrector::RunRange runRange = checkRunRange(runNumber);
  
  return correctForNoise(iso, p.isEB(), runRange, isData);
}

float EcalIsolationCorrector::correctForNoise	(	reco::Photon	p,
		std::string	runName,
		bool	isData = false
	)

Definition at line 180 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                          {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  EcalIsolationCorrector::RunRange runRange = RunAB;
  
  if (runName == "RunAB") 
    runRange = RunAB;
  else if (runName == "RunC")
    runRange = RunC;
  else if (runName == "RunD")
    runRange = RunD;
  else {
    std::cerr << "Error: Unknown run range " << runName << std::endl;
    abort();
  }
  
  return correctForNoise(iso, p.isEB(), runRange, isData);
}

float EcalIsolationCorrector::correctForNoise	(	pat::Electron	e,
		bool	isData = false,
		float	intL_AB = 5.5,
		float	intL_C = 6.7,
		float	intL_D = 7.3
	)

Definition at line 155 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                                                   {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
  
  float iso = e.dr03EcalRecHitSumEt();
  float combination = (intL_AB * correctForNoise(iso, e.isEB(), RunAB, isData) +
               intL_C  * correctForNoise(iso, e.isEB(), RunC,  isData) +
               intL_D  * correctForNoise(iso, e.isEB(), RunD,  isData))/(intL_AB + intL_C + intL_D);
  
  return combination;
}

float EcalIsolationCorrector::correctForNoise	(	pat::Electron	e,
		int	runNumber,
		bool	isData = false
	)

Definition at line 122 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), and isElectron_.

                                                                                       {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;
  
  float iso = e.dr03EcalRecHitSumEt();
  EcalIsolationCorrector::RunRange runRange = checkRunRange(runNumber);
  
  return correctForNoise(iso, e.isEB(), runRange, isData);
}

float EcalIsolationCorrector::correctForNoise	(	pat::Electron	e,
		std::string	runName,
		bool	isData = false
	)

Definition at line 133 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), reco::GsfElectron::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                           {
  
  if (!isElectron_)
    std::cerr << "Warning: this corrector is setup for photons and you are passing an electron !" << std::endl;  
  
  float iso = e.dr03EcalRecHitSumEt();
  EcalIsolationCorrector::RunRange runRange = RunAB;
  
  if (runName == "RunAB") 
    runRange = RunAB;
  else if (runName == "RunC")
    runRange = RunC;
  else if (runName == "RunD")
    runRange = RunD;
  else {
    std::cerr << "Error: Unknown run range " << runName << std::endl;
    abort();
  }
  
  return correctForNoise(iso, e.isEB(), runRange, isData);
}

float EcalIsolationCorrector::correctForNoise	(	pat::Photon	p,
		bool	isData = false,
		float	intL_AB = 5.5,
		float	intL_C = 6.7,
		float	intL_D = 7.3
	)

Definition at line 249 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                                                 {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  float combination = (intL_AB * correctForNoise(iso, p.isEB(), RunAB, isData) +
               intL_C  * correctForNoise(iso, p.isEB(), RunC,  isData) +
               intL_D  * correctForNoise(iso, p.isEB(), RunD,  isData))/(intL_AB + intL_C + intL_D);
  
  return combination;
}

float EcalIsolationCorrector::correctForNoise	(	pat::Photon	p,
		int	runNumber,
		bool	isData = false
	)

Definition at line 216 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), and isElectron_.

                                                                                     {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  EcalIsolationCorrector::RunRange runRange = checkRunRange(runNumber);
  
  return correctForNoise(iso, p.isEB(), runRange, isData);
}

float EcalIsolationCorrector::correctForNoise	(	pat::Photon	p,
		std::string	runName,
		bool	isData = false
	)

Definition at line 227 of file EcalIsolationCorrector.cc.

References ecal_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), reco::Photon::isEB(), isElectron_, RunAB, RunC, and RunD.

                                                                                         {
  
  if (isElectron_)
    std::cerr << "Warning: this corrector is setup for electrons and you are passing a photon !" << std::endl;
  
  float iso = p.ecalRecHitSumEtConeDR03();
  EcalIsolationCorrector::RunRange runRange = RunAB;
  
  if (runName == "RunAB") 
    runRange = RunAB;
  else if (runName == "RunC")
    runRange = RunC;
  else if (runName == "RunD")
    runRange = RunD;
  else {
    std::cerr << "Error: Unknown run range " << runName << std::endl;
    abort();
  }
  
  return correctForNoise(iso, p.isEB(), runRange, isData);
}

float EcalIsolationCorrector::correctForNoise ( float  iso, bool  isBarrel, EcalIsolationCorrector::RunRange  runRange, bool  isData  )

protected

Definition at line 22 of file EcalIsolationCorrector.cc.

References gather_cfg::cout, isElectron_, mps_fire::result, RunAB, RunC, and RunD.

                                                                                                                            {
  
  float result = iso;
  
  if (!isElectron_) {
    if (!isData) {
      if (runRange == RunAB) {
    if (!isBarrel)
      result = (iso-0.2827)/1.0949;
    else
      result = (iso-0.0931)/1.0738;
      } else if (runRange == RunC) {
    if (!isBarrel)
      result = (iso-0.5690)/0.9217;
    else
      result = (iso-0.1824)/0.9279;
      } else if (runRange == RunD) {
    if (!isBarrel) 
      result = (iso-0.9997)/0.8781;
    else
      result = (iso-0.0944)/0.8140;
      }
    } else {
      std::cout << "Warning: you should correct MC to data" << std::endl;
    }
  } else {
    if (!isData) {
      if (runRange == RunAB) {
    if (!isBarrel)
      result = (iso+0.1174)/1.0012;
    else
      result = (iso+0.2736)/0.9948;
      } else if (runRange == RunC) {
    if (!isBarrel)
      result = (iso+0.2271)/0.9684;
    else
      result = (iso+0.5962)/0.9568;
      } else if (runRange == RunD) {
    if (!isBarrel) 
      result = (iso+0.2907)/1.0005;
    else
      result = (iso+0.9098)/0.9395;
      }
    } else {
      std::cout << "Warning: you should correct MC to data" << std::endl;
    } 
  }
  
  return result;
}

Member Data Documentation

bool EcalIsolationCorrector::isElectron_

private

Definition at line 93 of file EcalIsolationCorrector.h.

Referenced by correctForHLTDefinition(), and correctForNoise().