Message: Re: Double first step in the volume Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Feedback Re: Double first step in the volume 

Forum: Event and Track Management
Re: None Double first step in the volume
Re: Feedback Re: Double first step in the volume (Gumplinger Peter)
Re: None Re: Double first step in the volume
Date: 19 Nov, 2015
From: Gumplinger Peter <Gumplinger Peter>

> Step# X Y Z KineE dEStep StepLeng TrakLeng Volume Process
> 0 -2.35 mm 8.03 mm -10 cm 8 MeV 0 eV 0 fm 0 fm World initStep
> 1 -2.35 mm 8.03 mm -2.25 mm 8 MeV 0 eV 9.78 cm 9.78 cm World
> Transportation
> 2 -2.35 mm 8.03 mm 2.25 mm 8 MeV 0 eV 4.5 mm 10.2 cm physicTarget
> Transportation
> 3 -2.35 mm 8.03 mm 3.81 cm 8 MeV 0 eV 3.58 cm 13.8 cm World
> Transportation
> Stepnumber: 3 Particle: gamma ProcessName: Transportation

So, you print for the step that arrives at the Crystall. It is NOT THE FIRST STEP INSIDE CRYSTALL! The step is inside World. The preStepPoint has fGeomBoundary because it so happens that is also started where the track entered World.

> 4 -2.35 mm 8.03 mm 4.85 cm 8 MeV 0 eV 1.04 cm 14.9 cm Crystall
> Transportation

> 2.
> Step# X Y Z KineE dEStep StepLeng TrakLeng Volume Process
> 0 3.21 mm 6.07 mm -10 cm 8 MeV 0 eV 0 fm 0 fm World initStep
> 1 3.21 mm 6.07 mm -2.25 mm 8 MeV 0 eV 9.78 cm 9.78 cm World
> Transportation
> 2 3.21 mm 6.07 mm 2.25 mm 8 MeV 0 eV 4.5 mm 10.2 cm physicTarget
> Transportation
> 3 3.21 mm 6.07 mm 4.15 cm 8 MeV 0 eV 3.92 cm 14.1 cm World
> Transportation
> Stepnumber: 3 Particle: gamma ProcessName: Transportation

Same as above.

> 4 3.21 mm 6.07 mm 4.26 cm 533 keV 0 eV 1.15 mm 14.3 cm Crystall compt
> Stepnumber: 4 Particle: gamma ProcessName: compt

Yes, fGeomBoundary is where preStepPoint entered Crystall and since the step ends inside Crystall the postname == "Crystall" is also true.

> 5 8.71 mm 1.24 cm 4.35 cm 486 keV 0 eV 8.42 mm 15.1 cm Crystall compt

> Looks like if interaction is happening in the Crystall step will still
> satisfy condition
> 
>  if (thePrePoint ->GetStepStatus()==fGeomBoundary && postname ==
> "Crystall")

Yes, so long as this is the first step inside Crystall.

> Does the prestep Status is still fGeomBoundary after transpostation has
> occured from another volume?

Your logic prints for the step arriving at Crystall if it crossed another volume and for the first step inside Crystall in case the step ends inside the Crystall.

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

1 Agree: Re: Double first step in the volume   (Dmytro - 20 Nov, 2015)
 Add Message Add Message
to: "Re: Double first step in the volume"

 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 ]