JetSubStructureUtils::Pull Class Reference

#include <Pull.h>

Inheritance diagram for JetSubStructureUtils::Pull:

Collaboration diagram for JetSubStructureUtils::Pull:

Public Member Functions
	Pull (bool useEtaInsteadOfY=false)
virtual std::map< std::string, double >	result (const fastjet::PseudoJet &jet) const
virtual TOut	result (const xAOD::Jet &jet) const

Private Attributes
bool	m_useEtaInsteadOfY

Detailed Description

Definition at line 11 of file Pull.h.

Constructor & Destructor Documentation

Pull()

JetSubStructureUtils::Pull::Pull ( bool useEtaInsteadOfY = false )

inline

Definition at line 13 of file Pull.h.

: m_useEtaInsteadOfY(useEtaInsteadOfY) {};

Member Function Documentation

result() [1/2]

map< string, double > Pull::result ( const fastjet::PseudoJet & jet ) const

virtual

Definition at line 12 of file Pull.cxx.

{
  vector<fastjet::PseudoJet> constit_pseudojets = jet.constituents();
  map<string, double> Variables;
  Variables["PullMag"] =  -999.*1000.;
  Variables["PullPhi"] =  -999.*1000.;
  Variables["Pull_C00"] =  -999.*1000.;
  Variables["Pull_C01"] =  -999.*1000.;
  Variables["Pull_C10"] =  -999.*1000.;
  Variables["Pull_C11"] =  -999.*1000.;
 
  if (constit_pseudojets.size() < 2) {
    return Variables;
  }
 
  double T[6];
  for (int i=0; i<6; i++){
    T[i] = 0;
  }
 
  TVector2 r_pull;//pull (times jet pt)
  double jeta = m_useEtaInsteadOfY ? jet.eta() : jet.rapidity();
 
  for(vector<fastjet::PseudoJet>::iterator cit = constit_pseudojets.begin(); cit != constit_pseudojets.end(); ++cit) {
    const fastjet::PseudoJet &constit = *cit;
 
    double dphi = jet.delta_phi_to(constit); 
    double ceta = m_useEtaInsteadOfY ? constit.eta() : constit.rapidity();
 
    TVector2 r_i(ceta-jeta, dphi);  //position on eta-phi plane relative to weighted centroid
    T[2] += (ceta*ceta) - (jeta*jeta);
    T[3] += r_i.X() * r_i.Y();
    T[4] += r_i.Y() * r_i.X();
    T[5] += (constit.phi_std() * constit.phi_std()) - (jet.phi_std() * jet.phi_std());
    r_i *= (constit.pt() * r_i.Mod());
    r_pull += r_i;
  }
 
  if(jet.pt() != 0) {
    T[0] = r_pull.Mod()/jet.pt();
  }
 
  T[1] = TVector2::Phi_mpi_pi(r_pull.Phi()); //this means something quite different than statshape package at ANL!!
 
  Variables["PullMag"]  = T[0];
  Variables["PullPhi"]  = T[1];
  Variables["Pull_C00"] = T[2];
  Variables["Pull_C01"] = T[3];
  Variables["Pull_C10"] = T[4];
  Variables["Pull_C11"] = T[5];
 
  return Variables; 
}

result() [2/2]

virtual TOut JetSubStructureUtils::SubstructureCalculator< TOut >::result ( const xAOD::Jet & jet ) const

inlinevirtual

Reimplemented from JetSubStructureUtils::SubstructureCalculator< std::map< std::string, double > >.

Definition at line 25 of file SubstructureCalculator.h.

                                                    {
        // PS 4/18 master developent
        // std::vector<fastjet::PseudoJet> constit_pseudojets =
        //   jet::JetConstituentFiller::constituentPseudoJets(jet);
 
        std::vector<fastjet::PseudoJet> constit_pseudojets;
        std::transform(jet.getConstituents().begin(),
                       jet.getConstituents().end(),
                       std::back_inserter(constit_pseudojets),
                       [](const auto& c){
                         const xAOD::IParticle* ip = c->rawConstituent();
                         return
                           // fastjet::PseudoJet((c->rawConstituent())->p4());
                           fastjet::PseudoJet(ip->p4());
                       });
 
        fastjet::PseudoJet pjet = fastjet::join(constit_pseudojets);
 
        return result(pjet);
      }

Member Data Documentation

m_useEtaInsteadOfY

bool JetSubStructureUtils::Pull::m_useEtaInsteadOfY

private

Definition at line 19 of file Pull.h.

---

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

• Pull.h
• Pull.cxx