本文整理汇总了C++中JetCollection::end方法的典型用法代码示例。如果您正苦于以下问题：C++ JetCollection::end方法的具体用法？C++ JetCollection::end怎么用？C++ JetCollection::end使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类JetCollection的用法示例。
在下文中一共展示了JetCollection::end方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: analyseAndReturn
BAT::TtbarHypothesis HitFitAnalyser::analyseAndReturn(const EventPtr event, const JetCollection jets, const JetCollection bjets, const LeptonPointer selectedLepton ) {
weight_ = event->weight() * prescale_ * scale_;
treeMan_->setCurrentFolder(histogramFolder_);
treeMan_->Fill("EventWeight", weight_ );
// const JetCollection jets(event->getCleanedJets( SelectionCriteria::ElectronPlusJetsReference ));
// const JetCollection bJets(event->getCleanedBJets( SelectionCriteria::ElectronPlusJetsReference ));
// const LeptonPointer selectedLepton = event->getSignalLepton( SelectionCriteria::ElectronPlusJetsReference );
string metPrefix = METAlgorithm::names.at(0);
const METPointer met(event->MET((METAlgorithm::value) 0));
// Get cleaned jets that aren't b tagged
JetCollection leadingLightJets;
JetCollection leadingBJets;
unsigned int maxNJet = std::min(5, int(jets.size()));
for ( unsigned int jetIndex=0; jetIndex < maxNJet; ++jetIndex ) {
bool isBJet = false;
JetPointer thisJet = jets[jetIndex];
for ( unsigned int bJetIndex=0; bJetIndex < bjets.size(); ++bJetIndex ) {
JetPointer thisBJet = bjets[bJetIndex];
if ( thisJet == thisBJet ) {
isBJet = true;
break;
}
}
if ( !isBJet ) leadingLightJets.push_back( thisJet );
else leadingBJets.push_back( thisJet );
}
//set MC matching flag
if (event->getDataType() == DataType::TTJets_amcatnloFXFX)
do_MC_matching = true;
else
do_MC_matching = false;
//prepare the jets collection
// Copy jets into an array
JetCollection jetCopy;
for (JetCollection::const_iterator j = leadingLightJets.begin(); j != leadingLightJets.end(); ++j) {
jetCopy.push_back(*j);
}
JetCollection bJetCopy;
for (JetCollection::const_iterator j = leadingBJets.begin(); j != leadingBJets.end(); ++j) {
bJetCopy.push_back(*j);
}
std::sort(jetCopy.begin(), jetCopy.end(), jetPtComp);
std::sort(bJetCopy.begin(), bJetCopy.end(), jetPtComp);
jetsForFitting.clear();
unsigned numJetsToFit = jetCopy.size();
if (jetCopy.size() >= 2) {
if (numJetsToFit > 5)
numJetsToFit = 5;
jetsForFitting.insert(jetsForFitting.begin(), jetCopy.begin(), jetCopy.begin() + numJetsToFit);
}
bJetsForFitting.clear();
numJetsToFit = bJetCopy.size();
if (bJetCopy.size() >= 2) {
if (numJetsToFit > 5)
numJetsToFit = 5;
bJetsForFitting.insert(bJetsForFitting.begin(), bJetCopy.begin(), bJetCopy.begin() + numJetsToFit);
}
BatHitFit hhFitter(electronTranslator_, muonTranslator_, jetTranslator_, metTranslator_, hitfitDefault_,
hitfitLepWMass_, hitfitHadWMass_, hitfitTopMass_);
// Clear the internal state
hhFitter.clear();
// Add lepton into hitfit
if ( isElectronChannel_ ) {
const ElectronPointer signalElectron(boost::static_pointer_cast<Electron>(selectedLepton));
//TODO: fix the fitter to accept lepton class OR particle!!
hhFitter.AddLepton(*signalElectron);
}
else {
const MuonPointer signalMuon(boost::static_pointer_cast<Muon>(selectedLepton));
//TODO: fix the fitter to accept lepton class OR particle!!
hhFitter.AddLepton(*signalMuon);
}
// Add jets into HitFit
// Also check if jets matched to ttbar partons are in the jets passed to the fit
bool quarkInCollection = false, quarkBarInCollection = false, lebBInCollection = false, hadBInCollection = false;
for (size_t jet = 0; jet != jetsForFitting.size(); ++jet) {
hhFitter.AddJet(*jetsForFitting.at(jet));
// cout << "Adding light jet with csv : " << jetsForFitting.at(jet)->getBTagDiscriminator(BAT::BtagAlgorithm::value::CombinedSecondaryVertexV2) << endl;
if ( jetsForFitting.at(jet)->ttbar_decay_parton() ) {
int partonPdg = jetsForFitting.at(jet)->ttbar_decay_parton();
if ( partonPdg == 3 ) quarkInCollection = true;
else if ( partonPdg == 4 ) quarkBarInCollection = true;
else if ( partonPdg == 5 ) lebBInCollection = true;
else if ( partonPdg == 6 ) hadBInCollection = true;
}
//.........这里部分代码省略.........