/scratch0/jsnuveri/BDSIM/BDSIMgit/bdsim/src/BDSDecapole.cc

/* BDSIM code.    Version 1.0
   Author: Grahame A. Blair, Royal Holloway, Univ. of London.
   Last modified 24.7.2002
   Copyright (c) 2002 by G.A.Blair.  ALL RIGHTS RESERVED. 

   Modified 22.03.05 by J.C.Carter, Royal Holloway, Univ. of London.
   Changed StringFromInt to be the BDSGlobal version
*/

#include "BDSGlobalConstants.hh" 

#include "BDSDecapole.hh"
#include "BDSDecMagField.hh"
#include "BDSDecStepper.hh"
#include "BDSMagnet.hh"
#include "BDSMagnetOuterInfo.hh"

#include "G4FieldManager.hh"
#include "G4LogicalVolume.hh"
#include "G4Tubs.hh"
#include "G4UserLimits.hh"
#include "G4VisAttributes.hh"
#include "G4VPhysicalVolume.hh"

class BDSTiltOffset;

BDSDecapole::BDSDecapole(G4String           name,
                         G4double           length,
                         G4double           bQuadPrime,
                         BDSBeamPipeInfo*   beamPipeInfo,
                         BDSMagnetOuterInfo magnetOuterInfo,
                         BDSTiltOffset      tiltOffset):
  BDSMagnet(BDSMagnetType::decapole, name, length,
            beamPipeInfo,magnetOuterInfo, tiltOffset),
  itsBQuadPrime(bQuadPrime)
{;}

void BDSDecapole::Build() {
  BDSMagnet::Build();
  if(BDSGlobalConstants::Instance()->GetIncludeIronMagFields())
    {
      G4double polePos[4];
      G4double Bfield[3];
      
      //coordinate in GetFieldValue
      polePos[0]=-BDSGlobalConstants::Instance()->GetMagnetPoleRadius()*sin(CLHEP::pi/10);
      polePos[1]=BDSGlobalConstants::Instance()->GetMagnetPoleRadius()*cos(CLHEP::pi/10);
      polePos[2]=0.;
      polePos[3]=-999.;//flag to use polePos rather than local track
      
      itsMagField->GetFieldValue(polePos,Bfield);
      G4double BFldIron=
        sqrt(Bfield[0]*Bfield[0]+Bfield[1]*Bfield[1])*
        BDSGlobalConstants::Instance()->GetMagnetPoleSize()/
        (BDSGlobalConstants::Instance()->GetComponentBoxSize()/2-
         BDSGlobalConstants::Instance()->GetMagnetPoleRadius());
      
      // Magnetic flux from a pole is divided in two directions
      BFldIron/=2.;
      
      BuildOuterFieldManager(10, BFldIron,CLHEP::pi/10);
    }
}

void BDSDecapole::BuildBPFieldAndStepper()
{
  // set up the magnetic field and stepper
  itsMagField=new BDSDecMagField(itsBQuadPrime);
  itsEqRhs=new G4Mag_UsualEqRhs(itsMagField);
  
  BDSDecStepper* decStepper=new BDSDecStepper(itsEqRhs);
  decStepper->SetBQuadPrime(itsBQuadPrime);
  itsStepper = decStepper;
}

---