BcBsPiAnalysis

Ntuple Production (2012Data)

All needed files are in

ntucms1:/disk1/home/ctao/BcBspi_ntuple/CMSSW_5_3_7_patch5/src/BcBspi_ntuple/BcBspi_ntuple

STEP 1: Environment setup and ntuple production

Prepare the environment at CERN

source /afs/cern.ch/cms/LCG/LCG-2/UI/cms_ui_env.csh

cd YOUR_DIR
For Run2012A, B, C
cmsrel CMSSW_5_3_7_patch5
cd CMSSW_5_3_7_patch5/src
For Run2012D
cmsrel CMSSW_5_3_7_patch6
cd CMSSW_5_3_7_patch6/src
cmsenv
create BcBspi_ntuple/BcBspi_ntuple and go to /disk1/home/ctao/BcBspi_ntuple/CMSSW_5_3_7_patch5/src/BcBspi_ntuple/BcBspi_ntuple
cp -a BuildFile.xml bcbspi_ntuple_cfg_HLT.py src/ interface/ YOUR_DIR/CMSSW_5_3_7_patch5/src/BcBspi_ntuple/BcBspi_ntuple
scramv1 b

To set up crab environment,

source /afs/cern.ch/cms/LCG/LCG-2/UI/cms_ui_env.sh
source /afs/cern.ch/cms/ccs/wm/scripts/Crab/crab.sh
voms-proxy-init -voms cms:/cms/twcms -valid 192:0

To submit crab jobs, create test/Run2012*-22Jan2013-v1, where * stands for A/B/C/D. And go to the directory.

