XrayTelAnalysis Class Reference

#include <XrayTelAnalysis.hh>

Collaboration diagram for XrayTelAnalysis:

Public Member Functions
	~XrayTelAnalysis ()
void	book (G4bool isMaster)
void	finish (G4bool isMaster)
void	analyseStepping (const G4Track &track, G4bool entering)
void	Update (G4double energy, G4int threadID)

Static Public Member Functions
static XrayTelAnalysis *	getInstance ()

Private Member Functions
	XrayTelAnalysis ()

Private Attributes
G4double	eKin
G4double	x
G4double	y
G4double	z
G4double	dirX
G4double	dirY
G4double	dirZ
G4String	asciiFileName
G4String	histFileName
std::ofstream *	asciiFile
std::map< G4int, G4int > *	nEnteringTracks
std::map< G4int, G4double > *	totEnteringEnergy

Static Private Attributes
static XrayTelAnalysis *	instance = 0

Detailed Description

Definition at line 58 of file XrayTelAnalysis.hh.

Constructor & Destructor Documentation

~XrayTelAnalysis()

XrayTelAnalysis::~XrayTelAnalysis

(

)

Definition at line 75 of file XrayTelAnalysis.cc.

{
  if (asciiFile)
    delete asciiFile;
  if (nEnteringTracks)
    delete nEnteringTracks;
  if (totEnteringEnergy)
    delete totEnteringEnergy;
}

XrayTelAnalysis()

XrayTelAnalysis::XrayTelAnalysis ( )

private

Definition at line 62 of file XrayTelAnalysis.cc.

                                 : 
  asciiFile(0),nEnteringTracks(0), totEnteringEnergy(0)
{
  histFileName = "xraytel";


  asciiFileName="xraytel.out";
  asciiFile = new std::ofstream(asciiFileName);

  if(asciiFile->is_open()) 
    (*asciiFile) << "Energy (keV)  x (mm)    y (mm)    z (mm)" << G4endl << G4endl;  
}

Here is the caller graph for this function:

Member Function Documentation

analyseStepping()

void XrayTelAnalysis::analyseStepping	(	const G4Track &	track,
		G4bool	entering
	)

Definition at line 233 of file XrayTelAnalysis.cc.

