ShipCompField.cxx

Go to the documentation of this file.

// SPDX-License-Identifier: LGPL-3.0-or-later
// SPDX-FileCopyrightText: Copyright CERN for the benefit of the SHiP
// Collaboration
 
#include "ShipCompField.h"
 
#include <iostream>
 
ShipCompField::ShipCompField(const std::string& label,
                             TVirtualMagField* firstField)
    : TVirtualMagField(label.c_str()), theFields_() {
  theFields_.push_back(firstField);
}
 
ShipCompField::ShipCompField(const std::string& label,
                             TVirtualMagField* firstField,
                             TVirtualMagField* secondField)
    : TVirtualMagField(label.c_str()), theFields_() {
  theFields_.push_back(firstField);
  theFields_.push_back(secondField);
}
 
ShipCompField::ShipCompField(const std::string& label,
                             const std::vector<TVirtualMagField*>& theFields)
    : TVirtualMagField(label.c_str()), theFields_(theFields) {}
 
ShipCompField::~ShipCompField() {
  // The destructor does nothing since this class does not own
  // the various TVirtualMagField pointers
}
 
void ShipCompField::Field(const Double_t* position, Double_t* B) {
  // Loop over the fields and do a simple linear superposition
 
  // First initialise the field components to zero
  B[0] = 0.0, B[1] = 0.0, B[2] = 0.0;
 
  for (auto* theField : theFields_) {
    if (theField) {
      // std::cout<<"Finding field for "<<theField->GetName()<<std::endl;
 
      // Find the magnetic field components for this part
      Double_t BVect[3] = {0.0, 0.0, 0.0};
      theField->Field(position, BVect);
 
      // Simple linear superposition of the B field components
      B[0] += BVect[0];
      B[1] += BVect[1];
      B[2] += BVect[2];
 
      // std::cout<<"B = "<<BVect[0]<<", "<<BVect[1]<<", "<<BVect[2]<<std::endl;
      // std::cout<<"BSum = "<<B[0]<<", "<<B[1]<<", "<<B[2]<<std::endl;
    }
  }
}