#include <HGCalTB16SD01.h>

Inheritance diagram for HGCalTB16SD01:

Public Member Functions
double	getEnergyDeposit (G4Step *) override

	HGCalTB16SD01 (const std::string &, const DDCompactView &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)

uint32_t	setDetUnitId (const G4Step *step) override

	~HGCalTB16SD01 () override

Public Member Functions inherited from CaloSD
	CaloSD (const std::string &aSDname, const DDCompactView &cpv, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, float timeSlice=1., bool ignoreTkID=false)

void	clear () override

void	clearHits () override

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

void	fillHits (edm::PCaloHitContainer &, const std::string &) override

void	Initialize (G4HCofThisEvent *HCE) override

void	PrintAll () override

bool	ProcessHits (G4Step step, G4TouchableHistory tHistory) override

bool	ProcessHits (G4GFlashSpot aSpot, G4TouchableHistory ) override

	~CaloSD () override

Public Member Functions inherited from SensitiveCaloDetector
	SensitiveCaloDetector (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 BeginOfRun * >
	Observer ()

void	slotForUpdate (const BeginOfRun * iT)

virtual	~Observer ()

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 EndOfTrack * >
	Observer ()

void	slotForUpdate (const EndOfTrack * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfEvent * >
	Observer ()

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

Static Public Member Functions
static uint32_t	packIndex (int det, int lay, int x, int y)

static void	unpackIndex (const uint32_t &idx, int &det, int &lay, int &x, int &y)

Private Member Functions
void	initialize (const G4StepPoint *point)

Private Attributes
double	birk1_

double	birk2_

double	birk3_

bool	initialize_

std::string	matName_

G4Material *	matScin_

bool	useBirk_

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

Protected Member Functions inherited from CaloSD
G4bool	checkHit ()

CaloG4Hit *	createNewHit ()

virtual bool	filterHit (CaloG4Hit *, double)

double	getAttenuation (const G4Step *aStep, double birk1, double birk2, double birk3)

virtual uint16_t	getDepth (const G4Step *)

int	getNumberOfHits ()

double	getResponseWt (const G4Track *)

virtual G4bool	getStepInfo (G4Step *aStep)

virtual int	getTrackID (const G4Track *)

G4bool	hitExists ()

virtual void	initRun ()

void	resetForNewPrimary (const G4ThreeVector &, double)

G4ThreeVector	setToGlobal (const G4ThreeVector &, const G4VTouchable *)

G4ThreeVector	setToLocal (const G4ThreeVector &, const G4VTouchable *)

void	update (const BeginOfRun *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfEvent *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfTrack *trk) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const EndOfTrack *trk) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const ::EndOfEvent *) override

void	updateHit (CaloG4Hit *)

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 > &)

Protected Member Functions inherited from Observer< const EndOfEvent * >
virtual void	update (const EndOfEvent *)=0
	This routine will be called when the appropriate signal arrives. More...

Protected Attributes inherited from CaloSD
int	checkHits

double	correctT

bool	corrTOFBeam

CaloG4Hit *	currentHit

CaloHitID	currentID

float	edepositEM

float	edepositHAD

double	eminHit

double	eminHitD

G4int	emPDG

double	energyCut

G4ThreeVector	entranceLocal

G4ThreeVector	entrancePoint

G4int	epPDG

bool	forceSave

G4int	gammaPDG

float	incidentEnergy

double	kmaxIon

double	kmaxNeutron

double	kmaxProton

const SimTrackManager *	m_trackManager

G4ThreeVector	posGlobal

G4StepPoint *	preStepPoint

CaloHitID	previousID

int	primIDSaved

bool	runInit

bool	suppressHeavy

G4Track *	theTrack

double	tmaxHit

bool	useMap

Detailed Description

Definition at line 19 of file HGCalTB16SD01.h.

Constructor & Destructor Documentation

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

Definition at line 24 of file HGCalTB16SD01.cc.

References birk1_, birk2_, birk3_, g, edm::ParameterSet::getParameter(), matName_, matScin_, MeV, AlCaHLTBitMon_QueryRunRegistry::string, and useBirk_.

                                                  : 
   CaloSD(name, cpv, clg, p, manager), initialize_(true) {
   
   // Values from NIM 80 (1970) 239-244: as implemented in Geant3
   edm::ParameterSet m_HC = p.getParameter<edm::ParameterSet>("HGCalTestBeamSD");
   matName_   = m_HC.getParameter<std::string>("Material");
   useBirk_   = m_HC.getParameter<bool>("UseBirkLaw");
   birk1_     = m_HC.getParameter<double>("BirkC1")*(g/(MeV*cm2));
   birk2_     = m_HC.getParameter<double>("BirkC2");
   birk3_     = m_HC.getParameter<double>("BirkC3");
   matScin_   = nullptr;
 
   edm::LogInfo("HGCSim") << "HGCalTB16SD01:: Use of Birks law is set to " 
              << useBirk_ << " for " << matName_ 
              << " with three constants kB = " << birk1_ 
              << ", C1 = " << birk2_ << ", C2 = " << birk3_;
 }

