EDepSim::NuMIRockKinematicsGenerator Class Reference

#include <EDepSimNuMIRockKinematicsGenerator.hh>

Inheritance diagram for EDepSim::NuMIRockKinematicsGenerator:

Public Member Functions
	NuMIRockKinematicsGenerator (const G4String &name, const G4String &fluxName)
virtual	~NuMIRockKinematicsGenerator ()
virtual GeneratorStatus	GeneratePrimaryVertex (G4Event *evt, const G4LorentzVector &position)
	Add a primary vertex to the event. More...
Public Member Functions inherited from EDepSim::VKinematicsGenerator
	VKinematicsGenerator (const G4String &name)
virtual	~VKinematicsGenerator ()
G4String	GetName () const
	Return the name of the generator. More...

Private Attributes
std::string	fFluxName
	The static part of the file name field. More...

Additional Inherited Members
Public Types inherited from EDepSim::VKinematicsGenerator
enum	GeneratorStatus { kFail = 0, kSuccess, kLastVertex, kEndEvent }
	A status value that can be returned by GeneratePrimaryVertex. More...

Detailed Description

Definition at line 12 of file EDepSimNuMIRockKinematicsGenerator.hh.

Constructor & Destructor Documentation

EDepSim::NuMIRockKinematicsGenerator::NuMIRockKinematicsGenerator	(	const G4String &	name,
		const G4String &	fluxName
	)

Construct a new generator. The name is the name of the generator (e.g. GENIE, rock, &c). The fluxName is the name of the neutrino spectrum.

Definition at line 23 of file EDepSimNuMIRockKinematicsGenerator.cc.

    : EDepSim::VKinematicsGenerator(name), fFluxName(fluxName) { }

EDepSim::NuMIRockKinematicsGenerator::~NuMIRockKinematicsGenerator ( )

virtual

Definition at line 27 of file EDepSimNuMIRockKinematicsGenerator.cc.

{ }

Member Function Documentation

EDepSim::VKinematicsGenerator::GeneratorStatus EDepSim::NuMIRockKinematicsGenerator::GeneratePrimaryVertex ( G4Event *  evt, const G4LorentzVector &  position  )

virtual

Add a primary vertex to the event.

Implements EDepSim::VKinematicsGenerator.

Definition at line 30 of file EDepSimNuMIRockKinematicsGenerator.cc.

                            {

    // Make a position vertex...
    double boxSize = 300*cm;
    G4LorentzVector eventVertex;
    for (int i=0; i<3; ++i) {
        eventVertex[i] = 2*boxSize*G4UniformRand() - boxSize;
    }
    eventVertex[0] = -3.0*meter;
    eventVertex[2] -= 0.5*meter;
    eventVertex[3] = 10.0*microsecond*G4UniformRand();
    
    // Create a new vertex to add the new particles, and add the vertex to the
    // event.
    G4PrimaryVertex* theVertex 
        = new G4PrimaryVertex(G4ThreeVector(eventVertex.x(),
                                            eventVertex.y(),
                                            eventVertex.z()),
                              eventVertex.t());
    anEvent->AddPrimaryVertex(theVertex);

    // Fill the particles to be tracked (status ==1).  These particles are
    // attached to the primary vertex.  Also save the incident neutrino
    // particle and the incident target nucleus; these particles are attached
    // to informational vertex.
    G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
    G4ParticleDefinition* particleDef = particleTable->FindParticle("mu-");

    // Get the momentum.
    double momentum;;
    do {
        momentum = G4RandGauss::shoot(4.0*GeV,2*GeV);
    } while (momentum<2*GeV);
    momentum *= G4UniformRand();

    G4ThreeVector dir(1.0,
                      G4RandGauss::shoot(0.0,0.15),
                      G4RandGauss::shoot(0.0,0.15));
    dir = dir.unit();

    std::cout << "Add rock muon at "
              << " " << eventVertex.x()
              << " " << eventVertex.y()
              << " " << eventVertex.z()
              << " " << eventVertex.t()
              << " going " << dir.x()
              << " " << dir.y()
              << " " << dir.z()
              << std::endl;
        
    // Create the particle.
    G4PrimaryParticle* theParticle
        = new G4PrimaryParticle(particleDef,
                                momentum*dir.x(),
                                momentum*dir.y(),
                                momentum*dir.z());
    theVertex->SetPrimary(theParticle);

    return kSuccess;
}

Member Data Documentation

std::string EDepSim::NuMIRockKinematicsGenerator::fFluxName

private

The static part of the file name field.

Definition at line 27 of file EDepSimNuMIRockKinematicsGenerator.hh.

---

The documentation for this class was generated from the following files:

• edep-sim/src/kinem/EDepSimNuMIRockKinematicsGenerator.hh
• edep-sim/src/kinem/EDepSimNuMIRockKinematicsGenerator.cc