#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

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

Detailed Description

Definition at line 11 of file Angularity.h.

Member Function Documentation

◆ result() [1/3]

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/3]

virtual TOut JetSubStructureUtils::SubstructureCalculator< TOut >::result

inline

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

◆ result() [3/3]

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

inlinevirtualinherited