{
  G4AutoLock l(&dataManipulationMutex);
  eKin = track.GetKineticEnergy()/keV;
  G4ThreeVector pos = track.GetPosition()/mm;
  y = pos.y();
  z = pos.z();
  G4ThreeVector dir = track.GetMomentumDirection();
  dirX = dir.x();
  dirY = dir.y();
  dirZ = dir.z();

  // Fill histograms
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillH1(1,eKin);
  man->FillH2(1,y,z);
  
  // Fill histograms and ntuple, tracks entering the detector
  if (entering) {
    // Fill and plot histograms
    man->FillH1(2,eKin);
    man->FillH2(2,y,z);

    man->FillNtupleDColumn(0,eKin);
    man->FillNtupleDColumn(1,x);
    man->FillNtupleDColumn(2,y);
    man->FillNtupleDColumn(3,z);
    man->FillNtupleDColumn(4,dirX);
    man->FillNtupleDColumn(5,dirY);
    man->FillNtupleDColumn(6,dirZ);
    man->AddNtupleRow();  
  }


  // Write to file
  if (entering) {
    if(asciiFile->is_open()) {
      (*asciiFile) << std::setiosflags(std::ios::fixed)
           << std::setprecision(3)
           << std::setiosflags(std::ios::right)
           << std::setw(10);
      (*asciiFile) << eKin;
      (*asciiFile) << std::setiosflags(std::ios::fixed)
           << std::setprecision(3)
           << std::setiosflags(std::ios::right)
           << std::setw(10);
      (*asciiFile) << x;
      (*asciiFile) << std::setiosflags(std::ios::fixed)
           << std::setprecision(3)
           << std::setiosflags(std::ios::right)
           << std::setw(10);
      (*asciiFile) << y;
      (*asciiFile) << std::setiosflags(std::ios::fixed)
           << std::setprecision(3)
           << std::setiosflags(std::ios::right)
           << std::setw(10);
      (*asciiFile) << z
           << G4endl;
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

book()

void XrayTelAnalysis::book

(

G4bool

isMaster

)

Definition at line 95 of file XrayTelAnalysis.cc.

{
  G4AutoLock l(&dataManipulationMutex);

  //reset counters: do be done only once, by the master
  if (isMaster)
    {
      if (nEnteringTracks)    
    {
      delete nEnteringTracks;
      nEnteringTracks = 0;
    }
      nEnteringTracks = new std::map<G4int,G4int>;
  
      if (totEnteringEnergy)
    {
      delete totEnteringEnergy;
      totEnteringEnergy = 0;
    }
      totEnteringEnergy = new std::map<G4int,G4double>;
    }

  // Get/create analysis manager
  G4AnalysisManager* man = G4AnalysisManager::Instance();

  // Open an output file: it is done in master and threads. The 
  // printout is done only by the master, for tidyness
  if (isMaster)
    G4cout << "Opening output file " << histFileName << " ... ";
  man->OpenFile(histFileName);
  man->SetFirstHistoId(1);
  if (isMaster)
    G4cout << " done" << G4endl;

  // Book 1D histograms
  man->CreateH1("h1","Energy, all /keV",  100,0.,100.);
  man->CreateH1("h2","Energy, entering detector /keV", 500,0.,500.);

  // Book 2D histograms (notice: the numbering is independent)
  man->CreateH2("d1","y-z, all /mm", 100,-500.,500.,100,-500.,500.); 
  man->CreateH2("d2","y-z, entering detector /mm", 200,-50.,50.,200,-50.,50.);
  
  // Book ntuples
  man->CreateNtuple("tree", "Track ntuple");
  man->CreateNtupleDColumn("energy");
  man->CreateNtupleDColumn("x");
  man->CreateNtupleDColumn("y");
  man->CreateNtupleDColumn("z");
  man->CreateNtupleDColumn("dirx");
  man->CreateNtupleDColumn("diry");
  man->CreateNtupleDColumn("dirz");
  man->FinishNtuple();
}

Here is the call graph for this function:

Here is the caller graph for this function:

finish()

void XrayTelAnalysis::finish

(

G4bool

isMaster

)

Definition at line 150 of file XrayTelAnalysis.cc.

{
  G4AutoLock l(&dataManipulationMutex);
  // Save histograms
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->Write();
  man->CloseFile();
  // Complete clean-up
  delete G4AnalysisManager::Instance();

  if (!isMaster)
    return;

  //only master performs these operations
  asciiFile->close();
 
  //Sequential run: output results!
  if (nEnteringTracks->count(-2))
    {
      G4cout << "End of Run summary (sequential)" << G4endl << G4endl;
      G4cout << "Total Entering Detector : " << nEnteringTracks->find(-2)->second  
         << G4endl;
      G4cout << "Total Entering Detector Energy : " 
         << (totEnteringEnergy->find(-2)->second)/MeV  
         << " MeV"
         << G4endl;
      return;
    }
  

  //MT run: sum results
 

  //MT build, but sequential run  
  if (nEnteringTracks->count(-1))
    {
      G4cout << "End of Run summary (sequential with MT build)" << G4endl << G4endl;
      G4cout << "Total Entering Detector : " << nEnteringTracks->find(-1)->second  
         << G4endl;
      G4cout << "Total Entering Detector Energy : " 
         << (totEnteringEnergy->find(-1)->second)/MeV  
         << " MeV"
         << G4endl;
      G4cout << "##########################################" << G4endl;
      return;
    }

  G4bool loopAgain = true;
  G4int totEntries = 0;
  G4double totEnergy = 0.;

  G4cout << "##########################################" << G4endl;
  for (size_t i=0; loopAgain; i++)
    {
      //ok, this thread was found
      if (nEnteringTracks->count(i))
    {
      G4cout << "End of Run summary (thread= " << i << ")" << G4endl;
      G4int part = nEnteringTracks->find(i)->second;
      G4cout << "Total Entering Detector : " << part << G4endl;
      G4double ene = totEnteringEnergy->find(i)->second;
      G4cout << "Total Entering Detector Energy : " 
         << ene/MeV  
         << " MeV" << G4endl;
      totEntries += part;
      totEnergy += ene;
      G4cout << "##########################################" << G4endl;
    }
      else
    loopAgain = false;
    }
  //Report global numbers
  if (totEntries)
    {
      G4cout << "End of Run summary (MT) global" << G4endl << G4endl;
      G4cout << "Total Entering Detector : " << totEntries << G4endl;
      G4cout << "Total Entering Detector Energy : " 
         << totEnergy/MeV  
         << " MeV" << G4endl;
      G4cout << "##########################################" << G4endl;
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

getInstance()

XrayTelAnalysis * XrayTelAnalysis::getInstance ( )

static

Definition at line 86 of file XrayTelAnalysis.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

Update()

void XrayTelAnalysis::Update	(	G4double	energy,
		G4int	threadID
	)

Definition at line 295 of file XrayTelAnalysis.cc.

{
  G4AutoLock l(&dataManipulationMutex);
  //It already exists: increase the counter
  if (nEnteringTracks->count(threadID))
    {
      (nEnteringTracks->find(threadID)->second)++;     
    }
  else //enter a new one
    {
      G4int tracks = 1;
      nEnteringTracks->insert(std::make_pair(threadID,tracks));
    }

  //It already exists: increase the counter
  if (totEnteringEnergy->count(threadID))
    (totEnteringEnergy->find(threadID)->second) += energy;
  else //enter a new one    
    totEnteringEnergy->insert(std::make_pair(threadID,energy));
    
  return;
}

Here is the caller graph for this function:

Member Data Documentation

asciiFile

std::ofstream* XrayTelAnalysis::asciiFile

private

Definition at line 92 of file XrayTelAnalysis.hh.

asciiFileName

G4String XrayTelAnalysis::asciiFileName

private

Definition at line 89 of file XrayTelAnalysis.hh.

dirX

G4double XrayTelAnalysis::dirX

private

Definition at line 85 of file XrayTelAnalysis.hh.

dirY

G4double XrayTelAnalysis::dirY

private

Definition at line 86 of file XrayTelAnalysis.hh.

dirZ

G4double XrayTelAnalysis::dirZ

private

Definition at line 87 of file XrayTelAnalysis.hh.

eKin

G4double XrayTelAnalysis::eKin

private

Definition at line 81 of file XrayTelAnalysis.hh.

histFileName

G4String XrayTelAnalysis::histFileName

private

Definition at line 90 of file XrayTelAnalysis.hh.

instance

XrayTelAnalysis * XrayTelAnalysis::instance = 0

staticprivate

Definition at line 78 of file XrayTelAnalysis.hh.

nEnteringTracks

std::map<G4int,G4int>* XrayTelAnalysis::nEnteringTracks

private

Definition at line 95 of file XrayTelAnalysis.hh.

totEnteringEnergy

std::map<G4int,G4double>* XrayTelAnalysis::totEnteringEnergy

private

Definition at line 96 of file XrayTelAnalysis.hh.

x

G4double XrayTelAnalysis::x

private

Definition at line 82 of file XrayTelAnalysis.hh.

y

G4double XrayTelAnalysis::y

private

Definition at line 83 of file XrayTelAnalysis.hh.

z

G4double XrayTelAnalysis::z

private

Definition at line 84 of file XrayTelAnalysis.hh.

---

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

• XrayTelAnalysis.hh
• XrayTelAnalysis.cc