d3/d6f/EicRcEvent_8h_source.html

//

// AYK (ayk@bnl.gov), 2014/07/09

//

//  EIC code reconstructed event class; consider to reuse eic-smear codes;

//

//  A good fraction of this code is cut'n'paste from EventSmear.h (eic-smear

//  package);


#include <eicsmear/erhic/EventMC.h>


#include <EicRcParticle.h>


#ifndef _EIC_RC_EVENT_

#define _EIC_RC_EVENT_


class EicRcEvent : public erhic::EventDis {

  friend class EicEventAssembler;


 public:


 EicRcEvent(): mScatteredLeptonID(-1)/*, femc(0.0), fhac(0.0)*/ {};


  virtual ~EicRcEvent() { ClearParticles(); };


#if _LATER_


  virtual void Reset();

#endif


  virtual void ClearParticles();


  virtual UInt_t GetNTracks() const { return mParticles.size(); };


  virtual const EicRcParticle* GetTrack(UInt_t u) const {

    return (u < mParticles.size() ? mParticles.at(u) : NULL);

  };

  virtual       EicRcParticle* GetTrack(UInt_t u)       {

    return (u < mParticles.size() ? mParticles.at(u) : NULL);

  };


#if _LATER_

  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; }

#endif


  // Basically wrappers for EicRootManager calls; yet it is convenient to have

  // them here to easily access original generator event(track) classes (for

  // instance to see which of the simulated particles have no reconstructed

  // counterparts -> so were not detected/reconstructed);

  const erhic::EventMC*    GetGenMcEvent()                            const;

  // NB: numbering scheme here and in the below call follows ERHIC convention

  // rather than PANDA one (so indices 0,1 correspond to beam particles);

  const erhic::ParticleMC* GetGenMcTrack(int genMcIndex)              const;

  // Find this MC track in array of reconstructed particles;

  EicRcParticle* GetRcParticleMatchingGenMcTrack(unsigned genMcIndex) const;


  // The only source of information for these two guys is original EventMC;

  // assume that beam lepton is #0 and beam hadron is #1 in respective

  // EventMC internal array; 2016/11/15: use erhic::EventMC-inherited class methods instead;

  //const erhic::ParticleMC* BeamLepton()    const { return GetGenMcTrack(0); };

  const erhic::ParticleMC* BeamLepton()    const { return GetGenMcEvent()->BeamLepton(); };

  //const erhic::ParticleMC* BeamHadron()    const { return GetGenMcTrack(1); };

  const erhic::ParticleMC* BeamHadron()    const { return GetGenMcEvent()->BeamHadron(); };

  const erhic::ParticleMC* ExchangeBoson() const {return NULL; };


  virtual const EicRcParticle* ScatteredLepton() const {

    // Just return the stored pointer; its assignment is a different story;

    return mScatteredLeptonID == -1 ? 0 : mParticles[mScatteredLeptonID];

  };


  void HadronicFinalState(ParticlePtrList&) const;


#if _LATER_


  virtual void AddLast(ParticleMCS*);


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


  virtual void SetScattered(int);


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

#endif


  private:

  Int_t mScatteredLeptonID;               // scattered lepton ID in mParticles[]

  std::vector<EicRcParticle*> mParticles;


 public:

  // Clearly a hack for now;

  //double femc, fhac;


  ClassDef(EicRcEvent, 7);

};


#endif