G4ParticleHPField Class Reference

#include <G4ParticleHPField.hh>

Collaboration diagram for G4ParticleHPField:

Public Member Functions
	G4ParticleHPField ()
	~G4ParticleHPField ()
void	InitY (G4int i, G4int n)
void	SetData (G4int i, G4double x, G4int j, G4double y)
void	SetEnergy (G4int i, G4double e)
void	SetX (G4int i, G4double e)
void	SetY (G4int i, G4int j, G4double x)
G4double	GetEnergy (G4int i)
G4double	GetX (G4int i)
G4double	GetY (G4int i, G4int j)
G4ParticleHPFieldPoint &	GetPoint (G4int i)
G4double	GetY (G4double e, G4int j)
G4int	GetFieldLength ()
void	Dump ()

Private Member Functions
void	Check (G4int i)

Private Attributes
G4ParticleHPFieldPoint *	theData
G4int	nEntries
G4int	nPoints

Detailed Description

Definition at line 35 of file G4ParticleHPField.hh.

Constructor & Destructor Documentation

G4ParticleHPField()

G4ParticleHPField::G4ParticleHPField

(

)

Definition at line 37 of file G4ParticleHPField.cc.

  {
    theData = new G4ParticleHPFieldPoint[100]; 
    nPoints=100;
    nEntries=0;
    theData->SetData(0,0,0);
  }

Here is the call graph for this function:

~G4ParticleHPField()

G4ParticleHPField::~G4ParticleHPField

(

)

Definition at line 45 of file G4ParticleHPField.cc.

{ delete [] theData;}

Member Function Documentation

Check()

void G4ParticleHPField::Check ( G4int  i )

private

Definition at line 104 of file G4ParticleHPField.cc.

  {
    if(i>nEntries) throw G4HadronicException(__FILE__, __LINE__, "Skipped some index numbers in G4ParticleHPField");
    if(i==nPoints)
    {
      nPoints += 50;
      G4ParticleHPFieldPoint * buff = new G4ParticleHPFieldPoint[nPoints];
//      G4cout << "copying 1"<<G4endl;
      for (G4int j=0; j<nEntries; j++) 
      {
        buff[j] = theData[j];
      }
//      G4cout << "copying 2"<<G4endl;
      delete [] theData;
      theData = buff;
    }
    if(i==nEntries) nEntries=i+1;
  }

Here is the caller graph for this function:

Dump()

void G4ParticleHPField::Dump

(

)

Definition at line 90 of file G4ParticleHPField.cc.

  {
    G4cout << nEntries<<G4endl;
    for(G4int i=0; i<nEntries; i++)
    {
      G4cout << theData[i].GetX()<<" ";
      for(G4int j=0; j<theData[i].GetDepth(); j++)
      {
        G4cout << theData[i].GetY(j)<<" ";
      }
      G4cout << G4endl;
    }
  }

Here is the call graph for this function:

Here is the caller graph for this function:

GetEnergy()

G4double G4ParticleHPField::GetEnergy ( G4int  i )

inline

Definition at line 68 of file G4ParticleHPField.hh.

{ return theData[i].GetX(); }

Here is the call graph for this function:

GetFieldLength()

G4int G4ParticleHPField::GetFieldLength ( )

inline

Definition at line 75 of file G4ParticleHPField.hh.

{return nEntries;}

Here is the call graph for this function:

GetPoint()

G4ParticleHPFieldPoint& G4ParticleHPField::GetPoint ( G4int  i )

inline

Definition at line 71 of file G4ParticleHPField.hh.

{ return theData[i]; }

Here is the call graph for this function:

GetX()

G4double G4ParticleHPField::GetX ( G4int  i )

inline

Definition at line 69 of file G4ParticleHPField.hh.

{ return theData[i].GetX(); }

Here is the call graph for this function:

Here is the caller graph for this function:

GetY() [1/2]

G4double G4ParticleHPField::GetY ( G4int  i, G4int  j  )

inline

Definition at line 70 of file G4ParticleHPField.hh.

{ return theData[i].GetY(j); }

Here is the call graph for this function:

Here is the caller graph for this function:

GetY() [2/2]

G4double G4ParticleHPField::GetY	(	G4double	e,
		G4int	j
	)

Definition at line 47 of file G4ParticleHPField.cc.

  {
    G4int low   = 0;
    G4int high  = 0;
    G4int i;
    for (i=1; i<nEntries/10; i++)
    {
      if(theData[10*i].GetX()>e) break;
    }
    if(i==(nEntries/10))
    {
      i=10*i;
      while (i<nEntries) // Loop checking, 11.05.2015, T. Koi
      {
        if(theData[i++].GetX()>e) break;
      } 
      if (i==nEntries)
      {
        low  = nEntries-1;
        high = nEntries-2;
      }else{
        low = i-1;
        high = i;
      }
    }else{
      for (G4int jj=0; jj<10; jj++)
      {
        if(theData[i].GetX()<e) break;
        i--;
      }
      low = i;
      high = i+1;
    }
    G4double x1, x2, y1, y2, x, y;
    x = e;
    x1 = theData[low] .GetX();
    x2 = theData[high].GetX();
    y1 = theData[low] .GetY(j);
    y2 = theData[high].GetY(j);
    y = x*(y2-y1)/(x2-x1);
    return y += y2-x2*(y2-y1)/(x2-x1);
  }

Here is the call graph for this function:

InitY()

void G4ParticleHPField::InitY ( G4int  i, G4int  n  )

inline

Definition at line 43 of file G4ParticleHPField.hh.

  {
    Check(i);
    theData[i].InitY(n);
  }

Here is the call graph for this function:

SetData()

void G4ParticleHPField::SetData ( G4int  i, G4double  x, G4int  j, G4double  y  )

inline

Definition at line 48 of file G4ParticleHPField.hh.

  { 
    Check(i);
    theData[i].SetData(x, j, y);
  }

Here is the call graph for this function:

SetEnergy()

void G4ParticleHPField::SetEnergy ( G4int  i, G4double  e  )

inline

Definition at line 53 of file G4ParticleHPField.hh.

  {
    Check(i);
    theData[i].SetX(e);
  }

Here is the call graph for this function:

SetX()

void G4ParticleHPField::SetX ( G4int  i, G4double  e  )

inline

Definition at line 58 of file G4ParticleHPField.hh.

  {
    Check(i);
    theData[i].SetX(e);
  }

Here is the call graph for this function:

SetY()

void G4ParticleHPField::SetY ( G4int  i, G4int  j, G4double  x  )

inline

Definition at line 63 of file G4ParticleHPField.hh.

  {
    Check(i);
    theData[i].SetY(j, x);
  }

Here is the call graph for this function:

Member Data Documentation

nEntries

G4int G4ParticleHPField::nEntries

private

Definition at line 84 of file G4ParticleHPField.hh.

nPoints

G4int G4ParticleHPField::nPoints

private

Definition at line 85 of file G4ParticleHPField.hh.

theData

G4ParticleHPFieldPoint* G4ParticleHPField::theData

private

Definition at line 83 of file G4ParticleHPField.hh.

---

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

• G4ParticleHPField.hh
• G4ParticleHPField.cc