Dilepton PInG MC requests

Requests

Quarkonia_QM_Feb2011

Signals to embed: PYTHIA+EvtGen(PHOTOS), combined with a Selector

• prompt JPsi->mumu
  • cfg: PyquenEvtGen_jpsiMuMu_cfi.py

• non-prompt JPsi->mumu (B->JPsi)
  • cfg: PyquenEvtGen_bJpsiMuMu_cfi.py

• Upsilon->mumu
  • cfg: PyquenEvtGen_upsilon1sMuMu_cfi.py

Technical (bkg, cmssw, tag, datatier, eventcontent)

• Background: HYDJET, HI MC sample 2010, tune BASS (hydjetSetting with pyquenSettings)
• CMSSW: 3_11_1 (problem solving incompatibilities with smaller versions, like 3_9_9_patch1)
• TAG: START311_V1::All (for compatibility with START39_V7HI::All, that has a different BS, we overwrite on-the-fly the BS)
• Datatier: GEN-SIM-DIGI-RAW-RECO
• Eventcontent: FEVTDEBUGHLT
• Note: HLT:HIon is not available in the GT used (the HLT will have to be run separately)

Details about the settings:

• Need to

           cvs co -r HEAD GeneratorInterface/Core
           cvs co -r HEAD GeneratorInterface/HiGenCommon
           cvs co -r HEAD GeneratorInterface/PyquenInterface
           cvs co -r HEAD Configuration/GenProduction/python
   

• Need to add in the mixing file:
  • • 1.

       process.hiGenParticles.abortOnUnknownPDGCode = cms.untracked.bool(False)  

    • 2.

             def overrideBeamSpot(process):
             process.GlobalTag.toGet = cms.VPSet(
                          cms.PSet(record = cms.string("BeamSpotObjectsRcd"),
                                          tag = cms.string("Realistic2.76ATeVCollisions_STARTUP_v0_mc"),
                                          connect = cms.untracked.string("frontier://FrontierProd/CMS_COND_31X_BEAMSPOT")
                                            ),
                                          )
          return process
           

• • • 3.

    cms.PSet(
        record  = cms.string( "L1GtTriggerMenuRcd" ),
        tag     = cms.string( "L1GtTriggerMenu_L1Menu_CollisionsHeavyIons2010_v2_mc" ),
        connect = cms.untracked.string( process.GlobalTag.connect.value().replace('CMS_COND_31X_GLOBALTAG', 'CMS_COND_31X_L1T') )
      )

• To change the pT range, the only thing that needs to be done is to modify the value for hadronPtMax and hadronPtMin

  p_T JPsi bin	[0,3]	[3,6]	[6,9]	[9,12]	[12,15]	[15,30]
  setting	hadronPtMax = cms.vdouble(3.)	hadronPtMax = cms.vdouble(6.)	hadronPtMax = cms.vdouble(9.)	hadronPtMax = cms.vdouble(12.)	hadronPtMax = cms.vdouble(15.)	hadronPtMax = cms.vdouble(30.)
  setting	hadronPtMin = cms.vdouble(0.)	hadronPtMin = cms.vdouble(3.)	hadronPtMin = cms.vdouble(6.)	hadronPtMin = cms.vdouble(9.)	hadronPtMin = cms.vdouble(12.)	hadronPtMin = cms.vdouble(15.)

