JetSubStructureUtils::Angularity Class Reference

#include <Angularity.h>

Inheritance diagram for JetSubStructureUtils::Angularity:

Collaboration diagram for JetSubStructureUtils::Angularity:

Public Member Functions
virtual double	result (const fastjet::PseudoJet &jet) const
virtual TOut	result (const xAOD::Jet &jet) const

Detailed Description

Definition at line 11 of file Angularity.h.

Member Function Documentation

result() [1/2]

double Angularity::result ( const fastjet::PseudoJet & jet ) const

virtual

Definition at line 11 of file Angularity.cxx.

{
  if(jet.constituents().empty()) return -1;
  if(jet.m()==0.0) return -1;
 
  vector<fastjet::PseudoJet> constit_pseudojets = jet.constituents();
  TLorentzVector jet_p4(jet.px(), jet.py(), jet.pz(), jet.e());
 
  double Angularity2=-1.;
  const double a2=-2.;
  double sum_a2=0.;
 
  for(unsigned int iConstit=0; iConstit < constit_pseudojets.size(); iConstit++) {
    TLorentzVector tclus = TLorentzVector(constit_pseudojets[iConstit].px(),constit_pseudojets[iConstit].py(),constit_pseudojets[iConstit].pz(),constit_pseudojets[iConstit].e());
    double theta_i = jet_p4.Angle(tclus.Vect());
    double sintheta_i = sin(theta_i);
    if( sintheta_i == 0 ) continue; // avoid FPE
    double e_theta_i_a2 = constit_pseudojets[iConstit].E()*pow(sintheta_i,a2)*pow(1-cos(theta_i),1-a2);
    sum_a2 += e_theta_i_a2;
  }
 
  if ( jet.m() < 1.e-20 ) return -1.0;
  Angularity2 = sum_a2/jet.m(); 
  return Angularity2;
}

result() [2/2]

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

inlinevirtual

Reimplemented from JetSubStructureUtils::SubstructureCalculator< 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);
      }

---

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

• Angularity.h
• Angularity.cxx