FDsimG4CorrectorRingLowerLimit.cc

Go to the documentation of this file.

#include "FDsimG4CorrectorRingLowerLimit.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 ;

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

  Init() ;

}

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

}
FDsimG4CorrectorRingLowerLimit::FDsimG4CorrectorRingLowerLimit(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......

FDsimG4CorrectorRingLowerLimit::~FDsimG4CorrectorRingLowerLimit(){;}

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

void FDsimG4CorrectorRingLowerLimit::Init(){

  SetProfileParameters() ;
  MakeLogicalVolume2() ;

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

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

  DumpInfo() ;

}

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

void FDsimG4CorrectorRingLowerLimit::SetProfileParameters()
{

  fa1 =  -22.25/1000./m ;
  fa2 =   85.20/1000./m2 ;
  fa3 =  -87.51/1000./m3 ;
}

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

G4double FDsimG4CorrectorRingLowerLimit::GetRadius(G4double z){

  return 0.*z ;
                             
}

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

G4double FDsimG4CorrectorRingLowerLimit::GetSagita(G4double radius){

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

  return Sagita_value ;

                             
}

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