MOTDen_dplyr: MOT Density Calculation for use with dplyr

In bgrich/JLstartup: Start-up Functions for Analysis

Description

MOTDen_dplyr calculates the density of Rydberg atoms in the MOT.

Usage

MOTDen_dplyr(Time, Signal, show_plot = FALSE, density_only = TRUE)

Arguments

Time	a vector. This is the the time axis (x-axis)
Signal	a vector. This is the voltage/signal axis (y-axis)
show_plot	a logical. This option produces a plot of the original data and the new gaussian fit if set to TRUE.
density_only	a logical. This option changes the output of the function. When TRUE, the function only outputs the density. When FALSE, the function outputs a list containing the density, the total number of atoms, and the half-width-half-max of the MOT in cm.

Details

This function calculates the density of Rydberg atoms in the MOT for a given MOT trace. It takes a data frame and, using trapz and fitG, computes the integral of the MOT fluorescence and it's full-width at half-maximum (FWHM). To get the total integral of the MOT, it is assumed that the MOT is a Gaussian in all dimensions. It then calculates the total number of atoms by converting the integral of the MOT fluorescence into the power and dividing that by the power per photon emitted. The number of Rydberg atoms takes the total number of atoms and divides it by 3.5. This is a compromise between cases where 1/3 of the atoms enter the Rydberg state and 1/4 of the atoms enter the Rydberg state.

The function returns the density in atoms/cm^{-1}.

This version of the MOT density code is for use with a tibble and dplyr.

Note: This is only for excitation to a single Rydberg state! For excitation to multiple states, the number of Rydberg atoms compared to the total number of atoms needs to be adjusted.

This is an R implementation of the pcamera function/script written by Mary Kutteruf for Matlab.

bgrich/JLstartup documentation built on May 12, 2019, 8:21 p.m.