exGPSHistoManager Class Reference

#include <exGPSHistoManager.hh>

Public Member Functions
	exGPSHistoManager ()
	~exGPSHistoManager ()
void	book ()
void	save ()
void	Fill (G4int PDGid, G4double e, G4double x, G4double y, G4double z, G4double t, G4double p, G4double w)
void	PrintStatistic ()
void	SetEMin (G4double emin)
void	SetEMax (G4double emax)
void	SetPosMin (G4double posmin)
void	SetPosMax (G4double posmax)
void	SetFileName (G4String name)

Detailed Description

Definition at line 50 of file exGPSHistoManager.hh.

Constructor & Destructor Documentation

exGPSHistoManager::exGPSHistoManager

(

)

Definition at line 42 of file exGPSHistoManager.cc.

                                    :
fMinpos(-10.),fMaxpos(10),fMineng(0.),fMaxeng(1000.),
fEnerHisto(0),fPosiXY(0),fPosiZX(0),fPosiYZ(0),fAnglCTP(0),fAnglTP(0)
{
  fFileName[0] = "expGPS";
  fFactoryOn = false;
  fMessenger = new exGPSHistoMessenger(this);
}

exGPSHistoManager::~exGPSHistoManager

(

)

Definition at line 53 of file exGPSHistoManager.cc.

{ }

Member Function Documentation

void exGPSHistoManager::book

(

)

Definition at line 58 of file exGPSHistoManager.cc.

References G4VAnalysisManager::CreateH1(), G4VAnalysisManager::CreateH2(), G4VAnalysisManager::CreateNtuple(), G4VAnalysisManager::CreateNtupleDColumn(), G4VAnalysisManager::CreateNtupleIColumn(), G4VAnalysisManager::FinishNtuple(), G4cout, G4endl, G4VAnalysisManager::GetFileType(), G4RootAnalysisManager::GetH1(), G4RootAnalysisManager::GetH2(), G4VAnalysisManager::OpenFile(), G4VAnalysisManager::SetFirstHistoId(), G4VAnalysisManager::SetFirstNtupleId(), G4VAnalysisManager::SetHistoDirectoryName(), G4VAnalysisManager::SetNtupleDirectoryName(), and G4VAnalysisManager::SetVerboseLevel().

Referenced by exGPSRunAction::BeginOfRunAction().

{
  // Create or get analysis manager
  // The choice of analysis technology is done via selection of a namespace
  // in exGPSHistoManager.hh
  G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
  analysisManager->SetVerboseLevel(2);
  G4String extension = analysisManager->GetFileType();
  fFileName[1] = fFileName[0] + "." + extension;
      
  // Create directories 
  analysisManager->SetHistoDirectoryName("histo");
  analysisManager->SetNtupleDirectoryName("ntuple");
    
  // Open an output file
  //
  G4bool fileOpen = analysisManager->OpenFile(fFileName[0]);
  if (!fileOpen) {
    G4cout << "\n---> exGPSHistoManager::book(): cannot open " << fFileName[1]
           << G4endl;
    return;
  }
  
  // create selected histograms
  //
  analysisManager->SetFirstHistoId(1);
  analysisManager->SetFirstNtupleId(1);

  G4int id = analysisManager->CreateH1("1","Ekin primary",100,fMineng,fMaxeng);
  fEnerHisto= analysisManager->GetH1(id);

  id = analysisManager->CreateH2("2","Position XY",100,fMinpos,fMaxpos,
                                                   100,fMinpos,fMaxpos);
  fPosiXY= analysisManager->GetH2(id);

  id = analysisManager->CreateH2("3","Position YZ",100,fMinpos,fMaxpos,
                                                     100,fMinpos,fMaxpos);
  fPosiYZ= analysisManager->GetH2(id);

  id = analysisManager->CreateH2("4","Position ZX",100,fMinpos,fMaxpos,
                                                       100,fMinpos,fMaxpos);
  fPosiZX= analysisManager->GetH2(id);

  id =analysisManager->CreateH2("5","Source phi-std::cos(theta) distribution",
          360,0,360,100, -1, 1);
  fAnglCTP =analysisManager->GetH2(id);

  id =analysisManager->CreateH2("6","Source phi-theta distribution",
          360,0,360,180,0,180);
  fAnglTP = analysisManager->GetH2(id);

  // Create 1st ntuple (id = 1)
  //    
  analysisManager->CreateNtuple("MyTuple", "Primary Particle Tuple");
  id = analysisManager->CreateNtupleIColumn("particleID");
  id = analysisManager->CreateNtupleDColumn("Ekin");
  id = analysisManager->CreateNtupleDColumn("posX");
  id = analysisManager->CreateNtupleDColumn("posY");
  id = analysisManager->CreateNtupleDColumn("posZ");
  id = analysisManager->CreateNtupleDColumn("dirTheta");
  id = analysisManager->CreateNtupleDColumn("dirPhi");
  id = analysisManager->CreateNtupleDColumn("weight");
  analysisManager->FinishNtuple();
  
  fFactoryOn = true;       
  G4cout << "\n----> Histogram Tree is opened in " << fFileName[1] << G4endl;
}

