JetSoftDrop.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
// JetSoftDrop.h
 
#ifndef JetSoftDrop_H
#define JetSoftDrop_H
 
// Joe Taenzer
// October 2016
//
// Tool to groom jets with softdrop.
// SoftDrop is described in this paper: arXiv:1402.2657
// See also: http://fastjet.hepforge.org/svn/contrib/contribs/RecursiveTools/tags/1.0.0/SoftDrop.hh
 
#include "AsgTools/AsgTool.h"
#include "JetInterface/IJetGroomer.h"
#include "JetInterface/IJetFromPseudojet.h"
#include "AsgTools/ToolHandle.h"
 
#include "fastjet/ClusterSequence.hh"
#include "fastjet/contrib/RecursiveSymmetryCutBase.hh"
#include "fastjet/contrib/SoftDrop.hh"
 
class JetSoftDrop
: public asg::AsgTool,
  virtual public IJetGroomer {
ASG_TOOL_CLASS(JetSoftDrop, IJetGroomer)
 
public:
 
  // Ctor.
  JetSoftDrop(const std::string& name);
 
  // Dtor.
  ~JetSoftDrop();
 
  // Initilization.
  StatusCode initialize();
 
  // Groom a jet and add result to a container.
  int groom(const xAOD::Jet& jin,
            const PseudoJetContainer&,
            xAOD::JetContainer& jout) const;
 
  // Dump to log.
  void print() const;
 
private:  // data
 
  // Job options.
  // SoftDrop algorithm:
  // z > zcut * (dR12/R0)^beta
  // z = min(pT1, pT2)/(pT1+pT2)
  // R0 = characteristic jet radius
  float m_zcut;                        // pT fraction for retaining subjets
  float m_beta;                        // How much to consider angular dependence
  float m_R0;                            // Normalization of angular distance, usually the characteristic jet radius (default R0 = 1)
  ToolHandle<IJetFromPseudojet> m_bld;  // Tool to build jets.
 
};
 
#endif