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

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.

References cuy::ii, 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 48 of file fastjetfortran_madfks.cc.

References mps_fire::i, 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 65 of file fastjetfortran_madfks.cc.

References mps_fire::i.

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

unique_ptr<ClusterSequence> fwrapper::cs

Definition at line 45 of file fastjetfortran_madfks.cc.

Referenced by PedestalTask::_dump(), RawTask::_process(), DigiPhase1Task::_process(), DigiTask::_process(), cond::DataProxyWrapperBase::addInfo(), DigiTask::bookHistograms(), CSCGeometryBuilder::build(), CmsShowCommonPopup::changeGeomColor(), CmsShowCommonPopup::changeSelectionColorSet(), CmsShowModelPopup::clicked(), XrdAdaptor::XrdStatisticsService::condorUpdate(), edmNew::copyDetSetRange(), HcalTimeSlewSim::delay(), edm::CosMuoGenProducer::endRunProduce(), CSCStripTopology::equationOfStrip(), fastjetarea_(), fastjetconstituents_(), fastjetdmerge_(), fastjetdmergemax_(), fastjetglobalrhoandsigma_(), FFTJetPatRecoProducer::FFTJetPatRecoProducer(), LHEJetFilter::filter(), EcalDeadCellBoundaryEnergyFilter::filter(), gen::Pythia6Hadronizer::finalizeEvent(), KinematicConstrainedVertexFitter::fit(), KinematicConstrainedVertexFitterT< nTrk, nConstraint >::fit(), BPHRecoBuilder::BPHSpecificCollection< T >::get(), AnomalousECALVariables::isDeadEcalCluster(), CSCRadialStripTopology::localError(), TkRadialStripTopology::localError(), CSCRadialStripTopology::measurementError(), TkRadialStripTopology::measurementError(), l1t::ClusterShapes::operator+(), HBHEPhase1Reconstructor::processData(), edm::DataFrameContainer::push_back(), Qjets::Qjets(), HFPreRecAlgo::reconstruct(), HcalSimpleRecAlgoImpl::removePileup(), SubjetFilterAlgorithm::run(), BPHRecoBuilder::sameTrack(), edm::DataFrameContainer::set_back(), cscdqm::Configuration::setChamberCounterValue(), Model::setCocoaStatus(), ESElectronicsSim::setMIPToGeV(), HBHENegativeFlagSetter::setPulseShapeFlags(), QjetsPlugin::SetRNEngine(), gen::ExhumeHadronizer::statistics(), gen::Pythia6Hadronizer::statistics(), L1RCTProducer::updateFedVector(), LagrangeParentParticleFitter::~LagrangeParentParticleFitter(), and ParentParticleFitter::~ParentParticleFitter().

vector<PseudoJet> fwrapper::input_particles

Definition at line 42 of file fastjetfortran_madfks.cc.

JetDefinition fwrapper::jet_def

Definition at line 44 of file fastjetfortran_madfks.cc.

Referenced by fastjetppgenkt_(), fastjetppgenktwitharea_(), fastjetsiscone_(), fastjetsisconewitharea_(), and heppy::ReclusterJets::ReclusterJets().

vector<PseudoJet> fwrapper::jets

Definition at line 42 of file fastjetfortran_madfks.cc.

