CML2SDWithParticle Class Reference

#include <ML2SDWithParticle.hh>

Inheritance diagram for CML2SDWithParticle:

Collaboration diagram for CML2SDWithParticle:

Public Member Functions
	CML2SDWithParticle ()
	CML2SDWithParticle (G4int idType, G4int max_N_particles_in_PhSp_File, G4int seed, G4int nMaxParticlesInRamPhaseSpace, G4String name, G4String PhaseSpaceOutFile, G4bool bSavePhaseSpace, G4bool bStopAtVolatilePhaseSpace, SPrimaryParticle *primaryParticleData, G4double accTargetZPosition)
	~CML2SDWithParticle (void)
G4bool	ProcessHits (G4Step *aStep, G4TouchableHistory *ROHist)
G4int	getTotalNumberOfParticles ()
CML2SDWithParticle *	getCML2SensitiveDetectorParticle ()
Sparticle	getParticle (int i)
void	setActive (G4bool act)
void	save ()
Public Member Functions inherited from G4VSensitiveDetector
	G4VSensitiveDetector (G4String name)
	G4VSensitiveDetector (const G4VSensitiveDetector &right)
virtual	~G4VSensitiveDetector ()
G4VSensitiveDetector &	operator= (const G4VSensitiveDetector &right)
G4int	operator== (const G4VSensitiveDetector &right) const
G4int	operator!= (const G4VSensitiveDetector &right) const
virtual void	Initialize (G4HCofThisEvent *)
virtual void	EndOfEvent (G4HCofThisEvent *)
virtual void	clear ()
virtual void	DrawAll ()
virtual void	PrintAll ()
G4bool	Hit (G4Step *aStep)
void	SetROgeometry (G4VReadOutGeometry *value)
void	SetFilter (G4VSDFilter *value)
G4int	GetNumberOfCollections () const
G4String	GetCollectionName (G4int id) const
void	SetVerboseLevel (G4int vl)
void	Activate (G4bool activeFlag)
G4bool	isActive () const
G4String	GetName () const
G4String	GetPathName () const
G4String	GetFullPathName () const
G4VReadOutGeometry *	GetROgeometry () const
G4VSDFilter *	GetFilter () const
virtual G4VSensitiveDetector *	Clone () const

Private Member Functions
void	saveDataParticles (G4int nParticle)
void	saveHeaderParticles ()

Private Attributes
G4ThreeVector	halfSize
G4ThreeVector	pos
SPrimaryParticle *	primaryParticleData
Sparticle *	particles
G4String	fullOutFileData
G4int	nTotalParticles
G4int	nParticle
G4int	idType
G4int	nMaxParticlesInRamPhaseSpace
G4int	max_N_particles_in_PhSp_File
G4bool	bStopAtPhaseSpace
G4bool	bSavePhaseSpace
G4bool	bActive
G4double	accTargetZPosition

Additional Inherited Members
Protected Member Functions inherited from G4VSensitiveDetector
virtual G4int	GetCollectionID (G4int i)
Protected Attributes inherited from G4VSensitiveDetector
G4CollectionNameVector	collectionName
G4String	SensitiveDetectorName
G4String	thePathName
G4String	fullPathName
G4int	verboseLevel
G4bool	active
G4VReadOutGeometry *	ROgeometry
G4VSDFilter *	filter

Detailed Description

Definition at line 55 of file ML2SDWithParticle.hh.

Constructor & Destructor Documentation

CML2SDWithParticle() [1/2]

CML2SDWithParticle::CML2SDWithParticle

(

)

Definition at line 48 of file ML2SDWithParticle.cc.

: G4VSensitiveDetector("killer_plane"),particles(0)
{
    idType=idSD_KillerPlane;
    bStopAtPhaseSpace=true;
    nTotalParticles=0;
    bActive=true;
}

CML2SDWithParticle() [2/2]

CML2SDWithParticle::CML2SDWithParticle	(	G4int	idType,
		G4int	max_N_particles_in_PhSp_File,
		G4int	seed,
		G4int	nMaxParticlesInRamPhaseSpace,
		G4String	name,
		G4String	PhaseSpaceOutFile,
		G4bool	bSavePhaseSpace,
		G4bool	bStopAtVolatilePhaseSpace,
		SPrimaryParticle *	primaryParticleData,
		G4double	accTargetZPosition
	)

Definition at line 57 of file ML2SDWithParticle.cc.

