00001 // -*- C++ -*-
00002 #include "Rivet/Projections/UnstableFinalState.hh"
00003 #include "Rivet/Tools/Logging.hh"
00004 #include "Rivet/Cmp.hh"
00005 
00006 #define IS_PARTON_PDGID(id) ( abs(id) <= 100 && abs(id) != 22 && (abs(id) < 11 || abs(id) > 18) )
00007 
00008 namespace Rivet {
00009 
00010 
00011   int UnstableFinalState::compare(const Projection& p) const {
00012     const UnstableFinalState& other = dynamic_cast<const UnstableFinalState&>(p);
00013     return \
00014       cmp(_etamin, other._etamin) ||
00015       cmp(_etamax, other._etamax) ||
00016       cmp(_ptmin, other._ptmin);
00017   }
00018 
00019 
00020   void UnstableFinalState::project(const Event& e) {
00021     Log& log = getLog();
00022     _theParticles.clear();
00023 
00024     for (GenEvent::particle_const_iterator p = e.genEvent().particles_begin();
00025          p != e.genEvent().particles_end(); ++p) {
00026       const int st = (*p)->status();
00027       const double pT = (*p)->momentum().perp();
00028       const double eta = (*p)->momentum().eta();
00029       const GenVertex* pv = (*p)->production_vertex();
00030       const GenVertex* dv = (*p)->end_vertex();
00031       bool passed = ( st == 1 || (st == 2 && abs((*p)->pdg_id()) != 22) ) &&
00032                     !isZero(pT) &&
00033                     pT >= _ptmin &&
00034                     eta > _etamin &&
00035                     eta < _etamax &&
00036                     !IS_PARTON_PDGID((*p)->pdg_id());
00037       if (passed) {
00038         if (pv!=NULL) {
00039           for (GenVertex::particles_in_const_iterator pp = pv->particles_in_const_begin() ;
00040               pp != pv->particles_in_const_end() ; ++pp) {
00041             // Avoid double counting if particle == parent
00042             if ( (*p)->pdg_id() == (*pp)->pdg_id() )
00043               passed = false;
00044           }
00045         }
00046       }
00047 
00048       if (log.isActive(Log::TRACE)) {
00049         log << Log::TRACE << std::boolalpha
00050             << "ID = " << (*p)->pdg_id() << ", status = " << st << ", pT = " << pT
00051             << ", eta = " << eta << ": result = " << passed << endl;
00052         if (pv!=NULL) {
00053           for (GenVertex::particles_in_const_iterator pp = pv->particles_in_const_begin() ;
00054               pp != pv->particles_in_const_end() ; ++pp) {
00055             log << Log::TRACE << std::boolalpha
00056                 << "     parent ID = " << (*pp)->pdg_id() << endl;
00057           }
00058         }
00059         if (dv!=NULL) {
00060           for (GenVertex::particles_out_const_iterator pp = dv->particles_out_const_begin() ;
00061               pp != dv->particles_out_const_end() ; ++pp) {
00062             log << Log::TRACE << std::boolalpha
00063                 << "     child ID  = " << (*pp)->pdg_id() << endl;
00064           }
00065         }
00066       }
00067       if (passed) _theParticles.push_back(Particle(**p));
00068     }
00069     log << Log::DEBUG << "Number of final-state particles = "
00070         << _theParticles.size() << endl;
00071   }
00072 
00073 }
UnstableFinalState.cc
