12/25/2010

windows 7 自带刻录机刻录安装盘

1 下载iso 文件
2  右键单击文件,选择打开方式为
   windows 7 光盘映像刻录机

3 选择刻录后验证光盘
4 刻录

12/20/2010

Statistic box display or hide in root

One can delete the statistics box for a histogram TH1* h with:
      h->SetStats(0)
and activate it again with:
      h->SetStats(1).

reference:
http://root.cern.ch/root/html/TPaveStats.html

but do not forget to include the library


#include "TPaveStats.h"

12/19/2010

Research 12/19/2010

  1.  Do the vem trace integration
  2. Apply the energy events weighting of the QGSjetII events

12/17/2010

Previous Work Summary month by month for reference later

  • 2010
  1. Jan       To identify the muon peaks in the simulated ADST files
  2. Feb      Energy histogram of surface detector
  3. Mar     Energy histogram change with years
  4. Apr      Different method for jumps on simulation files
  5. May     Jump method fitting
  6. Jun       Hybrid Curvature Xmax relations
  7. Jul        Hybrid Curvature Shower age correlations
  8. Aug      Plot muon and EM traces separately
  9. Sep      Calculate muon/em ratio
  10. Oct      Calculate average muon jumps and compare data and simulation
  11. Nov     Calculate average muon jumps and compare data and simulation
  12. Dec     Calculate average muon jumps and compare data and simulation

  • 2009
  1. Jan
  2. Feb
  3. Mar
  4. Apr
  5. May
  6. Jun
  7. Jul
  8. Aug
  9. Sep
  10. Oct
  11. Nov
  12. Dec

Acer 5740 touchpad not working under ubuntu 10.04

Got solution from this link:
http://ubuntuforums.org/showthread.php?t=1482801
Open the terminal and type:
$gconf-editor
Then look into
/desktop/gnome/peripherals/touchpad
if touchpad_enabled is set to true. If not set it to true.

12/14/2010

Delete old emails in gmail

  1. First you need to look at the dates of those email, and then do a search based on date. Use this search string in the search box:

    before:2008-01-01

    and press Enter. 
  2. After that, click on "Select: All" on the blue area above the first email in the list.
  3. Then click on "All 20 conversations on this page are selected. Select all conversations that match this search".  at the top of the list of emails  Now press Delete.
  4.  
     Official Gmail Cheat sheet

12/13/2010

research 12/14/2010

editting prog1212, get the ratio.txt first

12/06/2010

远程查看网络状态

$bwm-ng
 or
 $ sudo iftop -p -i eth1

download multiple files from command line

got a lot help from this webpage
 1 creat a text file "download"
$ gedit download

2  put the download links in that file
3 download them from terminal
$ wget -i download

11/29/2010

Rename files in linux

for file in `ls $DATA/*.dat`

do
   mv "$file" "`basename $file .dat`.QGSjetII.dat"
done

where  basename command removing directories and suffix from file names

for example
 basename include/stdio.h .h
              Output "stdio".

11/18/2010

set environment variable in shell scripts

  1. put export commands in example.sh
  2. source example.sh

11/15/2010

Uncertainty calculations

1 adding or subtracting 
a=10+/-1
b=20+/-2
c=a+b=30+/-3
 d=b-a=10+/-3
2 multiplying or dividing

a=10+/-1
b=20+/-2
c=a*b=200+/-(1/10+2/20)*200=200+/-40
d=b/a=2+/-(1/10+2/20)*2=2+/-0.4

11/08/2010

qsub submit a job to jekyll

qsub basesub.job
文件内容(basesub.job)
#PBS -l walltime=01:00:00
#PBS -l nodes=1:ppn=1
#PBS -N PROCESSNAME
#PBS -S /bin/bash
#PBS -j oe

WD=/home/guofeng/Desktop/ADST_backup/prog2010/Nov/prog1101_jumps_simulation/

cd ${WD}
./Process_all_data.sh

11/05/2010

Root graph drawing options

gr->Draw("ACP");
gr->Draw("AB");

Google image index, looks cool

http://www.google.com/logos/index.html

Research 11/05/2010

Get average number of muons from two primaries, to see if there is difference

11/02/2010

Get mean and rms out of histogram

Assuming TH1* h, you can do
  h->GetMean()
  h->GetRMS()
  h->GetEntries()

10/20/2010

10/15/2010

Research 10/15/2010

1 Look at what condition is needed for a photon to be considered to be a photon, also for electron and muon.

The minimum energy required is 5 MeV.

2 Separate the muon and em wave to find good functions of it.

10/08/2010

Lyon shower location

