CCalAnalysis Class Reference

Alternative here is a XML file. More...

#include <CCalAnalysis.hh>

Collaboration diagram for CCalAnalysis:

Public Member Functions
virtual	~CCalAnalysis ()
void	BeginOfRun (G4int n)
void	EndOfRun (G4int n)
void	EndOfEvent (G4int flag)
void	Init ()
void	Finish ()
int	maxbin ()
void	InsertEnergyHcal (float *)
void	InsertEnergyEcal (float *)
void	InsertEnergy (float v)
void	InsertLateralProfile (float *)
void	InsertTime (float *)
void	InsertTimeProfile (int, double, double)
void	setNtuple (float *hcalE, float *ecalE, float elab, float x, float y, float z, float edep, float edec, float edhc)

Static Public Member Functions
static CCalAnalysis *	getInstance ()

Private Member Functions
	CCalAnalysis ()

Private Attributes
G4int	fVerbosity
G4int	numberOfTimeSlices
G4int	energy
G4int	hcalE
G4int	ecalE
G4int	timeHist
G4int	lateralProfile
G4int	timeProfile

Static Private Attributes
static CCalAnalysis *	instance = 0

Detailed Description

Alternative here is a XML file.

Defines the format which is used for the output file default is a ROOT file. Comment the g4root and un-comment one of the others, to change the output format

Definition at line 44 of file CCalAnalysis.hh.

Constructor & Destructor Documentation

~CCalAnalysis()

CCalAnalysis::~CCalAnalysis ( )

virtual

Definition at line 140 of file CCalAnalysis.cc.

                            {  
  Finish();
  delete G4AnalysisManager::Instance();
}

Here is the call graph for this function:

CCalAnalysis()

CCalAnalysis::CCalAnalysis ( )

private

