TkAccumulatingSensitiveDetector Class Reference

#include <TkAccumulatingSensitiveDetector.h>

Inheritance diagram for TkAccumulatingSensitiveDetector:

Public Member Functions
void	clearHits () override
void	EndOfEvent (G4HCofThisEvent *) override
void	fillHits (edm::PSimHitContainer &, const std::string &) override
bool	ProcessHits (G4Step *, G4TouchableHistory *) override
uint32_t	setDetUnitId (const G4Step *) override
	TkAccumulatingSensitiveDetector (const std::string &, const edm::EventSetup &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
	~TkAccumulatingSensitiveDetector () override
Public Member Functions inherited from SensitiveTkDetector
	SensitiveTkDetector (const std::string &iname, const edm::EventSetup &es, 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
bool	isCaloSD () const
	SensitiveDetector (const std::string &iname, const edm::EventSetup &es, const SensitiveDetectorCatalog &, edm::ParameterSet const &p, bool calo)
	~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 ()

Protected Member Functions
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...
Protected Member Functions inherited from SensitiveDetector
TrackInformation *	cmsTrackInformation (const G4Track *aTrack)
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 > &)

Private Member Functions
bool	closeHit (const G4Step *)
void	createHit (const G4Step *)
bool	newHit (const G4Step *)
void	sendHit ()
void	updateHit (const G4Step *)

Private Attributes
bool	allowZeroEnergyLoss
float	energyCut
float	energyHistoryCut
int	eventno
Local3DPoint	globalEntryPoint
Local3DPoint	globalExitPoint
uint32_t	lastId
int	lastTrack
UpdatablePSimHit *	mySimHit
bool	neverAccumulate
const G4VPhysicalVolume *	oldVolume
std::string	pname
bool	printHits
float	px
float	py
float	pz
double	rTracker
double	rTracker2
std::unique_ptr< TrackingSlaveSD >	slaveHighTof
std::unique_ptr< TrackingSlaveSD >	slaveLowTof
std::unique_ptr< const G4ProcessTypeEnumerator >	theG4ProcTypeEnumerator
const SimTrackManager *	theManager
TrackerG4SimHitNumberingScheme *	theNumberingScheme
std::unique_ptr< FrameRotation >	theRotation
float	theTofLimit
double	zTracker

Additional Inherited Members
Protected Types inherited from SensitiveDetector
enum	coordinates { WorldCoordinates, LocalCoordinates }

Detailed Description

Definition at line 28 of file TkAccumulatingSensitiveDetector.h.

Constructor & Destructor Documentation

TkAccumulatingSensitiveDetector::TkAccumulatingSensitiveDetector	(	const std::string &	name,
		const edm::EventSetup &	es,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

Definition at line 44 of file TkAccumulatingSensitiveDetector.cc.

References allowZeroEnergyLoss, energyCut, energyHistoryCut, edm::ParameterSet::getParameter(), GeV, neverAccumulate, printHits, rTracker, rTracker2, SensitiveDetector::setNames(), slaveHighTof, slaveLowTof, AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, theG4ProcTypeEnumerator, theNumberingScheme, theRotation, theTofLimit, and zTracker.

    : SensitiveTkDetector(name, es, clg, p),
      theManager(manager),
      rTracker(1200. * CLHEP::mm),
      zTracker(3000. * CLHEP::mm),
      mySimHit(nullptr),
      lastId(0),
      lastTrack(0),
      oldVolume(nullptr),
      px(0.0f),
      py(0.0f),
      pz(0.0f),
      eventno(0),
      pname("") {
  edm::ParameterSet m_TrackerSD = p.getParameter<edm::ParameterSet>("TrackerSD");
  allowZeroEnergyLoss = m_TrackerSD.getParameter<bool>("ZeroEnergyLoss");
  neverAccumulate = m_TrackerSD.getParameter<bool>("NeverAccumulate");
  printHits = m_TrackerSD.getParameter<bool>("PrintHits");
  theTofLimit = m_TrackerSD.getParameter<double>("ElectronicSigmaInNanoSeconds") * 3 * CLHEP::ns;  // 3 sigma
  energyCut =
      m_TrackerSD.getParameter<double>("EnergyThresholdForPersistencyInGeV") * CLHEP::GeV;  //default must be 0.5
  energyHistoryCut =
      m_TrackerSD.getParameter<double>("EnergyThresholdForHistoryInGeV") * CLHEP::GeV;  //default must be 0.05
  rTracker2 = rTracker * rTracker;

  // No Rotation given in input, automagically choose one based upon the name
  std::string rotType;
  theRotation.reset(new TrackerFrameRotation());
  rotType = "TrackerFrameRotation";

#ifdef FAKEFRAMEROTATION
  theRotation.reset(new FakeFrameRotation());
  rotType = "FakeFrameRotation";
#endif

  edm::LogInfo("TrackerSimInfo") << " TkAccumulatingSensitiveDetector: "
                                 << " Criteria for Saving Tracker SimTracks: \n"
                                 << " History: " << energyHistoryCut << " MeV; Persistency: " << energyCut
                                 << " MeV;  TofLimit: " << theTofLimit << " ns"
                                 << "\n FrameRotation type " << rotType << " rTracker(cm)= " << rTracker / CLHEP::cm
                                 << " zTracker(cm)= " << zTracker / CLHEP::cm
                                 << " allowZeroEnergyLoss: " << allowZeroEnergyLoss
                                 << " neverAccumulate: " << neverAccumulate << " printHits: " << printHits;

  slaveLowTof.reset(new TrackingSlaveSD(name + "LowTof"));
  slaveHighTof.reset(new TrackingSlaveSD(name + "HighTof"));

  std::vector<std::string> temp;
  temp.push_back(slaveLowTof.get()->name());
  temp.push_back(slaveHighTof.get()->name());
  setNames(temp);

  theG4ProcTypeEnumerator.reset(new G4ProcessTypeEnumerator);
  theNumberingScheme = nullptr;
}

TkAccumulatingSensitiveDetector::~TkAccumulatingSensitiveDetector ( )

override

Definition at line 104 of file TkAccumulatingSensitiveDetector.cc.

{}

Member Function Documentation

void TkAccumulatingSensitiveDetector::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 362 of file TkAccumulatingSensitiveDetector.cc.

References slaveHighTof, and slaveLowTof.

Referenced by update().

                                                {
  slaveLowTof.get()->Initialize();
  slaveHighTof.get()->Initialize();
}

bool TkAccumulatingSensitiveDetector::closeHit ( const G4Step *  aStep )

private

Definition at line 333 of file TkAccumulatingSensitiveDetector.cc.

References PSimHit::exitPoint(), SensitiveDetector::LocalPreStepPosition(), LogDebug, mag(), mag2(), mySimHit, and theRotation.

Referenced by newHit().

                                                                  {
  const float tolerance2 = 0.0025f;  // (0.5 mm)^2 are allowed between entry and exit
  Local3DPoint theEntryPoint = theRotation.get()->transformPoint(LocalPreStepPosition(aStep));
  LogDebug("TrackerSimDebug") << " closeHit: distance = " << (mySimHit->exitPoint() - theEntryPoint).mag();

  return ((mySimHit->exitPoint() - theEntryPoint).mag2() < tolerance2) ? true : false;
}

void TkAccumulatingSensitiveDetector::createHit ( const G4Step *  aStep )

private

Definition at line 206 of file TkAccumulatingSensitiveDetector.cc.

References SensitiveDetector::cmsTrackInformation(), SensitiveDetector::ConvertToLocal3DPoint(), PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), Exception, PSimHit::exitPoint(), GeV, globalEntryPoint, globalExitPoint, lastId, lastTrack, SensitiveDetector::LocalPostStepPosition(), SensitiveDetector::LocalPreStepPosition(), LogDebug, mySimHit, oldVolume, G4TrackToParticleID::particleID(), PV3DBase< T, PVType, FrameType >::phi(), pname, printHits, px, py, pz, setDetUnitId(), TrackInformation::storeTrack(), groupFilesInBlocks::temp, theG4ProcTypeEnumerator, theRotation, PV3DBase< T, PVType, FrameType >::theta(), and PSimHit::trackId().

Referenced by ProcessHits().

                                                                   {
  // VI: previous hit should be already deleted
  //     in past here was a check if a hit is inside a sensitive detector,
  //     this is not needed, because call to senstive detector happens
  //     only inside the volume
  const G4Track* theTrack = aStep->GetTrack();
  Local3DPoint theExitPoint = theRotation.get()->transformPoint(LocalPostStepPosition(aStep));
  Local3DPoint theEntryPoint;
  //
  //  Check particle type - for gamma and neutral hadrons energy deposition
  //  should be local (VI)
  //
  if (0.0 == theTrack->GetDefinition()->GetPDGCharge()) {
    theEntryPoint = theExitPoint;
  } else {
    theEntryPoint = theRotation.get()->transformPoint(LocalPreStepPosition(aStep));
  }

  //
  //    This allows to send he skipEvent if it is outside!
  //
  const G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
  float thePabs = preStepPoint->GetMomentum().mag() / GeV;
  float theTof = preStepPoint->GetGlobalTime() / nanosecond;
  float theEnergyLoss = aStep->GetTotalEnergyDeposit() / GeV;
  int theParticleType = G4TrackToParticleID::particleID(theTrack);
  uint32_t theDetUnitId = setDetUnitId(aStep);
  int theTrackID = theTrack->GetTrackID();
  if (theDetUnitId == 0) {
    edm::LogWarning("TkAccumulatingSensitiveDetector::createHit") << " theDetUnitId is not valid for " << GetName();
    throw cms::Exception("TkAccumulatingSensitiveDetector::createHit")
        << "cannot get theDetUnitId for G4Track " << theTrackID;
  }

  // To whom assign the Hit?
  // First iteration: if the track is to be stored, use the current number;
  // otherwise, get to the mother
  unsigned int theTrackIDInsideTheSimHit = theTrackID;

  const TrackInformation* temp = cmsTrackInformation(theTrack);
  if (!temp->storeTrack()) {
    // Go to the mother!
    theTrackIDInsideTheSimHit = theTrack->GetParentID();
    LogDebug("TrackerSimDebug") << " TkAccumulatingSensitiveDetector::createHit(): setting the TrackID from "
                                << theTrackIDInsideTheSimHit << " to the mother one " << theTrackIDInsideTheSimHit
                                << " " << theEnergyLoss;
  } else {
    LogDebug("TrackerSimDebug") << " TkAccumulatingSensitiveDetector:createHit(): leaving the current TrackID "
                                << theTrackIDInsideTheSimHit;
  }

  const G4ThreeVector& gmd = preStepPoint->GetMomentumDirection();
  // convert it to local frame
  G4ThreeVector lmd =
      ((G4TouchableHistory*)(preStepPoint->GetTouchable()))->GetHistory()->GetTopTransform().TransformAxis(gmd);
  Local3DPoint lnmd = theRotation.get()->transformPoint(ConvertToLocal3DPoint(lmd));
  float theThetaAtEntry = lnmd.theta();
  float thePhiAtEntry = lnmd.phi();

  mySimHit = new UpdatablePSimHit(theEntryPoint,
                                  theExitPoint,
                                  thePabs,
                                  theTof,
                                  theEnergyLoss,
                                  theParticleType,
                                  theDetUnitId,
                                  theTrackIDInsideTheSimHit,
                                  theThetaAtEntry,
                                  thePhiAtEntry,
                                  theG4ProcTypeEnumerator.get()->processId(theTrack->GetCreatorProcess()));
  lastId = theDetUnitId;
  lastTrack = theTrackID;

  // only for debugging
  if (printHits) {
    // point on Geant4 unit (mm)
    globalEntryPoint = ConvertToLocal3DPoint(preStepPoint->GetPosition());
    globalExitPoint = ConvertToLocal3DPoint(aStep->GetPostStepPoint()->GetPosition());
    // in CMS unit (GeV)
    px = preStepPoint->GetMomentum().x() / CLHEP::GeV;
    py = preStepPoint->GetMomentum().y() / CLHEP::GeV;
    pz = preStepPoint->GetMomentum().z() / CLHEP::GeV;
    oldVolume = preStepPoint->GetPhysicalVolume();
    pname = theTrack->GetDefinition()->GetParticleName();
    LogDebug("TrackerSimDebug") << " Created PSimHit: " << pname << " " << mySimHit->detUnitId() << " "
                                << mySimHit->trackId() << " " << theTrackID
                                << " p= " << aStep->GetPreStepPoint()->GetMomentum().mag() << " "
                                << mySimHit->energyLoss() << " " << mySimHit->entryPoint() << " "
                                << mySimHit->exitPoint();
  }
}

void TkAccumulatingSensitiveDetector::EndOfEvent ( G4HCofThisEvent *  )

override

Definition at line 341 of file TkAccumulatingSensitiveDetector.cc.

References LogDebug, mySimHit, and sendHit().

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

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

overridevirtual

Implements SensitiveTkDetector.

Definition at line 367 of file TkAccumulatingSensitiveDetector.cc.

References slaveHighTof, and slaveLowTof.

                                                                                              {
  if (slaveLowTof.get()->name() == hname) {
    cc = slaveLowTof.get()->hits();
  } else if (slaveHighTof.get()->name() == hname) {
    cc = slaveHighTof.get()->hits();
  }
}

bool TkAccumulatingSensitiveDetector::newHit ( const G4Step *  aStep )

private

Definition at line 316 of file TkAccumulatingSensitiveDetector.cc.

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

Referenced by ProcessHits().

                                                                {
  const G4Track* theTrack = aStep->GetTrack();

  // for neutral particles do not merge hits (V.I.)
  if (0.0 == theTrack->GetDefinition()->GetPDGCharge())
    return true;

  uint32_t theDetUnitId = setDetUnitId(aStep);
  int theTrackID = theTrack->GetTrackID();

  LogDebug("TrackerSimDebug") << "newHit: OLD(detID,trID) = (" << lastId << "," << lastTrack << "), NEW = ("
                              << theDetUnitId << "," << theTrackID << ") Step length(mm)= " << aStep->GetStepLength()
                              << " Edep= " << aStep->GetTotalEnergyDeposit()
                              << " p= " << aStep->GetPreStepPoint()->GetMomentum().mag();
  return ((theTrackID == lastTrack) && (lastId == theDetUnitId) && closeHit(aStep)) ? false : true;
}

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

overridevirtual

Implements SensitiveDetector.

Definition at line 106 of file TkAccumulatingSensitiveDetector.cc.

References allowZeroEnergyLoss, createHit(), LogDebug, mySimHit, neverAccumulate, newHit(), sendHit(), and updateHit().

Referenced by LaserAlignmentSimulation::update().

                                                                                    {
  LogDebug("TrackerSimDebug") << " Entering a new Step " << aStep->GetTotalEnergyDeposit() << " "
                              << aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName();

  if (aStep->GetTotalEnergyDeposit() > 0. || allowZeroEnergyLoss) {
    if (!mySimHit) {
      createHit(aStep);
    } else if (neverAccumulate || newHit(aStep)) {
      sendHit();
      createHit(aStep);
    } else {
      updateHit(aStep);
    }
    return true;
  }
  return false;
}

void TkAccumulatingSensitiveDetector::sendHit ( )

private

Definition at line 173 of file TkAccumulatingSensitiveDetector.cc.

References PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), eventno, PSimHit::exitPoint(), globalEntryPoint, globalExitPoint, lastId, lastTrack, LogDebug, mySimHit, oldVolume, PSimHit::pabs(), pname, TkSimHitPrinter::printGlobal(), TkSimHitPrinter::printGlobalMomentum(), TkSimHitPrinter::printHitData(), printHits, TkSimHitPrinter::printLocal(), px, py, pz, slaveHighTof, slaveLowTof, TkSimHitPrinter::startNewSimHit(), theTofLimit, PSimHit::timeOfFlight(), and PSimHit::trackId().