• cmsDrive.py commands:

  • cmsDriver.py Configuration/GenProduction/python/PyquenEvtGen_jpsiMuMu_cfi.py -n 2000 -s GEN,SIM,DIGI,L1,DIGI2RAW,HLT:HIon,RAW2DIGI,RECO --scenario HeavyIons --himix --geometry Extended --conditions START311_V1::All --datatier GEN-SIM-DIGI-RAW-RECO --eventcontent=FEVTDEBUGHLT --filein=XXXX --fileout=Hydjet_MinBias_2760GeV_pt6_GEN-SIM-RAW-HLT-RECO_EmbeddingJpsiMuMu_ptJPsi03.root --processName HISIGNAL --python_filename Hydjet_MinBias_2760_pt6_EmbeddingJpsiMuMu_ptJPsi03_GEN-SIM-L1-DIGI-RAW-HLT-RECO.py --no_exec 

  •  cmsDriver.py Configuration/GenProduction/python/PyquenEvtGen_bJpsiMuMu_cfi.py -n 2000 -s GEN,SIM,DIGI,L1,DIGI2RAW,HLT:HIon,RAW2DIGI,RECO --scenario HeavyIons --himix --geometry Extended --conditions START311_V1::All --datatier GEN-SIM-DIGI-RAW-RECO --eventcontent=FEVTDEBUGHLT --filein=XXXX --fileout=Hydjet_MinBias_2760GeV_pt6_GEN-SIM-RAW-HLT-RECO_EmbeddingBJPsiMuMu_ptJPsi03.root --processName HISIGNAL --python_filename Hydjet_MinBias_2760_pt6_EmbeddingBJPsiMuMu_ptJPsi03_GEN-SIM-L1-DIGI-RAW-HLT-RECO.py --no_exec 

  • cmsDriver.py Configuration/GenProduction/python/PyquenEvtGen_upsilon1sMuMu_cfi.py -n 2000 -s GEN,SIM,DIGI,L1,DIGI2RAW,HLT:HIon,RAW2DIGI,RECO --scenario HeavyIons --himix --geometry Extended --conditions START311_V1::All --datatier GEN-SIM-DIGI-RAW-HLTDEBUG-RECO --eventcontent=FEVTDEBUGHLT --filein=XXXX --fileout=Hydjet_MinBias_2760GeV_pt6_GEN-SIM-RAW-HLT-RECO_EmbeddingUpsilonMuMu_ptUps03.root --processName HISIGNAL --python_filename Hydjet_MinBias_2760_pt6_EmbeddingUpsilonMuMu_ptUps03_GEN-SIM-L1-DIGI-RAW-HLT-RECO.py --no_exec 

W embedding

• Request 20k of W+ and 20k of W- embedded in HYDJET, flat gun Wpt [0,30].

• Need to commit
  • a new version of PyquenDefaultSettings_cff.py that includes WtoMuNu decays : updated PyquenDefaultSettings_cff.py
  • 2 cfi for the generation of flat W- and W+
    • WMinusMuNu_FlatPt0to30_cfi.py
    • WPlusMuNu_FlatPt0to30_cfi.py

• cmssw 399patch1
• global tag : START39_V7HI::All
• Eventcontent: FEVTDEBUGHLT
• Background: HYDJET, HI MC sample 2010, tune BASS hydjetSetting
• Datatier: GEN-SIM-DIGI-RAW-RECO
• cmsDriver command

  Configuration/GenProduction/python/WMinusMuNu_FlatPt0to30_cfi.py -n 2 -s GEN:hiSignal,SIM,DIGI,L1,DIGI2RAW,HLT:HIon,RAW2DIGI,RECO --scenario HeavyIons --himix --geometry DB --conditions  START39_V7HI::All --datatier GEN-SIM-DIGI-RAW-RECO --eventcontent=FEVTDEBUGHLT --filein=rfio:/castor/cern.ch/user/d/dmoon/cms394/HydjetMB_Reco/Z0_2.76Tev_Hydjet_MB_Reco_90.root --fileout=file:/tmp/silvest/Hydjet_MinBias_2760GeV_pt6_GEN-SIM-RAW-RECO_EmbeddingWMinusMuNu.root --processName HISIGNAL --python_filename Hydjet_MinBias_2760_pt6_EmbeddingWMinusMuNu_GEN-SIM-L1-DIGI-RAW-RECO.py --no_exec 

Trigger studies

• Signals to embed:
  • 5K Z0
  • 5K Jpsi,
  • 5K Upsilon
