#include <PerJetFlavourUncertaintyComponent.h>

Inheritance diagram for jet::PerJetFlavourUncertaintyComponent:

Collaboration diagram for jet::PerJetFlavourUncertaintyComponent:

Public Member Functions
	PerJetFlavourUncertaintyComponent (const ComponentHelper &component)

	PerJetFlavourUncertaintyComponent (const PerJetFlavourUncertaintyComponent &toCopy)

virtual PerJetFlavourUncertaintyComponent *	clone () const

virtual	~PerJetFlavourUncertaintyComponent ()

virtual StatusCode	initialize (TFile *histFile)

virtual std::vector< int >	getRelevantLabels () const

virtual TString	getName () const

virtual TString	getValidName () const

virtual CompScaleVar::TypeEnum	getScaleVar () const

virtual JetTopology::TypeEnum	getTopology () const

virtual bool	isAlwaysZero () const

virtual bool	getValidity (const xAOD::Jet &jet, const xAOD::EventInfo &eInfo) const

virtual double	getUncertainty (const xAOD::Jet &jet, const xAOD::EventInfo &eInfo) const

virtual bool	getValidUncertainty (double &unc, const xAOD::Jet &jet, const xAOD::EventInfo &eInfo) const

void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

Protected Member Functions
virtual bool	getValidityImpl (const xAOD::Jet &jet, const xAOD::EventInfo &eInfo) const

virtual double	getUncertaintyImpl (const xAOD::Jet &jet, const xAOD::EventInfo &eInfo) const

virtual bool	getValidBool (const double validity) const

virtual double	getSplitFactor (const xAOD::Jet &jet) const

virtual double	getAbsMass (const xAOD::Jet &jet, const CompMassDef::TypeEnum massDef) const

virtual double	getMassOverPt (const xAOD::Jet &jet, const CompMassDef::TypeEnum massDef) const

virtual double	getMassOverE (const xAOD::Jet &jet, const CompMassDef::TypeEnum massDef) const

Protected Attributes
bool	m_isInit {}

const TString	m_uncHistName

const TString	m_validHistName

const CompScaleVar::TypeEnum	m_scaleVar

const JetTopology::TypeEnum	m_topology

const float	m_energyScale {}

const Interpolate::TypeEnum	m_interpolate

const int	m_splitNumber {}

int	m_numExpectedHist {}

UncertaintyHistogram *	m_uncHist {}

UncertaintyHistogram *	m_validHist {}

Private Member Functions
	PerJetFlavourUncertaintyComponent (const std::string &name="")

double	getFlavourResponseUncertainty (const xAOD::Jet &jet, const xAOD::EventInfo &eInfo) const

bool	isSupportedLabel (const int label) const

bool	checkTruthLabel (const xAOD::Jet &jet) const

void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Private Attributes
const bool	m_absEta {}

const std::vector< int >	m_labels

const FlavourComp::TypeEnum	m_flavourType {}

const bool	m_constrainZresponse {}

const TString	m_constrainZresponseFunc {}

const TF1 *	m_ZjetQuarkFrac {}

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...

Detailed Description

Definition at line 15 of file PerJetFlavourUncertaintyComponent.h.

Constructor & Destructor Documentation

◆ PerJetFlavourUncertaintyComponent() [1/3]

jet::PerJetFlavourUncertaintyComponent::PerJetFlavourUncertaintyComponent ( const ComponentHelper & component )

Definition at line 29 of file PerJetFlavourUncertaintyComponent.cxx.

     : UncertaintyComponent(component,1)
     , m_absEta(CompParametrization::isAbsEta(component.parametrization))
     , m_labels(component.truthLabels)
     , m_flavourType(component.flavourType)
     , m_constrainZresponse(component.constrainZresponse)
     , m_constrainZresponseFunc(component.constrainZresponseFunc)
 {
     ATH_MSG_DEBUG("Created PerJetFlavourUncertaintyComponent named " << m_uncHistName.Data());
 }

◆ PerJetFlavourUncertaintyComponent() [2/3]

jet::PerJetFlavourUncertaintyComponent::PerJetFlavourUncertaintyComponent ( const PerJetFlavourUncertaintyComponent & toCopy )