Referenced by EndOfEvent(), and ProcessHits().

                                              {
  if (mySimHit == nullptr)
    return;
  if (printHits) {
    TkSimHitPrinter thePrinter("TkHitPositionOSCAR.dat");
    thePrinter.startNewSimHit(GetName(),
                              oldVolume->GetLogicalVolume()->GetName(),
                              mySimHit->detUnitId(),
                              mySimHit->trackId(),
                              lastTrack,
                              eventno);
    thePrinter.printLocal(mySimHit->entryPoint(), mySimHit->exitPoint());
    thePrinter.printGlobal(globalEntryPoint, globalExitPoint);
    thePrinter.printHitData(pname, mySimHit->pabs(), mySimHit->energyLoss(), mySimHit->timeOfFlight());
    thePrinter.printGlobalMomentum(px, py, pz);
    LogDebug("TrackerSimDebug") << " Storing PSimHit: " << mySimHit->detUnitId() << " " << mySimHit->trackId() << " "
                                << mySimHit->energyLoss() << " " << mySimHit->entryPoint() << " "
                                << mySimHit->exitPoint();
  }

  if (mySimHit->timeOfFlight() < theTofLimit) {
    slaveLowTof.get()->processHits(*mySimHit);  // implicit conversion (slicing) to PSimHit!!!
  } else {
    slaveHighTof.get()->processHits(*mySimHit);  // implicit conversion (slicing) to PSimHit!!!
  }
  //
  // clean up
  delete mySimHit;
  mySimHit = nullptr;
  lastTrack = 0;
  lastId = 0;
}

