Sensitive Detector implementation of the CSI detector. More...

#include <SensitiveDetector.h>

Inheritance diagram for SensitiveDetector:

Collaboration diagram for SensitiveDetector:

[legend]

Public Member Functions
	SensitiveDetector ()
	Constructor.

	~SensitiveDetector ()
	Destructor.

int	saveSimHit (const G4int cellId, const G4int trackID, const G4int pid, const G4double tof, const G4double edep, G4ThreeVector mom, G4ThreeVector pos, G4ThreeVector PosCell, G4ThreeVector VecCell)
	Save ECLSimHits of the event into datastore. More...

Static Public Member Functions
static const std::map< std::string, RelationArray::EConsolidationAction > &	getMCParticleRelations ()
	Return a list of all registered Relations with MCParticles.

static void	setActive (bool activeStatus)
	Enable/Disable all Sensitive Detectors. More...

static void	registerMCParticleRelation (const std::string &name, RelationArray::EConsolidationAction ignoreAction=RelationArray::c_negativeWeight)
	Register an relation involving MCParticles. More...

static void	registerMCParticleRelation (const RelationArray &relation, RelationArray::EConsolidationAction ignoreAction=RelationArray::c_negativeWeight)
	Overload to make it easer to register MCParticle relations. More...

Protected Member Functions
bool	step (G4Step aStep, G4TouchableHistory ) override
	Step processing method. More...

Private Member Functions
virtual bool	ProcessHits (G4Step aStep, G4TouchableHistory aROhist)
	Check if recording hits is enabled and if so call step() and set the correct MCParticle flag. More...

Private Attributes
int	m_hitNum
	members of SensitiveDetector

int	m_EvnetNumber
	The current number of created hits in an event. More...

int	m_oldEvnetNumber
	The current number of created hits in an event. More...

int	m_trackID
	The current number of created hits in an event. More...

G4ThreeVector	m_startPos
	track id

G4ThreeVector	m_endPos
	Position of prestep.

G4ThreeVector	m_WightedPos
	Position of poststep.

G4ThreeVector	m_momentum
	Wighted step Position.

double	m_startTime
	momentum of track

double	m_endTime
	global time

double	m_WightedTime
	global time

double	m_startEnergy
	global time

double	m_energyDeposit
	particle energy at the entrance in volume

double	m_trackLength
	energy deposited in volume

int	CsiHitIndex [18][80]
	length of the track in the volume

int	iECLCell
	Hit index of StoreArray.

int	TimeIndex
	Hit Energy of StoreArray.

double	local_pos
	Hit Time of StoreArray.

double	T_ave
	position alongthe vector of crystal axis

int	firstcall
	flight time to diode sensor

int	m_phiID
	flag of first call

int	m_thetaID
	The current phi ID in an event. More...

int	m_cellID
	The current theta ID in an event. More...

Const::EDetector	m_subdetector
	Subdetector the class belongs to.

Static Private Attributes
static std::map< std::string, RelationArray::EConsolidationAction >	s_mcRelations
	Static set holding all relations which have to be updated at the end of the Event.

static bool	s_active
	Static bool which indicates wether recording of hits is enabled.

Detailed Description

Sensitive Detector implementation of the CSI detector.

Definition at line 22 of file SensitiveDetector.h.

Member Function Documentation

◆ ProcessHits()

bool ProcessHits	(	G4Step *	aStep,
		G4TouchableHistory *	aROhist
	)

inlineprivatevirtualinherited

Check if recording hits is enabled and if so call step() and set the correct MCParticle flag.

Called by Geant4 for each step inside the sensitive volumes attached

Definition at line 94 of file SensitiveDetectorBase.h.

◆ registerMCParticleRelation() [1/2]

static void registerMCParticleRelation	(	const RelationArray &	relation,
		RelationArray::EConsolidationAction	ignoreAction = `RelationArray::c_negativeWeight`
	)

inlinestaticinherited

Overload to make it easer to register MCParticle relations.

Parameters

relation	RelationArray to register
ignoreAction

Definition at line 66 of file SensitiveDetectorBase.h.

◆ registerMCParticleRelation() [2/2]

void registerMCParticleRelation	(	const std::string &	name,
		RelationArray::EConsolidationAction	ignoreAction = `RelationArray::c_negativeWeight`
	)

