#include <TRTProcessingOfEndCapHits.h>

Collaboration diagram for TRTProcessingOfEndCapHits:

Public Member Functions
	TRTProcessingOfEndCapHits (TRTSensitiveDetector *)
	~TRTProcessingOfEndCapHits ()
bool	ProcessHit (G4Step *)

Private Member Functions
	TRTProcessingOfEndCapHits (const TRTProcessingOfEndCapHits &)
TRTProcessingOfEndCapHits &	operator= (const TRTProcessingOfEndCapHits &)
void	Initialize ()
void	DeleteArrays ()

Private Attributes
int	m_printMessages
int	m_sectorsABC
int	m_testLocalCoordinatesOfHits
int	m_numberOfStrawsInPlanesAB
int	m_numberOfStrawsInPlaneC
int	m_numberOfWheelsAB
int	m_numberOfStrawsInSectorsAB
int	m_numberOfStrawsInSectorC
int	m_initialStrawIDInSectorsAB
int	m_initialStrawIDInSectorC
int	m_numberOfStrawsInIDSectorsAB
int	m_numberOfStrawsInIDSectorC
const TRTParameters *	m_pParameters
TRTSensitiveDetector *	m_pSensitiveDetector
int	m_verboseLevel

Friends
class	TRTParametersForEndCapHits

Detailed Description

Definition at line 15 of file TRTProcessingOfEndCapHits.h.

Constructor & Destructor Documentation

◆ TRTProcessingOfEndCapHits() [1/2]

TRTProcessingOfEndCapHits::TRTProcessingOfEndCapHits ( TRTSensitiveDetector * pSensitiveDet )

Definition at line 31 of file TRTProcessingOfEndCapHits.cxx.

                                     :
  m_printMessages(0), //FIXME obsolete?
  m_sectorsABC(0),
  m_testLocalCoordinatesOfHits(0),
  m_numberOfStrawsInPlanesAB(0),
  m_numberOfStrawsInPlaneC(0),
  m_numberOfWheelsAB(0),
  m_numberOfStrawsInSectorsAB(0),
  m_numberOfStrawsInSectorC(0),
  m_initialStrawIDInSectorsAB(0),
  m_initialStrawIDInSectorC(0),
  m_numberOfStrawsInIDSectorsAB(0),
  m_numberOfStrawsInIDSectorC(0),
  m_pParameters(nullptr),
  m_pSensitiveDetector(pSensitiveDet),
  m_verboseLevel(pSensitiveDet->verboseLevel)
{
  m_pParameters = TRTParameters::GetPointer();
 
  m_printMessages = m_pParameters->GetInteger("PrintMessages"); //FIXME Obsolete?
 
  if (m_verboseLevel>5) { G4cout << "##### Constructor TRTProcessingOfEndCapHits" << G4endl; }
 
  this->Initialize();
 
  m_pParameters = nullptr;
 
  if (m_verboseLevel>5) { G4cout << "##### Constructor TRTProcessingOfEndCapHits done" << G4endl; }
}

◆ ~TRTProcessingOfEndCapHits()

TRTProcessingOfEndCapHits::~TRTProcessingOfEndCapHits ( )

Definition at line 65 of file TRTProcessingOfEndCapHits.cxx.

66{

67}

◆ TRTProcessingOfEndCapHits() [2/2]

TRTProcessingOfEndCapHits::TRTProcessingOfEndCapHits ( const TRTProcessingOfEndCapHits & )

private

Member Function Documentation

◆ DeleteArrays()

void TRTProcessingOfEndCapHits::DeleteArrays ( )

private

◆ Initialize()

void TRTProcessingOfEndCapHits::Initialize ( )

private

Definition at line 72 of file TRTProcessingOfEndCapHits.cxx.

{
  if (m_verboseLevel>5) { G4cout << "######### Method TRTProcessingOfEndCapHits::Initialize" << G4endl; }
 
  m_sectorsABC = m_pParameters->GetInteger("SectorsABC");
 
  m_testLocalCoordinatesOfHits =
    m_pParameters->GetInteger("TestLocalCoordinatesOfHits");
 
  if (m_verboseLevel>5) { G4cout << "######### Method TRTProcessingOfEndCapHits::Initialize done" << G4endl; }
 
}

