#include <CCalG4Hcal.hh>

Inheritance diagram for CCalG4Hcal:

Public Member Functions
	CCalG4Hcal (const G4String &name)

virtual	~CCalG4Hcal ()

Public Member Functions inherited from CCalHcal
	CCalHcal (const G4String &name)

virtual	~CCalHcal ()

G4String	getGenMat () const

double	getDy_2Cal () const

double	getDx_2Cal () const

double	getXposCal () const

G4String	getBoxMat () const

int	getNBox () const

double	getDy_2Box () const

double	getDx_2Box () const

double	getWallThickBox () const

double	getXposBox (unsigned int i) const

int	getNLayerScnt () const

int	getTypeScnt (unsigned int i) const

int	getMotherScnt (unsigned int i) const

double	getXposScnt (unsigned int i) const

int	getNLayerAbs () const

int	getTypeAbs (unsigned int i) const

int	getMotherAbs (unsigned int i) const

double	getXposAbs (unsigned int i) const

G4String	getAbsMat () const

int	getNAbsorber () const

double	getDy_2Abs () const

double	getDx_2Abs (unsigned int i) const

G4String	getScntMat () const

G4String	getWrapMat () const

G4String	getPlasMat () const

int	getNScintillator () const

double	getDy_2ScntLay (unsigned int i) const

double	getDx_2ScntLay (unsigned int i) const

double	getDx_2Wrap (unsigned int i) const

double	getDx_2FrontP (unsigned int i) const

double	getDx_2BackP (unsigned int i) const

double	getDx_2Scnt (unsigned int i) const

Public Member Functions inherited from CCalDetector
	CCalDetector (const G4String &name)

virtual	~CCalDetector ()

void	constructHierarchy ()

void	construct ()

void	addDetector (CCalDetector *)

G4String	Name () const

G4String	baseFileName () const

G4String	File () const

CCalDetector *	getDaughter (int i) const

int	getNDaughters () const

G4bool	operator== (const CCalDetector &left) const

G4bool	operator!= (const CCalDetector &left) const

Public Member Functions inherited from CCalG4Able
	CCalG4Able (G4String name)

virtual	~CCalG4Able ()

G4VPhysicalVolume *	PhysicalVolume (G4VPhysicalVolume *)

void	setVisType (CCalVisualisable::visType, G4LogicalVolume *)

void	setSensitivity (G4bool sens=true)

G4bool	isSensitive () const

G4String	G4Name () const

void	setName (const G4String &name)

G4bool	operator== (const CCalG4Able &right) const

Protected Member Functions
virtual G4VPhysicalVolume *	constructIn (G4VPhysicalVolume *)

virtual void	constructDaughters ()

G4LogicalVolume *	constructScintillatorLayer (G4int)

G4LogicalVolume *	constructAbsorberLayer (G4int)

virtual void	constructSensitive ()

Protected Member Functions inherited from CCalHcal
virtual int	readFile ()

Protected Member Functions inherited from CCalDetector
int	buildFromFile ()

Protected Member Functions inherited from CCalG4Able
void	AddCCalG4Able (CCalG4Able *)

Additional Inherited Members
Protected Attributes inherited from CCalDetector
G4String	detectorName

G4String	fileName

CCalDetectorTable	theDetectorsInside

int	constructFlag

Protected Attributes inherited from CCalG4Able
G4VPhysicalVolume *	detPhysicalVolume

CCalG4AbleTable	theG4DetectorsInside

G4String	g4ableName

G4bool	sensitivity

CCalVisualisable	visProperties

G4VisAttributes *	g4VisAtt [CCalVisualisable::TotalVisTypes]

Static Protected Attributes inherited from CCalDetector
static G4String	pathName = getenv("CCAL_GEOMPATH")

Detailed Description

Definition at line 40 of file CCalG4Hcal.hh.

Constructor & Destructor Documentation

CCalG4Hcal::CCalG4Hcal ( const G4String & name )

Definition at line 57 of file CCalG4Hcal.cc.

57 :

58 CCalHcal(name), CCalG4Able(name), sclLog(0), absLog(0) {}

CCalG4Able::CCalG4Able

CCalG4Able(G4String name)

Definition: CCalG4Able.cc:49

CCalHcal::CCalHcal

CCalHcal(const G4String &name)

Definition: CCalHcal.cc:37

CCalG4Hcal::~CCalG4Hcal ( )

virtual

Definition at line 60 of file CCalG4Hcal.cc.

                        {
   if (sclLog)
     delete[] sclLog;
   if (absLog)
     delete[] absLog;
 }

