#include <AHCalSD.h>

Inheritance diagram for AHCalSD:

Public Member Functions
	AHCalSD (G4String, const DDCompactView &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)

double	getEnergyDeposit (G4Step *) override

uint32_t	setDetUnitId (G4Step *step) override

bool	unpackIndex (const uint32_t &idx, int &row, int &col, int &depth)

	~AHCalSD () override

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

void	clear () override

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

void	fillHits (edm::PCaloHitContainer &, std::string n) 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 (std::string &iname, const DDCompactView &cpv, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)

Public Member Functions inherited from SensitiveDetector
virtual void	AssignSD (const std::string &vname)

Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point)

void	EndOfEvent (G4HCofThisEvent *eventHC) override

Local3DPoint	FinalStepPosition (G4Step *s, coordinates)

virtual std::vector< std::string >	getNames ()

void	Initialize (G4HCofThisEvent *eventHC) override

Local3DPoint	InitialStepPosition (G4Step *s, coordinates)

std::string	nameOfSD ()

void	NaNTrap (G4Step *step)

void	Register ()

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

Protected Member Functions
bool	filterHit (CaloG4Hit *, double) override

Protected Member Functions inherited from CaloSD
G4bool	checkHit ()

void	clearHits () override

CaloG4Hit *	createNewHit ()

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

virtual uint16_t	getDepth (G4Step *)

int	getNumberOfHits ()

double	getResponseWt (G4Track *)

virtual G4bool	getStepInfo (G4Step *aStep)

virtual int	getTrackID (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 Observer< const EndOfEvent * >
virtual void	update (const EndOfEvent *)=0
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
double	betaThr

double	birk1

double	birk2

double	birk3

double	eminHit

bool	useBirk

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

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 12 of file AHCalSD.h.

Constructor & Destructor Documentation

AHCalSD::AHCalSD	(	G4String	name,
		const DDCompactView &	cpv,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

Definition at line 20 of file AHCalSD.cc.

References birk1, birk2, birk3, eminHit, g, edm::ParameterSet::getParameter(), MeV, and useBirk.

                                                                     : 
   CaloSD(name, cpv, clg, p, manager,
          (float)(p.getParameter<edm::ParameterSet>("AHCalSD").getParameter<double>("TimeSliceUnit")),
          p.getParameter<edm::ParameterSet>("AHCalSD").getParameter<bool>("IgnoreTrackID")) {
 
   edm::ParameterSet m_HC = p.getParameter<edm::ParameterSet>("AHCalSD");
   useBirk          = m_HC.getParameter<bool>("UseBirkLaw");
   birk1            = m_HC.getParameter<double>("BirkC1")*(CLHEP::g/(CLHEP::MeV*CLHEP::cm2));
   birk2            = m_HC.getParameter<double>("BirkC2");
   birk3            = m_HC.getParameter<double>("BirkC3");
   eminHit          = m_HC.getParameter<double>("EminHit")*CLHEP::MeV;
 
   edm::LogInfo("HcalSim") << "AHCalSD::  Use of Birks law is set to      " 
                           << useBirk << "  with three constants kB = "
                           << birk1 << ", C1 = " << birk2 << ", C2 = " << birk3
               << "\nAHCalSD:: Threshold for storing hits: "
                           << eminHit << std::endl;
 }

AHCalSD::~AHCalSD ( )

override

Definition at line 41 of file AHCalSD.cc.

41 { }

Member Function Documentation

bool AHCalSD::filterHit	(	CaloG4Hit *	aHit,
		double	time
	)

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 106 of file AHCalSD.cc.

References eminHit, CaloG4Hit::getEnergyDeposit(), and CaloSD::tmaxHit.

                                                     {
   return ((time <= tmaxHit) && (aHit->getEnergyDeposit() > eminHit));
 }

double AHCalSD::getEnergyDeposit ( G4Step * aStep )

overridevirtual

Reimplemented from CaloSD.

Definition at line 43 of file AHCalSD.cc.

References birk1, birk2, birk3, CaloSD::getAttenuation(), useBirk, and mps_merge::weight.

                                               {
 
   double destep = aStep->GetTotalEnergyDeposit();
   double wt2    = aStep->GetTrack()->GetWeight();
   double weight = (wt2 > 0.0) ? wt2 : 1.0;
 #ifdef EDM_ML_DEBUG
   double weight0 = weight;
 #endif
   if (useBirk) weight *= getAttenuation(aStep, birk1, birk2, birk3);
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HcalSim") << "AHCalSD: weight " << weight0 << " " << weight 
               << std::endl;
 #endif
   double edep = weight*destep;
   return edep;
 }

uint32_t AHCalSD::setDetUnitId ( G4Step * step )

overridevirtual

Implements CaloSD.

Definition at line 60 of file AHCalSD.cc.

References cuy::col, TauDecayModes::dec, particleFlowClusterECALTimeSelected_cfi::depth, RemoveAddSevLevel::flag, HcalOther, diffTreeTool::index, CaloSD::preStepPoint, DetId::rawId(), and unpackIndex().

                                              { 
 
   G4StepPoint* preStepPoint = aStep->GetPreStepPoint(); 
   const G4VTouchable* touch = preStepPoint->GetTouchable();
 
   int depth = (touch->GetReplicaNumber(1));
   int incol = ((touch->GetReplicaNumber(0))%10);
   int inrow = ((touch->GetReplicaNumber(0))/10)%10;
   int jncol = ((touch->GetReplicaNumber(0))/100)%10;
   int jnrow = ((touch->GetReplicaNumber(0))/1000)%10;
   int col   = (jncol == 0) ? incol : -incol;
   int row   = (jnrow == 0) ? inrow : -inrow;
   uint32_t index = AHCalDetId(row,col,depth).rawId();
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HcalSim") << "AHCalSD: det = " << HcalOther 
                           << " depth = " << depth << " row = " << row 
                           << " column = " << col << " packed index = 0x" 
               << std::hex << index << std::dec << std::endl;
   bool flag = unpackIndex(index, row, col, depth);
   edm::LogInfo("HcalSim") << "Results from unpacker for 0x" << std::hex 
               << index << std::dec << " Flag " << flag << " Row " 
               << row << " Col " << col << " Depth " << depth 
               << std::endl;
 #endif
   return index;
 }

bool AHCalSD::unpackIndex	(	const uint32_t &	idx,
		int &	row,
		int &	col,
		int &	depth
	)

Definition at line 87 of file AHCalSD.cc.

References TauDecayModes::dec, AHCalDetId::depth(), DetId::det(), relval_steps::gen(), DetId::Hcal, HcalOther, AHCalDetId::icol(), AHCalDetId::irow(), and DetId::subdetId().

Referenced by setDetUnitId().

                                                                              {
 
   DetId gen(idx);
   HcalSubdetector subdet = (HcalSubdetector(gen.subdetId()));
   bool rcode = (gen.det()==DetId::Hcal && subdet!=HcalOther);
   row = col = depth = 0;
   if (rcode) {
     row   = AHCalDetId(idx).irow();
     col   = AHCalDetId(idx).icol();
     depth = AHCalDetId(idx).depth();
   }
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HcalSim") << "AHCalSD: packed index = 0x" << std::hex << idx 
               << std::dec << " Row " << row << " Column " << col
               << " Depth " << depth << " OK " << rcode << std::endl;
 #endif
   return rcode;
 }