ExtBesCrystalParameterisation Class Reference

#include <ExtBesCrystalParameterisation.h>

Inheritance diagram for ExtBesCrystalParameterisation:

Public Member Functions
	ExtBesCrystalParameterisation (const G4int, const G4int, const G4int, ExtBesEmcGeometry *, const G4int)
virtual	~ExtBesCrystalParameterisation ()
void	ComputeTransformation (const G4int, G4VPhysicalVolume *) const
void	ComputeDimensions (G4Trap &, const G4int, const G4VPhysicalVolume *) const
G4Material *	ComputeMaterial (const G4int, G4VPhysicalVolume *)
G4VSolid *	ComputeSolid (const G4int, G4VPhysicalVolume *)
void	ComputeIDAndSide (const G4VPhysicalVolume *) const
void	SetVerboseLevel (G4int val)

Detailed Description

Definition at line 23 of file ExtBesCrystalParameterisation.h.

Constructor & Destructor Documentation

ExtBesCrystalParameterisation()

ExtBesCrystalParameterisation::ExtBesCrystalParameterisation	(	const G4int	startID,
		const G4int	endID,
		const G4int	thetaAllCrystals,
		ExtBesEmcGeometry *	besEMCGeometry,
		const G4int	verboseLevel )

Definition at line 23 of file ExtBesCrystalParameterisation.cxx.

         :fVerboseLevel(0)
{
  //for debug
  //G4Exception("ExtBesCrystalParameterisation::ExtBesCrystalParameterisation() starting........");
  fBesEmcGeometry   =  besEMCGeometry;
  fVerboseLevel     =  verboseLevel; //if need to debug this class,open it
  fStartID          =  startID;
  fAllCrystals      =  thetaAllCrystals;
  fCrystalID        =  0;
  fFlagLeft         =  true;
  if (startID>endID||endID>thetaAllCrystals||startID<1)
    {
      G4Exception("ExtBesCrystalParameterisation construction: ID of crystal>No avaible!");
    }
 
}

~ExtBesCrystalParameterisation()

ExtBesCrystalParameterisation::~ExtBesCrystalParameterisation ( )

virtual

Definition at line 48 of file ExtBesCrystalParameterisation.cxx.

{}

Member Function Documentation

ComputeDimensions()

void ExtBesCrystalParameterisation::ComputeDimensions	(	G4Trap &	trackerChamber,
		const G4int	copyNo,
		const G4VPhysicalVolume *	physVol ) const

Definition at line 110 of file ExtBesCrystalParameterisation.cxx.

{
  //for debug
  //G4Exception("ExtBesCrystalParameterisation::ComputeDimensions() starting........");
  ComputeIDAndSide(physVol);
 
  G4double zHalfLength  = fBesEmcGeometry->GetZHalfLength(fCrystalID);
  G4double thetaAxis    = fBesEmcGeometry->GetThetaAxis(fCrystalID);
  G4double phiAxis      = fBesEmcGeometry->GetPhiAxis(fCrystalID);
  G4double yHalfLength1 = fBesEmcGeometry->GetYHalfLength1(fCrystalID);
  G4double xHalfLength1 = fBesEmcGeometry->GetXHalfLength1(fCrystalID);
  G4double xHalfLength2 = fBesEmcGeometry->GetXHalfLength2(fCrystalID);
  G4double tanAlpha1    = fBesEmcGeometry->GetTanAlpha1(fCrystalID);
  G4double yHalfLength2 = fBesEmcGeometry->GetYHalfLength2(fCrystalID);
  G4double xHalfLength3 = fBesEmcGeometry->GetXHalfLength3(fCrystalID);
  G4double xHalfLength4 = fBesEmcGeometry->GetXHalfLength4(fCrystalID);
  G4double tanAlpha2    = fBesEmcGeometry->GetTanAlpha2(fCrystalID);
 
  if(!fFlagLeft)
    {
      phiAxis=-phiAxis;
      tanAlpha1=-tanAlpha1;
      tanAlpha2=-tanAlpha2;
      G4double tmp;
      tmp=xHalfLength1;
      xHalfLength1=xHalfLength2;
      xHalfLength2=tmp;
      tmp=xHalfLength3;
      xHalfLength3=xHalfLength4;
      xHalfLength4=tmp;
    }
  if(fVerboseLevel>5)
    G4cout << "The size of No." << copyNo << " crystal(placed) are:" << G4endl
       << "zHalfLength =" << zHalfLength/cm  << "(cm), "     << G4endl
       << "thetaAxis   =" << thetaAxis/deg   << "(degree), " << G4endl
       << "phiAxis     =" << phiAxis/deg     << "(degree), " << G4endl
       << "yHalfLength1=" << yHalfLength1/cm << "(cm), "     << G4endl
       << "xHalfLength1=" << xHalfLength1/cm << "(cm), "     << G4endl
       << "xHalfLength2=" << xHalfLength2/cm << "(cm), "     << G4endl
       << "tanAlpha1   =" << tanAlpha1       << "(), "       << G4endl
       << "yHalfLength2=" << yHalfLength2/cm << "(cm), "     << G4endl
       << "xHalfLength3=" << xHalfLength3/cm << "(cm), "     << G4endl
       << "xHalfLength4=" << xHalfLength4/cm << "(cm)."      << G4endl
       << "tanAlpha2   =" << tanAlpha2       << "(), "       << G4endl
       << "(x4-x3)*y1/(x2-x1)/y2="
       << (xHalfLength4-xHalfLength3)*yHalfLength1/(xHalfLength2-xHalfLength1)/yHalfLength2 << G4endl
       << "tanAlpha2/tanAlpha1=" << tanAlpha2/tanAlpha1 <<G4endl;
  //G4cout << fStartID << "-->This is an output sentence for debug." << G4endl;
 
  trackerChamber.SetAllParameters(
                  zHalfLength,
                  thetaAxis,
                  phiAxis,
                  yHalfLength1,
                  xHalfLength1,
                  xHalfLength2,
                  tanAlpha1,
                  yHalfLength2,
                  xHalfLength3,
                  xHalfLength4,
                  tanAlpha2);
 
}

