PTGuiRollG3 Class Reference

#include <PTGuiRollG3.h>

List of all members.

Public Member Functions
	PTGuiRollG3 ()
virtual	~PTGuiRollG3 ()
TClonesArray *	GetMaterials () const
TClonesArray *	GetMedia () const
Private Member Functions
	PTGuiRollG3 (const PTGuiRollG3 &that)
PTGuiRollG3 &	operator= (const PTGuiRollG3 &that)
void	RollMaterials ()
void	RollMedia ()
Double_t	Cut (Int_t imed, Int_t icut)
Private Attributes
TClonesArray *	fMaterials
TClonesArray *	fMedia
TGeant3 *	fGeant3

---

Detailed Description

Definition at line 28 of file PTGuiRollG3.h.

---

Constructor & Destructor Documentation

PTGuiRollG3::PTGuiRollG3

(

)

Definition at line 36 of file PTGuiRollG3.cxx.

                         : fMaterials(0), fMedia(0), fGeant3(0) {
  // Normal Constructor

    fMaterials = new TClonesArray("PTGuiMaterial",100);
    fMedia     = new TClonesArray("PTGuiMedium",100);
    
    if ( gMC ) fGeant3 = (TGeant3*) gMC;
    if ( !fGeant3 ) {
      cerr << "PTGuiRollG3 ctor. Null TGeant3 global ptr. Abort." << endl;
      abort();
    }

  // Convert Geant3 objects to PTGui objects
    RollMaterials();
    RollMedia();
    
}

PTGuiRollG3::~PTGuiRollG3

(

)

[virtual]

Definition at line 55 of file PTGuiRollG3.cxx.

References fMaterials, and fMedia.

                          {
  // Destructor

  // Delete TClonesArrays and contained objects
  if ( fMaterials ) { fMaterials->Delete(); delete fMaterials; }
  if ( fMedia ) { fMedia -> Delete(); delete fMedia; }
  
  // fGeant3 is not owned

}

PTGuiRollG3::PTGuiRollG3

(

const PTGuiRollG3 &

that

)

[private]

---

Member Function Documentation

Double_t PTGuiRollG3::Cut	(	Int_t	imed,
		Int_t	icut
	)			[private]

Definition at line 169 of file PTGuiRollG3.cxx.

References fGeant3, TGeant3::Gclink(), Gclink_t::jtmed, TGeant3::Lq(), and TGeant3::Q().

Referenced by RollMedia().

                                                {
  // Return cut value for cut icut, for medium idmed
  // "Cut" energy thresholds and physics process flags can be customized 
  // by medium.  Check here for such customization or otherwise apply defaults.
  // As described in the G3 GSTPAR documentation, valid icuts are in the
  // range 1-22,26-29,and 31-33

  Gclink_t* gclink = fGeant3->Gclink(); // pointer to GCLINK common block
  Int_t* zlq = fGeant3->Lq(); // pointer to Zebra bank lq
  Float_t* zq = fGeant3->Q();  // pointer to Zebra bank q
  
  Int_t jtm = zlq[gclink->jtmed-imed];
  Int_t jtmn = zlq[jtm];
  if ( jtmn > 0 ) {
    // the cuts for this medium have been specialized
    jtm = jtmn;
  }
  else {
    jtm = gclink -> jtmed;
  }

  return Double_t(zq[jtm+icut]);
  
}

TClonesArray* PTGuiRollG3::GetMaterials

(

)

const [inline]

Definition at line 36 of file PTGuiRollG3.h.

References fMaterials.

Referenced by PTGuiMainFrame::UpdateMaterialList().

{return fMaterials;}    

TClonesArray* PTGuiRollG3::GetMedia

(

)

const [inline]

Definition at line 37 of file PTGuiRollG3.h.

References fMedia.

Referenced by PTGuiMainFrame::AddComboMedia(), and PTGuiMainFrame::ProcessMessage().

{return fMedia;}    

PTGuiRollG3& PTGuiRollG3::operator=

(

const PTGuiRollG3 &

that

)

[private]

void PTGuiRollG3::RollMaterials

(

)

[private]

Definition at line 67 of file PTGuiRollG3.cxx.

