EditWYSIWYGAttachPrintablePDF
r2 - 22 Dec 2011 - 06:33:39 - GunmoRyuYou are here: KoreaCmsWiki >  Top_SM Web > DoublyChargedHiggsToEMuAnalysis
Search for Pair produced H++/H-- decaying to electron-muon channel

Introduction

Monte-Carlo production

Signal samples

We produced signal samples using PYTHIA MC generator. The parameters are inherited from the default values of CMSSW configuration. In this MC samples, we assumed masses of HL and HR are degenerated. We consider pair production of H±± via Drell-Yan but not consider WW fusion. We turned on the decay to e-mu channel only to simplify the analysis. Thus the number of signal entries will be the maximum number of expected signals and the upper limit we get will be the conservative one.

  • To be understood
    • MSEL=0
    • CKIN(45)=5
    • MSTP(25)=2 : HZZ4f angular decay? why here?

The following code is the fragment of PYTHIA configuration file used in the MC generation.

More...Close

    PythiaParameters = cms.PSet(
        pythiaUESettingsBlock,
        processParameters = cms.vstring(
            'PMAS(353,1)=%s          !mass of H_L' % hppMass,
            'PMAS(354,1)=%s          !mass of H_R' % hppMass,
            'MSEL=0                  !(D=1) to select between full user control (0, then use MSUB) and some preprogrammed alternative: QCD hight pT processe
s (1, then ISUB=11, 12, 13, 28, 53, 68), QCD low pT processes (2, then ISUB=11, 12, 13, 28, 53, 68, 91, 92, 94, 95)',
            'MSUB(349)=1             !H_L H_L',
            'MSUB(350)=1             !H_R H_R',
            'MSUB(351)=0             !H_L ff (WW fusion)',
            'MSUB(352)=0             !H_R ff (WW fusion)',
            'CKIN(45)=5.             !high mass cut on m2 in 2 to 2 process Registered by Chris.Seez@cern.ch',
            'MSTP(25)=2              !Angular decay correlations in H->ZZ->4fermions Registered by Alexandre.Nikitenko@cern.ch',
            'PARP(181)= 0.0          !H++ coupling with ee',
            'PARP(182)= 0.1          !H++ coupling with emu',
            'PARP(183)= 0.0          !H++ coupling with etau',
            'PARP(184)= 0.0          !H++ coupling with mue',
            'PARP(185)= 0.0          !H++ coupling with mumu',
            'PARP(186)= 0.0          !H++ coupling with mutau',
            'PARP(187)= 0.0          !H++ coupling with taue',
            'PARP(188)= 0.0          !H++ coupling with taumu',
            'PARP(189)= 0.0          !H++ coupling with taumutau',
            'PARP(190)= 0.0          !H++_L coupling with W',
            'PARP(191)= 0.0          !H++_R coupling with W'
         ),        # This is a vector of ParameterSet names to be read, in this order
        parameterSets = cms.vstring('pythiaUESettings',
            'processParameters')
    )

The CMSSW Default is defined like this:

More...Close

            'MSTJ(11)=3     ! Choice of the fragmentation function',
            'MSTJ(22)=2     ! Decay those unstable particles',
            'PARJ(71)=10 .  ! for which ctau  10 mm',
            'MSTP(2)=1      ! which order running alphaS',
            'MSTP(33)=0     ! no K factors in hard cross sections',
            'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)',
            'MSTP(52)=2     ! work with LHAPDF',
            'MSTP(81)=1     ! multiple parton interactions 1 is Pythia default',
            'MSTP(82)=4     ! Defines the multi-parton model',
            'MSTU(21)=1     ! Check on possible errors during program execution',
            'PARP(82)=1.8387   ! pt cutoff for multiparton interactions',
            'PARP(89)=1960. ! sqrts for which PARP82 is set',
            'PARP(83)=0.5   ! Multiple interactions: matter distrbn parameter',
            'PARP(84)=0.4   ! Multiple interactions: matter distribution parameter',
            'PARP(90)=0.16  ! Multiple interactions: rescaling power',
            'PARP(67)=2.5    ! amount of initial-state radiation',
            'PARP(85)=1.0  ! gluon prod. mechanism in MI',
            'PARP(86)=1.0  ! gluon prod. mechanism in MI',
            'PARP(62)=1.25   ! ',
            'PARP(64)=0.2    ! ',
            'MSTP(91)=1      !',
            'PARP(91)=2.1   ! kt distribution',
            'PARP(93)=15.0  ! '

The configuration fragment files for various mass hypothesis are already uploaded in the CVS. 

cvs co -d HiggsAnalysis/DoublyChargedHiggs/python UserCode/JHGoh/HiggsAnalysis/DoublyChargedHiggs/python
# If you are using CMSSW_2_2_X versions, use the branch for 22X releases
cvs co -r B22X -d HiggsAnalysis/DoublyChargedHiggs/python UserCode/JHGoh/HiggsAnalysis/DoublyChargedHiggs/python

With this configuration, we produced MC samples at the KNU T2 or IN2P3? T2 using CRAB grid tool.

You can use the following bash script to generate CMSSW,CRAB configuration files.

More...Close

#!/bin/bash

Reconstruction and preselection

There are 2 electrons and 2 muons in the final state of this channel. This final state is already considered in the Higgs to ZZ to 4 leptons study. We collect electrons and muons using the PAT(Physics Analysis Toolkit) interface, which allows to remove

-- JunghwanGoh - 06 Jan 2010

Edit | WYSIWYG | Attach | Printable | Raw View | Backlinks: Web, All Webs | History: r2 < r1 | More topic actions
Top_SM.DoublyChargedHiggsToEMuAnalysis moved from Top_SM.DoublyChargedHiggsToEMuAnalysis on 22 Dec 2011 - 06:33 by GunmoRyu - put it back
 
KOREA-CMS
This site is powered by the TWiki collaboration platformCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding KoreaCmsWiki? Send feedback