Definition at line 40 of file PerJetFlavourUncertaintyComponent.cxx.

     : UncertaintyComponent(toCopy)
     , m_absEta(toCopy.m_absEta)
     , m_labels(toCopy.m_labels)
     , m_flavourType(toCopy.m_flavourType)
     , m_constrainZresponse(toCopy.m_constrainZresponse)
     , m_constrainZresponseFunc(toCopy.m_constrainZresponseFunc)
 {
     ATH_MSG_DEBUG("Creating copy of PerJetFlavourUncertaintyComponent named " << m_uncHistName.Data());
 }

◆ ~PerJetFlavourUncertaintyComponent()

virtual jet::PerJetFlavourUncertaintyComponent::~PerJetFlavourUncertaintyComponent ( )

inlinevirtual

Definition at line 22 of file PerJetFlavourUncertaintyComponent.h.

22 {}

◆ PerJetFlavourUncertaintyComponent() [3/3]

jet::PerJetFlavourUncertaintyComponent::PerJetFlavourUncertaintyComponent ( const std::string & name = "" )

private

Definition at line 17 of file PerJetFlavourUncertaintyComponent.cxx.

     : UncertaintyComponent(ComponentHelper(name),0)
     , m_absEta(false)
     , m_labels()
     , m_flavourType(FlavourComp::UNKNOWN)
     , m_constrainZresponse(false)
     , m_constrainZresponseFunc("Zjet_qFrac_AntiKt4EMPFlow")
     , m_ZjetQuarkFrac(nullptr)
 {
     JESUNC_NO_DEFAULT_CONSTRUCTOR;
 }

Member Function Documentation

◆ checkTruthLabel()

bool jet::PerJetFlavourUncertaintyComponent::checkTruthLabel ( const xAOD::Jet & jet ) const

private

Definition at line 173 of file PerJetFlavourUncertaintyComponent.cxx.

 {
     static const SG::AuxElement::ConstAccessor<int> truthLabelAccessor ("PartonTruthLabelID");
     if (!truthLabelAccessor.isAvailable(jet))
     {
         ATH_MSG_ERROR("Unable to find PartonTruthLabelID on the jet");
         return JESUNC_ERROR_CODE;
     }
     const int truthHighestEparton = abs(truthLabelAccessor(jet));
  
     for (const int label : m_labels)
         if (label == truthHighestEparton)
             return true;
     return false;
 }

◆ clone()

PerJetFlavourUncertaintyComponent * jet::PerJetFlavourUncertaintyComponent::clone ( ) const

virtual

Implements jet::UncertaintyComponent.

Definition at line 51 of file PerJetFlavourUncertaintyComponent.cxx.

 {
     return new PerJetFlavourUncertaintyComponent(*this);
 }

◆ getAbsMass()

double jet::UncertaintyComponent::getAbsMass	(	const xAOD::Jet &	jet,
		const CompMassDef::TypeEnum	massDef
	)		const

protectedvirtualinherited

Definition at line 309 of file UncertaintyComponent.cxx.

 {
     // Check for the simple case (where we want the default four-vector itself)
     if (massDef == CompMassDef::UNKNOWN || massDef == CompMassDef::FourVecMass)
         return jet.m();
     
     // Not the simple case, check for the specified four-vector instead and return it if it is available
     JetFourMomAccessor scale(CompMassDef::getJetScaleString(massDef).Data());
     if (scale.isAvailable(jet))
         return scale(jet).M();
  
     // Fall-back on the TA moment as a float if applicable (legacy support)
     SG::AuxElement::ConstAccessor<float> scaleTAMoment("JetTrackAssistedMassCalibrated");
     if (massDef == CompMassDef::TAMass && scaleTAMoment.isAvailable(jet))
         return scaleTAMoment(jet);
  
     // Fall-back on the calo mass as the 4-vec if applicable (legacy support)
     if (massDef == CompMassDef::CaloMass)
         return jet.m();
  
     // Specified scale is not available, error
     ATH_MSG_ERROR("Failed to retrieve the " << CompMassDef::enumToString(massDef).Data() << " mass from the jet");
     return JESUNC_ERROR_CODE;
 }

◆ getFlavourResponseUncertainty()

double jet::PerJetFlavourUncertaintyComponent::getFlavourResponseUncertainty	(	const xAOD::Jet &	jet,
		const xAOD::EventInfo &	eInfo
	)		const

private

