src/Pulses/PulGauss.cc
// ------ Gaussian Pulse Strength To Obtain On Resonance Rotation Angle -------
Gaussian Pulse Strength To Obtain On Resonance Rotation Angle
Gaussian Pulse Strength To Obtain On Resonance Rotation Angle
double GgamB1(double angle, double tau, int N, double fact)
   {
   double den, Z;				
   Gprep(N, fact, den, Z);			
   row_vector GNs = GNvect(N, fact);		
   double normang = Re(GNs.sum())*tau/double(N);
   return (angle/360.0)/normang;		
   }
double GgamB1 ( double angle , double tau , int N , double fact ) { double den , Z ; Gprep ( N , fact , den , Z ) ; row_vector GNs = GNvect ( N , fact ) ; double normang = Re ( GNs . sum ( ) ) * tau / double ( N ) ; return ( angle / 360.0 ) / normang ; }
Gaussian Pulse Strength To Obtain On Resonance Rotation Angle
Gaussian Pulse Strength To Obtain On Resonance Rotation Angle
// Input        angle   : The pulse rotation angle (degrees) //              tau     : Gaussian pulse length (sec) //              N       : Number of Gaussian steps //		fact    : Cutoff factor [0.0,1.0] //			  defaults to 2.5 % //              fact    : Cutoff factor // Output       gamB1   : Gaussian pulse stength to obtain //			  specified rotation angle on resonance //			  for indicated pulse length and cutoff // Denominator & numerator // Prep. Gaussian values // Normalized Gaussian // Normalized rot. angle // Return gamB1 in Hz