Message: Re: calculate particles within a cone Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Re: calculate particles within a cone 

Forum: Particles
Re: Question calculate particles within a cone (Summit)
Date: 21 May, 2014
From: Summit <Summit>

Hi

Sorry for posting in old thread!!!
I have a confusion related to this. 
Incident particles are moving in Z-direction. I have to score particles in cone whose half angle is 7.5 deg and having azimuth symmetry. In order to do this, 
      G4double Z = pd->GetPDGCharge()/eplus;
      G4double dirY = track->GetMomentumDirection().y(); 
      G4double dirZ = track->GetMomentumDirection().z(); 
      G4double theta = std::atan(dirY/dirZ);
      G4double Theta2 = track->GetMomentumDirection().theta();
      if(Z>0)
      if(-7.5<= theta/deg && theta/deg <= 7.5)
      histo->fill(2,Z,1.0);

I am not sure that am I doing right to calculate the (secondary) particles emitting from a desired cone?

I would also like to know what is the difference between theta and Theta2? And how track->GetMomentumDirection().phi() is different from this? From where can I get the complete documentation related to G4ThreeVector so that I can complete understand about all these.

Thanking you.

With kind regards, Summit.

 Add Message Add Message
to: "Re: calculate particles within a cone"

 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 ]