XrayFluoAnalysisManager Class Reference

#include <XrayFluoAnalysisManager.hh>

Collaboration diagram for XrayFluoAnalysisManager:

Public Member Functions
virtual	~XrayFluoAnalysisManager ()
void	book ()
void	finish ()
void	analyseStepping (const G4Step *aStep)
void	analyseEnergyDep (G4double eDep)
void	analysePrimaryGenerator (G4double energy)
void	PhaseSpaceOn ()
void	PhaseSpaceOff ()
void	SetOutputFileName (G4String)
const std::pair< G4double, G4String >	GetEmittedParticleEnergyAndType ()
void	LoadGunData (G4String, G4bool)
void	SetPhysicFlag (G4bool)

Static Public Member Functions
static XrayFluoAnalysisManager *	getInstance ()

Private Member Functions
	XrayFluoAnalysisManager ()

Private Attributes
G4String	outputFileName
G4bool	phaseSpaceFlag
G4bool	physicFlag
std::vector< G4double > *	gunParticleEnergies
std::vector< G4String > *	gunParticleTypes
XrayFluoAnalysisMessenger *	analisysMessenger
G4bool	dataLoaded
G4int	fParticleEnergyAndTypeIndex

Static Private Attributes
static XrayFluoAnalysisManager *	instance = 0

Detailed Description

Definition at line 55 of file XrayFluoAnalysisManager.hh.

Constructor & Destructor Documentation

~XrayFluoAnalysisManager()

XrayFluoAnalysisManager::~XrayFluoAnalysisManager ( )

virtual

Definition at line 88 of file XrayFluoAnalysisManager.cc.

{
  if ( gunParticleEnergies ) delete gunParticleEnergies;
  gunParticleEnergies = 0;
  if ( gunParticleTypes ) delete gunParticleTypes;
  gunParticleTypes = 0;

  delete instance;
  instance = 0;
  
  delete G4AnalysisManager::Instance();  
}

XrayFluoAnalysisManager()

XrayFluoAnalysisManager::XrayFluoAnalysisManager ( )

private

Definition at line 70 of file XrayFluoAnalysisManager.cc.

  :outputFileName("xrayfluo"), phaseSpaceFlag(false), physicFlag (false), 
   gunParticleEnergies(0), gunParticleTypes(0)
{
  dataLoaded = false;
  fParticleEnergyAndTypeIndex = 0;

  //creating the messenger
  analisysMessenger = new XrayFluoAnalysisMessenger(this);
  
  //Instantiate the analysis manager
  G4AnalysisManager::Instance();

  G4cout << "XrayFluoAnalysisManager created" << G4endl;
}

Here is the caller graph for this function:

Member Function Documentation

analyseEnergyDep()

void XrayFluoAnalysisManager::analyseEnergyDep

(

G4double

eDep

)

Definition at line 430 of file XrayFluoAnalysisManager.cc.

{
  G4AutoLock l(&dataManipulationMutex);
  // Filling of Energy Deposition, called by XrayFluoEventAction  
  G4AnalysisManager* man = G4AnalysisManager::Instance();

  if (!phaseSpaceFlag)
    man->FillH1(1,energyDep/keV);
  
}

Here is the call graph for this function:

Here is the caller graph for this function:

analysePrimaryGenerator()

void XrayFluoAnalysisManager::analysePrimaryGenerator

(

G4double

energy

)

Definition at line 443 of file XrayFluoAnalysisManager.cc.

{
  G4AutoLock l(&dataManipulationMutex);
  // Filling of energy spectrum histogram of the primary generator
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  if (!phaseSpaceFlag)
    man->FillH1(7,energy/keV);
  
}

Here is the call graph for this function:

Here is the caller graph for this function:

analyseStepping()

void XrayFluoAnalysisManager::analyseStepping

(

const G4Step *

aStep

)

Definition at line 243 of file XrayFluoAnalysisManager.cc.

