fwrapper Namespace Reference

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
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, GetRecoTauVFromDQM_MC_cff::kk, 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++;
      } 
    }
  }

Variable Documentation

auto_ptr<ClusterSequence> fwrapper::cs

Definition at line 46 of file fastjetfortran_madfks.cc.

Referenced by cond::DataProxyWrapperBase::addInfo(), CSCGeometryBuilder::build(), CmsShowCommonPopup::changeGeomColor(), CmsShowCommonPopup::changeSelectionColorSet(), CmsShowModelPopup::clicked(), edmNew::copyDetSetRange(), edm::CosMuoGenSource::endRun(), edm::CosMuoGenProducer::endRun(), edm::ExhumeProducer::endRun(), CSCStripTopology::equationOfStrip(), fastjetarea_(), fastjetconstituents_(), fastjetdmerge_(), fastjetdmergemax_(), fastjetglobalrhoandsigma_(), FFTJetPatRecoProducer::FFTJetPatRecoProducer(), EcalDeadCellBoundaryEnergyFilter::filter(), AnomalousECALVariables::isDeadEcalCluster(), RadialStripTopology::localError(), edm::service::MessageLoggerScribe::log(), RadialStripTopology::measurementError(), ato< std::string >::operator()(), edm::DataFrameContainer::push_back(), SubjetFilterAlgorithm::run(), edm::DataFrameContainer::set_back(), cscdqm::Configuration::setChamberCounterValue(), Model::setCocoaStatus(), gen::Pythia6Hadronizer::statistics(), gen::ExhumeHadronizer::statistics(), transfer_cluster_transfer(), and L1RCTProducer::updateFedVector().

vector<PseudoJet> fwrapper::input_particles

Definition at line 43 of file fastjetfortran_madfks.cc.

Referenced by transfer_cluster_transfer(), and transfer_input_particles().

JetDefinition fwrapper::jet_def

Definition at line 45 of file fastjetfortran_madfks.cc.

Referenced by fastjetppgenkt_(), fastjetppgenktwitharea_(), fastjetsiscone_(), and fastjetsisconewitharea_().

vector<PseudoJet> fwrapper::jets

Definition at line 43 of file fastjetfortran_madfks.cc.