uint32_t TkAccumulatingSensitiveDetector::setDetUnitId ( const G4Step *  step )

overridevirtual

Implements SensitiveDetector.

Definition at line 124 of file TkAccumulatingSensitiveDetector.cc.

References TrackerG4SimHitNumberingScheme::g4ToNumberingScheme(), and theNumberingScheme.

Referenced by createHit(), and newHit().

                                                                         {
  return theNumberingScheme->g4ToNumberingScheme(step->GetPreStepPoint()->GetTouchable());
}

void TkAccumulatingSensitiveDetector::update ( const BeginOfEvent *  )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 347 of file TkAccumulatingSensitiveDetector.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();
  delete mySimHit;
  mySimHit = nullptr;
}

void TkAccumulatingSensitiveDetector::update ( const BeginOfTrack *  )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 128 of file TkAccumulatingSensitiveDetector.cc.

References funct::abs(), SensitiveDetector::cmsTrackInformation(), energyCut, energyHistoryCut, info(), LogDebug, TrackInformation::putInHistory(), rTracker2, TrackInformation::storeTrack(), and zTracker.

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)();

#ifdef DUMPPROCESSES
  if (gTrack->GetCreatorProcess()) {
    edm::LogVerbatim("TrackerSimInfo") << " -> PROCESS CREATOR : " << gTrack->GetCreatorProcess()->GetProcessName();
  } else {
    edm::LogVerbatim("TrackerSimInfo") << " -> No Creator process";
  }
