12/09/2011

latex

\renewcommand*\thesection{\thechapter.\arabic{section}}
使得section 以分数开始

chapter1
section 1.1

12/01/2011

More is different

Reason to do condensed matter physics

11/29/2011

不加数字编码

\thispagestyle{empty}

11/02/2011

SD station time requirement

delta t= t- t0
require
delta t between -1000 ns and 2000 ns
otherwise it is considered accidental

11/01/2011

Total Number of Events

HAS(2004-09/2011):  3370117

SD(2004-09/2009):2023206

10/16/2011

latex 左双引号的输入法

在键盘上左上角第二排,Esc键下方,数字键1左侧,输入两个就是左双引号了。

10/11/2011

[Latex]no title page after \maketitle

Usually in latex
\maketitle
will create a new page after the commands, for example, if you want to put a title page, you can use that command.
But if you want your article start right after the abstract instead of putting it
on another page, you can do this in the document class.

\documentclass[12pt,aps,prb,preprint,notitlepage]{revtex4}

10/09/2011

L/W的意思在neutrino 文章里


L是指trigger stations 的长度,W是trigger station 的宽度

9/21/2011

9/11/2011

Research 09/12/2011

读文件,找出那些eventid,满足大角度和多个peak.

9/07/2011

惹不起还躲不起

root tfile top close too slowly,百般折腾不知原因,后来发现不定义他就好了,全部删除,没有办法之下的办法。

8/30/2011

Get one single page from pdf documents

Create one page test.pdf from example.pdf page 17

$pdftk A=example.pdf cat A17 output test.pdf

8/29/2011

Research 08/30/2011

按照Rubtsov的条件来cut,
看结果如何。

8/26/2011

读博

黑手语:很长,但是对有志读博做研究的同学很有现实意义,原文为英文,不乐意自己看英文的,已经表示你没有那种为了做学术静下心的潜力了。

8/24/2011

Research 08/24/2011

Shower front curvature can be used as a composition-sensitive parameter

7/19/2011

Research 07/19/2011

See if you can improve the correlation accuracy in this picture.

7/15/2011

Research 07/15/2011

最近把过去一年所做的总结一下给老板看看,让老板修改修改,可作为论文的基础材料。

7/13/2011

root fit 注意定义方程的区域

TF1 *f_q_iron_Particle_Num = new TF1("f_q_iron_Particle_Num","[0]*x+[1]",18.3,19.4);
这个区域18.3到19.4非常重要,若不对无法进行拟合

7/02/2011

New thought

Put Electron and Photon Number together to replace Electron Number

S1000 on ADST

it is the signal at 1000 m, it is not S38

6/27/2011

What I did for the last year

2010/05
2010/06  curvature xmax
2010/07  Shower age curvature
2010/09   Fit Jump method
2010/10  Continues Jump Method
2010/11 Jump method as proton and iron nuclei  primary study
2010/12 Risetime studies  and Jump method fit studies
2011/01 Different Method Correlation with Muon Number
2011/02 Different Method Correlation with Muon Number
2011/03 Different Method explore
2011/04 Signal and Variation Studies
2011/05 Continues
2011/06 Find Energy different definitions
2011/07 ICRC not much research
2011/08 ICRC zenith and rise time
2011/09 zenith and rise time
2011/10 Not much
2011/11 Not much
2011/12. write that gap notes

6/26/2011

6/13/2011

发现问题

energy 似乎就是这个了,但是这两个有何不同?
distance
都不同

Ratio(Peak/Charge) used by the simulation

0.353671

6/03/2011

Difference found on average energy deposit

1 later program applied start slot and end slot 排除
2 later program does not use artificial stations
3 later program does not use any threshold cut 排除

Programs running:
prog 0435, 0436, 0445, 0446, 0455, 0456
prog 0522, 0523

5/23/2011

5/10/2011

ntp time sync ubuntu

you will get this error

gyuan@zen:~$ sudo ntpdate
10 May 14:17:31 ntpdate[14675]: no servers can be used, exiting
The solution


gyuan@zen:~$ sudo /etc/init.d/ntp stop
 * Stopping NTP server ntpd                                              [ OK ]
gyuan@zen:~$ sudo ntpdate-debian
10 May 14:15:10 ntpdate[14624]: adjust time server 91.189.94.4 offset 0.000043 sec

5/05/2011

Research 05/05/2011

1 Change Energy deposit to the standard signal calculation
2 See if there is any way to get the electron content working. and see if there is any relationship between the number of particles and energy

E=a(Ne+bNmu)

4/18/2011

Auger SD T1 trigger

1,simple threshold trigger.
coincidence of the three PMTs each above 1.75 IVem peak


2,time over threshold trigger.

