Skip to content

updated CMB magnitude velocity and treatment of stars velocity disper… #1

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
fcneirad wants to merge 1 commit into
base: master
Choose a base branch
from
Open

updated CMB magnitude velocity and treatment of stars velocity disper… #1

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion include/velocities.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -39,7 +39,7 @@ namespace gerlumph {
const double l0_deg = 122.9320; // in deg
const double two_pi = 2.0*3.14159265358979323846;

const double v_apex = 387; // km/s
const double v_apex = 371; // km/s, updated magnitude from Fixsen et al. 1996
// We move towards (l,b)=(264.4,48.4) (Kogut et al. 1993), so everything is moving towards the opposite direction
const double l_dip_deg = 264.4 - 180; // in deg
const double b_dip_deg = -48.4; // in deg
Expand Down
8 changes: 4 additions & 4 deletions src/normal_source/velocities.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -42,12 +42,12 @@ void velocityComponents::createVelocitiesK04(int seed,double ra,double dec,doubl
this->pec[i].phi = getUniform(-180,180);

this->disp[i].v = velDisp(sigma_disp,epsilon,z_l,D_s,D_l);
this->disp[i].phi = getUniform(-180,180);
//this->disp[i].phi = getUniform(-180,180); // no longer needed


double vtot_x = this->pec[i].v*cos(this->pec[i].phi*this->d2r) + this->disp[i].v*cos(this->disp[i].phi*this->d2r) + this->cmb[i].v*cos(this->cmb[i].phi*this->d2r);
double vtot_y = this->pec[i].v*sin(this->pec[i].phi*this->d2r) + this->disp[i].v*sin(this->disp[i].phi*this->d2r) + this->cmb[i].v*sin(this->cmb[i].phi*this->d2r);
this->tot[i].v = sqrt( pow(vtot_x,2) + pow(vtot_y,2) );
double vtot_x = this->pec[i].v*cos(this->pec[i].phi*this->d2r) + this->cmb[i].v*cos(this->cmb[i].phi*this->d2r);
double vtot_y = this->pec[i].v*sin(this->pec[i].phi*this->d2r) + this->cmb[i].v*sin(this->cmb[i].phi*this->d2r);
this->tot[i].v = sqrt( pow(vtot_x,2) + pow(vtot_y,2) ) + this->disp[i].v; // the vel disp of the stars adds directly to the magnitude of the effective transverse velocity
this->tot[i].phi = atan2(vtot_y,vtot_x)*this->r2d;
}
}
Expand Down