TWiki> Main Web>VBFWWJetsMCGeneration (2012-07-14, unknown) EditAttachPDF

Private Production for VBF MC WW Background Samples with Summer12 Configurations

People

Qiang Li, Wei Zou, Zijun Xu (Peking University)

STEP 1: MC LHE Files 

    some description about using the MadGraph/MadEvent to generate the LHE files

Sample description # of evts submitted/finished Dataset Name Responsible person Site
VBF W+ W- > l nu j j 339359 WmWpEW_v1_unweighted_events.lhe.gz Zijun Xu ftp://162.105.246.56/pub/CMS/xuzj/VBF/WmWp_EW/
VBF W Z > l nu j j 400000 WZEW_v2_unweighted_events.lhe.gz Zijun Xu ftp://hepfarm02.phy.pku.edu.cn/pub/CMS/xuzj/VBF/WZ_EW/
QCD W+ W- plus 1Jet, W+ W- > l nu j j 500000 WWplus1Jet_v2_unweighted_events.lhe.gz Zijun Xu ftp://hepfarm02.phy.pku.edu.cn/pub/CMS/xuzj/VBF/WWplus1Jet_QCD/
QCD W Z plus 1Jet, W Z > l nu j j 500000 WZplus1Jet_QCD_v2_unweighted_events.lhe.gz Zijun Xu ftp://hepfarm02.phy.pku.edu.cn/pub/CMS/xuzj/VBF/WZplus1Jet_QCD/
QCD W+ W- plus 2Jets, W+ W- > l nu j j 135253（is reproducing now） WWplus2Jet_QCD_v1_unweighted_events.lhe.gz Zijun Xu ftp://hepfarm02.phy.pku.edu.cn/pub/CMS/xuzj/VBF/WWplus2Jets_QCD/
QCD W Z plus 2Jets, W Z > l nu j j 400000 WZplus2Jet_QCD_v1_unweighted_events.lhe.gz Zijun Xu ftp://hepfarm02.phy.pku.edu.cn/pub/CMS/xuzj/VBF/WZplus2Jets_QCD/
QCD W plus 4Jets, W > l nu 438082 Wplus4Jets_QCD_v1_unweighted_events.lhe.gz Zijun Xu ftp://hepfarm02.phy.pku.edu.cn/pub/CMS/xuzj/VBF/Wplus4Jets_QCD

STEP 2: GEN

Prepare the environment at LPC 

(cd to your working area at LPC)

setenv SCRAM_ARCH slc5_amd64_gcc462
source /uscmst1/prod/sw/cms/cshrc prod
cmsrel CMSSW_5_2_4
cd CMSSW_5_2_4/src
cmsenv
source /uscmst1/prod/grid/CRAB/crab.csh
source /uscmst1/prod/grid/gLite_SL5.csh
voms-proxy-init

Create the GEN-SIM python configuration (please replace crab_vbf_WpWmEW with a different name, change the lhe files to the corresponding files): 

mkdir crab_vbf_WpWmEW
cd crab_vbf_WpWmEW
cmsDriver.py -s GEN,SIM --beamspot Realistic8TeVCollision --conditions START52_V9::All --pileup NoPileUp --datamix NODATAMIXER --eventcontent RAWSIM --datatier GEN-SIM -n 10 --no_exec --mc --fileout zprime_GENSIM.root --filein run_01_unweighted_events.lhe --filetype LHE --python_filename qiangli-WpWmEW_GENSIM_cfg.py

Manually add the madgraph-pythia matching code in the qiangli-WpWmEW_GENSIM_cfg.py 

import FWCore.ParameterSet.Config as cms
from Configuration.Generator.PythiaUEZ2starSettings_cfi import *

