#include <Puppi.h>

Collaboration diagram for Puppi:

Public Member Functions
	Puppi (double R0, double Rmin, double beta, double centralPTCutOffset, double centralPTCutSlope, double forwardPTCutOffset, double forwardPTCutSlope, double etaBoundary)
void	setParticles (const std::vector< fastjet::PseudoJet > &chargedHS, const std::vector< fastjet::PseudoJet > &chargedPU, const std::vector< fastjet::PseudoJet > &neutral, const std::vector< fastjet::PseudoJet > &forward, int nPU)
double	getWeight (const fastjet::PseudoJet &pfo)
double	getAlpha (const fastjet::PseudoJet &pfo)
double	getMedian ()
double	getRMS ()

Private Member Functions
double	getChi2 (const fastjet::PseudoJet &pfo)
void	findAlphaMedianAndRMS ()

Private Attributes
std::vector< fastjet::PseudoJet >	m_chargedHS
std::vector< fastjet::PseudoJet >	m_chargedPU
std::vector< fastjet::PseudoJet >	m_neutral
std::vector< fastjet::PseudoJet >	m_forward
std::vector< std::vector< fastjet::PseudoJet > * >	m_allParticles
double	m_R0 {}
double	m_Rmin {}
double	m_beta {}
double	m_centralPTCutOffset {}
double	m_centralPTCutSlope {}
double	m_forwardPTCutOffset {}
double	m_forwardPTCutSlope {}
double	m_etaBoundary {}
bool	m_includeCentralNeutralsInAlpha {}
double	m_median {}
double	m_rms {}
double	m_nPV {}

Detailed Description

Definition at line 22 of file Puppi.h.

Constructor & Destructor Documentation

◆ Puppi()

Puppi::Puppi	(	double	R0,
		double	Rmin,
		double	beta,
		double	centralPTCutOffset,
		double	centralPTCutSlope,
		double	forwardPTCutOffset,
		double	forwardPTCutSlope,
		double	etaBoundary )

Definition at line 19 of file Puppi.cxx.

                                                                                                                                                                             : 
  m_R0(R0), m_Rmin(Rmin), m_beta(beta), m_centralPTCutOffset(centralPTCutOffset), m_centralPTCutSlope(centralPTCutSlope), m_forwardPTCutOffset(forwardPTCutOffset), m_forwardPTCutSlope(forwardPTCutSlope), m_etaBoundary(etaBoundary)
{}

Member Function Documentation

◆ findAlphaMedianAndRMS()

void Puppi::findAlphaMedianAndRMS ( )

private

Definition at line 133 of file Puppi.cxx.

                                 {
  vector<double> values;
 
  for(const auto& p: m_chargedPU){
 
    // Don't want to include particles on the boundary
    if( fabs(p.eta()) > m_etaBoundary-m_R0) continue;
 
    double value = getAlpha(p);
    if(value > -999) values.push_back(value);  //-9999 is the value assigned to pfo with zero neighbors
  }
 
  std::sort(values.begin(),values.end());
            
  if(!values.empty()) m_median=values[int(values.size()*0.5)];
  else m_median=-9999;
 
  // now compute the LHS RMS
  double sum=0;
  int n = 0;
  for(auto value: values){
    if (value - m_median > 0) continue;
    sum += (value - m_median)*(value - m_median);
    ++n;
  }
 
  if(n > 0) m_rms = TMath::Sqrt(sum/n);
  else m_rms=-9999;
}

◆ getAlpha()

double Puppi::getAlpha ( const fastjet::PseudoJet & pfo )

Definition at line 83 of file Puppi.cxx.

                                                 {
  fastjet::Selector sel = fastjet::SelectorCircle(m_R0);
  sel.set_reference(pfo);
 
  double sum=0;
  int nNeighbors=0;
 
  if (fabs(pfo.eta())<m_etaBoundary+m_R0){
    vector<fastjet::PseudoJet> chargedHSNeighbors = sel(m_chargedHS);
    for (const auto& p: chargedHSNeighbors){
      float dR=pfo.delta_R(p);
      if (dR>m_Rmin){
    sum+=p.pt()/pow(dR, m_beta);
    nNeighbors+=1;
      }
    }
  }
 
  if (m_includeCentralNeutralsInAlpha){
    if (fabs(pfo.eta())<m_etaBoundary+m_R0){
      vector<fastjet::PseudoJet> neutralNeighbors = sel(m_neutral);
      for (const auto& p: neutralNeighbors){
    float dR=pfo.delta_R(p);
    if (dR>m_Rmin){
      sum+=pow(p.pt()/dR, m_beta);
      nNeighbors+=1;
    }
      }
    }
  }
 
  if (fabs(pfo.eta())>m_etaBoundary-m_R0){
    vector<fastjet::PseudoJet> forwardNeighbors = sel(m_forward);
    for (const auto& p: forwardNeighbors){
      float dR=pfo.delta_R(p);
      if (dR>m_Rmin){
    sum+=pow(p.pt()/dR, m_beta);
    nNeighbors+=1;
      }
    }
  }
 
  if (sum<=FLT_MIN) return -99999;
  if (nNeighbors!=0) return log(sum);
  return -9999;
}

◆ getChi2()

double Puppi::getChi2 ( const fastjet::PseudoJet & pfo )

private

