Message: Re: GPS bias command Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Re: GPS bias command 

Forum: Event and Track Management
Re: Question GPS bias command (JungEun Elena Suh)
Re: More Re: GPS bias command (JungEun Elena Suh)
Date: 11 Sep, 2015
From: JungEun Elena Suh <JungEun Elena Suh>

So for example
If I saw GPS Command Example 24:

  Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2015/09/11/07.29-95592-test24.gif

    /gps/particle geantino
    /gps/pos/type Volume
    /gps/pos/shape Sphere
    /gps/pos/centre 1. 2. 1. cm
    /gps/pos/radius 2. cm
    /gps/ang/type iso
    /gps/ene/type Arb
    /gps/ene/diffspec 0
    /gps/hist/type arb
    /gps/hist/point 0.0 11.
    /gps/hist/point 1.0 10.
    /gps/hist/point 2.0 9.
    /gps/hist/point 3.0 8.
    /gps/hist/point 4.0 7.
    /gps/hist/point 7.0 4.
    /gps/hist/point 8.0 3.
    /gps/hist/point 9.0 2.
    /gps/hist/point 10.0 1.
    /gps/hist/point 11.0 0.
    /gps/hist/inter Lin
    /gps/hist/type biasz
    /gps/hist/point 0. 0.
    /gps/hist/point 0.4 0.5
    /gps/hist/point 0.6 1.
    /gps/hist/point 1. 0.2
    /gps/hist/type biast
    /gps/hist/point 0. 0.
    /gps/hist/point 0.1 1.
    /gps/hist/point 0.5 0.1
    /gps/hist/point 1. 1.
    /gps/hist/type biasp
    /gps/hist/point 0. 0.
    /gps/hist/point 0.125 1.
    /gps/hist/point 0.375 4.
    /gps/hist/point 0.625 1.
    /gps/hist/point 0.875 4.
    /gps/hist/point 1. 1.



Here, as I understood,

1. Energy

/gps/hist/type arb
/gps/hist/point 0.0 11.  : lower boundary of the first bin? why weight 11?? (0 MeV)
/gps/hist/point 1.0 10.  : 0~1 MeV, weight 10
/gps/hist/point 2.0 9.   : 1~2 MeV, weight 9 ...
/gps/hist/point 3.0 8.


2. Theta

/gps/hist/type biast
/gps/hist/point 0. 0.   : lower boundary of the first bin (0 degree)
/gps/hist/point 0.1 1.  : 0~18 deg, weight 1
/gps/hist/point 0.5 0.1 : 18~90 deg, weight 0.1
/gps/hist/point 1. 1.   : 90~180 deg, weight 1

3. BiasZ

/gps/hist/type biasz
/gps/hist/point 0. 0.   : lower boundary of the first bean (-1 cm)
/gps/hist/point 0.4 0.5 : -1 ~ 0.6 cm (z), weight 0.5  (because center of the sphere is at 1 2 1 cm and the radius is 2 cm)
/gps/hist/point 0.6 1.  : 0.6 ~ 1.4 cm (z), weight 1
/gps/hist/point 1. 0.2  : 1.4 ~ 3 cm (z), weight 0.2


Is that right?

I think(hope) it is working like this,
so I put biasy in my gps macro.

This is my gps macro

/gps/particle ion
/gps/ion 4 9

/gps/ene/type Mono
/gps/ene/mono 20 GeV

/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/pos/centre 0 0 0 cm
/gps/pos/radius 1 m

/gps/ang/type iso

/gps/ang/minphi 0 deg
/gps/ang/maxphi 0 deg
/gps/ang/mintheta 0 deg
/gps/ang/maxtheta 0 deg

(I'll expand maxphi to 360 deg and maxtheta to certain angle later. I wish it doesn't cause any problem now.)

And I added 4 lines,

/gps/hist/type biasy
/gps/hist/point 0. 0. 
/gps/hist/point 0.9 0.
/gps/hist/point 1. 1

And I hope it doesn't initiate any ions between -1m to 0.8m along y axis and only initiate ions between 0.8 to 1 m along y axis.
Am I wrong?
It doesn't show me any difference if I put those 4 lines or not.

I draw initial x y z position of primary particle and judging from this.
I got initial x y z position of primary particle by

  G4PrimaryVertex* primaryVertex = evt->GetPrimaryVertex();
  G4PrimaryParticle* primaryParticle = primaryVertex->GetPrimary();
  G4double vertx = primaryVertex->GetX0();
  G4double verty = primaryVertex->GetY0();
  G4double vertz = primaryVertex->GetZ0();

in EventAction.cc


Am I wrong in some point?
or completely wrong?


Ah...and when I load this macro after I run my project,
it shows me a message
-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : PART11117
      issued by : G4ParticleTable::FindIon()
This method is obsolete and will be dropped from v10.0. Use G4IonTable::FindIon().
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------

But I never used FindIon() anywhere...
Can it be the problem...?


Thank you so much for the reading.
I'm going to investigate this problem again... :(




 

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

1 None: Re: GPS bias command   (JungEun Elena Suh - 06 Oct, 2015)
 Add Message Add Message
to: "Re: GPS bias command"

 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 ]