#include <EDepSimSpillTimeGenerator.hh>

Inheritance diagram for EDepSim::SpillTimeGenerator:

Public Member Functions
	SpillTimeGenerator (G4String name, double spillTime, double bunchSeparation, double bunchLength, const std::vector< double > &bunchPower)

virtual	~SpillTimeGenerator ()

double	GetTime (const G4LorentzVector &vertex)
	Return the time of the event to be generated. More...

virtual bool	ForceTime ()

Public Member Functions inherited from EDepSim::VTimeGenerator
	VTimeGenerator (const G4String &name)

virtual	~VTimeGenerator ()

G4String	GetName () const
	Return the name of the generator. More...

Private Attributes
double	fSpillTime
	The time that the neutrinos are crossing (0,0,0). More...

double	fBunchSeparation
	The time between bunch starts. More...

double	fBunchLength
	The length of a bunch. This is the gaussian width of the bunch. More...

std::vector< double >	fBunchPower
	The power in each bunch. More...

double	fMaxPower
	The maximum power in a bunch. More...

Detailed Description

Definition at line 13 of file EDepSimSpillTimeGenerator.hh.

Constructor & Destructor Documentation

EDepSim::SpillTimeGenerator::SpillTimeGenerator	(	G4String	name,
		double	spillTime,
		double	bunchSeparation,
		double	bunchLength,
		const std::vector< double > &	bunchPower
	)

Create the generatory. The name is for documentation, the spillTime gives the start time of the first bunch at the center of the EDepSim:: hall, the bunchSeparation is the time between bunch starts, and the bunchLength is the time length of a bunch with respect to it's start time.

Definition at line 9 of file EDepSimSpillTimeGenerator.cc.

     : EDepSim::VTimeGenerator(name), fSpillTime(spillTime),
       fBunchSeparation(bunchSeparation), fBunchLength(bunchLength),
       fBunchPower(bunchPower) {
     if (fBunchPower.size()<1) fBunchPower.push_back(1.0);
     fMaxPower = 0.0;
     for (std::vector<double>::iterator b = fBunchPower.begin();
          b != fBunchPower.end();
          ++b) {
         if (fMaxPower<(*b)) fMaxPower = *b;
     }
 
     if (fMaxPower<0.001) {
         EDepSimError("Spill created without any power in bunchs");
     }
 }

EDepSim::SpillTimeGenerator::~SpillTimeGenerator ( )

virtual

Definition at line 28 of file EDepSimSpillTimeGenerator.cc.

28 {}

Member Function Documentation

bool EDepSim::SpillTimeGenerator::ForceTime ( )

virtual

Flag if the time should be forced. This tells the calling routine that the time returned by this generator should override any times found in the original vertex. This returns true, so that the generated spill time overrides the interaction time.

Implements EDepSim::VTimeGenerator.

Definition at line 54 of file EDepSimSpillTimeGenerator.cc.

                                           {
     return true;
 }

double EDepSim::SpillTimeGenerator::GetTime ( const G4LorentzVector & vertex )

virtual

Return the time of the event to be generated.

Implements EDepSim::VTimeGenerator.

Definition at line 30 of file EDepSimSpillTimeGenerator.cc.

                                                                     {
     // Choose a spill;
     int spills = fBunchPower.size();
     int spill = 0;
     do {
         spill = int(spills*G4UniformRand());
         if (spill>=spills) spill=spills-1;
         if (spill<0) spill = 0;
     } while (fBunchPower[spill]<fMaxPower*G4UniformRand());
 
     // Get the spill start time.  This assumes the neutrino wave is traveling
     // along the positive Z axis.
     double time = fSpillTime + vtx.z()/c_light + spill*fBunchSeparation;
 
     // Adjust the time in the bunch.
     double truncTime = G4RandGauss::shoot(0.0, 1.0);
     while (truncTime < -2.0 || truncTime > 2.0) {
         truncTime = G4RandGauss::shoot(0.0, 1.0);
     }
     time += fBunchLength*truncTime;
 
     return time;
 }