Definition at line 42 of file Puppi.cxx.

                                                {
  double chi=(getAlpha(pfo)-m_median)/m_rms;
  return chi*fabs(chi);
}

◆ getMedian()

double Puppi::getMedian ( )

Definition at line 163 of file Puppi.cxx.

                       {
  return m_median;
}

◆ getRMS()

double Puppi::getRMS ( )

Definition at line 166 of file Puppi.cxx.

                    {
  return m_rms;
}

◆ getWeight()

double Puppi::getWeight ( const fastjet::PseudoJet & pfo )

Definition at line 49 of file Puppi.cxx.

                                                  {
 
  double chi2Total=getChi2(pfo);
  int nDF=1;
 
  if(pfo.has_user_info<PuppiUserInfo>()){
    //This is untested!!!!
    //Need to think about how to handle sign of chi2 - for now just throw it away
    chi2Total=fabs(chi2Total);
    for(auto v: pfo.user_info<PuppiUserInfo>().otherChi2Vec) {
      chi2Total+=v;
      nDF+=1;
    }
  }
 
  double w = ROOT::Math::chisquared_cdf(chi2Total, nDF);
 
  double offset,slope;
  if( fabs(pfo.eta()) < m_etaBoundary) {
    offset=m_centralPTCutOffset;
    slope=m_centralPTCutSlope;
  }
  else {
    offset=m_forwardPTCutOffset;
    slope=m_forwardPTCutSlope;
  }
 
  w*=(pfo.pt()>offset+m_nPV*slope);
 
  return w;
}

◆ setParticles()

void Puppi::setParticles	(	const std::vector< fastjet::PseudoJet > &	chargedHS,
		const std::vector< fastjet::PseudoJet > &	chargedPU,
		const std::vector< fastjet::PseudoJet > &	neutral,
		const std::vector< fastjet::PseudoJet > &	forward,
		int	nPU )

Definition at line 26 of file Puppi.cxx.

                                                                                                                                                                                                                   {
 
  m_chargedHS=chargedHS;
  m_chargedPU=chargedPU;
  m_neutral=neutral;
  m_forward=forward;
    
  m_nPV=nPV;
 
  m_allParticles = {&m_chargedHS, &m_chargedPU, &m_neutral, &m_forward};
    
  findAlphaMedianAndRMS();
}

Member Data Documentation

◆ m_allParticles

std::vector< std::vector< fastjet::PseudoJet >* > Puppi::m_allParticles

private

Definition at line 49 of file Puppi.h.

◆ m_beta

double Puppi::m_beta {}

private

Definition at line 53 of file Puppi.h.

53{};

◆ m_centralPTCutOffset

double Puppi::m_centralPTCutOffset {}

private

Definition at line 54 of file Puppi.h.

54{};

◆ m_centralPTCutSlope

double Puppi::m_centralPTCutSlope {}

private

Definition at line 55 of file Puppi.h.

55{};

◆ m_chargedHS

std::vector<fastjet::PseudoJet> Puppi::m_chargedHS

private

Definition at line 44 of file Puppi.h.

◆ m_chargedPU

std::vector<fastjet::PseudoJet> Puppi::m_chargedPU

private

Definition at line 45 of file Puppi.h.

◆ m_etaBoundary

double Puppi::m_etaBoundary {}

private

Definition at line 58 of file Puppi.h.

58{};

◆ m_forward

std::vector<fastjet::PseudoJet> Puppi::m_forward

private

Definition at line 47 of file Puppi.h.

◆ m_forwardPTCutOffset

double Puppi::m_forwardPTCutOffset {}

private

Definition at line 56 of file Puppi.h.

56{};

◆ m_forwardPTCutSlope

double Puppi::m_forwardPTCutSlope {}

private

Definition at line 57 of file Puppi.h.

57{};

◆ m_includeCentralNeutralsInAlpha

bool Puppi::m_includeCentralNeutralsInAlpha {}

private

Definition at line 59 of file Puppi.h.

59{};

◆ m_median

double Puppi::m_median {}

private

Definition at line 61 of file Puppi.h.

61{};

◆ m_neutral

std::vector<fastjet::PseudoJet> Puppi::m_neutral

private

Definition at line 46 of file Puppi.h.

◆ m_nPV

double Puppi::m_nPV {}

private

Definition at line 64 of file Puppi.h.

64{};

◆ m_R0

double Puppi::m_R0 {}

private

Definition at line 51 of file Puppi.h.

51{};

◆ m_Rmin

double Puppi::m_Rmin {}

private

Definition at line 52 of file Puppi.h.

52{};

◆ m_rms

double Puppi::m_rms {}

private

Definition at line 62 of file Puppi.h.

62{};

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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Puppi()

Member Function Documentation

◆ findAlphaMedianAndRMS()

◆ getAlpha()

◆ getChi2()

◆ getMedian()

◆ getRMS()

◆ getWeight()

◆ setParticles()

Member Data Documentation

◆ m_allParticles

◆ m_beta

◆ m_centralPTCutOffset

◆ m_centralPTCutSlope

◆ m_chargedHS

◆ m_chargedPU

◆ m_etaBoundary

◆ m_forward

◆ m_forwardPTCutOffset

◆ m_forwardPTCutSlope

◆ m_includeCentralNeutralsInAlpha

◆ m_median

◆ m_neutral

◆ m_nPV

◆ m_R0

◆ m_Rmin

◆ m_rms