Changeset 4069 in Sophya for trunk


Ignore:
Timestamp:
Apr 27, 2012, 11:47:27 AM (13 years ago)
Author:
ansari
Message:

dernières corrections (proofreading) du papier avant publication par A&A, 26 Mars 2012, 27/04/2012

Location:
trunk/Cosmo/RadioBeam
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/Cosmo/RadioBeam/interfconfigs.cc

    r3947 r4069  
    1818  cout << ">>>CreateFilledSqConfig(" << nd << "," << Ddish << "," << Eta << ") ---> NDishes=" << vd.size() << endl;
    1919 
     20  return vd;
     21}
     22
     23/* --Fonction -- */
     24vector<Dish> CreateSparseSqConfig(int nd, double LSep, double Ddish, double Eta)
     25{
     26  vector<Dish>  vd;
     27  int cnt=0;
     28  for(int i=0; i<nd; i++)
     29    for(int j=0; j<nd; j++) {
     30      cnt++;
     31      vd.push_back(Dish(cnt, i*LSep, j*LSep, Eta*Ddish));
     32    }
     33  cout << ">>>CreateSparseSqConfig(" << nd << ", LSep=" << LSep
     34       << ", Ddish=" << Ddish << "," << Eta << ") ---> NDishes=" << vd.size() << endl;
    2035  return vd;
    2136}

  • trunk/Cosmo/RadioBeam/interfconfigs.h

    r3947 r4069  
    1818//  Filled square array of ndxnd dishes
    1919vector<Dish> CreateFilledSqConfig(int nd, double Ddish=5., double Eta=0.9);
     20//  Sparse regular square array of ndxnd dishes, each with diameter Ddish, separated at Lsep
     21vector<Dish> CreateSparseSqConfig(int nd, double LSep=5., double Ddish=1., double Eta=0.9);
    2022//  Semi filled square array of ndxnd dishes
    2123vector<Dish> CreateSemiFilledSqConfig(int nd, double Ddish=5., double Eta=0.9);

  • trunk/Cosmo/RadioBeam/repicon.cc

    r3948 r4069  
    4343  cout << " Usage: repicon [-parname Value] configId OutPPFName \n"
    4444       << " configIds: f4x4,f8x8,f11x11,f20x20, confA,confB,confC,confD, hex12,cross11, \n"
    45        <<  "           f4cyl,f8cyl,f4cylp,f4cylp, nan12,nan24,nan36,nan40,nan128 \n"
     45       <<  "           f4cyl,f8cyl,f4cylp,f4cylp, s4x4,s11x11,s20x20  nan12,nan24,nan36,nan40,nan128 \n"
    4646       << "   f4x4 , f8x8 , f11x11 , f20x20 Filled array of nxn dishes  \n"
     47       << "   s4x4 , s11x11 , s20x20  sparse/semi-Filled regular array of nxn dishes  \n"
    4748       << "   confA , confB, confC, confD : semi-filled array of dishes \n"
    4849       << "   hex12,cross11 : ASKAP like double hexagonal (12xD=12m), cross config (11xD=12m) \n"
    … …  
    5657       << "    -z redshift (default=0.7) --> determines Lambda \n"
    5758       << "    -D DishDiameter (default=5 m) \n"
     59       << "    -lsep Lsep (default=5 m) dish separation for sparse/semi-filled array \n"
    5860       << "    -eta fill_factor (default=0.90) \n"
    5961       << "    -lmax array extension (default=100 m ) for response calculation kmax \n"
    … …  
    8890 
    8991  double Ddish=5.;
     92  double Lsep=5.;
    9093  double Eta=0.9;
    9194  bool fgDfixed=false;
    … …  
    114117      Ddish=atof(arg[ka+1]);  fgDfixed=true;  ka+=2;
    115118    }
     119    else if (strcmp(arg[ka],"-lsep")==0) {
     120      Lsep=atof(arg[ka+1]);  ka+=2;
     121    }
    116122    else if (strcmp(arg[ka],"-eta")==0) {
    117123      Eta=atof(arg[ka+1]);  ka+=2;
    … …  
    178184      vdishes=CreateFilledSqConfig(20,Ddish, Eta);
    179185    }
     186
     187    else if (config=="s4x4") {
     188      vdishes=CreateSparseSqConfig(4, Lsep, Ddish, Eta);
     189    }
     190    else if (config=="s11x11") {
     191      vdishes=CreateSparseSqConfig(11, Lsep, Ddish, Eta);
     192    }
     193    else if (config=="s20x20") {
     194      vdishes=CreateSparseSqConfig(20, Lsep, Ddish, Eta);
     195    }
     196
    180197    else if (config=="f4cyl") {
    181198      cylW=12.;   cylRL=2*LAMBDA;
    … …  
    318335void SaveDTVecDishesH2Resp(POutPersist& po, vector<Dish>& vdishes, Four2DRespTable& mdresp)
    319336{
    320   char* names[5]={"did","posx","posy","diam","diamy"};
     337  const char* names[5]={"did","posx","posy","diam","diamy"};
    321338  NTuple ntvd(5,names,64,false);
    322339  r_4 xnt[10];

  • trunk/Cosmo/RadioBeam/sensfgnd21cm.tex

    r4050 r4069  
    219219$1 \lesssim z \lesssim 2$, and possibly up to the reionization redshift \citep{wyithe.08}.
    220220
    221 In section 2, we  discuss the intensity mapping and its potential for measuring of the
     221In section 2, we  discuss the intensity mapping and its potential for measuring the
    222222\HI mass distribution power spectrum. The method used in this paper to characterize
    223223a radio instrument response and sensitivity for $P_{\mathrm{H_I}}(k)$ is presented in section 3.
    … …  
    398398compared to its current value $\gHI(z=1.5) \sim 0.025$.
    399399The 21 cm brightness temperature and the corresponding power spectrum can be written as
    400 (\cite{madau.97}; \cite{zaldarriaga.04}); \cite{barkana.07})
     400(\cite{madau.97}; \cite{zaldarriaga.04}; \cite{barkana.07})
    401401\begin{eqnarray}
    402402 P_{T_{21}}(k) & = & \left( \bar{T}_{21}(z)  \right)^2 \, P(k)    \label{eq:pk21z} \\
    … …  
    628628{ \changemark
    629629\begin{eqnarray}
    630 \alpha , \beta & \rightarrow &  \ell_\perp = l_x, l_y = (1+z) \, \dang(z) \,  \alpha,\beta  \\
     630\alpha , \beta & \rightarrow &  \ell_\perp = \ell_x, \ell_y = (1+z) \, \dang(z) \,  \alpha,\beta  \\
    631631\uv & \rightarrow & k_\perp = k_x, k_y = 2 \pi \frac{ \uvu , \uvv }{  (1+z) \, \dang(z)  }   \label{eq:uvkxky} \\
    632632\delta \nu & \rightarrow & \delta \ell_\parallel = (1+z) \frac{c}{H(z)} \frac{\delta \nu}{\nu}
    … …  
    672672to a 3D white noise, with a uniform noise spectral density:}
    673673\begin{equation}
    674 P_{noise}(k_\perp, l_\parallel(\nu) ) = P_{noise} = 2 \, \frac{\Tsys^2}{t_{int} \, \nu_{21} } \, \frac{\lambda^2}{D^2}  \, \dang^2(z) \frac{c}{H(z)} \, (1+z)^4 
     674P_{noise}(k_\perp, \ell_\parallel(\nu) ) = P_{noise} = 2 \, \frac{\Tsys^2}{t_{int} \, \nu_{21} } \, \frac{\lambda^2}{D^2}  \, \dang^2(z) \frac{c}{H(z)} \, (1+z)^4 
    675675\label{ctepnoisek}
    676676\end{equation}
    … …  
    980980 
    981981\subsection{ Synchrotron and radio sources }
    982 We modeled the radio sky in the form of three 3D maps (data cubes) of sky temperature
     982We modeled the radio sky in the form of 3D maps (data cubes) of sky temperature
    983983brightness $T(\alpha, \delta, \nu)$ as a function of two equatorial angular coordinates $(\alpha, \delta)$
    984984and the frequency $\nu$. Unless otherwise specified, the results presented here are based on simulations of
    … …  
    15791579The \HI power spectrum is divided by an envelop curve $P(k)_{ref}$
    15801580corresponding to the power spectrum without baryonic oscillations.
    1581 The dots represents one simulation for a "packed" array of cylinders 
     1581The dots represents one simulation for a "packed" array of dishes 
    15821582with a system temperature,$T_{sys}=50$K, an observation time,
    15831583$T_{obs}=$ 1 year,
    … …  
    16091609\includegraphics[width=8.5cm]{Figs/AveragedPk.pdf}
    16101610\caption{1D projection of the power spectrum averaged over 100 simulations
    1611 of the packed cylinder array $b$.
     1611of the packed dish array.
    16121612The simulations are performed for the following conditions: a system
    16131613temperature $T_{sys}=50$K, an observation time $T_{obs}=1$ year,
    … …  
    18091809
    18101810\section{Conclusions}
    1811 The 3D mapping of redshifted 21 cm emission though {\it intensity mapping} is a novel and complementary
     1811The 3D mapping of redshifted 21 cm emission through {\it intensity mapping} is a novel and complementary
    18121812approach to optical surveys for studying the statistical properties of the LSS in the universe
    18131813up to redshifts $z \lesssim 3$. A radio instrument with a large instantaneous field of view

  • trunk/Cosmo/RadioBeam/specpk.cc

    r4030 r4069  
    421421  Histo& hmcntok=(*hmcntok_p_);
    422422  Histo fracmodok=hmcntok/hmcnt;
    423   char* nomcol[5] = {"k","pnoise","nmode","nmodok","fracmodok"};
     423  const char* nomcol[5] = {"k","pnoise","nmode","nmodok","fracmodok"};
    424424  dt.Clear();
    425425  dt.AddDoubleColumn(nomcol[0]);   
Note: See TracChangeset for help on using the changeset viewer.