Definition at line 39 of file CCalAnalysis.cc.

                           :
  energy(0), hcalE(0), ecalE(0), timeHist(0), lateralProfile(0),
  timeProfile(0)
{

#ifdef debug
  fVerbosity = 1;
#else
  fVerbosity = 0;  
#endif

  numberOfTimeSlices = 200;

  // Create analysis manager
  G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
  analysisManager->SetVerboseLevel(1);
  analysisManager->SetFirstHistoId(1);
  analysisManager->SetFirstNtupleId(1);

  // Open an output file
  analysisManager->OpenFile("ccal");
  G4cout << "********************************************" << G4endl
     << "* o/p file ccal"  << G4endl
     << "********************************************" << G4endl 
     << G4endl;
  

  // Create a tuple :
  // Create ntuple
  analysisManager->CreateNtuple("ntuple1", "Event info");
  char tupleid[6];
  for (int i=0;i<28;i++)
    {
      sprintf(tupleid,"hcal%d",i);
      analysisManager->CreateNtupleFColumn(tupleid);
    }
  for (int i=0; i<49; i++) 
    {
      sprintf(tupleid,"ecal%d",i);
      analysisManager->CreateNtupleFColumn(tupleid);
    }
  analysisManager->CreateNtupleFColumn("ELAB");
  analysisManager->CreateNtupleFColumn("XPOS");
  analysisManager->CreateNtupleFColumn("YPOS");
  analysisManager->CreateNtupleFColumn("ZPOS");
  analysisManager->CreateNtupleFColumn("EDEP");
  analysisManager->CreateNtupleFColumn("EDEC");
  analysisManager->CreateNtupleFColumn("EHDC");
  analysisManager->FinishNtuple();

  //
  //Create histograms. Save the ID of the first histogram in each block
  //
  char id[5], ntupletag[50];
  //Energy deposit in Hcal layers 
  for (int i = 0; i<28; i++) {
    sprintf(id,"h%d",i+100);
    sprintf(ntupletag, "Energy Deposit in Hcal Layer%d   in GeV",i);
    G4int histoID = analysisManager->CreateH1(id,ntupletag, 100, 0., 1.0);
    if (!i)
      hcalE = histoID;
  }
  // Energy deposits in Ecal towers
  for (int i = 0; i<49; i++) 
    {
      sprintf(id, "h%d",i+200);
      sprintf(ntupletag, "Energy Deposit in Ecal Tower%d   in GeV",i);
      G4int histoID = analysisManager->CreateH1(id,ntupletag, 100, 0., 1.0);
      if (!i)
    ecalE = histoID;
    }
  // Total energy deposit
  energy = analysisManager->CreateH1("h4000", "Total energy deposited in GeV", 
                     100, 0., 100.0);

  // Time slices      
  for (int i=0; i<numberOfTimeSlices; i++){
    sprintf(id, "h%d",i+300);
    sprintf(ntupletag, "Time slice %d nsec energy profile   in GeV",i);
    G4int histoID = analysisManager->CreateH1(id,ntupletag,100, 0., 100.0);
    if (!i)
      timeHist = histoID;
  }

  // Profile of lateral energy deposit in Hcal
  for (int i = 0; i<70; i++) {
    sprintf(id, "h%d",i+500);
    sprintf(ntupletag, "Lateral energy profile at %d cm  in GeV",i);
    G4int histoID = analysisManager->CreateH1(id, ntupletag, 100, 0., 10.0);
    if (!i)
      lateralProfile = histoID;
  }

  // Time profile 
  timeProfile = analysisManager->CreateH1("h901", "Time Profile in Sensitive Detector", 
                        200, 0., 200.);
  analysisManager->CreateH1("h902", "Time Profile in Sensitive+Passive", 
                200, 0., 200.);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Function Documentation

BeginOfRun()

void CCalAnalysis::BeginOfRun

(

G4int

n

)

Definition at line 312 of file CCalAnalysis.cc.

{ 
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  for (int i=0; i<numberOfTimeSlices; i++){
    man->GetH1(timeHist+1)->reset();
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

EndOfEvent()

void CCalAnalysis::EndOfEvent

(

G4int

flag

)

Definition at line 333 of file CCalAnalysis.cc.

                                        {
  // The plotter is updated only if there is some
  // hits in the event
  if (!flag) return;
}

Here is the caller graph for this function:

EndOfRun()

void CCalAnalysis::EndOfRun

(

G4int

n

)

Definition at line 323 of file CCalAnalysis.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 CCalAnalysis::Finish

(

)

Definition at line 149 of file CCalAnalysis.cc.

{;}             

Here is the caller graph for this function:

getInstance()

CCalAnalysis * CCalAnalysis::getInstance ( )

static

Definition at line 152 of file CCalAnalysis.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

Init()

void CCalAnalysis::Init

(

)

Definition at line 146 of file CCalAnalysis.cc.

                        {
}                       

InsertEnergy()

void CCalAnalysis::InsertEnergy

(

float

v

)

Definition at line 225 of file CCalAnalysis.cc.

{
 
  G4double x = v;
  G4AnalysisManager::Instance()->FillH1(energy,x);
  if (fVerbosity)
    G4cout << "CCalAnalysis::InsertEnergy: Fill Total energy Hcal histo with " 
       << x << G4endl;
}

Here is the caller graph for this function:

InsertEnergyEcal()

void CCalAnalysis::InsertEnergyEcal

(

float *

v

)

Definition at line 181 of file CCalAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  
  G4double totalFilledEnergyEcal = 0.0;

  for (G4int i=0; i<49; i++) {    
    G4double x = v[i];
    man->FillH1(ecalE+i,x);
    if (fVerbosity)
      {
    G4cout << "Fill Ecal histo " << i << " with " << x << G4endl;
    totalFilledEnergyEcal += x;
      }
  }
  if (fVerbosity)
    G4cout << 
      "CCalAnalysis::InsertEnergyEcal: Total filled Energy Ecal histo " 
       << totalFilledEnergyEcal << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

InsertEnergyHcal()

void CCalAnalysis::InsertEnergyHcal

(

float *

v

)

Definition at line 159 of file CCalAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  G4double totalFilledEnergyHcal = 0.0;
  for (int i=0; i<28; i++) {
    G4double x = v[i];
    man->FillH1(hcalE+i,x);
    if (fVerbosity)
      {
    G4cout << "Fill Hcal histo " << i << " with " << x << G4endl;
    totalFilledEnergyHcal += x;
      }
  }  

  if (fVerbosity)
    G4cout << 
      "CCalAnalysis::InsertEnergyHcal: Total filled Energy Hcal histo " 
       << totalFilledEnergyHcal << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

InsertLateralProfile()

void CCalAnalysis::InsertLateralProfile

(

float *

v

)

Definition at line 204 of file CCalAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  G4double totalFilledProfileHcal = 0.0;

  for (G4int i=0; i<70; i++) {    
      G4double x = v[i];
      man->FillH1(lateralProfile+1,x);
      if (fVerbosity)
    {
      G4cout << "Fill Profile Hcal histo " << i << " with " << x << G4endl;
      totalFilledProfileHcal += x;
    }    
  }
  if (fVerbosity)
    G4cout << "CCalAnalysis::InsertLateralProfile: Total filled Profile Hcal"
       << " histo " << totalFilledProfileHcal << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

InsertTime()

void CCalAnalysis::InsertTime

(

float *

v

)

Definition at line 237 of file CCalAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  
  G4double totalFilledTimeProfile = 0.0;
  for (G4int j=0; j<numberOfTimeSlices; j++) 
    {
      G4double x = v[j];
      man->FillH1(timeHist+j,x);
      if (fVerbosity)
    {
      G4cout << "Fill Time slice histo " << j << " with " << x << G4endl;
      totalFilledTimeProfile += x;
    }    

      G4double t = j + 0.5;
      man->FillH1(timeProfile+1,t,x);
      if (fVerbosity)
    G4cout << "Fill Time profile histo 1 with " << t << " " << x << G4endl;
    }
  if (fVerbosity)
    G4cout << "CCalAnalysis::InsertTime: Total filled Time profile histo " 
       << totalFilledTimeProfile << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

InsertTimeProfile()

void CCalAnalysis::InsertTimeProfile	(	int	hit,
		double	time,
		double	edep
	)

Definition at line 264 of file CCalAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();
  man->FillH1(timeProfile,time,edep);

  if (fVerbosity)
    G4cout << "CCalAnalysis:: Fill Time Profile with Hit " << hit
       << " Edeposit " << edep << " Gev at " << time << " ns" << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

maxbin()

int CCalAnalysis::maxbin ( )

inline

Definition at line 56 of file CCalAnalysis.hh.

{return numberOfTimeSlices;}

Here is the call graph for this function:

Here is the caller graph for this function:

setNtuple()

void CCalAnalysis::setNtuple	(	float *	hcalE,
		float *	ecalE,
		float	elab,
		float	x,
		float	y,
		float	z,
		float	edep,
		float	edec,
		float	edhc
	)

Definition at line 276 of file CCalAnalysis.cc.

{
  G4AnalysisManager* man = G4AnalysisManager::Instance();  
  G4int counter=0;
  for (int i=0; i<28; i++) 
    {
      man->FillNtupleFColumn(counter,HCalE[i]);
      counter++;
    }
  for (int i=0; i<49; i++) 
    {
      man->FillNtupleFColumn(counter,ECalE[i]);
      counter++;
    }
  man->FillNtupleFColumn(counter,elab);
  man->FillNtupleFColumn(counter+1,x);
  man->FillNtupleFColumn(counter+2,y);
  man->FillNtupleFColumn(counter+3,z);
  man->FillNtupleFColumn(counter+4,edep);
  man->FillNtupleFColumn(counter+5,edec);
  man->FillNtupleFColumn(counter+6,edhc);

  man->AddNtupleRow();  

  if (fVerbosity)
    G4cout << "CCalAnalysis:: Fill Ntuple " << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

ecalE

G4int CCalAnalysis::ecalE

private

Definition at line 87 of file CCalAnalysis.hh.

energy

G4int CCalAnalysis::energy

private

Definition at line 85 of file CCalAnalysis.hh.

fVerbosity

G4int CCalAnalysis::fVerbosity

private

Definition at line 76 of file CCalAnalysis.hh.

hcalE

G4int CCalAnalysis::hcalE

private

Definition at line 86 of file CCalAnalysis.hh.

instance

CCalAnalysis * CCalAnalysis::instance = 0

staticprivate

Definition at line 74 of file CCalAnalysis.hh.

lateralProfile

G4int CCalAnalysis::lateralProfile

private

Definition at line 89 of file CCalAnalysis.hh.

numberOfTimeSlices

G4int CCalAnalysis::numberOfTimeSlices

private

Definition at line 78 of file CCalAnalysis.hh.

timeHist

G4int CCalAnalysis::timeHist

private

Definition at line 88 of file CCalAnalysis.hh.

timeProfile

G4int CCalAnalysis::timeProfile

private

Definition at line 91 of file CCalAnalysis.hh.

---

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

• CCalAnalysis.hh
• CCalAnalysis.cc