TauolaEvent.cxx

00001 #include "TauolaEvent.h"
00002 #include "Plots.h"
00003 
00004 using namespace std;
00005 
00006 namespace Tauolapp
00007 {
00008 
00009 void TauolaEvent::undecayTaus(){
00010   
00011   std::vector<TauolaParticle*> particle_list;
00012   particle_list = findParticles(Tauola::getDecayingParticle());
00013 
00014   for(int p=0; p < (int) particle_list.size(); p++)
00015     particle_list.at(p)->findLastSelf()->undecay();
00016 
00017 }
00018 
00019 void TauolaEvent::decayTaus(){
00020 
00021   std::vector<TauolaParticle*> particle_list;
00022   particle_list = findStableParticles(Tauola::getDecayingParticle());
00023 
00024   while(particle_list.size()!=0){
00025 
00026     // tau and its matching tau-like partner is removed from the list here:
00027     TauolaParticlePair t_pair(particle_list);
00028 
00029     //t_pair.print();
00030     t_pair.decayTauPair();
00031     t_pair.checkMomentumConservation();
00032   }
00033 
00034   // Final event record modifications
00035   eventEndgame();
00036 }
00037 
00038 } // namespace Tauolapp