The collision era is here!

fireworks_pat_123596.jpg

This first data twiki refers to the 2009 Data collisions running. If you want the 2010 7 TeV collisions please see WorkBookPATOnData.

Data samples

CMSSW_3_3_5

Run 122314

The first run of first data can be accessed with this snippet:

# Add the files for run 122314
readFiles = cms.untracked.vstring()
secFiles = cms.untracked.vstring()

readFiles.extend( [
       '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/122/314/F62B040F-6CD8-DE11-9007-001D09F24664.root',
       '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/122/314/F4387297-74D8-DE11-996C-001D09F24F1F.root',
       '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/122/314/EE7B1AC4-6CD8-DE11-97BB-0030487A1FEC.root',
       '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/122/314/9264A28F-87D8-DE11-83A1-001D09F24763.root',
       '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/122/314/7C7495C2-71D8-DE11-ACF2-001D09F248F8.root',
       '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/122/314/4CAB3B6C-6BD8-DE11-845C-000423D9890C.root',
       '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/122/314/02E8544C-70D8-DE11-85CF-001617C3B66C.root' ] );
process.source.fileNames = readFiles
process.source.firstRun = cms.untracked.uint32(122314)

Runs 123575 and 123596

These two runs can be accessed with this snippet:

readFiles = cms.untracked.vstring()
secFiles = cms.untracked.vstring()

readFiles.extend( [
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/575/102A5CB2-B3E1-DE11-A45E-001617E30CE8.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/575/3891C1EB-A3E1-DE11-AB9A-0030487C6090.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/575/46EFF8F0-AAE1-DE11-9232-0016177CA7A0.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/575/82355A91-A4E1-DE11-A5AA-003048D37456.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/575/D0CB515A-A0E1-DE11-8592-001617C3B6FE.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/575/F04F8FC7-A6E1-DE11-AD29-000423D951D4.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/0A71AE7F-4DE2-DE11-8B2F-001D09F251CC.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/14D22A92-62E2-DE11-9B14-000423D990CC.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/263E80C6-41E2-DE11-A194-001617C3B66C.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/26EC6965-4CE2-DE11-ABCF-003048D373AE.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/2A6903C0-68E2-DE11-B6BA-001D09F28D54.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/2EC732CC-5CE2-DE11-A781-001D09F290BF.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/4693BF9A-40E2-DE11-BDBD-000423D944F8.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/54A3A373-4CE2-DE11-8658-000423D99AAA.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/5C4B3A8E-63E2-DE11-A02A-000423D99E46.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/68ED95B5-50E2-DE11-B4C8-001D09F27003.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/7E34865F-45E2-DE11-896A-000423D98F98.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/7EFA67BE-66E2-DE11-AE17-001617C3B79A.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/903D2066-61E2-DE11-9F6E-0019B9F704D6.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/A0FC9BDF-65E2-DE11-A6A1-000423D174FE.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/A4D9D21B-58E2-DE11-8F7A-000423D986A8.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/B60CC58F-5CE2-DE11-9FC7-001D09F24DDF.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/D0586F90-5FE2-DE11-8976-001D09F24691.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/E432BCD7-55E2-DE11-B670-001617C3B6CC.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/123/596/F67BCF17-48E2-DE11-98B1-000423D94534.root'

        ] );

If desired, you can select the bsc "physics bits" for this run like this:

# configure HLT
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('0 AND (40 OR 41)')

# let it run
process.p = cms.Path(
    process.hltLevel1GTSeed*
    process.patDefaultSequence
    )

Run 124120: 2360 GeV collisions

This run can be accessed with this snippet:

readFiles = cms.untracked.vstring()
secFiles = cms.untracked.vstring()

readFiles.extend( [
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/124/120/F08F782B-77E8-DE11-B1FC-0019B9F72BFF.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/124/120/EE9412FD-80E8-DE11-9FDD-000423D94908.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/124/120/7C9741F5-78E8-DE11-8E69-001D09F2AD84.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/124/120/44255E49-80E8-DE11-B6DB-000423D991F0.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/124/120/3C02A810-7CE8-DE11-BB51-003048D375AA.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/124/120/04F15557-7BE8-DE11-8A41-003048D2C1C4.root',
        '/store/data/BeamCommissioning09/MinimumBias/RECO/v2/000/124/120/04092AB7-75E8-DE11-958F-000423D98750.root'
        ] );

