JetSubStructureUtils::Dipolarity Class Reference

#include <Dipolarity.h>

Inheritance diagram for JetSubStructureUtils::Dipolarity:

Collaboration diagram for JetSubStructureUtils::Dipolarity:

Public Member Functions
	Dipolarity (unsigned int n1, unsigned int n2, bool exclusive, double subjetR)
virtual double	result (const fastjet::PseudoJet &jet) const
virtual TOut	result (const xAOD::Jet &jet) const

Private Member Functions
double	dipolarity (std::vector< fastjet::PseudoJet > &constit_pseudojets, const fastjet::PseudoJet &jet1, const fastjet::PseudoJet &jet2) const

Private Attributes
bool	m_exclusive
unsigned int	m_n1
unsigned int	m_n2
double	m_subjetR

Detailed Description

Definition at line 12 of file Dipolarity.h.

Constructor & Destructor Documentation

Dipolarity()

Dipolarity::Dipolarity	(	unsigned int	n1,
		unsigned int	n2,
		bool	exclusive,
		double	subjetR )

Definition at line 15 of file Dipolarity.cxx.

                                                                                       :
  m_exclusive(exclusive),
  m_n1(n1),
  m_n2(n2),
  m_subjetR(subjetR) {
 
    if(m_n2 <= m_n1 || m_n1 == 0) {
    cout << "dipolarity was called with invalid arguments" << endl;
    }
}

Member Function Documentation

dipolarity()

double Dipolarity::dipolarity ( std::vector< fastjet::PseudoJet > & constit_pseudojets, const fastjet::PseudoJet & jet1, const fastjet::PseudoJet & jet2 ) const

private

Definition at line 62 of file Dipolarity.cxx.

{
  float dipolarity = 0;
  float sumpt = 0;
 
  TVector2 v12(jet2.eta() - jet1.eta(), jet1.delta_phi_to(jet2));
  if(v12.Mod2() < 0.001) return -1;
 
  for(unsigned int iConstit = 0; iConstit < constit_pseudojets.size(); iConstit++) {
    fastjet::PseudoJet constituent = constit_pseudojets[iConstit];
 
    sumpt += constituent.perp();
 
    TVector2 v;
 
    v.Set(constituent.eta() - jet1.eta(), jet1.delta_phi_to(constituent));
    float test = v * v12.Unit();
 
    if(test < 0)
      dipolarity += constituent.perp()*v.Mod2(); 
    else if(test < v12.Mod())
      dipolarity += constituent.perp()*(v.Mod2() - pow(v * v12.Unit(), 2));                                              
    else {
      v.Set(constituent.eta() - jet2.eta(), jet2.delta_phi_to(constituent));
      dipolarity += constituent.perp()*v.Mod2();                                                                           
    }
  }
 
  if(sumpt < 0.001) return -1;
  return dipolarity/(sumpt*v12.Mod2());
}

result() [1/2]

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

virtual

Definition at line 26 of file Dipolarity.cxx.

{
  vector<fastjet::PseudoJet> constit_pseudojets = jet.constituents();
  if(constit_pseudojets.empty()) return -1;
 
  vector<fastjet::PseudoJet> kt_subjets;
  fastjet::PseudoJet jet1, jet2;
 
  if(!m_exclusive) {
    fastjet::JetDefinition jet_def = fastjet::JetDefinition(fastjet::kt_algorithm, m_subjetR,
                                                            fastjet::E_scheme, fastjet::Best);
    fastjet::ClusterSequence kt_clust_seq(constit_pseudojets, jet_def);
    kt_subjets = fastjet::sorted_by_pt(kt_clust_seq.inclusive_jets(5000.0));
 
    if(kt_subjets.size() < m_n2) {
      //Not enough subjets
      return 0.0;
    }
  }
  else {
    if(m_n2 > constit_pseudojets.size()) {
      //We were asked to calculate exclusiveDipolarity, but there are not enough constituents
      return 0.0;
    }
 
    fastjet::JetDefinition jet_def = fastjet::JetDefinition(fastjet::kt_algorithm, 1.5,
                                                            fastjet::E_scheme, fastjet::Best);
    fastjet::ClusterSequence kt_clust_seq(constit_pseudojets, jet_def);
    kt_subjets = fastjet::sorted_by_pt(kt_clust_seq.exclusive_jets((int)m_n2));
  }
 
  jet1 = kt_subjets[m_n1-1];
  jet2 = kt_subjets[m_n2-1];
  return dipolarity(constit_pseudojets, jet1, jet2);
}

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

Member Data Documentation

m_exclusive

bool JetSubStructureUtils::Dipolarity::m_exclusive

private

Definition at line 23 of file Dipolarity.h.

m_n1

unsigned int JetSubStructureUtils::Dipolarity::m_n1

private

Definition at line 24 of file Dipolarity.h.

m_n2

unsigned int JetSubStructureUtils::Dipolarity::m_n2

private

Definition at line 25 of file Dipolarity.h.

m_subjetR

double JetSubStructureUtils::Dipolarity::m_subjetR

private

Definition at line 26 of file Dipolarity.h.

---

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

• Dipolarity.h
• Dipolarity.cxx