#include <HFWedgeSD.h>

Inheritance diagram for HFWedgeSD:

Public Member Functions
virtual void	clear ()

virtual void	DrawAll ()

virtual void	EndOfEvent (G4HCofThisEvent *eventHC)

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

virtual void	Initialize (G4HCofThisEvent *HCE)

virtual void	PrintAll ()

virtual bool	ProcessHits (G4Step step, G4TouchableHistory tHistory)

virtual	~HFWedgeSD ()

Public Member Functions inherited from SensitiveCaloDetector
	SensitiveCaloDetector (std::string &iname, const DDCompactView &cpv, SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)

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

Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point)

Local3DPoint	FinalStepPosition (G4Step *s, coordinates)

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

Local3DPoint	InitialStepPosition (G4Step *s, coordinates)

std::string	nameOfSD ()

void	NaNTrap (G4Step *step)

void	Register ()

	SensitiveDetector (std::string &iname, const DDCompactView &cpv, SensitiveDetectorCatalog &, edm::ParameterSet const &p)

virtual	~SensitiveDetector ()

Protected Member Functions
virtual void	clearHits ()

HFShowerG4Hit *	createNewHit ()

virtual void	fillHits (edm::PCaloHitContainer &, std::string)

G4bool	hitExists ()

virtual uint32_t	setDetUnitId (G4Step *)

void	updateHit (HFShowerG4Hit *)

Private Attributes
HFShowerG4Hit *	currentHit

int	currentID

double	edep

G4ThreeVector	globalPos

int	hcID

std::map< int, HFShowerG4Hit * >	hitMap

G4ThreeVector	localPos

const SimTrackManager *	m_trackManager

G4ThreeVector	momDir

int	previousID

HFShowerG4HitsCollection *	theHC

std::string	theName

double	time

int	trackID

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

Detailed Description

Definition at line 22 of file HFWedgeSD.h.

Constructor & Destructor Documentation

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

Definition at line 19 of file HFWedgeSD.cc.

References SensitiveDetector::AssignSD(), ecaldqm::collectionName, LogDebug, SensitiveDetectorCatalog::logicalNames(), and SensitiveDetector::Register().

                                          :
   SensitiveCaloDetector(name, cpv, clg, p), theName(name),
   m_trackManager(manager), hcID(-1), theHC(0), currentHit(0) {
 
   collectionName.insert(name);
   LogDebug("FiberSim") << "***************************************************"
                << "\n"
                << "*                                                 *"
                << "\n"
                << "* Constructing a HFWedgeSD  with name " << GetName()
                << "\n"
                << "*                                                 *"
                << "\n"
                << "***************************************************";
   //
   // Now attach the right detectors (LogicalVolumes) to me
   //
   std::vector<std::string> lvNames = clg.logicalNames(name);
   this->Register();
   for (std::vector<std::string>::iterator it=lvNames.begin();
        it !=lvNames.end(); it++){
     this->AssignSD(*it);
     LogDebug("FiberSim") << "HFWedgeSD : Assigns SD to LV " << (*it);
   }
 }

HFWedgeSD::~HFWedgeSD ( )

virtual

Definition at line 47 of file HFWedgeSD.cc.

References theHC.

                       {
   if (theHC)    delete theHC;
 }

Member Function Documentation

void HFWedgeSD::clear ( void )

virtual

Definition at line 88 of file HFWedgeSD.cc.

Referenced by EndOfEvent().

88 {}

void HFWedgeSD::clearHits ( )

protectedvirtual

Implements SensitiveDetector.

Definition at line 143 of file HFWedgeSD.cc.

References hitMap, and previousID.

Referenced by Initialize().

                           {
   hitMap.erase (hitMap.begin(), hitMap.end());
   previousID = -1;
 }

HFShowerG4Hit * HFWedgeSD::createNewHit ( )

protected

Definition at line 111 of file HFWedgeSD.cc.

References currentID, edep, globalPos, hitMap, localPos, LogDebug, MeV, momDir, previousID, HFShowerG4Hit::setGlobalPos(), HFShowerG4Hit::setHitId(), HFShowerG4Hit::setLocalPos(), HFShowerG4Hit::setPrimMomDir(), HFShowerG4Hit::setTime(), HFShowerG4Hit::setTrackId(), theHC, time, trackID, and updateHit().

