QweakSimTungstenPlug_DetectorSD.cc

Go to the documentation of this file.

// QweakSimTungstenPlug_DetectorSD.cc
// Martin McHugh
// 2015-05-04

/////// --------------------------------------------------------------------
#include "QweakSimTungstenPlug_DetectorSD.hh"

//--- user classes
#include "QweakSimTungstenPlug_DetectorHit.hh"
#include "QweakSimTrackInformation.hh"
#include "QweakSimSensitiveDetectorDefinition.hh"

/////// --------------------------------------------------------------------
QweakSimTungstenPlug_DetectorSD::QweakSimTungstenPlug_DetectorSD(G4String name)
:G4VSensitiveDetector(name)
{
    G4cout << name << G4endl;
    collectionName.insert("TungstenPlugCollection"); 
    TungstenPlug_CollectionID = -1;
}



/////// --------------------------------------------------------------------
QweakSimTungstenPlug_DetectorSD::~QweakSimTungstenPlug_DetectorSD()
{
    //delete TungstenPlug_DetectorHitsCollection; 
}



/////// --------------------------------------------------------------------

void QweakSimTungstenPlug_DetectorSD::Initialize(G4HCofThisEvent* HCE)
{ 
    TungstenPlug_DetectorHitsCollection = new QweakSimTungstenPlug_DetectorHitsCollection(SensitiveDetectorName,collectionName[0]);
    //if(TungstenPlug_CollectionID<0)
    TungstenPlug_CollectionID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]); 
    HCE->AddHitsCollection(TungstenPlug_CollectionID , TungstenPlug_DetectorHitsCollection);

}


/////// --------------------------------------------------------------------

G4bool QweakSimTungstenPlug_DetectorSD::ProcessHits(G4Step* aStep, G4TouchableHistory* )
{ 
    //--- Not Entering Geometry 
    if (aStep->GetPreStepPoint()->GetStepStatus() != fGeomBoundary ) return false;
  
    G4TouchableHandle theTouchable = aStep->GetPreStepPoint()->GetTouchableHandle();
  
    G4int         motherCopyNo = theTouchable->GetVolume(1)->GetCopyNo();   // one Mother Volume
  
    //G4double    parentID     = aStep->GetTrack()->GetParentID();
    G4double      trackID      = aStep->GetTrack()->GetTrackID();
  
    G4ParticleDefinition*  fpParticleDefinition = aStep->GetTrack()->GetDefinition();
    G4String  particleName = fpParticleDefinition->GetParticleName();
    //G4double  pdgCharge    = fpParticleDefinition->GetPDGCharge();
    G4int pdgEncoding  = fpParticleDefinition->GetPDGEncoding();
  
    //GlobalTimeOfHit
    G4double      globalTime   = aStep->GetPreStepPoint()->GetGlobalTime();
  
    G4ThreeVector worldPos             = aStep->GetPreStepPoint()->GetPosition();
    G4ThreeVector localPos             = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldPos);

    // Local Origin (where did the hit come from)
    G4ThreeVector originVertexPosition = aStep->GetTrack()->GetVertexPosition();

    G4ThreeVector worldMomentum                 = aStep->GetPreStepPoint()->GetMomentum();
    G4ThreeVector localMomentum                 = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldMomentum);
    G4ThreeVector currentMomentumDirection      = aStep->GetTrack()->GetMomentumDirection();
    G4ThreeVector originVertexMomentumDirection = aStep->GetTrack()->GetVertexMomentumDirection();

    G4double      originVertexKineticEnergy = aStep->GetTrack()->GetVertexKineticEnergy();
    G4double      originVertexTotalEnergy   = aStep->GetTrack()->GetVertexKineticEnergy();//GetVertexTotalEnergy();
  
    G4double      currentKineticEnergy      = aStep->GetTrack()->GetKineticEnergy();
    G4double      currentTotalEnergy        = aStep->GetTrack()->GetTotalEnergy();
  
    //--- TungstenPlug deposited energy
    G4double      depositedEnergy      = aStep->GetTotalEnergyDeposit();
  
    //----------------- Store Hit information
    QweakSimTungstenPlug_DetectorHit* aHit = new QweakSimTungstenPlug_DetectorHit(motherCopyNo);
  
    aHit -> StoreTrackID(trackID);
  
    aHit -> StoreParticleName(particleName);
    aHit -> StoreParticleType(pdgEncoding);
  
    aHit -> StoreGlobalTime(globalTime);
    
    aHit -> StoreWorldPosition(worldPos);
    aHit -> StoreLocalPosition(localPos);
    aHit -> StoreOriginVertexPosition(originVertexPosition);
  
    aHit -> StoreWorldMomentum(worldMomentum);
    aHit -> StoreLocalMomentum(localMomentum);
    aHit -> StoreMomentumDirection(currentMomentumDirection);
    aHit -> StoreOriginVertexMomentumDirection(originVertexMomentumDirection);
  
    aHit -> StoreOriginVertexKineticEnergy(originVertexKineticEnergy);
    aHit -> StoreOriginVertexTotalEnergy(originVertexTotalEnergy);  // kinetic energy
  
    aHit -> StoreKineticEnergy(currentKineticEnergy);
    aHit -> StoreTotalEnergy(currentTotalEnergy);
  
    //--- TungstenPlug deposited energy
    aHit -> StoreDepositedEnergy(depositedEnergy);
  
    TungstenPlug_DetectorHitsCollection->insert(aHit);  
    
    return true;
}



/////// --------------------------------------------------------------------

void QweakSimTungstenPlug_DetectorSD::EndOfEvent(G4HCofThisEvent* )
{
    //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();
  
}



/////// --------------------------------------------------------------------