Producing a Sample

to be edited!

** don't forget to include RECO-SIM ** validating releases:

https://cms-pdmv.cern.ch/valdb/

Producer

The Producer code produces a single muon gun sample through a custom EDProducer and crab for use in muon alignement. It requires updating with new releases. Here, I will describe how to make a sample and update it for a new release.

Makeing a sample:

To create a new sample, edit these locations:

crabConfig.py

Name your event here:

config.General.requestName = 'singleMuonGun_10_X_13TeV_realistic_45_phase1_2018_design' #name here'

Naming convention: singleMuonGun + release + type of sample + number of events + Global Tag or other conditions

For type of sample, realistic = simulating pt and eta distribution of muons in current run. Could also be a single pt e.g. "50GeV"

Number of events "totalUnits":

config.Data.unitsPerJob = 500 NJOBS = 9000 config.Data.totalUnits = config.Data.unitsPerJob * NJOBS

output location:

config.Data.outLFNDirBase = '/store/group/alca_muonalign/'

Config file to run:

config.JobType.psetName = 'singleMuonGun_10X.py'

In that config file 'singleMuonGun_10X', edit these lines:

generator controls the EDProcuer "SingleMuonGun":

process.generator = cms.EDProducer("SingleMuonGun", Verbosity = cms.untracked.int32(30), # 1 - print begin and end of event # 10 - print type of gun (constant or spectrum) # 20 - print muon parameters (q, pt, eta, phi) # 30 - print CMSSW event info

IMPORTANT! *****************************************************************

ConstPt_eq_MinPt = cms.bool(False), # if TRUE then generate muons with CONSTANT pT = MinPt # if FALSE then generate muons with pT spectrum of muons as in 2017 data (hardcoded in SingleMuonGun/plugins/SingleMuonGun.cc)

****************************************************************************

MinPt = cms.double(1000.0), MaxPt = cms.double(200.0), MinEta = cms.double(-2.5), MaxEta = cms.double(2.5), MinPhi = cms.double(-3.14159265359), MaxPhi = cms.double(3.14159265359) )

This controls the GT:

process.GlobalTag = GlobalTag(process.GlobalTag, '100X_upgrade2018_design_IdealBS_v7', '')

check autoCond for updated GTs https://github.com/cms-sw/cmssw/blob/CMSSW_10_0_X/Configuration/AlCa/python/autoCond.py

Once you've updated all these things, you are ready to submit. Just one last thing: check that the group space has enough free space! eos quota | grep -B 0 -A 5 "muonalign"

cmsDriver.py SingleMuPt10_pythia8_cfi --conditions auto:phase1_2017_design -n 10 --eventcontent RECOSIM -s GEN,SIM,DIGI,L1,DIGI2RAW,RAW2DIGI,RECO --datatier GEN-SIM-RAW-RECO --era Run2_2017 --beamspot GaussSigmaZ4cm --geometry DB:Extended --fileout file:step1.root --no_exec

cmsDriver.py SingleMu_cfi --conditions auto:phase1_2018_design -n 10 --eventcontent RECOSIM -s GEN,SIM,DIGI,L1,DIGI2RAW,RAW2DIGI,RECO --datatier GEN-SIM-RAW-RECO --era Run2_2018 --beamspot GaussSigmaZ4cm --geometry DB:Extended --fileout file:step1.root --no_exec