staticinherited

Register an relation involving MCParticles.

All Relations which point from an MCParticle to something have to be registered with addMCParticleRelation() because the index of the MCParticles might change at the end of the event. During simulation, the TrackID should be used as index of the MCParticle

Parameters

name	Name of the relation to register
ignoreAction

Definition at line 22 of file SensitiveDetectorBase.cc.

◆ saveSimHit()

int saveSimHit	(	const G4int	cellId,
		const G4int	trackID,
		const G4int	pid,
		const G4double	tof,
		const G4double	edep,
		G4ThreeVector	mom,
		G4ThreeVector	pos,
		G4ThreeVector	PosCell,
		G4ThreeVector	VecCell
	)

Save ECLSimHits of the event into datastore.

Parameters

cellId	the channel identifier
trackID	the track identifier
pid	the pdg code of the particle
tof	the time of flight of the particle
edep	the energy deposited in the crystal by the particle
mom	the momentum of the particle
pos	TO BE COMPLETED
PosCell	TO BE COMPLETED
VecCell	TO BE COMPLETED

Returns: int the number of hits saved in this event

Definition at line 153 of file SensitiveDetector.cc.

     {
  
       //Get the datastore arrays
       StoreArray<MCParticle>   mcParticles;
       StoreArray<CsiSimHit>    simHits;
  
       RelationArray simHitRel(mcParticles, simHits);
  
       // convert Hep3Vectors to TVector3
       TVector3 momentum(mom.getX(), mom.getY(), mom.getZ());
       TVector3 position(pos.getX(), pos.getY(), pos.getZ());
  
       //Append
       simHits.appendNew(cellId, trackID, pid, tof / CLHEP::ns, edep / CLHEP::GeV,
                         momentum * (1.0 / CLHEP::GeV), position * (1.0 / CLHEP::cm));
  
       int simhitNumber = simHits.getEntries() - 1;
       simHitRel.add(trackID, simhitNumber);
       B2DEBUG(150, "Saving CsiSimHit Number: " << simhitNumber);
  
       // Store CsiHit: total energy deposited per hit, with average time.
  
       // Old names:
       /*
       StoreArray<CsiHit> eclHitArray;
       RelationArray eclHitRel(mcParticles, eclHitArray);
       */
  
       // New names:
  
       StoreArray<CsiHit> csiHits;
  
       RelationArray csiHitRel(mcParticles, csiHits);
  
  
       StoreObjPtr<EventMetaData> eventMetaDataPtr;
       int m_currentEvnetNumber = eventMetaDataPtr->getEvent();
  
  
       if (firstcall == 0 || m_currentEvnetNumber != m_oldEvnetNumber) {
         m_oldEvnetNumber = m_currentEvnetNumber;
         for (int iCSICell = 0; iCSICell < 18; iCSICell++) {
           for (int  TimeInd = 0; TimeInd < 80; TimeInd++) {
             CsiHitIndex[iCSICell][TimeInd] = -1;
           }
         }
         firstcall++;
       }
  
  
       if (m_currentEvnetNumber == m_oldEvnetNumber) {
         if ((tof / CLHEP::ns) < 8000) {
           TimeIndex = (int)((tof / CLHEP::ns) / 100);
  
           // local_pos is the distance between the hit and the PIN-diode end of the crystal.
           local_pos = (15. - (pos  - PosCell) * VecCell);
  
           // From the ECL package T_ave =  6.05 + 0.0749 * local_pos - 0.00112 * local_pos * local_pos + (tof / CLHEP::ns)  ;
  
           //For CsI: we remove the offset (I hypothetize it if from photo-sensor delay)
           //         we remove the acceleration of the light
           //         we correct the velocity term using n_CsI=1.80
           T_ave = 0.0600 * local_pos + (tof / CLHEP::ns)  ;
  
  
           double E_cell = (edep / CLHEP::GeV);
           if (CsiHitIndex[cellId][TimeIndex] == -1) {
             //m_hitNum = csiHits->GetLast() + 1;
             //new(csiHits->AddrAt(m_hitNum)) CSIHit();
             csiHits.appendNew();
             m_hitNum = csiHits.getEntries() - 1;
             CsiHitIndex[cellId][TimeIndex] = m_hitNum;
             csiHits[m_hitNum]->setCellId(cellId);
             csiHits[m_hitNum]->setEnergyDep(E_cell);
             csiHits[m_hitNum]->setTimeAve(T_ave);
             csiHits[m_hitNum]->setTimeVar(0);
             csiHitRel.add(trackID, m_hitNum);
  
           } else {
             m_hitNum = CsiHitIndex[cellId][TimeIndex];
             double old_edep = csiHits[m_hitNum]->getEnergyDep();
             double old_TimeAve = csiHits[m_hitNum]->getTimeAve();
             double old_TimeVar = csiHits[m_hitNum]->getTimeVar();
  
             //double new_edep = E_cell + old_edep;
             //double new_TimeAve = (old_edep * old_TimeAve + E_cell * T_ave) / new_edep;
  
             double temp = E_cell + old_edep;
             double delta = T_ave - old_TimeAve;
             double R = delta * E_cell / temp;
             double new_TimeAve = old_TimeAve + R;
             double new_TimeVar = old_TimeVar + old_edep * delta * R;
             double new_edep = temp;
  
             //Algo from Wiki:
             /*
               for x, weight in dataWeightPairs:  # Alternatively "for x, weight in zip(data, weights):"
                   temp = weight + sumweight
                 delta = x - mean
                 R = delta * weight / temp
                 mean = mean + R
                 M2 = M2 + sumweight * delta * R  # Alternatively, "M2 = M2 + weight * delta * (x−mean)"
                 sumweight = temp
  
               variance_n = M2/sumweight
             */
  
             csiHits[m_hitNum]->setEnergyDep(new_edep);
             csiHits[m_hitNum]->setTimeAve(new_TimeAve);
             csiHits[m_hitNum]->setTimeVar(new_TimeVar);
             csiHitRel.add(trackID, m_hitNum);
  
           }
         }
       } else { B2ERROR("m_currentEvnetNumber  ERROR: m_oldEvnetNumber==m_oldEvnetNumber"); }
  
       return (simhitNumber);
     }//saveSimHit

