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