GN_Data.cc

Go to the documentation of this file.

//! \file
//!
//! Source file for the GN_Data class.
//!
//! Defines the GN_Data class and the member routines which
//! are used in collecting the information for an event.
//!
//! \author D.K. Hasell
//! \version 1.0
//! \date 2010-10-31
//!
//! \ingroup detector

// Include the Data header file.

#include "GN_Data.h"

// Include the required GEANT header files.

#include "CLHEP/Units/PhysicalConstants.h"

// Include the STL header files referenced in this file.

#include <vector>

// Include the C++ header files referenced in this file.

#include <cmath>

// Use the STD namespace.

using namespace std;

// MCData class 

ClassImp( GN_Data )

// Constructor.

GN_Data::GN_Data() {}

// Destructor.

GN_Data::~GN_Data() {}

// Member functions.

// Routine to clear the event vectors.

void GN_Data::Reset() {

   nGN = 0;

   id.clear();
   q.clear();
   tr.clear();

   x.clear();
   y.clear();
   z.clear();

   px.clear();
   py.clear();
   pz.clear();

}

// Function to calculate the total momentum for a given particle.

double GN_Data::P( int i ) {
   return sqrt( px[i] * px[i] + py[i] * py[i] + pz[i] * pz[i] );
}

// Function to calculate the polar angle for a given particle.

double GN_Data::Theta( int i ) {
   return acos( pz[i] / P(i) );
}

// Function to calculate the azimuthal angle for a given particle.

double GN_Data::Phi( int i ) {
   double ph = atan2( py[i], px[i] );
   return ( ph > -CLHEP::halfpi ) ? ph : ph + CLHEP::twopi;
}