◆ operator=()

TRTProcessingOfEndCapHits & TRTProcessingOfEndCapHits::operator= ( const TRTProcessingOfEndCapHits & )

private

◆ ProcessHit()

bool TRTProcessingOfEndCapHits::ProcessHit ( G4Step * pStep )

Definition at line 88 of file TRTProcessingOfEndCapHits.cxx.

{
  G4Track* pTrack = pStep->GetTrack();
  // get the HepMC barcode using the track helper
  TrackHelper trHelp(pTrack);
 
  G4StepPoint* pPreStepPoint = pStep->GetPreStepPoint();
  G4StepPoint* pPostStepPoint = pStep->GetPostStepPoint();
 
  G4TouchableHandle pTouchableHandle = pPreStepPoint->GetTouchableHandle();
  int depth = 1;
 
  int strawID = pTouchableHandle->GetReplicaNumber(depth++);
 
  int sectorID = 0;
 
  if (m_sectorsABC)
    sectorID = pTouchableHandle->GetReplicaNumber(depth++);
 
  int planeID = pTouchableHandle->GetReplicaNumber(depth++);
 
  int wheelID = pTouchableHandle->GetReplicaNumber(depth);
 
  G4ThreeVector globalPreStepPoint = pPreStepPoint->GetPosition();
 
  int endCapID;
 
  if (globalPreStepPoint.z() > 0.)
    endCapID = 0;
  else
    endCapID = 1;
 
  if (m_sectorsABC)
    {
      if (wheelID < m_numberOfWheelsAB)
        {
          strawID += m_numberOfStrawsInSectorsAB * sectorID -
            m_initialStrawIDInSectorsAB;
          if (strawID < 0)
            strawID += m_numberOfStrawsInPlanesAB;
        }
      else
        {
          strawID += m_numberOfStrawsInSectorC * sectorID - m_initialStrawIDInSectorC;
          if (strawID < 0)
            strawID += m_numberOfStrawsInPlaneC;
        }
    }
 
  G4ThreeVector globalPostStepPoint = pPostStepPoint->GetPosition();
 
  const G4TouchableHistory* pTouchableHistory =
    static_cast<const G4TouchableHistory*>(pPreStepPoint->GetTouchable());
 
  const G4AffineTransform& topTransform = pTouchableHistory->GetHistory()->
    GetTopTransform();
 
  G4ThreeVector localPreStepPoint = topTransform.TransformPoint(globalPreStepPoint);
  G4ThreeVector localPostStepPoint = topTransform.TransformPoint(globalPostStepPoint);
 
  double preStepX = localPreStepPoint.x();
  double preStepY = localPreStepPoint.y();
  double preStepZ = localPreStepPoint.z();
  double postStepX = localPostStepPoint.x();
  double postStepY = localPostStepPoint.y();
  double postStepZ = localPostStepPoint.z();
 
 
  if (m_testLocalCoordinatesOfHits)
    {
      double preStepR = std::sqrt(preStepX * preStepX + preStepY * preStepY);
      double postStepR = std::sqrt(postStepX * postStepX + postStepY * postStepY);
 
      if (preStepR > 2.0000001 || postStepR > 2.0000001)
        {
          G4int particleEncoding = m_pSensitiveDetector->m_particleEncoding;
          G4double kineticEnergy = m_pSensitiveDetector->m_kineticEnergy;
          G4double energyDeposit = m_pSensitiveDetector->m_energyDeposit;
 
          std::cout << "!!!!! End-caps. Error in local coordinates of hits!" << std::endl;
          std::cout << "  endCapID=" << endCapID << "  wheelID=" << wheelID
                    << "  planeID=" << planeID << "  strawID=" << strawID
                    << "  trackID=" << trHelp.GetUniqueID() << std::endl;
          std::cout << "  particleEncoding=" << particleEncoding;
 
          if (kineticEnergy < 0.0001)
            std::cout << "  kineticEnergy=" << kineticEnergy * 1000000. << " eV";
          else if (kineticEnergy < 0.1)
            std::cout << "  kineticEnergy=" << kineticEnergy * 1000. << " keV";
          else if (kineticEnergy >= 100.)
            std::cout << "  kineticEnergy=" << kineticEnergy / 1000. << " GeV";
          else if (kineticEnergy >= 100000.)
            std::cout << "  kineticEnergy=" << kineticEnergy / 1000000. << " TeV";
          else
            std::cout << "  kineticEnergy=" << kineticEnergy << " MeV";
 
          if (energyDeposit < 0.1)
            std::cout << "  energyDeposit=" << energyDeposit * 1000. << " eV"
                      << std::endl;
          else if (energyDeposit >= 100.)
            std::cout << "  energyDeposit=" << energyDeposit / 1000. << " MeV"
                      << std::endl;
          else
            std::cout << "  energyDeposit=" << energyDeposit << " keV" << std::endl;
 
          std::cout << "  preStepX=" << preStepX << " mm  preStepY=" << preStepY
                    << " mm  preStepR=" << preStepR << " mm" << std::endl;
          std::cout << "  postStepX=" << postStepX << " mm  postStepY=" << postStepY
                    << " mm  postStepR=" << postStepR << " mm" << std::endl << std::endl;
        }
    }
 
  double preStepGlobalTime = pPreStepPoint->GetGlobalTime();
  double postStepGlobalTime = pPostStepPoint->GetGlobalTime();
  double globalTime = (preStepGlobalTime + postStepGlobalTime) / 2.;
 
 
  int numberOfStrawsInPlane;
  int numberOfStrawsInIDSector;
 
  if (wheelID < m_numberOfWheelsAB) {
    numberOfStrawsInPlane = m_numberOfStrawsInPlanesAB;
    numberOfStrawsInIDSector = m_numberOfStrawsInIDSectorsAB;
  } else {
    numberOfStrawsInPlane = m_numberOfStrawsInPlaneC;
    numberOfStrawsInIDSector = m_numberOfStrawsInIDSectorC;
  }
 
  // Mapping for negative endcap is
  //    nSectors = 32
  //    nStraws = num straws in sector
  //    straw =  straw number within sector
  //    sector -> (nSectors + nSectors/2 - sector - 1) % nSectors
  //    straw  -> nStraws - 1 - straw
  //
  // Since we calculate sector number below from strawID (where
  // strawID is straw number in plane), its equivalent to map strawID
  // first and then calculate sector number and straw number within sector.
 
  // Mapping straw number in plane for negative endcap
  if (endCapID > 0) {
    strawID = (3*numberOfStrawsInPlane/2 - strawID - 1) % numberOfStrawsInPlane;
  }
 
  // Calculate sector number and straw number within sector.
  sectorID = strawID / numberOfStrawsInIDSector;
  strawID %= numberOfStrawsInIDSector;
 
 
  int hitID;
  if (endCapID == 0)
    hitID = 0x00200000;
  else
    hitID = 0x00300000;
 
  hitID += (wheelID << 15);
  hitID += (sectorID << 10);
  hitID += (planeID << 5);
  hitID += strawID;
 
  m_pSensitiveDetector->m_hitID = hitID;
  auto* eventInfo = m_pSensitiveDetector->m_g4UserEventInfo;
  m_pSensitiveDetector->m_partLink =  trHelp.GenerateParticleLink(eventInfo ? eventInfo->GetEventStore() : nullptr);
  m_pSensitiveDetector->m_preStepX = preStepX;
  m_pSensitiveDetector->m_preStepY = preStepY;
  m_pSensitiveDetector->m_preStepZ = preStepZ;
  m_pSensitiveDetector->m_postStepX = postStepX;
  m_pSensitiveDetector->m_postStepY = postStepY;
  m_pSensitiveDetector->m_postStepZ = postStepZ;
  m_pSensitiveDetector->m_globalTime = globalTime;
 
  return true;
}