CCalG4Ecal Class Reference

#include <CCalG4Ecal.hh>

Inheritance diagram for CCalG4Ecal:

Collaboration diagram for CCalG4Ecal:

Public Types
enum	CMType { module1, module2 }

Public Member Functions
	CCalG4Ecal (const G4String &name)
virtual	~CCalG4Ecal ()
void	setType (CMType ty)
Public Member Functions inherited from CCalEcal
	CCalEcal (const G4String &name)
virtual	~CCalEcal ()
G4String	getGenMat () const
double	getWidBox () const
double	getLengBox () const
double	getXpos () const
double	getYpos () const
double	getZpos () const
double	getThetaX () const
double	getPhiX () const
double	getThetaY () const
double	getPhiY () const
double	getThetaZ () const
double	getPhiZ () const
G4String	getLayMat () const
int	getLayNum () const
double	getLayRadius () const
double	getLayAngle () const
double	getLengFront () const
double	getLayPar (unsigned int i) const
G4String	getCrystMat () const
int	getCrystNum () const
double	getCrystLength () const
double	getCrystTol () const
double	getCrystPar (unsigned int i) const
G4String	getSuppMat () const
double	getDxSupp () const
double	getDySupp () const
double	getDzSupp () const
double	getDistSupp () 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

Static Public Attributes
static G4String	idName = "CrystalMatrix"

Protected Member Functions
virtual G4VPhysicalVolume *	constructIn (G4VPhysicalVolume *)
virtual void	constructSensitive ()
Protected Member Functions inherited from CCalEcal
virtual int	readFile ()
virtual void	constructDaughters ()
Protected Member Functions inherited from CCalDetector
int	buildFromFile ()
Protected Member Functions inherited from CCalG4Able
void	AddCCalG4Able (CCalG4Able *)

Private Member Functions
G4LogicalVolume *	constructGlobal ()

Private Attributes
CMType	type
std::vector< ptrG4Log >	sensitiveLogs

Static Private Attributes
static G4LogicalVolume *	crystalmatrixLog = 0

Additional Inherited Members
Protected Attributes inherited from CCalDetector
G4String	detectorName
G4String	fileName
G4String	pathName
CCalDetectorTable	theDetectorsInside
int	constructFlag
Protected Attributes inherited from CCalG4Able
G4VPhysicalVolume *	detPhysicalVolume
CCalG4AbleTable	theG4DetectorsInside
G4String	g4ableName
G4bool	sensitivity
CCalVisualisable	visProperties
G4VisAttributes *	g4VisAtt [CCalVisualisable::TotalVisTypes]

Detailed Description

Definition at line 39 of file CCalG4Ecal.hh.

Member Enumeration Documentation

CMType

enum CCalG4Ecal::CMType

Enumerator
module1
module2

Definition at line 42 of file CCalG4Ecal.hh.

{module1, module2};

Constructor & Destructor Documentation

CCalG4Ecal()

CCalG4Ecal::CCalG4Ecal

(

const G4String &

name

)

Definition at line 65 of file CCalG4Ecal.cc.

                                          :
  CCalEcal(name), CCalG4Able(name), type(module1) {}

~CCalG4Ecal()

CCalG4Ecal::~CCalG4Ecal ( )

virtual

Definition at line 68 of file CCalG4Ecal.cc.

{}

Member Function Documentation

constructGlobal()

G4LogicalVolume * CCalG4Ecal::constructGlobal ( )

private

