SubjetFilterAlgorithm Class Reference

#include <SubjetFilterAlgorithm.h>

Public Member Functions
void	run (const std::vector< fastjet::PseudoJet > &inputs, std::vector< CompoundPseudoJet > &fatJets, const edm::EventSetup &iSetup)
	SubjetFilterAlgorithm (const std::string &moduleLabel, const std::string &jetAlgorithm, unsigned nFatMax, double rParam, double rFilt, double jetPtMin, double massDropCut, double asymmCut, bool asymmCutLater, bool doAreaFastjet, double ghostEtaMax, int activeAreaRepeats, double ghostArea, bool verbose)
std::string	summary () const
virtual	~SubjetFilterAlgorithm ()

Private Attributes
int	activeAreaRepeats_
double	asymmCut2_
bool	asymmCutLater_
bool	doAreaFastjet_
fastjet::AreaDefinition *	fjAreaDef_
fastjet::JetDefinition *	fjJetDef_
double	ghostArea_
double	ghostEtaMax_
std::string	jetAlgorithm_
double	jetPtMin_
double	massDropCut_
std::string	moduleLabel_
unsigned	nevents_
unsigned	nFatMax_
unsigned	nfound_
unsigned	ntotal_
double	rFilt_
double	rParam_
bool	verbose_

Detailed Description

Definition at line 29 of file SubjetFilterAlgorithm.h.

Constructor & Destructor Documentation

SubjetFilterAlgorithm::SubjetFilterAlgorithm	(	const std::string &	moduleLabel,
		const std::string &	jetAlgorithm,
		unsigned	nFatMax,
		double	rParam,
		double	rFilt,
		double	jetPtMin,
		double	massDropCut,
		double	asymmCut,
		bool	asymmCutLater,
		bool	doAreaFastjet,
		double	ghostEtaMax,
		int	activeAreaRepeats,
		double	ghostArea,
		bool	verbose
	)

Definition at line 40 of file SubjetFilterAlgorithm.cc.

References activeAreaRepeats_, doAreaFastjet_, Exception, fjAreaDef_, fjJetDef_, ghostArea_, ghostEtaMax_, jetAlgorithm_, and rParam_.

  : moduleLabel_(moduleLabel)
  , jetAlgorithm_(jetAlgorithm)
  , nFatMax_(nFatMax)
  , rParam_(rParam)
  , rFilt_(rFilt)
  , jetPtMin_(jetPtMin)
  , massDropCut_(massDropCut)
  , asymmCut2_(asymmCut*asymmCut)
  , asymmCutLater_(asymmCutLater)
  , doAreaFastjet_(doAreaFastjet)
  , ghostEtaMax_(ghostEtaMax)
  , activeAreaRepeats_(activeAreaRepeats)
  , ghostArea_(ghostArea)
  , verbose_(verbose)
  , nevents_(0)
  , ntotal_(0)
  , nfound_(0)
  , fjJetDef_(nullptr)
  , fjAreaDef_(nullptr)
{
  // FASTJET JET DEFINITION
  if (jetAlgorithm=="CambridgeAachen"||jetAlgorithm_=="ca")
    fjJetDef_=new fastjet::JetDefinition(fastjet::cambridge_algorithm,rParam_);
  else if (jetAlgorithm=="AntiKt"||jetAlgorithm_=="ak")
    fjJetDef_=new fastjet::JetDefinition(fastjet::antikt_algorithm,rParam_);
  else if (jetAlgorithm=="Kt"||jetAlgorithm_=="kt")
    fjJetDef_=new fastjet::JetDefinition(fastjet::kt_algorithm,rParam_);
  else
    throw cms::Exception("InvalidJetAlgo")
      <<"Jet Algorithm for SubjetFilterAlgorithm is invalid: "
      <<jetAlgorithm_<<", use (ca|CambridgeAachen)|(Kt|kt)|(AntiKt|ak)"<<endl;

  // FASTJET AREA DEFINITION
  if (doAreaFastjet_) {
    fastjet::GhostedAreaSpec ghostSpec(ghostEtaMax_,activeAreaRepeats_,ghostArea_);
    fjAreaDef_=new fastjet::AreaDefinition(fastjet::active_area_explicit_ghosts,
                       ghostSpec); 
  }
}

