1 下载iso 文件
2 右键单击文件,选择打开方式为
windows 7 光盘映像刻录机
3 选择刻录后验证光盘
4 刻录
12/25/2010
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
- Do the vem trace integration
- Apply the energy events weighting of the QGSjetII events
12/17/2010
Previous Work Summary month by month for reference later
- 2010
- Jan To identify the muon peaks in the simulated ADST files
- Feb Energy histogram of surface detector
- Mar Energy histogram change with years
- Apr Different method for jumps on simulation files
- May Jump method fitting
- Jun Hybrid Curvature Xmax relations
- Jul Hybrid Curvature Shower age correlations
- Aug Plot muon and EM traces separately
- Sep Calculate muon/em ratio
- Oct Calculate average muon jumps and compare data and simulation
- Nov Calculate average muon jumps and compare data and simulation
- Dec Calculate average muon jumps and compare data and simulation
- 2009
- Jan
- Feb
- Mar
- Apr
- May
- Jun
- Jul
- Aug
- Sep
- Oct
- Nov
- 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
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
- 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. - After that, click on "Select: All" on the blue area above the first email in the list.
- 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.
12/13/2010
12/10/2010
12/06/2010
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
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
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
11/15/2010
Uncertainty calculations
1 adding or subtracting
a=10+/-1b=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
文件内容(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
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.
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
a massive amount of complete simulation files under
/sps/pauger/Simulations/libraries/augermgr.ccin2p3
9/29/2010
PowerPoint 2007: Delete all slide notes
- Click the Office button, point to Prepare, and then click Inspect Document.
- Click Inspect.
- 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.
- 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.
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
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.
Since e.g. indicates a partial list, it is redundant(多余的) to add “etc.” at the end of a list introduced by this abbreviation.
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
"\\" 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.
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
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)
8/08/2010
7/14/2010
6/29/2010
Segmentation fault
When you see that, first check to see if you define the canvas before you define the graph.
6/28/2010
6/25/2010
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:
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
5/25/2010
5/17/2010
5/14/2010
First wikipedia article
Take two people for a few days, but looks nice now.
http://en.wikipedia.org/wiki/Shear_wave_splitting
http://en.wikipedia.org/wiki/Shear_wave_splitting
5/11/2010
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!
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
4/28/2010
Grep ignore case
- By default, grep is case-sensitive (use -i to ignore case).
- By default, grep ignores the context of a string (use -w to match words only).
- 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();
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.
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
4/19/2010
Overflow and Underflow display in root
gStyle->SetOptStat(111111);
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
$ 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();
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
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);
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
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.
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
3/24/2010
一些概念
AGN:
Active Galactic Nucleus. The core of an active galaxy responsible for the enormous luminosity. Many contains massive black holes.
The air fluorescence yield
is defined as the number of photons produced by a charged particle
per meter of travel.
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/16/2010
3/03/2010
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
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.
2 Look into the program to see how did it do the number of muons.
1/25/2010
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
订阅:
评论 (Atom)





