ExtBesCrystalParameterisation.cxx

Go to the documentation of this file.

//---------------------------------------------------------------------------//
//      BOOST --- BESIII Object_Oreiented Simulation Tool                    //
//---------------------------------------------------------------------------//
//Descpirtion: EMC detector 
//Author: Fu Chengdong
//Created: Sep 4, 2003
//Comment:
//---------------------------------------------------------------------------//
//
#include "TrkExtAlg/ExtBesCrystalParameterisation.h"
 
#include "TrkExtAlg/ExtBesEmcConstruction.h"
#include "TrkExtAlg/ExtBesEmcGeometry.h"
//#include "BesDetectorConstruction.hh"
 
#include "G4VPhysicalVolume.hh"
#include "G4LogicalVolume.hh"
#include "G4ThreeVector.hh"
#include "G4Trap.hh"
#include "G4SystemOfUnits.hh"
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
ExtBesCrystalParameterisation::ExtBesCrystalParameterisation(  
        const G4int    startID,
    const G4int    endID,
    const G4int    thetaAllCrystals,
    ExtBesEmcGeometry* besEMCGeometry,
    const G4int    verboseLevel
        ):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)
    {
      G4cout<<"ExtBesCrystalParameterisation construction: ID of crystal>No avaible!"<<G4endl; exit(-1);
    }
 
}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
ExtBesCrystalParameterisation::~ExtBesCrystalParameterisation()
{}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
void ExtBesCrystalParameterisation::ComputeTransformation
(const G4int copyNo, G4VPhysicalVolume* physVol) const
{
  //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);
  
}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
void ExtBesCrystalParameterisation::ComputeDimensions
(G4Trap& trackerChamber, const G4int copyNo,
 const G4VPhysicalVolume* physVol) const
{
  //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);
 
}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
G4VSolid* ExtBesCrystalParameterisation::ComputeSolid(const G4int copyNo,
                                       G4VPhysicalVolume *pPhysicalVol) 
{
  ComputeIDAndSide(pPhysicalVol);
  return pPhysicalVol->GetLogicalVolume()->GetSolid();
}
 
G4Material* ExtBesCrystalParameterisation::ComputeMaterial(const G4int copyNo,
                                       G4VPhysicalVolume *pPhysicalVol) 
{
  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();
}
 
void ExtBesCrystalParameterisation::ComputeIDAndSide(
                    const G4VPhysicalVolume* pPhyVol)const
{
  //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)
  {
    G4cout<<"The point of PhysicCrystal error!!!!!!!!!!!"<<G4endl;
    exit(-1);
  }
 
  //*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;
}