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