References MuELoss::a, fGeant3, TGeant3::Gclink(), TGeant3::Gcnum(), TGeant3::Gfmate(), TGeant3::Iq(), Gclink_t::jmate, TGeant3::Lq(), nentries, Gcnum_t::nmate, and Lit::Print().

                                {
  // Read materials from G3 banks using TGeant3 interface, and store
  // in fMaterials array of PTGuiMaterial objects.
  //  

  char namate[21];
  Float_t a, z, dens, radl, absl;
  Int_t npar = 0;
  Float_t* par = 0;
        
  // Loop over materials
  Int_t nmat = fGeant3->Gcnum()->nmate;
    
  Int_t nentries = 0; // number of array entries
  // Start at 1 instead of 0 because G3 fortran indexes from 1
  for ( Int_t imat = 1; imat < nmat; imat++ ) {
    // Get material parameters
    fGeant3->Gfmate(imat,namate,a,z,dens,radl,absl,par,npar);
    if ( a < 0. ) continue; // this is used to indicate non-existent material
    // How odd that this is necessary, but namate returned by Gfmate is
    // incomplete
    Int_t* zlq = fGeant3->Lq();
    Int_t* ziq = fGeant3->Iq();
    Int_t jma = Int_t(zlq[fGeant3->Gclink()->jmate-imat]);
    strncpy(namate,(char*)&ziq[jma+1],20);
    namate[20] = '\0';

    // Fuss with the material name to remove trailing spaces
    std::string matname = namate;
    
    int pos(matname.size());
    for ( ; pos && matname[pos-1]==' '; --pos );
    matname.erase(pos,matname.size()-pos);

    // Create new material object
    PTGuiMaterial* material=new PTGuiMaterial(matname,imat,a,z,dens,radl,absl);
    (*fMaterials)[nentries++] = material;
    material -> Print();
  }

}

void PTGuiRollG3::RollMedia

(

)

[private]

Definition at line 110 of file PTGuiRollG3.cxx.

References Cut(), fGeant3, TGeant3::Gcnum(), TGeant3::Gftmed(), nentries, Gcnum_t::ntmed, and Lit::Print().

                            {
  // Read media from G3 banks using TGeant3 interface, and store
  // in fMedia array of PTGuiMedium objects.
  // Cuts (e.g. CUTGAM) and physics process flags (e.g. LOSS), can
  // be specified on a medium by medium basis (although the default
  // is to apply one set to all media), so are store here as
  // part of the PTGuiMedium object.   

  char named[21];
  Int_t imat, isvol, ifield;
  Float_t fieldm, tmaxfd, stemax, deemax, epsil, stmin;
  Int_t nwbuf = 0;
  Float_t* ubuf = 0;
    
  // Loop over media
  Int_t nmed = fGeant3->Gcnum()->ntmed;
    
  Int_t nentries = 0; // number of array entries

  // Start at 1 instead of 0 because G3 fortran indexes from 1
  for ( Int_t imed = 1; imed < nmed; imed++ ) {

    fGeant3->Gftmed (imed,named,imat,isvol,ifield,fieldm,tmaxfd,
                      stemax,deemax,epsil,stmin,ubuf,&nwbuf);

    named[20] = '\0';

    if ( imat <= 0 ) continue; // this is used to indicate non-existent medium

    // Fuss with the medium name to remove trailing spaces
    std::string medname = named;
    int pos(medname.size());
    for ( ; pos && medname[pos-1]==' '; --pos );
    medname.erase(pos,medname.size()-pos);

    // Create new medium object
    PTGuiMedium* medium = new PTGuiMedium(medname,imed,imat,isvol,ifield,
                                   fieldm,tmaxfd,stemax,deemax,epsil,stmin);
    (*fMedia)[nentries++] = medium;

    // Store cut/physic process flags associated with medium in
    // medium parameters array

    // Valid medium associated cuts are described in the GSTPAR documentation
    for ( int icut = 0; icut < 10; icut++ ) 
      medium -> SetCut(icut,Cut(imed,icut+1));
     
    // Valid medium associated physics process flags are described in the
    // GSTPar documention
    for ( int ip = 0; ip < 12; ip++ )
      medium -> SetProcess(ip,(Int_t)Cut(imed,ip+11));

    medium -> Print();
    
  }
  
}

---

Member Data Documentation

TGeant3* PTGuiRollG3::fGeant3 [private]

Definition at line 56 of file PTGuiRollG3.h.

Referenced by Cut(), RollMaterials(), and RollMedia().

TClonesArray* PTGuiRollG3::fMaterials [private]

Definition at line 51 of file PTGuiRollG3.h.

Referenced by GetMaterials(), and ~PTGuiRollG3().

TClonesArray* PTGuiRollG3::fMedia [private]

Definition at line 52 of file PTGuiRollG3.h.

Referenced by GetMedia(), and ~PTGuiRollG3().

---

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

• PTGuiRollG3.h
• PTGuiRollG3.cxx

---