src/Pulses/PulGauss.cc
// ------- Gaussian Pulse Length To Obtain On Resonance Rotation Angle --------
Gaussian Pulse Length To Obtain On Resonance Rotation Angle
Gaussian Pulse Length To Obtain On Resonance Rotation Angle
double Gtime(double angle, double gamB1, int N, double fact)
   {
   double den, Z;				
   Gprep(N, fact, den, Z);			
   row_vector GNs = GNvect(N, fact);		
   double normang = gamB1*Re(GNs.sum());	
   return double(N)*angle/(360.0*normang);	
   }
double Gtime ( double angle , double gamB1 , int N , double fact ) { double den , Z ; Gprep ( N , fact , den , Z ) ; row_vector GNs = GNvect ( N , fact ) ; double normang = gamB1 * Re ( GNs . sum ( ) ) ; return double ( N ) * angle / ( 360.0 * normang ) ; }
Gaussian Pulse Length To Obtain On Resonance Rotation Angle
Gaussian Pulse Length To Obtain On Resonance Rotation Angle
// Input        angle   : The pulse rotation angle (degrees) //              gamB1   : Gaussian pulse rf-field strength (Hz) //              N       : Number of Gaussian steps //		fact    : Cutoff factor [0.0,1.0] //			  defaults to 2.5 % // Output       time	: Gaussian pulse length to obtain //			  specified rotation angle on resonance //			  for indicated pulse strength & cutoff // Denominator & numerator // Prep. Gaussian values // Normalized Gaussian // Get Gaussian integral (Hz) // Scale to seconds