2009 Summer

• Main works
  • Service work : CMSSW integration build renewal with SDT group
  • Physics : MC generation with PYTHIA8 on CMSSW_3_1_0 and Jet analysis
  • Meeting : High Pt QCD group, SDT group, DQM

• Sub works
  • Computing Lectures: c++ lecture
  • Tutorials : CMSSW_3_1_0 tutorial
  • TDR study
  • LHC Physics study
  • Research paper study

Service work

• I got work for CMSSW Integration Builds
• In the Q/A page, there are many contents making by text. It is difficult to understand information for user or development group.
• I made some plot(rem using rate, cpu time, file size, python script errors rate) for text data. it was successful!
• I added new interface for package dependencies.
• I made Twiki page for Q/A page.
• http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/hykim/

Physics

• PYTHIA8 on CMSSW_3_1_0 :
  • make config file

...MY CONFIG

import ParameterSet?.Config as cms

source = cms.Source("EmptySource")

generator = cms.EDFilter("Pythia8Source",
pythiaPylistVerbosity = cms.untracked.int32(1),
pythiaHepMCVerbosity = cms.untracked.bool(False),
maxEventsToPrint = cms.untracked.int32(1),
crossSection = cms.untracked.double(-1.0),
filterEfficiency = cms.untracked.double(-1.0),
comEnergy = cms.untracked.double(14000.0),
PythiaParameters? = cms.PSet(
myset = cms.vstring('HardQCD:all = on','PhaseSpace:pTHatMin = 20.'),
parameterSets = cms.vstring('myset')
)
)
ProductionFilterSequence? = cms.Sequence(generator)

...END

Using cmsDriver.py , I made the GEN-SIM-DIGI-L1-DIGI to RAW-HLT chain cfg file

...Argument

cmsDriver.py MY_CONFIG.py -s GEN,SIM,DIGI,L1,DIGI2RAW,HLT --eventcontent FEVT --datatier GEN-SIM-RAW --conditions FrontierConditions?_GlobalTag,MC_31X_V3::All --number 100 --no_exec --mc

...END

using summer09 globalTag

cmsRun MY_CFG.py

I got bug result

cms::Exception caught in cmsRun ----
PluginNotFound? BEGIN Unable to find plugin 'Pythia8Source'. Please check spelling of name.
---- PluginNotFound? END

this problems came from CMSSW_3 version does not support PYTHIA8 yet. CMSSW_3 little big has been changed. first of all, generator interface uses EDFilter on CMSSW_3. no more use source. so, if use CMSSW_3 you have to use emptysource for mc data.

Changing Generator interface.

I had tested on CMSSW_3_2_0. It also does not support PYTHIA8 interface.

CMSSW_3_2_0 is designed for CRAFT09

I got information about PYTHIA8 interface on CMSSW_3_1_0 at hypernew.

I changed my code.

...MY CONFIG

import ParameterSet?.Config as cms

source = cms.Source("EmptySource")

generator = cms.EDFilter("Pythia8GeneratorFliter",
pythiaPylistVerbosity = cms.untracked.int32(1),
pythiaHepMCVerbosity = cms.untracked.bool(False),
maxEventsToPrint = cms.untracked.int32(1),
crossSection = cms.untracked.double(-1.0),
filterEfficiency = cms.untracked.double(-1.0),
comEnergy = cms.double(14000.0),
PythiaParameters? = cms.PSet(
myset = cms.vstring('HardQCD:all = on','PhaseSpace:pTHatMin = 20.'),
parameterSets = cms.vstring('myset')
)
)
ProductionFilterSequence? = cms.Sequence(generator)

...END

on CMSSW_3_1_0, must use "Pythia8GeneratorFilter"

I got 10 raw events.

next step. Reconstruction.

using cmsDriver I made reco events.

...Argument cmsDriver.py step2 -s RECO --filein file:My_cfg_py_GEN_SIM_DIGI_L1_DIGI2RAW_HLT.root --datatier GEN-SIM-DIGI-RAW-HLT-RECO --eventcontent RECO --conditions FrontierConditions?_GlobalTag,MC_31X_V3::All --no_exec --mc