Member Function Documentation

G4LogicalVolume * CCalG4Hcal::constructAbsorberLayer ( G4int lay )

protected

Definition at line 320 of file CCalG4Hcal.cc.

References CCalVisualisable::Absorber, CCalMaterialFactory::findMaterial(), G4cout, G4endl, CCalHcal::getAbsMat(), CCalHcal::getDx_2Abs(), CCalHcal::getDy_2Abs(), CCalMaterialFactory::getInstance(), python.hepunit::mm, CCalDetector::Name(), CCalG4Able::setVisType(), and tab().

Referenced by constructIn().

                                                              {
 
   //Pointers to the Materials
   CCalMaterialFactory* matfact       = CCalMaterialFactory::getInstance();
 
   //Now the absorber layer
   G4Material* matter = matfact->findMaterial(getAbsMat());
   G4String    name   = Name() + "Absorber" + lay;
   G4VSolid*   solid  = new G4Box (name, getDx_2Abs(lay)*mm, getDy_2Abs()*mm,
                   getDy_2Abs()*mm);
   G4LogicalVolume* log = new G4LogicalVolume(solid, matter, name);
   setVisType(CCalVisualisable::Absorber,log);
 #ifdef debug
   G4cout << tab << name << " Box made of " << getAbsMat() << " of dimension " 
        << getDx_2Abs(lay)*mm << " " << getDy_2Abs()*mm << " " 
        << getDy_2Abs()*mm << G4endl;
 #endif
 
   return log;
 }

void CCalG4Hcal::constructDaughters ( )

protectedvirtual

Reimplemented from CCalHcal.

Definition at line 361 of file CCalG4Hcal.cc.

361 {}

G4VPhysicalVolume * CCalG4Hcal::constructIn ( G4VPhysicalVolume * mother )

protectedvirtual

Implements CCalG4Able.

Definition at line 71 of file CCalG4Hcal.cc.

