fastjetfortran_madfks.cc File Reference

#include <iostream>
#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 386 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 358 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 409 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 426 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 443 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 257 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 320 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 156 of file fastjetfortran_madfks.cc.

                                                                    {
    
    // prepare jet def
    plugin.reset(new 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 205 of file fastjetfortran_madfks.cc.

                                                                    {
    
    // prepare jet def
    plugin.reset(new 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.