jet::PerJetFlavourUncertaintyComponent Class Reference

#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.
Functions providing the same interface as AthMessaging
bool	msgLvl (const MSG::Level lvl) const
	Test the output level of the object.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.

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.

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.
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels)
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level.
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging)

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.

{}

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.

{ return m_uncHistName;   }

getRelevantLabels()

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

inlinevirtual

Definition at line 26 of file PerJetFlavourUncertaintyComponent.h.

{ return m_labels; }

getScaleVar()

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

inlinevirtualinherited

Definition at line 37 of file UncertaintyComponent.h.

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

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

{ 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();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    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.

{ nullptr };

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.

{ MSG::NIL };

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:

• PerJetFlavourUncertaintyComponent.h
• PerJetFlavourUncertaintyComponent.cxx