#endif

  //
  //Position
  //
  const G4ThreeVector& pos = gTrack->GetPosition();
  LogDebug("TrackerSimDebug") << " update(..) of " << gTrack->GetDefinition()->GetParticleName()
                              << " trackID= " << gTrack->GetTrackID() << " E(MeV)= " << gTrack->GetKineticEnergy()
                              << " Ecut= " << energyCut << " R(mm)= " << pos.perp() << " Z(mm)= " << pos.z();

  //
  // Check if in Tracker Volume
  //
  if (pos.x() * pos.x() + pos.y() * pos.y() < rTracker2 && std::abs(pos.z()) < zTracker) {
    //
    // inside the Tracker
    //
    TrackInformation* info = nullptr;
    if (gTrack->GetKineticEnergy() > energyCut) {
      info = cmsTrackInformation(gTrack);
      info->storeTrack(true);
    }
    //
    // Save History?
    //
    if (gTrack->GetKineticEnergy() > energyHistoryCut) {
      if (!info) {
        info = cmsTrackInformation(gTrack);
      }
      info->putInHistory();
      LogDebug("TrackerSimDebug") << " Track inside the tracker selected for HISTORY"
                                  << " Track ID= " << gTrack->GetTrackID();
    }
  }
}

void TkAccumulatingSensitiveDetector::update ( const BeginOfJob *  )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 354 of file TkAccumulatingSensitiveDetector.cc.

