GFluxWindow.cxx

Go to the documentation of this file.

// ROOT includes
#include "TVector3.h"

// GENIE includes
#include "GENIE/Framework/Messenger/Messenger.h"
#include "GENIE/Framework/Numerical/RandomGen.h"
#include "GENIE/Framework/Conventions/Units.h"

// NOvASoft includes
#include "GenieSNova/src/GFluxWindow.h"

using namespace genie::supernova;


//.......................................................................
GFluxWindow::GFluxWindow(TVector3 pos, double size):
            fPosition(pos),
            fDirection(-pos.Unit()),
            fSize(size)
{
    CalculateBasis();
}


//.......................................................................
void GFluxWindow::GenerateNext(GenParticle& particle)
{
    particle.SetDirection(fDirection);
    particle.fPos.SetVect(GenPosition() * 0.01);
}


//.......................................................................
TVector3 GFluxWindow::Direction()
{
    return fDirection;
}


//.......................................................................
double GFluxWindow::Area()
{
    return fSize * fSize;
}


//.......................................................................
TVector3 GFluxWindow::GenPosition()
{
    genie::RandomGen * rnd = genie::RandomGen::Instance();
    double x = rnd->RndFlux().Uniform(-0.5, 0.5);
    double y = rnd->RndFlux().Uniform(-0.5, 0.5);
    return (fBasisX * x) + (fBasisY * y) + fPosition;
}


//.......................................................................
void GFluxWindow::CalculateBasis()
{
    // Calculate Basis vectors
    fBasisX = fDirection.Cross(TVector3(0, 1, 0));
    if (fBasisX.Mag2() == 0)
        fBasisX = fDirection.Cross(TVector3(0, 0, 1));

    fBasisY = fBasisX.Cross(fDirection);

    //normalize basis
    fBasisX = fBasisX.Unit() * fSize;
    fBasisY = fBasisY.Unit() * fSize;
}