FP420G4Hit.cc

Go to the documentation of this file.

// File: FP420G4Hit.cc
// Date: 02.2006
// Description: Transient Hit class for the FP420
#include "SimG4CMS/FP420/interface/FP420G4Hit.h"
#include <iostream>

FP420G4Hit::FP420G4Hit():entry(0) {

  entrylp(0);
  exitlp(0);
  elem     = 0.;
  hadr     = 0.;
  theIncidentEnergy = 0.;
  theTimeSlice = 0.;
  theTrackID = -1;
  theUnitID  =  0;
  thePabs =0.;
  theTof=0. ;
  theEnergyLoss=0.   ;
  theParticleType=0 ;
  theUnitID=0;
  theTrackID=-1;
  theThetaAtEntry=-10000. ;
  thePhiAtEntry=-10000. ;
  theParentId=0;
  
  theX = 0.;
  theY = 0.;
  theZ = 0.;
  theVx = 0.;
  theVy = 0.;
  theVz = 0.;
}


FP420G4Hit::~FP420G4Hit(){}


FP420G4Hit::FP420G4Hit(const FP420G4Hit &right) {
  theUnitID         = right.theUnitID;
  
  theTrackID        = right.theTrackID;
  theTof            = right.theTof ;
  theEnergyLoss      = right.theEnergyLoss   ;
  theParticleType    = right.theParticleType ;
  thePabs            = right.thePabs;
  elem               = right.elem;
  hadr               = right.hadr;
  theIncidentEnergy  = right.theIncidentEnergy;
  theTimeSlice       = right.theTimeSlice;
  entry              = right.entry;
  entrylp            = right.entrylp;
  exitlp             = right.exitlp;
  theThetaAtEntry    = right.theThetaAtEntry;
  thePhiAtEntry      = right.thePhiAtEntry;
  theParentId        = right.theParentId;
  
  theX = right.theX;
  theY = right.theY;
  theZ = right.theZ;
  
  theVx = right.theVx;
  theVy = right.theVy;
  theVz = right.theVz;
  
  
}


const FP420G4Hit& FP420G4Hit::operator=(const FP420G4Hit &right) {
  theUnitID         = right.theUnitID;
  
  theTrackID        = right.theTrackID;
  theTof            = right.theTof ;
  theEnergyLoss      = right.theEnergyLoss   ;
  theParticleType    = right.theParticleType ;
  thePabs            = right.thePabs;
  elem               = right.elem;
  hadr               = right.hadr;
  theIncidentEnergy  = right.theIncidentEnergy;
  theTimeSlice       = right.theTimeSlice;
  entry              = right.entry;
  entrylp            = right.entrylp;
  exitlp             = right.exitlp;
  theThetaAtEntry    = right.theThetaAtEntry;
  thePhiAtEntry      = right.thePhiAtEntry;
  theParentId        = right.theParentId;
  
  theX = right.theX;
  theY = right.theY;
  theZ = right.theZ;
  
  theVx = right.theVx;
  theVy = right.theVy;
  theVz = right.theVz;
  
  
  return *this;
}


void FP420G4Hit::addEnergyDeposit(const FP420G4Hit& aHit) {

  elem += aHit.getEM();
  hadr += aHit.getHadr();
}


void FP420G4Hit::Print() {
  std::cout << (*this);
}

G4ThreeVector   FP420G4Hit::getEntry() const           {return entry;}
void         FP420G4Hit::setEntry(const G4ThreeVector& xyz)   { entry    = xyz; }

G4ThreeVector    FP420G4Hit::getEntryLocalP() const           {return entrylp;}
void         FP420G4Hit::setEntryLocalP(const G4ThreeVector& xyz1)   { entrylp    = xyz1; }

G4ThreeVector     FP420G4Hit::getExitLocalP() const           {return exitlp;}
void         FP420G4Hit::setExitLocalP(const G4ThreeVector& xyz1)   { exitlp    = xyz1; }

double       FP420G4Hit::getEM() const              {return elem; }
void         FP420G4Hit::setEM (double e)           { elem     = e; }
      
double       FP420G4Hit::getHadr() const            {return hadr; }
void         FP420G4Hit::setHadr (double e)         { hadr     = e; }
      
double       FP420G4Hit::getIncidentEnergy() const  {return theIncidentEnergy; }
void         FP420G4Hit::setIncidentEnergy (double e){theIncidentEnergy  = e; }

unsigned int  FP420G4Hit::getTrackID() const         {return theTrackID; }
void         FP420G4Hit::setTrackID (int i)         { theTrackID = i; }

