MaGeGeneratorCosmicRayMuons.hh

Go to the documentation of this file.

//
// ********************************************************************
// * DISCLAIMER                                                       *
// *                                                                  *
// * 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.                                                             *
// *                                                                  *
// ********************************************************************
//
// This is supposed to be a generator for cosmic ray muons (with the 
// correct spectrum) in the tunnel. It is work in progress
//
// History:
// --------
// 28 Oct 2004   L.Pandola    First implementation (not working)

#ifndef MaGeGeneratorCosmicRayMuons_h
#define MaGeGeneratorCosmicRayMuons_h 1

#include "G4VPrimaryGenerator.hh"
#include "G4Navigator.hh"
#include "G4ParticleMomentum.hh"
#include "G4ParticleDefinition.hh"
#include "globals.hh"
#include <fstream>

#include "generators/MaGeGeneratorCosmicRayMuonsMessenger.hh"
class G4DataVector;

class MaGeGeneratorCosmicRayMuons : public G4VPrimaryGenerator {

public:
  MaGeGeneratorCosmicRayMuons (G4String,G4String); 
  ~MaGeGeneratorCosmicRayMuons ();
  void GeneratePrimaryVertex(G4Event *evt);

public:

  // position distribution  
  void SetHalfZ(G4double);
  void SetRadius(G4double);

private: 
  // angular distribution
  void GenerateAngularSpectrum();

  // energy distribution 
  void GenerateEnergySpectrum();

  // initial position
  void SampleInitialPosition();

  void OpenEnergyFile(G4String);
  void OpenAngularFile(G4String);

public:
  //Get methods
  inline G4double GetParticleEnergy() {return particle_energy;}
  inline G4String GetParticleName() {return particle_definition->GetParticleName();};

  // verbosity
  void SetVerbosity(G4int);
  
  //Set methods
  void SetParticleDefinition(G4ParticleDefinition * aParticleDefinition);
  inline void SetParticleCharge(G4double aCharge)
  { particle_charge = aCharge; }
  
private:
  // position distribution
  G4double halfz;
  G4double Radius;
  
  // particle properties 
  G4int                  NumberOfParticlesToBeGenerated;
  G4ParticleDefinition*  particle_definition;
  G4ParticleMomentum     particle_momentum_direction;
  G4double               particle_energy;
  G4double               particle_charge;
  G4ThreeVector          particle_position;
  G4double               particle_time;
  G4ThreeVector          particle_polarization;

  //files to be read
  G4String fFileEnergyName;
  G4String fFileAngularName;
  std::ifstream fEnergyInputFile;
  std::ifstream fAngularInputFile;

  // Verbose
  G4int verbosityLevel;

private:
  G4DataVector* fEnergy;
  G4DataVector* fEnergySpectrum;
  MaGeGeneratorCosmicRayMuonsMessenger *theMessenger;
  G4Navigator *gNavigator;
  
};


#endif

---