generator = cms.EDFilter("Pythia6HadronizerFilter",
	pythiaHepMCVerbosity = cms.untracked.bool(True),
	maxEventsToPrint = cms.untracked.int32(0),
	pythiaPylistVerbosity = cms.untracked.int32(1),
	comEnergy = cms.double(8000.0),
	UseExternalGenerators = cms.untracked.bool(True),
	PythiaParameters = cms.PSet(
    	pythiaUESettingsBlock,
    	processParameters = cms.vstring('MSEL=0     	! User defined processes',
                    	'PMAS(5,1)=4.8   ! b quark mass',
                    	'PMAS(6,1)=172.5 ! t quark mass',
     	'MSTJ(1)=1   	! Fragmentation/hadronization on or off',
     	'MSTP(61)=1  	! Parton showering on or off'),
    	# This is a vector of ParameterSet names to be read, in this order
    	parameterSets = cms.vstring('pythiaUESettings',
        	'processParameters')
	),
)

Create a cb_qiangli-WpWmEW_GENSIM.cfg with the contents below: 

[CRAB]
jobtype  = CMSSW
scheduler = condor

[CMSSW]
datasetpath = None
get_edm_output = 1
first_run = 1
total_number_of_events = 550000
events_per_job = 500
pset = qiangli-WpWmEW_GENSIM_cfg.py 
generator = lhe

[USER]
return_data = 0
copy_data = 1
publish_data = 1
storage_element = cmssrm.fnal.gov
storage_path = /srm/managerv2?SFN=11
user_remote_dir = /store/user/lnujj
dbs_url_for_publication =https://cmsdbsprod.cern.ch:8443/cms_dbs_ph_analysis_01_writer/servlet/DBSServlet
publish_data_name = qiangli-WpWmEW_5_2_4_SIM_Full_NEW_SEED
ui_working_dir = qiangli-WpWmEW_5_2_4_SIM_Full_NewSeed
check_user_remote_dir   = 0
additional_input_files = input01.lhe, input02.lhe, input03.lhe, input04.lhe, input05.lhe, input06.lhe, input07.lhe, input08.lhe, input09.lhe, input10.lhe, input11.lhe

And then you need replace the source lhe file name in your qiangli-WpWmEW_GENSIM_cfg.py : 

process.source = cms.Source("LHESource",
    fileNames = cms.untracked.vstring(
   'file:input01.lhe',
   'file:input02.lhe',
   'file:input03.lhe',
   'file:input04.lhe',
   'file:input05.lhe',
   'file:input06.lhe',
   'file:input07.lhe',
   'file:input08.lhe',
   'file:input09.lhe',
   'file:input10.lhe',
   'file:input11.lhe'
   )   
)

And there is a simple script to rename the files name to the desired format: 

ls genWpWmEW/WpWmEW_* | awk 'BEGIN {i = 0} {if(/genWpWmEW/) {i =i + 1;if(i < 10) {print "ln " $0" input0"i".lhe "};if(i >= 10){print "ln " $0" input"i".lhe "} }}

Then, just create the job and submit: 

crab -create
crab -submit

Once the job is done, please publish it to db: 

crab -getoutput
crab -report
crab -publish

Here are the table for book keeping:

STEP2 Summary Table

Sample description # of evts submitted/finished Dataset Name Responsible person Site

STEP 3: SIM

Make sure you are still in the crab_vbf_WpWmEW working directory.

Create the python configuration for "step2" (DIGI,L1,DIGI2RAW,HLT:GRun), again replace the crab_vbf_WpWmEW with a different directory name: 

cmsDriver.py REDIGI --step DIGI,L1,DIGI2RAW,HLT:GRun --conditions START50_V15::All --pileup 2012_Startup_50ns_PoissonOOTPU --customise Configuration/StandardSequences/DigiToRecoNoPU.py --datamix NODATAMIXER --eventcontent RAWSIM --datatier GEN-SIM-RAW --no_exec --mc --filein tc_GENSIM.root --fileout tc_RAWSIM.root --python_filename qiangli-WpWmEW_RAWSIM_cfg.py

Manually add the pile-up samples into the configuration (insert the following block into the qiangli-WpWmEW_RAWSIM_cfg.py): 