◆ setActive()

static void setActive ( bool activeStatus )

inlinestaticinherited

Enable/Disable all Sensitive Detectors.

By default, all sensitive detectors won't create hits to make it possible to use the Geant4 Navigator for non-simulation purposes. Only during simulation the sensitive detectors will be enabled to record hits

Parameters

activeStatus bool to indicate wether hits should be recorded

Definition at line 50 of file SensitiveDetectorBase.h.

◆ step()

bool step	(	G4Step *	aStep,
		G4TouchableHistory *
	)

overrideprotectedvirtual

Step processing method.

Parameters

aStep the G4Step with the current step information

Returns: true if a Hit has been created, false if the hit was ignored

Implements SensitiveDetectorBase.

Definition at line 85 of file SensitiveDetector.cc.

Member Data Documentation

◆ m_cellID

int m_cellID

private

The current theta ID in an event.

Used to fill the DataStore ECL array

Definition at line 93 of file SensitiveDetector.h.

◆ m_EvnetNumber

int m_EvnetNumber

private

The current number of created hits in an event.

Used to fill the DataStore ECLHit.

Definition at line 69 of file SensitiveDetector.h.

◆ m_oldEvnetNumber

int m_oldEvnetNumber

private

The current number of created hits in an event.

Used to fill the DataStore ECL EB array.

Definition at line 70 of file SensitiveDetector.h.

◆ m_thetaID

int m_thetaID

private

The current phi ID in an event.

Used to fill the DataStore ECL array

Definition at line 92 of file SensitiveDetector.h.

◆ m_trackID

int m_trackID

private

The current number of created hits in an event.

Used to fill the DataStore

Definition at line 71 of file SensitiveDetector.h.

The documentation for this class was generated from the following files:

beast/csi/simulation/include/SensitiveDetector.h
beast/csi/simulation/src/SensitiveDetector.cc

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Function Documentation

◆ ProcessHits()

◆ registerMCParticleRelation() [1/2]

◆ registerMCParticleRelation() [2/2]

◆ saveSimHit()

◆ setActive()

◆ step()

Member Data Documentation

◆ m_cellID

◆ m_EvnetNumber

◆ m_oldEvnetNumber

◆ m_thetaID

◆ m_trackID