: G4VSensitiveDetector(name),particles(0)
{
    accTargetZPosition=ZPos;
    max_N_particles_in_PhSp_File=maxPartFile;
    nMaxParticlesInRamPhaseSpace = nMaxPart;
    idType=id;
    primaryParticleData=pData;
    bActive=true;
    nParticle=0;
    nTotalParticles=0;
    bSavePhaseSpace=bSave;
    bStopAtPhaseSpace=bStop;
    if (bSavePhaseSpace)
    {particles=new Sparticle[nMaxPart];}
    G4String seedName;
    char a[10];
    sprintf(a,"%d", seed);
    seedName=(G4String)a;

    fullOutFileData=PhaseSpaceOutFile+"_"+seedName+".txt";
    fullOutFileData=PhaseSpaceOutFile+"_"+seedName+".txt";
}

~CML2SDWithParticle()

CML2SDWithParticle::~CML2SDWithParticle

(

void

)

Definition at line 83 of file ML2SDWithParticle.cc.

{
        if (particles!=0)
    {
        delete [] particles;
    }
}

Member Function Documentation

getCML2SensitiveDetectorParticle()

CML2SDWithParticle* CML2SDWithParticle::getCML2SensitiveDetectorParticle ( )

inline

Definition at line 63 of file ML2SDWithParticle.hh.

{return this;}

Here is the caller graph for this function:

getParticle()

Sparticle CML2SDWithParticle::getParticle ( int  i )

inline

Definition at line 64 of file ML2SDWithParticle.hh.

{return particles[i];}

getTotalNumberOfParticles()

G4int CML2SDWithParticle::getTotalNumberOfParticles ( )

inline

Definition at line 62 of file ML2SDWithParticle.hh.

{return nTotalParticles;}

Here is the caller graph for this function:

ProcessHits()

G4bool CML2SDWithParticle::ProcessHits ( G4Step *  aStep, G4TouchableHistory *  ROHist  )

virtual

Implements G4VSensitiveDetector.

Definition at line 121 of file ML2SDWithParticle.cc.

{
    if (bActive &&  (CML2AcceleratorConstruction::GetInstance()->getPhysicalVolume()->GetRotation()->isIdentity()))
    {
        G4double energyKin= aStep->GetTrack()->GetKineticEnergy();
        static bool bFirstTime=true;
        if (idType==idSD_KillerPlane)
        {
            nTotalParticles++;
            aStep->GetTrack()->SetTrackStatus(fStopAndKill);
        }
        else
        {
            if (energyKin>0.)
            {
                particles[nParticle].volumeName="";
                particles[nParticle].pos=aStep->GetPreStepPoint()->GetPosition();
                particles[nParticle].dir=aStep->GetPreStepPoint()->GetMomentumDirection();
                particles[nParticle].kinEnergy=energyKin;
                particles[nParticle].partPDGE=aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
                particles[nParticle].primaryParticlePDGE=primaryParticleData->partPDGE;
                particles[nParticle].nPrimaryPart=primaryParticleData->nPrimaryParticle;
                nParticle++;
                nTotalParticles++;
                if (nTotalParticles==max_N_particles_in_PhSp_File)
                {
                    if (bFirstTime)
                    {
                        bFirstTime=false;
                        saveHeaderParticles();
                    }
                    saveDataParticles(nParticle);
                    nParticle=0;
                    bActive =false;// to stop the phase space creation
                }
                if (nParticle==nMaxParticlesInRamPhaseSpace)
                {
                    if (bFirstTime)
                    {
                        bFirstTime=false;
                        saveHeaderParticles();
                    }
                    saveDataParticles(nParticle);
                    nParticle=0;
                }

                Sparticle *particle=new Sparticle;
                particle->dir=aStep->GetPreStepPoint()->GetMomentumDirection();
                particle->pos=aStep->GetPreStepPoint()->GetPosition();
                particle->kinEnergy=energyKin; 
                particle->nPrimaryPart=nTotalParticles; // to pass the id of this phase space particle 
                particle->partPDGE=aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
                particle->primaryParticlePDGE=primaryParticleData->partPDGE;
                particle->volumeId=-1;
                particle->volumeName="-1";
            }
            if (bStopAtPhaseSpace)
            {aStep->GetTrack()->SetTrackStatus(fStopAndKill);}
        }
    }
    else
    {
        if (bStopAtPhaseSpace)
        {aStep->GetTrack()->SetTrackStatus(fStopAndKill);}
    }

    return true;
}

Here is the call graph for this function:

save()

void CML2SDWithParticle::save

(

)

Definition at line 190 of file ML2SDWithParticle.cc.

