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
void	transfer_jets (double *f77jets, int &njets)
	helper routine to help transfer jets -> f77jets[4*ijet+0..3]
Variables
auto_ptr< ClusterSequence >	cs
vector< PseudoJet >	input_particles
JetDefinition	jet_def
vector< PseudoJet >	jets
auto_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

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 80 of file fastjetfortran_madfks.cc.

References cs, input_particles, jets, npart, transfer_input_particles(), and transfer_jets().

Referenced by fastjetppgenkt_().

                                                                                           {

    // 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.reset(new 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.reset(new 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;
    }
    
  }

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

helper routine to transfer fortran input particles into

Definition at line 49 of file fastjetfortran_madfks.cc.

References i, input_particles, j, and npart.

Referenced by transfer_cluster_transfer().

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

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

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

Definition at line 66 of file fastjetfortran_madfks.cc.

References i, j, and jets.

Referenced by transfer_cluster_transfer().

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