Definition at line 117 of file PerJetFlavourUncertaintyComponent.cxx.

 {
     // sigma_x = l_x * dR_x(x response modelling uncertainty)
     //      sigma_x is the flavour response uncertainty for jet type x
     //      l_x is the jet flavour label, and is 0 or 1
     //      dR_x is the response modelling uncertainty for jets of type x
     //      dR_x is typically taken as the difference in x response between Pythia and Herwig
     
     // Check the jet label (is this uncertainty relevant?)
     if (!checkTruthLabel(jet))
         return 0;
  
     // This is a jet of the relevant type
     // Now calculate and return the uncertainty
     const double pT  = jet.pt()*m_energyScale;
     const double eta = m_absEta ? fabs(jet.eta()) : jet.eta();
  
     FlavourComp::TypeEnum ThisJetFlavourType = FlavourComp::UNKNOWN;
     if(m_labels.at(0)==21 || m_labels.at(0)==0)   ThisJetFlavourType = FlavourComp::PerJetResponse_Gluon;
     else if(m_labels.at(0)==1 || m_labels.at(0)==2 || m_labels.at(0)==3) ThisJetFlavourType = FlavourComp::PerJetResponse_LQ;
     else if(m_labels.at(0)==5)  ThisJetFlavourType = FlavourComp::PerJetResponse_B;
     else if(m_labels.at(0)==4 )  ThisJetFlavourType = FlavourComp::PerJetResponse_C;
     
     // bool DoesItPass = false;
     if(m_flavourType == ThisJetFlavourType){
         double unc = m_uncHist->getValue(pT,eta);
         if (m_constrainZresponse){
           if (m_flavourType==FlavourComp::PerJetResponse_Gluon) unc *= m_ZjetQuarkFrac->Eval(pT);
           if (m_flavourType==FlavourComp::PerJetResponse_LQ) unc *= (m_ZjetQuarkFrac->Eval(pT)-1.0);
         }
         // Return the uncertainty
         return unc;
  
     }else{ return 0;}
  
 }

◆ getMassOverE()

double jet::UncertaintyComponent::getMassOverE	(	const xAOD::Jet &	jet,
		const CompMassDef::TypeEnum	massDef
	)		const

protectedvirtualinherited

Definition at line 360 of file UncertaintyComponent.cxx.

 {
     // Check for the simple case (where we want the default four-vector itself)
     if (massDef == CompMassDef::UNKNOWN || massDef == CompMassDef::FourVecMass)
         return jet.m()/jet.e();
  
     // Not the simple case, check for the specified four-vector instead and return it if it is available
     JetFourMomAccessor scale(CompMassDef::getJetScaleString(massDef).Data());
     if (scale.isAvailable(jet))
         return scale(jet).M()/scale(jet).E();
     
     // Fall-back on the TA moment as a float if applicable (legacy support)
     SG::AuxElement::ConstAccessor<float> scaleTAMoment("JetTrackAssistedMassCalibrated");
     if (massDef == CompMassDef::TAMass && scaleTAMoment.isAvailable(jet))
         return scaleTAMoment(jet)/jet.e();
     
     // Fall-back on the calo mass as the 4-vec if applicable (legacy support)
     if (massDef == CompMassDef::CaloMass)
         return jet.m()/jet.e();
  
     // Specified scale is not available, error
     ATH_MSG_ERROR("Failed to retrieve the " << CompMassDef::enumToString(massDef).Data() << " mass from the jet");
     return JESUNC_ERROR_CODE;
  
 }

◆ getMassOverPt()

double jet::UncertaintyComponent::getMassOverPt	(	const xAOD::Jet &	jet,
		const CompMassDef::TypeEnum	massDef
	)		const

protectedvirtualinherited