HGCalTB16SD01::~HGCalTB16SD01 ( )

override

Definition at line 45 of file HGCalTB16SD01.cc.

45 {}

Member Function Documentation

double HGCalTB16SD01::getEnergyDeposit ( G4Step * aStep )

overridevirtual

Reimplemented from CaloSD.

Definition at line 47 of file HGCalTB16SD01.cc.

References birk1_, birk2_, birk3_, gather_cfg::cout, CaloSD::getAttenuation(), initialize(), initialize_, matScin_, point, and useBirk_.

                                                     {
 
   G4StepPoint* point = aStep->GetPreStepPoint();
   if (initialize_) initialize(point);
   double destep = aStep->GetTotalEnergyDeposit();
   double wt2    = aStep->GetTrack()->GetWeight();
   double weight = (wt2 > 0.0) ? wt2 : 1.0;
   if (useBirk_ && matScin_ == point->GetMaterial()) {
     weight *= getAttenuation(aStep, birk1_, birk2_, birk3_);
   }
 #ifdef EDM_ML_DEBUG
   std::cout << "HGCalTB16SD01: Detector " 
         << point->GetTouchable()->GetVolume()->GetName() << " with "
         << point->GetMaterial()->GetName() << " weight " << weight 
         << ":" << wt2 << std::endl;
 #endif
   return weight*destep;
 }

void HGCalTB16SD01::initialize ( const G4StepPoint * point )

private

Definition at line 109 of file HGCalTB16SD01.cc.

References gather_cfg::cout, initialize_, matName_, and matScin_.

Referenced by getEnergyDeposit().

                                                        {
   if (matName_ == point->GetMaterial()->GetName()) {
     matScin_    =  point->GetMaterial();
     initialize_ = false;
   }
 #ifdef EDM_ML_DEBUG
   std::cout << "HGCalTB16SD01: Material pointer for " << matName_
         << " is initialized to : " << matScin_ << std::endl;
 #endif
 }

uint32_t HGCalTB16SD01::packIndex	(	int	det,
		int	lay,
		int	x,
		int	y
	)

static

Definition at line 78 of file HGCalTB16SD01.cc.

References gather_cfg::cout, TauDecayModes::dec, training_settings::idx, ixx, iyy, x, and y.

Referenced by setDetUnitId().

                                                                 {
 
   int ix = 0, ixx = x;
   if (x < 0) { ix = 1; ixx =-x;}
   int iy = 0, iyy = y;
   if (y < 0) { iy = 1; iyy =-y;}
   uint32_t idx = (det&15)<<28;      //bits 28-31
   idx         += (lay&127)<<21;     //bits 21-27
   idx         += (iy&1)<<19;        //bit  19
   idx         += (iyy&511)<<10;     //bits 10-18
   idx         += (ix&1)<<9;         //bit   9
   idx         += (ixx&511);         //bits  0-8
 
 #ifdef EDM_ML_DEBUG
   std::cout << "HGCalTB16SD01: Detector " << det << " Layer "  << lay << " x "
         << x << " " << ix << " " << ixx << " y " << y << " " << iy << " " 
         << iyy << " ID " << std::hex << idx << std::dec << std::endl;
 #endif
   return idx;
 }

uint32_t HGCalTB16SD01::setDetUnitId ( const G4Step * step )

overridevirtual

Implements CaloSD.

Definition at line 66 of file HGCalTB16SD01.cc.

References dataset::name, packIndex(), CaloSD::preStepPoint, x, and y.

                                                          { 
 
   G4StepPoint* preStepPoint = aStep->GetPreStepPoint(); 
   const G4VTouchable* touch = preStepPoint->GetTouchable();
   G4String name             = preStepPoint->GetPhysicalVolume()->GetName();
 
   int det(1), x(0), y(0);
   int lay = (touch->GetReplicaNumber(0));
 
   return packIndex (det, lay, x, y);
 }

void HGCalTB16SD01::unpackIndex	(	const uint32_t &	idx,
		int &	det,
		int &	lay,
		int &	x,
		int &	y
	)

static

Definition at line 99 of file HGCalTB16SD01.cc.

References training_settings::idx, x, and y.

                                  {
 
   det  = (idx>>28)&15;
   lay  = (idx>>21)&127;
   y    = (idx>>10)&511; if (((idx>>19)&1) == 1) y = -y;
   x    = (idx)&511;     if (((idx>>9)&1)  == 1) x = -x;
 
 }