void exGPSHistoManager::Fill	(	G4int	PDGid,
		G4double	e,
		G4double	x,
		G4double	y,
		G4double	z,
		G4double	t,
		G4double	p,
		G4double	w
	)

Definition at line 141 of file exGPSHistoManager.cc.

References G4VAnalysisManager::AddNtupleRow(), python.hepunit::cm, python.hepunit::deg, G4VAnalysisManager::FillNtupleDColumn(), G4VAnalysisManager::FillNtupleIColumn(), and python.hepunit::MeV.

Referenced by exGPSEventAction::EndOfEventAction().

{
   fEnerHisto->fill(e, w);
   fEnerHisto->fill(e/MeV,w);
   fPosiXY->fill(x/cm,y/cm,w);
   fPosiZX->fill(z/cm,x/cm,w);
   fPosiYZ->fill(y/cm,z/cm,w);
   fAnglCTP->fill(p/deg,std::cos(t),w);
   fAnglTP->fill(p/deg,t/deg,w);

   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
   analysisManager->FillNtupleIColumn(1,0,PDGid);
   analysisManager->FillNtupleDColumn(1,1,e);
   analysisManager->FillNtupleDColumn(1,2,x);
   analysisManager->FillNtupleDColumn(1,3,y);
   analysisManager->FillNtupleDColumn(1,4,z);
   analysisManager->FillNtupleDColumn(1,5,t);
   analysisManager->FillNtupleDColumn(1,6,p);
   analysisManager->FillNtupleDColumn(1,7,w);
   analysisManager->AddNtupleRow(1);
}

void exGPSHistoManager::PrintStatistic

(

)

Definition at line 167 of file exGPSHistoManager.cc.

References G4cout, and G4endl.

{
  if(fFactoryOn) {
    G4cout << "\n ----> print histograms statistic \n" << G4endl;
    
    /*G4cout
       << " EAbs : mean = " << G4BestUnit(fHistPt[1]->mean(), "Energy") 
               << " rms = " << G4BestUnit(fHistPt[1]->rms(),  "Energy") 
               << G4endl;
    G4cout                
       << " EGap : mean = " << G4BestUnit(fHistPt[2]->mean(), "Energy") 
               << " rms = " << G4BestUnit(fHistPt[2]->rms(),  "Energy") 
               << G4endl;
    G4cout 
       << " LAbs : mean = " << G4BestUnit(fHistPt[3]->mean(), "Length") 
               << " rms = " << G4BestUnit(fHistPt[3]->rms(),  "Length") 
               << G4endl;
    G4cout 
       << " LGap : mean = " << G4BestUnit(fHistPt[4]->mean(), "Length") 
               << " rms = " << G4BestUnit(fHistPt[4]->rms(),  "Length") 
               << G4endl;
     */
  }
}

void exGPSHistoManager::save

(

)

Definition at line 128 of file exGPSHistoManager.cc.

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

Referenced by exGPSRunAction::EndOfRunAction().

{
  if (fFactoryOn) {
    G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();    
    analysisManager->Write();
    analysisManager->CloseFile();  
    G4cout << "\n----> Histogram Tree is saved in " << fFileName[1] << G4endl;
      
    delete G4AnalysisManager::Instance();
    fFactoryOn = false;
  }                    
}

void exGPSHistoManager::SetEMax ( G4double  emax )

inline

Definition at line 77 of file exGPSHistoManager.hh.

Referenced by exGPSHistoMessenger::SetNewValue().

{fMaxeng = emax;}

void exGPSHistoManager::SetEMin ( G4double  emin )

inline

Definition at line 76 of file exGPSHistoManager.hh.

Referenced by exGPSHistoMessenger::SetNewValue().

{fMineng = emin;}

void exGPSHistoManager::SetFileName ( G4String  name )

inline

Definition at line 80 of file exGPSHistoManager.hh.

Referenced by exGPSHistoMessenger::SetNewValue().

{fFileName[0] = name;}

void exGPSHistoManager::SetPosMax ( G4double  posmax )

inline

Definition at line 79 of file exGPSHistoManager.hh.

Referenced by exGPSHistoMessenger::SetNewValue().

{fMaxpos = posmax;}

void exGPSHistoManager::SetPosMin ( G4double  posmin )

inline

Definition at line 78 of file exGPSHistoManager.hh.

Referenced by exGPSHistoMessenger::SetNewValue().

{fMinpos = posmin;}

---

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

• exGPSHistoManager.hh
• exGPSHistoManager.cc