QweakSimLumi_DetectorSD Class Reference

#include <QweakSimLumi_DetectorSD.hh>

Inherits G4VSensitiveDetector.

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

Private Attributes
QweakSimLumi_DetectorHitsCollection *	Lumi_DetectorHitsCollection
G4int	Lumi_CollectionID

Detailed Description

Definition at line 7 of file QweakSimLumi_DetectorSD.hh.

Constructor & Destructor Documentation

QweakSimLumi_DetectorSD::QweakSimLumi_DetectorSD

(

G4String

name

)

Definition at line 6 of file QweakSimLumi_DetectorSD.cc.

References Lumi_CollectionID.

                                                             :G4VSensitiveDetector(name) {
    collectionName.insert("LumiCollection");
    Lumi_CollectionID = -1;
}

QweakSimLumi_DetectorSD::~QweakSimLumi_DetectorSD

(

)

Definition at line 11 of file QweakSimLumi_DetectorSD.cc.

                                                  {
    //delete Lumi_HitCollection;
}

Member Function Documentation

void QweakSimLumi_DetectorSD::EndOfEvent

(

G4HCofThisEvent *

HCE

)

Definition at line 115 of file QweakSimLumi_DetectorSD.cc.

                                                         {
    // placeholder 
}

void QweakSimLumi_DetectorSD::Initialize

(

G4HCofThisEvent *

HCE

)

Definition at line 15 of file QweakSimLumi_DetectorSD.cc.

References Lumi_CollectionID, and Lumi_DetectorHitsCollection.

                                                             {
    Lumi_DetectorHitsCollection = new QweakSimLumi_DetectorHitsCollection(SensitiveDetectorName, collectionName[0]);
    Lumi_CollectionID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
    HCE->AddHitsCollection(Lumi_CollectionID, Lumi_DetectorHitsCollection);
}

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

Definition at line 21 of file QweakSimLumi_DetectorSD.cc.

References Lumi_DetectorHitsCollection.

                                                                              {
    // Kill chargeless particles that don't deposit much energy
    G4double charge = aStep->GetTrack()->GetDefinition()->GetPDGCharge();
    if(charge==0. && aStep->GetTrack()->GetTotalEnergy()/MeV < 0.1) {
        return false;
    }

    // Kill optical photons (don't need to save these tracks
    if( aStep->GetTrack()->GetDefinition()->GetParticleName() == "opticalphoton") {
        return false;
    }
    // begin experimental
    if(aStep->GetPreStepPoint()->GetStepStatus() != fGeomBoundary)
        return false;
    // end experimental
    G4TouchableHandle theTouchable = aStep->GetPreStepPoint()->GetTouchableHandle();
    G4int motherCopyNo = theTouchable->GetVolume(1)->GetCopyNo();

    G4double trackID = aStep->GetTrack()->GetTrackID();

    G4ParticleDefinition* fpParticleDefinition = aStep->GetTrack()->GetDefinition();
    G4String particleName = fpParticleDefinition->GetParticleName();

    G4int pdgEncoding = fpParticleDefinition->GetPDGEncoding();

    G4double      globalTime   = aStep->GetPreStepPoint()->GetGlobalTime();

    G4ThreeVector worldPos             = aStep->GetPreStepPoint()->GetPosition();
    G4ThreeVector localPos             = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldPos);
    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();

    G4double      depositedEnergy      = aStep->GetTotalEnergyDeposit();

    QweakSimLumi_DetectorHit* aHit = new QweakSimLumi_DetectorHit(motherCopyNo);

    aHit -> StoreTrackID(trackID);

    aHit -> StoreParticleName(particleName);
    aHit -> StoreParticleType(pdgEncoding);

    aHit -> StoreGlobalTime(globalTime);

    //aHit -> StoreHasBeenHit(hasBeenHit);
    //aHit -> StoreEdgeEventFlag(edgeEventFlag);
    //aHit -> SotreNbOfHits(nbOfHits);

    aHit -> StoreWorldPosition(worldPos);
    aHit -> StoreLocalPosition(localPos);
    //aHit -> StoreLocalExitPosition(localExitPos);
    aHit -> StoreOriginVertexPosition(originVertexPosition);

    aHit -> StoreWorldMomentum(worldMomentum);
    aHit -> StoreLocalMomentum(localMomentum);
    aHit -> StoreMomentumDirection(currentMomentumDirection);
    aHit -> StoreOriginVertexMomentumDirection(originVertexMomentumDirection);

    //aHit -> StoreOriginVertexThetaAngle; //calculated in EventAction
    //aHit -> StoreOriginVertexPhiAngle;

    aHit -> StoreOriginVertexKineticEnergy(originVertexKineticEnergy);
    aHit -> StoreOriginVertexTotalEnergy(originVertexTotalEnergy);  // kinetic energy

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

    aHit -> StoreHitDepositedEnergy(depositedEnergy);

    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());
    }

    Lumi_DetectorHitsCollection->insert(aHit);  

    return true;
}

Field Documentation

G4int QweakSimLumi_DetectorSD::Lumi_CollectionID

private

Definition at line 18 of file QweakSimLumi_DetectorSD.hh.

Referenced by Initialize(), and QweakSimLumi_DetectorSD().

QweakSimLumi_DetectorHitsCollection* QweakSimLumi_DetectorSD::Lumi_DetectorHitsCollection

private

Definition at line 17 of file QweakSimLumi_DetectorSD.hh.

Referenced by Initialize(), and ProcessHits().

---

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

• QweakSimLumi_DetectorSD.hh
• QweakSimLumi_DetectorSD.cc