PltSD Class Reference

#include <PltSD.h>

Inheritance diagram for PltSD:

Public Member Functions
	PltSD (const std::string &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
uint32_t	setDetUnitId (const G4Step *) override
	~PltSD () override
Public Member Functions inherited from TimingSD
void	clearHits () override
void	EndOfEvent (G4HCofThisEvent *eventHC) override
void	fillHits (edm::PSimHitContainer &, const std::string &) override
void	Initialize (G4HCofThisEvent *HCE) override
void	PrintAll () override
bool	ProcessHits (G4Step *, G4TouchableHistory *) override
	TimingSD (const std::string &, const SensitiveDetectorCatalog &, const SimTrackManager *)
	~TimingSD () override
Public Member Functions inherited from SensitiveTkDetector
	SensitiveTkDetector (const std::string &iname, const SensitiveDetectorCatalog &clg)
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 SensitiveDetectorCatalog &, bool calo, const std::string &newcollname="")
	~SensitiveDetector () override
Public Member Functions inherited from Observer< const BeginOfEvent *>
	Observer ()
void	slotForUpdate (const BeginOfEvent * iT)
virtual	~Observer ()

Protected Member Functions
bool	checkHit (const G4Step *, BscG4Hit *) override
Protected Member Functions inherited from TimingSD
const G4ThreeVector &	getGlobalEntryPoint () const
const G4ThreeVector &	getLocalEntryPoint () const
virtual int	getTrackID (const G4Track *)
void	setCuts (double eCut, double historyCut)
void	setTimeFactor (double)
void	setToLocal (const G4StepPoint *stepPoint, const G4ThreeVector &globalPoint, G4ThreeVector &localPoint)
void	update (const BeginOfEvent *) 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 Attributes
double	energyCut
double	energyHistoryCut

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

Detailed Description

Definition at line 13 of file PltSD.h.

Constructor & Destructor Documentation

PltSD()

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

Definition at line 19 of file PltSD.cc.

References energyCut, energyHistoryCut, edm::ParameterSet::getParameter(), AlCaHLTBitMon_ParallelJobs::p, and TimingSD::setCuts().

    : TimingSD(name, clg, manager) {
  edm::ParameterSet m_TrackerSD = p.getParameter<edm::ParameterSet>("PltSD");
  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

  setCuts(energyCut, energyHistoryCut);
}

~PltSD()

PltSD::~PltSD ( )

override

Definition at line 33 of file PltSD.cc.

{}

Member Function Documentation

checkHit()

bool PltSD::checkHit ( const G4Step *  , BscG4Hit *  hit  )

overrideprotectedvirtual

Reimplemented from TimingSD.

Definition at line 126 of file PltSD.cc.

References nano_mu_digi_cff::float, TimingSD::getLocalEntryPoint(), and mag2().

                                                 {
  // 50 micron are allowed between the exit
  // point of the current hit and the entry point of the new hit
  static const float tolerance2 = (float)(0.0025 * CLHEP::mm * CLHEP::mm);
  return ((hit->getExitLocalP() - getLocalEntryPoint()).mag2() < tolerance2);
}

setDetUnitId()

uint32_t PltSD::setDetUnitId ( const G4Step *  aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 35 of file PltSD.cc.

References hcalRecHitTable_cff::detId, ForwardName::getName(), personalPlayback::level, AlCaHLTBitMon_QueryRunRegistry::string, and groupFilesInBlocks::temp.

                                                {
  unsigned int detId = 0;

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("PltSD") << " DetID = " << detId;
#endif
  //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
    std::string sensorName = ForwardName::getName(touch->GetVolume(2)->GetName());
    std::string telName = ForwardName::getName(touch->GetVolume(3)->GetName());
    std::string volumeName = ForwardName::getName(touch->GetVolume(4)->GetName());
    if (sensorName != "PLTSensorPlane")
      edm::LogVerbatim("PltSD") << " PltSD::setDetUnitId -w- Sensor name not PLTSensorPlane ";
    if (telName != "Telescope")
      edm::LogVerbatim("PltSD") << " PltSD::setDetUnitId -w- Telescope name not Telescope ";
    if (volumeName != "PLT")
      edm::LogVerbatim("PltSD") << " PltSD::setDetUnitId -w- Volume name not PLT ";

    //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;
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("PltSD") << "Hit Recorded at "
                              << "plt:" << pltNum << " hc:" << halfCarriageNum << " tel:" << telNum
                              << " plane:" << sensorNum;
#endif
  }
  return detId;
}

Member Data Documentation

energyCut

double PltSD::energyCut

private

Definition at line 24 of file PltSD.h.

Referenced by PltSD().

energyHistoryCut

double PltSD::energyHistoryCut

private

Definition at line 25 of file PltSD.h.

Referenced by PltSD().