ssh username@ccali.in2p3.fr
a massive amount of complete simulation files under

 /sps/pauger/Simulations/libraries/augermgr.ccin2p3

9/29/2010

PowerPoint 2007: Delete all slide notes

  1. Click the Office button, point to Prepare, and then click Inspect Document.
  2. Click Inspect.
  3. Click Remove All next to Presentation Notes. You can also choose to remove invisible on-slide or off-slide content, personal information, comments, and more.
  4. Save this version of the presentation separately (be sure not to overwrite your full version), and share it safely.

From:
CTRL+ALT+GEEK: PowerPoint 2007: Delete all slide notes:

9/27/2010

Research 09/27/2010

1, Try to fit the signal without the 0.1~0.2 portion
2, for the muon portion, fix p2 to 1, see what happens.

9/24/2010

Research 09/24/2010

1 see the jump distribution of muon and electron(completed)
2 count the photon signal and photon numbers

8/27/2010

e.g./i.e.

When you mean “for example,” use e.g. It is an abbreviation for the Latin phrase exempli gratia. When you mean “that is,” use “i.e.” It is an abbreviation for the Latin phrase id est. Either can be used to clarify a preceding statement, the first by example, the second by restating the idea more clearly or expanding upon it. Because these uses are so similar, the two abbreviations are easily confused. If you just stick with good old English “for example” and “that is” you won’t give anyone a chance to sneer at you. If you insist on using the abbreviation, perhaps “example given” will remind you to use “e.g.,” while “in effect” suggests “I.E.”

Since e.g. indicates a partial list, it is redundant(多余的) to add “etc.” at the end of a list introduced by this abbreviation.

8/23/2010

\bibitem

\bibitem[label]{cite_key}

when needed, \cite

8/20/2010

latex equation super or subscript size

If you wanted to keep the size consistent, you could declare each fraction to use the display style instead.

\begin{equation}
  x = a_0 + \frac{1}{\displaystyle a_1 
          + \frac{1}{\displaystyle a_2 
          + \frac{1}{\displaystyle a_3 + a_4}}}
\end{equation}
 

Latex pragraph break and line break

"\par" is paragraph break which have same function as a "blank line"

 "\\" is a line break

8/11/2010

Mean free path calculation

http://en.wikipedia.org/wiki/Mean_free_path
it is about the reciprocal 倒数 of the product of the particle density and  cross section.

8/09/2010

General Exam Written Part

Totally about 30 pages,
about 15 pages background
about 15 pages my own staff

Some Units(eV)

  • 1,000 Gev = 1 Tev (Tera-electron-volt)
  • 1,000,000 Gev = 1 Pev (Peta-electron volt)
  • 1,000,000,000 Gev = 1 Eev (Exa-electron-volt)

7/14/2010

Editing Menu list with Ubuntu 10.04

/boot/grub/grub.cfg

have already replaced menu.lst

6/29/2010

Segmentation fault

When you see that, first check to see if you define the canvas before you define the graph.

6/25/2010

Try to focus on your work

World Cup is nice to watch but not during your work!!!

Focus!

6/09/2010

Tcanvas can only start from one

If you define  Tcanvas *cas[0], there will be error pump out,

This is the entire stack trace of all threads:
what a shame for root!!! since it is based on c++ for a long time and c++ is always start an array at zero!

6/01/2010

Import text files to spreed sheet in OpenOffice

Need to change the file extension to ".csv", then it just rolls!!!

5/27/2010

Offline Examples

/nfs/augersoft/ApeInstalled/offline/2.6.2-AsterixObelix/share/auger-offline/doc

5/25/2010

5/17/2010

Errors when define TGraph before TCanvas

Remember that.
Always define TCanvas before TGraph.

5/14/2010

5/11/2010

How do I insert a new line in wikipedia?

http://en.wikipedia.org/wiki/Wikipedia:FAQ/Editing#How_do_I_insert_a_new_line.3F

Try a different aproach

The way to use em and muon particles have been failed,
may be I should try the area of muon and total signal,
give it a try!

5/09/2010

shear wave splitting_matlab_code


% animation.m
%

%
% Author : Pieter Coppens
% Date : 8-11-2008
% Modified by:Guofeng Yuan
% Date:5-08-2010
% close all windows and clear all variables



clear;
clc;
sr=0.2;% sampling time interval
fr=1/sr; %sampling rate
max=15; %total time
split=5; %the time when the wave splits
fs = 100;
t_end =1; % end of the time interval
t = 0:1/fs:t_end;
 j = length(t);

  z1=zeros(1,101);
  z2=z1;