• Background: HYDJET minbias

• CMSSW: 3_9_0_pre5
• TAG: START39_V1::All
• Datatier: GEN-SIM-DIGI-RAW-HLTDEBUG-RECO
• Eventcontent: FEVTDEBUGHLT
• driver command:

  • cmsDriver.py Configuration/GenProduction/python/SingleJPsimumu_FlatPt0to20_cfi.py -n 5000 -s GEN:hiSignal,SIM,DIGI,L1,DIGI2RAW,HLT:HIon,RAW2DIGI,RECO --scenario HeavyIons --himix --geometry DB --conditions START39_V1::All --datatier GEN-SIM-DIGI-RAW-HLTDEBUG-RECO --eventcontent=FEVTDEBUGHLT --filein=XXXX --fileout=Hydjet_MinBias_2760GeV_pt6_GEN-SIM-RAW-HLT-RECO_EmbeddingJpsiMuMu.root --processName HISIGNAL --python_filename Hydjet_MinBias_2760_pt6_EmbeddingJpsiMuMu_GEN-SIM-L1-DIGI-RAW-HLT-RECO.py --no_exec 

  • cmsDriver.py Configuration/GenProduction/python/SingleUpsilonmumu_FlatPt0to20_cfi.py -n 5000 -s GEN:hiSignal,SIM,DIGI,L1,DIGI2RAW,HLT:HIon,RAW2DIGI,RECO --scenario HeavyIons --himix --geometry DB --conditions START39_V1::All --datatier GEN-SIM-DIGI-RAW-HLTDEBUG-RECO --eventcontent=FEVTDEBUGHLT --filein=XXXX --fileout=Hydjet_MinBias_2760GeV_pt6_GEN-SIM-RAW-HLT-RECO_EmbeddingUpsilonMuMu.root --processName HISIGNAL --python_filename Hydjet_MinBias_2760_pt6_EmbeddingUpsilonMuMu_GEN-SIM-L1-DIGI-RAW-HLT-RECO.py --no_exec 

  •  cmsDriver.py Configuration/GenProduction/python/SingleZmumu_FlatPt0to20_cfi.py -n 5000 -s GEN:hiSignal,SIM,DIGI,L1,DIGI2RAW,HLT:HIon,RAW2DIGI,RECO --scenario HeavyIons --himix --geometry DB --conditions START39_V1::All --datatier GEN-SIM-DIGI-RAW-HLTDEBUG-RECO --eventcontent=FEVTDEBUGHLT --filein=XXXX --fileout=Hydjet_MinBias_2760GeV_pt6_GEN-SIM-RAW-HLT-RECO_EmbeddingZMuMu.root --processName HISIGNAL --python_filename Hydjet_MinBias_2760_pt6_EmbeddingZMuMu_GEN-SIM-L1-DIGI-RAW-HLT-RECO.py --no_exec 

Efficiency studies

• Signals to embed:
  • 10K Z0
  • 10K Jpsi,
  • 10K Upsilon
• Background: HYDJET minbias

• CMSSW: 3_X_X
• TAG: XXX
• Datatier: GEN-SIM-DIGI-RAW-RECO
• Eventcontent: FEVTDEBUGHLT

• Validation plots
  • usingembedConfigs.pdf: embedConfigs.pdf

• Configuration files ( *HEAD version* )
  • Settings: flat pT and Y (Pythia6PtYDistGun):
    • 0->20 GeV/c
    • |y|<2.5
  • python files:
    • JPsi: Configuration/GenProduction/python/SingleJPsimumu_FlatPt0to20_cfi.py
    • Upsilon: Configuration/GenProduction/python/SingleUpsilonmumu_FlatPt0to20_cfi.py
    • Z0: Configuration/GenProduction/python/SingleZmumu_FlatPt0to20_cfi.py

-- CameliaMironov - 30-Sep-2010