References constructAbsorberLayer(), constructScintillatorLayer(), CCalMaterialFactory::findMaterial(), G4cout, G4endl, CCalHcal::getBoxMat(), CCalHcal::getDx_2Box(), CCalHcal::getDx_2Cal(), CCalHcal::getDy_2Box(), CCalHcal::getDy_2Cal(), CCalHcal::getGenMat(), CCalMaterialFactory::getInstance(), CCalHcal::getMotherAbs(), CCalHcal::getMotherScnt(), CCalHcal::getNAbsorber(), G4VPhysicalVolume::GetName(), G4LogicalVolume::GetName(), CCalHcal::getNBox(), CCalHcal::getNLayerAbs(), CCalHcal::getNLayerScnt(), CCalHcal::getNScintillator(), CCalHcal::getTypeAbs(), CCalHcal::getTypeScnt(), CCalHcal::getWallThickBox(), CCalHcal::getXposAbs(), CCalHcal::getXposBox(), CCalHcal::getXposCal(), CCalHcal::getXposScnt(), python.hepunit::mm, CCalDetector::Name(), CCalVisualisable::PseudoVolumes, CCalG4Able::setVisType(), CCalVisualisable::Support, and tab().

                                                                     {
   G4cout << "==>> Constructing CCalG4Hcal..." << G4endl;
 
   //Common logical volumes between methods.
 #ifdef debug
   G4cout << tab << "Common logical volumes initialization: " 
        << getNScintillator() << " scintillaor and " << getNAbsorber()
        << " absorber layers." << G4endl;
 #endif
   G4int i = 0;
   sclLog  = new ptrG4Log[getNScintillator()];
   absLog  = new ptrG4Log[getNAbsorber()];
   for (i=0; i < getNScintillator(); i++)
     sclLog[i]  = 0;
   for (i=0; i < getNAbsorber(); i++)
     absLog[i] = 0;
 
   //Pointers to the Materials
   CCalMaterialFactory* matfact       = CCalMaterialFactory::getInstance();
 
   //Mother volume
   G4Material* matter = matfact->findMaterial(getGenMat());
   G4VSolid*   solid  = new G4Box (Name(), getDx_2Cal()*mm, getDy_2Cal()*mm,
                   getDy_2Cal()*mm);
   G4LogicalVolume* logh = new G4LogicalVolume(solid, matter, Name());
   setVisType(CCalVisualisable::PseudoVolumes,logh);
 #ifdef debug
     G4cout << tab << Name() << " Box made of " << getGenMat()
          << " of dimension " << getDx_2Cal()*mm << " " << getDy_2Cal()*mm
          << " " << getDy_2Cal()*mm << G4endl;
 #endif
 
   G4PVPlacement* hcal = new G4PVPlacement(0,G4ThreeVector(getXposCal()*mm,0,0),
                       Name(), logh, mother, false, 1);
   G4String name("Null");
 #ifdef pdebug
   if (mother != 0) name = mother->GetName();
   G4cout << Name() << " Number 1 positioned in " << name << " at ("
        << getXposCal()*mm << ",0,0) with no rotation" << G4endl;
 #endif
 
   //Wall of the Boxes
   solid  = new G4Box (name, 0.5*getWallThickBox()*mm, getDy_2Box()*mm, 
               getDy_2Box()*mm);
   matter = matfact->findMaterial(getBoxMat());
   name   = Name() + "Wall";
   G4LogicalVolume* logw = new G4LogicalVolume(solid, matter, name);
   setVisType(CCalVisualisable::Support,logw);
 #ifdef debug
   G4cout << tab << name << " Box made of " << getBoxMat()
        << " of dimension " << 0.5*getWallThickBox()*mm << " " 
        << getDy_2Box()*mm << " " << getDy_2Box()*mm << G4endl;
 #endif
 
   //Now the boxes
   ptrG4Log* logb = new ptrG4Log[getNBox()];
   matter = matfact->findMaterial(getGenMat());
   for (i=0; i<getNBox(); i++) {
     name   = Name() + "Box" + i;
     solid  = new G4Box (name, getDx_2Box()*mm, getDy_2Box()*mm, 
             getDy_2Box()*mm);
     logb[i]= new G4LogicalVolume(solid, matter, name);
     setVisType(CCalVisualisable::PseudoVolumes,logb[i]);
 #ifdef debug
     G4cout << tab << name << " Box made of " << getGenMat()
      << " of dimension " << getDx_2Box()*mm << " " << getDy_2Box()*mm
          << " " << getDy_2Box()*mm << G4endl;
 #endif
 
     G4double xpos = -(getDx_2Box() - 0.5*getWallThickBox());
     new G4PVPlacement (0, G4ThreeVector(xpos*mm,0,0), logw, logw->GetName(), 
                logb[i], false, 1);
 #ifdef pdebug
     G4cout << logw->GetName() << " Number 1 positioned in " << name
      << " at (" << xpos*mm << ",0,0) with no rotation" << G4endl;
 #endif
     xpos = (getDx_2Box() - 0.5*getWallThickBox());
     new G4PVPlacement (0, G4ThreeVector(xpos*mm,0,0), logw, logw->GetName(), 
                logb[i], false, 2);
 #ifdef pdebug
     G4cout << logw->GetName() << " Number 2 positioned in " << name
      << " at (" << xpos*mm << ",0,0) with no rotation" << G4endl;
 #endif
 
     new G4PVPlacement (0, G4ThreeVector(getXposBox(i)*mm,0,0), logb[i], name, 
                logh, false, i+1);
 #ifdef pdebug
     G4cout << name << " Number " << i+1 << " positioned in " << logh->GetName()
      << " at (" << getXposBox(i)*mm << ",0,0) with no rotation" << G4endl;
 #endif
   }
 
   //Loop over scintillator layers
   for (i=0; i<getNLayerScnt(); i++) {
     G4int lay = getTypeScnt(i);
     if (!sclLog[lay])
       sclLog[lay] = constructScintillatorLayer(lay);
     if (getMotherScnt(i) < 0 || getMotherScnt(i) >= getNScintillator()) {
       logw = logh;
     } else {
       logw = logb[getMotherScnt(i)];
     }
     G4double xpos = getXposScnt(i);
     new G4PVPlacement (0, G4ThreeVector(xpos*mm,0,0), sclLog[lay], 
                sclLog[lay]->GetName(), logw, false, i+1);
 #ifdef pdebug
     G4cout << sclLog[lay]->GetName() << " Number " << i+1 << " positioned in " 
      << logw->GetName() << " at (" << xpos*mm << ",0,0) with no rotation" 
      << G4endl;
 #endif
   }
 
   //Loop over absorber layers
   for (i=0; i<getNLayerAbs(); i++) {
     G4int lay = getTypeAbs(i);
     if (!absLog[lay])
       absLog[lay] = constructAbsorberLayer(lay);
     if (getMotherAbs(i) < 0 || getMotherAbs(i) >= getNAbsorber()) {
       logw = logh;
     } else {
       logw = logb[getMotherAbs(i)];
     }
     G4double xpos = getXposAbs(i);
     new G4PVPlacement (0, G4ThreeVector(xpos*mm,0,0), absLog[lay], 
                absLog[lay]->GetName(), logw, false, i+1);
 #ifdef pdebug
     G4cout << absLog[lay]->GetName() << " Number " << i+1 << " positioned in " 
      << logw->GetName() << " at (" << xpos*mm << ",0,0) with no rotation" 
      << G4endl;
 #endif
   }
 
   delete [] logb;
 
   G4cout << "<<== End of CCalG4Hcal construction ..." << G4endl;
 
   return hcal;
 }

