#include <CMSBoostedTauSeedingAlgorithm.h>

Inheritance diagram for contrib::CMSBoostedTauSeedingAlgorithm:

Public Types
typedef CMSBoostedTauSeedingAlgorithmStructure	StructureType

Public Member Functions
	CMSBoostedTauSeedingAlgorithm (double ptMin, double muMin, double muMax, double yMin, double yMax, double dRMin, double dRMax, int maxDepth, int verbosity=0)

std::string	description () const override

PseudoJet	result (const PseudoJet &jet) const override

	~CMSBoostedTauSeedingAlgorithm () override

Protected Member Functions
void	dumpSubJetStructure (const fastjet::PseudoJet &jet, int depth, int maxDepth, const std::string &depth_and_idx_string) const

std::pair< PseudoJet, PseudoJet >	findSubjets (const PseudoJet &jet, int depth, bool &subjetsFound) const

Private Attributes
double	dRMax_
	the max value of the dR parameter More...

double	dRMin_
	the min value of the dR parameter More...

int	maxDepth_
	the max depth for descending into clustering sequence More...

double	muMax_
	the max value of the mass-drop parameter More...

double	muMin_
	the min value of the mass-drop parameter More...

double	ptMin_
	minimum sub-jet pt More...

int	verbosity_
	flag to enable/disable debug output More...

double	yMax_
	the max value of the asymmetry parameter More...

double	yMin_
	the min value of the asymmetry parameter More...

Detailed Description

This class implements the CMS boosted tau algorithm

Definition at line 52 of file CMSBoostedTauSeedingAlgorithm.h.

Member Typedef Documentation

◆ StructureType

typedef CMSBoostedTauSeedingAlgorithmStructure contrib::CMSBoostedTauSeedingAlgorithm::StructureType

Definition at line 74 of file CMSBoostedTauSeedingAlgorithm.h.

Constructor & Destructor Documentation

◆ CMSBoostedTauSeedingAlgorithm()

contrib::CMSBoostedTauSeedingAlgorithm::CMSBoostedTauSeedingAlgorithm	(	double	ptMin,
		double	muMin,
		double	muMax,
		double	yMin,
		double	yMax,
		double	dRMin,
		double	dRMax,
		int	maxDepth,
		int	verbosity = `0`
	)

Definition at line 30 of file CMSBoostedTauSeedingAlgorithm.cc.

       : ptMin_(iminPt),
         muMin_(iminMassDrop),
         muMax_(imaxMassDrop),
         yMin_(iminY),
         yMax_(imaxY),
         dRMin_(iminDeltaR),
         dRMax_(imaxDeltaR),
         maxDepth_(maxDepth),
         verbosity_(verbosity) {}

◆ ~CMSBoostedTauSeedingAlgorithm()

contrib::CMSBoostedTauSeedingAlgorithm::~CMSBoostedTauSeedingAlgorithm ( )

inlineoverride

Definition at line 66 of file CMSBoostedTauSeedingAlgorithm.h.

66 {}

Member Function Documentation

◆ description()

std::string contrib::CMSBoostedTauSeedingAlgorithm::description ( ) const

override

Definition at line 51 of file CMSBoostedTauSeedingAlgorithm.cc.

                                                              {
     std::ostringstream oss;
     oss << "CMSBoostedTauSeedingAlgorithm algorithm";
     return oss.str();
   }

◆ dumpSubJetStructure()

void contrib::CMSBoostedTauSeedingAlgorithm::dumpSubJetStructure	(	const fastjet::PseudoJet &	jet,
		int	depth,
		int	maxDepth,
		const std::string &	depth_and_idx_string
	)		const

protected