SubjetFilterAlgorithm::~SubjetFilterAlgorithm ( )

virtual

Definition at line 96 of file SubjetFilterAlgorithm.cc.

References fjAreaDef_, and fjJetDef_.

{
  if (nullptr!=fjJetDef_)  delete fjJetDef_;
  if (nullptr!=fjAreaDef_) delete fjAreaDef_;
}

Member Function Documentation

void SubjetFilterAlgorithm::run	(	const std::vector< fastjet::PseudoJet > &	inputs,
		std::vector< CompoundPseudoJet > &	fatJets,
		const edm::EventSetup &	iSetup
	)

Definition at line 108 of file SubjetFilterAlgorithm.cc.

References custom_jme_cff::area, asymmCut2_, asymmCutLater_, gather_cfg::cout, fwrapper::cs, doAreaFastjet_, fjAreaDef_, fjJetDef_, mps_fire::i, jetPtMin_, massDropCut_, min(), nevents_, nFatMax_, nfound_, ntotal_, rFilt_, rParam_, mathSSE::sqrt(), edm::swap(), and verbose_.

Referenced by SubjetFilterJetProducer::runAlgorithm().

{
  nevents_++;
  
  if (verbose_) cout<<endl<<nevents_<<". EVENT"<<endl;
  
  fastjet::ClusterSequence* cs = (doAreaFastjet_) ?
    new fastjet::ClusterSequenceArea(fjInputs,*fjJetDef_,*fjAreaDef_) :
    new fastjet::ClusterSequence    (fjInputs,*fjJetDef_);
  
  vector<fastjet::PseudoJet> fjFatJets =
    fastjet::sorted_by_pt(cs->inclusive_jets(jetPtMin_));
  
  size_t nFat =
    (nFatMax_==0) ? fjFatJets.size() : std::min(fjFatJets.size(),(size_t)nFatMax_);
  
  for (size_t iFat=0;iFat<nFat;iFat++) {
    
    if (verbose_) cout<<endl<<iFat<<". FATJET: "<<fjFatJets[iFat]<<endl;
    
    fastjet::PseudoJet fjFatJet = fjFatJets[iFat];
    fastjet::PseudoJet fjCurrentJet(fjFatJet);
    fastjet::PseudoJet fjSubJet1,fjSubJet2;
    bool hadSubJets;
    
    vector<CompoundPseudoSubJet> subJets;
    
    
    // FIND SUBJETS PASSING MASSDROP [AND ASYMMETRY] CUT(S)
    while ((hadSubJets = cs->has_parents(fjCurrentJet,fjSubJet1,fjSubJet2))) {
      
      if (fjSubJet1.m() < fjSubJet2.m()) swap(fjSubJet1,fjSubJet2);
      
      if (verbose_) cout<<"SUBJET CANDIDATES: "<<fjSubJet1<<endl
            <<"                   "<<fjSubJet2<<endl;
      if (verbose_) cout<<"md="<<fjSubJet1.m()/fjCurrentJet.m()
            <<" y="<<fjSubJet1.kt_distance(fjSubJet2)/fjCurrentJet.m2()
            <<endl;
      
      if (fjSubJet1.m()<massDropCut_*fjCurrentJet.m() &&
      (asymmCutLater_||
       fjSubJet1.kt_distance(fjSubJet2)>asymmCut2_*fjCurrentJet.m2())) {
    break;
      }
      else {
    fjCurrentJet = fjSubJet1;
      }
    }
    
    if (!hadSubJets) {
      if (verbose_) cout<<"FAILED TO SELECT SUBJETS"<<endl;
    }
    // FOUND TWO GOOD SUBJETS PASSING MASSDROP CUT
    else {
      if (verbose_) cout<<"SUBJETS selected"<<endl;
      
      if (asymmCutLater_&&
      fjSubJet1.kt_distance(fjSubJet2)<=asymmCut2_*fjCurrentJet.m2()) {
    if (verbose_) cout<<"FAILED y-cut"<<endl;
      }
      // PASSED ASYMMETRY (Y) CUT
      else {
    if (verbose_) cout<<"PASSED y cut"<<endl;
    
    vector<fastjet::PseudoJet> fjFilterJets;
    double       Rbb   = std::sqrt(fjSubJet1.squared_distance(fjSubJet2));
    double       Rfilt = std::min(0.5*Rbb,rFilt_);
    double       dcut  = Rfilt*Rfilt/rParam_/rParam_;
    fjFilterJets=fastjet::sorted_by_pt(cs->exclusive_subjets(fjCurrentJet,dcut));
    
    if (verbose_) {
      cout<<"Rbb="<<Rbb<<", Rfilt="<<Rfilt<<endl;
      cout<<"FILTER JETS: "<<flush;
      for (size_t i=0;i<fjFilterJets.size();i++) {
        if (i>0) cout<<"             "<<flush; 
        cout<<fjFilterJets[i]<<endl;
      }
    }
    
    vector<fastjet::PseudoJet> fjSubJets;
    fjSubJets.push_back(fjSubJet1);
    fjSubJets.push_back(fjSubJet2);
    size_t nFilter = fjFilterJets.size();
    for (size_t iFilter=0;iFilter<nFilter;iFilter++)
      fjSubJets.push_back(fjFilterJets[iFilter]);
    
    for (size_t iSub=0;iSub<fjSubJets.size();iSub++) {
      vector<fastjet::PseudoJet> fjConstituents=
        cs->constituents(fjSubJets[iSub]);
      vector<int> constituents;
      for (size_t iConst=0;iConst<fjConstituents.size();iConst++) {
        int userIndex = fjConstituents[iConst].user_index();
        if (userIndex>=0)  constituents.push_back(userIndex);
      }

      double subJetArea=(doAreaFastjet_) ?
        ((fastjet::ClusterSequenceArea*)cs)->area(fjSubJets[iSub]) : 0.0;
      
      if (iSub<2||!constituents.empty())
        subJets.push_back(CompoundPseudoSubJet(fjSubJets[iSub],
                           subJetArea,
                           constituents));
    }
    
      } // PASSED Y CUT
      
    } // PASSED MASSDROP CUT
    
    if (verbose_) cout<<"write fatjet with "<<subJets.size()<<" sub+filter jets"
              <<endl;
    
    double fatJetArea = (doAreaFastjet_) ?
      ((fastjet::ClusterSequenceArea*)cs)->area(fjFatJet) : 0.0;

    fjJets.push_back(CompoundPseudoJet(fjFatJet,fatJetArea,subJets));
    
    ntotal_++;
    if (subJets.size()>3) nfound_++;
    
  } // LOOP OVER FATJETS
  
  if (verbose_) cout<<endl<<fjJets.size()<<" FATJETS written\n"<<endl;
  
  delete cs;
  
  return;
}

