PltSD Class Reference

#include <PltSD.h>

Inheritance diagram for PltSD:

Public Member Functions
void	clearHits () override
void	EndOfEvent (G4HCofThisEvent *) override
void	fillHits (edm::PSimHitContainer &, const std::string &) override
	PltSD (const std::string &, const DDCompactView &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
bool	ProcessHits (G4Step *, G4TouchableHistory *) override
uint32_t	setDetUnitId (const G4Step *) override
	~PltSD () override
Public Member Functions inherited from SensitiveTkDetector
	SensitiveTkDetector (const std::string &iname, const DDCompactView &cpv, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)
Public Member Functions inherited from SensitiveDetector
void	EndOfEvent (G4HCofThisEvent *eventHC) override
const std::vector< std::string > &	getNames () const
void	Initialize (G4HCofThisEvent *eventHC) override
	SensitiveDetector (const std::string &iname, const DDCompactView &cpv, const SensitiveDetectorCatalog &, edm::ParameterSet const &p)
	~SensitiveDetector () override
Public Member Functions inherited from Observer< const BeginOfEvent * >
	Observer ()
void	slotForUpdate (const BeginOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfTrack * >
	Observer ()
void	slotForUpdate (const BeginOfTrack * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfJob * >
	Observer ()
void	slotForUpdate (const BeginOfJob * iT)
virtual	~Observer ()

Private Member Functions
virtual bool	closeHit (G4Step *)
virtual void	createHit (G4Step *)
TrackInformation *	getOrCreateTrackInformation (const G4Track *)
virtual bool	newHit (G4Step *)
virtual void	sendHit ()
void	update (const BeginOfEvent *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfTrack *) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfJob *) override
	This routine will be called when the appropriate signal arrives. More...
virtual void	updateHit (G4Step *)

Private Attributes
float	energyCut
float	energyHistoryCut
int	eventno
Local3DPoint	globalEntryPoint
Local3DPoint	globalExitPoint
uint32_t	lastId
unsigned int	lastTrack
G4TrackToParticleID *	myG4TrackToParticleID
UpdatablePSimHit *	mySimHit
G4VPhysicalVolume *	oldVolume
std::string	pname
TrackingSlaveSD *	slave
G4ProcessTypeEnumerator *	theG4ProcessTypeEnumerator

Additional Inherited Members
Protected Types inherited from SensitiveDetector
enum	coordinates { WorldCoordinates, LocalCoordinates }
Protected Member Functions inherited from SensitiveDetector
Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point) const
Local3DPoint	FinalStepPosition (const G4Step *step, coordinates) const
Local3DPoint	InitialStepPosition (const G4Step *step, coordinates) const
Local3DPoint	LocalPostStepPosition (const G4Step *step) const
Local3DPoint	LocalPreStepPosition (const G4Step *step) const
void	NaNTrap (const G4Step *step) const
void	setNames (const std::vector< std::string > &)

Detailed Description

Definition at line 30 of file PltSD.h.

Constructor & Destructor Documentation

PltSD::PltSD	(	const std::string &	name,
		const DDCompactView &	cpv,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

Definition at line 35 of file PltSD.cc.

References energyCut, energyHistoryCut, edm::ParameterSet::getParameter(), GeV, myG4TrackToParticleID, slave, and theG4ProcessTypeEnumerator.

                                             :
SensitiveTkDetector(name, cpv, clg, p), mySimHit(nullptr),
oldVolume(nullptr), lastId(0), lastTrack(0), eventno(0) {
    
    edm::ParameterSet m_TrackerSD = p.getParameter<edm::ParameterSet>("PltSD");
    energyCut           = m_TrackerSD.getParameter<double>("EnergyThresholdForPersistencyInGeV")*GeV; //default must be 0.5 (?)
    energyHistoryCut    = m_TrackerSD.getParameter<double>("EnergyThresholdForHistoryInGeV")*GeV;//default must be 0.05 (?)
    
    edm::LogInfo("PltSD") <<"Criteria for Saving Tracker SimTracks: \n "
    <<" History: "<<energyHistoryCut<< " MeV ; Persistency: "<< energyCut<<" MeV\n"
    <<" Constructing a PltSD";
    
    slave  = new TrackingSlaveSD(name);
        
    theG4ProcessTypeEnumerator = new G4ProcessTypeEnumerator;
    myG4TrackToParticleID = new G4TrackToParticleID;
}

PltSD::~PltSD ( )

override

Definition at line 57 of file PltSD.cc.

References myG4TrackToParticleID, slave, and theG4ProcessTypeEnumerator.

              {
    delete slave;
    delete theG4ProcessTypeEnumerator;
    delete myG4TrackToParticleID;
}

Member Function Documentation

void PltSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 319 of file PltSD.cc.

References TrackingSlaveSD::Initialize(), and slave.

Referenced by update().

                      {
    slave->Initialize();
}

bool PltSD::closeHit ( G4Step *  aStep )

privatevirtual

Definition at line 246 of file PltSD.cc.

References PSimHit::exitPoint(), SensitiveDetector::InitialStepPosition(), SensitiveDetector::LocalCoordinates, LogDebug, mag(), mySimHit, and RecoTauCleanerPlugins::tolerance.

Referenced by newHit().

                                   {
    
    if (mySimHit == nullptr) return false;
    const float tolerance = 0.05 * mm; // 50 micron are allowed between the exit
    // point of the current hit and the entry point of the new hit
    Local3DPoint theEntryPoint = SensitiveDetector::InitialStepPosition(aStep,LocalCoordinates);
    LogDebug("PltSD")<< " closeHit: distance = " << (mySimHit->exitPoint()-theEntryPoint).mag();
    
    if ((mySimHit->exitPoint()-theEntryPoint).mag()<tolerance) return true;
    return false;
}

void PltSD::createHit ( G4Step *  aStep )

privatevirtual

Definition at line 258 of file PltSD.cc.

References SensitiveDetector::ConvertToLocal3DPoint(), PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), SensitiveDetector::FinalStepPosition(), GeV, globalEntryPoint, globalExitPoint, SensitiveDetector::InitialStepPosition(), lastId, lastTrack, SensitiveDetector::LocalCoordinates, LogDebug, myG4TrackToParticleID, mySimHit, oldVolume, G4TrackToParticleID::particleID(), PV3DBase< T, PVType, FrameType >::phi(), pname, G4ProcessTypeEnumerator::processId(), setDetUnitId(), theG4ProcessTypeEnumerator, PV3DBase< T, PVType, FrameType >::theta(), PSimHit::trackId(), findQualityFiles::v, and SensitiveDetector::WorldCoordinates.

Referenced by ProcessHits().

                                    {
    
    if (mySimHit != nullptr) {
        delete mySimHit;
        mySimHit=nullptr;
    }
    
    G4Track * theTrack  = aStep->GetTrack();
    G4VPhysicalVolume * v = aStep->GetPreStepPoint()->GetPhysicalVolume();
    
    Local3DPoint theEntryPoint = SensitiveDetector::InitialStepPosition(aStep,LocalCoordinates);
    Local3DPoint theExitPoint  = SensitiveDetector::FinalStepPosition(aStep,LocalCoordinates);
    
    float thePabs             = aStep->GetPreStepPoint()->GetMomentum().mag()/GeV;
    float theTof              = aStep->GetPreStepPoint()->GetGlobalTime()/nanosecond;
    float theEnergyLoss       = aStep->GetTotalEnergyDeposit()/GeV;
    int theParticleType       = myG4TrackToParticleID->particleID(theTrack);
    uint32_t theDetUnitId     = setDetUnitId(aStep);
    
    globalEntryPoint = SensitiveDetector::InitialStepPosition(aStep,WorldCoordinates);
    globalExitPoint = SensitiveDetector::FinalStepPosition(aStep,WorldCoordinates);
    pname = theTrack->GetDynamicParticle()->GetDefinition()->GetParticleName();
    
    unsigned int theTrackID = theTrack->GetTrackID();
    
    G4ThreeVector gmd  = aStep->GetPreStepPoint()->GetMomentumDirection();
    // convert it to local frame
    G4ThreeVector lmd = ((G4TouchableHistory *)(aStep->GetPreStepPoint()->GetTouchable()))->GetHistory()->GetTopTransform().TransformAxis(gmd);
    Local3DPoint lnmd = ConvertToLocal3DPoint(lmd);
    float theThetaAtEntry = lnmd.theta();
    float thePhiAtEntry = lnmd.phi();
    
    mySimHit = new UpdatablePSimHit(theEntryPoint,theExitPoint,thePabs,theTof,
                                    theEnergyLoss,theParticleType,theDetUnitId,
                                    theTrackID,theThetaAtEntry,thePhiAtEntry,
                                    theG4ProcessTypeEnumerator->processId(theTrack->GetCreatorProcess()));
    
    LogDebug("PltSD") << " Created PSimHit: " << pname << " "
    << mySimHit->detUnitId() << " " << mySimHit->trackId()
    << " " << mySimHit->energyLoss() << " "
    << mySimHit->entryPoint() << " " << mySimHit->exitPoint();
    lastId = theDetUnitId;
    lastTrack = theTrackID;
    oldVolume = v;
}

void PltSD::EndOfEvent ( G4HCofThisEvent *  )

override

Definition at line 188 of file PltSD.cc.

References LogDebug, mySimHit, and sendHit().

                                        {
    
    LogDebug("PltSD")<< " Saving the last hit in a ROU " << GetName();
    
    if (mySimHit == nullptr) return;
    sendHit();
}

void PltSD::fillHits ( edm::PSimHitContainer &  cc, const std::string &  hname  )

overridevirtual

Implements SensitiveTkDetector.

Definition at line 196 of file PltSD.cc.

References TrackingSlaveSD::hits(), TrackingSlaveSD::name(), and slave.

                                                                   {
  if (slave->name() == hname) { cc=slave->hits(); }
}

TrackInformation * PltSD::getOrCreateTrackInformation ( const G4Track *  gTrack )

private

Definition at line 323 of file PltSD.cc.

References info(), and groupFilesInBlocks::temp.

Referenced by ProcessHits().

                                                                           {
    G4VUserTrackInformation* temp = gTrack->GetUserInformation();
    if (temp == nullptr){
        edm::LogError("PltSD") <<" ERROR: no G4VUserTrackInformation available";
        abort();
    }else{
        TrackInformation* info = dynamic_cast<TrackInformation*>(temp);
        if (info == nullptr){
            edm::LogError("PltSD") <<" ERROR: TkSimTrackSelection: the UserInformation does not appear to be a TrackInformation";
            abort();
        }
        return info;
    }
}

bool PltSD::newHit ( G4Step *  aStep )

privatevirtual

Definition at line 232 of file PltSD.cc.

References closeHit(), lastId, lastTrack, LogDebug, mySimHit, and setDetUnitId().

Referenced by ProcessHits().

                                 {
    
    G4Track * theTrack = aStep->GetTrack();
    uint32_t theDetUnitId = setDetUnitId(aStep);
    unsigned int theTrackID = theTrack->GetTrackID();
    
    LogDebug("PltSD") << " OLD (d,t) = (" << lastId << "," << lastTrack
    << "), new = (" << theDetUnitId << "," << theTrackID << ") return "
    << ((theTrackID == lastTrack) && (lastId == theDetUnitId));
    if ((mySimHit != nullptr) && (theTrackID == lastTrack) && (lastId == theDetUnitId) && closeHit(aStep))
    return false;
    return true;
}

bool PltSD::ProcessHits ( G4Step *  aStep, G4TouchableHistory *    )

overridevirtual

Implements SensitiveDetector.

Definition at line 63 of file PltSD.cc.

References createHit(), energyCut, energyHistoryCut, getOrCreateTrackInformation(), info(), lastTrack, LogDebug, newHit(), TrackInformation::putInHistory(), sendHit(), TrackInformation::storeTrack(), and updateHit().

                                                             {
    
    LogDebug("PltSD") << " Entering a new Step "
    << aStep->GetTotalEnergyDeposit() << " "
    << aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName();
    
    G4Track * gTrack  = aStep->GetTrack();
    if ((unsigned int)(gTrack->GetTrackID()) != lastTrack) {
        
        if (gTrack->GetKineticEnergy() > energyCut){
            TrackInformation* info = getOrCreateTrackInformation(gTrack);
            info->storeTrack(true);
        }
        if (gTrack->GetKineticEnergy() > energyHistoryCut){
            TrackInformation* info = getOrCreateTrackInformation(gTrack);
            info->putInHistory();
        }
    }
    
    if (aStep->GetTotalEnergyDeposit()>0.) {
        if (newHit(aStep) == true) {
            sendHit();
            createHit(aStep);
        } else {
            updateHit(aStep);
        }
        return true;
    }
    return false;
}

void PltSD::sendHit ( )

privatevirtual

Definition at line 200 of file PltSD.cc.

References PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), lastId, lastTrack, LogDebug, mySimHit, pname, TrackingSlaveSD::processHits(), slave, and PSimHit::trackId().

