DMoon's Monthly Logs : 2007 Dec.

2007 12 01 (Sat)

L2 를 하기 위해 L1을 HLT에서 하는 중.. 문제는 그냥 L1을 하면.. filtering 이 안되는 듯 하다.. L1정보가 실제의 muon의 개수보다 많기 때문이다.

그런데 DBS 샘플을 보면 그렇지 않게 나온다. 따라서 내 L1은 문제가 있다.. ㅡㅡ; 그것을 찾기 위해 gen 파일을 보내고 있다.. 거기서 pdigi, raw를 할 생각임.. 그런다음.. level 2를 하려고함..

cmsRun RelVal?_Digi_Digi2Raw.cfg Unable to find plugin 'RCTConfigProducers'. Please check spelling of name.

lowqualityFlag 관련 정보 https://twiki.cern.ch/twiki/bin/view/CMS/GMTEmulator

2007 12 03 (Mon)

DBS raw data 옮기는 중 오늘 할일 DBS로 L2 시도.. CandFrommumu? 모듈 조사 L3와 globla muon 비교.. efficiency 구하기..

CandFromMumu?

module CandFromMumu? = CandCombiner? { string decay = "muTracks@+ muTracks@-" string cut = "1.00 < mass < 6.00 & pt > 4.0" }

그럼 이제 muTrack 만들기 module allMuTracks = allTracks from "PhysicsTools/RecoCandAlgos/data/allTracks.cfi" replace allMuTracks.src = ctfWithMaterialTracksMumu replace allMuTracks.particleType = "mu-"

module muTracks = CandSelector? { InputTag? src = allMuTracks string cut = "pt > 4." bool filter = false }

PhysicsTools?/RecoCandAlgos/data/allTracks.cfi allTracks ctfWithMaterialTracksMumu

L1Muon?.cfg -> jpsi path에서 double muon type 3가 필요하기 때문에 L1Muon?.cfg가 필요하다.. 해보는 중.. Error occured while creating HLTLevel1Seed? with label L1MuLevel1Seed? MissingParameter?: Parameter 'L1ExtraParticleMap' not found.

이 에러가 나는 이유는 include "HLTrigger/Muon/data/CommonModules_2x1033.cff" 에서 다음을 호출하게 되는데 include "HLTrigger/Configuration/data/common/HLTSetupCase2.cff 이 안에는 module hlt2GetRaw = hltGetRaw from "HLTrigger/HLTanalyzers/data/hltGetRaw.cfi" module l1SeedBegin = hltLevel1Seed from "HLTrigger/HLTfilters/data/hltLevel1Seed.cfi" sequence hltBegin = { hlt2GetRaw & l1SeedBegin } l1SeedBegin 에서 이미 hltLevel1Seed 를 호출하면서 그 안에 있는 InputTag? L1ExtraParticleMap? = l1extraParticleMap 에러가 발생하고 있다.

HLTrigger/Muon/data/CommonModules_1032_HLTExercise.cff

sequence l1muonreco = { hltBegin } sequence l2muonrecoNocand = { doLocalMuon, L2MuonSeeds?, L2Muons? } sequence l2muonreco = { l2muonrecoNocand, L2MuonCandidates? } sequence l2muonisoreco = { doRegionalCaloForMuons, L2MuonIsolations? } sequence l3muonrecoNocand = { (doLocalPixel & doLocalStrip), L3muonAltReco? } sequence l3muonreco = { l3muonrecoNocand, L3MuonCandidates? } sequence l3muonisoreco = { regionalTrackingForMuIso2, L3MuonIsolations? }

l1muonreco -> ????

계속 에러가 발생

hlt 말고 L2MuonProducer? 이용 중

L2MuonSeedGenerator? 까지는 잘 되고 있으나 L2Muons?L2MuonCandidates?가 아직 안되고 있음..

- Tip - N/A : 상관없다. BX : Bunch Corssing


2007 12 04 (Tue)

CMSSW/ HLTrigger/ Muon/ src/ HLTMuonDimuonFilter?.cc 코드를 응용해서 L2 part를 만들수 있었고 현재 L2 정보를 불러 올 수 있게 되었다..

