GammaRayTelAnalysis Class Reference

#include <GammaRayTelAnalysis.hh>

Public Member Functions
virtual	~GammaRayTelAnalysis ()
void	BeginOfRun ()
void	EndOfRun ()
void	EndOfEvent (G4int flag)
void	Init ()
void	Finish ()
void	SetHisto2DMode (G4String str)
G4String	GetHisto2DMode ()
void	InsertPositionXZ (double x, double z)
void	InsertPositionYZ (double y, double z)
void	InsertEnergy (double en)
void	InsertHits (int nplane)
void	setNtuple (float E, float p, float x, float y, float z)

Static Public Member Functions
static GammaRayTelAnalysis *	getInstance ()

Detailed Description

Definition at line 71 of file GammaRayTelAnalysis.hh.

Constructor & Destructor Documentation

GammaRayTelAnalysis::~GammaRayTelAnalysis ( )

virtual

Definition at line 81 of file GammaRayTelAnalysis.cc.

                                          {
  Finish();
  // Complete clean-up
  delete G4AnalysisManager::Instance();
}

Here is the call graph for this function:

Member Function Documentation

void GammaRayTelAnalysis::BeginOfRun

(

)

Definition at line 150 of file GammaRayTelAnalysis.cc.

{ 
  G4AnalysisManager* man = G4AnalysisManager::Instance();

  // Open an output file

  G4cout << "Opening output file " << histoFileName << " ... ";
  man->OpenFile(histoFileName);
  man->SetFirstHistoId(1);
  G4cout << " done" << G4endl;
  

  int Nplane = GammaRayTelDetector->GetNbOfTKRLayers();
  int Nstrip = GammaRayTelDetector->GetNbOfTKRStrips();
  int Ntile = GammaRayTelDetector->GetNbOfTKRTiles();
  double sizexy = GammaRayTelDetector->GetTKRSizeXY();
  double sizez = GammaRayTelDetector->GetTKRSizeZ();
  int N = Nstrip*Ntile;      

  // Book1D histograms 
  //------------------

  // 1D histogram that store the energy deposition of the
  // particle in the last (number 0) TKR X-plane
  man->CreateH1("1","Edep in the last X plane (keV)", 100, 50, 200);

  // 1D histogram that store the hits distribution along the TKR X-planes
  man->CreateH1("2","Hits dist in TKR X planes",Nplane, 0, Nplane-1);

  // Book 2D histograms 
  //-------------------

  // 2D histogram that store the position (mm) of the hits (XZ projection)

  if (histo2DMode == "strip")
    {
      man->CreateH2("1","Tracker Hits XZ (strip,plane)", 
            N, 0, N-1, 
            2*Nplane, 0, Nplane-1); 
    }
  else
    {
      man->CreateH2("1","Tracker Hits XZ (x,z) in mm", 
            int(sizexy/5), -sizexy/2, sizexy/2, 
            int(sizez/5), -sizez/2, sizez/2);
    }  

  // 2D histogram that store the position (mm) of the hits (YZ projection)  
  if (histo2DMode == "strip")
    {
      man->CreateH2("2","Tracker Hits YZ (strip,plane)", 
            N, 0, N-1, 
            2*Nplane, 0, Nplane-1); 
    }
  else
    {
      man->CreateH2("2","Tracker Hits YZ (x,z) in mm", 
            int(sizexy/5), -sizexy/2, sizexy/2, 
            int(sizez/5), -sizez/2, sizez/2);
    }  
  
    
  // Book Ntuples (energy / plane/ x / y / z)
  //------------------------------------------  
  man->CreateNtuple("1","Track ntuple");
  man->CreateNtupleDColumn("energy");
  man->CreateNtupleDColumn("plane"); // can I use Int values? 
  man->CreateNtupleDColumn("x"); // can I use Int values?
  man->CreateNtupleDColumn("y"); // can I use Int values?
  man->CreateNtupleDColumn("z"); // can I use Int values?
  man->FinishNtuple();

}

Here is the call graph for this function:

Here is the caller graph for this function:

void GammaRayTelAnalysis::EndOfEvent

(

G4int

flag

)

Definition at line 239 of file GammaRayTelAnalysis.cc.

{;}

Here is the caller graph for this function:

void GammaRayTelAnalysis::EndOfRun

(

)

Definition at line 229 of file GammaRayTelAnalysis.cc.

{
  //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:

void GammaRayTelAnalysis::Finish

(

)

Definition at line 91 of file GammaRayTelAnalysis.cc.

{
  delete analysisMessenger;
  analysisMessenger = 0;
}             

Here is the caller graph for this function:

G4String GammaRayTelAnalysis::GetHisto2DMode ( )

inline

Definition at line 88 of file GammaRayTelAnalysis.hh.

{return histo2DMode;};

Here is the caller graph for this function:

GammaRayTelAnalysis * GammaRayTelAnalysis::getInstance ( )

static

Definition at line 97 of file GammaRayTelAnalysis.cc.

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

Here is the caller graph for this function:

void GammaRayTelAnalysis::Init

(

)

Definition at line 88 of file GammaRayTelAnalysis.cc.

{;}                       

void GammaRayTelAnalysis::InsertEnergy

(

double

en

)

Definition at line 118 of file GammaRayTelAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillH1(1,en);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void GammaRayTelAnalysis::InsertHits

(

int

nplane

)

Definition at line 125 of file GammaRayTelAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillH1(2,nplane);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void GammaRayTelAnalysis::InsertPositionXZ	(	double	x,
		double	z
	)

Definition at line 104 of file GammaRayTelAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillH2(1,x,z);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void GammaRayTelAnalysis::InsertPositionYZ	(	double	y,
		double	z
	)

Definition at line 111 of file GammaRayTelAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillH2(2,y,z);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void GammaRayTelAnalysis::SetHisto2DMode ( G4String  str )

inline

Definition at line 87 of file GammaRayTelAnalysis.hh.

{histo2DMode = str;};

Here is the caller graph for this function:

void GammaRayTelAnalysis::setNtuple	(	float	E,
		float	p,
		float	x,
		float	y,
		float	z
	)

Definition at line 131 of file GammaRayTelAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillNtupleDColumn(0,E);
  man->FillNtupleDColumn(1,p);
  man->FillNtupleDColumn(2,x);
  man->FillNtupleDColumn(3,y);
  man->FillNtupleDColumn(4,z);
  man->AddNtupleRow();
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p02/examples/advanced/gammaray_telescope/include/GammaRayTelAnalysis.hh
• source/geant4.10.03.p02/examples/advanced/gammaray_telescope/src/GammaRayTelAnalysis.cc