unsigned int FP420G4Hit::getUnitID() const          {return theUnitID; }
void         FP420G4Hit::setUnitID (unsigned int i) { theUnitID = i; }

double       FP420G4Hit::getTimeSlice() const       {return theTimeSlice; }
void         FP420G4Hit::setTimeSlice (double d)    { theTimeSlice = d; }
int          FP420G4Hit::getTimeSliceID() const     {return (int)theTimeSlice;}

void         FP420G4Hit::addEnergyDeposit(double em, double hd)
                                                   {elem  += em ; hadr += hd;}

double       FP420G4Hit::getEnergyDeposit() const   {return elem+hadr;}

float FP420G4Hit::getPabs() const {return thePabs;}
float FP420G4Hit::getTof() const {return theTof;}
float FP420G4Hit::getEnergyLoss() const {return theEnergyLoss;}
int FP420G4Hit::getParticleType() const {return theParticleType;}

void FP420G4Hit::setPabs(float e) {thePabs = e;}
void FP420G4Hit::setTof(float e) {theTof = e;}
void FP420G4Hit::addEnergyLoss(float e) {theEnergyLoss += e;}
void FP420G4Hit::setEnergyLoss(float e) {theEnergyLoss = e;}
void FP420G4Hit::setParticleType(short i) {theParticleType = i;}

float FP420G4Hit::getThetaAtEntry() const {return theThetaAtEntry;}   
float FP420G4Hit::getPhiAtEntry() const{ return thePhiAtEntry;}

void FP420G4Hit::setThetaAtEntry(float t){theThetaAtEntry = t;}
void FP420G4Hit::setPhiAtEntry(float f) {thePhiAtEntry = f ;}

float FP420G4Hit::getX() const{ return theX;}
void FP420G4Hit::setX(float t){theX = t;}

float FP420G4Hit::getY() const{ return theY;}
void FP420G4Hit::setY(float t){theY = t;}

float FP420G4Hit::getZ() const{ return theZ;}
void FP420G4Hit::setZ(float t){theZ = t;}

int FP420G4Hit::getParentId() const {return theParentId;}
void FP420G4Hit::setParentId(int p){theParentId = p;}

float FP420G4Hit::getVx() const{ return theVx;}
void FP420G4Hit::setVx(float t){theVx = t;}

float FP420G4Hit::getVy() const{ return theVy;}
void FP420G4Hit::setVy(float t){theVy = t;}

float FP420G4Hit::getVz() const{ return theVz;}
void FP420G4Hit::setVz(float t){theVz = t;}





std::ostream& operator<<(std::ostream& os, const FP420G4Hit& hit) {
  os << " Data of this FP420G4Hit are:" << std::endl
     << " hitEntryLocalP: " << hit.getEntryLocalP() << std::endl
     << " hitExitLocalP: " << hit.getExitLocalP() << std::endl
     << " Time slice ID: " << hit.getTimeSliceID() << std::endl
     << " Time slice : " << hit.getTimeSlice() << std::endl
     << " Tof : " << hit.getTof() << std::endl
     << " EnergyDeposit = " << hit.getEnergyDeposit() << std::endl
     << " elmenergy = " << hit.getEM() << std::endl
     << " hadrenergy = " << hit.getHadr() << std::endl
     << " EnergyLoss = " << hit.getEnergyLoss() << std::endl
     << " ParticleType = " << hit.getParticleType() << std::endl
     << " ParticleType = " << hit.getParticleType() << std::endl
     << " Theta at Entry = " << hit.getThetaAtEntry() << std::endl
     << " Phi at Entry = " << hit.getPhiAtEntry() << std::endl
     << " X at Entry = " << hit.getX() << std::endl
     << " Y at Entry = " << hit.getY() << std::endl
     << " Z at Entry = " << hit.getZ() << std::endl
     << " VtxX = " << hit.getVx() << std::endl
     << " VtxY = " << hit.getVy() << std::endl
     << " VtxZ = " << hit.getVz() << std::endl
     << " parentID = " << hit.getParentId() << std::endl
     << " Pabs = " << hit.getPabs() << std::endl
     << " Energy of primary particle (ID = " << hit.getTrackID()
     << ") = " << hit.getIncidentEnergy() << " (MeV)"<<std::endl
     << " Entry point in FP420 unit number " << hit.getUnitID()
     << " is: " << hit.getEntry() << " (mm)" << std::endl;
  os << "~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"
     << std::endl;
  return os;
}