MJGeneratorTUNLFEL.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: MJGeneratorTUNLFEL.cc,v 1.2 2004/11/09 13:42:39 xliu Exp $ 
//      
// CLASS IMPLEMENTATION:  MJGeneratorTUNLFEL.cc
//
//---------------------------------------------------------------------------//
// 
//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//
//

#include <CLHEP/Random/RandGauss.h>

#include "G4ios.hh"
#include "G4Event.hh"
#include "G4Gamma.hh"
#include "G4ParticleGun.hh"
#include "G4Run.hh"
#include "G4ThreeVector.hh"

#include "generators/MJGeneratorTUNLFELMessenger.hh"
#include "io/MJLogger.hh"

//---------------------------------------------------------------------------//

#include "generators/MJGeneratorTUNLFEL.hh" 

//---------------------------------------------------------------------------//

MJGeneratorTUNLFEL::MJGeneratorTUNLFEL()
{
  fGeneratorName = "TUNLFEL";
  fG4Messenger = new MJGeneratorTUNLFELMessenger(this);
  fParticleGun = new G4ParticleGun(1);
  SetDefaults();
}

//---------------------------------------------------------------------------//

MJGeneratorTUNLFEL::MJGeneratorTUNLFEL(const MJGeneratorTUNLFEL & other)
{;}

//---------------------------------------------------------------------------//

MJGeneratorTUNLFEL::~MJGeneratorTUNLFEL()
{
  delete fG4Messenger;
  delete fParticleGun;
}

//---------------------------------------------------------------------------//

void MJGeneratorTUNLFEL::BeginOfRunAction(G4Run const *run)
{
  if(fMajorSigma < fMinorSigma){
    G4double temp = fMajorSigma;
    fMajorSigma = fMinorSigma;
    fMinorSigma = temp;
    MJLog(warning) << "Minor sigma > Major Sigma, values swapped." << endlog;
  } 
  if(MJLogger::GetSeverity() <= MJLogger::routine)
    Dump();
}  

//---------------------------------------------------------------------------//

void MJGeneratorTUNLFEL::Dump()
{
  G4cout << "  TUNL FEL Parameters :" << G4endl;
  G4cout << "  -----------------------------------------" << G4endl;
  G4cout << "  Rho (degrees): " << fRho/deg << G4endl;
  G4cout << "  Major sigma (cm): " << fMajorSigma/cm << G4endl;
  G4cout << "  Minor sigma (cm): " << fMinorSigma/cm << G4endl;
  G4cout << "  Mean Energy (MeV): " << fMeanEnergy/MeV << G4endl;
  G4cout << "  Energy sigma (keV): " << fEnergySigma/keV << G4endl;
  G4cout << "  Origin (cm) : " << fOrigin/cm << G4endl;
  G4cout << "-------------------------------------------" << G4endl << G4endl;
}

//---------------------------------------------------------------------------//

void MJGeneratorTUNLFEL::EndOfRunAction(G4Run const *run)
{;}

//---------------------------------------------------------------------------//

void MJGeneratorTUNLFEL::GeneratePrimaryVertex(G4Event *event)
{
  fParticleGun->SetParticleDefinition(G4Gamma::GammaDefinition());
  fParticleGun->SetParticleMomentumDirection(fDirection);
  fCurrentPosition.set(RandGauss::shoot(0.0, fMajorSigma), 
                       RandGauss::shoot(0.0, fMinorSigma), 0.0);
  fCurrentPosition.rotateZ(fRho);
  fCurrentPosition += fOrigin;
  fParticleGun->SetParticlePosition(fCurrentPosition);
  fCurrentEnergy = RandGauss::shoot(fMeanEnergy, fEnergySigma);
  fParticleGun->SetParticleEnergy(fCurrentEnergy);
  fParticleGun->GeneratePrimaryVertex(event);
}
  
//---------------------------------------------------------------------------//

void MJGeneratorTUNLFEL::SetDefaults()
{
  fDirection.set(0.0, 0.0, -1.0);
  fMajorSigma = fMinorSigma = 0.6 * cm;
  fOrigin.set(0.0, 0.0, 10.0 * m);
  fRho = 0.0;
  fMeanEnergy = 2.2 * MeV;
  fEnergySigma = 22.0 * keV;
}

//---------------------------------------------------------------------------//
//---------------------------------------------------------------------------//

---