Functions
void	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	transfer_input_particles (const double *p, const int &npart)
	helper routine to transfer fortran input particles into More...

void	transfer_jets (double *f77jets, int &njets)
	helper routine to help transfer jets -> f77jets[4*ijet+0..3] More...

Variables
unique_ptr< ClusterSequence >	cs

vector< PseudoJet >	input_particles

JetDefinition	jet_def

vector< PseudoJet >	jets

unique_ptr< JetDefinition::Plugin >	plugin

Detailed Description

a namespace for the fortran-wrapper which contains commonly-used structures and means to transfer fortran <-> C++

Function Documentation

◆ transfer_cluster_transfer()

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

helper routine packaging the transfers, the clustering and the extraction of the jets

Definition at line 81 of file fastjetfortran_madfks.cc.

                                                                  {
     // transfer p[4*ipart+0..3] -> input_particles[i]
     transfer_input_particles(p, npart);
  
     // perform the clustering
     if (ghost_maxrap == 0.0) {
       // cluster without areas
       cs = std::make_unique<ClusterSequence>(input_particles, jet_def);
     } else {
       // cluster with areas
       GhostedAreaSpec area_spec(ghost_maxrap, nrepeat, ghost_area);
       AreaDefinition area_def(active_area, area_spec);
       cs = std::make_unique<ClusterSequenceArea>(input_particles, jet_def, area_def);
     }
     // extract jets (pt-ordered)
     jets = sorted_by_pt(cs->inclusive_jets(ptmin));
  
     // transfer jets -> f77jets[4*ijet+0..3]
     transfer_jets(f77jets, njets);
  
     // Determine which parton/particle ended-up in which jet
     // set all jet entrie to zero first
     for (int ii = 0; ii < npart; ++ii)
       whichjet[ii] = 0;
  
     // Loop over jets and find constituents
     for (int kk = 0; kk < njets; ++kk) {
       vector<PseudoJet> constit = cs->constituents(jets[kk]);
       for (unsigned int ll = 0; ll < constit.size(); ++ll)
         whichjet[constit[ll].user_index()] = kk + 1;
     }
   }

References cs, cuy::ii, input_particles, jet_def, jets, GetRecoTauVFromDQM_MC_cff::kk, BTaggingMonitoring_cff::njets, npart, AlCaHLTBitMon_ParallelJobs::p, ptmin, transfer_input_particles(), and transfer_jets().

Referenced by fastjetppgenkt_().

◆ transfer_input_particles()

void fwrapper::transfer_input_particles	(	const double *	p,
		const int &	npart
	)

helper routine to transfer fortran input particles into

Definition at line 50 of file fastjetfortran_madfks.cc.

                                                                    {
     input_particles.resize(0);
     input_particles.reserve(npart);
     for (int i = 0; i < npart; i++) {
       valarray<double> mom(4);  // mom[0..3]
       for (int j = 0; j <= 3; j++) {
         mom[j] = *(p++);
         // RF-MZ: reorder the arguments because in madfks energy goes first; in fastjet last.
         // mom[(j+3) % 4] = *(p++);
       }
       PseudoJet psjet(mom);
       psjet.set_user_index(i);
       input_particles.push_back(psjet);
     }
   }

References mps_fire::i, input_particles, dqmiolumiharvest::j, npart, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by transfer_cluster_transfer().

◆ transfer_jets()

void fwrapper::transfer_jets	(	double *	f77jets,
		int &	njets
	)

helper routine to help transfer jets -> f77jets[4*ijet+0..3]

Definition at line 67 of file fastjetfortran_madfks.cc.

                                                   {
     njets = jets.size();
     for (int i = 0; i < njets; i++) {
       for (int j = 0; j <= 3; j++) {
         *f77jets = jets[i][j];
         // RF-MZ: reorder the arguments because in madfks energy goes first; in fastjet last.
         //*f77jets = jets[i][(j+3) % 4];
         f77jets++;
       }
     }
   }