CCaloSD Class Reference

#include <CCaloSD.hh>

Inheritance diagram for CCaloSD:

Collaboration diagram for CCaloSD:

Public Member Functions
	CCaloSD (G4String aSDname, CCalVOrganization *numberingScheme)
virtual	~CCaloSD ()
void	Initialize (G4HCofThisEvent *HCE)
G4bool	ProcessHits (G4Step *aStep, G4TouchableHistory *ROhist)
void	EndOfEvent (G4HCofThisEvent *HCE)
void	clear ()
void	DrawAll ()
void	PrintAll ()
void	SetPrimaryID (int i)
int	GetPrimaryID ()
void	SetOrganization (CCalVOrganization *org)
Public Member Functions inherited from G4VSensitiveDetector
	G4VSensitiveDetector (G4String name)
	G4VSensitiveDetector (const G4VSensitiveDetector &right)
virtual	~G4VSensitiveDetector ()
G4VSensitiveDetector &	operator= (const G4VSensitiveDetector &right)
G4int	operator== (const G4VSensitiveDetector &right) const
G4int	operator!= (const G4VSensitiveDetector &right) const
G4bool	Hit (G4Step *aStep)
void	SetROgeometry (G4VReadOutGeometry *value)
void	SetFilter (G4VSDFilter *value)
G4int	GetNumberOfCollections () const
G4String	GetCollectionName (G4int id) const
void	SetVerboseLevel (G4int vl)
void	Activate (G4bool activeFlag)
G4bool	isActive () const
G4String	GetName () const
G4String	GetPathName () const
G4String	GetFullPathName () const
G4VReadOutGeometry *	GetROgeometry () const
G4VSDFilter *	GetFilter () const
virtual G4VSensitiveDetector *	Clone () const

Private Member Functions
G4ThreeVector	SetToLocal (G4ThreeVector globalPoint)
void	getStepInfo (G4Step *aStep)
G4bool	hitExists ()
void	createNewHit ()
void	updateHit ()
void	StoreHit (CCalG4Hit *ahit)
void	ResetForNewPrimary ()
void	summarize ()
G4double	curve_LY (G4StepPoint *stepPoint)

Private Attributes
G4ThreeVector	EntrancePoint
float	IncidentEnergy
G4int	PrimID
G4int	HCID
G4String	SDname
CCalG4HitCollection *	theHC
G4int	TSID
CCalG4Hit *	CurrentHit
G4Track *	theTrack
G4VPhysicalVolume *	CurrentPV
G4VPhysicalVolume *	PreviousPV
unsigned int	UnitID
unsigned int	PreviousUnitID
G4int	PrimaryID
G4int	TSliceID
G4double	TSlice
G4StepPoint *	PreStepPoint
G4StepPoint *	PostStepPoint
float	EdepositEM
float	EdepositEHAD
G4ThreeVector	HitPoint
CCalVOrganization *	theDescription

Additional Inherited Members
Protected Member Functions inherited from G4VSensitiveDetector
virtual G4int	GetCollectionID (G4int i)
Protected Attributes inherited from G4VSensitiveDetector
G4CollectionNameVector	collectionName
G4String	SensitiveDetectorName
G4String	thePathName
G4String	fullPathName
G4int	verboseLevel
G4bool	active
G4VReadOutGeometry *	ROgeometry
G4VSDFilter *	filter

Detailed Description

Definition at line 59 of file CCaloSD.hh.

Constructor & Destructor Documentation

CCaloSD()

CCaloSD::CCaloSD	(	G4String	aSDname,
		CCalVOrganization *	numberingScheme
	)

Definition at line 42 of file CCaloSD.cc.

                                                                 :
  G4VSensitiveDetector(name), HCID(-1), SDname(name), theHC(0),
  CurrentHit(0), theTrack(0), CurrentPV(0), PreviousPV(0), UnitID(0), 
  PreviousUnitID(0), PreStepPoint(0), PostStepPoint(0), 
  theDescription(numberingScheme) {
  
  collectionName.insert(name);
  
  G4cout << "*******************************************************" << G4endl;
  G4cout << "*                                                     *" << G4endl;
  G4cout << "* Constructing a CCaloSD  with name " << name            << G4endl;  
  G4cout << "*                                                     *" << G4endl;
  G4cout << "*******************************************************" << G4endl;

  CCalSDList::getInstance()->addCalo(name);
}

Here is the call graph for this function:

~CCaloSD()

CCaloSD::~CCaloSD ( )

virtual

