Kinematics.h

Go to the documentation of this file.

#ifndef INCLUDE_EICSMEAR_ERHIC_KINEMATICS_H_
#define INCLUDE_EICSMEAR_ERHIC_KINEMATICS_H_

#include <list>
#include <vector>

#include <Rtypes.h>
#include <TLorentzVector.h>

#include "eicsmear/erhic/ParticleMC.h"
#include "eicsmear/erhic/ParticleIdentifier.h"
#include "eicsmear/erhic/VirtualParticle.h"

namespace erhic {

class EventDis;
class VirtualParticle;

struct DisKinematics : public TObject {
  DisKinematics();
  DisKinematics(double x, double y, double nu, double Q2, double W2);
  Double32_t mX;
  Double32_t mQ2;
  Double32_t mW2;
  Double32_t mNu;
  Double32_t mY;

  ClassDef(erhic::DisKinematics, 1)
};

class KinematicsComputer {
 public:
  virtual ~KinematicsComputer() { }
  virtual TObject* Calculate() = 0;
  ClassDef(erhic::KinematicsComputer, 1)
};

class LeptonKinematicsComputer : public KinematicsComputer {
 public:
  virtual ~LeptonKinematicsComputer() { }
  explicit LeptonKinematicsComputer(const EventDis&);
  virtual DisKinematics* Calculate();

 protected:
  std::vector<const VirtualParticle*> mBeams;

  ClassDef(erhic::LeptonKinematicsComputer, 1)
};

class JacquetBlondelComputer : public KinematicsComputer {
 public:
  virtual ~JacquetBlondelComputer();
  explicit JacquetBlondelComputer(const EventDis&);
  virtual DisKinematics* Calculate();

 protected:
  virtual Double_t ComputeY() const;
  virtual Double_t ComputeQSquared() const;
  virtual Double_t ComputeX() const;
  const EventDis& mEvent;
  std::vector<const VirtualParticle*> mParticles;

  ClassDef(erhic::JacquetBlondelComputer, 1)
};

class DoubleAngleComputer : public KinematicsComputer {
 public:
  virtual ~DoubleAngleComputer();
  explicit DoubleAngleComputer(const EventDis&);
  virtual DisKinematics* Calculate();

 protected:
  const EventDis& mEvent;
  virtual Double_t ComputeQuarkAngle() const;
  virtual Double_t ComputeY() const;
  virtual Double_t ComputeQSquared() const;
  virtual Double_t ComputeX() const;
  mutable Bool_t mHasChanged;
  mutable Double_t mAngle;
  std::vector<const VirtualParticle*> mParticles;

  ClassDef(erhic::DoubleAngleComputer, 1)
};

}  // namespace erhic

#endif  // INCLUDE_EICSMEAR_ERHIC_KINEMATICS_H_