{
    if ((bActive) && (nParticle>0) && (CML2AcceleratorConstruction::GetInstance()->getPhysicalVolume()->GetRotation()->isIdentity()))
    {saveDataParticles(nParticle);nParticle=0;}
}

Here is the call graph for this function:

Here is the caller graph for this function:

saveDataParticles()

void CML2SDWithParticle::saveDataParticles ( G4int  nParticle )

private

Definition at line 98 of file ML2SDWithParticle.cc.

{
    std::ofstream out;
    out.open(fullOutFileData, std::ios::app);
    static G4int nTotParticles=0;
    for (int i=0; i< nPart; i++)
    {
        out << nTotParticles++ << '\t';
//      out << particles[i].volumeName << '\t';
        out << particles[i].pos.getX()/mm  << '\t';
        out << particles[i].pos.getY()/mm  << '\t';
        out << (accTargetZPosition + particles[i].pos.getZ())/mm  << '\t'; // it translates the current z value in global coordinates to the accelerator local coordinates (only z)
        out << particles[i].dir.getX() << '\t';
        out << particles[i].dir.getY() << '\t';
        out << particles[i].dir.getZ() << '\t';
        out << particles[i].kinEnergy/MeV << '\t';
        out << particles[i].partPDGE<< '\t';
        out << particles[i].primaryParticlePDGE<< '\t';
        out << particles[i].nPrimaryPart<< G4endl;
    }
    out.close();
}

Here is the call graph for this function:

Here is the caller graph for this function:

saveHeaderParticles()

void CML2SDWithParticle::saveHeaderParticles ( )

private

Definition at line 90 of file ML2SDWithParticle.cc.

{
    std::ofstream out;
    out.open(fullOutFileData, std::ios::out);
    out << "Sensitive Detector-Particles"<<G4endl;
    out << "n Total Events,\t x [mm],\t y [mm],\t z [mm],\t dirX,\t dirY,\t dirZ,\t KinEnergy [MeV],\t part Type,\t primary part type,\t nPrimaryPart" << G4endl;
    out.close();
}

Here is the caller graph for this function:

setActive()

void CML2SDWithParticle::setActive ( G4bool  act )

inline

Definition at line 65 of file ML2SDWithParticle.hh.

{bActive=act;}

Here is the call graph for this function:

Member Data Documentation

accTargetZPosition

G4double CML2SDWithParticle::accTargetZPosition

private

Definition at line 80 of file ML2SDWithParticle.hh.

bActive

G4bool CML2SDWithParticle::bActive

private

Definition at line 78 of file ML2SDWithParticle.hh.

bSavePhaseSpace

G4bool CML2SDWithParticle::bSavePhaseSpace

private

Definition at line 78 of file ML2SDWithParticle.hh.

bStopAtPhaseSpace

G4bool CML2SDWithParticle::bStopAtPhaseSpace

private

Definition at line 78 of file ML2SDWithParticle.hh.

fullOutFileData

G4String CML2SDWithParticle::fullOutFileData

private

Definition at line 75 of file ML2SDWithParticle.hh.

halfSize

G4ThreeVector CML2SDWithParticle::halfSize

private

Definition at line 71 of file ML2SDWithParticle.hh.

idType

G4int CML2SDWithParticle::idType

private

Definition at line 77 of file ML2SDWithParticle.hh.

max_N_particles_in_PhSp_File

G4int CML2SDWithParticle::max_N_particles_in_PhSp_File

private

Definition at line 77 of file ML2SDWithParticle.hh.

nMaxParticlesInRamPhaseSpace

G4int CML2SDWithParticle::nMaxParticlesInRamPhaseSpace

private

Definition at line 77 of file ML2SDWithParticle.hh.

nParticle

G4int CML2SDWithParticle::nParticle

private

Definition at line 76 of file ML2SDWithParticle.hh.

nTotalParticles

G4int CML2SDWithParticle::nTotalParticles

private

Definition at line 76 of file ML2SDWithParticle.hh.

particles

Sparticle* CML2SDWithParticle::particles

private

Definition at line 74 of file ML2SDWithParticle.hh.

pos

G4ThreeVector CML2SDWithParticle::pos

private

Definition at line 72 of file ML2SDWithParticle.hh.

primaryParticleData

SPrimaryParticle* CML2SDWithParticle::primaryParticleData

private

Definition at line 73 of file ML2SDWithParticle.hh.

---

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

• ML2SDWithParticle.hh
• ML2SDWithParticle.cc