Definition at line 60 of file CCaloSD.cc.

                  {
  G4cout << "CCaloSD: " << SDname << " deleted" << G4endl;
  if (theDescription) 
    delete theDescription;
}

Member Function Documentation

clear()

void CCaloSD::clear ( void  )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 290 of file CCaloSD.cc.

                    {
} 

createNewHit()

void CCaloSD::createNewHit ( )

private

Definition at line 202 of file CCaloSD.cc.

                           {

#ifdef debug
  G4int currentCopyNo = -999;
  G4int motherCopyNo  = -999;
  G4TouchableHistory* theTouchable
    = (G4TouchableHistory*)( theTrack->GetTouchable() );
  if ( theTouchable ) {
    currentCopyNo = theTouchable->GetReplicaNumber( 0 );
    if ( theTouchable->GetHistoryDepth() > 0 ) {
      motherCopyNo = theTouchable->GetReplicaNumber( 1 );
    }
  }
  G4cout << "CCaloSD createNewHit for"
         << " PV "     << CurrentPV->GetName()
         << " PVid = " << currentCopyNo
         << " MVid = " << motherCopyNo
         << " Unit "   << UnitID <<G4endl;
  G4cout << " primary "    << PrimaryID
         << " time slice " << TSliceID 
         << " For Track  " << theTrack->GetTrackID()
         << " which is a " <<  theTrack->GetDefinition()->GetParticleName();
       
  if (theTrack->GetTrackID()==1) {
    G4cout << " of energy "     << theTrack->GetTotalEnergy();
  } else {
    G4cout << " daughter of part. " << theTrack->GetParentID();
  }

  G4cout  << " and created by " ;
  if (theTrack->GetCreatorProcess()!=NULL)
    G4cout << theTrack->GetCreatorProcess()->GetProcessName() ;
  else 
    G4cout << "NO process";
  G4cout << G4endl;
#endif          
    

  CurrentHit = new CCalG4Hit;
  CurrentHit->setTrackID(PrimaryID);
  CurrentHit->setTimeSlice(TSlice);
  CurrentHit->setUnitID(UnitID);
  CurrentHit->setEntry(EntrancePoint);
  CurrentHit->setIncidentEnergy(IncidentEnergy);
  updateHit();
  
  StoreHit(CurrentHit);

}    

Here is the call graph for this function:

Here is the caller graph for this function:

curve_LY()

G4double CCaloSD::curve_LY ( G4StepPoint *  stepPoint )

private

Definition at line 312 of file CCaloSD.cc.

                                                 {

  G4double weight = 1.;
  G4ThreeVector  localPoint = SetToLocal(stepPoint->GetPosition());
  G4double crlength = 230.;
  G4double dapd = 0.5 * crlength - localPoint.z();
  if (dapd >= -0.1 || dapd <= crlength+0.1) {
    if (dapd <= 100.)
      weight = 1.05 - dapd * 0.0005;
  } else {
    G4cout << "CCaloSD, light coll curve : wrong distance to APD " << dapd 
       << " crlength = " << crlength
       << " z of localPoint = " << localPoint.z() 
       << " take weight = " << weight << G4endl;
  }
#ifdef ddebug
  G4cout << "CCaloSD, light coll curve : " << dapd 
     << " crlength = " << crlength
     << " z of localPoint = " << localPoint.z() 
     << " take weight = " << weight << G4endl;
#endif
  return weight;
}

Here is the call graph for this function:

Here is the caller graph for this function:

DrawAll()

void CCaloSD::DrawAll ( void  )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 294 of file CCaloSD.cc.

                      {
} 

EndOfEvent()

void CCaloSD::EndOfEvent ( G4HCofThisEvent *  HCE )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 281 of file CCaloSD.cc.

                                         {
  summarize();
}

Here is the call graph for this function:

GetPrimaryID()

int CCaloSD::GetPrimaryID ( )

inline

Definition at line 79 of file CCaloSD.hh.

{return PrimaryID;}

Here is the call graph for this function:

getStepInfo()

void CCaloSD::getStepInfo ( G4Step *  aStep )

private

