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

RunRange

enum EcalIsolationCorrector::RunRange

Enumerator
RunAB
RunC
RunD

Definition at line 19 of file EcalIsolationCorrector.h.

{ RunAB, RunC, RunD };

Constructor & Destructor Documentation

EcalIsolationCorrector()

EcalIsolationCorrector::EcalIsolationCorrector

(

bool

forElectrons

)

Definition at line 7 of file EcalIsolationCorrector.cc.

: isElectron_(forElectrons) {}

~EcalIsolationCorrector()

EcalIsolationCorrector::~EcalIsolationCorrector ( )

inline

Definition at line 22 of file EcalIsolationCorrector.h.

{};

Member Function Documentation

checkRunRange()

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

protected

Definition at line 9 of file EcalIsolationCorrector.cc.

References RunAB, RunC, RunD, convertSQLiteXML::runNumber, and crab_NZS2018C_1::runRange.

Referenced by correctForHLTDefinition(), and correctForNoise().

                                                                                  {
  EcalIsolationCorrector::RunRange runRange = RunAB;

  if (runNumber <= 203755 && runNumber > 197770)
    runRange = RunC;
  else if (runNumber > 203755)
    runRange = RunD;

  return runRange;
}

correctForHLTDefinition() [1/13]

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 EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, 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;
}

correctForHLTDefinition() [2/13]

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

Definition at line 335 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, checkRunRange(), correctForHLTDefinition(), correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, isElectron_, convertSQLiteXML::runNumber, and crab_NZS2018C_1::runRange.

                                                                                                   {
  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);
}

correctForHLTDefinition() [3/13]

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

Definition at line 348 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForHLTDefinition(), correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, isElectron_, RunAB, RunC, RunD, and crab_NZS2018C_1::runRange.

                                                                                                       {
  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);
}

correctForHLTDefinition() [4/13]

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 EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForHLTDefinition(), correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, 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;
}

correctForHLTDefinition() [5/13]

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

Definition at line 443 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, checkRunRange(), correctForHLTDefinition(), correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, convertSQLiteXML::runNumber, and crab_NZS2018C_1::runRange.

                                                                                              {
  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);
}

correctForHLTDefinition() [6/13]

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

Definition at line 456 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForHLTDefinition(), correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, RunAB, RunC, RunD, and crab_NZS2018C_1::runRange.

                                                                                                  {
  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);
}

correctForHLTDefinition() [7/13]

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 EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForHLTDefinition(), correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, 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;
}

correctForHLTDefinition() [8/13]

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

Definition at line 389 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, checkRunRange(), correctForHLTDefinition(), correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, isElectron_, convertSQLiteXML::runNumber, and crab_NZS2018C_1::runRange.

                                                                                               {
  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);
}

correctForHLTDefinition() [9/13]

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

Definition at line 402 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForHLTDefinition(), correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, isElectron_, RunAB, RunC, RunD, and crab_NZS2018C_1::runRange.

                                                                                                   {
  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);
}

correctForHLTDefinition() [10/13]

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 EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForHLTDefinition(), correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, 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;
}

correctForHLTDefinition() [11/13]

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

Definition at line 499 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, checkRunRange(), correctForHLTDefinition(), correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, convertSQLiteXML::runNumber, and crab_NZS2018C_1::runRange.

                                                                                             {
  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);
}

correctForHLTDefinition() [12/13]

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

Definition at line 512 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForHLTDefinition(), correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, RunAB, RunC, RunD, and crab_NZS2018C_1::runRange.

                                                                                                 {
  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);
}

correctForHLTDefinition() [13/13]

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

protected

Definition at line 289 of file EcalIsolationCorrector.cc.

References PixelPluginsPhase0_cfi::isBarrel, isElectron_, mps_fire::result, RunAB, RunC, RunD, and crab_NZS2018C_1::runRange.

                                                                                               {
  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;
}

correctForNoise() [1/13]

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 EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, 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;
}

correctForNoise() [2/13]

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

Definition at line 75 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, checkRunRange(), correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, isElectron_, convertSQLiteXML::runNumber, and crab_NZS2018C_1::runRange.

                                                                                           {
  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);
}

correctForNoise() [3/13]

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

Definition at line 85 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, isElectron_, RunAB, RunC, RunD, and crab_NZS2018C_1::runRange.

                                                                                               {
  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);
}

correctForNoise() [4/13]

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 EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, 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;
}

correctForNoise() [5/13]

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

Definition at line 166 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, checkRunRange(), correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, convertSQLiteXML::runNumber, and crab_NZS2018C_1::runRange.

                                                                                      {
  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);
}

correctForNoise() [6/13]

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

Definition at line 176 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, RunAB, RunC, RunD, and crab_NZS2018C_1::runRange.

                                                                                          {
  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);
}

correctForNoise() [7/13]

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 EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, 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;
}

correctForNoise() [8/13]

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

Definition at line 121 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, checkRunRange(), correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, isElectron_, convertSQLiteXML::runNumber, and crab_NZS2018C_1::runRange.

                                                                                       {
  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);
}

correctForNoise() [9/13]

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

Definition at line 131 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForNoise(), MillePedeFileConverter_cfg::e, HLTBitAnalyser_cfi::isData, isElectron_, RunAB, RunC, RunD, and crab_NZS2018C_1::runRange.

                                                                                           {
  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);
}

correctForNoise() [10/13]

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 EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, 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;
}

correctForNoise() [11/13]

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

Definition at line 211 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, checkRunRange(), correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, convertSQLiteXML::runNumber, and crab_NZS2018C_1::runRange.

                                                                                     {
  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);
}

correctForNoise() [12/13]

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

Definition at line 221 of file EcalIsolationCorrector.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, correctForNoise(), HLTBitAnalyser_cfi::isData, isElectron_, AlCaHLTBitMon_ParallelJobs::p, RunAB, RunC, RunD, and crab_NZS2018C_1::runRange.

                                                                                         {
  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);
}

correctForNoise() [13/13]

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, PixelPluginsPhase0_cfi::isBarrel, HLTBitAnalyser_cfi::isData, isElectron_, mps_fire::result, RunAB, RunC, RunD, and crab_NZS2018C_1::runRange.

                                                           {
  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

isElectron_

bool EcalIsolationCorrector::isElectron_

private

Definition at line 105 of file EcalIsolationCorrector.h.

Referenced by correctForHLTDefinition(), and correctForNoise().