G4LogicalVolume * CCalG4Hcal::constructScintillatorLayer ( G4int lay )

protected

Definition at line 211 of file CCalG4Hcal.cc.

References CCalVisualisable::Cable, CCalMaterialFactory::findMaterial(), G4cout, G4endl, CCalHcal::getDx_2BackP(), CCalHcal::getDx_2FrontP(), CCalHcal::getDx_2Scnt(), CCalHcal::getDx_2ScntLay(), CCalHcal::getDx_2Wrap(), CCalHcal::getDy_2ScntLay(), CCalHcal::getGenMat(), CCalMaterialFactory::getInstance(), G4LogicalVolume::GetName(), CCalHcal::getPlasMat(), CCalHcal::getScntMat(), CCalHcal::getWrapMat(), python.hepunit::mm, CCalDetector::Name(), CCalVisualisable::PseudoVolumes, CCalVisualisable::Sensitive, CCalG4Able::setVisType(), CCalVisualisable::Support, and tab().

Referenced by constructIn().

                                                                  {
 
   //Pointers to the Materials
   CCalMaterialFactory* matfact       = CCalMaterialFactory::getInstance();
 
   //The scintillator layer
   G4Material* matter = matfact->findMaterial(getGenMat());
   G4String    name   = Name() + "ScntLayer" + lay;
   G4VSolid*   solid  = new G4Box (name, getDx_2ScntLay(lay)*mm, 
                   getDy_2ScntLay(lay)*mm,
                   getDy_2ScntLay(lay)*mm);
   G4LogicalVolume* log = new G4LogicalVolume(solid, matter, name);
   setVisType(CCalVisualisable::PseudoVolumes,log);
 #ifdef debug
   G4cout << tab << name << " Box made of " << getGenMat() << " of dimension " 
        << getDx_2ScntLay(lay)*mm << " " << getDy_2ScntLay(lay)*mm << " " 
        << getDy_2ScntLay(lay)*mm << G4endl;
 #endif
 
   G4LogicalVolume* logd;
   G4double         xpos;
   //Wrappers if any
   if (getDx_2Wrap(lay) > 0) {
     name   = Name() + "ScntWrapper" + lay;
     matter = matfact->findMaterial(getWrapMat());
     solid  = new G4Box (name, getDx_2Wrap(lay)*mm, 
             getDy_2ScntLay(lay)*mm, getDy_2ScntLay(lay)*mm);
     logd   = new G4LogicalVolume(solid, matter, name);
     setVisType(CCalVisualisable::Support,logd);
 #ifdef debug
     G4cout << tab << name << " Box made of " << getWrapMat() << " of dimension " 
      << getDx_2Wrap(lay)*mm << " " << getDy_2ScntLay(lay)*mm << " " 
      << getDy_2ScntLay(lay)*mm << G4endl;
 #endif
     xpos   =-(getDx_2ScntLay(lay)-getDx_2Wrap(lay));
     new G4PVPlacement(0, G4ThreeVector(xpos*mm,0,0), logd, name, log, false,1);
 #ifdef pdebug
     G4cout << logd->GetName() << " Number 1 positioned in " << log->GetName() 
      << " at (" << xpos*mm << ",0,0) with no rotation" << G4endl;
 #endif
     xpos   = (getDx_2ScntLay(lay)-getDx_2Wrap(lay));
     new G4PVPlacement(0, G4ThreeVector(xpos*mm,0,0), logd, name, log, false,2);
 #ifdef pdebug
     G4cout << logd->GetName() << " Number 2 positioned in " << log->GetName() 
      << " at (" << xpos*mm << ",0,0) with no rotation" << G4endl;
 #endif
   }
 
   //Plastic covers
   matter = matfact->findMaterial(getPlasMat());
   name   = Name() + "FrontPlastic" + lay;
   solid  = new G4Box (name, getDx_2FrontP(lay)*mm, getDy_2ScntLay(lay)*mm, 
               getDy_2ScntLay(lay)*mm);
   logd   = new G4LogicalVolume(solid, matter, name);
   setVisType(CCalVisualisable::Cable,logd);
 #ifdef debug
   G4cout << tab << name << " Box made of " << getPlasMat() << " of dimension " 
        << getDx_2FrontP(lay)*mm << " " << getDy_2ScntLay(lay)*mm << " " 
        << getDy_2ScntLay(lay)*mm << G4endl;
 #endif
   xpos   =-getDx_2ScntLay(lay)+2.*getDx_2Wrap(lay)+getDx_2FrontP(lay);
   new G4PVPlacement(0, G4ThreeVector(xpos*mm,0,0), logd, name, log, false,1);
 #ifdef pdebug
   G4cout << logd->GetName() << " Number 1 positioned in " << log->GetName() 
        << " at (" << xpos*mm << ",0,0) with no rotation" << G4endl;
 #endif
   name   = Name() + "BackPlastic" + lay;
   solid  = new G4Box (name, getDx_2BackP(lay)*mm, getDy_2ScntLay(lay)*mm, 
               getDy_2ScntLay(lay)*mm);
   logd   = new G4LogicalVolume(solid, matter, name);
   setVisType(CCalVisualisable::Cable,logd);
 #ifdef debug
   G4cout << tab << name << " Box made of " << getPlasMat() << " of dimension " 
        << getDx_2BackP(lay)*mm << " " << getDy_2ScntLay(lay)*mm << " " 
        << getDy_2ScntLay(lay)*mm << G4endl;
 #endif
   xpos   =(-getDx_2ScntLay(lay)+2.*getDx_2Wrap(lay)+2.*getDx_2FrontP(lay)+
            2.*getDx_2Scnt(lay)+getDx_2BackP(lay));
   new G4PVPlacement(0, G4ThreeVector(xpos*mm,0,0), logd, name, log, false,1);
 #ifdef pdebug
   G4cout << logd->GetName() << " Number 1 positioned in " << log->GetName() 
        << " at (" << xpos*mm << ",0,0) with no rotation" << G4endl;
 #endif
 
   //Now the scintillators
   matter = matfact->findMaterial(getScntMat());
   name   = Name() + "Scintillator" + lay;
   solid  = new G4Box (name, getDx_2Scnt(lay)*mm, getDy_2ScntLay(lay)*mm, 
               getDy_2ScntLay(lay)*mm);
   logd   = new G4LogicalVolume(solid, matter, name);
   setVisType(CCalVisualisable::Sensitive,logd);
   allSensitiveLogs.push_back(logd);
 #ifdef debug
   G4cout << tab << name << " Box made of " << getScntMat() << " of dimension " 
        << getDx_2Scnt(lay)*mm << " " << getDy_2ScntLay(lay)*mm << " " 
        << getDy_2ScntLay(lay)*mm << G4endl;
 #endif
   xpos   =(-getDx_2ScntLay(lay)+2.*getDx_2Wrap(lay)+2.*getDx_2FrontP(lay)+
        getDx_2Scnt(lay));
   new G4PVPlacement(0, G4ThreeVector(xpos*mm,0,0), logd, name, log, false,1);
 #ifdef pdebug
   G4cout << logd->GetName() << " Number 1 positioned in " << log->GetName() 
        << " at (" << xpos*mm << ",0,0) with no rotation" << G4endl;
 #endif
 
   return log;
 }

void CCalG4Hcal::constructSensitive ( )

protectedvirtual

Reimplemented from CCalG4Able.

Definition at line 342 of file CCalG4Hcal.cc.

References G4cerr, G4cout, G4endl, CCalSensitiveDetectors::getInstance(), CCalDetector::Name(), and CCalSensitiveDetectors::registerVolume().

                                    {
 
   if (allSensitiveLogs.size()>0) {
     CCalSensitiveDetectors* sensDets = CCalSensitiveDetectors::getInstance();
     G4String SDname = Name();
     for (std::vector<ptrG4Log>::iterator iter=allSensitiveLogs.begin(); 
          iter<allSensitiveLogs.end(); iter++) {
       sensDets->registerVolume(SDname, (*iter));
 #ifdef sdebug
       G4cout << "Register volume " << (*iter)->GetName() << " for" << SDname 
        << G4endl;
 #endif
     }
   } else {
     G4cerr << "CCalG4Hcal ERROR: Could not construct Sensitive Detector" 
        << G4endl;
   }
 }

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation