Message: Re: Spallation events identification Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Spallation events identification 

Forum: Hadronic Processes
Re: Question Spallation events identification (Stefania stucci)
Re: None Re: Spallation events identification (Michael H. Kelsey)
Re: None Re: Spallation events identification (Stefania stucci)
Re: None Re: Spallation events identification (Michael H. Kelsey)
Re: None Re: Spallation events identification (Stefania stucci)
Date: Mar 04, 03:00
From: michel maire <michel maire>

On Sun, 03 Mar 2019 05:36:22 GMT, Stefania Antonia Stucci wrote:
> Hi, still me. 
> I noticed that by using this code: 
> 
> if (aStep->GetTrack()->GetParentID()!=0) {
>         particle= aStep->GetTrack()->GetDefinition();
>         if (particle->GetParticleType() == "nucleus") {     
>             G4int A = particle->GetAtomicMass();
>             analysisManager->FillH1(1, A);
>         }        
>     }
> 
> instead of 
> 
> for (size_t lp=0; lp<(*secondary).size(); lp++) {
>             particle = (*secondary)[lp]->GetDefinition();
>             G4String name   = particle->GetParticleName();
>             if (particle->GetParticleType() == "nucleus") {
>        G4int A = particle->GetAtomicMass();
>        analysisManager->FillNtupleDColumn(1, A);       
>           }       
>     }
> 
> gives me two different numbers, in particular the first method gives me a larger number of ions. Following event by event I noticed that the first method is counting a larger number of deuterons than actually geant4 produces as secondary. Do you know why there is this double counting? What’s wrong with the first method? 
> 

  With the second method you scan the list of secondaries created when a nuclear interaction occurs (therefore, only once).

 With the first method, you check the current particle being tracked, step by step. Usually, a particle makes several steps in a given volume.

                         Michel

 Add Message Add Message
to: "Re: Spallation events identification"

 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 ]