Message: Re: Counting Number of Rkt Li6(n,t)alpha Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Counting Number of Rkt Li6(n,t)alpha 

Keywords: reaktion, secondaries
Forum: Hits, Digitization and Pileup
Re: None Counting Number of Rkt Li6(n,t)alpha (Sylvia Studeny)
Date: 01 Aug, 2006
From: Sylvia Studeny <sylvia.studeny@gsf.de>

Hi,

I finally solved the problem in another way (I found another listing to the topic: 35). I count the reaction in UserSteppingAction, where it is possible to access secondary particle information via the SteppingManager:

void RE02SteppingAction::UserSteppingAction(const G4Step * theStep) {

G4Track * theTrack = theStep->GetTrack();

  // check if it is alive
  if(theTrack->GetTrackStatus()!=fAlive) { return; }

  // check if it is primary
  if(theTrack->GetParentID()!=0) { return; }

 fSecondary = fpSteppingManager->GetfSecondary();
 anzahlSekTeilchen = fpSteppingManager-> GetfN2ndariesPostStepDoIt();
 if(anzahlSekTeilchen!=2){return;}
  particleName0 = (*fSecondary)[0]->GetDefinition()->GetParticleName();
  particleName1 = (*fSecondary)[1]->GetDefinition()->GetParticleName();

 if(particleName0=="triton" && particleName1== "alpha"){
   if(theStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()
                                                          ->GetName()=="TLD"){
      copyNumber = 0;
      copyNumber = theStep->GetPostStepPoint()->GetTouchable()->GetReplicaNumber();
      anzahlRkt[copyNumber]++;
   }
 }
}

How this is possible in a primitive scorer I do not know. Sylvia

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

 Add Message Add Message
to: "Re: Counting Number of Rkt Li6(n,t)alpha"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]