{
  G4AutoLock l(&dataManipulationMutex);
  G4AnalysisManager* man = G4AnalysisManager::Instance();

  if (phaseSpaceFlag){       
    G4ParticleDefinition* particleType= 0;
    G4String parentProcess="";
    G4ThreeVector momentum(0.,0.,0.);
    G4double particleEnergy=0;
    G4String sampleMaterial="";
    G4double particleDepth=0;
    G4int isBornInTheSample=0;
    XrayFluoDetectorConstruction* detector = XrayFluoDetectorConstruction::GetInstance();    

    // Select volume from which the step starts
    if(aStep->GetPreStepPoint()->GetPhysicalVolume()->GetName()=="Sample"){
      G4ThreeVector creationPos = aStep->GetTrack()->GetVertexPosition();
      // qua serve ancora una selezione: deve essere interno al sample
      //codice "a mano" per controllare il volume di orgine della traccia
     
      G4VPhysicalVolume* creationPosVolume = detector->GetGeometryNavigator()->LocateGlobalPointAndSetup(creationPos);

      // if physicflag is on, I record all the data of all the particle born in the sample. 
      // If it is off (but phase space production is on) I record data 
      // only for particles creted inside the sample and exiting it.  
      if (physicFlag ^ (!physicFlag && 
            (aStep->GetPostStepPoint()->GetStepStatus() == fGeomBoundary) 
            )) 
    {   
      // extracting information needed
      particleType = aStep->GetTrack()->GetDynamicParticle()->GetDefinition();
      momentum = aStep->GetTrack()->GetDynamicParticle()->GetMomentum();
      particleEnergy = aStep->GetPreStepPoint()->GetKineticEnergy();
      if (creationPosVolume->GetName() == "Sample" ) {
        isBornInTheSample = 1;
        particleDepth = creationPosVolume->GetLogicalVolume()->GetSolid()
          ->DistanceToOut(creationPos, G4ThreeVector(0,0,-1));
      }
      else {
        particleDepth = -1;
      }
      // metodo per ottenere la profondita' "a mano"
      //    particleDepth = std::abs(creationPos.z() - detector->GetSampleIlluminatedFaceCoord()); 

      G4int parent=0;
      if(aStep->GetTrack()->GetCreatorProcess())
        {
          parentProcess = aStep->GetTrack()->GetCreatorProcess()->GetProcessName();
          parent = 5;
          // hack for HBK
          if (parentProcess == "") parent = 0;
          if (parentProcess == "ioni") parent = 1;
          if (parentProcess == "LowEnPhotoElec") parent = 2;
          if (parentProcess == "Transportation") parent = 3;// should never happen
          if (parentProcess == "initStep") parent = 4;// should never happen
        }   
      else {
        parentProcess = "Not Known -- (primary generator + Rayleigh)";
        parent = 5;
      }
      G4int sampleMat=0;
      if(aStep->GetTrack()){
        sampleMaterial = aStep->GetTrack()->GetMaterial()->GetName();
        if (sampleMaterial == ("Dolorite" || "Anorthosite" || "Mars1" || "IceBasalt" || "HPGe")) sampleMat=1;
      }
    
      G4int part = -1 ;
      if (particleType == G4Gamma::Definition()) part =1; 
      if (particleType == G4Electron::Definition()) part = 0;
      if (particleType == G4Proton::Definition()) part = 2;
      
      //Fill ntuple  
      man->FillNtupleIColumn(1,0, part);
      man->FillNtupleDColumn(1,1,particleEnergy);
      man->FillNtupleDColumn(1,2,momentum.theta());
      man->FillNtupleDColumn(1,3,momentum.phi());
      man->FillNtupleIColumn(1,4,parent); 
      man->FillNtupleIColumn(1,5,sampleMat);
      man->FillNtupleIColumn(1,6,isBornInTheSample);
      man->FillNtupleDColumn(1,7,particleDepth);
      man->AddNtupleRow(1);
      
    }
    }
  }
    
  // Normal behaviour, without creation of phase space
  else 
    {
   
      // Filling the histograms with data, passing thru stepping.

      // Select volume from wich the step starts
      if(aStep->GetPreStepPoint()->GetPhysicalVolume()->GetName()=="Sample"){
    // Select volume from wich the step ends
    if(aStep->GetTrack()->GetNextVolume()->GetName() == "World" ) { 
      // Select the particle type
      if ((aStep->GetTrack()->GetDynamicParticle()
           ->GetDefinition()) == G4Gamma::Definition() )
    
        {         
          G4double gammaLeavingSample = 
        aStep->GetPreStepPoint()->GetKineticEnergy();
          man->FillH1(4,gammaLeavingSample/keV);

        }   
      else if ((aStep->GetTrack()->GetDynamicParticle()
            ->GetDefinition()) == G4Electron::Definition() ) 
        {
          G4double eleLeavingSample = 
        aStep->GetPreStepPoint()->GetKineticEnergy();
          man->FillH1(5,eleLeavingSample/keV);
        }
      else if ((aStep->GetTrack()->GetDynamicParticle()
            ->GetDefinition()) == G4Proton::Definition() )
        {
          G4double 
        protonsLeavSam = aStep->GetPreStepPoint()->GetKineticEnergy();
          man->FillH1(9,protonsLeavSam/keV);
        }
    }   
      }
    
      if((aStep->GetTrack()->GetDynamicParticle()
      ->GetDefinition()) == G4Gamma::Definition() )
        {
      if(aStep->GetTrack()->GetCurrentStepNumber() == 1)
        {         
          if(aStep->GetTrack()->GetParentID() != 0)
        {   
          if(aStep->GetTrack()->GetVolume()->GetName() == "Sample")
            {
              G4double gammaBornInSample = 
            aStep->GetPreStepPoint()->GetKineticEnergy();
              man->FillH1(2,gammaBornInSample/keV);
            }
        }
        }
    }    
      else if ((aStep->GetTrack()->GetDynamicParticle()
        ->GetDefinition() ) == G4Electron::Definition() )
    {
      if(aStep->GetTrack()->GetCurrentStepNumber() == 1)
        {         
          if(aStep->GetTrack()->GetParentID() != 0)
        {   
          if(aStep->GetTrack()->GetVolume()->GetName() == "Sample")
            {
              G4double eleBornInSample = 
            aStep->GetPreStepPoint()->GetKineticEnergy();
              man->FillH1(3,eleBornInSample/keV);
            }
        }
        }
    }
      
  
      if(aStep->GetTrack()->GetNextVolume())
    {

      
      if(aStep->GetTrack()->GetNextVolume()->GetName() == "HPGeDetector")
        
        { 
          if ((aStep->GetTrack()->GetDynamicParticle()
           ->GetDefinition()) == G4Gamma::Definition() ) 
        {

          G4double gammaAtTheDetPre = 
            aStep->GetPreStepPoint()->GetKineticEnergy();
          man->FillH1(6,gammaAtTheDetPre/keV);
        }
          else if ((aStep->GetTrack()->GetDynamicParticle()
            ->GetDefinition() ) == G4Proton::Definition() ) 
        {
          G4double protonsAtTheDetPre = 
            aStep->GetPreStepPoint()->GetKineticEnergy();
          man->FillH1(8,protonsAtTheDetPre);
        }
        }
    }
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

book()

void XrayFluoAnalysisManager::book

(

)

Definition at line 113 of file XrayFluoAnalysisManager.cc.

{
  G4AutoLock l(&dataManipulationMutex);
  // Get analysis manager
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  // Open an output file
  man->OpenFile(outputFileName);
  man->SetVerboseLevel(1);
  man->SetFirstHistoId(1);
  man->SetFirstNtupleId(1);

  G4cout << "Open output file: " << outputFileName << G4endl;

  if (phaseSpaceFlag) 
    {   
      // Book output Tuple (ID = 1)
      man->CreateNtuple("101","OutputNTuple");
      man->CreateNtupleIColumn("particle");      
      man->CreateNtupleDColumn("energies");
      man->CreateNtupleDColumn("momentumTheta");
      man->CreateNtupleDColumn("momentumPhi");
      man->CreateNtupleIColumn("processes");
      man->CreateNtupleIColumn("material");
      man->CreateNtupleIColumn("origin");
      man->CreateNtupleDColumn("depth");
      man->FinishNtuple();
      G4cout << "Created ntuple for phase space" << G4endl;
    }
  else {
    // Book histograms
    man->CreateH1("h1","Energy Deposit", 500,0.,10.); //20eV def.
    man->CreateH1("h2","Gamma born in the sample", 100,0.,10.);
    man->CreateH1("h3","Electrons  born in the sample", 100,0.,10.);
    man->CreateH1("h4","Gammas leaving the sample", 300,0.,10.);
    man->CreateH1("h5","Electrons leaving the sample ",200000 ,0.,10.0); // .05 eV def.
    man->CreateH1("h6","Gammas reaching the detector", 100,0.,10.);
    man->CreateH1("h7","Spectrum of the incident particles", 100,0.,10.);
    man->CreateH1("h8","Protons reaching the detector", 100,0.,10.);
    man->CreateH1("h9","Protons leaving the sample", 100,0.,10.);
    G4cout << "Created histos" << G4endl;
  }
} 

Here is the call graph for this function:

Here is the caller graph for this function:

finish()

void XrayFluoAnalysisManager::finish

(

)

Definition at line 223 of file XrayFluoAnalysisManager.cc.

{ 
  G4AutoLock l(&dataManipulationMutex);
  G4cout << "Going to save histograms" << G4endl;
  // Save histograms
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->Write();
  man->CloseFile();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetEmittedParticleEnergyAndType()

const std::pair< G4double, G4String > XrayFluoAnalysisManager::GetEmittedParticleEnergyAndType

(

)

Definition at line 203 of file XrayFluoAnalysisManager.cc.

{
  G4AutoLock l(&dataManipulationMutex);
  std::pair<G4double,G4String> result;
  
  if(fParticleEnergyAndTypeIndex < (G4int) gunParticleEnergies->size())
    {
      G4double energy = gunParticleEnergies->at(fParticleEnergyAndTypeIndex);
      G4String name = gunParticleTypes->at(fParticleEnergyAndTypeIndex);
      result.first = energy;
      result.second = name;
    }

  fParticleEnergyAndTypeIndex++;
  return result;
}

Here is the caller graph for this function:

getInstance()

XrayFluoAnalysisManager * XrayFluoAnalysisManager::getInstance ( )

static

Definition at line 103 of file XrayFluoAnalysisManager.cc.

{
  G4AutoLock l(&instanceMutex);
  if (instance == 0) {instance = new XrayFluoAnalysisManager;}
  return instance;
}

Here is the call graph for this function:

Here is the caller graph for this function:

LoadGunData()

void XrayFluoAnalysisManager::LoadGunData	(	G4String	fileName,
		G4bool	raileighFlag
	)

Definition at line 158 of file XrayFluoAnalysisManager.cc.

{
  G4AutoLock l(&dataManipulationMutex);

  //Do not allow more than one thread to trigger the file reading
  if (dataLoaded)
    return;

  // Get analysis reader manager
  G4AnalysisReader* analysisReader = G4AnalysisReader::Instance();
  analysisReader->SetVerboseLevel(1);

  //This is for testing purposes
  G4int ntupleId = analysisReader->GetNtuple("101",fileName);
  G4cout << "Got ntupleId: " << ntupleId << G4endl;

  gunParticleEnergies = new std::vector<G4double>;
  gunParticleTypes = new std::vector<G4String>;

  G4int particleID, processID;
  G4double energy;
  analysisReader->SetNtupleIColumn("processes", processID);
  analysisReader->SetNtupleDColumn("energies", energy);
  analysisReader->SetNtupleIColumn("particles", particleID);

  while (analysisReader->GetNtupleRow() ) 
    {
      if (raileighFlag ^ (!raileighFlag && (processID == 1 || 
                        processID == 2)) )  
    {
      gunParticleEnergies->push_back(energy*MeV);
      if (particleID == 1)
        gunParticleTypes->push_back("gamma");
      else if (particleID == 0)
        gunParticleTypes->push_back("e-");
    }
    }

  G4cout << "Maximum mumber of events: "<< gunParticleEnergies->size() <<G4endl;  

  dataLoaded= true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

PhaseSpaceOff()

void XrayFluoAnalysisManager::PhaseSpaceOff ( )

inline

Definition at line 80 of file XrayFluoAnalysisManager.hh.

{phaseSpaceFlag = false;}

Here is the call graph for this function:

PhaseSpaceOn()

void XrayFluoAnalysisManager::PhaseSpaceOn ( )

inline

Definition at line 78 of file XrayFluoAnalysisManager.hh.

{phaseSpaceFlag = true;}

Here is the caller graph for this function:

SetOutputFileName()

void XrayFluoAnalysisManager::SetOutputFileName

(

G4String

newName

)

Definition at line 455 of file XrayFluoAnalysisManager.cc.

{
  
  outputFileName = newName;
}

Here is the caller graph for this function:

SetPhysicFlag()

void XrayFluoAnalysisManager::SetPhysicFlag

(

G4bool

val

)

Definition at line 236 of file XrayFluoAnalysisManager.cc.

{
  physicFlag = val;
}

Here is the caller graph for this function:

Member Data Documentation

analisysMessenger

XrayFluoAnalysisMessenger* XrayFluoAnalysisManager::analisysMessenger

private

Definition at line 108 of file XrayFluoAnalysisManager.hh.

dataLoaded

G4bool XrayFluoAnalysisManager::dataLoaded

private

Definition at line 110 of file XrayFluoAnalysisManager.hh.

fParticleEnergyAndTypeIndex

G4int XrayFluoAnalysisManager::fParticleEnergyAndTypeIndex

private

Definition at line 112 of file XrayFluoAnalysisManager.hh.

gunParticleEnergies

std::vector<G4double>* XrayFluoAnalysisManager::gunParticleEnergies

private

Definition at line 101 of file XrayFluoAnalysisManager.hh.

gunParticleTypes

std::vector<G4String>* XrayFluoAnalysisManager::gunParticleTypes

private

Definition at line 102 of file XrayFluoAnalysisManager.hh.

instance

XrayFluoAnalysisManager * XrayFluoAnalysisManager::instance = 0

staticprivate

Definition at line 105 of file XrayFluoAnalysisManager.hh.

outputFileName

G4String XrayFluoAnalysisManager::outputFileName

private

Definition at line 95 of file XrayFluoAnalysisManager.hh.

phaseSpaceFlag

G4bool XrayFluoAnalysisManager::phaseSpaceFlag

private

Definition at line 97 of file XrayFluoAnalysisManager.hh.

physicFlag

G4bool XrayFluoAnalysisManager::physicFlag

private

Definition at line 99 of file XrayFluoAnalysisManager.hh.

---

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

• XrayFluoAnalysisManager.hh
• XrayFluoAnalysisManager.cc