Definition at line 149 of file CCalG4Ecal.cc.

                                             {

  //Pointers to the Materials and Rotation Matrix factory
  CCalMaterialFactory* matfact       = CCalMaterialFactory::getInstance();
  CCalRotationMatrixFactory* rotfact = CCalRotationMatrixFactory::getInstance();
  
  G4Material* matter = matfact->findMaterial(getGenMat());
  G4VSolid* solid = new G4Box (idName, 0.5*getWidBox()*mm, 0.5*getWidBox()*mm,
                   0.5*getLengBox()*mm);
#ifdef debug
  G4cout << tab << idName << " Box made of " << getGenMat() << " of dimension " 
       << 0.5*getWidBox()*mm << ", " << 0.5*getWidBox()*mm << ", "
       << 0.5*getLengBox()*mm << G4endl;
#endif
  G4LogicalVolume* glog = new G4LogicalVolume (solid, matter, idName);
  setVisType(CCalVisualisable::PseudoVolumes,glog);

  //Now the layers
  G4String name = idName + "Layer";
  matter = matfact->findMaterial(getLayMat());
  solid  = new G4Trd(name, getLayPar(0)*mm, getLayPar(1)*mm, getLayPar(2)*mm, 
             getLayPar(3)*mm, getLayPar(4)*mm);
#ifdef debug
  G4cout << tab << name << " Trd made of " << getLayMat() << " of dimension " 
       << getLayPar(0)*mm << ", " << getLayPar(1)*mm << ", " << getLayPar(2)*mm
       << ", " << getLayPar(3)*mm << ", " << getLayPar(4)*mm << G4endl;
#endif
  G4LogicalVolume* laylog = new G4LogicalVolume (solid, matter, name);
  setVisType(CCalVisualisable::OtherServices,laylog);

  G4int i = 0;
  G4String rotstr;
  G4double xp, yp, zp, angle;
  G4double zshift = -0.5 * (getLengBox() - getCrystLength()) + getLengFront();
  G4RotationMatrix* rot = 0;
  for (i = 0; i < getLayNum(); i++) {
    angle  = 0.5 * getLayAngle() * (2*i + 1 - getLayNum());
    xp     = angle * (getLayPar(4) + getLayRadius()) * mm;
    zp     = (zshift + getLayPar(0)*std::abs(std::sin(angle))) * mm;
    rotstr = idName + "Layer" + i;
    rot    = rotfact->findMatrix(rotstr);
    if (!rot) {
#ifdef ddebug
      G4cout << "Creating a new rotation: " << rotstr << tab 
       << (90.0*deg+angle) << "," << 0.0*deg << "," << 90.0*deg << "," 
       << 90.0*deg << "," << angle << "," << 0.0*deg << G4endl;
#endif
      rot = rotfact->AddMatrix(rotstr, (90.0*deg+angle), 0.0*deg, 90.0*deg,
                   90.0*deg, angle, 0.0*deg);
    }
    new G4PVPlacement(rot, G4ThreeVector(xp,0.,zp), laylog, name, glog,
              false, i+1);
#ifdef pdebug
    G4cout << laylog->GetName() << " number " << i+1 << " positioned in " 
         << glog->GetName()  << " at (" << xp << ", 0," << zp
     << ") with rotation angle " << angle/deg << G4endl;
#endif
  }

  //Now the crystals
  name   = idName + "Crystal";
  matter = matfact->findMaterial(getCrystMat());
  solid  = new G4Trd(name, getCrystPar(0)*mm, getCrystPar(1)*mm, 
             getCrystPar(2)*mm, getCrystPar(3)*mm, getCrystPar(4)*mm);
#ifdef debug
  G4cout << tab << name << " Trd made of " << getCrystMat() << " of dimension " 
       << getCrystPar(0)*mm << ", " << getCrystPar(1)*mm << ", " 
       << getCrystPar(2)*mm << ", " << getCrystPar(3)*mm << ", " 
       << getCrystPar(4)*mm << G4endl;
#endif

  G4LogicalVolume* detLog = new G4LogicalVolume (solid, matter, name);
  setVisType(CCalVisualisable::Sensitive,detLog);
  sensitiveLogs.push_back(detLog);
  for (i = 0; i < getCrystNum(); i++) {
    angle  = 0.5 * getLayAngle() * (2*i + 1 - getCrystNum());
    yp     = angle * (getCrystPar(4) + getLayRadius()) * mm;
    zp     = (getCrystPar(0)*std::abs(std::sin(angle)) - getCrystTol()) * mm;
    rotstr = idName + "Crystal" + i;
    rot    = rotfact->findMatrix(rotstr);
    if (!rot) {
#ifdef ddebug
      G4cout << "Creating a new rotation: " << rotstr << tab << 90.0*deg << ","
       << 0.0*deg << "," << (90.0*deg+angle) << "," << 0.0*deg << "," 
       << angle << "," << 90.0*deg << G4endl;
#endif
      rot = rotfact->AddMatrix(rotstr, 90.0*deg, 0.0*deg, (90.0*deg+angle),
                   90.0*deg, angle, 90.0*deg);
    }
    new G4PVPlacement(rot, G4ThreeVector(0,yp,zp), detLog, name, laylog,
              false, i+1);
#ifdef pdebug
    G4cout << detLog->GetName() << " number " << i+1 << " positioned in " 
         << laylog->GetName()  << " at (0," << yp << "," << zp
     << ") with rotation angle " << angle/deg << G4endl;
#endif
  }

  //Support boxes
  name   = idName + "Support";
  matter = matfact->findMaterial(getSuppMat());
  solid  = new G4Box (name, 0.5*getDxSupp()*mm, 0.5*getDySupp()*mm, 
              0.5*getDzSupp()*mm);
#ifdef debug
  G4cout << tab << name << " Box made of " << getSuppMat() << " of dimension " 
       << 0.5*getDxSupp()*mm << ", " << 0.5*getDySupp()*mm << ", "
       << 0.5*getDzSupp()*mm << G4endl;
#endif
  G4LogicalVolume* slog = new G4LogicalVolume (solid, matter, name);
  setVisType(CCalVisualisable::Support,slog);

  zp   = (-0.5 * getLengBox() + getCrystLength() + getLengFront() +
      0.5 * getDzSupp() + getDistSupp()) * mm;
  for (i = 0; i < getCrystNum(); i++) {
    yp   = getLayPar(1) * (2*i + 1 - getCrystNum()) * mm;
    new G4PVPlacement(0, G4ThreeVector(0,yp,zp), slog, name, glog,
              false, i+1);
#ifdef pdebug
    G4cout << slog->GetName() << " number " << i+1 << " positioned in " 
         << glog->GetName()  << " at (0," << yp << "," << zp
     << ") with no rotation" << G4endl;
#endif
  }

  return glog;
}