Referenced by TtFullHadKinFitter::addKinFitInfo(), TtSemiLepKinFitter::addKinFitInfo(), PartonJetCorrectionExample::analyze(), FlavorJetCorrectionExample::analyze(), CaloMCTruthTreeProducer::analyze(), PFMCTruthTreeProducer::analyze(), myFastSimVal::analyze(), JetCorExample< Jet >::analyze(), JetPlotsExample< Jet >::analyze(), JetAnaPythia< Jet >::analyze(), PatZjetsJetAnalyzer::analyze(), TopJetAnalyzer::analyze(), QcdHighPtDQM::analyze(), PatTopSelectionAnalyzer::analyze(), JetResolutionDemo::analyze(), PatBasicAnalyzer::analyze(), PatBTagAnalyzer::analyze(), DQMHcalDiJetsAlCaReco::analyze(), TtFullHadSignalSelMVATrainer::analyze(), JetCorrectorOnTheFly< Jet >::analyze(), TtSemiLepSignalSelMVATrainer::analyze(), l1t::L1TStage2CaloAnalyzer::analyze(), TtSemiLepJetCombMVATrainer::analyze(), TkAlCaRecoMonitor::analyze(), ResolutionCreator::analyze(), cms::ProducerAnalyzer::analyze(), PatJetAnalyzer::analyze(), HLTInclusiveVBFSource::analyze(), B2GHadronicHLTValidation::analyze(), B2GSingleLeptonHLTValidation::analyze(), TopSingleLeptonHLTValidation::analyze(), l1t::GtInputDump::analyze(), METplusTrackMonitor::analyze(), TopDiLeptonHLTValidation::analyze(), TTbar_GenJetAnalyzer::analyze(), SusyDQM< Mu, Ele, Jet, Met >::analyze(), ObjMonitor::analyze(), l1t::L1TGlobalAnalyzer::analyze(), NoBPTXMonitor::analyze(), TopMonitor::analyze(), l1t::GtRecordDump::analyze(), PhotonMonitor::analyze(), RazorMonitor::analyze(), CastorMonitorModule::analyze(), METMonitor::analyze(), HTMonitor::analyze(), IsolatedTracksNxN::analyze(), FFTJetCorrectionProducer::applyCorrections(), GctFormatTranslateV35::blockToGctJetCandsAndCounts(), GctFormatTranslateV38::blockToGctJetCandsAndCounts(), GctFormatTranslateMCLegacy::blockToGctJetCandsAndCounts(), TtFullLepHypKinSolution::buildHypo(), TtSemiLepHypWMassDeltaTopMass::buildHypo(), TtFullHadHypGenMatch::buildHypo(), TtFullLepHypGenMatch::buildHypo(), TtSemiLepHypGenMatch::buildHypo(), TtSemiLepHypGeom::buildHypo(), TtSemiLepHypMaxSumPtWMass::buildHypo(), TtSemiLepHypMVADisc::buildHypo(), TtSemiLepHypWMassMaxSumPt::buildHypo(), TtSemiLepHypothesis::buildHypo(), TtFullLepHypKinSolution::buildKey(), TtFullHadHypGenMatch::buildKey(), TtFullLepHypGenMatch::buildKey(), HiFJGridEmptyAreaCalculator::calculateAreaFractionOfJets(), PrimaryVertexAssignment::chargedHadronVertex(), AlignmentGlobalTrackSelector::checkIsolation(), AlignmentGlobalTrackSelector::checkJetCount(), Rivet::HiggsTemplateCrossSections::classifyEvent(), L1TStage2CaloLayer2Comp::compareJets(), RazorComputer::compute(), TtSemiLepKinFitProducer< LeptonCollection >::doBTagging(), HLTTauRefProducer::doJets(), l1t::GtRecordDump::dumpTestVectors(), fastjetarea_(), fastjetconstituents_(), SingleTopTChannelLepton_miniAOD::MonitorEnsemble::fill(), TopSingleLepton_miniAOD::MonitorEnsemble::fill(), HLTOfflineDQMTopSingleLepton::MonitorSingleLepton::fill(), HLTOfflineDQMTopDiLepton::MonitorDiLepton::fill(), SingleTopTChannelLepton::MonitorEnsemble::fill(), TopSingleLepton::MonitorEnsemble::fill(), TopDiLeptonOffline::MonitorEnsemble::fill(), AdHocNTupler::fill(), JetIDFailureFilter::filter(), TrackingFailureFilter::filter(), LHEJetFilter::filter(), HLTRHemisphere::filter(), HeavyChHiggsToTauNuSkim::filter(), pat::PATMHTProducer::getJets(), HLTJetPairDzMatchFilter< T >::hltFilter(), HLTAlphaTFilter< T >::hltFilter(), HLTMinDPhiMETFilter::hltFilter(), HLTFiltCand::hltFilter(), HLTJetSortedVBFFilter< T >::hltFilter(), HLTJetHbbFilter< T >::hltFilter(), HLTTauRefProducer::HLTTauRefProducer(), reco::isodeposit::OtherJetConstituentsDeltaRVeto::initialize(), TtSemiLepJetCombGeom::isValid(), TtFullHadHypothesis::isValid(), TtFullLepHypothesis::isValid(), TtSemiLepHypothesis::isValid(), FFTJetProducer::iterateJetReconstruction(), MvaMEtUtilities::jetP4(), main(), IPProducerHelpers::FromJetAndCands::makeBaseVector(), HTXS::MakeTLVs(), TemplatedSecondaryVertexProducer< IPTI, VTX >::matchGroomedJets(), JetFlavourClustering::matchGroomedJets(), JetFlavourClustering::matchReclusteredJets(), TtDilepLRSignalSelObservables::operator()(), TtSemiLRJetCombObservables::operator()(), l1t::stage2::JetPacker::pack(), l1t::stage1::CentralJetPacker::pack(), l1t::stage1::ForwardJetPacker::pack(), egHLT::EgJetTagProbeCut< T >::pass(), egHLT::EgJetB2BCut< T >::pass(), TtSemiLeptonicEvent::print(), TtFullLeptonicEvent::print(), TtFullHadronicEvent::print(), TtSemiLepJetCombMaxSumPtWMass::produce(), PseudoTopProducer::produce(), ECFAdder::produce(), TtSemiLepJetCombWMassMaxSumPt::produce(), TtSemiLepJetCombGeom::produce(), StEvtSolutionMaker::produce(), TtDilepEvtSolutionMaker::produce(), TtSemiLepKinFitProducer< LeptonCollection >::produce(), ParticleLevelProducer::produce(), QGTagger::produce(), TtSemiLepJetCombWMassDeltaTopMass::produce(), cms::SubEventGenJetProducer::produce(), HLTJetL1MatchProducer< T >::produce(), AnyJetToCaloJetProducer::produce(), HLTJetL1TMatchProducer< T >::produce(), TtFullHadSignalSelMVAComputer::produce(), TtFullLepKinSolutionProducer::produce(), TtSemiLepHitFitProducer< LeptonCollection >::produce(), JetTracksAssociationToTrackRefs::produce(), PileupJPTJetIdProducer::produce(), PFJetsMaxInvMassModule::produce(), PATPuppiJetSpecificProducer::produce(), MVAJetPuIdProducer::produce(), TtSemiEvtSolutionMaker::produce(), TtSemiLepSignalSelMVAComputer::produce(), RecoTauPiZeroFlattener::produce(), TtSemiLepJetCombMVAComputer::produce(), TtHadEvtSolutionMaker::produce(), L2TauPixelIsoTagProducer::produce(), TtFullHadKinFitProducer::produce(), QjetsAdder::produce(), reco::CorrectedJetProducer< T >::produce(), JVFJetIdProducer::produce(), HiFJRhoProducer::produce(), cms::JetVertexAssociation::produce(), cms::JetCorrectionProducer< T >::produce(), HLTMETCleanerUsingJetID::produce(), HLTMhtProducer::produce(), TtFullHadHypothesis::produce(), NoPileUpPFMEtProducer::produce(), HLTHtMhtProducer::produce(), pat::PATJetUpdater::produce(), TtSemiLepHypothesis::produce(), DeepBoostedJetTagInfoProducer::produce(), TtFullLepHypothesis::produce(), JetEnergyShift::produce(), PrimaryVertexSorter< ParticlesCollection >::produce(), pat::__class__Producer::produce(), HLTJetCollectionsForBoostedLeptonPlusJets< jetType >::produce(), RecoTauGenericJetRegionProducer< JetType, CandType >::produce(), HLTTrackMETProducer::produce(), HiGenCleaner< T2 >::produce(), GenJetFlavourTableProducer::produce(), TestBXVectorRefProducer::produce(), NoPileUpPFMEtDataProducer::produce(), NjettinessAdder::produce(), ShiftedPFCandidateProducerForNoPileUpPFMEt::produce(), ShiftedPFCandidateProducerForPFNoPUMEt::produce(), cms::METSignificanceProducer::produce(), GenJetTauTaggerProducer::produce(), TtJetPartonMatch< C >::produce(), l1t::FakeInputProducer::produce(), pat::PATJetProducer::produce(), DeepFlavourTagInfoProducer::produce(), PATObjectCrossLinker::produce(), DeepDoubleXTagInfoProducer::produce(), HiL1Subtractor::produce(), NanoAODBaseCrossCleaner::produce(), PFRecoTauChargedHadronProducer::produce(), FastPrimaryVertexWithWeightsProducer::produce(), PixelJetPuId::produce(), L1TStage2Layer2Producer::produce(), JetIdSelector< T >::produce(), l1t::GenToInputProducer::produce(), l1t::BXVectorInputProducer::produce(), JetConstituentSelector< T, C >::produce(), pat::JetCorrFactorsProducer::produce(), L1TStage1Layer2Producer::produce(), FastPrimaryVertexProducer::produce(), HLTJetsCleanedFromLeadingLeptons< JetType >::produce(), BoostedTopProducer::produce(), SoftLepton::produce(), PileupJetIdProducer::produce(), CaloJetMETcorrInputProducerT< T, Textractor >::produce(), JetCleanerForType1METT< T, Textractor >::produce(), BTagSFProducer::produce(), JetFlavourClustering::produce(), PFJetMETcorrInputProducerT< T, Textractor >::produce(), SmearedJetProducerT< T >::produce(), FastjetJetProducer::produceTrackJets(), QGTagger::putInEvent(), HFJetShowerShape::putInEvent(), HITRegionalPixelSeedGenerator::regions(), PrimaryVertexSorting::score(), HLTL1TSeed::seedsL1TriggerObjectMaps(), JetCombinatorics::Verbose(), GctFormatTranslateMCLegacy::writeGctOutJetBlock(), VirtualJetProducer::writeJets(), FFTJetProducer::writeJets(), Calculate::~Calculate(), Calculate_miniAOD::~Calculate_miniAOD(), ParticleLevelProducer::~ParticleLevelProducer(), and PrimaryVertexAssignment::~PrimaryVertexAssignment().

unique_ptr<JetDefinition::Plugin> fwrapper::plugin

Definition at line 43 of file fastjetfortran_madfks.cc.

Referenced by cms::CATopJetProducer::CATopJetProducer(), fastjetsiscone_(), fastjetsisconewitharea_(), GEDGsfElectronFinalizer::GEDGsfElectronFinalizer(), Vispa.Main.MainWindow.MainWindow::newAnalysisDesignerSlot(), Vispa.Main.MainWindow.MainWindow::newPxlSlot(), Vispa.Main.MainWindow.MainWindow::openAnalysisFileSlot(), and Vispa.Main.MainWindow.MainWindow::openPxlFileSlot().