QweakSimTriggerScintillator_DetectorSD Class Reference

Storing hit information of the Trigger Scintillator detector. More...

#include <QweakSimTriggerScintillator_DetectorSD.hh>

Inherits G4VSensitiveDetector.

Public Member Functions
	QweakSimTriggerScintillator_DetectorSD (G4String name)
	~QweakSimTriggerScintillator_DetectorSD ()
void	Initialize (G4HCofThisEvent *HCE)
G4bool	ProcessHits (G4Step *aStep, G4TouchableHistory *ROhist)
void	EndOfEvent (G4HCofThisEvent *HCE)

Private Attributes
QweakSimTriggerScintillator_DetectorHitsCollection *	TriggerScintillator_DetectorHitsCollection
G4int	TriggerScintillator_CollectionID

Detailed Description

Storing hit information of the Trigger Scintillator detector.

Placeholder for a long explaination

Definition at line 57 of file QweakSimTriggerScintillator_DetectorSD.hh.

Constructor & Destructor Documentation

QweakSimTriggerScintillator_DetectorSD::QweakSimTriggerScintillator_DetectorSD

(

G4String

name

)

Definition at line 15 of file QweakSimTriggerScintillator_DetectorSD.cc.

References TriggerScintillator_CollectionID.

:G4VSensitiveDetector(name)
{
  //G4cout << G4endl << "###### Calling QweakSimTriggerScintillator_DetectorSD::QweakSimTriggerScintillator_DetectorSD() " << G4endl << G4endl; 

  collectionName.insert("TriggerScintillatorCollection"); 
  TriggerScintillator_CollectionID = -1;

  //G4cout << G4endl << "###### Leaving QweakSimTriggerScintillator_DetectorSD::QweakSimTriggerScintillator_DetectorSD() " << G4endl << G4endl; 
}

QweakSimTriggerScintillator_DetectorSD::~QweakSimTriggerScintillator_DetectorSD

(

)

Definition at line 27 of file QweakSimTriggerScintillator_DetectorSD.cc.

{
  //G4cout << G4endl << "###### Calling QweakSimTriggerScintillator_DetectorSD::~QweakSimTriggerScintillator_DetectorSD() " << G4endl << G4endl; 

  //delete TriggerScintillator_HitsCollection;

}

Member Function Documentation

void QweakSimTriggerScintillator_DetectorSD::EndOfEvent

(

G4HCofThisEvent *

HCE

)

Definition at line 142 of file QweakSimTriggerScintillator_DetectorSD.cc.

{
  //G4cout << G4endl << "###### Calling QweakSimTriggerScintillator_DetectorSD::EndOfEvent() " << G4endl << G4endl; 
  
//   G4int NbDCHits = DC_hitsCollection->entries();
  
//   G4cout << "\n-------->Hits Collection: in this event they are " << NbDCHits
//   << " hits in the Drift Cells : " << G4endl;
//   for (G4int i=0;i<NbDCHits;i++) (*DC_hitsCollection)[i]->Print();
  


   //G4cout << G4endl << "###### Leaving QweakSimTriggerScintillator_DetectorSD::EndOfEvent() " << G4endl << G4endl; 

}

void QweakSimTriggerScintillator_DetectorSD::Initialize

(

G4HCofThisEvent *

HCE

)

Definition at line 36 of file QweakSimTriggerScintillator_DetectorSD.cc.

References TriggerScintillator_CollectionID, and TriggerScintillator_DetectorHitsCollection.

{
   //G4cout << G4endl << "###### Calling QweakSimTriggerScintillator_DetectorSD::Initialize() " << G4endl << G4endl; 

   TriggerScintillator_DetectorHitsCollection = new QweakSimTriggerScintillator_DetectorHitsCollection(SensitiveDetectorName,collectionName[0]);

   //if(TriggerScintillator_CollectionID<0)

   TriggerScintillator_CollectionID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]); 

   HCE->AddHitsCollection(TriggerScintillator_CollectionID , TriggerScintillator_DetectorHitsCollection);
   
   //G4cout << G4endl << "###### Leaving QweakSimTriggerScintillator_DetectorSD::Initialize() " << G4endl << G4endl; 
}

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

Todo:

beware: total.neq.kinetic (testing only)

Definition at line 53 of file QweakSimTriggerScintillator_DetectorSD.cc.

References QweakSimTriggerScintillator_DetectorHit::GetLogVolume(), QweakSimTriggerScintillator_DetectorHit::StoreCellPosition(), QweakSimTriggerScintillator_DetectorHit::StoreCellRotation(), QweakSimTriggerScintillator_DetectorHit::StoreDepositedEnergy(), QweakSimTriggerScintillator_DetectorHit::StoreGlobalTime(), QweakSimTriggerScintillator_DetectorHit::StoreKineticEnergy(), QweakSimTriggerScintillator_DetectorHit::StoreLocalMomentum(), QweakSimTriggerScintillator_DetectorHit::StoreLocalPosition(), QweakSimTriggerScintillator_DetectorHit::StoreLogVolume(), QweakSimTriggerScintillator_DetectorHit::StoreMomentumDirection(), QweakSimTriggerScintillator_DetectorHit::StoreOriginVertexKineticEnergy(), QweakSimTriggerScintillator_DetectorHit::StoreOriginVertexMomentumDirection(), QweakSimTriggerScintillator_DetectorHit::StoreOriginVertexPosition(), QweakSimTriggerScintillator_DetectorHit::StoreOriginVertexTotalEnergy(), QweakSimTriggerScintillator_DetectorHit::StoreParticleName(), QweakSimTriggerScintillator_DetectorHit::StoreParticleType(), QweakSimTriggerScintillator_DetectorHit::StoreTotalEnergy(), QweakSimTriggerScintillator_DetectorHit::StoreTrackID(), QweakSimTriggerScintillator_DetectorHit::StoreWorldMomentum(), QweakSimTriggerScintillator_DetectorHit::StoreWorldPosition(), and TriggerScintillator_DetectorHitsCollection.

