egammaLayerRecalibTool Class Reference

#include <egammaLayerRecalibTool.h>

Inheritance diagram for egammaLayerRecalibTool:

Collaboration diagram for egammaLayerRecalibTool:

Public Types
typedef std::vector< std::pair< InputModifier *, GetAmountBase * > >	ModifiersList

Public Member Functions
	egammaLayerRecalibTool (const std::string &name, const std::string &tune, int SaccEnable=1)
	egammaLayerRecalibTool (const std::string &tune, int SaccEnable=1)
	~egammaLayerRecalibTool ()
CP::CorrectionCode	applyCorrection (xAOD::Egamma &, const xAOD::EventInfo &event_info) const
CP::CorrectionCode	scale_inputs (StdCalibrationInputs &inputs) const
	apply layer calibration to the More...
void	add_scale (InputModifier *modifier, GetAmountBase *amount)
	add custom layer scale correction. More...
void	add_scale (const std::string &scale)
	add scale correction from string. More...
void	clear_corrections ()
	remove all the scale corrections More...
void	fixForMissingCells (bool fix=true)
void	scaleMC (bool scaleMC=true)
void	disable_PSCorrections ()
void	disable_S12Corrections ()
void	disable_SaccCorrections ()
void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

Static Public Member Functions
static std::pair< std::string, egammaLayerRecalibTool * >	create (const std::string &type, const std::string &args)
	helper to create a tool from a string (useful for command line arguments) More...

Private Member Functions
const std::string	resolve_path (std::string filename) const
void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Static Private Member Functions
static std::string	resolve_alias (const std::string &tune)

Private Attributes
std::string	m_tune
bool	m_doPSCorrections = true
bool	m_doS12Corrections = true
bool	m_doSaccCorrections = true
ModifiersList	m_modifiers
corr_pileupShift *	m_pileup_tool = nullptr
bool	m_aodFixMissingCells = false
bool	m_scaleMC = false
std::string	m_nm
	Message source name. More...
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels) More...
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer. More...
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level. More...
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging) More...

Static Private Attributes
static const unsigned int	m_Run2Run3runNumberTransition = 400000

Detailed Description

Definition at line 282 of file egammaLayerRecalibTool.h.

Member Typedef Documentation

ModifiersList

typedef std::vector<std::pair<InputModifier*, GetAmountBase*> > egammaLayerRecalibTool::ModifiersList

Definition at line 285 of file egammaLayerRecalibTool.h.

Constructor & Destructor Documentation

egammaLayerRecalibTool() [1/2]

egammaLayerRecalibTool::egammaLayerRecalibTool	(	const std::string &	name,
		const std::string &	tune,
		int	SaccEnable = 1
	)

Parameters

tune	string to configure the tuning "" the tool has no effect as default it is "current_default" "test1" just for testing

Definition at line 980 of file egammaLayerRecalibTool.cxx.

  : asg::AsgMessaging(name), m_tune(tune), m_doSaccCorrections(SaccEnable)
{
  add_scale(tune);
}

egammaLayerRecalibTool() [2/2]

egammaLayerRecalibTool::egammaLayerRecalibTool	(	const std::string &	tune,
		int	SaccEnable = 1
	)

Definition at line 987 of file egammaLayerRecalibTool.cxx.

  : egammaLayerRecalibTool("egammaLayerRecalibTool", tune, SaccEnable) { }

~egammaLayerRecalibTool()

egammaLayerRecalibTool::~egammaLayerRecalibTool ( )

inline

Definition at line 294 of file egammaLayerRecalibTool.h.

{ clear_corrections(); delete m_pileup_tool; }

Member Function Documentation

add_scale() [1/2]

void egammaLayerRecalibTool::add_scale

(

const std::string &

scale

)

add scale correction from string.

Can be called multiple times. The list of valid values is on the twiki

Definition at line 245 of file egammaLayerRecalibTool.cxx.