Here is the call graph for this function:

Here is the caller graph for this function:

constructIn()

G4VPhysicalVolume * CCalG4Ecal::constructIn ( G4VPhysicalVolume *  mother )

protectedvirtual

Implements CCalG4Able.

Definition at line 74 of file CCalG4Ecal.cc.

                                                                    {
  G4cout << "==>> Constructing CCalG4Ecal..." << G4endl;

  // Construction of global volume as a Box

  if (!crystalmatrixLog) {
    crystalmatrixLog = constructGlobal();
  }
  CCalRotationMatrixFactory* rotfact = CCalRotationMatrixFactory::getInstance();

  G4double x, y, z;
  if (mother != 0) {
    x = getXpos()*mm;
    y = getYpos()*mm;
    z = getZpos()*mm;
  } else {
    x = y = z = 0;
  }
    
  int num;
  if (type == module2) {
    num = 2;
  } else {
    num = 1;
  }
#ifdef pdebug
  G4String name("Null");
  if (mother != 0) name = mother->GetName();
  G4cout << crystalmatrixLog->GetName() << " Number " << num << " positioned in "
       << name << " at (" << x << ", " << y << ", " << z << ")";
#endif

  G4RotationMatrix* cmrot = 0;
  if (mother != 0) {
    G4String rotstr = idName + num;
    cmrot  = rotfact->findMatrix(rotstr);
    if (!cmrot) {
#ifdef ddebug
      G4cout << "Creating a new rotation: " << rotstr << tab 
       << getThetaX()*deg << "," << getPhiX()*deg << "," 
       << getThetaY()*deg << "," << getPhiY()*deg << "," 
       << getThetaZ()*deg << "," << getPhiZ()*deg << G4endl;
#endif
      cmrot = rotfact->AddMatrix(rotstr, getThetaX()*deg, getPhiX()*deg, 
                 getThetaY()*deg, getPhiY()*deg,
                 getThetaZ()*deg, getPhiZ()*deg);
    } // if !cmrot
#ifdef pdebug
    G4cout << " rotation by (" <<  getThetaX() << ", " << getPhiX() << ", " 
     << getThetaY() << "," << getPhiY() << ", "  << getThetaZ() << ", " 
     << getPhiZ() << ")" << G4endl;
#endif
  } else {
#ifdef pdebug
    G4cout << " without rotation..." << G4endl;
#endif
  }

  G4PVPlacement* crystalmatrix;
  if (mother != 0) {
    crystalmatrix = new G4PVPlacement(cmrot, G4ThreeVector(x,y,z),
                      crystalmatrixLog, idName,
                      mother->GetLogicalVolume(), false, num);
  } else {
    crystalmatrix = new G4PVPlacement(cmrot, G4ThreeVector(x,y,z),
                      idName, crystalmatrixLog,
                      mother, false, num);
  }
  G4cout << "<<== End of CCalG4Ecal construction ..." << G4endl;

  return crystalmatrix;
}

Here is the call graph for this function:

constructSensitive()

void CCalG4Ecal::constructSensitive ( )

protectedvirtual

Reimplemented from CCalG4Able.

Definition at line 276 of file CCalG4Ecal.cc.

                                    {

#ifdef debug
  G4cout << "Now registering CrystalMatrix LogicalVolume's to SD's:" << G4endl;
#endif
  if (sensitiveLogs.size()>0) {
    CCalSensitiveDetectors* sensDets = CCalSensitiveDetectors::getInstance();
    G4String SDname = idName;
    for(std::vector<ptrG4Log>::iterator iter=sensitiveLogs.begin(); 
                               iter<sensitiveLogs.end(); iter++) {
      sensDets->registerVolume(SDname, (*iter));
#ifdef sdebug
      G4cout << "Register volume " << (*iter)->GetName() << " for" << SDname 
       << G4endl;
#endif
    }
  }

}

Here is the call graph for this function:

setType()

void CCalG4Ecal::setType ( CMType  ty )

inline

Definition at line 48 of file CCalG4Ecal.hh.

{type = ty;}

Here is the caller graph for this function:

Member Data Documentation

crystalmatrixLog

G4LogicalVolume * CCalG4Ecal::crystalmatrixLog = 0

staticprivate

Definition at line 70 of file CCalG4Ecal.hh.

idName

G4String CCalG4Ecal::idName = "CrystalMatrix"

static

Definition at line 51 of file CCalG4Ecal.hh.

sensitiveLogs

std::vector<ptrG4Log> CCalG4Ecal::sensitiveLogs

private

Definition at line 73 of file CCalG4Ecal.hh.

type

CMType CCalG4Ecal::type

private

Definition at line 67 of file CCalG4Ecal.hh.

---

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

• CCalG4Ecal.hh
• CCalG4Ecal.cc