process.mix.input.fileNames = cms.untracked.vstring(
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/00024E74-D96A-E111-9CE8-0030487D5DC3.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/003F446E-DC67-E111-A0F9-00266CFFA7A8.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/00DDA010-A769-E111-8A12-0030487D5EAF.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/02124719-B769-E111-8971-003048C66BBE.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/0261C34F-C969-E111-B0E6-003048C692E2.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/026B25B2-D967-E111-8561-00266CFFA5E0.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/0289AEA8-6165-E111-8512-0030487D5EB1.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/029CC6E3-9967-E111-95B2-00266CF32684.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/02C1D0AB-F267-E111-B434-003048F0E812.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/02E93ED0-6765-E111-885A-0030487F1F23.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/02F1619A-6E69-E111-A950-0030487E4EC7.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/0413ADCF-DD69-E111-85EC-003048C69414.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/04412411-4468-E111-9826-003048C66BBE.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/044C8BDF-516A-E111-AFFF-003048D436D2.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/044ED2FA-286B-E111-83D9-0030487EAFF9.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/04D7763A-2A6B-E111-A4C4-003048C66184.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/0646CB95-D567-E111-ACED-0025901D4C3C.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/067CF4F5-1568-E111-8B12-00266CF33288.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/06932BAC-BF69-E111-AE28-002481E14F2A.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0002/06972F6F-116B-E111-BC1A-0030487E0A2D.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0005E496-3661-E111-B31E-003048F0E426.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/003EEBD4-8061-E111-9A23-003048D437F2.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/005825F1-F260-E111-BD97-003048C692DA.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0065594C-B35E-E111-8B8C-003048C693EA.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0091FFD0-6B5E-E111-92FE-003048C693DA.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/00C69AE3-FE60-E111-BC48-0030487D8633.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02079692-AC61-E111-97BB-0025901D4D54.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/021BD915-2D61-E111-8BAD-002481E76052.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02386E4D-DC5E-E111-9413-00266CF1074C.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02446A08-515E-E111-82C7-00266CF330B8.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/024BBC60-CD5D-E111-A41A-0025901D4C32.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/025AB53D-F15E-E111-A7F4-003048C66BBE.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0274007A-A05E-E111-9227-003048D4DFB8.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02AC6853-E65D-E111-A515-0025901D4B04.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02B17485-B861-E111-84FC-003048C6930E.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/02CFF160-B85E-E111-B547-003048D3CD92.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0407D67D-F360-E111-A737-0030487F1A73.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0407DB67-2F5E-E111-8115-003048D47976.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/0464DA8A-4D61-E111-95F5-003048C64787.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0000/04670075-3561-E111-AD81-002481E0D50C.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/0015C0BD-0C64-E111-81F8-003048D43942.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/0060F354-3364-E111-A5C7-003048D43960.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/00852170-E061-E111-BA8F-003048D3CDE0.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/00B6AAA8-0D65-E111-966C-0025901D493A.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/00EE834E-3464-E111-9B62-003048D43960.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/00EEF21D-0762-E111-9A59-003048F0E1B2.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02163E93-0A65-E111-B2DE-00266CF32E78.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/021CC150-0A65-E111-90F2-00266CF33340.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/0220793B-9862-E111-A38C-003048CF6338.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02370007-2062-E111-98D4-0030487D5DC3.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/025EC033-EC63-E111-A9D4-0025901D4B22.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/026FA1DC-5F64-E111-8D64-003048D439B4.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/029C9303-0A65-E111-B8DE-0025901D4B22.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02B30A54-C561-E111-AE69-003048D43700.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02B60BCF-DF63-E111-906F-003048C693E6.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02BC9C04-BD64-E111-9B60-003048D436CA.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02CE813D-0A65-E111-9721-0025901D4A58.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/02D8803D-EE63-E111-ABE4-0030487F1657.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/0400B93B-E264-E111-B99D-003048C6930E.root',
'/store/mc/Summer12/MinBias_TuneZ2star_8TeV-pythia6/GEN-SIM/START50_V13-v3/0001/04067303-0A65-E111-BAD0-0025901D4ADE.root',
)

