TransportedMonopole.cxx

Go to the documentation of this file.

#include "AtlfastEvent/TransportedMonopole.h"
#include "AtlfastEvent/ChargeService.h"
#include "AtlfastEvent/MagField.h"

#include <iostream>

namespace Atlfast {
  
  void TransportedMonopole::deflect(){
    
    // std::abs accepts doubles and floats
    using std::abs;
    using std::cout;
    using std::endl;
    
    // Get undeflected eta value
    double eta=m_eta;
    
    // Use ChargeService to get particle's charge
    double charge;
    ChargeService cs;
    try{
      charge=cs(m_particle);
    }catch(std::string errMsg){
      throw errMsg;
    }
    
    
    MagField *bfield = MagField::Instance();
    
    if( bfield->is_on() && !(charge==0.) ) {
      
      assert(std::abs(charge)>0.5);
      
      double def=0.;
      double pt = m_particle->momentum().perp();
      
      
      if (abs(eta) < 1.7) {
        
        // The equation is a little complicated - broken up
        // into 3 parts
        
        double equ1 = sqrt(0.25*pt*(exp(2*eta)+1));
        
        double equ2 = 0.5*exp(eta)*sqrt(4.8e7-(pt*pt));

        double equ3 = sqrt((equ1-equ2)*(equ1-equ2) + 4.9e7);
        
        def = 0.007*(equ3-1.);
        
      }

      // Replace member eta with deflected version
      m_eta = eta + def*charge;
      
    }
    
    return;
    
  }
  
}

---