Definition at line 334 of file UncertaintyComponent.cxx.

 {
     // Check for the simple case (where we want the default four-vector itself)
     if (massDef == CompMassDef::UNKNOWN || massDef == CompMassDef::FourVecMass)
         return jet.m()/jet.pt();
     
     // Not the simple case, check for the specified four-vector instead and return it if it is available
     JetFourMomAccessor scale(CompMassDef::getJetScaleString(massDef).Data());
     if (scale.isAvailable(jet))
         return scale(jet).M()/scale(jet).Pt();
     
     // Fall-back on the TA moment as a float if applicable (legacy support)
     SG::AuxElement::ConstAccessor<float> scaleTAMoment("JetTrackAssistedMassCalibrated");
     if (massDef == CompMassDef::TAMass && scaleTAMoment.isAvailable(jet))
         return scaleTAMoment(jet)/jet.pt();
     
     // Fall-back on the calo mass as the 4-vec if applicable (legacy support)
     if (massDef == CompMassDef::CaloMass)
         return jet.m()/jet.pt();
  
     // Specified scale is not available, error
     ATH_MSG_ERROR("Failed to retrieve the " << CompMassDef::enumToString(massDef).Data() << " mass from the jet");
     return JESUNC_ERROR_CODE;
  
 }

◆ getName()

virtual TString jet::UncertaintyComponent::getName ( ) const

inlinevirtualinherited

Definition at line 35 of file UncertaintyComponent.h.

35 { return m_uncHistName; }

◆ getRelevantLabels()

virtual std::vector<int> jet::PerJetFlavourUncertaintyComponent::getRelevantLabels ( ) const

inlinevirtual

Definition at line 26 of file PerJetFlavourUncertaintyComponent.h.

26 { return m_labels; }

◆ getScaleVar()

virtual CompScaleVar::TypeEnum jet::UncertaintyComponent::getScaleVar ( ) const

inlinevirtualinherited

Definition at line 37 of file UncertaintyComponent.h.

37 { return m_scaleVar; }

◆ getSplitFactor()

double jet::UncertaintyComponent::getSplitFactor ( const xAOD::Jet & jet ) const

protectedvirtualinherited

Definition at line 196 of file UncertaintyComponent.cxx.

 {
     // Immediate return for the most common case
     if (!m_splitNumber)
         return 1;
  
     // SplitNumber was specified, we have to calculate the factor
     double splitFactor = 1;
     const TH1* histo = m_uncHist->getHisto();
  
     if (m_splitNumber == 1 || m_splitNumber == -1)
     {
         // Straight line in log(pT) from 0 to 1
         // y = mx+b
         // m = 1/[log(max)-log(min)]
         // x = log(min) --> b = -m*log(min)
  
         const double minPt = histo->GetXaxis()->GetBinLowEdge(1);
         const double maxPt = histo->GetXaxis()->GetBinLowEdge(histo->GetNbinsX()+1);
         const double valPt = jet.pt()*m_energyScale;
  
         const double slope = 1./(log(maxPt)-log(minPt));
         const double intercept = -slope*log(minPt);
  
         splitFactor = slope*log(valPt <= minPt ? minPt+1.e-3 : valPt >= maxPt ? maxPt-1.e-3 : valPt)+intercept;
     }
     else if (m_splitNumber == 2 || m_splitNumber == -2)
     {
         // Straight line in |eta| from 0 to 1
         // y = mx + b
         // m = 1/(max - min)
         // x = min --> b = -m*min
         const double minEta = 0;
         const double maxEta = 4.5;
         const double absEta = fabs(jet.eta());
  
         const double slope = 1./(maxEta - minEta);
         const double intercept = -slope*minEta;
  
         splitFactor = slope*(absEta <= minEta ? minEta+1.e-3 : absEta >= maxEta ? maxEta-1.e-3 : absEta)+intercept;
     }
     else if (m_splitNumber == 3 || m_splitNumber == -3)
     {
         // Increasing with log(pT) and increasing with |eta|
         // z = mx + ny + c
         // z(min,min) = 0
         // z(max,max) = 1
         // Linear in both dimensions means need factor of 1/2 in a single dimension
         // m = 0.5/[log(maxPt)-log(minPt)]
         // n = 0.5/(maxEta - minEta)
         // c = -minPt*m - minEta*n
         // 2*z = (logx-logxmin)/(logxmax-logxmin) + (y-ymin)/(ymax-ymin)
         const double minPt = histo->GetXaxis()->GetBinLowEdge(1);
         const double maxPt = histo->GetXaxis()->GetBinLowEdge(histo->GetNbinsX()+1);
         const double valPt = jet.pt()*m_energyScale;
  
         const double minEta = 0;
         const double maxEta = 4.5;
         const double absEta = fabs(jet.eta());
  
         const double slopePt = 1./(log(maxPt)-log(minPt));
         const double slopeEta = 1./(maxEta-minEta);
         
         const double fixedPt  = valPt <= minPt ? minPt+1.e-3 : valPt >= maxPt ? maxPt - 1.e-3 : valPt;
         const double fixedEta = absEta <= minEta ? minEta+1.e-3 : absEta >= maxEta ? maxEta-1.e-3 : absEta;
  
         splitFactor = (slopePt*(log(fixedPt)-log(minPt)) + slopeEta*(fixedEta-minEta))/2.;
     }
     else if (m_splitNumber == 4 || m_splitNumber == -4)
     {
         // Increasing with log(pT) and decreasing with |eta|
         // See description above, follows similarly
         // 2*z = (logx-logxmin)/(logxmax-logxmin) + (ymax-y)/(ymax-ymin)
         const double minPt = histo->GetXaxis()->GetBinLowEdge(1);
         const double maxPt = histo->GetXaxis()->GetBinLowEdge(histo->GetNbinsX()+1);
         const double valPt = jet.pt()*m_energyScale;
  
         const double minEta = 0;
         const double maxEta = 4.5;
         const double absEta = fabs(jet.eta());
  
         const double slopePt = 1./(log(maxPt)-log(minPt));
         const double slopeEta = 1./(maxEta-minEta);
         
         const double fixedPt  = valPt <= minPt ? minPt+1.e-3 : valPt >= maxPt ? maxPt - 1.e-3 : valPt;
         const double fixedEta = absEta <= minEta ? minEta+1.e-3 : absEta >= maxEta ? maxEta-1.e-3 : absEta;
  
         splitFactor = (slopePt*(log(fixedPt)-log(minPt)) + slopeEta*(maxEta-fixedEta))/2.;
     }
  
  
     // Now check if this is the functional part or the complementary part
     if (m_splitNumber < 0)
         splitFactor = sqrt(1-splitFactor*splitFactor);
  
     return splitFactor;
 }