If desired, you can select the bsc "physics bits" for this run like this:

# configure HLT
process.load('L1TriggerConfig.L1GtConfigProducers.L1GtTriggerMaskTechTrigConfig_cff')
process.load('HLTrigger/HLTfilters/hltLevel1GTSeed_cfi')
process.hltLevel1GTSeed.L1TechTriggerSeeding = cms.bool(True)
process.hltLevel1GTSeed.L1SeedsLogicalExpression = cms.string('0 AND (40 OR 41)')

# let it run
process.p = cms.Path(
    process.hltLevel1GTSeed*
    process.patDefaultSequence
    )

You can switch to the 2360 MC-derived JEC by this:

from CMS.PhysicsTools.PatAlgos.tools.jetTools import *
switchJECSet( process, "2360GeV")

CMSSW_3_4_1

The 341 re-reco of the minimum bias sample is done. This prescription is using a secondary dataset (skim) of this sample.

Runs 123575 and 123596

These two runs can be accessed with this snippet:

readFiles = cms.untracked.vstring()
secFiles = cms.untracked.vstring()

readFiles.extend( [
        '/store/data/BeamCommissioning09/MinimumBias/RAW-RECO/SD_InterestingEvents-Dec19thSkim_341_v1/0006/5A5773A3-B9ED-DE11-A99A-0026189437F8.root',
        '/store/data/BeamCommissioning09/MinimumBias/RAW-RECO/SD_InterestingEvents-Dec19thSkim_341_v1/0006/CA156CCC-B9ED-DE11-B022-00261894390B.root'
        ] );

The bsc "physics bits" are already required for this secondary dataset.

Run 124120: 2360 GeV collisions

This run can be accessed with this snippet:

readFiles = cms.untracked.vstring()
secFiles = cms.untracked.vstring()

readFiles.extend( [
        '/store/data/BeamCommissioning09/MinimumBias/RAW-RECO/SD_InterestingEvents-Dec19thSkim_341_v1/0005/B641315C-ACED-DE11-82E1-0030486792B6.root',
        '/store/data/BeamCommissioning09/MinimumBias/RAW-RECO/SD_InterestingEvents-Dec19thSkim_341_v1/0005/8EB877E1-A5ED-DE11-9680-002618943826.root',
        '/store/data/BeamCommissioning09/MinimumBias/RAW-RECO/SD_InterestingEvents-Dec19thSkim_341_v1/0005/52B38DBB-A1ED-DE11-AAB0-00248C0BE014.root',
        '/store/data/BeamCommissioning09/MinimumBias/RAW-RECO/SD_InterestingEvents-Dec19thSkim_341_v1/0005/404C0D27-A8ED-DE11-87B2-00261894397F.root',
        '/store/data/BeamCommissioning09/MinimumBias/RAW-RECO/SD_InterestingEvents-Dec19thSkim_341_v1/0005/3C4F0CB2-A1ED-DE11-9200-0026189438E8.root',
        '/store/data/BeamCommissioning09/MinimumBias/RAW-RECO/SD_InterestingEvents-Dec19thSkim_341_v1/0005/2AA408F8-A3ED-DE11-8B3B-002618943821.root'
        ] );

The bsc "physics bits" are already required for this secondary dataset.

You can switch to the 2360 MC-derived JEC by this:

from CMS.PhysicsTools.PatAlgos.tools.jetTools import *
switchJECSet( process, "2360GeV")

Setting up the environment

CMSSW_3_3_5

The prescription here is almost identical to CMSSW_3_3_5 "out of the box" except for:

  • The addition of a couple of config files to allow you to run over the first run of the first data
  • Updating to the JEC for 900 GeV collisions, more info here.
  • Updating to the JEC for 2360 GeV collisions, more info here.

scram p CMSSW CMSSW_3_3_5
cd CMSSW_3_3_5/src
cmsenv
addpkg PhysicsTools/PatAlgos V07-09-34-06
addpkg  CondFormats/JetMETObjects V01-08-08-09
addpkg JetMETCorrections/Configuration V01-08-23-05
scram b -j4
cd PhysicsTools/PatAlgos/test
cmsRun patLayer1_fromRECO_900GeV_firstdata_cfg.py

or

