CML2SDWithParticle Class Reference

#include <ML2SDWithParticle.hh>

Inheritance 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 ()
const 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

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::CML2SDWithParticle

(

)

Definition at line 48 of file ML2SDWithParticle.cc.

References idSD_KillerPlane.

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

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.

References test::a.

: 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

(

void

)

Definition at line 83 of file ML2SDWithParticle.cc.

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

Member Function Documentation

CML2SDWithParticle* CML2SDWithParticle::getCML2SensitiveDetectorParticle ( )

inline

Definition at line 63 of file ML2SDWithParticle.hh.

Referenced by CML2PhaseSpaces::getCML2SensitiveDetectorParticle().

{return this;}

Sparticle CML2SDWithParticle::getParticle ( int  i )

inline

Definition at line 64 of file ML2SDWithParticle.hh.

{return particles[i];}

G4int CML2SDWithParticle::getTotalNumberOfParticles ( )

inline

Definition at line 62 of file ML2SDWithParticle.hh.

Referenced by CML2PhaseSpaces::getCML2SensDetNParticle().

{return nTotalParticles;}

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

virtual

Implements G4VSensitiveDetector.

Definition at line 121 of file ML2SDWithParticle.cc.

References Sparticle::dir, fStopAndKill, G4Track::GetDefinition(), CML2AcceleratorConstruction::GetInstance(), G4Track::GetKineticEnergy(), G4StepPoint::GetMomentumDirection(), G4ParticleDefinition::GetPDGEncoding(), G4StepPoint::GetPosition(), G4Step::GetPreStepPoint(), pyG4VTouchable::GetRotation, G4Step::GetTrack(), idSD_KillerPlane, Sparticle::kinEnergy, Sparticle::nPrimaryPart, SPrimaryParticle::nPrimaryParticle, Sparticle::partPDGE, SPrimaryParticle::partPDGE, Sparticle::pos, Sparticle::primaryParticlePDGE, G4Track::SetTrackStatus(), Sparticle::volumeId, and Sparticle::volumeName.

{
    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;
}

void CML2SDWithParticle::save

(

)

Definition at line 190 of file ML2SDWithParticle.cc.

References CML2AcceleratorConstruction::GetInstance(), and pyG4VTouchable::GetRotation.

Referenced by CML2PhaseSpaces::save().

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

void CML2SDWithParticle::setActive ( G4bool  act )

inline

Definition at line 65 of file ML2SDWithParticle.hh.

{bActive=act;}

---

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

• ML2SDWithParticle.hh
• ML2SDWithParticle.cc