FDsimG4CorrectorRingUpperLimit.cc

Go to the documentation of this file.

#include "FDsimG4CorrectorRingUpperLimit.hh"
#include "FDsimG4Colours.hh"


#include "G4SDManager.hh"
#include "G4IntersectionSolid.hh"
#include "G4Box.hh"
#include "G4Polycone.hh"

#include "G4Material.hh"
#include "G4MaterialTable.hh"
#include "G4Tubs.hh"
#include "G4Cons.hh"
#include "G4LogicalVolume.hh"
#include "G4ThreeVector.hh"
#include "G4VisAttributes.hh"

using namespace TelescopeSimulatorLX ;

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

  Init() ;

}

FDsimG4CorrectorRingUpperLimit::FDsimG4CorrectorRingUpperLimit(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() ;

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

  Init() ;

}

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

FDsimG4CorrectorRingUpperLimit::~FDsimG4CorrectorRingUpperLimit(){;}

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

void FDsimG4CorrectorRingUpperLimit::Init(){

  SetProfileParameters() ;
  MakeLogicalVolume2() ;

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

  G4cerr << " ================================================================  " << G4endl ;
  G4cerr << " =                     Using Upper Limit Lens Profile              " << G4endl ;

  DumpInfo() ;

}

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

void FDsimG4CorrectorRingUpperLimit::SetProfileParameters()
{

  fa1 =  5.906/1000./m ;
  fa2 =  94.52/1000./m2 ;
  fa3 = -118.8/1000./m3 ;
}

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

G4double FDsimG4CorrectorRingUpperLimit::GetRadius(G4double z){

  return 0.*z ;
                             
}

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

G4double FDsimG4CorrectorRingUpperLimit::GetSagita(G4double radius){

  static G4double fz0 = (fa1/2.*pow(fMinRadius,2) + fa2/3.*pow(fMinRadius,3) + fa3/4.*pow(fMinRadius,4)) ;

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

  return Sagita_value ;

                             
}

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