FDsimG4CorrectorRing.cc

Go to the documentation of this file.

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

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

  Init() ;

}

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

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

FDsimG4CorrectorRing::~FDsimG4CorrectorRing(){;}

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

void FDsimG4CorrectorRing::Init(){

  SetProfileParameters() ;
  MakeLogicalVolume1() ;

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

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

  DumpInfo() ;

}

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

void FDsimG4CorrectorRing::SetProfileParameters()
{

  //  Parameters from GAP-2000-003

  fR = 8383.0*mm ;
  fy0 = 795.27*mm ;
  fz0 = 8380.9*mm ;
}

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

G4double FDsimG4CorrectorRing::GetRadius(G4double z){

  G4double Radius_value = fy0 + std::sqrt(std::pow(fR,2)-pow((fz0-z),2)) ;

  return Radius_value ;
                             
}

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

G4double FDsimG4CorrectorRing::GetSagita(G4double radius){

  G4double Sagita_value = fz0 - std::sqrt(std::pow(fR,2)-pow((radius-fy0),2)) ;

  return Sagita_value ;

                             
}

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