%creat the vedio file
mov = avifile('shear_wave_splitting.avi','quality',50,'compression','None','fps',5)

 for i=0:sr:max
  
 if i<=split   %before the splits
    x0=t+i;
    y0= sin(2*pi*t)*cos(pi/4);
    z0=sin(2*pi*t)*sin(pi/4);
 
    plot3(x0(1:j),y0(1:j),z0(1:j),'-r','LineWidth',1.5)
     title(['Shear wave splitting']) % title
     xlabel(['Time = ',num2str(i,'%.2f'),'s']) % the x-axis shows the time-instant in seconds.
    hold off
  
    axis([0 25 -1 1 -1 1 -1 1 ])
   grid on
   axis on

 pause(0.1)

  
 else  %after the splits
    x1=t+i; %the slow wave
    y1=0.5*sin(2*pi*t);
  
    x2=t+1.8*(i-split)+split; %the fast wave
    y2=0.87*sin(2*pi*t);
 
    plot3(x1(1:j),y1(1:j),z1(1:j),x2(1:j),z2(1:j),y2(1:j),'-m','LineWidth',1.2)
    hold off
     title(['Shear wave splitting']) % title
     xlabel(['Time = ',num2str(i,'%.2f'),'s']) % the x-axis shows the time-instant in seconds.
  
    axis([0 25 -1 1 -1 1 -1 1 ])
   grid on
   axis on

pause(0.1)
 
 end
  F = getframe(gcf); % save the current figure
    mov = addframe(mov,F); % add it as the next frame of the movie
 end

 % close the open windows and clear all variables
mov = close(mov);
close all
clear all

% Open the .avi-file in the external .avi-movie application of windows.
% Comment this if you're not running this file under a Windows environment.
winopen('shear_wave_splitting.avi') %open the movie

5/05/2010

科学家名字

Nicolaus Copernicus 哥白尼
Coriolis effect 科里奥利效应

4/28/2010

