ZXanalysis < Sandbox

Overview

Usefull Links

Instructions

Generating MC samples

Production Diffractive Z boson

Instructions for simulating diffractive Z bosons using Pythia8.2 or POMWIG MC generator using the 2017 standard run conditions. The code should be executed within the CMSSW software.

setup CMSSW

Hide setup

scram p CMSSW CMSSW_10_6_16
cd CMSSW_10_6_16/src
eval `scram r -sh`
scram b
cd -

Creating generator fragment: Create a new fragment file (SMP-RunII-00001-fragment.py) containing either Pythia or POMWIG generator settings.

pythia82-fragment.py

Hide fragment

import FWCore.ParameterSet.Config as cms
from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
generator = cms.EDFilter("Pythia8GeneratorFilter",
        comEnergy = cms.double(13000.0),
        crossSection = cms.untracked.double(45.5),
        filterEfficiency = cms.untracked.double(1),
        maxEventsToPrint = cms.untracked.int32(1),
        pythiaHepMCVerbosity = cms.untracked.bool(False),
        pythiaPylistVerbosity = cms.untracked.int32(1),
        PythiaParameters = cms.PSet(
               pythia8CommonSettingsBlock,
               pythia8CP5SettingsBlock,
               processParameters = cms.vstring(
                   'PhaseSpace:pTHatMin = 12.',
                   'WeakSingleBoson:ffbar2ffbar(s:gmZ) = on',
                   'Diffraction:doHard = on',
                   'Diffraction:sampleType = 4',
                   'SigmaDiffractive:mode = 2',
                   'SigmaDiffractive:PomFlux = 7', #The H1 Fit B parametrisation of the Pomeron flux
                   'PDF:PomSet = 4',                     # H1 2006 Fit B NLO
                   'PDF:pSet = 20',                        # NNPDF3.1 QCD+LUXQED NNLO
                   '23:onMode = off',                     # turn off all Z decays
                   '23:onIfAny = 11 13 15' ),          # turn ON Z->leplep
                parameterSets = cms.vstring('pythia8CommonSettings',
                                           'pythia8CP5Settings',
                                           'processParameters')
        )
)
ProductionFilterSequence = cms.Sequence(generator)

pomwig-fragment.py

Hide fragment

POMWIG is based on the Hadron Emission Reactions With Interfering Gluons (HERWIG) MC Generator. The CMSSW code is available here

import FWCore.ParameterSet.Config as cms
herwig6Parameters = cms.PSet(
        comEnergy = cms.double(13000.0),
        useJimmy = cms.bool(False),
        doMPInteraction = cms.bool(False),
        herwigHepMCVerbosity = cms.untracked.bool(False),
        herwigVerbosity = cms.untracked.int32(1),
        printCards = cms.untracked.bool(True),
        maxEventsToPrint = cms.untracked.int32(2),
        crossSection = cms.untracked.double(11.51),
        filterEfficiency = cms.untracked.double(1.0),
)
generator = cms.EDFilter("PomwigGeneratorFilter",
    herwig6Parameters,
    HerwigParameters = cms.PSet(
        parameterSets = cms.vstring('SDInclusiveWmunu'),
        SDInclusiveWmunu = cms.vstring('NSTRU      = 14         ! H1 Pomeron Fit B',
            'Q2WWMN     = 1E-6       ! Minimum |t|',
            'Q2WWMX     = 4.0         ! Maximum |t|',
            'EMMIN      = 50.0           ! minimum DY Mass',
            'YWWMIN     = 1E-2       ! Minimum xi',
            'YWWMAX     = 0.15       ! Maximum xi',
            'IPROC      = 11350         ! Process PomP -> Z -> Lep Lep ',
            'MODPDF(1)  = 10042    ! Set MODPDF CTEQ6L1',
            'MODPDF(2)  = -1           ! Set MODPDF')
    ),
    diffTopology = cms.int32(2),
    survivalProbability = cms.double(0.05),
    h1fit = cms.int32(2),
    doPDGConvert = cms.bool(False)
)
ProductionFilterSequence = cms.Sequence(generator)

Generating RAW file: The following sequence is used to generate RAW.root output w/o HLT as explained here.

cmsDriver.py SMP-RunII-00001-fragment.py --fileout file:stepRAW.root --mc --eventcontent PREMIXRAW --datatier GEN-SIM-RAW --step GEN,SIM,DIGI,DATAMIX,L1,DIGI2RAW,HLT:@fake2 --conditions 106X_mc2017_realistic_v6 --beamspot Realistic25ns13TeVEarly2017Collision --procModifiers premix_stage2 --datamix PreMix --era Run2_2017 --nThreads 8 --geometry DB:Extended --pileup_input "dbs:/Neutrino_E-10_gun/RunIISummer19ULPrePremix-UL17_106X_mc2017_realistic_v6-v1/PREMIX" --no_exec -n 100

Generating MINIAOD file:

cmsDriver.py step2 --filein file:stepRAW.root --fileout file:MiniAODSIM.root --mc --eventcontent MINIAODSIM --datatier MINIAODSIM --runUnscheduled --step RAW2DIGI,L1Reco,RECO,RECOSIM,PAT --conditions 106X_mc2017_realistic_v6  --nThreads 8 --geometry DB:Extended --era Run2_2017 --no_exec -n 100

In this example, 100 events will be generated. To change the final number of events, modify the -n 100 flag.

Production Exclusive ZZ and Zgamma pair production

First one should create Madgraph gridpacks. The most recent recommendations are available in QuickGuideMadGraph5aMCatNLO#Quick_tutorial_on_how_to_produce

To create gridpacks, execute the following lines:

git clone https://github.com/cms-sw/genproductions.git
cd genproductions/bin/MadGraph5_aMCatNLO/
cp -r ~mpitt/public/CEP/Madgraph/* cards/
./gridpack_generation.sh  ExclZZ_box cards/ExclZZ_box

For job submission in condor, use code in ~mpitt/public/TopWidth/condor

Analysis

run this and that

Plots

execute this and that

Topic revision: r6 - 2021-06-12 - MichaelPitt

Account
- Log In

- Cern Search
- TWiki Search
- Google Search
Sandbox All webs

Edit
Attach

Copyright &© 2008-2023 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