FDsimG4VCorrectorRing.cc

Go to the documentation of this file.

#include "FDsimG4VCorrectorRing.hh"

#include "G4Polycone.hh"

#include "G4LogicalVolume.hh"
#include "G4PVPlacement.hh"
#include "G4ThreeVector.hh"

using namespace TelescopeSimulatorLX ;

FDsimG4VCorrectorRing::FDsimG4VCorrectorRing(G4double rmin, G4double rmax, G4double minPhi, G4double deltaPhi, G4double thickness, G4Material* material, G4int npoints) 
{

  fMinRadius   =  rmin  ; 
  fMaxRadius   =  rmax  ; 

  fStartPhi    = minPhi ;
  fDeltaPhi    = deltaPhi ; 
  fThickness   = thickness  ;

  fNpoints     = npoints ;

  fMaterial    = material ;
  fLogicalVol  = NULL ;

}

FDsimG4VCorrectorRing::~FDsimG4VCorrectorRing(){

  delete [] fRminVec ;
  delete [] fRmaxVec ;
  delete [] fZVec ;

  fRminVec = NULL ; 
  fRmaxVec = NULL ; 
  fZVec = NULL ; 

}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
void FDsimG4VCorrectorRing::MakeLogicalVolume1(){


  static const  G4int NUMZPLANES  = fNpoints ; 
  static const G4double z0 = GetSagita(fMinRadius) ;

  fRminVec = new G4double[NUMZPLANES] ; 
  fRmaxVec = new G4double[NUMZPLANES] ; 
  fZVec = new G4double[NUMZPLANES] ; 

  G4double SagMin = GetSagita(fMinRadius) ;
  G4double SagMax = GetSagita(fMaxRadius) ;

  G4double Sagita ;

  for (G4int iz=0 ; iz < NUMZPLANES-1 ; iz++){
    fRmaxVec[iz] = fMaxRadius ; 
    Sagita  = SagMax -  (iz)*((SagMax-SagMin)/(NUMZPLANES-2)) ;
    fRminVec[iz] = GetRadius(Sagita) ;  
    fZVec[iz] = Sagita - z0 ; 

  }
 

  fRminVec[NUMZPLANES-1] = fMinRadius ; 
  fRmaxVec[NUMZPLANES-1] = fMaxRadius ;  
  fZVec[NUMZPLANES-1] = fZVec[NUMZPLANES-2] - fThickness ; 


  G4Polycone *ring_solid
    = new G4Polycone("CorrectorRing",fStartPhi,fDeltaPhi,NUMZPLANES,fZVec,fRminVec,fRmaxVec);

  fLogicalVol 
    = new G4LogicalVolume(ring_solid,fMaterial,"CorrectorRing",0,0,0); 


}


//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
void FDsimG4VCorrectorRing::MakeLogicalVolume2(){

  static const  G4int NUMZPLANES  = fNpoints ; 

  static const G4double z0 = GetSagita(fMinRadius) ;

  fRminVec = new G4double[NUMZPLANES] ; 
  fRmaxVec = new G4double[NUMZPLANES] ; 
  fZVec = new G4double[NUMZPLANES] ; 

  for (G4int iRad=0 ; iRad < NUMZPLANES-1 ; iRad++){
    fRmaxVec[iRad] = fMaxRadius ; 
    fRminVec[iRad] = fMaxRadius - (iRad)*((fMaxRadius-fMinRadius)/(NUMZPLANES-2)) ;
    fZVec[iRad]    = GetSagita(fRminVec[iRad]) - z0;  
  }
 

  fRminVec[NUMZPLANES-1] = fMinRadius ; 
  fRmaxVec[NUMZPLANES-1] = fMaxRadius ;  
  fZVec[NUMZPLANES-1]    = fZVec[NUMZPLANES-2]  - fThickness ; 


  G4Polycone *ring_solid
    = new G4Polycone("CorrectorRing",fStartPhi,fDeltaPhi,NUMZPLANES,fZVec,fRminVec,fRmaxVec);

  fLogicalVol 
    = new G4LogicalVolume(ring_solid,fMaterial,"CorrectorRing",0,0,0); 


}


//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void FDsimG4VCorrectorRing::DumpInfo(){

  G4cerr << " ================================================================  " << G4endl ;
  G4cerr << " =                   Corrector Lens  parameters                    " << G4endl ;
  G4cerr << " =                                                                 " << G4endl ;
  G4cerr << " =  Inner radius     :  "  << fMinRadius/mm << " (mm)"                 << G4endl; 
  G4cerr << " =  Outer radius     :  "  << fMaxRadius/mm << " (mm)"                 << G4endl; 
  G4cerr << " =                                                                 " << G4endl ;
  G4cerr << " =  Azimuthal span   :  "  << fStartPhi*180./pi << " < Phi < " << (fStartPhi+fDeltaPhi)*180.0/pi << " (deg.) " << G4endl ; 
  G4cerr << " =                                                                 " << G4endl ;
  G4cerr << " =  Thickness        :  " << fThickness/mm << " (mm) " << G4endl ;
  G4cerr << " =                                                                 " << G4endl ;
  G4cerr << " =  Material         :  " << fMaterial->GetName()                   << G4endl ;
  G4cerr << " =                                                                 " << G4endl ;    
  G4cerr << " ================================================================  " << G4endl ;

#ifdef VERBOSE_USE
  for (G4int iz=0 ; iz < 10 ; iz++){
    G4cerr << "Rmin,Rmax,z " << fRminVec[iz]/mm << " " <<  fRmaxVec[iz]/mm << " " << 
      fZVec[iz]/mm << G4endl ;
  }

  G4cerr << "." << G4endl <<  "." << G4endl <<  "." << G4endl <<  "." << G4endl ;

  for (G4int iz= (fNpoints<10) ? fNpoints : fNpoints-10  ; iz < fNpoints ; iz++){
    G4cerr << "Rmin,Rmax,z " << fRminVec[iz]/mm << " " <<  fRmaxVec[iz]/mm << " " << 
      fZVec[iz]/mm << G4endl ;
  }

#endif

}