DMoon's Monthly Logs : 2008 Jan.
2007 01 10 (Thu)
지현이랑 일하고 있음..
BBbar_2_Jpsi.cfg
필터가 제대로 작동하고 있는 것 같지 않음..
-6 6 까지 나옴..
다시 해보는 중..
gen, psim, pdigi, reconstruction
내경우.. L1 L2를.. 프리셀렉션..
2008 01 24 (Thu)
어제 CERN 에 도착..
오늘부터 일 시작.. HI
https://twiki.cern.ch/twiki/bin/view/CMS/CmsHiSoftware#J_psi_and_upsilon_generation_acc
jpsi & upsilon Heavy Ion test 시작
testjpsi.cfg - pure Pythia
testjpsifilter_gen.cfg - Pythia + possibility to filter muons
testjpsifilter_sim.cfg - Pythia+Filter+SIM-RECO
For upsilon generation one has to change in cfg:
replace
PythiaSource?.pdfFile = 'upsipbpb'
replace
ParticleID? = 553
replace
PythiaSource?.Etamin = -6.5
replace
PythiaSource?.Etamax = 6.5
2008 01 25 (Fri)
1. 헤비아이온 소프트 웨어 시도, L1 L2 트리거
2. L3 다이뮤온 트리거 만들기
3. lowqualityFlag 비교..
시립대 써버의 고장으로 하이레벨 트리거를 보지 못하고 있음..
testjpsifilter_gen.cfg
testjpsifilter_sim.cfg
testjpsifilter_digi.cfg
testjpsifilter_rec.cfg
로 나눠 진행..
gen은 5000개가 5분안에 만들어지지만.. 나머지 단계에서 시간이 많이 소요됨..
sim 5000개 진행 하고 들어감
2008 01 26 (Sat)
2560개에서 terminate 됨
2000개 gen 과 sim 하는 중..
끝나면.. rec 까지 해야함..
digi 하는 중..
reco 마무리..
PhysVal?-JPsiToMuMu-NoPU_SingleParticleFilter_gen_sim_dig_rec_e2000.root
/castor/cern.ch/user/d/dmoon/hidata_jpsi
밑에 저장..
cand 분석(generator level)
hi_jpsi_e2000_singleParticleFilter_JpsiAna.root
hi_jpsi_e2000_singleParticleFilter_JPsiCand.root
high level trigger 할 차례
생각과는 많이 다름..
2000개 생산..
eta, mass, pt 가 다 조금씩 이상함..
확인중..
생산 개수는 맞지만..
pt, eta, 의 모양이 굉장히 flat하다.
어쩌면.. gen과 sim을 같이해서 그럴지도 모른다.
gen, psim을 같이 해보자.
upsilon 도 테스트 중..
그것도 확인 요함..
개수의 문제가 있다.
2008 01 28 (Mon)
psim을 사용한 것과 simulation을 사용한 것과는 차이가 있다.
사실 차이가 많이 난다.
들여다 봐야 한다.
untracked를 쓰는 것과 쓰지 않는 것의 차이를 공부하자..
simulation을 사용하는 경우.. event가 0가 나온다.
* simulation 은 psim과 pdigi를 포함하고 있다.
sequence simulation = {psim,pdigi,genParticleCandidates}
sequence psim = {
VtxSmeared?, g4SimHits }
total chain 이용해서 simulation하고 있다.
-> 1460 에서 terminated 됨..
예전에 했던 파일을 이용해서 다시 hlt하는 중..
그결과를 메일 보내야 함..
내일은 upsilon 으로 하려고 함..
참고 : untracked
InputTag? GMTInputTag? = l1GmtEmulDigis
hi_jpsi_hlt_e2000.root 생성
TestMuL1L2?모듈 사용중..
2008 01 29 (Tue)
high eta region 에서 missing 이 있는 이유는
CSCTrackFinder?.cfi 때문이다.
다시 default로 하고 있다.
L2 pt 가 200 이상에서도 나타난다.
2008 01 30 (Wed)
제너래이션 된 정보가 flat하다..
다시 확인 하는 중.. 올가의 정보는 다르다.
Olga가 자신의 프로젝트 디렉토리를 알려주어서 거기서 해보는 중..
~kodolova/scratch0/CMSSW_TEST/cmssw167_test/src/GeneratorInterface/Pythia6Interface/test
analysis 는
~kodolova/scratch0/CMSSW_TEST/cmssw167_test/src/UserCode/OlgaKodolova/test
Quarkonia.C
Quarkonia.h
analysis_l1l2.cfg
이것을 ana_hi.cfg로 바꿔서 하는 중..
ana_hi.cfg 파일안에.. generated 된 파일...
PhysVal?-JPsiToMuMu-NoPU_SingleParticleFilter_gen.root을 넣고
run..
test.root가 생기고
그 다음 root -l test.root
하고 나서
다음과 같이 하면
Root > .L Quarkonia.C
// Root > Quarkonia t
// Root > t.GetEntry(12); // Fill t data members with entry number 12
// Root > t.Show(); // Show values of entry 12
// Root > t.Show(16); // Read and show values of entry 16
// Root > t.Loop(); // Loop on all entries
eta, pt 플롯을 얻을 수 있다.
실행 시킬때는
eta.Draw();
pt.Draw();
라고 하면 된다.
아무래도 cms167에서 프로그램을 설치할때 문제가 생겼던 것으로 생각된다.
현재 일하는 project folder
/afs/cern.ch/user/d/dmoon/scratch0/cms167/src/GeneratorInterface/Pythia6Interface/test
testjpsifilter_gen.cfg, event : 2000
sim, digi, reco chain.
server 213
tip: 링크 ln -s aaa.root bbb.root
>
bbb.root -> aaa.root
root 에서 branch entuple 불러서 fillup 할 수 있는 코드 만들기
.ls
InputClassName?(ex:Quarkonia)->MakeClass("OutputCodeName(ana)");
example :
ROOT>
Quarkonia->MakeClass("ana");
-> Olga가 준 코드 분석
추가 histogram 이름
TH1F?* etahist = new
TH1F?("eta"," ",100,-10., 10.);
TH1F?* pthist = new
TH1F?("pt"," ",100,0., 40.);
for(Int_t i = 0; i<NumPart;i++){
//if(abs(Code[i]) == 443)
if(abs(Code[i]) == 553)
//if(abs(Code[i]) == 13)
{
etahist->Fill(eta[i]);
pthist->Fill(pt[i]);
}
}
Quarkonia.C 를 이용해서 edmHepMCproduct 의 정보에서 빼옴..
같은 방식으로 다시 해서 jpsi는 gen이 되었으나
upsilon의 개수는 문제가 있음..
200000 개 생산 했지만.. 8000여개 정도의 upsilon이 나옴..
오늘 일 한 것 정리..
우선.. 2000개의 jpsi 생산..
4000개의 뮤온 만들어짐
프리셀렉션..
1800개
L1 트리거 개수 27개
L2 트리거 개수 26개
프리셀렉트 에피션시
1800/4000 = 0.45
L1 트리거 에피션시
27/1800 = 0.015
--
DongHoMoon - 21 Feb 2008