fastjetfortran_madfks.cc File Reference

#include <iostream>
#include <memory>
#include "fastjet/ClusterSequence.hh"
#include "fastjet/ClusterSequenceArea.hh"
#include "fastjet/Selector.hh"
#include "fastjet/SISConePlugin.hh"

Go to the source code of this file.

Namespaces
	fwrapper

Functions
double	fastjetarea_ (const int &ijet)
void	fastjetconstituents_ (const int &ijet, int *constituent_indices, int &nconstituents)
double	fastjetdmerge_ (const int &n)
double	fastjetdmergemax_ (const int &n)
void	fastjetglobalrhoandsigma_ (const double &rapmin, const double &rapmax, const double &phimin, const double &phimax, double &rho, double &sigma, double &meanarea)
void	fastjetppgenkt_ (const double *p, const int &npart, const double &R, const double &ptjetmin, const double &palg, double *f77jets, int &njets, int *whichjet)
void	fastjetppgenktwitharea_ (const double *p, const int &npart, const double &R, const double &palg, const double &ghost_rapmax, const int &nrepeat, const double &ghost_area, double *f77jets, int &njets, int *whichjet)
void	fastjetsiscone_ (const double *p, const int &npart, const double &R, const double &f, double *f77jets, int &njets, int *whichjet)
void	fastjetsisconewitharea_ (const double *p, const int &npart, const double &R, const double &f, const double &ghost_rapmax, const int &nrepeat, const double &ghost_area, double *f77jets, int &njets, int *whichjet)
void	fwrapper::transfer_cluster_transfer (const double *p, const int &npart, const JetDefinition &jet_def, const double &ptmin, double *f77jets, int &njets, int *whichjet, const double &ghost_maxrap=0.0, const int &nrepeat=0, const double &ghost_area=0.0)
void	fwrapper::transfer_input_particles (const double *p, const int &npart)
	helper routine to transfer fortran input particles into More...
void	fwrapper::transfer_jets (double *f77jets, int &njets)
	helper routine to help transfer jets -> f77jets[4*ijet+0..3] More...

Variables
unique_ptr< ClusterSequence >	fwrapper::cs
vector< PseudoJet >	fwrapper::input_particles
JetDefinition	fwrapper::jet_def
vector< PseudoJet >	fwrapper::jets
unique_ptr< JetDefinition::Plugin >	fwrapper::plugin

Function Documentation

fastjetarea_()

double fastjetarea_

(

const int &

ijet

)

f77 interface to provide access to the area of a jet found in the jet clustering with one of the above "...witharea" routines.

Given the index ijet of a jet (in the range 1...njets) obtained in the last call to jet clustering, return its area. If the jets have not been obtained with a "...witharea" soutine it returns 0.

Definition at line 392 of file fastjetfortran_madfks.cc.

                                     {
  assert(ijet > 0 && ijet <= int(jets.size()));
  const ClusterSequenceAreaBase *csab = dynamic_cast<const ClusterSequenceAreaBase *>(cs.get());
  if (csab != nullptr) {
    // we have areas and can use csab to access all the area-related info
    return csab->area(jets[ijet - 1]);
  } else {
    return 0.;
    //  Error("No area information associated to this jet.");
  }
}

References cms::cuda::assert(), fwrapper::cs, and fwrapper::jets.

fastjetconstituents_()

void fastjetconstituents_	(	const int &	ijet,
		int *	constituent_indices,
		int &	nconstituents
	)

f77 interface to provide access to the constituents of a jet found in the jet clustering with one of the above routines.

Given the index ijet of a jet (in the range 1...njets) obtained in the last call to jet clustering, fill the array constituent_indices, with nconstituents entries, with the indices of the constituents that belong to that jet (which will be in the range 1...npart)

Definition at line 366 of file fastjetfortran_madfks.cc.

                                                                                         {
  assert(cs.get() != nullptr);
  assert(ijet > 0 && ijet <= int(jets.size()));
 
  vector<PseudoJet> constituents = cs->constituents(jets[ijet - 1]);
 
  nconstituents = constituents.size();
  for (int i = 0; i < nconstituents; i++) {
    constituent_indices[i] = constituents[i].cluster_hist_index() + 1;
  }
}

References cms::cuda::assert(), fwrapper::cs, mps_fire::i, and fwrapper::jets.

fastjetdmerge_()

double fastjetdmerge_

(

const int &

n

)

return the dmin corresponding to the recombination that went from n+1 to n jets (sometimes known as d_{n n+1}).

Definition at line 413 of file fastjetfortran_madfks.cc.

                                    {
  assert(cs.get() != nullptr);
  return cs->exclusive_dmerge(n);
}

References cms::cuda::assert(), fwrapper::cs, and dqmiodumpmetadata::n.

fastjetdmergemax_()

double fastjetdmergemax_

(

const int &

n

)

return the maximum of the dmin encountered during all recombinations up to the one that led to an n-jet final state; identical to exclusive_dmerge, except in cases where the dmin do not increase monotonically.

Definition at line 429 of file fastjetfortran_madfks.cc.

                                       {
  assert(cs.get() != nullptr);
  return cs->exclusive_dmerge_max(n);
}

References cms::cuda::assert(), fwrapper::cs, and dqmiodumpmetadata::n.

fastjetglobalrhoandsigma_()

void fastjetglobalrhoandsigma_	(	const double &	rapmin,
		const double &	rapmax,
		const double &	phimin,
		const double &	phimax,
		double &	rho,
		double &	sigma,
		double &	meanarea
	)