int ilmt1 = 0; for (cand1=mucands->begin(); cand1!=mucands->end(); cand1++) { TrackRef? traj1 = cand1->get();

chargel2[ilmt1]=traj1->charge(); etal2[ilmt1]=traj1->eta(); ptl2[ilmt1]=traj1->pt(); phil2[ilmt1]=traj1->phi(); std::cout<<" L2 Muon Candidates "<<endl; std::cout<<" L2 Muons charge "<<chargel2[ilmt1]<<endl; std::cout<<" L2 Muons eta "<<etal2[ilmt1]<<endl; std::cout<<" L2 Muons pt "<<ptl2[ilmt1]<<endl; std::cout<<" L2 Muons phi "<<phil2[ilmt1]<<endl;

hL2Charge->Fill(chargel2[ilmt1]); hL2Eta->Fill(etal2[ilmt1]); hL2Pt->Fill(ptl2[ilmt1]); hL2Phi->Fill(phil2[ilmt1]); ilmt1++;

}

L1Analyzer?를 돌리고 있다..

2007 12 05 (Wed)

to do list 1.L2 Muon producer : L2 muon candidates 와 L2 muons간의 관계 L2 Muons를 먼저 하는 것인지 아니면.. L2 muons cand를 먼저 하는 것인지.. 확인 2. L2 muon candidates 의 결과가 L1의 결과와 상반 되는 부분이 있음 3. lowqualityFlag 공부 마무리 4. L1Anayzer? code 다시 보기


> lowqualityFlag code 이해 CMS IN 2004/022

Specification of the Interface Between the Regional Muon Triggers and the Global Muon Trigger D. Acosta 외..~~

muon candidates 는 pt, charge sign, eta, phi, quality code 등의 정보를 갖고 있다. 각각의 chamber 들은 (DT, CSC, RPC(barrel), RPC(endcap)), 4개의 muon을 candidates 하게 된다. 그러면 전체 16개의 muon 들 중 최종 4개의 뮤온이 GMT 로 들어가게 된다. 그리고 GMT bit 는 GCT(Globla Calorimeter Trigger) 로 부터 특정 threthhold의 에너지 deposit 된 것들을 통과 시킨다.(MIP signal)

GMT 와 CSC의 quality code 를 이해해야 한다. table 5 in CMS IN 2004/022 table 4 in CMS IN 2004/006 을 연결 시켜 생각해 봐야한다. https://twiki.cern.ch/twiki/bin/view/CMS/GMTEmulator 의 accept and rejection 이해

http://cmslxr.fnal.gov/lxr/source/DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTCand.h#082 /// get quality

/// Quality codes: /// /// 0 .. no muon /// 1 .. beam halo muon (CSC) /// 2 .. very low quality level 1 (e.g. ignore in single and di-muon trigger) /// 3 .. very low quality level 2 (e.g. ignore in single muon trigger use in di-muon trigger) /// 4 .. very low quality level 3 (e.g. ignore in di-muon trigger, use in single-muon trigger) /// 5 .. unmatched RPC /// 6 .. unmatched DT or CSC /// 7 .. matched DT-RPC or CSC-RPC /// /// attention: try not to rely on quality codes in analysis: they may change again ///

2007 12 06 (Thu)

L2 part modified I can get L2 information.

CandShwallowSelector? 이용한 L1 L2 jpsi reconstruction OL1 을 위한 Trigger 공부

2007 12 10 (Mon)

50000개 짜리 gen 하는 중.. 시작 시간 11:15 -> 20:50 거의 9시간 20분정도..

2007 12 11 (Tue)

10000개의 data를 생산 하였으나 1180개의 jpsi가 생산되었음.. 0.018 efficiency

digiToRaw를 위해 지금 digi 하는 중..(30만개) 기대 jpsi 수 : 4620 개

뮤온 : 9240

RCTConfigProducers? 늘 이것이 문제다..

CERN에서는 어떻게 되는지 살펴보기 위해 해보는 중..

hltBToJPsiTomumu.cfg

> Cannot find parameter block CkfBaseTrajectoryFilter?_block errot