I got root file with this message.

• Z -> qqbar evnet
  • Generating events on PYTHIA6 or PYTHIA8

source = cms.Source("EmptySource")

from PythiaUESettings?_cfi import *
generator = cms.EDFilter("Pythia6GeneratorFilter",
pythiaHepMCVerbosity = cms.untracked.bool(False),
maxEventsToPrint = cms.untracked.int32(0),
pythiaPylistVerbosity = cms.untracked.int32(0),
filterEfficiency = cms.untracked.double(1.0),
comEnergy = cms.double(14000.0),
PythiaParameters? = cms.PSet(
pythiaUESettingsBlock,
processParameters = cms.vstring('MSEL = 11 ',
'MDME( 174,1) = 1 Z decay into d dbar',
'MDME( 175,1) = 1 Z decay into u ubar',
'MDME( 176,1) = 1 Z decay into s sbar',
'MDME( 177,1) = 1 Z decay into c cbar',
'MDME( 178,1) = 0 Z decay into b bbar',
'MDME( 179,1) = 0 Z decay into t tbar',
#'MDME( 182,1) = 0 Z decay into e- e+',
#'MDME( 183,1) = 0 Z decay into nu_e nu_ebar',
#'MDME( 184,1) = 0 Z decay into mu- mu+',
#'MDME( 185,1) = 0 Z decay into nu_mu nu_mubar',
#'MDME( 186,1) = 0 Z decay into tau- tau+',
#'MDME( 187,1) = 0 Z decay into nu_tau nu_taubar',
'CKIN( 1) = 40. !(D=2. GeV?)',
'CKIN( 2) = -1. !(D=-1. GeV?)',
'MDME(89,1) = 0 ! no tau->electron',
'MDME(90,1) = 0 ! no tau->muon'),
# This is a vector of ParameterSet? names to be read, in this order
parameterSets = cms.vstring('pythiaUESettings',
'processParameters')
)
)

ProductionFilterSequence? = cms.Sequence(generator)

...end

RECO step report

• FireWorks? examples

• cmsDriver for GEN-SIM-DIGI-RECO chain

• Step1

cmsDriver.py My_cfi.py

-s GEN,SIM,DIGI,L1,DIGI2RAW,HLT #Generation chain step

--eventcontent FEVTDEBUG --datatier GEN-SIM-RAW-RECO

--conditions FrontierConditions?_GlobalTag,MC_31X_V3::All

--number 100 --no_exec --mc

• Step2

cmsDriver.py step2 -s RAW2DIGI?,RECO

--filein file:My_cfi_py_GEN_SIM_DIGI_L1_DIGI2RAW_HLT.root

--datatier GEN-SIM-DIGI-RAW-HLT-RECO --eventcontent RECO

--conditions FrontierConditions?_GlobalTag,MC_31X_V3::All

-n 100 --no_exec --mc

• QCD Z' Dijet event (my physics subtopic!)

PYTHIA6 Setting

This configuration make Z' boson to qqbar process excepting bbbar and ttbar process. In the processing, two out going partons will become Di-Jet.

Results on Fireworks

short comment

I expected Di-jet, but on results there are three jets. I have supposed two jets form two partons and third jet from gluon radiation.

Meeting

• 17/Jul/2009 meeting with SDT group
  • Got a introduction of CMSSW IB
  • Got a job of CMSSSW IB
    • make plot of rem state
    • make plot of pycheck

• 23/Jul/2009 QCD HighPT? meeting
  • Update on event shppes study
  • Validation of 31X MC samples

• 29/Jul/2009 All CMS meeting

• 19/Aug/2009 Biweekly Korea CMS group meeting
  • I had presentation about CMSSW IB
  • Introducing IB and using IB
  • http://indico.cern.ch/getFile.py/access?contribId=0&resId=6&materialId=slides&confId=66494

-- HyunYongKim - 17 Aug 2009