cmsRun patLayer1_fromRECO_2360GeV_firstdata_cfg.py

CMSSW_3_4_2

The prescription here is almost identical to CMSSW_3_4_1 "out of the box" except for:

  • The addition of a couple of config files to allow you to run over the first run of the first data
  • Updating to the JEC for 900 GeV collisions, more info here.
  • Updating to the JEC for 2360 GeV collisions, more info here.

scram p CMSSW CMSSW_3_4_2
cd CMSSW_3_4_1/src
cmsenv
addpkg PhysicsTools/PatAlgos 
cd PhysicsTools/PatAlgos/test
cmsRun patLayer1_fromRECO_900GeV_firstdata_cfg.py

or

cmsRun patLayer1_fromRECO_2360GeV_firstdata_cfg.py

How to run on Monte Carlo

It is possible to run on the latest Monte Carlo with the same configuration, it is sufficient to comment out the removal of MC matching:

#removeMCMatching(process, 'All')

Be advised that this file is keeping everything in the event so if you want to reduce your output size, adjust this accordingly.

The MinBias MC samples can be found from here.

Running Fireworks

To use standalone Fireworks with PAT objects in this release, please take the latest version (34x) which has an updated DataFormats for PAT candidates (which changed).

A configuration to run fireworks on the PAT objects is here.

Example FWLite Macro to plot everything

Make sure you followed the instructions at WorkBookFWLite for setting up the root logon.

Also make sure to change the name of the input root file in the macro plotAllPATObjects.C to the one you have made. The default name is given in the line TFile * file = new TFile("PATLayer1_Output.fromAOD_full.root"); in the macro plotAllPATObjects.C.

cd PhysicsTools/PatAlgos/test/
root -l
root [0] gSystem->CompileMacro("fwlite/plotAllPATObjects.C", "k");
root [1] plotAllPATObjects();                                     
root [2] .q

Example Data Plots from PAT

The PAT MET and Jet distributions are shown here:

http://indico.cern.ch/getFile.py/access?contribId=4&resId=0&materialId=slides&confId=75029

Note: There are a lot of other plots shown than the simple analysis macro above (which serves as a skeleton).

Related links

For up-to-date details on the Commissioning of first data see here:

https://twiki.cern.ch/twiki/bin/viewauth/CMS/CollisionsDec2009

https://twiki.cern.ch/twiki/bin/viewauth/CMS/FirstCollisionsAnalysis

Update history:

  • 9-Dec-09
    • Added Monte Carlo description.
    • Advised to update to latest Fireworks.
    • Renamed default file to reflect COM energy of 900 GeV.
    • Added MC-based JEC for 900 GeV.
    • Require bit "0" in addition to "40 or 41" to make sure the collision is in the appropriate BX.
  • 14-Dec-09
    • Added MC-based JEC for 2360 GeV, and added a new file for latest run at 2360 GeV.
  • 22-Dec-09
    • Added prescription for 341 re-reco.

-- BenediktHegner - 09-Dec-2009

Topic attachments
I Attachment History Action Size Date Who Comment
JPEGjpg CollisionEvent.jpg r1 manage 837.6 K 2009-11-24 - 16:28 SalvatoreRRappoccio  
JPEGjpg fireworks_pat_123596.jpg r1 manage 230.9 K 2009-12-06 - 23:01 SalvatoreRRappoccio  
JPEGjpg jetEtaVsPhi.jpg r1 manage 35.8 K 2009-11-24 - 17:39 SalvatoreRRappoccio  
JPEGjpg jetPt.jpg r1 manage 16.0 K 2009-11-24 - 17:39 SalvatoreRRappoccio  
JPEGjpg metPhi.jpg r1 manage 18.0 K 2009-11-24 - 17:41 SalvatoreRRappoccio  
JPEGjpg metPt.jpg r1 manage 16.3 K 2009-11-24 - 17:41 SalvatoreRRappoccio  
JPEGjpg nJet.jpg r1 manage 15.9 K 2009-11-24 - 17:40 SalvatoreRRappoccio  
Unknown file formatfwc pat.fwc r1 manage 42.0 K 2009-12-06 - 23:04 SalvatoreRRappoccio  
Edit | Attach | Watch | Print version | History: r14 < r13 < r12 < r11 < r10 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r14 - 2010-04-28 - RogerWolf
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    CMSPublic All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback