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 다시 보기

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