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 79 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.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;
     }
     
   }

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 48 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 65 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++;
       } 
     }
   }