Definition at line 98 of file CCaloSD.cc.

                                       {
  
  PreStepPoint = aStep->GetPreStepPoint(); 
  PostStepPoint= aStep->GetPostStepPoint(); 
  HitPoint     = PreStepPoint->GetPosition();   

  theTrack = aStep->GetTrack();   
  CurrentPV= PreStepPoint->GetPhysicalVolume();

  G4String     particleType =  theTrack->GetDefinition()->GetParticleName();

  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
  G4double weight;
  if (touch->GetVolume(0)->GetName() == "CrystalMatrixCrystal")
    weight = curve_LY(PreStepPoint);
  else
    weight = 1.;

  if (particleType == "e-" ||
      particleType == "e+" ||
      particleType == "gamma" ){
    EdepositEM   = weight*(aStep->GetTotalEnergyDeposit());
    EdepositEHAD = 0.;
  } else {
    EdepositEM   = 0.;
    EdepositEHAD = weight*(aStep->GetTotalEnergyDeposit());
  }

  TSlice = PostStepPoint->GetGlobalTime() / nanosecond;

  if ( TSlice > 1.0E9 ) TSliceID = 999999999;
  else                  TSliceID = (int) TSlice;

  if (theDescription!=0) 
    UnitID = theDescription->GetUnitID(aStep);
  else
    UnitID = 0;
   
}

Here is the call graph for this function:

Here is the caller graph for this function:

hitExists()

G4bool CCaloSD::hitExists ( )

private

Definition at line 139 of file CCaloSD.cc.

                          {
   
  if (PrimaryID<1) {
    G4cerr << "***** CCaloSD error: PrimaryID = " << PrimaryID
       << " Maybe detector name changed"
       << G4endl;
  }
   
      
  if ( CurrentPV==PreviousPV && PrimaryID == PrimID && TSliceID == TSID &&
       UnitID==PreviousUnitID) {
    updateHit();
    return true;
  }
   

  if (PrimaryID != PrimID)
    ResetForNewPrimary();
   

  //look in HC whether a hit with the same primID,UnitID,TSliceID exists already:
   
  G4bool found = false;
  for (G4int j=0; j<theHC->entries()&&!found; j++) {

    CCalG4Hit* aPreviousHit = (*theHC)[j];
    if (aPreviousHit->getTrackID()  == PrimaryID &&
    aPreviousHit->getTimeSliceID() == TSliceID  &&
    aPreviousHit->getUnitID()== UnitID       ) {
      CurrentHit = aPreviousHit;
      found = true;
    }
  }          

  if (found) {
    updateHit();
    return true;
  } else {
    return false;
  }    
}

Here is the call graph for this function:

Here is the caller graph for this function:

Initialize()

void CCaloSD::Initialize ( G4HCofThisEvent *  HCE )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 67 of file CCaloSD.cc.

                                            {

#ifdef debug
  G4cout << "CCaloSD : Initialize called for " << SDname << G4endl;
#endif
  //This initialization is performed at the beginning of an event
  //------------------------------------------------------------

  theHC = new CCalG4HitCollection(SDname, collectionName[0]);
  if (HCID<0) 
    HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
  HCE->AddHitsCollection( HCID, theHC );

  TSID = -2;
  PrimID = -2;
}

Here is the call graph for this function:

PrintAll()

void CCaloSD::PrintAll ( void  )

virtual

Reimplemented from G4VSensitiveDetector.

Definition at line 298 of file CCaloSD.cc.

                       {
  G4cout << "CCaloSD: Collection " << theHC->GetName() << G4endl;
  theHC->PrintAllHits();
} 

Here is the call graph for this function:

ProcessHits()

G4bool CCaloSD::ProcessHits ( G4Step *  aStep, G4TouchableHistory *  ROhist  )

virtual

Implements G4VSensitiveDetector.

Definition at line 86 of file CCaloSD.cc.

                                                            {

  if (aStep == NULL) return true;
   
  getStepInfo(aStep);
  if (hitExists() == false && EdepositEM+EdepositEHAD>0.) 
    createNewHit();

  return true;
} 

Here is the call graph for this function:

ResetForNewPrimary()

void CCaloSD::ResetForNewPrimary ( )

private

Definition at line 182 of file CCaloSD.cc.

                                 {
  
  EntrancePoint = SetToLocal(HitPoint);
  IncidentEnergy = PreStepPoint->GetKineticEnergy();

}

Here is the call graph for this function:

Here is the caller graph for this function:

SetOrganization()

void CCaloSD::SetOrganization

(

CCalVOrganization *

org

)

Definition at line 304 of file CCaloSD.cc.

                                                   {

  if (theDescription!=0) 
    delete theDescription;
  theDescription = org;
}

Here is the caller graph for this function:

SetPrimaryID()

void CCaloSD::SetPrimaryID ( int  i )

inline

Definition at line 74 of file CCaloSD.hh.

                           {PrimaryID = i;
#ifdef debug
  G4cout << "CCaloSD SetPrimaryID primID =" << i << G4endl;
#endif
  }

