DigiScintHit.cxx

Go to the documentation of this file.

// $Id: DigiScintHit.cxx,v 1.9 2009/05/19 15:05:12 kasahara Exp $
//
// Methods for simulated scintillator hits
//
// messier@indiana.edu
#include "Digitization/DigiScintHit.h"
#include "Conventions/Munits.h"

#include <TDatabasePDG.h>
#include <TString.h>

#include <iostream>
ClassImp(DigiScintHit)

//......................................................................

DigiScintHit::DigiScintHit() :
  fStripId(),
  fT1(0),
  fX1(0),
  fY1(0),
  fZ1(0),
  fT2(0),
  fX2(0),
  fY2(0),
  fZ2(0),
  fDS(0),
  fDE(0),
  fParticleEnergy(0),
  fParticlePX1(0),
  fParticlePY1(0),
  fParticlePZ1(0),
  fParticleId(0),
  fTrackId(0),
  fFailBits(0)
{ }

//......................................................................

DigiScintHit::DigiScintHit(const DigiScintHit& in) :
  TObject(in),
  fStripId(in.StripEndId()),
  fT1(in.T1()),
  fX1(in.X1()),
  fY1(in.Y1()),
  fZ1(in.Z1()),
  fT2(in.T2()),
  fX2(in.X2()),
  fY2(in.Y2()),
  fZ2(in.Z2()),
  fDS(in.DS()),
  fDE(in.DE()),
  fParticleEnergy(in.ParticleEnergy()),
  fParticlePX1(in.ParticlePX1()),
  fParticlePY1(in.ParticlePY1()),
  fParticlePZ1(in.ParticlePZ1()),
  fParticleId(in.ParticleId()),
  fTrackId(in.TrackId()),
  fFailBits(in.FailBits())
{
}

//......................................................................

DigiScintHit::DigiScintHit(Int_t trkId, Int_t trkPid, Float_t trkEnergy,
                           PlexStripEndId sid,
                           Float_t t1, Float_t x1, Float_t y1, Float_t z1,
                           Float_t t2, Float_t x2, Float_t y2, Float_t z2,
                           Float_t ds, Float_t de, 
                           Float_t trkPx1, Float_t trkPy1, Float_t trkPz1) :
  fStripId(sid),
  fT1(t1),
  fX1(x1),
  fY1(y1),
  fZ1(z1),
  fT2(t2),
  fX2(x2),
  fY2(y2),
  fZ2(z2),
  fDS(ds),
  fDE(de),
  fParticleEnergy(trkEnergy),
  fParticlePX1(trkPx1),
  fParticlePY1(trkPy1),
  fParticlePZ1(trkPz1),
  fParticleId(trkPid),
  fTrackId(trkId),
  fFailBits(0)
{ }

//......................................................................

Float_t DigiScintHit::ParticleKineticEnergy() const
{
  TParticlePDG* part = TDatabasePDG::Instance()->GetParticle(ParticleId());
  Float_t mass = 0;
  if(part) 
    mass = part->Mass()*Munits::GeV;
  return ParticleEnergy()-mass;
}


//......................................................................

ostream& operator<<(ostream& os, const DigiScintHit& h) 
{ return h.FormatToOStream(os); }

//_____________________________________________________________________________
void DigiScintHit::Print(Option_t *option) const
{
   std::cout << this->AsString(option) << std::endl;
}

//_____________________________________________________________________________
std::ostream& DigiScintHit::FormatToOStream(std::ostream& os,
                                            Option_t *option) const
{
   os << AsString(option);
   return os;
}

//_____________________________________________________________________________
const char* DigiScintHit::AsString(Option_t *option) const
{
  TParticlePDG* part = TDatabasePDG::Instance()->GetParticle(ParticleId());

  switch(option[0]) {
  case 'e': // (E)xtended printout.
  case 'E':
    return Form("DigiScintHit Track=%d  Particle=%s(%d) E=%.1f GeV K=%f GeV\n" 
                "             Strip: %s\n"
                "             (%5.2f, %4.3f, %4.3f), %.1f ns\n"
                "             dE/dx = %.2f MeV / %.1f cm\n",
                TrackId(),
                (part==0)?"unknown":part->GetName(),
                ParticleId(),
                ParticleEnergy()/Munits::GeV,
                ParticleKineticEnergy()/Munits::GeV,
                StripEndId().AsString(),
                X1(),Y1(),Z1(),T1()/Munits::ns,
                DE()/Munits::MeV,
                DS()/Munits::cm           
                );
  case 'o': // old
    return Form("DigiScintHit (Trk %d) (%s) (%.1f GeV) (%.1f MeV/%.1f cm) (Plane %d Strip %d)",
                TrackId(),
                (part==0)?"unknown":part->GetName(),
                ParticleEnergy()/Munits::GeV,
                DE()/Munits::MeV,
                DS()/Munits::cm,
                Plane(),
                Strip() 
                );                
  default:
    return Form("DigiScintHit %s (Trk %4d=%6s %.2f GeV) (%.2f MeV/%.2f cm) %f",
                StripEndId().AsString("c"),
                TrackId(),
                (part==0)?"unknown":part->GetName(),
                ParticleEnergy()/Munits::GeV,
                DE()/Munits::MeV,
                DS()/Munits::cm,
                T1()/Munits::ns
                );                
  }
  return ""; // Keep compiler happy.
}

---