PSimHit.h

Go to the documentation of this file.

#ifndef PSimHit_H
#define PSimHit_H

#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
#include "SimDataFormats/EncodedEventId/interface/EncodedEventId.h"

class TrackingSlaveSD;  // for friend declaration only

class PSimHit {
public:
  static constexpr unsigned int k_tidOffset = 200000000;

  PSimHit() : theDetUnitId(0) {}

  PSimHit(const Local3DPoint& entry,
          const Local3DPoint& exit,
          float pabs,
          float tof,
          float eloss,
          int particleType,
          unsigned int detId,
          unsigned int trackId,
          float theta,
          float phi,
          unsigned short processType = 0)
      : theEntryPoint(entry),
        theSegment(exit - entry),
        thePabs(pabs),
        theEnergyLoss(eloss),
        theThetaAtEntry(theta),
        thePhiAtEntry(phi),
        theTof(tof),
        theParticleType(particleType),
        theProcessType(processType),
        theDetUnitId(detId),
        theTrackId(trackId) {}

  Local3DPoint entryPoint() const { return theEntryPoint; }

  Local3DPoint exitPoint() const { return theEntryPoint + theSegment; }

  Local3DPoint localPosition() const { return theEntryPoint + 0.5 * theSegment; }

  LocalVector momentumAtEntry() const { return LocalVector(thetaAtEntry(), phiAtEntry(), pabs()); }

  LocalVector localDirection() const { return LocalVector(thetaAtEntry(), phiAtEntry(), 1.f); }

  Geom::Theta<float> thetaAtEntry() const { return Geom::Theta<float>(theThetaAtEntry); }

  Geom::Phi<float> phiAtEntry() const { return Geom::Phi<float>(thePhiAtEntry); }

  float pabs() const { return thePabs; }

  float timeOfFlight() const { return tof(); }

  float tof() const { return theTof; }

  float energyLoss() const { return theEnergyLoss; }

  int particleType() const { return theParticleType; }

  unsigned int detUnitId() const { return theDetUnitId; }

  unsigned int trackId() const { return theTrackId; }

  unsigned int originalTrackId() const { return theTrackId % k_tidOffset; }

  unsigned int offsetTrackId() const { return theTrackId / k_tidOffset; }

  static unsigned int addTrackIdOffset(unsigned int tId, unsigned int offset) { return offset * k_tidOffset + tId; }

  EncodedEventId eventId() const { return theEventId; }

  void setEventId(EncodedEventId e) { theEventId = e; }

  unsigned short processType() const { return theProcessType; }

  void setTof(float tof) { theTof = tof; }

protected:
  // properties
  Local3DPoint theEntryPoint;  // position at entry
  Local3DVector theSegment;    // exitPos - entryPos
  float thePabs;               // momentum
  float theEnergyLoss;         // Energy loss
  float theThetaAtEntry;
  float thePhiAtEntry;

  float theTof;  // Time Of Flight
  int theParticleType;
  unsigned short theProcessType;  // ID of the process which created the track
                                  // which created the PSimHit

  // association
  unsigned int theDetUnitId;
  unsigned int theTrackId;
  EncodedEventId theEventId;

  friend class TrackingSlaveSD;
};

std::ostream& operator<<(std::ostream& o, const PSimHit& hit);

#endif  // PSimHit_H