Create a cb_qiangli-WpWmEW_RAWSIM.cfg with the contents below: 

[CRAB]
jobtype  = CMSSW
scheduler = condor

[CMSSW]
datasetpath = /qiangli-WpWmEW_4_2_8_SIM_Full_NEW_SEED/weizou-qiangli-WpWmEW_4_2_8_SIM_Full_NEW_SEED-53a21e6f3c7b0deae641bf4fdf3b97bb/USER 
get_edm_output = 1
first_run = 1
total_number_of_events = -1
number_of_jobs = 500
pset = qiangli-WpWmEW_RAWSIM_cfg.py 

[USER]
return_data = 0
copy_data = 1
publish_data = 1
storage_element = cmssrm.fnal.gov
storage_path = /srm/managerv2?SFN=11
user_remote_dir = /store/user/lnujj
dbs_url_for_publication =https://cmsdbsprod.cern.ch:8443/cms_dbs_ph_analysis_01_writer/servlet/DBSServlet
publish_data_name = qiangli-WpWmEW_5_2_4_RAW_FULL_NEW_SEED
check_user_remote_dir   = 0
ui_working_dir = qiangli-WpWmEW_5_2_4_FULL_RAW_Seed

Create the job and submit: 

crab -create
crab -submit

Once the job is done, please publish it to db: 

crab -getoutput
crab -report
crab -publishNoInp

Here are the table for book keeping:

STEP3 Summary Table

Sample description # of evts submitted/finished Dataset Name Responsible person Site

STEP 4: RECO

Make sure you are still in the crab_vbf_WpWmEW working directory.

Create the python configuration for "step3" (RAW2DIGI,L1Reco,RECO), again replace the crab_vbf_WpWmEW with a different directory name: 

cmsDriver.py STEP2 --step RAW2DIGI,L1Reco,RECO,VALIDATION:validation_prod,DQM:DQMOfflinePOGMC --conditions START50_V15::All --pileup 2012_Startup_50ns_PoissonOOTPU --datamix NODATAMIXER --eventcontent AODSIM --datatier AODSIM -n 10 --no_exec --mc --filein tc_RAWSIM.root --fileout tc_AODSIM.root --python_filename qiangli-WpWmEW_RECO_cfg.py

Create a cb_qiangli-WpWmEW_RECO.cfg with the contents below: 

[CRAB]
jobtype  = CMSSW
scheduler = condor

[CMSSW]
datasetpath = /qiangli-WpWmEW_4_2_8_SIM_Full_NEW_SEED/weizou-qiangli-WpWmEW_4_2_8_RAW_FULL_NEW_SEED-fbfd28b452cfd6916cfaf24ef9bea208/USER 
dbs_url = http://cmsdbsprod.cern.ch/cms_dbs_ph_analysis_01/servlet/DBSServlet
get_edm_output = 1
first_run = 1
total_number_of_events = -1
number_of_jobs = 300
pset = qiangli-WpWmEW_RECO_cfg.py 

[USER]
return_data = 0
copy_data = 1
publish_data = 1
storage_element = cmssrm.fnal.gov
storage_path = /srm/managerv2?SFN=11
user_remote_dir = /store/user/lnujj
dbs_url_for_publication =https://cmsdbsprod.cern.ch:8443/cms_dbs_ph_analysis_01_writer/servlet/DBSServlet
publish_data_name = qiangli-WpWmEW_5_2_4_AODSIM_FULL_NEW_SEED
check_user_remote_dir   = 0
ui_working_dir = qiangli-WpWmEW_5_2_4_AODSIM_FULL_New_Seed

Create the job and submit: 

crab -create
crab -submit

Once the job is done, please publish it to db: 

crab -getoutput
crab -report
crab -publishNoInp

Here are the table for book keeping:

STEP4 Summary Table

Sample description # of evts submitted/finished Dataset Name Responsible person Site

-- WeiZou - 24-May-2012

Edit | Attach | Watch | Print version | History: r10 < r9 < r8 < r7 < r6 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r10 - 2012-07-14 - unknown
 
    • 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-2024 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