string SubjetFilterAlgorithm::summary

(

)

const

Definition at line 240 of file SubjetFilterAlgorithm.cc.

References moduleLabel_, nevents_, nfound_, and ntotal_.

Referenced by SubjetFilterJetProducer::endJob().

{
  double eff = (ntotal_>0) ? nfound_/(double)ntotal_ : 0;
  std::stringstream ss;
  ss<<"************************************************************\n"
    <<"* "<<moduleLabel_<<" (SubjetFilterAlgorithm) SUMMARY:\n"
    <<"************************************************************\n"
    <<"nevents = "<<nevents_<<endl
    <<"ntotal  = "<<ntotal_<<endl
    <<"nfound  = "<<nfound_<<endl
    <<"eff     = "<<eff<<endl
    <<"************************************************************\n";
  return ss.str();
}

Member Data Documentation

int SubjetFilterAlgorithm::activeAreaRepeats_

private

Definition at line 72 of file SubjetFilterAlgorithm.h.

Referenced by SubjetFilterAlgorithm().

double SubjetFilterAlgorithm::asymmCut2_

private

Definition at line 68 of file SubjetFilterAlgorithm.h.

Referenced by run().

bool SubjetFilterAlgorithm::asymmCutLater_

private

Definition at line 69 of file SubjetFilterAlgorithm.h.