ComputeIDAndSide()

void ExtBesCrystalParameterisation::ComputeIDAndSide

(

const G4VPhysicalVolume *

pPhyVol

)

const

Definition at line 200 of file ExtBesCrystalParameterisation.cxx.

{
  //for debug
  //G4Exception("ExtBesCrystalParameterisation::ComputeIDAndSide() starting......");
  G4int* pCrystalID=const_cast<G4int*>(&fCrystalID);
  G4bool* pFlagLeft=const_cast<G4bool*>(&fFlagLeft);
  G4int* pVerboseLevel=const_cast<G4int*>(&fVerboseLevel);
  
  *pVerboseLevel=
    ExtBesEmcConstruction::GetBesEmcConstruction()->GetVerboseLevel();
 
  //*pCrystalID=pPhyVol->GetMotherPhysical()->GetCopyNo(); 
  *pCrystalID=-1;
  for(G4int i=fStartID;i<=fAllCrystals;i++)
    {
      //G4cout << fBesEmcGeometry->GetPhysiBSCCrystal(i) << ".vs."
      //     << pPhyVol << G4endl;
      if(fBesEmcGeometry->GetPhysiBSCCrystal(i)==pPhyVol)
    {
      *pCrystalID=i;
    }
    }
  if(*pCrystalID==-1)
    G4Exception("The point of PhysicCrystal error!!!!!!!!!!!");
 
  //*pCrystalID=pPhyVol->GetCopyNo();
  //if(fVerboseLevel>6)
  //G4cout<<"******ExtBesCrystalParameterisation::ComputeIDAndSide******"<<G4endl
  //  <<"point of pPhyVol =" << pPhyVol << G4endl
  //  <<"point of mother  =" << pPhyVol->GetMother() << G4endl
  //  <<"CopyNo of pPhyVol=" << pPhyVol->GetCopyNo() << G4endl
  //  <<"CopyNo of mother =" << pPhyVol->GetMother()->GetCopyNo() << G4endl
  //  <<"********************************************************"<<G4endl;
  
  if(fCrystalID>fAllCrystals/2)
    {
      *pFlagLeft=false;
      *pCrystalID=fCrystalID-fAllCrystals/2;
    }
  else
    {
      *pFlagLeft=true;
      *pCrystalID=fAllCrystals/2-fCrystalID+1;
    }
  *pCrystalID=*pCrystalID-1;
}

Referenced by ComputeDimensions(), ComputeMaterial(), ComputeSolid(), and ComputeTransformation().

ComputeMaterial()

G4Material * ExtBesCrystalParameterisation::ComputeMaterial	(	const G4int	copyNo,
		G4VPhysicalVolume *	pPhysicalVol )

