#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)

virtual std::string	description () const

virtual PseudoJet	result (const PseudoJet &jet) const

virtual	~CMSBoostedTauSeedingAlgorithm ()

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 53 of file CMSBoostedTauSeedingAlgorithm.h.

Member Typedef Documentation

typedef CMSBoostedTauSeedingAlgorithmStructure contrib::CMSBoostedTauSeedingAlgorithm::StructureType

Definition at line 68 of file CMSBoostedTauSeedingAlgorithm.h.

Constructor & Destructor Documentation

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 32 of file CMSBoostedTauSeedingAlgorithm.cc.

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

virtual contrib::CMSBoostedTauSeedingAlgorithm::~CMSBoostedTauSeedingAlgorithm ( )

inlinevirtual

Definition at line 60 of file CMSBoostedTauSeedingAlgorithm.h.

References description(), metsig::jet, result(), and AlCaHLTBitMon_QueryRunRegistry::string.

60 {}

Member Function Documentation

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

virtual

Definition at line 48 of file CMSBoostedTauSeedingAlgorithm.cc.

Referenced by ~CMSBoostedTauSeedingAlgorithm().

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

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

protected

Definition at line 56 of file CMSBoostedTauSeedingAlgorithm.cc.

References gather_cfg::cout, particleFlowClusterECALTimeSelected_cfi::depth, PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, dumpSubJetStructure(), hpstanc_transforms::max, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dumpSubJetStructure(), and result().

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

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

protected

Definition at line 89 of file CMSBoostedTauSeedingAlgorithm.cc.

References gather_cfg::cout, PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, dRMin_, hpstanc_transforms::max, maxDepth_, RPCpg::mu, muMax_, muMin_, ptMin_, mathSSE::sqrt(), std::swap(), verbosity_, yMax_, and yMin_.

Referenced by result().

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

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

virtual

Definition at line 123 of file CMSBoostedTauSeedingAlgorithm.cc.

References contrib::CMSBoostedTauSeedingAlgorithmStructure::_dR, contrib::CMSBoostedTauSeedingAlgorithmStructure::_mu, contrib::CMSBoostedTauSeedingAlgorithmStructure::_pt, contrib::CMSBoostedTauSeedingAlgorithmStructure::_y, gather_cfg::cout, dumpSubJetStructure(), findSubjets(), join(), hpstanc_transforms::max, maxDepth_, alignCSCRings::s, mathSSE::sqrt(), and verbosity_.

Referenced by ~CMSBoostedTauSeedingAlgorithm().

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