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 19 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 22 of file EcalIsolationCorrector.h.

References checkRunRange(), correctForHLTDefinition(), correctForNoise(), MillePedeFileConverter_cfg::e, PixelPluginsPhase0_cfi::isBarrel, HLTBitAnalyser_cfi::isData, AlCaHLTBitMon_ParallelJobs::p, convertSQLiteXML::runNumber, and AlCaHLTBitMon_QueryRunRegistry::string.

{};

Member Function Documentation

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

protected

Definition at line 9 of file EcalIsolationCorrector.cc.

References RunAB, RunC, and RunD.

Referenced by correctForHLTDefinition(), correctForNoise(), and ~EcalIsolationCorrector().

                                                                                  {
  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 371 of file EcalIsolationCorrector.cc.

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

Referenced by correctForHLTDefinition(), and ~EcalIsolationCorrector().

                                                                               {
  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 335 of file EcalIsolationCorrector.cc.

References beam_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 348 of file EcalIsolationCorrector.cc.

References beam_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 480 of file EcalIsolationCorrector.cc.

References beam_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 443 of file EcalIsolationCorrector.cc.

References beam_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 456 of file EcalIsolationCorrector.cc.

References beam_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 425 of file EcalIsolationCorrector.cc.

References beam_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 389 of file EcalIsolationCorrector.cc.

References beam_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 402 of file EcalIsolationCorrector.cc.

References beam_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 536 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForHLTDefinition(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), PixelPluginsPhase0_cfi::isBarrel, HLTBitAnalyser_cfi::isData, reco::Photon::isEB(), isElectron_, RunAB, RunC, RunD, convertSQLiteXML::runNumber, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                         {
  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 499 of file EcalIsolationCorrector.cc.

References beam_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 512 of file EcalIsolationCorrector.cc.

References beam_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 289 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 106 of file EcalIsolationCorrector.cc.

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

Referenced by correctForHLTDefinition(), correctForNoise(), and ~EcalIsolationCorrector().

                                                                               {
  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 beam_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), HLTBitAnalyser_cfi::isData, 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 85 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), HLTBitAnalyser_cfi::isData, 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 197 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), HLTBitAnalyser_cfi::isData, 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 166 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), HLTBitAnalyser_cfi::isData, 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 176 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), HLTBitAnalyser_cfi::isData, 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 152 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), HLTBitAnalyser_cfi::isData, 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 121 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), HLTBitAnalyser_cfi::isData, 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 131 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::GsfElectron::dr03EcalRecHitSumEt(), HLTBitAnalyser_cfi::isData, 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 242 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), PixelPluginsPhase0_cfi::isBarrel, HLTBitAnalyser_cfi::isData, reco::Photon::isEB(), isElectron_, RunAB, RunC, RunD, convertSQLiteXML::runNumber, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                                 {
  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 211 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, checkRunRange(), correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), HLTBitAnalyser_cfi::isData, 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 221 of file EcalIsolationCorrector.cc.

References beam_dqm_sourceclient-live_cfg::cerr, correctForNoise(), reco::Photon::ecalRecHitSumEtConeDR03(), HLTBitAnalyser_cfi::isData, 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 20 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 105 of file EcalIsolationCorrector.h.

Referenced by correctForHLTDefinition(), and correctForNoise().