GammaRayTelAnalysis Class Reference

#include <GammaRayTelAnalysis.hh>

Collaboration diagram for GammaRayTelAnalysis:

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

Private Member Functions
	GammaRayTelAnalysis ()
void	Plot ()

Private Attributes
const GammaRayTelDetectorConstruction *	GammaRayTelDetector
G4String	histo2DMode
G4String	histoFileName
GammaRayTelAnalysisMessenger *	analysisMessenger

Static Private Attributes
static GammaRayTelAnalysis *	instance = 0

Detailed Description

Definition at line 71 of file GammaRayTelAnalysis.hh.

Constructor & Destructor Documentation

~GammaRayTelAnalysis()

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:

GammaRayTelAnalysis()

GammaRayTelAnalysis::GammaRayTelAnalysis ( )

private

Definition at line 68 of file GammaRayTelAnalysis.cc.

  :GammaRayTelDetector(0),histo2DMode("strip")
{
  GammaRayTelDetector =
    static_cast<const GammaRayTelDetectorConstruction*>
        (G4RunManager::GetRunManager()->GetUserDetectorConstruction());

  // Define the messenger and the analysis system
  analysisMessenger = new GammaRayTelAnalysisMessenger(this);  
  histoFileName = "gammaraytel";
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Function Documentation

BeginOfRun()

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:

EndOfEvent()

void GammaRayTelAnalysis::EndOfEvent

(

G4int

flag

)

Definition at line 239 of file GammaRayTelAnalysis.cc.

{;}

Here is the caller graph for this function:

EndOfRun()

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:

Finish()

void GammaRayTelAnalysis::Finish

(

)

Definition at line 91 of file GammaRayTelAnalysis.cc.

{
  delete analysisMessenger;
  analysisMessenger = 0;
}             

Here is the caller graph for this function:

GetHisto2DMode()

G4String GammaRayTelAnalysis::GetHisto2DMode ( )

inline

Definition at line 88 of file GammaRayTelAnalysis.hh.

{return histo2DMode;};

Here is the call graph for this function:

Here is the caller graph for this function:

getInstance()

GammaRayTelAnalysis * GammaRayTelAnalysis::getInstance ( )

static

Definition at line 97 of file GammaRayTelAnalysis.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

Init()

void GammaRayTelAnalysis::Init

(

)

Definition at line 88 of file GammaRayTelAnalysis.cc.

{;}                       

InsertEnergy()

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:

InsertHits()

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:

InsertPositionXZ()

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:

InsertPositionYZ()

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:

Plot()

void GammaRayTelAnalysis::Plot ( )

private

Here is the caller graph for this function:

SetHisto2DMode()

void GammaRayTelAnalysis::SetHisto2DMode ( G4String  str )

inline

Definition at line 87 of file GammaRayTelAnalysis.hh.

{histo2DMode = str;};

Here is the caller graph for this function:

setNtuple()

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:

Member Data Documentation

analysisMessenger

GammaRayTelAnalysisMessenger* GammaRayTelAnalysis::analysisMessenger

private

Definition at line 115 of file GammaRayTelAnalysis.hh.

GammaRayTelDetector

const GammaRayTelDetectorConstruction* GammaRayTelAnalysis::GammaRayTelDetector

private

Definition at line 110 of file GammaRayTelAnalysis.hh.

histo2DMode

G4String GammaRayTelAnalysis::histo2DMode

private

Definition at line 112 of file GammaRayTelAnalysis.hh.

histoFileName

G4String GammaRayTelAnalysis::histoFileName

private

Definition at line 113 of file GammaRayTelAnalysis.hh.

instance

GammaRayTelAnalysis * GammaRayTelAnalysis::instance = 0

staticprivate

Definition at line 108 of file GammaRayTelAnalysis.hh.

---

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

• GammaRayTelAnalysis.hh
• GammaRayTelAnalysis.cc