PTSimHit.cxx

Go to the documentation of this file.

//
// PTSimHit
//
// PTSimHit stores data from a single hit recorded during transport
// of a track through an active detector element.
//
// S. Kasahara  05/04 

#include <iostream>
using namespace std;

#include <TClonesArray.h>

#include "Digitization/DigiScintHit.h"

#include "ParticleTransportSim/PTSimHit.h"
#include "Conventions/Munits.h"
#include "MessageService/MsgService.h"

ClassImp(PTSimHit)

std::ostream& operator << (std::ostream& os, const PTSimHit& mh)
                                          { return mh.Print(os); }

CVSID("$Id: PTSimHit.cxx,v 1.10 2009/04/24 19:10:26 schubert Exp $");


//_________________________________________________________
PTSimHit::PTSimHit() : fUseLinearStep(false) {
   //  Default constructor
  this -> Clear(); // initialize data members
}

//__________________________________________________________
PTSimHit::~PTSimHit() {
// Destructor

}

//__________________________________________________________
void PTSimHit::Clear(Option_t* /* option */) {
// Resets hit data so that object can be reused

  fStripEndId.SetDetector(Detector::kFar);
  fStripEndId.SetPlane(999);
  fStripEndId.SetStrip(999);
  fPdg = 0;
  fTrkIndex = -1;
  fGlobalPos0.SetXYZT(0,0,0,0);
  fLocalPos0.SetXYZT(0,0,0,0);
  fP40.SetPxPyPzE(0,0,0,0);
  fGlobalPos.SetXYZT(0,0,0,0);
  fLocalPos.SetXYZT(0,0,0,0);
  fP4.SetPxPyPzE(0,0,0,0);
  fELoss = 0;
  fStep = 0;

}  

//__________________________________________________________
void PTSimHit::CreateDigiScintHit(TClonesArray& hitarray) const {
// Creates DigiScintHit from stored data and adds it to users
// TClonesArray at position entry

  MSG("PTSim",Msg::kVerbose) << "PTSimHit::CreateDigiScintHit\n" 
                               << *(this) << endl;

  if ( fStripEndId.GetPlane() == 999 ) {
    MSG("PTSim",Msg::kFatal) 
       << "CreateDigiScintHit called with invalid plane id " << endl;
    abort();
  }
  
  Int_t entry = hitarray.GetEntriesFast();
  Double_t step = fStep;
  if ( fUseLinearStep ) {
    // This is a bad idea but is useful for comparing to gminos
    double dx = fLocalPos.X() - fLocalPos0.X();
    double dy = fLocalPos.Y() - fLocalPos0.Y();
    double dz = fLocalPos.Z() - fLocalPos0.Z();
    step = sqrt(dx*dx+dy*dy+dz*dz);
  }
  
  new((hitarray)[entry]) DigiScintHit(fTrkIndex,fPdg,fP40.E(),
                           fStripEndId,
                           fLocalPos0.T(),
                           fLocalPos0.X()*Munits::cm,
                           fLocalPos0.Y()*Munits::cm,
                           fLocalPos0.Z()*Munits::cm,
                           fLocalPos.T(),
                           fLocalPos.X()*Munits::cm,
                           fLocalPos.Y()*Munits::cm,
                           fLocalPos.Z()*Munits::cm,
                           step*Munits::cm,fELoss,
                           fP40.Px(),fP40.Py(),fP40.Pz());
  return;
}  


//__________________________________________________________
void PTSimHit::Init(const PlexStripEndId& seid, Int_t pdg, Int_t trkind,
                      const TLorentzVector& globalpos,
                      const TLorentzVector& localpos, 
                      const TLorentzVector& p4) {
// Store entry point data for this hit
 
  fStripEndId = seid;
  fPdg = pdg;
  fTrkIndex = trkind;
  fGlobalPos0 = globalpos;
  fLocalPos0 = localpos;
  fP40 = p4;

  return;

}  

//__________________________________________________________
std::ostream& PTSimHit::Print(std::ostream& os) const {
  //  Print hit data on ostream.


  os << "PTSimHit:" << endl;
 
  os << "StripEndId " << fStripEndId.AsString() 
     << " particle id(pdg) " << fPdg 
     << " trk index " << fTrkIndex << endl;
  os << "Entry x,y,z (local)(cm): " << fLocalPos0.X() << " " << fLocalPos0.Y()
     << " " << fLocalPos0.Z() << endl;
  os << "      x,y,z(global)(cm): " << fGlobalPos0.X() << " "  
     << fGlobalPos0.Y() << " "
     << fGlobalPos0.Z() << endl;
  os << "      time(sec) " <<  fLocalPos0.T() << " p4(GeV) "
     << fP40.Px() << " " << fP40.Py() << " " 
     << fP40.Pz() << " " << fP40.E() << endl;
  os << "Exit x,y,z (local)(cm): " << fLocalPos.X() << " " << fLocalPos.Y()
     << " " << fLocalPos.Z() << endl;
  os << "     x,y,z(global)(cm): " << fGlobalPos.X() << " " 
     << fGlobalPos.Y() << " " 
     << fGlobalPos.Z() << endl;
  os << "     time(sec) " <<  fLocalPos.T() << " p4(GeV) "
     << fP4.Px() << " " << fP4.Py() << " " 
     << fP4.Pz() << " " << fP4.E() << endl;
  os << "Eloss(GeV) " << fELoss << " Step(cm) " << fStep << endl;

  return os;

}

//__________________________________________________________
void PTSimHit::Print(Option_t* /* option */) const {
  //  Print hit in form supported by TObject::Print.

  Print(std::cout);
  return;
}

//__________________________________________________________
void PTSimHit::SetCurrentState(const TLorentzVector& globalpos,
                                 const TLorentzVector& localpos, 
                                 const TLorentzVector& p4) {
// Store current data for this track

  fGlobalPos = globalpos;
  fLocalPos = localpos;
  fP4 = p4;

  return;

}  

---