◆ getTopology()

virtual JetTopology::TypeEnum jet::UncertaintyComponent::getTopology ( ) const

inlinevirtualinherited

Definition at line 38 of file UncertaintyComponent.h.

38 { return m_topology; }

◆ getUncertainty()

double jet::UncertaintyComponent::getUncertainty	(	const xAOD::Jet &	jet,
		const xAOD::EventInfo &	eInfo
	)		const

virtualinherited

Definition at line 169 of file UncertaintyComponent.cxx.

 {
     if (!m_isInit)
     {
         ATH_MSG_ERROR("Component hasn't been initialized: " << getName().Data());
         return JESUNC_ERROR_CODE;
     }
     return getUncertaintyImpl(jet,eInfo)*getSplitFactor(jet);
 }

◆ getUncertaintyImpl()

double jet::PerJetFlavourUncertaintyComponent::getUncertaintyImpl	(	const xAOD::Jet &	jet,
		const xAOD::EventInfo &	eInfo
	)		const

protectedvirtual

Implements jet::UncertaintyComponent.

Definition at line 112 of file PerJetFlavourUncertaintyComponent.cxx.

 {
     return getFlavourResponseUncertainty(jet,eInfo);
 }

◆ getValidBool()

bool jet::UncertaintyComponent::getValidBool ( const double validity ) const

protectedvirtualinherited

Definition at line 301 of file UncertaintyComponent.cxx.

 {
     if (validity < 1.e-5 && validity > -1.e-5) return false;
     if (validity < 1+1.e-5 && validity > 1-1.e-5) return true;
     ATH_MSG_ERROR(Form("Validity value not in expected range: %lf for histogram %s",validity,getValidName().Data()));
     return false;
 }

◆ getValidity()

bool jet::UncertaintyComponent::getValidity	(	const xAOD::Jet &	jet,
		const xAOD::EventInfo &	eInfo
	)		const

virtualinherited

Definition at line 159 of file UncertaintyComponent.cxx.

 {
     if (!m_isInit)
     {
         ATH_MSG_ERROR("Component hasn't been initialized: " << getName().Data());
         return false;
     }
     return getValidityImpl(jet,eInfo);
 }

◆ getValidityImpl()

bool jet::PerJetFlavourUncertaintyComponent::getValidityImpl	(	const xAOD::Jet &	jet,
		const xAOD::EventInfo &	eInfo
	)		const

protectedvirtual

Implements jet::UncertaintyComponent.