at least 13 bins in 120 FADC bins of a sliding window of 3 micro seconds are required to be above a threshold of 0.2 vem in coincidence of 2 out of 3 PMTs.

4/13/2011

baryon

baryon 重子
meson 介子
lepton 轻子

4/10/2011

New Idea

Have the Edeposit as a parameter, while only count the larger than 0.1 Vem bins,
Do the integral, the iron primary trace, because has more muon peaks, thus was less cut off, so will have a larger Edeposit compare to the one caused by proton primary trace.

问题

检查 prog0406 和 prog0303 对比一下,看怎么回事。
中间一个有问题。

4/07/2011

Research 04/07/2011

1 Investigate the plot of average energy deposit and plot the fluctuation vs the energy.

2 See if there is any way to get the electron content working. and see if there is any relationship between the number of particles and energy

E=a(Ne+bNmu)

4/03/2011

Research 04/04/2011

  • Adjust the y Axis of prog0310 plots
  • Change y Axis name
  • Add Epos Number

学习才是王道

不要淹没在信息世界的大海中

3/25/2011

时间 管理

  2.学会使用时间碎片和“死时间”。如果你做了上面的时间统计,你一定发现每天有很多时间流失掉了,例如等车、排队、走路、搭车等,可以用来背单 字、打电话、温习功课等。大家喜欢用MP3播放器听音乐,但是为什么不把一些学习的内容(英语会话、生字和录下来的广播新闻等)放入你的MP3播放器呢? 在国外,现在流行podcast,就是把一些可以用听来了解的信息,录入MP3播放器,在开车回家时听。我前一阵和同事一起出差,他们都很惊讶为什么我和 他们整天在一起,但是我的电子邮件都可以及时回答?后来,他们发现,当他们在飞机上和汽车上聊天、读杂志和发呆的时候,我就把电子邮件全回了。重点是,无 论自己忙还是不忙,你要把那些可以利用时间碎片做的事先准备好,到你有空闲的时候有计划地拿出来做。

    3. 每天一大早把一天该做的事排好优先次序,按照这个次序来做,并要求自己这天把最重要的三件事做完。我感到在工作和生活中每天都有干不完的事,唯一能够做的 就是分清轻重缓急。有的年轻人会说自己“没有时间学习”,其实,换个说法就是“学习没有被排上优先级次序”。曾经有一个教学生做时间管理的老师,他上课时 带来两个大玻璃缸和一堆大小不一的石头。他做了一个实验,在其中一个玻璃缸中先把小石、砂倒进去,最后大石头就放不下了。而另一个玻璃缸中先放大石头,其 它小石和砂却可以慢慢渗入。他以此为比喻说:“时间管理就是要找到自己的优先级,若颠倒顺序,一堆琐事占满了时间,重要的事情就没有空位了。”

Standard Deviation and error of the mean

Standard Deviation is defined as



While the error of mean is defined as

3/24/2011

Android sync problem___ Sync is currently experiencing problems. it will be back shortly

Had a search around and found a solution that worked for me (p.s. posted by someone else in another forum)
Go to "Settings".  Select "Applications".
Select "Manage Applications".
Scroll down to "Contacts Storage"
Click on "Clear Data"
Go back to to "Settings" again and in "Data Synchronization"
resync by pressing "Synchronize Now"
Good Luck
 

3/23/2011