Referenced by TtFullHadKinFitter::addKinFitInfo(), TtSemiLepKinFitter::addKinFitInfo(), CaloMCTruthTreeProducer::analyze(), FlavorJetCorrectionExample::analyze(), PartonJetCorrectionExample::analyze(), PFMCTruthTreeProducer::analyze(), myFastSimVal::analyze(), JetCorExample< Jet >::analyze(), JetPlotsExample< Jet >::analyze(), TopJetAnalyzer::analyze(), JetAnaPythia< Jet >::analyze(), HLTTauDQML1Plotter::analyze(), TtFullHadSignalSelMVATrainer::analyze(), TtSemiLepSignalSelMVATrainer::analyze(), TtSemiLepJetCombMVATrainer::analyze(), TkAlCaRecoMonitor::analyze(), DQMHcalDiJetsAlCaReco::analyze(), cms::ProducerAnalyzer::analyze(), JetCorrectorOnTheFly< Jet >::analyze(), QcdHighPtDQM::analyze(), LHEAnalyzer::analyze(), ResolutionCreator::analyze(), cms::DiJetAnalyzer::analyze(), HLTInclusiveVBFSource::analyze(), cms::GammaJetAnalysis::analyze(), TTbar_GenJetAnalyzer::analyze(), SusyDQM< Mu, Ele, Jet, Met >::analyze(), TopValidation::analyze(), IsolatedTracksNxN::analyze(), GctFormatTranslateV35::blockToGctJetCandsAndCounts(), GctFormatTranslateV38::blockToGctJetCandsAndCounts(), GctFormatTranslateMCLegacy::blockToGctJetCandsAndCounts(), TtFullLepHypKinSolution::buildHypo(), TtSemiLepHypMaxSumPtWMass::buildHypo(), TtFullHadHypGenMatch::buildHypo(), TtFullLepHypGenMatch::buildHypo(), TtSemiLepHypGenMatch::buildHypo(), TtSemiLepHypGeom::buildHypo(), TtSemiLepHypMVADisc::buildHypo(), TtSemiLepHypWMassDeltaTopMass::buildHypo(), TtSemiLepHypWMassMaxSumPt::buildHypo(), TtSemiLepHypothesis::buildHypo(), AlignmentGlobalTrackSelector::checkIsolation(), AlignmentGlobalTrackSelector::checkJetCount(), TtSemiLepKinFitProducer< LeptonCollection >::doBTagging(), HLTTauRefProducer::doJets(), fastjetarea_(), fastjetconstituents_(), TopSingleLepton::MonitorEnsemble::fill(), TopHLTSingleLepton::MonitorEnsemble::fill(), TopDiLeptonOffline::MonitorEnsemble::fill(), myFilter::filter(), JetIDFailureFilter::filter(), TrackingFailureFilter::filter(), HLTRHemisphere::filter(), HeavyChHiggsToTauNuSkim::filter(), EgammaProbeSelector::filter(), pat::PATMHTProducer::getJets(), HLTJetPairDzMatchFilter< T >::hltFilter(), HLTAlphaTFilter< T >::hltFilter(), HLTFiltCand::hltFilter(), HLTJetSortedVBFFilter< T >::hltFilter(), HLTTauRefProducer::HLTTauRefProducer(), reco::isodeposit::OtherJetConstituentsDeltaRVeto::initialize(), TtSemiLepJetCombGeom::isValid(), TtFullHadHypothesis::isValid(), TtFullLepHypothesis::isValid(), TtSemiLepHypothesis::isValid(), FFTJetProducer::iterateJetReconstruction(), main(), lhef::JetMatchingMLM::match(), TtDilepLRSignalSelObservables::operator()(), TtSemiLRJetCombObservables::operator()(), egHLT::EgJetTagProbeCut< T >::pass(), egHLT::EgJetB2BCut< T >::pass(), TtSemiLeptonicEvent::print(), TtFullLeptonicEvent::print(), TtFullHadronicEvent::print(), TtSemiLepJetCombWMassMaxSumPt::produce(), TtSemiLepJetCombMaxSumPtWMass::produce(), HLTJetL1MatchProducer< T >::produce(), TtSemiLepJetCombWMassDeltaTopMass::produce(), AnyJetToCaloJetProducer::produce(), TtSemiLepJetCombGeom::produce(), TtDilepEvtSolutionMaker::produce(), TtSemiLepKinFitProducer< LeptonCollection >::produce(), StEvtSolutionMaker::produce(), HLTMhtProducer::produce(), TtSemiLepJetCombMVAComputer::produce(), cms::SubEventGenJetProducer::produce(), TtFullHadSignalSelMVAComputer::produce(), TtSemiLepSignalSelMVAComputer::produce(), HLTHtMhtProducer::produce(), L2TauPixelIsoTagProducer::produce(), TtFullLepKinSolutionProducer::produce(), TtSemiLepHitFitProducer< LeptonCollection >::produce(), RecoTauJetRegionProducer::produce(), TtSemiEvtSolutionMaker::produce(), RecoTauPiZeroFlattener::produce(), TtFullHadKinFitProducer::produce(), cms::JetVertexAssociation::produce(), TtHadEvtSolutionMaker::produce(), cms::JetCorrectionProducer< T >::produce(), SoftPFElectronTagInfoProducer::produce(), TtFullHadHypothesis::produce(), TtSemiLepHypothesis::produce(), AlCaGammaJetProducer::produce(), SoftPFMuonTagInfoProducer::produce(), TtFullLepHypothesis::produce(), cms::AlCaDiJetsProducer::produce(), ShiftedPFCandidateProducerForNoPileUpPFMEt::produce(), RecoTauProducer::produce(), pat::xxxEventHypothesisProducer::produce(), pat::PATJetProducer::produce(), TtJetPartonMatch< C >::produce(), cms::METProducer::produce(), HiL1Subtractor::produce(), HiGenCleaner< T2 >::produce(), JetIdSelector< T >::produce(), pat::JetCorrFactorsProducer::produce(), JetCrystalsAssociator::produce(), FastPrimaryVertexProducer::produce(), SoftLepton::produce(), BoostedTopProducer::produce(), CaloJetMETcorrInputProducerT< T, Textractor >::produce(), PFJetMETcorrInputProducerT< T, Textractor >::produce(), SmearedJetProducerT< T, Textractor >::produce(), FastjetJetProducer::produceTrackJets(), HITRegionalPixelSeedGenerator::regions(), FastJetFWLiteWrapper::run(), PFSpecificAlgo::runSignificance(), transfer_cluster_transfer(), transfer_jets(), GctFormatTranslateMCLegacy::writeGctOutJetBlock(), VirtualJetProducer::writeJets(), and FFTJetProducer::writeJets().

auto_ptr<JetDefinition::Plugin> fwrapper::plugin

Definition at line 44 of file fastjetfortran_madfks.cc.