Definition at line 58 of file CMSBoostedTauSeedingAlgorithm.cc.

                                                                                                          {
     if (maxDepth != -1 && depth > maxDepth)
       return;
     fastjet::PseudoJet subjet1, subjet2;
     bool hasSubjets = jet.has_parents(subjet1, subjet2);
     if (!hasSubjets)
       return;
     std::string depth_and_idx_string_subjet1 = depth_and_idx_string;
     if (depth_and_idx_string_subjet1.length() > 0)
       depth_and_idx_string_subjet1.append(".");
     depth_and_idx_string_subjet1.append("0");
     for (int iSpace = 0; iSpace < depth; ++iSpace) {
       std::cout << " ";
     }
     std::cout << " jetConstituent #" << depth_and_idx_string_subjet1 << " (depth = " << depth
               << "): Pt = " << subjet1.pt() << ","
               << " eta = " << subjet1.eta() << ", phi = " << subjet1.phi() << ", mass = " << subjet1.m()
               << " (constituents = " << subjet1.constituents().size() << ")" << std::endl;
     dumpSubJetStructure(subjet1, depth + 1, maxDepth, depth_and_idx_string_subjet1);
     std::string depth_and_idx_string_subjet2 = depth_and_idx_string;
     if (depth_and_idx_string_subjet2.length() > 0)
       depth_and_idx_string_subjet2.append(".");
     depth_and_idx_string_subjet2.append("1");
     for (int iSpace = 0; iSpace < depth; ++iSpace) {
       std::cout << " ";
     }
     std::cout << " jetConstituent #" << depth_and_idx_string_subjet2 << " (depth = " << depth
               << "): Pt = " << subjet2.pt() << ","
               << " eta = " << subjet2.eta() << ", phi = " << subjet2.phi() << ", mass = " << subjet2.m()
               << " (constituents = " << subjet2.constituents().size() << ")" << std::endl;
     dumpSubJetStructure(subjet2, depth + 1, maxDepth, depth_and_idx_string_subjet2);
     for (int iSpace = 0; iSpace < depth; ++iSpace) {
       std::cout << " ";
     }
     double dR = subjet1.delta_R(subjet2);
     std::cout << " (mass-drop @ " << depth_and_idx_string << " = " << std::max(subjet1.m(), subjet2.m()) / jet.m()
               << ", dR = " << dR << ")" << std::endl;
   }

References gather_cfg::cout, LEDCalibrationChannels::depth, HGC3DClusterGenMatchSelector_cfi::dR, dumpSubJetStructure(), metsig::jet, SiStripPI::max, HLT_FULL_cff::maxDepth, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dumpSubJetStructure(), and result().

◆ findSubjets()

std::pair< PseudoJet, PseudoJet > contrib::CMSBoostedTauSeedingAlgorithm::findSubjets	(	const PseudoJet &	jet,
		int	depth,
		bool &	subjetsFound
	)		const

protected

Definition at line 98 of file CMSBoostedTauSeedingAlgorithm.cc.

                                                                  {
     const float etaMax_ = 3.0;  // cut on eta subjets <3. Anyway taus are cut at 2.3
     if (verbosity_ >= 2) {
       std::cout << "<CMSBoostedTauSeedingAlgorithm::findSubjets>:" << std::endl;
       std::cout << " jet: Pt = " << jet.pt() << ", eta = " << jet.eta() << ", phi = " << jet.phi()
                 << ", mass = " << jet.m() << std::endl;
     }
  
     PseudoJet subjet1, subjet2;
     bool hasSubjets = jet.has_parents(subjet1, subjet2);
     if (hasSubjets && (maxDepth_ == -1 || depth <= maxDepth_)) {
       // make subjet1 the more massive jet
       if (subjet1.m2() < subjet2.m2()) {
         std::swap(subjet1, subjet2);
       }
       double dR = subjet1.delta_R(subjet2);
       double kT = subjet1.kt_distance(subjet2);
       double mu = (jet.m() > 0.) ? sqrt(std::max(subjet1.m2(), subjet2.m2()) / jet.m2()) : -1.;
       // check if subjets pass selection required for seeding boosted tau reconstruction
       if (subjet1.pt() > ptMin_ && fabs(subjet1.eta()) < etaMax_ && subjet2.pt() > ptMin_ &&
           fabs(subjet2.eta()) < etaMax_ && dR > dRMin_ && dR < dRMax_ && mu > muMin_ && mu < muMax_ &&
           kT < (yMax_ * jet.m2()) && kT > (yMin_ * jet.m2())) {
         subjetsFound = true;
         return std::make_pair(subjet1, subjet2);
       } else if (subjet1.pt() > ptMin_) {
         return findSubjets(subjet1, depth + 1, subjetsFound);
       } else if (subjet2.pt() > ptMin_) {
         return findSubjets(subjet2, depth + 1, subjetsFound);
       }
     }
     subjetsFound = false;
     PseudoJet dummy_subjet1, dummy_subjet2;
     return std::make_pair(dummy_subjet1, dummy_subjet2);
   }