Referenced by EndOfEvent(), and ProcessHits().

                    {
    if (mySimHit == nullptr) return;
    LogDebug("PltSD") << " Storing PSimHit: " << pname << " " << mySimHit->detUnitId()
    << " " << mySimHit->trackId() << " " << mySimHit->energyLoss()
    << " " << mySimHit->entryPoint() << " " << mySimHit->exitPoint();
    
    slave->processHits(*mySimHit);
    
    // clean up
    delete mySimHit;
    mySimHit = nullptr;
    lastTrack = 0;
    lastId = 0;
}

uint32_t PltSD::setDetUnitId ( const G4Step *  aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 94 of file PltSD.cc.

References gather_cfg::cout, hcalDigis_cfi::level, LogDebug, and groupFilesInBlocks::temp.

Referenced by createHit(), and newHit().

                                                  {
    
    unsigned int detId = 0;
    
    LogDebug("PltSD")<< " DetID = "<<detId;
    
    //Find number of levels
    const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
    int level = 0;
    if (touch) level = ((touch->GetHistoryDepth())+1);
    
    //Get name and copy numbers
    if ( level > 1 ) {
        //some debugging with the names
        G4String sensorName   = touch->GetVolume(2)->GetName();
        G4String telName      = touch->GetVolume(3)->GetName();
        G4String volumeName   = touch->GetVolume(4)->GetName();
        if ( sensorName != "PLTSensorPlane" )
        std::cout << " PltSD::setDetUnitId -w- Sensor name not PLTSensorPlane " << std::endl;
        if ( telName != "Telescope" )
        std::cout << " PltSD::setDetUnitId -w- Telescope name not Telescope " << std::endl;
        if ( volumeName != "PLT" )
        std::cout << " PltSD::setDetUnitId -w- Volume name not PLT " << std::endl;
        
        //Get the information about which telescope, plane, row/column was hit
        int columnNum = touch->GetReplicaNumber(0);
        int rowNum = touch->GetReplicaNumber(1);
        int sensorNum = touch->GetReplicaNumber(2);
        int telNum  = touch->GetReplicaNumber(3);
        //temp stores the PLTBCM volume the hit occured in (i.e. was the hit on the + or -z side?)
        int temp = touch->GetReplicaNumber(5);
        //convert to the PLT hit id standard
        int pltNum;
        if (temp == 2) pltNum = 0;
        else pltNum = 1;

        //correct the telescope numbers on the -z side to have the same naming convention in phi as the +z side
        if (pltNum == 0){
            if (telNum == 0){
                telNum = 7;
            }
            else if (telNum == 1){
                telNum = 6;
            }
            else if (telNum == 2){
                telNum = 5;
            }
            else if (telNum == 3){
                telNum = 4;
            }
            else if (telNum == 4){
                telNum = 3;
            }
            else if (telNum == 5){
                telNum = 2;
            }
            else if (telNum == 6){
                telNum = 1;
            }
            else if (telNum == 7){
                telNum = 0;
            }
        }
        //the PLT is divided into sets of telescopes on the + and -x sides
        int halfCarriageNum = -1;

        //If the telescope is on the -x side of the carriage, halfCarriageNum=0.  If on the +x side, it is = 1.
        if(telNum == 0 || telNum == 1 || telNum == 2 || telNum == 3)
            halfCarriageNum = 0;
        else
            halfCarriageNum = 1;
        //correct the telescope numbers of the +x half-carriage to range from 0 to 3
        if(halfCarriageNum == 1){
            if(telNum == 4){
                telNum = 0;
            }
            else if (telNum == 5){
                telNum = 1;
            }
            else if (telNum == 6){
                telNum = 2;
            }
            else if (telNum == 7){
                telNum = 3;
            }
        }
        //Define unique detId for each pixel.  See https://twiki.cern.ch/twiki/bin/viewauth/CMS/PLTSimulationGuide for more information
        detId = 10000000*pltNum+1000000*halfCarriageNum+100000*telNum+10000*sensorNum+100*rowNum+columnNum;
        //std::cout <<  "Hit Recorded at " << "plt:" << pltNum << " hc:" << halfCarriageNum << " tel:" << telNum << " plane:" << sensorNum << std::endl;
    }
    
    return detId;
}

void PltSD::update ( const BeginOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 306 of file PltSD.cc.

References clearHits(), eventno, and mySimHit.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                         {
    
    clearHits();
    eventno = (*i)()->GetEventID();
    mySimHit = nullptr;
}

void PltSD::update ( const BeginOfTrack *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 313 of file PltSD.cc.

References pname.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                          {
    
    const G4Track* gTrack = (*bot)();
    pname = gTrack->GetDynamicParticle()->GetDefinition()->GetParticleName();
}

void PltSD::update ( const BeginOfJob *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 304 of file PltSD.cc.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{ }

void PltSD::updateHit ( G4Step *  aStep )

privatevirtual

Definition at line 215 of file PltSD.cc.

References UpdatablePSimHit::addEnergyLoss(), PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), SensitiveDetector::FinalStepPosition(), GeV, globalExitPoint, SensitiveDetector::LocalCoordinates, LogDebug, mySimHit, pname, UpdatablePSimHit::setExitPoint(), PSimHit::trackId(), and SensitiveDetector::WorldCoordinates.

Referenced by ProcessHits().

                                    {
    
    Local3DPoint theExitPoint = SensitiveDetector::FinalStepPosition(aStep,LocalCoordinates);
    float theEnergyLoss = aStep->GetTotalEnergyDeposit()/GeV;
    mySimHit->setExitPoint(theExitPoint);
    LogDebug("PltSD")<< " Before : " << mySimHit->energyLoss();
    mySimHit->addEnergyLoss(theEnergyLoss);
    globalExitPoint = SensitiveDetector::FinalStepPosition(aStep,WorldCoordinates);
    
    LogDebug("PltSD") << " Updating: new exitpoint " << pname << " "
    << theExitPoint << " new energy loss " << theEnergyLoss
    << "\n Updated PSimHit: " << mySimHit->detUnitId()
    << " " << mySimHit->trackId()
    << " " << mySimHit->energyLoss() << " "
    << mySimHit->entryPoint() << " " << mySimHit->exitPoint();
}

Member Data Documentation

float PltSD::energyCut

private

Definition at line 67 of file PltSD.h.

Referenced by PltSD(), and ProcessHits().

float PltSD::energyHistoryCut

private

Definition at line 68 of file PltSD.h.

Referenced by PltSD(), and ProcessHits().

int PltSD::eventno

private

Definition at line 75 of file PltSD.h.

Referenced by update().

Local3DPoint PltSD::globalEntryPoint

private

Definition at line 70 of file PltSD.h.

Referenced by createHit().

Local3DPoint PltSD::globalExitPoint

private

Definition at line 71 of file PltSD.h.

Referenced by createHit(), and updateHit().

uint32_t PltSD::lastId

private

Definition at line 73 of file PltSD.h.

Referenced by createHit(), newHit(), and sendHit().

unsigned int PltSD::lastTrack

private

Definition at line 74 of file PltSD.h.

Referenced by createHit(), newHit(), ProcessHits(), and sendHit().

G4TrackToParticleID* PltSD::myG4TrackToParticleID

private

Definition at line 65 of file PltSD.h.

Referenced by createHit(), PltSD(), and ~PltSD().

UpdatablePSimHit* PltSD::mySimHit

private

Definition at line 66 of file PltSD.h.

Referenced by closeHit(), createHit(), EndOfEvent(), newHit(), sendHit(), update(), and updateHit().

G4VPhysicalVolume* PltSD::oldVolume

private

Definition at line 72 of file PltSD.h.

Referenced by createHit().

std::string PltSD::pname

private

Definition at line 76 of file PltSD.h.

Referenced by createHit(), sendHit(), update(), and updateHit().

TrackingSlaveSD* PltSD::slave

private

Definition at line 63 of file PltSD.h.

Referenced by clearHits(), fillHits(), PltSD(), sendHit(), and ~PltSD().

G4ProcessTypeEnumerator* PltSD::theG4ProcessTypeEnumerator

private

Definition at line 64 of file PltSD.h.

Referenced by createHit(), PltSD(), and ~PltSD().