EventSmear.h

Go to the documentation of this file.

#ifndef INCLUDE_EICSMEAR_SMEAR_EVENTSMEAR_H_
#define INCLUDE_EICSMEAR_SMEAR_EVENTSMEAR_H_

#include <cmath>
#include <list>
#include <vector>

#include <TObject.h>

#include "eicsmear/erhic/EventDis.h"
#include "eicsmear/erhic/Kinematics.h"
#include "eicsmear/smear/ParticleMCS.h"
#include "eicsmear/erhic/VirtualParticle.h"

namespace Smear {

/*
 A generator-independent DIS event with smeared kinematics and particles.
 */
class Event : public erhic::EventDis {
 public:
  Event();

  virtual ~Event();

  virtual void Reset();

  virtual void ClearParticles();

  virtual UInt_t GetNTracks() const;

  virtual const ParticleMCS* GetTrack(UInt_t) const;

  virtual ParticleMCS* GetTrack(UInt_t);

  virtual void SetQ2(double Q2) { QSquared = Q2; }

  virtual void SetX(double xB) { x = xB; }

  virtual void SetY(double inelasticity) { y = inelasticity; }

  virtual void SetW2(double W2) { WSquared = W2; }

  virtual void SetNu(double Nu) { nu = Nu; }

  virtual const ParticleMCS* BeamLepton() const;

  virtual const ParticleMCS* BeamHadron() const;

  virtual const ParticleMCS* ExchangeBoson() const;

  virtual const ParticleMCS* ScatteredLepton() const;

  virtual void AddLast(ParticleMCS* particle);

  void HadronicFinalState(ParticlePtrList&) const;

  std::vector<const erhic::VirtualParticle*> GetTracks() const;

  virtual void SetScattered(int index);

  virtual void Print(Option_t* = "") const;

 protected:
  Int_t nTracks;  
  std::vector<ParticleMCS*> particles;  
  Int_t mScatteredIndex;

  ClassDef(Smear::Event, 1)
};

inline UInt_t Event::GetNTracks() const {
  return particles.size();
}

inline const Smear::ParticleMCS* Event::GetTrack(UInt_t u) const {
  return (u < particles.size() ? particles.at(u) : NULL);
}

inline Smear::ParticleMCS* Event::GetTrack(UInt_t u) {
  return (u < particles.size() ? particles.at(u) : NULL);
}

inline const ParticleMCS* Event::BeamLepton() const {
  return (particles.empty() ? NULL : particles.front());
}

inline const ParticleMCS* Event::BeamHadron() const {
  return (particles.size() > 1 ? particles.at(1) : NULL);
}

inline const ParticleMCS* Event::ExchangeBoson() const {
  return NULL;
}

}  // namespace Smear

typedef Smear::Event EventS;

#endif  // INCLUDE_EICSMEAR_SMEAR_EVENTSMEAR_H_