References edm::EventSetup::get(), numberingScheme(), and theNumberingScheme.

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().

                                                                {
  edm::ESHandle<GeometricDet> pDD;
  const edm::EventSetup* es = (*i)();
  es->get<IdealGeometryRecord>().get(pDD);

  theNumberingScheme = &(numberingScheme(*pDD));
}

void TkAccumulatingSensitiveDetector::updateHit ( const G4Step *  aStep )

private

Definition at line 298 of file TkAccumulatingSensitiveDetector.cc.

References UpdatablePSimHit::addEnergyLoss(), SensitiveDetector::ConvertToLocal3DPoint(), PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), GeV, globalExitPoint, SensitiveDetector::LocalPostStepPosition(), LogDebug, mySimHit, printHits, UpdatablePSimHit::setExitPoint(), theRotation, and PSimHit::trackId().

Referenced by ProcessHits().

                                                                   {
  // VI: in past here was a check if a hit is inside a sensitive detector,
  //     this is not needed, because call to senstive detector happens
  //     only inside the volume
  Local3DPoint theExitPoint = theRotation.get()->transformPoint(LocalPostStepPosition(aStep));
  float theEnergyLoss = aStep->GetTotalEnergyDeposit() / GeV;
  mySimHit->setExitPoint(theExitPoint);
  mySimHit->addEnergyLoss(theEnergyLoss);
  if (printHits) {
    globalExitPoint = ConvertToLocal3DPoint(aStep->GetPostStepPoint()->GetPosition());
    LogDebug("TrackerSimDebug") << " updateHit: for " << aStep->GetTrack()->GetDefinition()->GetParticleName()
                                << " trackID= " << aStep->GetTrack()->GetTrackID() << " deltaEloss= " << theEnergyLoss
                                << "\n Updated PSimHit: " << mySimHit->detUnitId() << " " << mySimHit->trackId() << " "
                                << mySimHit->energyLoss() << " " << mySimHit->entryPoint() << " "
                                << mySimHit->exitPoint();
  }
}