Referenced by ProcessHits().

                                        {
 
   LogDebug("FiberSim") << "HFWedgeSD::CreateNewHit for ID " << currentID
                << " Track " << trackID << " Edep: " << edep/MeV 
                << " MeV; Time: " << time << " ns; Position (local) " 
                << localPos << " (global ) " << globalPos 
                << " direction " << momDir;
   HFShowerG4Hit* aHit = new HFShowerG4Hit;
   aHit->setHitId(currentID);
   aHit->setTrackId(trackID);
   aHit->setTime(time);
   aHit->setLocalPos(localPos);
   aHit->setGlobalPos(globalPos);
   aHit->setPrimMomDir(momDir);
   updateHit(aHit);
 
   theHC->insert(aHit);
   hitMap.insert(std::pair<int,HFShowerG4Hit*>(previousID,aHit));
 
   return aHit;
 }    

void HFWedgeSD::DrawAll ( )

virtual

Definition at line 90 of file HFWedgeSD.cc.

90 {}

void HFWedgeSD::EndOfEvent ( G4HCofThisEvent * eventHC )

virtual

Reimplemented from SensitiveDetector.

Definition at line 82 of file HFWedgeSD.cc.

References clear(), LogDebug, and theHC.

                                                 {
  
   LogDebug("FiberSim") << "HFWedgeSD: Sees" << theHC->entries() << " hits";
   clear();
 }

void HFWedgeSD::fillHits	(	edm::PCaloHitContainer &	,
		std::string
	)

protectedvirtual

Implements SensitiveCaloDetector.

Definition at line 153 of file HFWedgeSD.cc.

153 {}

G4bool HFWedgeSD::hitExists ( )

protected

Definition at line 94 of file HFWedgeSD.cc.

References currentHit, currentID, hitMap, previousID, and updateHit().

Referenced by ProcessHits().

                             {
    
   // Update if in the same detector, time-slice and for same track   
   if (currentID == previousID) {
     updateHit(currentHit);
     return true;
   }
    
   std::map<int,HFShowerG4Hit*>::const_iterator it = hitMap.find(currentID);
   if (it != hitMap.end()) {
     updateHit(currentHit);
     return true;
   }
   
   return false;
 }

void HFWedgeSD::Initialize ( G4HCofThisEvent * HCE )

virtual

Reimplemented from SensitiveDetector.

Definition at line 51 of file HFWedgeSD.cc.

References clearHits(), ecaldqm::collectionName, hcID, LogDebug, and theHC.

                                                 {
 
   LogDebug("FiberSim") << "HFWedgeSD : Initialize called for " << GetName();
   theHC = new HFShowerG4HitsCollection(GetName(), collectionName[0]);
   if (hcID<0)
     hcID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
   HCE->AddHitsCollection(hcID, theHC);
 
   clearHits();
 }

void HFWedgeSD::PrintAll ( )

virtual

Definition at line 92 of file HFWedgeSD.cc.

92 {}

G4bool HFWedgeSD::ProcessHits	(	G4Step *	step,
		G4TouchableHistory *	tHistory
	)

virtual

Implements SensitiveDetector.

Definition at line 62 of file HFWedgeSD.cc.

References createNewHit(), currentHit, currentID, edep, globalPos, hitExists(), localPos, momDir, setDetUnitId(), time, and trackID.

                                                                  {
 
   G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
   const G4VTouchable* touch = preStepPoint->GetTouchable();
   currentID = setDetUnitId(aStep);
   trackID   = aStep->GetTrack()->GetTrackID();
   edep      = aStep->GetTotalEnergyDeposit();
   time      = (preStepPoint->GetGlobalTime())/ns;
 
   globalPos = preStepPoint->GetPosition();
   localPos  = touch->GetHistory()->GetTopTransform().TransformPoint(globalPos);
   const G4DynamicParticle* particle =  aStep->GetTrack()->GetDynamicParticle();
   momDir    = particle->GetMomentumDirection();
 
   if (hitExists() == false && edep>0.) 
     currentHit = createNewHit();
 
   return true;
 }

uint32_t HFWedgeSD::setDetUnitId ( G4Step * aStep )

protectedvirtual

Implements SensitiveDetector.

Definition at line 148 of file HFWedgeSD.cc.

Referenced by ProcessHits().

                                               {
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   return (touch->GetReplicaNumber(0));
 }

void HFWedgeSD::updateHit ( HFShowerG4Hit * aHit )

protected

Definition at line 133 of file HFWedgeSD.cc.

References currentID, edep, LogDebug, MeV, previousID, and HFShowerG4Hit::updateEnergy().

Referenced by createNewHit(), and hitExists().

                                              {
 
   if (edep != 0) {
     aHit->updateEnergy(edep);
     LogDebug("FiberSim") << "HFWedgeSD: Add energy deposit in " << currentID 
              << " edep " << edep/MeV << " MeV"; 
   }
   previousID = currentID;
 }