made changes to prog0218, please be careful about the numbers

    if(Num_muon>0 && Num_of_peaks>0 && Num_of_Jumps_Drops>0)
    {
    Group_Muon_Num[q]+=Num_muon;
     Group_Jump_Integral[q]+=Jump_Integral;
     Group_Jump_Drop_Integral[q]+=Jump_Drop_Integral;
     Group_Num_of_Jumps_Drops[q]+=Num_of_Jumps_Drops;
     Group_Num_of_peaks[q]+=Num_of_peaks;
     Group_Peak_integral[q]+=Peak_integral;
    j++;
    input>>Num_muon;
    }//if(Num_muon)
    else {m++;
    input>>Num_muon;
        continue;}

    if(j % Num_per_group == 0)
    {
    q++;
    output3<

c++ int operation 整数变换

c++整数变换是把小数部分去掉,而不是四舍五入。 
  • Example 1,

double b=3.4
int a=int(b);
输出
a=3;
  • Example 2,
double b=3.7
int a=int(b);
输出
a=3;

3/17/2011

第一个任务

Nmu vs Energy几种方法的图画出来

prog0217 的Method_comparison_QGSjetII.txt中包含了energy 信息,很好阿!

I am a little bit lost, please help me lord.

2/26/2011

Paring with Samsung WEP475

It's easy to pair the WEP475:
When you turn on your Samsung Bluetooth headset for the first time, it proceeds into pairing mode immediately.
If you have paired the headset with other Bluetooth devices, make sure that the headset is on. Press and hold the Talk button until the blue indicator light stays solid blue after blinking.
Set your Bluetooth phone to discover the headset by following your phone's guide. Typically, the steps involve going to a 'Setting', 'Connection', or 'Bluetooth' menu on your phone and then selecting the option to discover Bluetooth devices.
Select the headset from the list of devices found by your phone.
Enter a passkey or PIN, 0000 (4 zeros), then press the Yes or OK key.
If pairing was successful, the indicator light will flash in blue 10 times rapidly before going to flashing every 6 seconds. If unsuccessful, the light remains on, and you will need to re-attempt pairing.
When you turn on your headset after pairing with the phone, the headset is automatically connected to the phone.
When you turn on your phone after pairing with the headset, the phone is connected to the headset by pressing the Talk button

2/21/2011

加油

明天之前做出一些东西跟老板见面

第一件事,把epos和QGSjetII的文件合并,然后整体画图

2/18/2011

Research 2/18/2011

What should I do now?
1, create similar plot of the muon number vs energy on data and simulation

2/16/2011

Earth's atomosphere fact

  • 50% of the atmosphere by mass is below an altitude of 5.6 km (18,000 ft).
  • 90% of the atmosphere by mass is below an altitude of 16 km (52,000 ft). The common altitude of commercial airliners is about 10 km (33,000 ft) and Mt. Everest's summit is 8,848 m (29,029 ft) above sea level.
  • 99.99997% of the atmosphere by mass is below 100 km (62 mi; 330,000 ft), although in the rarefied region above this there are auroras and other atmospheric effects. The highest X-15 plane flight in 1963 reached an altitude of 108 km (354,300 ft).
  • 地球半径 6371 km

1/28/2011

Error 误差

delta x=2cm
delta y=2cm
  • if the two uncertainty are independent, the error of
z=x+y
will be
delta z=sqrt(2^2+2^2)=2.8 cm
the mean value of  x and y has an error of 
delta mean=2.8/2=1.4 cm
  • if you suspect the uncertainties may not be independent, the
error will be
delta z=2+2=4 cm
the mean value of  x and y has an error of 
delta mean=4/2=2 cm

1/27/2011

Creat a list of files

//define the variable
ostringstream hname,
                          htitle;
//put the number in the name
hname.str("");
        hname<< "Group_Num=" <
//use the name
c[l] = new TCanvas(hname.str().c_str(),htitle.str().c_str(),200,10,700,500);

Research 1/27/2011

Count the muons with peak searching method and also group them to see if there is correlation with the MC muon numbers

1/24/2011

install .bin file

chmod +x filename.bin
sudo ./filename.bin

1/23/2011

Research 01/23/2011

N_mu= P1/0.063

1 How well does it work in other core distance
2 in Monte Carlo , use the relationship above, plot the Fe primary and Proton primary curve of N_mu with Energy
3 in the future, counting peaks and compare them with real muons number

1/20/2011

Acer netbook windows xp can not boot

title        Windows XP Professional
root        (hd0,1)
savedefault
makeactive
chainloader    +1

The default number is 0, change it to 1 and give it a try

1/19/2011

Set Statistic box position in root

gStyle->SetStatX(0.5); 

gStyle->SetStatY(0.88); 

This is setting the middle point of the statistic box position on the graph, where the 0.5 and 0.88 was the percentage of the length and width of the canvas.

Edit the Places Menu on Ubuntu

Open the file browser, and from the menu it's Bookmarks/Edit Bookmarks

1/17/2011

Research 1/17/2010

1 Find the linear relation of NumMuons with the fitted parameter p1
2 Use Monte Carlo simulation to see what number of muon my relationship give me compare to the real number of muon.
3 Check NumMuon ~  p1 relation, to see if I can loose the restriction on the distance cut and angle cut or the number of stations needed to do jump fit.

1/14/2011

求余数 C++

英文中 余数 叫 reminder
求余数操作 叫 Modulo
Modulo is the operation that gives the remainder of a division of two values。
在c++中是百分号 %
例如:
a=11;
b=3;
c=a%b;
结果时c=2;

1/07/2011

Research 1/7/2011

Do the fit of stations by stations, then do a scatter plot of the number of muons and fitted parameter. give it a try!!!

When calculating the number of electrons, try no cut at all
Try both electron number and both electron and photon numbers.

1/06/2011

Research 1/6/2010

Try to improve jump method fitting

do a KS test to see which one works better

1 collect the jumps as we do before
2 get the muon and electron component out of it
3 try different method to see which one is working better.

1/04/2011

新的一年开始了,加油加油加油!!!

2011是很重要的一年,
将在我的一生中扮演重要的角色,
加油加油加油!!!