return the background transverse momentum density per unit scalar area rho, its fluctuation sigma, and the mean area of the jets used for the background estimation in a given event, as evaluated in the range [rapmin,rapmax] in rapidity and [phimin,phimax] in azimuth

Definition at line 445 of file fastjetfortran_madfks.cc.

                                                 {
  const ClusterSequenceAreaBase *csab = dynamic_cast<const ClusterSequenceAreaBase *>(cs.get());
  if (csab != nullptr) {
    // we have areas and can use csab to access all the area-related info
    Selector range = SelectorRapRange(rapmin, rapmax) && SelectorPhiRange(phimin, phimax);
    bool use_area_4vector = false;
    csab->get_median_rho_and_sigma(range, use_area_4vector, rho, sigma, meanarea);
  } else {
    Error("Clustering with area is necessary in order to be able to evaluate rho.");
  }
}

References fwrapper::cs, phimax, phimin, and FastTimerService_cff::range.

fastjetppgenkt_()

void fastjetppgenkt_	(	const double *	p,
		const int &	npart,
		const double &	R,
		const double &	ptjetmin,
		const double &	palg,
		double *	f77jets,
		int &	njets,
		int *	whichjet
	)

f77 interface to the pp generalised-kt (sequential recombination) algorithms, as defined in arXiv.org:0802.1189, which includes kt, Cambridge/Aachen and anti-kt as special cases.

Definition at line 259 of file fastjetfortran_madfks.cc.

                                    {
  // prepare jet def
  if (palg == 1.0) {
    jet_def = JetDefinition(kt_algorithm, R);
  } else if (palg == 0.0) {
    jet_def = JetDefinition(cambridge_algorithm, R);
  } else if (palg == -1.0) {
    jet_def = JetDefinition(antikt_algorithm, R);
  } else {
    jet_def = JetDefinition(genkt_algorithm, R, palg);
  }
 
  // do everything
  transfer_cluster_transfer(p, npart, jet_def, ptjetmin, f77jets, njets, whichjet);
}

References fwrapper::jet_def, BTaggingMonitoring_cff::njets, npart, AlCaHLTBitMon_ParallelJobs::p, dttmaxenums::R, and fwrapper::transfer_cluster_transfer().

fastjetppgenktwitharea_()

void fastjetppgenktwitharea_	(	const double *	p,
		const int &	npart,
		const double &	R,
		const double &	palg,
		const double &	ghost_rapmax,
		const int &	nrepeat,
		const double &	ghost_area,
		double *	f77jets,
		int &	njets,
		int *	whichjet
	)

f77 interface to the pp generalised-kt (sequential recombination) algorithms, as defined in arXiv.org:0802.1189, which includes kt, Cambridge/Aachen and anti-kt as special cases. Also calculates the active area of the jets, as defined in arXiv.org:0802.1188

Definition at line 324 of file fastjetfortran_madfks.cc.

                                            {
  // prepare jet def
  if (palg == 1.0) {
    jet_def = JetDefinition(kt_algorithm, R);
  } else if (palg == 0.0) {
    jet_def = JetDefinition(cambridge_algorithm, R);
  } else if (palg == -1.0) {
    jet_def = JetDefinition(antikt_algorithm, R);
  } else {
    jet_def = JetDefinition(genkt_algorithm, R, palg);
  }
 
  // do everything
  //    transfer_cluster_transfer(p,npart,jet_def,f77jets,njets,whichjet,ghost_rapmax,nrepeat,ghost_area);
}

References fwrapper::jet_def, and dttmaxenums::R.

fastjetsiscone_()

void fastjetsiscone_	(	const double *	p,
		const int &	npart,
		const double &	R,
		const double &	f,
		double *	f77jets,
		int &	njets,
		int *	whichjet
	)

f77 interface to SISCone (via fastjet), as defined in arXiv:0704.0292 [see below for the interface to kt, Cam/Aachen & kt]

Definition at line 159 of file fastjetfortran_madfks.cc.

                                                                                                                     {
  // prepare jet def
  plugin = std::make_unique<SISConePlugin>(R, f);
  jet_def = plugin.get();
 
  // do everything
  //    transfer_cluster_transfer(p,npart,jet_def,f77jets,njets,whichjet);
}

References f, fwrapper::jet_def, fwrapper::plugin, and dttmaxenums::R.

fastjetsisconewitharea_()

void fastjetsisconewitharea_	(	const double *	p,
		const int &	npart,
		const double &	R,
		const double &	f,
		const double &	ghost_rapmax,
		const int &	nrepeat,
		const double &	ghost_area,
		double *	f77jets,
		int &	njets,
		int *	whichjet
	)

f77 interface to SISCone (via fastjet), as defined in arXiv:0704.0292 [see below for the interface to kt, Cam/Aachen & kt] Also calculates the active area of the jets, as defined in arXiv.org:0802.1188

Definition at line 204 of file fastjetfortran_madfks.cc.

                                            {
  // prepare jet def
  plugin = std::make_unique<SISConePlugin>(R, f);
  jet_def = plugin.get();
 
  // do everything
  //    transfer_cluster_transfer(p,npart,jet_def,f77jets,njets,whichjet,ghost_rapmax,nrepeat,ghost_area);
}

References f, fwrapper::jet_def, fwrapper::plugin, and dttmaxenums::R.