HardSphereYukawaFormFactor.h

// -*- C++ -*-
#ifndef ThePEG_HardSphereYukawaFormFactor_H
#define ThePEG_HardSphereYukawaFormFactor_H
//
// This is the declaration of the HardSphereYukawaFormFactor class.
//
 
#include "NucleonFormFactor.h"
#include "ThePEG/PDT/ParticleData.h"
 
namespace ThePEG {
 
class HardSphereYukawaFormFactor: public NucleonFormFactor {
 
public:
 
  HardSphereYukawaFormFactor() : yukawaRange_(0.7*femtometer), aFact_(1.2*femtometer)
  {}
 
  Complex formFactor(tcPDPtr nucleon, Energy2 q2) const {
    pair<int,int> AZ = massAndZ(nucleon->id());
    Length R = aFact_*pow(AZ.first,1./3.);
    double Rq = sqrt(q2)*R/hbarc;
    return 3./pow(Rq,3)*(sin(Rq)-Rq*cos(Rq))/(1.+sqr(yukawaRange_/hbarc)*q2);
  }
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  pair<int,int> massAndZ(long pid) const {
    pair<int,int> output;
    output.first  = (abs(pid)%10000)/10;
    output.second = (abs(pid)%10000000)/10000;
    return output;
  }
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
private:
 
  HardSphereYukawaFormFactor & operator=(const HardSphereYukawaFormFactor &) = delete;
 
private :
 
  Length yukawaRange_;
 
  Length aFact_;
 
  
};
 
}
 
#endif /* ThePEG_HardSphereYukawaFormFactor_H */