Grep ignore case

  1. By default, grep is case-sensitive (use -i to ignore case).
  2. By default, grep ignores the context of a string (use -w to match words only).
  3. By default, grep shows the lines that match (use -v to show those that don't match).

TVector3 is very useful when calculating the distance

TVector3 coreTankVec   = stationPos - showerCore ;
    TVector3 distanceVector=coreTankVec.Cross(showerAxis);
    test_distance=distanceVector.Mag();

Two Important indications about shower primaries

Xmax:
As a shower induced by an heavy nuclei develops earlier and faster, the expected value of Xmax is smaller than in the case of a proton.

Muonic Fraction:
The muonic fraction is expected to be larger in the case of an heavy nuclei.

4/23/2010

缩写

WCD water cherenkov detectors
FADC Flash Analog-to-Digital Converters

4/19/2010

Research Focus

Get some thing proportional to muon/electron

Overflow and Underflow display in root

gStyle->SetOptStat(111111);
reference:
http://root.cern.ch/root/html/TStyle.html

but do not forget to include the library

#include "TStyle.h"

4/16/2010

Combining ps files and pdf files

  •  Ps files
To merge the files `slide1.ps', `slide2.ps', and `slide3.ps' into a new PostScript file, `slideshow.ps', type:
$ psmerge -oslideshow.ps slide1.ps slide2.ps slide3.ps
 
to install psmerge, in ubuntu type:
 
 $ sudo apt-get install psutils
 
  • Pdf files 
First of all you will need the program called pdftk. You can get it by opening terminal and typing:

$ sudo apt-get install pdftk

Now you can merge pdf files by putting them in the same directory and typing:

$ pdftk 1.pdf 2.pdf 3.pdf cat output merged.pdf

or if you wish to merge files in alphabetical order you can simply type:

$ pdftk *.pdf cat output merged.pdf

Make Scatter Plot, TGraph is better than TH2F

TGraph *gr1 = new TGraph(10);
   gr1->SetFillColor(6);
   gr1->SetFillStyle(3005);
   gr1->SetLineColor(6);
   gr1->SetLineWidth(603);
   gr1->SetPoint(0,140,0.5);
   gr1->SetPoint(1,130,2.9);
   gr1->SetPoint(2,124.677,3.83726);
   gr1->SetPoint(3,113.362,6.06903);
   gr1->SetPoint(4,108.513,8.00221);
   gr1->SetPoint(5,111.746,10.0272);
   gr1->SetPoint(6,119.828,12.8419);
   gr1->SetPoint(7,135.991,30.0872);
   gr1->SetPoint(8,140,40);
   gr1->SetPoint(9,135,60);
   gr1->Draw("C");
   TLatex *tex = new TLatex(140.841,37.9762,
      "#leftarrow t #rightarrow bH^{+}, H^{+} #rightarrow #tau#nu");
   tex->SetTextColor(6);
   tex->Draw();

4/11/2010

Use FindRoot in Mathematica

Sometimes Solve or NSolve in Mathematica does not work in Complicated
Equations, for this case, FindRoot is a good tool.

http://reference.wolfram.com/mathematica/ref/FindRoot.html

4/09/2010

Set X axis label

TH1I* particleDistr= new TH1I("particleDistr", "", 13, 1, 14);
  particleDistr->GetXaxis()->SetLabelSize(0.08);
 
  particleDistr->GetXaxis()->SetBinLabel(1,"#gamma");
  particleDistr->GetXaxis()->SetBinLabel(2,"e^{+}");
  particleDistr->GetXaxis()->SetBinLabel(3,"e^{-}");
  particleDistr->GetXaxis()->SetBinLabel(4,"#mu^{+}");
  particleDistr->GetXaxis()->SetBinLabel(5,"#mu^{-}");
  particleDistr->GetXaxis()->SetBinLabel(6,"#tau^{+}");
  particleDistr->GetXaxis()->SetBinLabel(7,"#tau^{-}");
  particleDistr->GetXaxis()->SetBinLabel(8,"#pi^{0}");
  particleDistr->GetXaxis()->SetBinLabel(9,"#pi^{+}");
  particleDistr->GetXaxis()->SetBinLabel(10,"#pi^{-}");
  particleDistr->GetXaxis()->SetBinLabel(11,"n");
  particleDistr->GetXaxis()->SetBinLabel(12,"p^{#pm}");
  particleDistr->GetXaxis()->SetBinLabel(13,"?");     
 
  fTracesObjects->Add(particleDistr);

3/26/2010

update unbutu from command line

sudo apt-get update &&  sudo apt-get upgrade

What I learned about the shower profile

A muonic component,which arrives early and is weakly spread in time ,as the muons propagate in straight line.

An electromagnetic component is more spread in time as its propagation is more diffusive.  

All airshower particles have a velocity of c however the muons travel in straight lines once created and therefore arrives first with a narrower time spread compared tothe electromagnetic component. 

In the WCD(water cherenkov detector),the electron and gammas are completely converted while the vertical muons deposit approximately 240 MeV of energy.

3/25/2010

New focus on my general

  1. Understand the physics of shower curvature
  2. Try to reproduce the jump method

3/24/2010

一些概念

AGN:
Active Galactic Nucleus. The core of an active galaxy responsible for the enormous luminosity. Many contains massive black holes.

Cerenkov Light:
Light emitted by a medium when traversed by a charged subatomic particle moving with a velocity greater than the speed of light in that medium. The light is emitted in the direction of motion of the inducing particle.

Extensive Air Shower
 A cascade of subatomic particles in the earth's atmosphere, initiated when an energetic cosmic ray collides with an atmospheric atom or molecule.

Flourescent Light
Light emitted through the excitation and ionization of atoms and molecule by high-speed charged subatomic particles. This light is emitted in all directions,unlike the beamed Cerenkov light.

Ionization
The process by which an tom loses electrons. May occur when a charged high-speed subatomic particle passes through a collection of atoms.

The air fluorescence yield
is defined as the number of photons produced by a charged particle
per meter of travel.

3/19/2010

elongation rate

if you plot the different Xmax data on y with the energy on x , the slope is the elongation rate

3/03/2010

Good Root Color Scheme

kBlue,
kGreen+4,
kRed,
kMagenta,
kCyan+2

2/19/2010

Purpose of the past works

2009
Energy_Hist_1208 Count the Number of Events above certain energy
for different time periods.


2010
prog0214 Get Number of Events and Energy information for all datas

2/11/2010

Ideas about research 02/11/2010

1Draw the plot of muon and total signal together, with different color.

2 Look into the program to see how did it do the number of muons.

1/25/2010

.bashrc

Put commands in .bashrc is better
than put that into .profile

It will load automatically you do not need to
do . ~/.profile

1/22/2010

Root Basic Colors and Basic Styles

kWhite, kBlack, kRed, kGreen,
kBlue,
kYellow, kMagenta, kCyan

http://root.cern.ch/root/html/TAttMarker.html#M1

gr4->SetMarkerStyle(20);
21,22,23
is good too

1/14/2010

激励

最后送给“蚁族”们一句《基督山伯爵》中的话:“永远不要忘记,在上帝揭露人的未来以前,人类的一切智慧是包含在这四个字里面的:‘等待’和‘希望’。”只有在等待中积蓄能量,充实自己,才不会在希望来临之时不知所措。

1/11/2010