cvs co 할때 cvs co -r $CMSSW_VERSION module_name

Olga's mail 2007 12 12(수) I suppose, we can estimate efficiency, as:

simu_eff * trig_eff = 0.0154 * trig_eff

What is your trig_eff over the generated sample?

-> DBS sample No 01

gen Muon : 12386

L1 muon : 6043

global muons & all muons: 5148

L1Trigger? eff : 6043/12386 = 0.487 0.0154르 곱하면 = 0.00751

지현이가 했던 버젼 1_6_6 으로 바꾼 결과 아주 잘 돌아감.. digi를 해서 바로 넣어도 됨.. 그대신.. HLTSetup.cff 에서 case1로 바꿔 줘야 함

한 것 정리

cvs co -r $CMSSW_VERSION HLTrigger/Muon cvs co -r $CMSSW_VERSION HLTrigger/btau cvs co -r $CMSSW_VERSION HLTrigger/Configuration/data/common cvs co -r $CMSSW_VERSION RecoMuon?/ cvs co -r $CMSSW_VERSION L1Trigger?/

cd /HLTrigger/Muon/test/ cmsRun Jpsimumu.cfg

2007 12 13 (Thu)

현재 시립대 써버에서 300000개짜리 jpsi 데이타를 트리거 하고 있다. 잘 돌아가고 있다.. 데이타들을 다시 정리할 필요가 있음..

MuonCloneSwallower?를 사용해서 jpsi의 수를 reco해봐야 함.. 그리고 L1Analyzer?도 해 볼 필요 있음.. Tirgger .. ㅡㅡ;

그렇지만.. 우선은.. Trigger 하는데 모든 source 가 사용되어 다른일은 하기 어려움..

2007 12 17 (Mon)

Jpsi HLT 정리 HLTrigger/Muon/test/Jpsimumu.cfg > include "HLTrigger/Muon/data/PathJpsimumu.cff"
path pJpsiMM = { jpsiMM }

> include "HLTrigger/Muon/data/CommonModules.cff"

> include "HLTrigger/Configuration/data/common/HLTSetup.cff"

> include "HLTrigger/Configuration/data/common/HLTSetupCase1.cff"

from Digis on input file, running L1Emulator? and L1Extra?

module hlt1GetDigi = hltGetDigi from "HLTrigger/HLTanalyzers/data/hltGetDigi.cfi" include "Configuration/StandardSequences/data/L1Emulator.cff" include "Configuration/StandardSequences/data/L1Extra.cff"

sequence hltBegin = { hlt1GetDigi & (L1Emulator? , L1Extra?) }

include "HLTrigger/Configuration/data/common/HLTSetupCommon.cff" include "HLTrigger/Configuration/data/common/HLTSetupFromDigis.cff"

lowqualityFlag : 1로 놓은 결과 619개 2.1까지.. 619 default lowqualityFlag : 4 2212

2007 12 19 (Wed)

lowqualityFlag 고치기 L1Trigger?/Configuration/data/L1MuonEmulator.cff include "L1Trigger/CSCTrackFinder/data/l1CscTfTrackEmulDigis.cfi" include "L1Trigger/CSCTrackFinder/data/l1CscTfEmulDigis.cfi"

여기서 고치면 된다..

지금 고쳐서 하는 중..

2007 12 26 (Wed)

Olga의 요청으로 pythia information을 알기 위해 eta : -10 ~ 10 pt : 0 이상 으로 필터를 조절.. 30000 개 생산 함.. jpsi_pt_0_20_e30000_gen_160_pt0_eta10.root

2007 12 27 (Thu)

gen 된 것을 digi-reco 중.. jpsi_pt_0_20_e30000_digi_reco_160_pt0_eta10.root cms160/src/DMoon/test/ 밑에서 작업중

-- DongHoMoon - 21 Feb 2008

Edit | Attach | Watch | Print version | History: r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r2 - 2008-02-21 - DongHoMoon
 
KOREA-CMS
This site is powered by the TWiki collaboration platform Powered by Perl This site is powered by the TWiki collaboration platformCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding KoreaCmsWiki? Send feedback