Referenced by run().

bool SubjetFilterAlgorithm::doAreaFastjet_

private

Definition at line 70 of file SubjetFilterAlgorithm.h.

Referenced by run(), and SubjetFilterAlgorithm().

fastjet::AreaDefinition* SubjetFilterAlgorithm::fjAreaDef_

private

Definition at line 81 of file SubjetFilterAlgorithm.h.

Referenced by run(), SubjetFilterAlgorithm(), and ~SubjetFilterAlgorithm().

fastjet::JetDefinition* SubjetFilterAlgorithm::fjJetDef_

private

Definition at line 80 of file SubjetFilterAlgorithm.h.

Referenced by run(), SubjetFilterAlgorithm(), and ~SubjetFilterAlgorithm().

double SubjetFilterAlgorithm::ghostArea_

private

Definition at line 73 of file SubjetFilterAlgorithm.h.

Referenced by SubjetFilterAlgorithm().

double SubjetFilterAlgorithm::ghostEtaMax_

private

Definition at line 71 of file SubjetFilterAlgorithm.h.

Referenced by SubjetFilterAlgorithm().

std::string SubjetFilterAlgorithm::jetAlgorithm_

private

Definition at line 62 of file SubjetFilterAlgorithm.h.

Referenced by SubjetFilterAlgorithm().

double SubjetFilterAlgorithm::jetPtMin_

private

Definition at line 66 of file SubjetFilterAlgorithm.h.

Referenced by run().

double SubjetFilterAlgorithm::massDropCut_

private

Definition at line 67 of file SubjetFilterAlgorithm.h.

Referenced by run().

std::string SubjetFilterAlgorithm::moduleLabel_

private

Definition at line 61 of file SubjetFilterAlgorithm.h.

Referenced by ExternalGeneratorFilter.ExternalGeneratorFilter::insertInto(), Modules.SwitchProducer::insertInto(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), Modules.ESProducer::nameInProcessDesc_(), and summary().

unsigned SubjetFilterAlgorithm::nevents_

private

Definition at line 76 of file SubjetFilterAlgorithm.h.

Referenced by run(), and summary().

unsigned SubjetFilterAlgorithm::nFatMax_

private

Definition at line 63 of file SubjetFilterAlgorithm.h.

Referenced by run().

unsigned SubjetFilterAlgorithm::nfound_

private

Definition at line 78 of file SubjetFilterAlgorithm.h.

Referenced by run(), and summary().

unsigned SubjetFilterAlgorithm::ntotal_

private

Definition at line 77 of file SubjetFilterAlgorithm.h.

Referenced by run(), and summary().

double SubjetFilterAlgorithm::rFilt_

private

Definition at line 65 of file SubjetFilterAlgorithm.h.

Referenced by run().

double SubjetFilterAlgorithm::rParam_

private

Definition at line 64 of file SubjetFilterAlgorithm.h.

Referenced by run(), and SubjetFilterAlgorithm().

bool SubjetFilterAlgorithm::verbose_

private

Definition at line 74 of file SubjetFilterAlgorithm.h.

Referenced by run().