smear/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_