FDsimG4CorrectorRingKG.cc

Go to the documentation of this file.

#include "FDsimG4CorrectorRingKG.hh"
#include "FDsimG4LensSD.hh"
#include "FDsimG4Colours.hh"


#include "G4SDManager.hh"
#include "G4Polycone.hh"

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

using namespace TelescopeSimulatorLX ;

FDsimG4CorrectorRingKG::FDsimG4CorrectorRingKG(G4double rmin, G4double rmax, G4double startPhi, G4double deltaPhi, G4double thick, G4Material* mat, G4int Npoints) : FDsimG4VCorrectorRing(rmin,rmax,startPhi,deltaPhi,thick,mat,Npoints){

  Init() ;

}

FDsimG4CorrectorRingKG::FDsimG4CorrectorRingKG(G4double startPhi, G4double deltaPhi, G4double thick, G4Material* mat, G4int Npoints) : 
  FDsimG4VCorrectorRing(1700.0*mm/2.0,2200.0*mm/2.,startPhi,deltaPhi,thick,mat,Npoints){

  Init() ;

}
FDsimG4CorrectorRingKG::FDsimG4CorrectorRingKG(G4double thick, G4Material* mat, G4int Npoints) : 
  FDsimG4VCorrectorRing(1700.0*mm/2.0,2200.0*mm/2.,0.,2.0*pi,thick,mat,Npoints){

  Init() ;

}

FDsimG4CorrectorRingKG::~FDsimG4CorrectorRingKG(){;}

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

void FDsimG4CorrectorRingKG::Init(){

  SetProfileParameters() ;
  MakeLogicalVolume1() ;

  fLogicalVol->SetVisAttributes(new G4VisAttributes(cyan)) ;

  G4cerr << " ================================================================  " << G4endl ;
  G4cerr << " =                        Using KG Lens Profile                    " << G4endl ;

  DumpInfo() ;

}

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

void FDsimG4CorrectorRingKG::SetProfileParameters()
{

  G4double f = 3.4*m - 1.743*m;
  G4double Rd = .85*m;
  G4double A = 3./2. * Rd*Rd;
  G4double n = 1.5;
  
  G4double denominator = 32.*(n-1.)*f*f*f;
  fa1 = 1./denominator; // r^4
  fa2 = A/denominator;  // r^2
  
  fz0 = 10.*mm;

}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

G4double FDsimG4CorrectorRingKG::GetRadius(G4double z){

  G4double Radius_value=sqrt(2*(sqrt(4*fa1*(z-fz0)+fa2*fa2)+fa2))/(2*sqrt(fa1));

  return Radius_value ;
                             
}

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

G4double FDsimG4CorrectorRingKG::GetSagita(G4double radius){

  G4double Sagita_value = fz0 + fa1 * pow(radius,4) - fa2 * pow(radius,2) ;


  return Sagita_value ;
                             
}

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