Member Data Documentation

bool TkAccumulatingSensitiveDetector::allowZeroEnergyLoss

private

Definition at line 66 of file TkAccumulatingSensitiveDetector.h.

Referenced by ProcessHits(), and TkAccumulatingSensitiveDetector().

float TkAccumulatingSensitiveDetector::energyCut

private

Definition at line 73 of file TkAccumulatingSensitiveDetector.h.

Referenced by TkAccumulatingSensitiveDetector(), and update().

float TkAccumulatingSensitiveDetector::energyHistoryCut

private

Definition at line 74 of file TkAccumulatingSensitiveDetector.h.

Referenced by TkAccumulatingSensitiveDetector(), and update().

int TkAccumulatingSensitiveDetector::eventno

private

Definition at line 86 of file TkAccumulatingSensitiveDetector.h.

Referenced by sendHit(), and update().

Local3DPoint TkAccumulatingSensitiveDetector::globalEntryPoint

private

Definition at line 82 of file TkAccumulatingSensitiveDetector.h.

Referenced by createHit(), and sendHit().

Local3DPoint TkAccumulatingSensitiveDetector::globalExitPoint

private

Definition at line 83 of file TkAccumulatingSensitiveDetector.h.

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

uint32_t TkAccumulatingSensitiveDetector::lastId

private

Definition at line 78 of file TkAccumulatingSensitiveDetector.h.

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

int TkAccumulatingSensitiveDetector::lastTrack

private

Definition at line 79 of file TkAccumulatingSensitiveDetector.h.

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

UpdatablePSimHit* TkAccumulatingSensitiveDetector::mySimHit

private

Definition at line 77 of file TkAccumulatingSensitiveDetector.h.

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

bool TkAccumulatingSensitiveDetector::neverAccumulate

private

Definition at line 68 of file TkAccumulatingSensitiveDetector.h.