cp /disk1/home/ctao/BcBspi_ntuple/CMSSW_5_3_7_patch5/src/BcBspi_ntuple/BcBspi_ntuple/test/Run2012C-22Jan2013-v1/* .
Change the datasetpath and user_remote_dir in crab.cfg (Changing C into B or D is enough.)
crab -create
crab -submit

STEP 2: AnalysisKit

MC Production

Configuration

/eos/cms/store/lhe/8049
/eos/cms/store/lhe/8079
/eos/cms/store/lhe/8047

  • Gen Fragment
import FWCore.ParameterSet.Config as cms
from Configuration.Generator.PythiaUEZ2Settings_cfi import *

generator = cms.EDFilter("Pythia6HadronizerFilter",
    maxEventsToPrint = cms.untracked.int32(2),
    pythiaHepMCVerbosity = cms.untracked.bool(False),
    pythiaPylistVerbosity = cms.untracked.int32(1),
    comEnergy = cms.double(8000.),
    PythiaParameters = cms.PSet(
        pythiaUESettingsBlock,
        processParameters = cms.vstring(
            'MSTP(51) = 10042',
            'MSTP(52) = 2',
            'MSTP(61)=0             ! Hadronization of the initial protons initial-state QCD and QED radiation OFF',
            'MDME(997,2) = 0        ! PHASE SPACE',
            'KFDP(997,1) = 531      ! Bs',
            'KFDP(997,2) = 211      ! pi+',
            'KFDP(997,3) = 0        ! nada',
            'KFDP(997,4) = 0        ! nada',
            'KFDP(997,5) = 0        ! nada',
            'MDME(953,1) = 0',
            'MDME(954,1) = 0',
            'MDME(955,1) = 0',
            'MDME(956,1) = 0',
            'MDME(957,1) = 0',
            'MDME(958,1) = 0',
            'MDME(959,1) = 0',
            'MDME(960,1) = 0',
            'MDME(961,1) = 0',
            'MDME(962,1) = 0',
            'MDME(963,1) = 0',
            'MDME(964,1) = 0',
            'MDME(965,1) = 0',
            'MDME(966,1) = 0',
            'MDME(967,1) = 0',
            'MDME(968,1) = 0',
            'MDME(969,1) = 0',
            'MDME(970,1) = 0',
            'MDME(971,1) = 0',
            'MDME(972,1) = 0',
            'MDME(973,1) = 0',
            'MDME(974,1) = 0',
            'MDME(975,1) = 0',
            'MDME(976,1) = 0',
            'MDME(977,1) = 0',
            'MDME(978,1) = 0',
            'MDME(979,1) = 0',
            'MDME(980,1) = 0',
            'MDME(981,1) = 0',
            'MDME(982,1) = 1    ! J/psi phi',
            'MDME(983,1) = 0',
            'MDME(984,1) = 0',
            'MDME(985,1) = 0',
            'MDME(986,1) = 0',
            'MDME(987,1) = 0',
            'MDME(988,1) = 0',
            'MDME(989,1) = 0',
            'MDME(990,1) = 0',
            'MDME(991,1) = 0',
            'MDME(858,1) = 0         ! J/psi->e+e-',
            'MDME(859,1) = 1         ! J/psi->mumu',
            'MDME(860,1) = 0         ! J/psi->quarks',
            'MDME(656,1) = 1         ! phi->K+K-',
            'MDME(657,1) = 0',
            'MDME(658,1) = 0',
            'MDME(659,1) = 0',
            'MDME(660,1) = 0',
            'MDME(661,1) = 0',
            'MDME(662,1) = 0',
            'MDME(663,1) = 0',
            'MDME(664,1) = 0',
            'MDME(665,1) = 0',
            'MDME(666,1) = 0',
            'MDME(997,1) = 2         ! on for Bc- only, antiparticle',
            'MDME(998,1) = 3',
            'MDME(999,1) = 3',
            'MDME(1000,1) = 3',
            'MDME(1001,1) = 3',
            'MDME(1002,1) = 3',
            'MDME(1003,1) = 3',
            'MDME(1004,1) = 3',
            'MDME(1005,1) = 3',
            'MDME(1006,1) = 3',
            'MDME(1007,1) = 3',
            'MDME(1008,1) = 3',
            'MDME(1009,1) = 3',
            'MDME(1010,1) = 3',
            'MDME(1011,1) = 3',
            'MDME(1012,1) = 3',
            'MDME(1013,1) = 3',
            'MDME(1014,1) = 3',
            'MDME(1015,1) = 3',
            'MDME(1016,1) = 3',
            'MDME(1017,1) = 3',
            'MDME(1018,1) = 3',
            'MDME(1019,1) = 3',
            'MDME(1020,1) = 3',
            'MDME(1021,1) = 3',
            'MDME(1022,1) = 3',
            'MDME(1023,1) = 3',
            'MDME(1024,1) = 3',
            'MDME(1025,1) = 3',
            'MDME(1026,1) = 3',
            'MDME(1027,1) = 3'),
        parameterSets = cms.vstring('pythiaUESettings','processParameters')
    )
)
bcfilter = cms.EDFilter("PythiaFilter",
    ParticleID = cms.untracked.int32(541)
)
mumugenfilter = cms.EDFilter("MCParticlePairFilter",
    Status = cms.untracked.vint32(1, 1),
    MaxEta = cms.untracked.vdouble(2.4, 2.4),
    MinEta = cms.untracked.vdouble(-2.4, -2.4),
    ParticleCharge = cms.untracked.int32(-1),
    ParticleID1 = cms.untracked.vint32(13),
    ParticleID2 = cms.untracked.vint32(13)
)

ProductionFilterSequence = cms.Sequence(generator*bcfilter*mumugenfilter)

Status

STEP0: GEN Efficiency

  • Command line options:
cmsDriver.py Configuration/GenProduction/python/bcbspi_genfragment.py --filein=rfio:/dpm/phys.ntu.edu.tw/home/cms/store/user/ctao/bcvegpy_00.lhe --python_filename BcBspi_GEN.py --fileout BcBspi_GEN.root --step GEN --beamspot Realistic8TeVCollision --conditions START53_V7C::All --pileup NoPileUp --datamix NODATAMIXER --eventcontent RAWSIM --datatier GEN --no_exec -n 100000
cmsRun BcBspi_GEN.py
  • Result:
BcBspi_GEN.root (1 runs, 1 lumis, 62545 events, 3516348328 bytes)
Efficiency = 62.545% +- 0.153%
Total number of event = 100000/0.62545 ~ 160000
Create 400 jobs. Each has 400 events.

STEP1: GEN-SIM

  • Command line options:
cmsDriver.py Configuration/GenProduction/python/bcbspi_genfragment.py --filein=rfio:/dpm/phys.ntu.edu.tw/home/cms/store/user/ctao/bcvegpy_00.lhe --python_filename BcBspi_GEN_SIM.py --fileout BcBspi_GEN_SIM.root --step GEN,SIM --beamspot Realistic8TeVCollision --conditions START53_V7C::All --pileup NoPileUp --datamix NODATAMIXER --eventcontent RAWSIM --datatier GEN-SIM --no_exec -n 10

  • Crab configuration:
[CRAB]
jobtype = cmssw
scheduler = glite
use_server = 0

[CMSSW]
datasetpath = None
generator = pythia
pset = BcBspi_GEN_SIM.py
output_file = BcBspi_GEN_SIM.root
total_number_of_events=160000
events_per_job=400

[USER]
return_data = 0
copy_data = 1
storage_element = T3_TW_NTU_HEP 
user_remote_dir = Bc_Bspi_8TeV_v1 

publish_data = 1
publish_data_name = Bc_Bspi_8TeV_GenFilter_v1
dbs_url_for_publication = https://cmsdbsprod.cern.ch:8443/cms_dbs_ph_analysis_01_writer/servlet/DBSServlet
check_user_remote_dir = 0

[GRID]
ce_white_list = ntugrid2.phys.ntu.edu.tw 
se_white_list = ntugrid4.phys.ntu.edu.tw 
virtual_organization = cms

  • Result:
Filter efficiency 62.545% +- 0.153%
# of events 102946
Dataset name /Bc_Bspi_8TeV_GenFilter_v1/ctao-Bc_Bspi_8TeV_GenFilter_v1-ef5ccca7d96eec95b1d60a9a6da4a974/USER

STEP2: DIGI + PU-mixing

  • command line options:
cmsDriver.py STEP2 --step DIGI,L1,DIGI2RAW,HLT:7E33v2 --python_filename BcBspi_PUMIX.py --fileout BcBspi_PUMIX.root --conditions START53_V7C::All --pileup 2012_Summer_50ns_PoissonOOTPU --datamix NODATAMIXER --eventcontent RAWSIM --datatier GEN-SIM-RAW --no_exec -n 10 
Add the section below in the configuration file:
process.mix.input.fileNames = cms.untracked.vstring(
'file:/wk1/kfjack/pileup_samples/summer12/00EEF21D-0762-E111-9A59-003048F0E1B2.root',
'file:/wk1/kfjack/pileup_samples/summer12/0220793B-9862-E111-A38C-003048CF6338.root',
'file:/wk1/kfjack/pileup_samples/summer12/02370007-2062-E111-98D4-0030487D5DC3.root',
'file:/wk1/kfjack/pileup_samples/summer12/048594F0-7E62-E111-ADCD-0030487F1797.root',
'file:/wk1/kfjack/pileup_samples/summer12/06B7EB97-6562-E111-AB9A-003048C693EA.root',
'file:/wk1/kfjack/pileup_samples/summer12/0A3D9816-6562-E111-95DD-0030487D8541.root',
'file:/wk1/kfjack/pileup_samples/summer12/0C25C4D8-0562-E111-8C61-003048D436CA.root',
'file:/wk1/kfjack/pileup_samples/summer12/0E699140-0362-E111-B364-003048D436CA.root',
'file:/wk1/kfjack/pileup_samples/summer12/0E9AAA91-1C62-E111-B296-003048D43700.root',
'file:/wk1/kfjack/pileup_samples/summer12/100F5F3E-0362-E111-8874-003048D436CA.root',
'file:/wk1/kfjack/pileup_samples/summer12/1034F28B-1B62-E111-A2DC-003048D43944.root',
'file:/wk1/kfjack/pileup_samples/summer12/10439A0E-9162-E111-9C71-002481E0D974.root',
'file:/wk1/kfjack/pileup_samples/summer12/10B43512-2262-E111-963E-003048C692B4.root',
'file:/wk1/kfjack/pileup_samples/summer12/1205F429-6262-E111-9F30-0030487D83B9.root',
'file:/wk1/kfjack/pileup_samples/summer12/160BE53C-6862-E111-AC84-003048D4365C.root',
'file:/wk1/kfjack/pileup_samples/summer12/16DCB5B2-5562-E111-AF22-003048C692DA.root',
'file:/wk1/kfjack/pileup_samples/summer12/188B3433-7D62-E111-A4C6-0030487F938F.root',
'file:/wk1/kfjack/pileup_samples/summer12/188C4639-1462-E111-BAD3-0030487F6A8B.root',
'file:/wk1/kfjack/pileup_samples/summer12/18BCBC23-7462-E111-83FC-002481E0D974.root',
'file:/wk1/kfjack/pileup_samples/summer12/1A4B2357-9E62-E111-9852-002481E15176.root',
'file:/wk1/kfjack/pileup_samples/summer12/1C8EF33E-1062-E111-8DE0-003048C662C8.root',
'file:/wk1/kfjack/pileup_samples/summer12/1CD51127-1D62-E111-905A-003048C69314.root',
'file:/wk1/kfjack/pileup_samples/summer12/1EB056A0-0962-E111-99E2-003048F02CB4.root',
'file:/wk1/kfjack/pileup_samples/summer12/223DAE52-0562-E111-AC14-003048CF6346.root',
'file:/wk1/kfjack/pileup_samples/summer12/224FA372-3362-E111-9583-003048F0E1AE.root',
'file:/wk1/kfjack/pileup_samples/summer12/22509669-8962-E111-836D-0030487F16FB.root',
'file:/wk1/kfjack/pileup_samples/summer12/2251AAAE-7562-E111-86DD-003048D479C0.root',
'file:/wk1/kfjack/pileup_samples/summer12/247A47A9-9462-E111-94C8-002481E946EE.root',
'file:/wk1/kfjack/pileup_samples/summer12/24A2C488-1F62-E111-A035-003048D462FE.root',
'file:/wk1/kfjack/pileup_samples/summer12/26474D1C-3E62-E111-AEAA-003048D4797E.root',
'file:/wk1/kfjack/pileup_samples/summer12/28144A44-5462-E111-8629-003048F0E3AE.root',
'file:/wk1/kfjack/pileup_samples/summer12/283BFB28-5B62-E111-BB68-003048C692D8.root',
'file:/wk1/kfjack/pileup_samples/summer12/2865B0E6-9A62-E111-8FC4-002481E15204.root',
'file:/wk1/kfjack/pileup_samples/summer12/2A16F88E-1B62-E111-8B75-003048D4DEA6.root',
'file:/wk1/kfjack/pileup_samples/summer12/2C055754-5662-E111-AA22-003048C69292.root',
'file:/wk1/kfjack/pileup_samples/summer12/2C55E61E-4862-E111-B608-0030487E4EC5.root',
'file:/wk1/kfjack/pileup_samples/summer12/2EB00C4C-2E62-E111-BA62-003048F02CB6.root',
'file:/wk1/kfjack/pileup_samples/summer12/303F20E6-2B62-E111-A742-003048C6763A.root',
'file:/wk1/kfjack/pileup_samples/summer12/3061CB45-6362-E111-B978-0030487D43E1.root',
'file:/wk1/kfjack/pileup_samples/summer12/3091E778-0C62-E111-A7FB-003048C662C8.root',
'file:/wk1/kfjack/pileup_samples/summer12/30AF0710-5A62-E111-A8AC-002481E0D6EE.root',
'file:/wk1/kfjack/pileup_samples/summer12/30E11B0D-0862-E111-A56B-003048F02CB4.root',
'file:/wk1/kfjack/pileup_samples/summer12/30EA9B2F-6B62-E111-BA50-0030487E4EB7.root',
'file:/wk1/kfjack/pileup_samples/summer12/323258DD-9862-E111-AAA6-003048CF6338.root',
'file:/wk1/kfjack/pileup_samples/summer12/32B47ABC-0B62-E111-9C9A-003048F0E1B2.root',
'file:/wk1/kfjack/pileup_samples/summer12/32EB056C-0D62-E111-850A-003048C662C8.root',
'file:/wk1/kfjack/pileup_samples/summer12/344A5882-1662-E111-92D6-0030487F6A8B.root',
'file:/wk1/kfjack/pileup_samples/summer12/349797D4-5F62-E111-817F-0030487D83B9.root',
'file:/wk1/kfjack/pileup_samples/summer12/34CB55FB-8C62-E111-B71E-003048C693B8.root',
'file:/wk1/kfjack/pileup_samples/summer12/3610A5F4-0762-E111-A013-003048F0E1B2.root'
)

  • Result:

STEP3: RECO

    • command line options:
cmsDriver.py STEP3 --step RAW2DIGI,L1Reco,RECO --python_filename BcBspi_RECO.py --fileout BcBspi_RECO.root --conditions START53_V7C::All --eventcontent RECOSIM --datatier GEN-SIM-RECO --no_exec -n 10
    • Result:
# of event: 102946 # of Bc at gen level: 20611 Have to regenerate MC with total # of event~800000
Edit | Attach | Watch | Print version | History: r10 < r9 < r8 < r7 < r6 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r10 - 2013-07-04 - ChiaoYuTao
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

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