Definition at line 107 of file PerJetFlavourUncertaintyComponent.cxx.

 {
     return !m_validHist ? true : getValidBool(m_validHist->getValue(jet.pt()*m_energyScale,m_absEta ? fabs(jet.eta()) : jet.eta()));
 }

◆ getValidName()

virtual TString jet::UncertaintyComponent::getValidName ( ) const

inlinevirtualinherited

Definition at line 36 of file UncertaintyComponent.h.

36 { return m_validHistName; }

◆ getValidUncertainty()

bool jet::UncertaintyComponent::getValidUncertainty	(	double &	unc,
		const xAOD::Jet &	jet,
		const xAOD::EventInfo &	eInfo
	)		const

virtualinherited

Definition at line 179 of file UncertaintyComponent.cxx.

 {
     if (getValidity(jet,eInfo))
     {
         unc = getUncertainty(jet,eInfo);
         return true;
     }
     return false;
 }

◆ initialize()

StatusCode jet::PerJetFlavourUncertaintyComponent::initialize ( TFile * histFile )

virtual

Reimplemented from jet::UncertaintyComponent.

Definition at line 56 of file PerJetFlavourUncertaintyComponent.cxx.

 {
     // Call the base class first
     if (UncertaintyComponent::initialize(histFile).isFailure())
         return StatusCode::FAILURE;
  
     // Ensure that the labels are sane
     if (m_labels.empty())
     {
         ATH_MSG_ERROR("Specified a PerJetFlavourUncertaintyComponent without the associated labels to use: " << getName().Data());
         return StatusCode::FAILURE;
     }
     for (int aLabel : m_labels)
     {
         if (!isSupportedLabel(aLabel))
         {
             ATH_MSG_ERROR("Unsupported label ID of " << aLabel << " in " << getName().Data());
             return StatusCode::FAILURE;
         }
     }
  
     // Load TF1 for Z-jet gluon fraction
     if (m_constrainZresponse)
     {
         // Find the histogram by name
         TObject* tmp_tf1 = histFile->Get(m_constrainZresponseFunc);
         if (!tmp_tf1)
         {
             ATH_MSG_ERROR(Form("Histogram file does not contain a histogram named %s",m_constrainZresponseFunc.Data()));
             return StatusCode::FAILURE;
         }
         // Ensure the object is a TF1
         m_ZjetQuarkFrac = dynamic_cast<const TF1*>(tmp_tf1);
         if (!m_ZjetQuarkFrac)
         {
             ATH_MSG_ERROR(Form("Histogram file contains the expected key, but it's not a TF1* (%s)",m_constrainZresponseFunc.Data()));
             return StatusCode::FAILURE;
         }
     }
  
     // Done checking the configuration
     return StatusCode::SUCCESS;
 }

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

◆ isAlwaysZero()

bool jet::UncertaintyComponent::isAlwaysZero ( ) const

virtualinherited

Reimplemented in jet::CombinedMassUncertaintyComponent.

Definition at line 140 of file UncertaintyComponent.cxx.

 {
     if (!m_isInit)
     {
         ATH_MSG_ERROR("Cannot call method before initialization, component: " << getName().Data());
         return false;
     }
  
     const TH1* histo = m_uncHist->getHisto();
     return !(fabs(histo->GetMinimum()) > 1.e-8 || fabs(histo->GetMaximum()) > 1.e-8);
 }

◆ isSupportedLabel()

bool jet::PerJetFlavourUncertaintyComponent::isSupportedLabel ( const int label ) const

private

Definition at line 156 of file PerJetFlavourUncertaintyComponent.cxx.

 {
     switch (abs(label))
     {
         case 1:  // parton-label, up
         case 2:  // parton-label, down
         case 3:  // parton-label, strange
         case 4:  // parton-label, charm
         case 5:  // parton-label, bottom
         case 21: // parton-label, gluon
         case 0: // parton-label, PileUp/non-identified
             return true;
         default:
             return false;
     }
 }

◆ 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
    }

◆ 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_absEta

const bool jet::PerJetFlavourUncertaintyComponent::m_absEta {}

private

Definition at line 39 of file PerJetFlavourUncertaintyComponent.h.

◆ m_constrainZresponse

const bool jet::PerJetFlavourUncertaintyComponent::m_constrainZresponse {}

private