Referenced by ProcessHits(), and TkAccumulatingSensitiveDetector().

const G4VPhysicalVolume* TkAccumulatingSensitiveDetector::oldVolume

private

Definition at line 84 of file TkAccumulatingSensitiveDetector.h.

Referenced by createHit(), and sendHit().

std::string TkAccumulatingSensitiveDetector::pname

private

Definition at line 87 of file TkAccumulatingSensitiveDetector.h.

Referenced by createHit(), and sendHit().

bool TkAccumulatingSensitiveDetector::printHits

private

Definition at line 67 of file TkAccumulatingSensitiveDetector.h.

Referenced by ntuplePrintersDiff.TrackingParticlePrinter::__call__(), createHit(), ntuplePrintersDiff.SeedPrinter::diff(), ntuplePrintersDiff.TrackPrinter::diff(), ntuplePrintersDiff.TrackingParticlePrinter::diff(), ntuplePrintersDiff.SeedPrinter::printSeed(), ntuplePrintersDiff.TrackPrinter::printTrack(), sendHit(), TkAccumulatingSensitiveDetector(), and updateHit().

float TkAccumulatingSensitiveDetector::px

private

Definition at line 85 of file TkAccumulatingSensitiveDetector.h.

Referenced by createHit(), and sendHit().

float TkAccumulatingSensitiveDetector::py

private

Definition at line 85 of file TkAccumulatingSensitiveDetector.h.

Referenced by createHit(), and sendHit().

float TkAccumulatingSensitiveDetector::pz

private

Definition at line 85 of file TkAccumulatingSensitiveDetector.h.

Referenced by createHit(), and sendHit().

double TkAccumulatingSensitiveDetector::rTracker

private

Definition at line 70 of file TkAccumulatingSensitiveDetector.h.

Referenced by TkAccumulatingSensitiveDetector().

double TkAccumulatingSensitiveDetector::rTracker2

private

Definition at line 69 of file TkAccumulatingSensitiveDetector.h.

Referenced by TkAccumulatingSensitiveDetector(), and update().

std::unique_ptr<TrackingSlaveSD> TkAccumulatingSensitiveDetector::slaveHighTof

private

Definition at line 62 of file TkAccumulatingSensitiveDetector.h.

Referenced by clearHits(), fillHits(), sendHit(), and TkAccumulatingSensitiveDetector().

std::unique_ptr<TrackingSlaveSD> TkAccumulatingSensitiveDetector::slaveLowTof

private

Definition at line 61 of file TkAccumulatingSensitiveDetector.h.

Referenced by clearHits(), fillHits(), sendHit(), and TkAccumulatingSensitiveDetector().

std::unique_ptr<const G4ProcessTypeEnumerator> TkAccumulatingSensitiveDetector::theG4ProcTypeEnumerator

private

Definition at line 64 of file TkAccumulatingSensitiveDetector.h.

Referenced by createHit(), and TkAccumulatingSensitiveDetector().

const SimTrackManager* TkAccumulatingSensitiveDetector::theManager

private

Definition at line 60 of file TkAccumulatingSensitiveDetector.h.

TrackerG4SimHitNumberingScheme* TkAccumulatingSensitiveDetector::theNumberingScheme

private

Definition at line 65 of file TkAccumulatingSensitiveDetector.h.

Referenced by setDetUnitId(), TkAccumulatingSensitiveDetector(), and update().

std::unique_ptr<FrameRotation> TkAccumulatingSensitiveDetector::theRotation

private

Definition at line 63 of file TkAccumulatingSensitiveDetector.h.

Referenced by closeHit(), createHit(), TkAccumulatingSensitiveDetector(), and updateHit().

float TkAccumulatingSensitiveDetector::theTofLimit

private

Definition at line 72 of file TkAccumulatingSensitiveDetector.h.

Referenced by sendHit(), and TkAccumulatingSensitiveDetector().

double TkAccumulatingSensitiveDetector::zTracker

private

Definition at line 71 of file TkAccumulatingSensitiveDetector.h.

Referenced by TkAccumulatingSensitiveDetector(), and update().