TauolaHepMC3Event.cxx

#include "TauolaHepMC3Event.h"
#include "Log.h"



namespace Tauolapp
{
using namespace std;
TauolaHepMC3Event::TauolaHepMC3Event(GenEvent * event){
  m_event=event;
/* NO UNITS YET
  // Default units
  m_momentum_unit = "GEV";
  m_length_unit   = "MM";

  if(m_event->momentum_unit() != Units::GEV) m_momentum_unit = "MEV";
  if(m_event->length_unit()   != Units::MM ) m_length_unit   = "CM";

  // If needed - change units used by HepMC to GEV and MM
  if( m_event->momentum_unit() != Units::GEV ||
      m_event->length_unit()   != Units::MM     )
  {
    m_event->use_units(Units::GEV,Units::MM);
  }
*/
}

TauolaHepMC3Event::~TauolaHepMC3Event(){

  while(m_tau_list.size()!=0){
    TauolaParticle * temp = m_tau_list.back();
    m_tau_list.pop_back();
    delete temp;
  }

}

GenEvent * TauolaHepMC3Event::getEvent(){
  return m_event;
}

std::vector<TauolaParticle*> TauolaHepMC3Event::findParticles(int pdg_id){

  if(m_tau_list.size()==0){

    //loop over all particle in the event looking for taus (or other)
    for( unsigned int i=0; i<m_event->particles().size(); ++i) {
      if(abs((m_event->particles()[i])->pid())==pdg_id)
        m_tau_list.push_back(new TauolaHepMC3Particle(m_event->particles()[i]));
    }
  }
  return m_tau_list;
}

std::vector<TauolaParticle*> TauolaHepMC3Event::findStableParticles(int pdg_id){

  /**  GenEvent::particle_const_iterator part_itr = m_event->particles_begin();
  //loop over all particle in the event looking for taus (or other)
  for( ; part_itr!=m_event->particles_end(); part_itr++){
    if(fabs((*part_itr)->pdg_id())==pdg_id){
      if((*part_itr)->end_vertex()){
        cout << "WARNING: Particle with pdg code " << (*part_itr)->pdg_id()
             << " has end vertex" <<endl;
      }
      else
        list.push_back(new TauolaHepMC3Particle(*part_itr));
    }
    }**/

  std::vector<TauolaParticle*> tau_list = findParticles(pdg_id);
  std::vector<TauolaParticle*> stable_tau_list;

  for(int i=0; i<(int) tau_list.size(); i++){

    if(!tau_list.at(i)->hasDaughters())
      stable_tau_list.push_back(tau_list.at(i));
    else
    {
      std::vector<TauolaParticle*> t = tau_list.at(i)->getDaughters();
      //Ignore taus that we won't be decaying anyway
      if(t.size()==1) continue;
      if(t.size()==2 && (abs(t[0]->getPdgID())==15 || abs(t[1]->getPdgID())==15) ) continue;
      Log::Warning()<<"Particle with pdg code "<<tau_list.at(i)->getPdgID()
                    <<" already has daughters" <<endl;
    }
  }

  return stable_tau_list;

}

void TauolaHepMC3Event::eventEndgame(){
/* NO UNITS YET
  //Set output units for the event
  string momentum("GEV"),length("MM");

  switch(Tauola::momentumUnit)
  {
    case Tauola::GEV:
      momentum = "GEV";
      break;
    case Tauola::MEV:
      momentum = "MEV";
      break;
    default:
      momentum = m_momentum_unit;
  }

  switch(Tauola::lengthUnit)
  {
    case Tauola::MM:
      length = "MM";
      break;
    case Tauola::CM:
      length = "CM";
      break;
    default:
      length = m_length_unit;
  }

  m_event->use_units(momentum,length);
*/
}

} // namespace Tauolapp