MJGeneratorUtil.cc

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu//
//                                                                           //
//                         MAJORANA Simulation                               //
//                                                                           //
//      This code implementation is the intellectual property of the         //
//      MAJORANA Collaboration. It is based on Geant4, an intellectual       //
//      property of the RD44 GEANT4 collaboration.                           //
//                                                                           //
//                        *********************                              //
//                                                                           //
//    Neither the authors of this software system, nor their employing       //
//    institutes, nor the agencies providing financial support for this      //
//    work  make  any representation or  warranty, express or implied,       //
//    regarding this software system or assume any liability for its use.    //
//    By copying, distributing or modifying the Program (or any work based   //
//    on on the Program) you indicate your acceptance of this statement,     //
//    and all its terms.                                                     //
//                                                                           //
//bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu//
//---------------------------------------------------------------------------//
//
// $Id: MJGeneratorUtil.cc,v 1.1 2004/11/22 13:07:41 gast Exp $
//
// CLASS IMPLEMENTATION:  MJGeneratorUtil.cc
//---------------------------------------------------------------------------//
// --------------------------------------------------------------------------//

#include "generators/MJGeneratorUtil.hh"

MJGeneratorUtil::MJGeneratorUtil()
{;}

MJGeneratorUtil::~MJGeneratorUtil()
{;}

G4ThreeVector MJGeneratorUtil::pick_isotropic() 
{
  G4double phi       = twopi * G4UniformRand();
  G4double cos_theta = -1.0 + 2.0*G4UniformRand();
  G4double theta     = acos(cos_theta);
  G4double sin_theta = sin(theta);
  G4double cos_x     = sin_theta*cos(phi);
  G4double cos_y     = sin_theta*sin(phi);
  G4double cos_z     = cos_theta;

  G4ThreeVector direction(cos_x,cos_y,cos_z);

  return direction;
}

G4ThreeVector MJGeneratorUtil::pick_point_in_box(G4double x_lo, G4double x_hi,
                                G4double y_lo, G4double y_hi,
                                G4double z_lo, G4double z_hi) {
  G4double dx = x_hi - x_lo;
  G4double dy = y_hi - y_lo;
  G4double dz = z_hi - z_lo;

  G4double x = x_lo + dx*G4UniformRand();
  G4double y = y_lo + dy*G4UniformRand();
  G4double z = z_lo + dz*G4UniformRand();

  G4ThreeVector pos(x,y,z);

  return pos;
}
 
//*************************************************************************
void MJGeneratorUtil::pick_point_in_circle(G4double Radius, G4double &x, G4double &y) {

  G4double b     = Radius*Radius;
  G4double x_rnd = G4UniformRand();
  G4double r     = sqrt(b*x_rnd);
  G4double phi   = twopi * G4UniformRand();

  x = r*cos(phi);
  y = r*sin(phi);
}

//*************************************************************************
G4ThreeVector MJGeneratorUtil::pick_point_in_annulus(G4double r1, G4double r2, G4double h) 
{
  G4double a = r1*r1;
  G4double b = r2*r2 - a;

  G4double x_rnd = G4UniformRand();
  G4double r = sqrt(b*x_rnd + a);

  G4double phi = twopi * G4UniformRand();

  G4double x = r * cos(phi);
  G4double y = r * sin(phi);

  G4double z = h * (-0.5 + G4UniformRand());

  G4ThreeVector pos(x,y,z);

  return pos;

}
//*************************************************************************
G4ThreeVector MJGeneratorUtil::pick_point_in_annulus(G4double r1, G4double r2, G4double h,
                                    G4double theta0, G4double dtheta) {

  G4double a = r1*r1;
  G4double b = r2*r2 - a;

  G4double x_rnd = G4UniformRand();
  G4double r = sqrt(b*x_rnd + a);

  G4double phi = theta0 + dtheta * G4UniformRand();

  G4double x = r * cos(phi);
  G4double y = r * sin(phi);

  G4double z = h * (-0.5 + G4UniformRand());

  G4ThreeVector pos(x,y,z);

  return pos;

}

//*************************************************************************
G4ThreeVector MJGeneratorUtil::pick_point_in_cylinder(G4double R, G4double L) {

  return pick_point_in_annulus(0.0, R, L);

}

//*************************************************************************
G4ThreeVector MJGeneratorUtil::pick_point_on_cylinder(G4double R, G4double L) 
{
  G4double phi = twopi * G4UniformRand();
  G4double x = R*cos(phi);
  G4double y = R*sin(phi);

  G4double z = L * (-0.5 + G4UniformRand());

  G4ThreeVector pos(x,y,z);

  return pos;
}


//*************************************************************************
G4ThreeVector MJGeneratorUtil::pick_point_in_shell(G4double r1, G4double r2) 
{
  const G4double ONE_THIRD = 1.0/3.0;

  G4double r1_3 = pow(r1,3);
  G4double r2_3 = pow(r2,3);
  G4double d    = r2_3 - r1_3;

  G4double x_rnd = r1_3 + d * G4UniformRand();
  G4double r = pow(x_rnd,ONE_THIRD);

  return r*pick_isotropic();
}

//*************************************************************************
G4ThreeVector MJGeneratorUtil::pick_point_in_sphere(G4double r) 
{
  return pick_point_in_shell(0.0, r);
}

//*************************************************************************
G4ThreeVector MJGeneratorUtil::pick_point_on_sphere(G4double r) 
{
  return r*pick_isotropic();
}

---