Here is the caller graph for this function:

SetToLocal()

G4ThreeVector CCaloSD::SetToLocal ( G4ThreeVector  globalPoint )

private

Definition at line 270 of file CCaloSD.cc.

                                                     {

  G4ThreeVector localPoint;
  const G4VTouchable*   touch= PreStepPoint->GetTouchable();
  localPoint=touch->GetHistory()->GetTopTransform().TransformPoint(global);
  
  return localPoint;  

}

Here is the call graph for this function:

Here is the caller graph for this function:

StoreHit()

void CCaloSD::StoreHit ( CCalG4Hit *  ahit )

private

Definition at line 190 of file CCaloSD.cc.

                                    {

  if (PrimID<0) return;
  if (hit == 0 ) {
    G4cout << "CCaloSD: hit to be stored is NULL !!" <<G4endl;
    return;
  }

  theHC->insert( hit );
}

Here is the call graph for this function:

Here is the caller graph for this function:

summarize()

void CCaloSD::summarize ( )

private

Definition at line 286 of file CCaloSD.cc.

                        {
}

Here is the caller graph for this function:

updateHit()

void CCaloSD::updateHit ( )

private

Definition at line 253 of file CCaloSD.cc.

                        {
  if (EdepositEM+EdepositEHAD != 0) {
    CurrentHit->addEnergyDeposit(EdepositEM,EdepositEHAD);
#ifdef debug
    G4cout << "Energy deposit in Unit " << UnitID << " em " << EdepositEM/MeV
     << " hadronic " << EdepositEHAD/MeV << " MeV" << G4endl;
#endif
  }

  // buffer for next steps:
  TSID = TSliceID;
  PrimID = PrimaryID;
  PreviousPV = CurrentPV;
  PreviousUnitID = UnitID;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

CurrentHit

CCalG4Hit* CCaloSD::CurrentHit

private

Definition at line 100 of file CCaloSD.hh.

CurrentPV

G4VPhysicalVolume* CCaloSD::CurrentPV

private

Definition at line 102 of file CCaloSD.hh.

EdepositEHAD

float CCaloSD::EdepositEHAD

private

Definition at line 110 of file CCaloSD.hh.

EdepositEM

float CCaloSD::EdepositEM

private

Definition at line 110 of file CCaloSD.hh.

EntrancePoint

G4ThreeVector CCaloSD::EntrancePoint

private

Definition at line 89 of file CCaloSD.hh.

HCID

G4int CCaloSD::HCID

private

Definition at line 95 of file CCaloSD.hh.

HitPoint

G4ThreeVector CCaloSD::HitPoint

private

Definition at line 111 of file CCaloSD.hh.

IncidentEnergy

float CCaloSD::IncidentEnergy

private

Definition at line 90 of file CCaloSD.hh.

PostStepPoint

G4StepPoint* CCaloSD::PostStepPoint

private

Definition at line 109 of file CCaloSD.hh.

PreStepPoint

G4StepPoint* CCaloSD::PreStepPoint

private

Definition at line 108 of file CCaloSD.hh.

PreviousPV

G4VPhysicalVolume* CCaloSD::PreviousPV

private

Definition at line 103 of file CCaloSD.hh.

PreviousUnitID

unsigned int CCaloSD::PreviousUnitID

private

Definition at line 104 of file CCaloSD.hh.

PrimaryID

G4int CCaloSD::PrimaryID

private

Definition at line 105 of file CCaloSD.hh.

PrimID

G4int CCaloSD::PrimID

private

Definition at line 91 of file CCaloSD.hh.

SDname

G4String CCaloSD::SDname

private

Definition at line 96 of file CCaloSD.hh.

theDescription

CCalVOrganization* CCaloSD::theDescription

private

Definition at line 126 of file CCaloSD.hh.

theHC

CCalG4HitCollection* CCaloSD::theHC

private

Definition at line 97 of file CCaloSD.hh.

theTrack

G4Track* CCaloSD::theTrack

private

Definition at line 101 of file CCaloSD.hh.

TSID

G4int CCaloSD::TSID

private

Definition at line 99 of file CCaloSD.hh.

TSlice

G4double CCaloSD::TSlice

private

Definition at line 106 of file CCaloSD.hh.

TSliceID

G4int CCaloSD::TSliceID

private

Definition at line 105 of file CCaloSD.hh.

UnitID

unsigned int CCaloSD::UnitID

private

Definition at line 104 of file CCaloSD.hh.

---

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

• CCaloSD.hh
• CCaloSD.cc