{
//   G4cout << G4endl << "###### Calling QweakSimTriggerScintillator_DetectorSD::ProcessHits() " << G4endl << G4endl; 

  G4Track *track = aStep->GetTrack();    
  G4String  particlename = track->GetDefinition()->GetParticleName();
  G4int     particletype = track->GetDefinition()->GetPDGEncoding();

  G4ThreeVector worldPos;
  G4ThreeVector localPos;
  G4ThreeVector worldMomentum;
  
  G4TouchableHandle theTouchable = aStep->GetPreStepPoint()->GetTouchableHandle();
 
  if (aStep->GetPreStepPoint()->GetStepStatus() == fGeomBoundary ) {
    // Entering Geometry
    worldPos      = aStep->GetPreStepPoint()->GetPosition();
    localPos      = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldPos);
    worldMomentum = aStep->GetPreStepPoint()->GetMomentum();
  } else { 
    return false;
  }

  G4ThreeVector localMomentum = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldMomentum);
 
  // Vertex: where this track was generated
  G4ThreeVector originVertexPosition          = aStep->GetTrack()->GetVertexPosition();
  G4ThreeVector originVertexMomentumDirection = aStep->GetTrack()->GetVertexMomentumDirection();
  G4double      originVertexKineticEnergy     = aStep->GetTrack()->GetVertexKineticEnergy();

  G4double      currentKineticEnergy     = aStep->GetTrack()->GetKineticEnergy();
  G4double      currentTotalEnergy       = aStep->GetTrack()->GetTotalEnergy();
  G4ThreeVector currentMomentumDirection = aStep->GetTrack()->GetMomentumDirection();

  G4int trackID = aStep->GetTrack()->GetTrackID();
  
  G4int    MotherCopyNo    = theTouchable->GetVolume(1)->GetCopyNo();   // one Mother Volume

  QweakSimTriggerScintillator_DetectorHit* aHit = new QweakSimTriggerScintillator_DetectorHit(MotherCopyNo);
  
  // Energy deposited in TS 
  G4double depositedEnergy = aStep->GetTotalEnergyDeposit();
  aHit->StoreDepositedEnergy(depositedEnergy);

  aHit->StoreTrackID(trackID);

  aHit->StoreGlobalTime(aStep->GetPreStepPoint()->GetGlobalTime());

  aHit->StoreParticleName(particlename);
  aHit->StoreParticleType(particletype);
    
  aHit->StoreWorldPosition(worldPos);
  aHit->StoreLocalPosition(localPos);

  aHit->StoreWorldMomentum(worldMomentum);
  aHit->StoreLocalMomentum(localMomentum);

  aHit->StoreOriginVertexPosition(originVertexPosition);
  aHit->StoreOriginVertexKineticEnergy(originVertexKineticEnergy);
  aHit->StoreOriginVertexTotalEnergy(originVertexKineticEnergy);   /// \todo beware: total.neq.kinetic (testing only)
  aHit->StoreOriginVertexMomentumDirection(originVertexMomentumDirection);

  aHit->StoreMomentumDirection(currentMomentumDirection);
  aHit->StoreKineticEnergy(currentKineticEnergy);
  aHit->StoreTotalEnergy(currentTotalEnergy);

  // check if it is first touch
  if(!(aHit->GetLogVolume()))
    {
      // store translation and rotation matrix of the drift cell 
      // for the sake of drawing the hit 
      aHit->StoreLogVolume(theTouchable->GetVolume()->GetLogicalVolume());
      G4AffineTransform aTrans = theTouchable->GetHistory()->GetTopTransform();
      aTrans.Invert();
      aHit->StoreCellRotation(aTrans.NetRotation());
      aHit->StoreCellPosition(aTrans.NetTranslation());
    }
 
//   G4cout << "%%%%%%%%%%%%%%%%%%%%%%%  before inserting the hit" << G4endl;

  TriggerScintillator_DetectorHitsCollection->insert(aHit); 
  
//   G4cout << G4endl << "###### Leaving QweakSimTriggerScintillatorSD::ProcessHits() " << G4endl << G4endl; 
  
  return true;
}

Here is the call graph for this function:

Field Documentation

G4int QweakSimTriggerScintillator_DetectorSD::TriggerScintillator_CollectionID

private

Definition at line 72 of file QweakSimTriggerScintillator_DetectorSD.hh.

Referenced by Initialize(), and QweakSimTriggerScintillator_DetectorSD().

QweakSimTriggerScintillator_DetectorHitsCollection* QweakSimTriggerScintillator_DetectorSD::TriggerScintillator_DetectorHitsCollection

private

Definition at line 70 of file QweakSimTriggerScintillator_DetectorSD.hh.

Referenced by Initialize(), and ProcessHits().

---

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

• QweakSimTriggerScintillator_DetectorSD.hh
• QweakSimTriggerScintillator_DetectorSD.cc