Definition at line 184 of file ExtBesCrystalParameterisation.cxx.

{
  ComputeIDAndSide(pPhysicalVol);
  G4int nBSCTheta=fBesEmcGeometry->GetBSCNbTheta();
  //if(fVerboseLevel>5)
  //G4cout <<"******ExtBesCrystalParameterisation::ComputeMaterial******"<<G4endl
  //   <<"Number of all crystals=" << nBSCTheta << G4endl
  //   <<"No of current crystal =" << pPhysicalVol->GetCopyNo() << G4endl
  //   <<"*******************************************************"<<G4endl;
  if(pPhysicalVol->GetCopyNo()>2*nBSCTheta)
    return ExtBesEmcConstruction::GetBesEmcConstruction()->GetCasingMaterial();
  else
    return ExtBesEmcConstruction::GetBesEmcConstruction()->GetCrystalMaterial();
}

ComputeSolid()

G4VSolid * ExtBesCrystalParameterisation::ComputeSolid	(	const G4int	copyNo,
		G4VPhysicalVolume *	pPhysicalVol )

Definition at line 177 of file ExtBesCrystalParameterisation.cxx.

{
  ComputeIDAndSide(pPhysicalVol);
  return pPhysicalVol->GetLogicalVolume()->GetSolid();
}

ComputeTransformation()

void ExtBesCrystalParameterisation::ComputeTransformation	(	const G4int	copyNo,
		G4VPhysicalVolume *	physVol ) const

Definition at line 53 of file ExtBesCrystalParameterisation.cxx.

{
  //for debug
  //G4Exception("ExtBesCrystalParameterisation::ComputeTransformation() starting........");
  ComputeIDAndSide(physVol);
 
  if(fVerboseLevel>3)
    {
      G4cout<<"*ExtBesCrystalParameterisation::ComputeTransformation()*"<<G4endl;
      G4cout << "copyNo(transfered)=" << copyNo << G4endl
         << "copyNo(gotten)    =" << physVol->GetCopyNo() << G4endl
         << "fStartID          =" << fStartID << G4endl
         << "fCrystalID        =" << fCrystalID << G4endl;
      G4cout << "point of fBesEmcGeometry=" << fBesEmcGeometry << G4endl;
    }
  
  G4double xPosition= fBesEmcGeometry->GetXPosition(fCrystalID);
  G4double yPosition= fBesEmcGeometry->GetYPosition(fCrystalID);
  G4double zPosition= fBesEmcGeometry->GetZPosition(fCrystalID);
  //G4Exception("ExtBesCrystalParameterisation::ComputeTransformation() running........");
  //G4double xPosition=0,yPosition=0,zPosition=0;
  if(fFlagLeft) zPosition=-zPosition;
  G4ThreeVector origin(xPosition,yPosition,zPosition);
  physVol->SetTranslation(origin);
 
  //G4double thetaPosition=fBesEmcGeometry->GetThetaPosition(fCrystalID);
  //G4RotationMatrix *rotateMatrix;
  //rotateMatrix = new G4RotationMatrix();
  //if(fFlagLeft)
  //  {
  //    rotateMatrix->rotateZ(-90*deg);
  //    rotateMatrix->rotateX(-180*deg+thetaPosition);
  //  }
  //else
  //  {
  //    rotateMatrix->rotateZ(-90*deg);
  //    rotateMatrix->rotateX(-thetaPosition);
  //  }
  //G4double delta = 0*deg;
  //G4ThreeVector axis = G4ThreeVector(0,0,0);
  //rotateMatrix->getAngleAxis(delta, axis);
  if(fVerboseLevel>5)
    {
    G4cout << "The crystals are at the position:" << G4endl
       << origin/cm <<"(cm)" << G4endl;
    //G4cout << thetaPosition/deg <<"(deg)." << G4endl
    //<< "The axis of crystals in the world system is: " << G4endl
    //<< delta/deg << "(deg)(delta) " << G4endl
    //<< axis  << "(Z axis)" << G4endl;
    }
  //physVol->SetRotation(rotateMatrix);
  
}

SetVerboseLevel()

void ExtBesCrystalParameterisation::SetVerboseLevel ( G4int val )

inline

Definition at line 41 of file ExtBesCrystalParameterisation.h.

{fVerboseLevel=val;};

---

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

• ExtBesCrystalParameterisation.h
• ExtBesCrystalParameterisation.cxx