00001 // -*- C++ -*-
00002 #include "Rivet/Rivet.hh"
00003 #include "Rivet/Projections/MergedFinalState.hh"
00004 #include "Rivet/Tools/ParticleIdUtils.hh"
00005 #include "Rivet/Cmp.hh"
00006 #include <algorithm>
00007 
00008 namespace Rivet {
00009 
00010 
00011   int MergedFinalState::compare(const Projection& p) const {
00012     /// @todo: Currently A+B is not recognised to be the same as B+A.
00013     return mkNamedPCmp(p, "FSA") || mkNamedPCmp(p, "FSB");
00014   }
00015 
00016 
00017   void MergedFinalState::project(const Event& e) {
00018     const FinalState& fsa = applyProjection<FinalState>(e, "FSA");
00019     const FinalState& fsb = applyProjection<FinalState>(e, "FSB");
00020     _theParticles.clear();
00021     foreach (const Particle& pa, fsa.particles()){
00022       _theParticles.push_back(pa);
00023     }
00024     foreach (const Particle& pb, fsb.particles()){
00025       const GenParticle* originalb = &(pb.genParticle());
00026       bool notfound = true;
00027       foreach (const Particle& pa, fsa.particles()){
00028         const GenParticle* originala = &(pa.genParticle());
00029         if (originala == originalb) {
00030           notfound = false;
00031           break;
00032         }
00033       }
00034       if (notfound) {
00035         _theParticles.push_back(pb);
00036       }
00037     }
00038     getLog() << Log::DEBUG << "Number of particles in the two final states to be merged: = \n"
00039              << "   1st final state = " << fsa.particles().size() << endl
00040              << "   2nd final state = " << fsb.particles().size() << endl;
00041     getLog() << Log::DEBUG << "Number of merged final-state particles = "
00042              << _theParticles.size() << endl;
00043   }
00044 
00045 
00046 }
MergedFinalState.cc
