TauolaParticle.h

#ifndef _TauolaParticle_h_included_
#define _TauolaParticle_h_included_

#include <iostream>
#include <math.h>
#include <vector>

#include "DecayList.h"
#include "Tauola.h"
#include "f_Decay.h"

namespace Tauolapp
{

class TauolaParticle{

 public:

  virtual ~TauolaParticle(){};

  static const int SAME_SIGN=1;

  static const int OPPOSITE_SIGN=2;

  static const int NA_SIGN=3;

  static const int STABLE=1;

  static const int DECAYED=2;

  static const int HISTORY=3;  

  static const int X_AXIS=1;

  static const int Y_AXIS=2;

  static const int Z_AXIS=3;

  static const int Z0 = 23;

  static const int HIGGS = 25;

  static const int HIGGS_A = 36;

  static const int HIGGS_PLUS = 37;

  static const int HIGGS_MINUS = -37;


  static const int W_PLUS = 24;

  static const int W_MINUS = -24;

  static const int GAMMA = 22;

  static const int TAU_PLUS = -15;

  static const int TAU_MINUS = 15;

  static const int TAU_NEUTRINO = 16;

  static const int TAU_ANTINEUTRINO = -16;


  static const int MUON_PLUS = -13;

  static const int MUON_MINUS = 13;

  static const int MUON_NEUTRINO = 14;

  static const int MUON_ANTINEUTRINO = -14;


  static const int POSITRON = -11;

  static const int ELECTRON = 11;

  static const int ELECTRON_NEUTRINO = 12;

  static const int ELECTRON_ANTINEUTRINO = -12;

  static const int UP = 2;

  static const int ANTIUP = -2;

  static const int DOWN = 1;

  static const int ANTIDOWN = -1;

  static const int GLUON = 21;
  static const int CHARM = 4;
  static const int TOP = 6;
  static const int STRANGE = 3;
  static const int BOTTOM = 5;

  static const int OTHER = 0;

  TauolaParticle * clone();

  double getAngle(TauolaParticle *);

  void add(TauolaParticle *);

  void subtract(TauolaParticle *);

  void decay();

  void addDecayToEventRecord();

  int  getSign();

  double getPolarimetricX();

  double getPolarimetricY();

  double getPolarimetricZ(); 

  bool hasDaughters();

  TauolaParticle * findLastSelf();

  std::vector<TauolaParticle *> findProductionMothers();

  void boostToRestFrame(TauolaParticle * boost);

  void boostDaughtersToRestFrame(TauolaParticle * boost);


  void boostFromRestFrame(TauolaParticle * boost);

  void boostDaughtersFromRestFrame(TauolaParticle * boost);

  void boostAlongZ(double pz, double e);

  void rotate(int axis, double phi, int second_axis=Z_AXIS);

  void rotateDaughters(int axis, double phi, int second_axis=Z_AXIS);

  double getRotationAngle(int axis, int second_axis=Z_AXIS);

  double getP();

  double getP(int axis);

  void  setP(int axis, double p_component);

  double getMass();



  /********************************************** 
      Beginning of virtual methods 

  ********************************************/

  virtual void undecay(){};

  virtual void checkMomentumConservation(){};

  virtual void decayEndgame(){};

  virtual double getPx()=0;

  virtual double getPy()=0;

  virtual double getPz()=0;

  virtual double getE()=0;

  virtual void setPx( double px )=0;

  virtual void setPy( double py )=0;

  virtual void setPz( double pz )=0;

  virtual void setE( double e )=0;

  virtual void setMothers(std::vector<TauolaParticle*> mothers)=0;

  virtual void setDaughters(std::vector<TauolaParticle*> daughters)=0;

  virtual std::vector<TauolaParticle*> getMothers()=0;

  virtual std::vector<TauolaParticle*> getDaughters()=0;

  virtual void setPdgID(int pdg_id)=0;

  virtual void setMass(double mass)=0;

  virtual void setStatus(int status)=0;

  virtual int getPdgID()=0;

  virtual int getStatus()=0;

  virtual int getBarcode()=0;

  virtual TauolaParticle * createNewParticle(int pdg_id, int status, 
                                             double mass, double px,
                                             double py, double pz,
                                             double e)=0;

  virtual void print()=0;

 private:

  double m_pol_x;

  double m_pol_y;

  double m_pol_z;

  double m_pol_n;
};

} // namespace Tauolapp
#endif