#include <PltSD.h>

Inheritance diagram for PltSD:

Public Member Functions
	PltSD (const std::string &, const DDCompactView &, 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 DDCompactView &, const SensitiveDetectorCatalog &, const edm::ParameterSet &, const SimTrackManager *)

	~TimingSD () 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 ()

Protected Member Functions
bool	checkHit (const G4Step , BscG4Hit ) override

Protected Member Functions inherited from TimingSD
const G4ThreeVector &	getGlobalEntryPoint () const

const G4ThreeVector &	getLocalEntryPoint () const

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 11 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 18 of file PltSD.cc.

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

   : TimingSD(name, cpv, clg, p, 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 ( )

override

Definition at line 30 of file PltSD.cc.

30 {

31 }

Member Function Documentation

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

overrideprotectedvirtual

Reimplemented from TimingSD.

Definition at line 126 of file PltSD.cc.

References objects.autophobj::float, BscG4Hit::getExitLocalP(), 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);
 }

uint32_t PltSD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 33 of file PltSD.cc.

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

                                                   {
     
     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;
 }

Member Data Documentation

double PltSD::energyCut

private

Definition at line 28 of file PltSD.h.

Referenced by PltSD().

double PltSD::energyHistoryCut

private

Definition at line 29 of file PltSD.h.

Referenced by PltSD().

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation