XrayTelAnalysis Class Reference

#include <XrayTelAnalysis.hh>

Public Member Functions
	~XrayTelAnalysis ()
void	book ()
void	finish ()
void	analyseStepping (const G4Track &track, G4bool entering)

Static Public Member Functions
static XrayTelAnalysis *	getInstance ()

Detailed Description

Definition at line 58 of file XrayTelAnalysis.hh.

Constructor & Destructor Documentation

XrayTelAnalysis::~XrayTelAnalysis

(

)

Definition at line 66 of file XrayTelAnalysis.cc.

{
  if (asciiFile)
    delete asciiFile;
}

Member Function Documentation

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

Definition at line 124 of file XrayTelAnalysis.cc.

References G4VAnalysisManager::AddNtupleRow(), G4VAnalysisManager::FillH1(), G4VAnalysisManager::FillH2(), G4VAnalysisManager::FillNtupleDColumn(), G4endl, G4Track::GetKineticEnergy(), G4Track::GetMomentumDirection(), G4Track::GetPosition(), python.hepunit::keV, python.hepunit::mm, right, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by XrayTelSteppingAction::UserSteppingAction().

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

void XrayTelAnalysis::book

(

)

Definition at line 80 of file XrayTelAnalysis.cc.

References G4VAnalysisManager::CreateH1(), G4VAnalysisManager::CreateH2(), G4VAnalysisManager::CreateNtuple(), G4VAnalysisManager::CreateNtupleDColumn(), G4VAnalysisManager::FinishNtuple(), G4cout, G4endl, G4VAnalysisManager::OpenFile(), and G4VAnalysisManager::SetFirstHistoId().

Referenced by XrayTelRunAction::BeginOfRunAction().

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

  // Open an output file
  G4cout << "Opening output file " << histFileName << " ... ";
  man->OpenFile(histFileName);
  man->SetFirstHistoId(1);
  G4cout << " done" << G4endl;

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

  // Book 2D histograms (notice: the numbering is independent)
  man->CreateH2("1","y-z, all /mm", 100,-500.,500.,100,-500.,500.); 
  man->CreateH2("2","y-z, entering detector /mm", 200,-50.,50.,200,-50.,50.);
  
  // Book ntuples
  man->CreateNtuple("10", "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();
}

void XrayTelAnalysis::finish

(

)

Definition at line 112 of file XrayTelAnalysis.cc.

References G4VAnalysisManager::CloseFile(), and G4VAnalysisManager::Write().

Referenced by XrayTelRunAction::EndOfRunAction().

{
  asciiFile->close();

  // Save histograms
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->Write();
  man->CloseFile();
  // Complete clean-up
  delete G4AnalysisManager::Instance();
}

XrayTelAnalysis * XrayTelAnalysis::getInstance ( )

static

Definition at line 73 of file XrayTelAnalysis.cc.

Referenced by XrayTelRunAction::BeginOfRunAction(), XrayTelRunAction::EndOfRunAction(), and XrayTelSteppingAction::UserSteppingAction().

{
  if (instance == 0) instance = new XrayTelAnalysis;
  return instance;
}

---

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

• XrayTelAnalysis.hh
• XrayTelAnalysis.cc