References gather_cfg::cout, LEDCalibrationChannels::depth, HGC3DClusterGenMatchSelector_cfi::dR, dRMin_, metsig::jet, SiStripPI::max, maxDepth_, amptDefaultParameters_cff::mu, muMax_, muMin_, ptMin_, mathSSE::sqrt(), std::swap(), verbosity_, yMax_, and yMin_.

Referenced by result().

◆ result()

PseudoJet contrib::CMSBoostedTauSeedingAlgorithm::result ( const PseudoJet & jet ) const

override

Definition at line 134 of file CMSBoostedTauSeedingAlgorithm.cc.

                                                                             {
     if (verbosity_ >= 1) {
       std::cout << "<CMSBoostedTauSeedingAlgorithm::findSubjets>:" << std::endl;
       std::cout << " jet: Pt = " << jet.pt() << ", eta = " << jet.eta() << ", phi = " << jet.phi()
                 << ", mass = " << jet.m() << std::endl;
     }
  
     if (verbosity_ >= 2) {
       dumpSubJetStructure(jet, 0, maxDepth_, "");
     }
  
     bool subjetsFound = false;
     std::pair<PseudoJet, PseudoJet> subjets = findSubjets(jet, 0, subjetsFound);
     if (subjetsFound) {
       // fill structure for returning result
       PseudoJet subjet1 = subjets.first;
       PseudoJet subjet2 = subjets.second;
       if (verbosity_ >= 1) {
         std::cout << "before recombination:" << std::endl;
         std::cout << " subjet #1: Pt = " << subjet1.pt() << ", eta = " << subjet1.eta() << ", phi = " << subjet1.phi()
                   << ", mass = " << subjet1.m() << std::endl;
         std::cout << " subjet #2: Pt = " << subjet2.pt() << ", eta = " << subjet2.eta() << ", phi = " << subjet2.phi()
                   << ", mass = " << subjet2.m() << std::endl;
       }
  
       const JetDefinition::Recombiner* rec = jet.associated_cluster_sequence()->jet_def().recombiner();
       PseudoJet result_local = join(subjet1, subjet2, *rec);
       if (verbosity_ >= 1) {
         std::cout << "after recombination:" << std::endl;
         std::vector<fastjet::PseudoJet> subjets = result_local.pieces();
         int idx_subjet = 0;
         for (std::vector<fastjet::PseudoJet>::const_iterator subjet = subjets.begin(); subjet != subjets.end();
              ++subjet) {
           std::cout << " subjet #" << idx_subjet << ": Pt = " << subjet->pt() << ", eta = " << subjet->eta()
                     << ", phi = " << subjet->phi() << ", mass = " << subjet->m()
                     << " (#constituents = " << subjet->constituents().size() << ")" << std::endl;
           std::vector<fastjet::PseudoJet> constituents = subjet->constituents();
           int idx_constituent = 0;
           for (std::vector<fastjet::PseudoJet>::const_iterator constituent = constituents.begin();
                constituent != constituents.end();
                ++constituent) {
             if (constituent->pt() < 1.e-3)
               continue;  // CV: skip ghosts
             std::cout << "  constituent #" << idx_constituent << ": Pt = " << constituent->pt()
                       << ", eta = " << constituent->eta() << ", phi = " << constituent->phi() << ","
                       << " mass = " << constituent->m() << std::endl;
             ++idx_constituent;
           }
           ++idx_subjet;
         }
       }
  
       CMSBoostedTauSeedingAlgorithmStructure* s = new CMSBoostedTauSeedingAlgorithmStructure(result_local);
       //s->_original_jet = jet;
       s->_mu = (jet.m2() > 0.) ? sqrt(std::max(subjet1.m2(), subjet2.m2()) / jet.m2()) : 0.;
       s->_y = (jet.m2() > 0.) ? subjet1.kt_distance(subjet2) / jet.m2() : 0.;
       s->_dR = subjet1.delta_R(subjet2);
       s->_pt = subjet2.pt();
  
       result_local.set_structure_shared_ptr(SharedPtr<PseudoJetStructureBase>(s));
  
       return result_local;
     } else {
       // no subjets for seeding boosted tau reconstruction found, return an empty PseudoJet
       if (verbosity_ >= 1) {
         std::cout << "No subjets found." << std::endl;
       }
       return PseudoJet();
     }
   }