BesCrystalParameterisation.cc

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 "EmcSim/BesCrystalParameterisation.hh"
 
#include "EmcSim/BesEmcConstruction.hh"
#include "EmcSim/BesEmcGeometry.hh"
// #include "BesDetectorConstruction.hh"
 
#include "G4LogicalVolume.hh"
#include "G4ThreeVector.hh"
#include "G4Trap.hh"
#include "G4VPhysicalVolume.hh"
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
BesCrystalParameterisation::BesCrystalParameterisation( const G4int startID, const G4int endID,
                                                        const G4int     thetaAllCrystals,
                                                        BesEmcGeometry* besEMCGeometry,
                                                        const G4int     verboseLevel )
    : fVerboseLevel( 0 ) {
  // for debug
  // G4Exception("BesCrystalParameterisation::BesCrystalParameterisation() 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 << "BesCrystalParameterisation construction: ID of crystal>No avaible!" << G4endl;
    exit( -1 );
  }
}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
BesCrystalParameterisation::~BesCrystalParameterisation() {}
 
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
void BesCrystalParameterisation::ComputeTransformation( const G4int        copyNo,
                                                        G4VPhysicalVolume* physVol ) const {
  // for debug
  // G4Exception("BesCrystalParameterisation::ComputeTransformation() starting........");
  ComputeIDAndSide( physVol );
 
  if ( fVerboseLevel > 3 )
  {
    G4cout << "*BesCrystalParameterisation::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("BesCrystalParameterisation::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 BesCrystalParameterisation::ComputeDimensions( G4Trap& trackerChamber, const G4int copyNo,
                                                    const G4VPhysicalVolume* physVol ) const {
  // for debug
  // G4Exception("BesCrystalParameterisation::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* BesCrystalParameterisation::ComputeSolid( const G4int        copyNo,
                                                    G4VPhysicalVolume* pPhysicalVol ) {
  if ( 0 ) G4cout << copyNo << G4endl;
  ComputeIDAndSide( pPhysicalVol );
  return pPhysicalVol->GetLogicalVolume()->GetSolid();
}
 
G4Material* BesCrystalParameterisation::ComputeMaterial( const G4int         copyNo,
                                                         G4VPhysicalVolume*  pPhysicalVol,
                                                         const G4VTouchable* parentTouch ) {
  if ( 0 ) G4cout << copyNo << parentTouch << G4endl;
  ComputeIDAndSide( pPhysicalVol );
  G4int nBSCTheta = fBesEmcGeometry->GetBSCNbTheta();
  // if(fVerboseLevel>5)
  // G4cout <<"******BesCrystalParameterisation::ComputeMaterial******"<<G4endl
  //    <<"Number of all crystals=" << nBSCTheta << G4endl
  //    <<"No of current crystal =" << pPhysicalVol->GetCopyNo() << G4endl
  //    <<"*******************************************************"<<G4endl;
  if ( pPhysicalVol->GetCopyNo() > 2 * nBSCTheta )
    return BesEmcConstruction::GetBesEmcConstruction()->GetCasingMaterial();
  else return BesEmcConstruction::GetBesEmcConstruction()->GetCrystalMaterial();
}
 
void BesCrystalParameterisation::ComputeIDAndSide( const G4VPhysicalVolume* pPhyVol ) const {
  // for debug
  // G4Exception("BesCrystalParameterisation::ComputeIDAndSide() starting......");
  G4int*  pCrystalID    = const_cast<G4int*>( &fCrystalID );
  G4bool* pFlagLeft     = const_cast<G4bool*>( &fFlagLeft );
  G4int*  pVerboseLevel = const_cast<G4int*>( &fVerboseLevel );
 
  *pVerboseLevel = BesEmcConstruction::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<<"******BesCrystalParameterisation::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;
}