{
  ATH_MSG_INFO("using scale " << tuneIn);
  std::string tune = resolve_alias(tuneIn);
 
  if (tune.empty()) { }
  else if ("es2024_run3_extrapolate_v0" == tune) {
    add_scale("run3_partial_ofc_extrapolate_v0");
  }
  // R22 layer tune with fixed E1E2 and repeated acc
  else if ("es2022_22.0_Precision_v1" == tune) {
    add_scale("run2_alt_with_layer2_r22_Precision_v1");
  }
  else if ("es2022_22.0_Precision" == tune) {
    add_scale("run2_alt_with_layer2_r22_Precision");
  }
  else if ("es2018_21.0_v0" == tune) {
    add_scale("run2_alt_with_layer2_r21_v1");
  }
  else if ("es2017_21.0_v0" == tune) {
    add_scale("run2_alt_with_layer2_r21_v0");
  }
  else if ("es2017_20.7_final" == tune) {
    add_scale("pileup_20.7");
    add_scale("run2_alt_with_layer2_modif");
  }
  else if ("es2017_20.7_improved" == tune) {
    add_scale("pileup_20.7");                // new pileup correction Guillaume for 20.7
    //TEMPORARY HACK REMOVED (two intermediate tags with this ES model use different layer corrections)
    add_scale("2012_alt_with_layer2_modif"); // temporary old corrections from run1 + EMECPS HV
  }
  else if ("pileup_20.7" == tune) {
    m_pileup_tool = new corr_pileupShift();
    add_scale(new ScaleE0(InputModifier::SUBTRACT), new GetAmountPileupE0(m_pileup_tool));
    add_scale(new ScaleE1(InputModifier::SUBTRACT), new GetAmountPileupE1(m_pileup_tool));
    add_scale(new ScaleE2(InputModifier::SUBTRACT), new GetAmountPileupE2(m_pileup_tool));
    add_scale(new ScaleE3(InputModifier::SUBTRACT), new GetAmountPileupE3(m_pileup_tool));
  }
  // Run3 2022+2023
  else if ("run3_partial_ofc_extrapolate_v0" == tune) {
    add_scale("layer2_run3_ofc_extrapolate_v0");
    add_scale("ps_run3_ofc_extrapolate_v0");
    if(m_doSaccCorrections) add_scale("acc_zee_run3_v0");
  }
  //Run 2 release 22 with fixed E1E2 and repeated acc 
  else if ("run2_alt_with_layer2_r22_Precision_v1"==tune) {
    add_scale("layer2_alt_el_mu_comb_r21_v0_fix");
    add_scale("ps_mu_r21_v0");
    if(m_doSaccCorrections) add_scale("acc_zee_r22_v1");
  }
  //Run 2 release 22
  else if ("run2_alt_with_layer2_r22_Precision"==tune) {
    add_scale("layer2_alt_el_mu_comb_r21_v0");
    add_scale("ps_mu_r21_v0");
    if(m_doSaccCorrections) add_scale("acc_zee_r22_v0");
  }
  else if ("run2_alt_with_layer2_r21_v1"==tune) {
    add_scale("layer2_alt_run2_r21_v1");
    add_scale("ps_2016_r21_v0");
  }
  else if ("run2_alt_with_layer2_r21_v0"==tune) {
    add_scale("layer2_alt_run2_r21_v0");
    add_scale("ps_2016_r21_v0");
  }
  else if("run2_alt_with_layer2_modif" == tune) {
    add_scale("ps_EMECHV1");
    add_scale("layer2_alt_run2_v1");
    add_scale("ps_2016");
  }
  // 2012
  else if ("2012" == tune) {
    add_scale("ps_HV1");
    add_scale("layer1_2012");
    add_scale("ps_2012");
  }
  else if("2012_with_layer2" == tune) {
    add_scale("ps_HV1");
    add_scale("layer2_2012_v5");
    add_scale("ps_2012");
  }
  else if ("2012_alt" == tune) {
    add_scale("ps_HV1");
    add_scale("layer1_alt_2012");
    add_scale("ps_2012");
  }
  else if("2012_alt_with_layer2" == tune) {
    add_scale("ps_HV1");
    add_scale("layer2_alt_2012_v5");
    add_scale("ps_2012");
  }
  else if("2012_alt_with_layer2_modif" == tune) {
    add_scale("ps_HV1");
    add_scale("ps_EMECHV1");
    add_scale("layer2_alt_2012_v5");
    add_scale("ps_2012");
  }
  else if("2010_with_layer2" == tune) {
    add_scale("layer2_2010_v5");
    add_scale("ps_2010");
  }
  else if ("2012_layer1_up" == tune) {
    add_scale("ps_HV1");
    add_scale("layer1_2012_up");
    add_scale("ps_2012");
  }
  else if ("2012_layer1_down" == tune) {
    add_scale("ps_HV1");
    add_scale("layer1_2012_down");
    add_scale("ps_2012");
  }
  else if ("2012_layer1_errup" == tune) {
    add_scale("layer1_2012_errup");
  }
  else if ("2012_layer1_errdown" == tune) {
    add_scale("layer1_2012_errdown");
  }
  else if ("2012_ps_down" == tune) {
    add_scale("ps_HV1");
    add_scale("layer1_2012");
    add_scale("ps_2012_down");
  }
  else if ("2012_ps_up" == tune) {
    add_scale("ps_HV1");
    add_scale("layer1_2012");
    add_scale("ps_2012_up");
  }
  else if ("2012_ps_errdown" == tune) {
    add_scale("ps_2012_errdown");
  }
  else if ("2012_ps_errup" == tune) {
    add_scale("ps_2012_errup");
  }
  else if ("2012_up" == tune) {
    add_scale("ps_HV1");
    add_scale("layer1_2012_up");
    add_scale("ps_2012_up");
  }
  else if ("2012_down" == tune) {
    add_scale("ps_HV1");
    add_scale("layer1_2012_down");
    add_scale("ps_2012_down");
  }
  else if ("2012_errup" == tune) {
    add_scale("layer1_2012_errup");
    add_scale("ps_2012_errup");
  }
  else if ("2012_errdown" == tune) {
    add_scale("layer1_2012_errdown");
    add_scale("ps_2012_errdown");
  }
  // 2011
  else if ("2011" == tune) {
    add_scale("layer1_2011");
    add_scale("ps_2011");
  }
  else if("2011_with_layer2" == tune) {
    add_scale("layer2_2011_v5");
    add_scale("ps_2011");
  }
  else if ("2011_alt" == tune) {
    add_scale("layer1_alt_2011");
    add_scale("ps_2011");
  }
  else if("2011_alt_with_layer2" == tune) {
    add_scale("layer2_alt_2011_v5");
    add_scale("ps_2011");
  }
  else if ("2011_layer1_up" == tune) {
    add_scale("layer1_2011_up");
    add_scale("ps_2011");
  }
  else if ("2011_layer1_down" == tune) {
    add_scale("layer1_2011_down");
    add_scale("ps_2011");
  }
  else if ("2011_layer1_errup" == tune) {
    add_scale("layer1_2011_errup");
  }
  else if ("2011_layer1_errdown" == tune) {
    add_scale("layer1_2011_errdown");
  }
  else if ("2011_ps_down" == tune) {
    add_scale("layer1_2011");
    add_scale("ps_2011_down");
  }
  else if ("2011_ps_up" == tune) {
    add_scale("layer1_2011");
    add_scale("ps_2011_up");
  }
  else if ("2011_ps_errdown" == tune) {
    add_scale("ps_2011_errdown");
  }
  else if ("2011_ps_errup" == tune) {
    add_scale("ps_2011_errup");
  }
  else if ("2011_up" == tune) {
    add_scale("layer1_2011_up");
    add_scale("ps_2011_up");
  }
  else if ("2011_down" == tune) {
    add_scale("layer1_2011_down");
    add_scale("ps_2011_down");
  }
  else if ("2011_errup" == tune) {
    add_scale("layer1_2011_errup");
    add_scale("ps_2011_errup");
  }
  else if ("2011_errdown" == tune) {
    add_scale("layer1_2011_errdown");
    add_scale("ps_2011_errdown");
  }
  // 2010
  else if ("2010" == tune) {
    add_scale("layer1_2010");
    add_scale("ps_2010");
  }
  else if ("2010_layer1_up" == tune) {
    add_scale("layer1_2010_up");
    add_scale("ps_2010");
  }
  else if ("2010_layer1_down" == tune) {
    add_scale("layer1_2010_down");
    add_scale("ps_2010");
  }
  else if ("2010_layer1_errup" == tune) {
    add_scale("layer1_2010_errup");
  }
  else if ("2010_layer1_errdown" == tune) {
    add_scale("layer1_2010_errdown");
  }
  else if ("2010_ps_down" == tune) {
    add_scale("layer1_2010");
    add_scale("ps_2010_down");
  }
  else if ("2010_ps_up" == tune) {
    add_scale("layer1_2010");
    add_scale("ps_2010_up");
  }
  else if ("2010_ps_errdown" == tune) {
    add_scale("ps_2010_errdown");
  }
  else if ("2010_ps_errup" == tune) {
    add_scale("ps_2010_errup");
  }
  else if ("2010_up" == tune) {
    add_scale("layer1_2010_up");
    add_scale("ps_2010_up");
  }
  else if ("2010_down" == tune) {
    add_scale("layer1_2010_down");
    add_scale("ps_2010_down");
  }
  else if ("2010_errup" == tune) {
    add_scale("layer1_2010_errup");
    add_scale("ps_2010_errup");
  }
  else if ("2010_errdown" == tune) {
    add_scale("layer1_2010_errdown");
    add_scale("ps_2010_errdown");
  }
  else if ("ps_HV1" == tune) {
    add_scale(new ScaleE0(InputModifier::ONEBASED), new GetAmountHVPSGuillaume());
  }
  else if ("ps_EMECHV1" == tune) {
    add_scale(new ScaleE0(InputModifier::ONEBASED), new GetAmountHVEMECPS207());
  }
  else if ("test1" == tune) {
    TH1F h_presampler("h_presampler", "h_presampler", 10, -2.5, 2.5);
    // just as an example, correct E0 by 0.1 * sign(eta)
    // and E1 by 1%
    for (int ibin = 1; ibin <= 5; ++ibin) {
      h_presampler.SetBinContent(ibin, -0.1);
      h_presampler.SetBinContent(ibin + 5, 0.1);
    }
    add_scale(new ScaleE0(InputModifier::ZEROBASED), new GetAmountHisto1D(h_presampler));
    add_scale(new ScaleE1(InputModifier::ZEROBASED), new GetAmountFixed(0.01));
  }
  else if ("acc_zee_run3_v0" == tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v13/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH2F* histo_acc = static_cast<TH2F*>(f.Get("hACC_Zee_rel23"));
    assert(histo_acc);
    add_scale(new ScaleEaccordion(InputModifier::ZEROBASED_ALPHA),
              new GetAmountHisto2DEtaCaloRunNumber(*histo_acc));    
  }
  // repeated acc scale based on layer2_alt_el_mu_comb_r21_v0_fix 
  else if ("acc_zee_r22_v1" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v12/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH2F* histo_acc = static_cast<TH2F*>(f.Get("hACC_Zee_rel22"));
    assert(histo_acc);
    add_scale(new ScaleEaccordion(InputModifier::ZEROBASED_ALPHA),
              new GetAmountHisto2DEtaCaloRunNumber(*histo_acc));
  }
  else if ("acc_zee_r22_v0" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v11/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH2F* histo_acc = static_cast<TH2F*>(f.Get("hACC_Zee_rel22"));
    assert(histo_acc);
    add_scale(new ScaleEaccordion(InputModifier::ZEROBASED_ALPHA),
              new GetAmountHisto2DEtaCaloRunNumber(*histo_acc));
  }
  else if ("layer1_1" == tune) {
    TFormula f("formula_layer1_1", "(abs(x)<1.425) ? 0.97 : 1");
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA), new GetAmountFormula(f));
  }
  else if ("layer1_2" == tune) {
    TFormula f("formula_layer1_2", "(abs(x)<1.425) ? 0.97 : 1.05");
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA), new GetAmountFormula(f));
  }
  else if ("layer1_alt_2012_v5" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_alt_2012"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1D(*histo));
  }
  else if ("layer1_2012_v5" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2012"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1D(*histo));
  }
  else if ("layer1_2012_v5_down" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2012"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DUp(*histo));
  }
  else if ("layer1_2012_v5_up" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2012"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DDown(*histo));
  }
  else if ("layer1_2012_v5_errdown" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2012"));
    add_scale(new ScaleE1(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorUp(*histo));
  }
  else if ("layer1_2012_v5_errup" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2012"));
    add_scale(new ScaleE1(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorDown(*histo));
  }
  else if ("layer1_alt_2011_v5" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_alt_2011"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1D(*histo));
  }
  else if ("layer1_2011_v5" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2011"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1D(*histo));
  }
  else if ("layer1_2011_v5_down" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2011"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DUp(*histo));
  }
  else if ("layer1_2011_v5_up" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2011"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DDown(*histo));
  }
  else if ("layer1_2011_v5_errdown" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2011"));
    add_scale(new ScaleE1(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorUp(*histo));
  }
  else if ("layer1_2011_v5_errup" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2011"));
    add_scale(new ScaleE1(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorDown(*histo));
  }
  else if ("layer1_2010_v5" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2010"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1D(*histo));
  }
  else if ("layer1_2010_v5_down" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2010"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DUp(*histo));
  }
  else if ("layer1_2010_v5_up" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2010"));
    add_scale(new ScaleE1(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DDown(*histo));
  }
  else if ("layer1_2010_v5_errdown" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2010"));
    add_scale(new ScaleE1(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorUp(*histo));
  }
  else if ("layer1_2010_v5_errup" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2010"));
    add_scale(new ScaleE1(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorDown(*histo));
  }
  else if ("layer2_run3_ofc_extrapolate_v0"==tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v12/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1D* histo = static_cast<TH1D*>(f.Get("hE1E2_emu_run2_rel21_v1_run3ofc"));
    assert(histo);
    add_scale(new ScaleE2(InputModifier::ONEBASED),
         new GetAmountHisto1D(*histo));    
  }
  // fix E1E2 scale from R21 precision model
  else if("layer2_alt_el_mu_comb_r21_v0_fix"==tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v12/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1D* histo = static_cast<TH1D*>(f.Get("hE1E2_emu_run2_rel21_v0_fix"));
    assert(histo);
    add_scale(new ScaleE2(InputModifier::ONEBASED),
         new GetAmountHisto1D(*histo));
  }
  else if("layer2_alt_el_mu_comb_r21_v0"==tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v11/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1D* histo = static_cast<TH1D*>(f.Get("hE1E2_emu_run2_rel21_v0"));
    assert(histo);
    add_scale(new ScaleE2(InputModifier::ONEBASED),
         new GetAmountHisto1D(*histo));
  }
  else if("layer2_alt_run2_r21_v1"==tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v6/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2_mu_run2_rel21_v1"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
              new GetAmountHisto1D(*histo));
  }
  else if("layer2_alt_run2_r21_v0"==tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v5/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2mu_2016_rel21_v1"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
              new GetAmountHisto1D(*histo));
  }
  else if("layer2_alt_run2_v1" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v3/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2mu_2016_v1"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
              new GetAmountHisto1D(*histo));
  }
  else if("layer2_alt_2012_v5" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_alt_2012"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1D(*histo));
  }
  else if("layer2_2012_v5" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2012"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1D(*histo));
  }
  else if("layer2_2012_v5_down" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2012"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DDown(*histo));
  }
  else if("layer2_2012_v5_up" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2012"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DUp(*histo));
  }
  else if ("layer2_2012_v5_errdown" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2012"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DErrorDown(*histo));
  }
  else if ("layer2_2012_v5_errup" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2012"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DErrorUp(*histo));
  }
  else if("layer2_alt_2011_v5" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_alt_2011"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1D(*histo));
  }
  else if("layer2_2011_v5" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2011"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1D(*histo));
  }
  else if("layer2_2011_v5_down" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2011"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DDown(*histo));
  }
  else if("layer2_2011_v5_up" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2011"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DUp(*histo));
  }
  else if ("layer2_2011_v5_errdown" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2011"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DErrorDown(*histo));
  }
  else if ("layer2_2011_v5_errup" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2011"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DErrorUp(*histo));
  }
  else if("layer2_2010_v5" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2010"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1D(*histo));
  }
  else if("layer2_2010_v5_down" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2010"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DDown(*histo));
  }
  else if("layer2_2010_v5_up" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2010"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DUp(*histo));
  }
  else if ("layer2_2010_v5_errdown" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2010"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DErrorDown(*histo));
  }
  else if ("layer2_2010_v5_errup" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo = checked_cast<TH1*>(f.Get("hE1E2ave_2010"));
    add_scale(new ScaleE2(InputModifier::ONEBASED),
          new GetAmountHisto1DErrorUp(*histo));
  }
  else if ("ps_2016_r21_v0" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v5/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2016_rel21"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
              new GetAmountHisto1D(*histo_ps_tot_error));
  }
  else if ("ps_run3_ofc_extrapolate_v0" == tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v12/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1F* histo_ps_tot_error = static_cast<TH1F*>(f.Get("hPS_MuonLowMu_rel21_run3ofc"));
    assert(histo_ps_tot_error);
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
              new GetAmountHisto1D(*histo_ps_tot_error));
  }
  else if ("ps_mu_r21_v0" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v11/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1F* histo_ps_tot_error = static_cast<TH1F*>(f.Get("hPS_MuonLowMu_rel21"));
    assert(histo_ps_tot_error);
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
              new GetAmountHisto1D(*histo_ps_tot_error));
  }
  else if ("ps_2016_v1" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v4/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2016"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1D(*histo_ps_tot_error));
  }
  else if ("ps_2012_v3" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2012"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1D(*histo_ps_tot_error));
  }
  else if ("ps_2012_v3_down" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2012"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DUp(*histo_ps_tot_error));
  }
  else if ("ps_2012_v3_up" == tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2012"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DDown(*histo_ps_tot_error));
  }
  else if ("ps_2012_v3_errdown" == tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2012"));
    add_scale(new ScaleE0(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorUp(*histo_ps_tot_error));
  }
  else if ("ps_2012_v3_errup" == tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2012"));
    add_scale(new ScaleE0(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorDown(*histo_ps_tot_error));
  }
  else if ("ps_2011_v3" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2011"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1D(*histo_ps_tot_error));
  }
  else if ("ps_2011_v3_down" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2011"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DUp(*histo_ps_tot_error));
  }
  else if ("ps_2011_v3_up" == tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2011"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DDown(*histo_ps_tot_error));
  }
  else if ("ps_2011_v3_errdown" == tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2011"));
    add_scale(new ScaleE0(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorUp(*histo_ps_tot_error));
  }
  else if ("ps_2011_v3_errup" == tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2011"));
    add_scale(new ScaleE0(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorDown(*histo_ps_tot_error));
  }
  // 2010
  else if ("ps_2010_v3" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2010"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1D(*histo_ps_tot_error));
  }
  else if ("ps_2010_v3_down" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2010"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DUp(*histo_ps_tot_error));
  }
  else if ("ps_2010_v3_up" == tune) {
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2010"));
    add_scale(new ScaleE0(InputModifier::ONEBASED_ALPHA),
          new GetAmountHisto1DDown(*histo_ps_tot_error));
  }
  else if ("ps_2010_v3_errdown" == tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2010"));
    add_scale(new ScaleE0(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorUp(*histo_ps_tot_error));
  }
  else if ("ps_2010_v3_errup" == tune){
    const std::string file = PathResolverFindCalibFile("egammaLayerRecalibTool/v1/egammaLayerRecalibTunes.root");
    TFile f(file.c_str());
    TH1* histo_ps_tot_error = checked_cast<TH1*>(f.Get("hPS_2010"));
    add_scale(new ScaleE0(InputModifier::ZEROBASED_ALPHA),
          new GetAmountHisto1DErrorDown(*histo_ps_tot_error));
  }
  else {
    throw std::runtime_error(tune+" is not a valid tune");
  }
}

add_scale() [2/2]

void egammaLayerRecalibTool::add_scale	(	InputModifier *	modifier,
		GetAmountBase *	amount
	)

add custom layer scale correction.

Can be called multiple times.

Definition at line 991 of file egammaLayerRecalibTool.cxx.

{
  m_modifiers.emplace_back(modifier, amount);
}

applyCorrection()

CP::CorrectionCode egammaLayerRecalibTool::applyCorrection	(	xAOD::Egamma &	particle,
		const xAOD::EventInfo &	event_info
	)		const

Definition at line 1010 of file egammaLayerRecalibTool.cxx.

{
  const xAOD::CaloCluster* cluster = particle.caloCluster();
  if (!cluster) {
    ATH_MSG_ERROR("egamma particle without CaloCluster");
    return CP::CorrectionCode::Error;
  }
 
  std::string fixT = "";
  double addE2 = 0, addE3 = 0;
  if (m_aodFixMissingCells &&
      event_info.runNumber() > m_Run2Run3runNumberTransition) {
    fixT = "_egFixForTopoTimingCut";
    unsigned short stat =
      xAOD::EgammaHelpers::energyInMissingCells(particle,addE2,addE3);
    if (stat) {
      ATH_MSG_WARNING("Fix for missing cells required"
              " but some layer info is not available,"
              " from L2 : " << stat%2 << " from L3 : " << stat/2);
    }
  }
 
  double eta_calo;
  static const SG::AuxElement::Accessor<float> accEtaCalo("etaCalo");
  if(particle.author() == xAOD::EgammaParameters::AuthorFwdElectron){
    eta_calo = cluster->eta();
  }
  else if (cluster->retrieveMoment(xAOD::CaloCluster::ETACALOFRAME, eta_calo)){
 
  }
  else if (accEtaCalo.isAvailable(*cluster)) {
    eta_calo = accEtaCalo(*cluster);
  }
  else{
    ATH_MSG_ERROR("etaCalo not available as auxilliary variable,"
          " using cluster eta as eta calo!");
    eta_calo=cluster->eta();
  }
 
  StdCalibrationInputs inputs {
    event_info.averageInteractionsPerCrossing(),
      event_info.runNumber(),
      cluster->eta(),
      cluster->phi(),
      cluster->energyBE(0),
      cluster->energyBE(1),
      cluster->energyBE(2) + addE2,
      cluster->energyBE(3) + addE3,
      eta_calo };
 
  bool isData = !event_info.eventType(xAOD::EventInfo::IS_SIMULATION);
  CP::CorrectionCode status = CP::CorrectionCode::Ok;
  if (isData || m_scaleMC)
    status = scale_inputs(inputs);
 
  static const SG::AuxElement::Decorator<double> deco_E0("correctedcl_Es0");
  static const SG::AuxElement::Decorator<double> deco_E1("correctedcl_Es1");
  static const SG::AuxElement::Decorator<double> deco_E2("correctedcl_Es2");
  static const SG::AuxElement::Decorator<double> deco_E3("correctedcl_Es3");
  static const SG::AuxElement::Decorator<std::string>
    deco_layer_correction("layer_correction");
 
  if (status == CP::CorrectionCode::Ok) {
    ATH_MSG_DEBUG("decorating cluster with corrected layer energies");
    deco_E0(*cluster) = m_doPSCorrections  ?
      inputs.E0raw : cluster->energyBE(0);
    deco_E1(*cluster) = m_doS12Corrections or m_doSaccCorrections ?
      inputs.E1raw : cluster->energyBE(1);
    deco_E2(*cluster) = m_doS12Corrections or m_doSaccCorrections ?
      inputs.E2raw : cluster->energyBE(2) + addE2;
    deco_E3(*cluster) = m_doSaccCorrections ?
      inputs.E3raw : cluster->energyBE(3) + addE3;
    deco_layer_correction(*cluster) = isData ? m_tune+fixT : fixT;
    return status;
  }
 
  ATH_MSG_DEBUG("cannot correct layer energies:"
        " decorating particle with non-corrected layer energies");
  // this is done for safety, since when a particle is decorated
  // all the particle in the container are decorated
  // it is not possible to distinguish between decorated / non-decorated
  // since all are decorated
  deco_E0(*cluster) = cluster->energyBE(0);
  deco_E1(*cluster) = cluster->energyBE(1);
  deco_E2(*cluster) = cluster->energyBE(2) + addE2;
  deco_E3(*cluster) = cluster->energyBE(3) + addE3;
  deco_layer_correction(*cluster) = isData ? m_tune + "_Err" + fixT : fixT;
  return status;
 
}

clear_corrections()

void egammaLayerRecalibTool::clear_corrections

(

)

remove all the scale corrections

Definition at line 1102 of file egammaLayerRecalibTool.cxx.

{
  for (auto modifier :  m_modifiers) {
    delete modifier.first;
    delete modifier.second;
  }
  m_modifiers.clear();
}

create()

std::pair< std::string, egammaLayerRecalibTool * > egammaLayerRecalibTool::create ( const std::string &  type, const std::string &  args  )

static

helper to create a tool from a string (useful for command line arguments)

Definition at line 1147 of file egammaLayerRecalibTool.cxx.

{
  std::map<std::string, std::string> args_map = parse(args);
  egammaLayerRecalibTool* tool = new egammaLayerRecalibTool("egammaLayerRecalibTool", "");
  std::string name = "";
  std::string amount_name = "";
  std::string type_name = "";
 
  GetAmountBase* amount_getter = nullptr;
  InputModifier* modifier = nullptr;
 
  if (args_map.find("amount") != args_map.end()) {
    std::string amount_str = args_map["amount"];
    bool perc = false;
    if (amount_str.back()=='%') {
      perc = true;
      amount_str.pop_back();
    }
    const float amount = TString(amount_str).Atof() * (perc ? 0.01 : 1);
 
    amount_getter = new GetAmountFixed(amount);
    std::stringstream amount_stream;
    amount_stream << amount;
    amount_name = amount_stream.str();
    std::replace(amount_name.begin(), amount_name.end(), '-', 'n');
    std::replace(amount_name.begin(), amount_name.end(), '.', 'p');
  }
  else if (args_map.find("name") != args_map.end()) {
    name = args_map["name"];
  }
  else if (args_map.find("histo") != args_map.end()) {
    int dim = 0;
    if (args_map.find("file") == args_map.end()) {
      std::cerr << "with histo you must specify file" << std::endl;
      assert(false);
    }
    if (args_map.find("formulax") != args_map.end()) dim = 1;
 
    if (dim == 0)
      {
    std::cerr << "with histo you must specify formulax" << std::endl;
    assert(false);
      }
    if (dim == 1) {
      TFile f(args_map["file"].c_str());
      std::cout << "opening histo " << args_map["histo"] << " from file " << args_map["file"] << std::endl;
      TH1F* histo = dynamic_cast<TH1F*>(f.Get(args_map["histo"].c_str()));
 
      if(histo){
    histo->SetDirectory(nullptr);
    amount_getter = new GetAmountHisto1D(*histo);
      }
      else{assert(false); }
    }
    else { assert(false); }
  }
  else {
    std::cerr << "cannot understand argument " << args << std::endl;
    assert(false);
  }
 
  if ("bias-E0" == type)                { modifier = new ScaleE0(InputModifier::ZEROBASED); type_name = "E0"; }
  else if ("bias-E1" == type)           { modifier = new ScaleE1(InputModifier::ZEROBASED); type_name = "E1"; }
  else if ("bias-E2" == type)           { modifier = new ScaleE2(InputModifier::ZEROBASED); type_name = "E2"; }
  else if ("bias-E3" == type)           { modifier = new ScaleE3(InputModifier::ZEROBASED); type_name = "E3"; }
  else if ("bias-E1overE2" == type)     { modifier = new ScaleE1overE2(InputModifier::ZEROBASED); type_name = "E1overE2"; }
  else if ("bias-Eaccordion" == type)   { modifier = new ScaleEaccordion(InputModifier::ZEROBASED); type_name = "Eaccordion"; }
  else if ("bias-Ecalorimeter" == type) { modifier = new ScaleEcalorimeter(InputModifier::ZEROBASED); type_name = "Ecalorimeter"; }
 
  if (not type_name.empty() and not amount_name.empty()) {
    name = type_name + "_" + amount_name;
  }
 
  if (name.empty()) {
    std::cerr << "you need to specify a name for the bias with type " << type << std::endl;
  }
 
  if (modifier and amount_getter) {
    tool->add_scale(modifier, amount_getter);
  }
  else{
    tool->add_scale(type);
    //release resources, if modifier false need to release amount_getter and vice versa
    //since they are not passed to the tool
    if(modifier) delete modifier;
    if(amount_getter) delete  amount_getter;
  }
 
  return {name, tool};
}

disable_PSCorrections()

void egammaLayerRecalibTool::disable_PSCorrections ( )

inline

Definition at line 324 of file egammaLayerRecalibTool.h.

{m_doPSCorrections=false;}

disable_S12Corrections()

void egammaLayerRecalibTool::disable_S12Corrections ( )

inline

Definition at line 325 of file egammaLayerRecalibTool.h.

{m_doS12Corrections=false;}

disable_SaccCorrections()

void egammaLayerRecalibTool::disable_SaccCorrections ( )

inline

Definition at line 326 of file egammaLayerRecalibTool.h.

{m_doSaccCorrections=false;}

fixForMissingCells()

void egammaLayerRecalibTool::fixForMissingCells ( bool  fix = true )

inline

Definition at line 322 of file egammaLayerRecalibTool.h.

{ m_aodFixMissingCells = fix; }

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  m_lvl = m_imsg ?
    static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
    MSG::INFO;
}

msg() [1/2]

MsgStream & asg::AsgMessaging::msg ( ) const

inherited

The standard message stream.

Returns

A reference to the default message stream of this object.

Definition at line 49 of file AsgMessaging.cxx.

                                      {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msg();
#else // not XAOD_STANDALONE
      return m_msg;
#endif // not XAOD_STANDALONE
   }

msg() [2/2]

MsgStream & asg::AsgMessaging::msg ( const MSG::Level  lvl ) const

inherited

The standard message stream.

Parameters

lvl	The message level to set the stream to

Returns

A reference to the default message stream, set to level "lvl"

Definition at line 57 of file AsgMessaging.cxx.

                                                          {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msg( lvl );
#else // not XAOD_STANDALONE
      m_msg << lvl;
      return m_msg;
#endif // not XAOD_STANDALONE
   }

msgLvl()

bool asg::AsgMessaging::msgLvl ( const MSG::Level  lvl ) const

inherited

Test the output level of the object.

Parameters

lvl	The message level to test against

Returns

boolean Indicting if messages at given level will be printed

true If messages at level "lvl" will be printed

Definition at line 41 of file AsgMessaging.cxx.

                                                       {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msgLvl( lvl );
#else // not XAOD_STANDALONE
      return m_msg.msgLevel( lvl );
#endif // not XAOD_STANDALONE
   }

resolve_alias()

std::string egammaLayerRecalibTool::resolve_alias ( const std::string &  tune )

staticprivate

Definition at line 207 of file egammaLayerRecalibTool.cxx.

                                                                     {
 
  if ("layer1_2012" == tune)         return "layer1_2012_v5";
  if ("layer1_alt_2012" == tune)     return "layer1_alt_2012_v5";
  if ("layer1_2011" == tune)         return "layer1_2011_v5";
  if ("layer1_alt_2011" == tune)     return "layer1_alt_2011_v5";
  if ("layer1_2010" == tune)         return "layer1_2010_v5";
  if ("ps_2012" == tune)             return "ps_2012_v3";
  if ("ps_2011" == tune)             return "ps_2011_v3";
  if ("ps_2010" == tune)             return "ps_2010_v3";
  if ("layer1_2012_up" == tune)      return "layer1_2012_v5_up";
  if ("layer1_2012_down" == tune)    return "layer1_2012_v5_down";
  if ("layer1_2012_errup" == tune)   return "layer1_2012_v5_errup";
  if ("layer1_2012_errdown" == tune) return "layer1_2012_v5_errdown";
  if ("layer1_2011_up" == tune)      return "layer1_2011_v5_up";
  if ("layer1_2011_down" == tune)    return "layer1_2011_v5_down";
  if ("layer1_2011_errup" == tune)   return "layer1_2011_v5_errup";
  if ("layer1_2011_errdown" == tune) return "layer1_2011_v5_errdown";
  if ("layer1_2010_up" == tune)      return "layer1_2010_v5_up";
  if ("layer1_2010_down" == tune)    return "layer1_2010_v5_down";
  if ("layer1_2010_errup" == tune)   return "layer1_2010_v5_errup";
  if ("layer1_2010_errdown" == tune) return "layer1_2010_v5_errdown";
  if ("ps_2012_up" == tune)          return "ps_2012_v3_up";
  if ("ps_2012_down" == tune)        return "ps_2012_v3_down";
  if ("ps_2012_errup" == tune)       return "ps_2012_v3_errup";
  if ("ps_2012_errdown" == tune)     return "ps_2012_v3_errdown";
  if ("ps_2011_up" == tune)          return "ps_2011_v3_up";
  if ("ps_2011_down" == tune)        return "ps_2011_v3_down";
  if ("ps_2011_errup" == tune)       return "ps_2011_v3_errup";
  if ("ps_2011_errdown" == tune)     return "ps_2011_v3_errdown";
  if ("ps_2010_up" == tune)          return "ps_2010_v3_up";
  if ("ps_2010_down" == tune)        return "ps_2010_v3_down";
  if ("ps_2010_errup" == tune)       return "ps_2010_v3_errup";
  if ("ps_2010_errdown" == tune)     return "ps_2010_v3_errdown";
 
  return tune;
}

resolve_path()

const std::string egammaLayerRecalibTool::resolve_path ( std::string  filename ) const

private

scale_inputs()

CP::CorrectionCode egammaLayerRecalibTool::scale_inputs

(

StdCalibrationInputs &

inputs

)

const

apply layer calibration to the

Parameters

inputs

Definition at line 996 of file egammaLayerRecalibTool.cxx.

{
  CP::CorrectionCode status = CP::CorrectionCode::Ok;
  for (const auto& modifier : m_modifiers)  {
    const float amount = (*modifier.second)(inputs);
    const auto s = (*modifier.first)(inputs, amount);
    ATH_MSG_DEBUG("  after E0|E1|E2|E3 = " << inputs.E0raw << "|" << inputs.E1raw << "|" << inputs.E2raw << "|" << inputs.E3raw);
    if (s != CP::CorrectionCode::Ok) {
      if (status != CP::CorrectionCode::Error) { status = s; }
    }
  }
  return status;
}

scaleMC()

void egammaLayerRecalibTool::scaleMC ( bool  scaleMC = true )

inline

Definition at line 323 of file egammaLayerRecalibTool.h.

{ m_scaleMC = scaleMC; }

setLevel()

void AthMessaging::setLevel ( MSG::Level  lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

m_aodFixMissingCells

bool egammaLayerRecalibTool::m_aodFixMissingCells = false

private

Definition at line 342 of file egammaLayerRecalibTool.h.

m_doPSCorrections

bool egammaLayerRecalibTool::m_doPSCorrections = true

private

Definition at line 333 of file egammaLayerRecalibTool.h.

m_doS12Corrections

bool egammaLayerRecalibTool::m_doS12Corrections = true

private

Definition at line 334 of file egammaLayerRecalibTool.h.

m_doSaccCorrections

bool egammaLayerRecalibTool::m_doSaccCorrections = true

private

Definition at line 335 of file egammaLayerRecalibTool.h.

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

m_modifiers

ModifiersList egammaLayerRecalibTool::m_modifiers

private

Definition at line 338 of file egammaLayerRecalibTool.h.

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_pileup_tool

corr_pileupShift* egammaLayerRecalibTool::m_pileup_tool = nullptr

private

Definition at line 340 of file egammaLayerRecalibTool.h.

m_Run2Run3runNumberTransition

const unsigned int egammaLayerRecalibTool::m_Run2Run3runNumberTransition = 400000

staticprivate

Definition at line 330 of file egammaLayerRecalibTool.h.

m_scaleMC

bool egammaLayerRecalibTool::m_scaleMC = false

private

Definition at line 343 of file egammaLayerRecalibTool.h.

m_tune

std::string egammaLayerRecalibTool::m_tune

private

Definition at line 332 of file egammaLayerRecalibTool.h.

---

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

• egammaLayerRecalibTool.h
• egammaLayerRecalibTool.cxx