Definition at line 42 of file PerJetFlavourUncertaintyComponent.h.

◆ m_constrainZresponseFunc

const TString jet::PerJetFlavourUncertaintyComponent::m_constrainZresponseFunc {}

private

Definition at line 43 of file PerJetFlavourUncertaintyComponent.h.

◆ m_energyScale

const float jet::UncertaintyComponent::m_energyScale {}

protectedinherited

Definition at line 55 of file UncertaintyComponent.h.

◆ m_flavourType

const FlavourComp::TypeEnum jet::PerJetFlavourUncertaintyComponent::m_flavourType {}

private

Definition at line 41 of file PerJetFlavourUncertaintyComponent.h.

◆ m_imsg

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

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

◆ m_interpolate

const Interpolate::TypeEnum jet::UncertaintyComponent::m_interpolate

protectedinherited

Definition at line 56 of file UncertaintyComponent.h.

◆ m_isInit

bool jet::UncertaintyComponent::m_isInit {}

protectedinherited

Definition at line 50 of file UncertaintyComponent.h.

◆ m_labels

const std::vector<int> jet::PerJetFlavourUncertaintyComponent::m_labels

private

Definition at line 40 of file PerJetFlavourUncertaintyComponent.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_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_numExpectedHist

int jet::UncertaintyComponent::m_numExpectedHist {}

protectedinherited

Definition at line 59 of file UncertaintyComponent.h.

◆ m_scaleVar

const CompScaleVar::TypeEnum jet::UncertaintyComponent::m_scaleVar

protectedinherited

Definition at line 53 of file UncertaintyComponent.h.

◆ m_splitNumber

const int jet::UncertaintyComponent::m_splitNumber {}

protectedinherited

Definition at line 57 of file UncertaintyComponent.h.

◆ m_topology

const JetTopology::TypeEnum jet::UncertaintyComponent::m_topology

protectedinherited

Definition at line 54 of file UncertaintyComponent.h.

◆ m_uncHist

UncertaintyHistogram* jet::UncertaintyComponent::m_uncHist {}

protectedinherited

Definition at line 60 of file UncertaintyComponent.h.

◆ m_uncHistName

const TString jet::UncertaintyComponent::m_uncHistName

protectedinherited

Definition at line 51 of file UncertaintyComponent.h.

◆ m_validHist

UncertaintyHistogram* jet::UncertaintyComponent::m_validHist {}

protectedinherited

Definition at line 61 of file UncertaintyComponent.h.

◆ m_validHistName

const TString jet::UncertaintyComponent::m_validHistName

protectedinherited

Definition at line 52 of file UncertaintyComponent.h.

◆ m_ZjetQuarkFrac

const TF1* jet::PerJetFlavourUncertaintyComponent::m_ZjetQuarkFrac {}

private

Definition at line 44 of file PerJetFlavourUncertaintyComponent.h.

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

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ PerJetFlavourUncertaintyComponent() [1/3]

◆ PerJetFlavourUncertaintyComponent() [2/3]

◆ ~PerJetFlavourUncertaintyComponent()

◆ PerJetFlavourUncertaintyComponent() [3/3]

Member Function Documentation

◆ checkTruthLabel()

◆ clone()

◆ getAbsMass()

◆ getFlavourResponseUncertainty()

◆ getMassOverE()

◆ getMassOverPt()

◆ getName()

◆ getRelevantLabels()

◆ getScaleVar()

◆ getSplitFactor()

◆ getTopology()

◆ getUncertainty()

◆ getUncertaintyImpl()

◆ getValidBool()

◆ getValidity()

◆ getValidityImpl()

◆ getValidName()

◆ getValidUncertainty()

◆ initialize()

◆ initMessaging()

◆ isAlwaysZero()

◆ isSupportedLabel()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ setLevel()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_absEta

◆ m_constrainZresponse

◆ m_constrainZresponseFunc

◆ m_energyScale

◆ m_flavourType

◆ m_imsg

◆ m_interpolate

◆ m_isInit

◆ m_labels

◆ m_lvl

◆ m_msg_tls

◆ m_nm

◆ m_numExpectedHist

◆ m_scaleVar

◆ m_splitNumber

◆ m_topology

◆ m_uncHist

◆ m_uncHistName

◆ m_validHist

◆ m_validHistName

◆ m_ZjetQuarkFrac