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

#include "globals.hh" //G4 global constants & types
#include "BDSTeleporter.hh"
#include "BDSAcceleratorComponent.hh"
#include "BDSGlobalConstants.hh"
#include "G4VisAttributes.hh"
#include "G4LogicalVolume.hh"
#include "G4Box.hh" 
#include "BDSMaterials.hh"
#include "G4ChordFinder.hh"
#include "G4FieldManager.hh"
#include "BDSDebug.hh"
#include "BDSMagField.hh"
#include "G4Mag_UsualEqRhs.hh"
#include "BDSTeleporterStepper.hh"
#include "G4ThreeVector.hh"
#include "BDSBeamline.hh"
#include <cmath>
#include "parser/enums.h"

BDSTeleporter::BDSTeleporter(G4String name,
                             G4double length):
  BDSAcceleratorComponent(name, length, 0, "teleporter"),
  itsChordFinder(NULL),itsFieldManager(NULL),itsStepper(NULL),itsMagField(NULL),itsEqRhs(NULL)
{
#ifdef BDSDEBUG
  G4cout << __METHOD_NAME__ << " Constructing Teleporter of length: " 
         << length/CLHEP::m << " m" << G4endl;
#endif
}

void BDSTeleporter::Build()
{
  BuildBPFieldAndStepper();   // create custom stepper
  BuildBPFieldMgr(itsStepper,itsMagField);  // register it in a manager
  BDSAcceleratorComponent::Build(); // create container and attach stepper
}

void BDSTeleporter::BuildContainerLogicalVolume()
{
  G4double radius = BDSGlobalConstants::Instance()->GetSamplerDiameter() * 0.5;
  containerSolid = new G4Box(name+"_container_solid",
                             radius,
                             radius,
                             chordLength*0.5);
  containerLogicalVolume = new G4LogicalVolume(containerSolid,
                                               emptyMaterial,
                                               name + "_container_lv");
  containerLogicalVolume->SetFieldManager(itsFieldManager,false); // modelled from BDSMagnet.cc

  // register extents with BDSGeometryComponent base class
  SetExtentX(-radius,radius);
  SetExtentY(-radius,radius);
  SetExtentZ(-chordLength*0.5, chordLength*0.5);
}
  
void BDSTeleporter::BuildBPFieldAndStepper()
{
#ifdef BDSDEBUG
  G4cout << "BDSTeleporter Build Stepper & Field " << G4endl;
#endif
  // set up the magnetic field and stepper
  itsMagField = new BDSMagField(); //Zero magnetic field.
  itsEqRhs    = new G4Mag_UsualEqRhs(itsMagField);
  itsStepper  = new BDSTeleporterStepper(itsEqRhs);
}

void BDSTeleporter::BuildBPFieldMgr( G4MagIntegratorStepper* stepper,
    G4MagneticField* field)
{
  //this is all copied from BDSMagnet.cc although names tidied a bit
  itsChordFinder = 
    new G4ChordFinder(field,
    chordLength*0.5/CLHEP::m,
                      stepper);

  itsChordFinder->SetDeltaChord(BDSGlobalConstants::Instance()->GetDeltaChord());
  itsFieldManager = new G4FieldManager();
  itsFieldManager->SetDetectorField(field);
  itsFieldManager->SetChordFinder(itsChordFinder);
  if(BDSGlobalConstants::Instance()->GetDeltaIntersection()>0){
    itsFieldManager->SetDeltaIntersection(BDSGlobalConstants::Instance()->GetDeltaIntersection());
  }
  if(BDSGlobalConstants::Instance()->GetMinimumEpsilonStep()>0)
    itsFieldManager->SetMinimumEpsilonStep(BDSGlobalConstants::Instance()->GetMinimumEpsilonStep());
  if(BDSGlobalConstants::Instance()->GetMaximumEpsilonStep()>0)
    itsFieldManager->SetMaximumEpsilonStep(BDSGlobalConstants::Instance()->GetMaximumEpsilonStep());
  if(BDSGlobalConstants::Instance()->GetDeltaOneStep()>0)
    itsFieldManager->SetDeltaOneStep(BDSGlobalConstants::Instance()->GetDeltaOneStep());
}

void BDS::CalculateAndSetTeleporterDelta(BDSBeamline* thebeamline)
{
  // get position of last item in beamline
  // and then calculate necessary offset teleporter should apply
  G4ThreeVector lastitemposition   = thebeamline->back()->GetReferencePositionEnd();
  G4ThreeVector firstitemposition  = thebeamline->front()->GetReferencePositionStart();
  G4ThreeVector  delta             = lastitemposition - firstitemposition;
#ifdef BDSDEBUG
  G4cout << "Calculating Teleporter delta" << G4endl;
  G4cout << "last item position  : " << lastitemposition  << " mm" << G4endl;
  G4cout << "first item position : " << firstitemposition << " mm" << G4endl;
#endif
  G4cout << "Teleport delta      : " << delta << " mm" << G4endl;
  BDSGlobalConstants::Instance()->SetTeleporterDelta(delta);
  
  // calculate length of teleporter
  // beamline is built along z and sbend deflects in x
  // setting length here ensures that length is always the z difference
  G4double teleporterlength       = fabs(delta.z()) - 1e-8;
  G4cout << "Calculated teleporter length : " << teleporterlength << " mm" << G4endl;
  BDSGlobalConstants::Instance()->SetTeleporterLength(teleporterlength);
}

BDSTeleporter::~BDSTeleporter()
{
  delete itsMagField;
  delete itsEqRhs;
  delete itsStepper;
}

---