Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

XeXe collision era and relval workflow #20749

Merged
merged 18 commits into from Oct 9, 2017
Merged
Show file tree
Hide file tree
Changes from 16 commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
Expand Up @@ -26,7 +26,8 @@
from RecoLocalCalo.EcalRecProducers.ecalGlobalUncalibRecHit_cfi import *
ecalUncalibRecHitSequence53X = cms.Sequence(ecalGlobalUncalibRecHit * ecalDetIdToBeRecovered)

uncalibRecHitSeq = cms.Sequence(ecalUncalibRecHitSequence, cms.Task(ecalDigis, ecalPreshowerDigis))
ecalAndPreshawerDigisForUncalibRecHitSeqTask = cms.Task(ecalDigis, ecalPreshowerDigis)
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

ecalAndPreshawerDigis -> ecalAndPreshowerDigis

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@perrotta Since this doesn't change the functionality, I will let the ongoing test finish rather than push a fix for this typo now. As a reviewer, I believe you can fix it directly if you'd like. Otherwise, I will wait till morning.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@mandrenguyen : I can't fix myself this PR. In any case, it is fine for me waiting for the test to finish, of course.

uncalibRecHitSeq = cms.Sequence(ecalUncalibRecHitSequence, ecalAndPreshawerDigisForUncalibRecHitSeqTask)

ALCARECOEcalUncalElectronECALSeq = cms.Sequence( uncalibRecHitSeq )

Expand Down
15 changes: 14 additions & 1 deletion Configuration/DataProcessing/python/Impl/pp.py
Expand Up @@ -19,6 +19,7 @@ def __init__(self):
self.recoSeq=''
self.cbSc='pp'
self.addEI=True
self.isRepacked=False
self.promptCustoms= [ 'Configuration/DataProcessing/RecoTLR.customisePrompt' ]
self.expressCustoms=[ ]
self.alcaHarvCustoms=[]
Expand All @@ -33,7 +34,10 @@ def __init__(self):

"""


def _setRepackedFlag(self,args):
if not 'repacked' in args:
args['repacked']= True

def promptReco(self, globalTag, **args):
"""
_promptReco_
Expand All @@ -50,6 +54,9 @@ def promptReco(self, globalTag, **args):
for c in self.promptCustoms:
args['customs'].append(c)

if self.isRepacked:
self._setRepackedFlag(args)

process = Reco.promptReco(self,globalTag, **args)

return process
Expand All @@ -70,6 +77,9 @@ def expressProcessing(self, globalTag, **args):
for c in self.expressCustoms:
args['customs'].append(c)

if self.isRepacked:
self._setRepackedFlag(args)

process = Reco.expressProcessing(self,globalTag, **args)

return process
Expand All @@ -86,6 +96,9 @@ def visualizationProcessing(self, globalTag, **args):

for c in self.visCustoms:
args['customs'].append(c)

if self.isRepacked:
self._setRepackedFlag(args)

process = Reco.visualizationProcessing(self,globalTag, **args)

Expand Down
@@ -0,0 +1,35 @@
#!/usr/bin/env python
"""
_ppEra_Run2_2017_pp_on_XeXe_

Scenario supporting proton collisions

"""

import os
import sys

from Configuration.DataProcessing.Reco import Reco
import FWCore.ParameterSet.Config as cms
from Configuration.Eras.Era_Run2_2017_pp_on_XeXe_cff import Run2_2017_pp_on_XeXe

from Configuration.DataProcessing.Impl.pp import pp

class ppEra_Run2_2017_pp_on_XeXe(pp):
def __init__(self):
pp.__init__(self)
self.recoSeq=''
self.cbSc='pp'
self.addEI=True
self.isRepacked=True
self.eras=Run2_2017_pp_on_XeXe
self.promptCustoms += [ 'Configuration/DataProcessing/RecoTLR.customisePostEra_Run2_2017_pp_on_XeXe' ]
self.expressCustoms += [ 'Configuration/DataProcessing/RecoTLR.customisePostEra_Run2_2017_pp_on_XeXe' ]
self.visCustoms += [ 'Configuration/DataProcessing/RecoTLR.customisePostEra_Run2_2017_pp_on_XeXe' ]
"""
_ppEra_Run2_2017_pp_on_XeXe_

Implement configuration building for data processing for proton
collision data taking for Run2

"""
3 changes: 3 additions & 0 deletions Configuration/DataProcessing/python/RecoTLR.py
Expand Up @@ -60,6 +60,9 @@ def customisePostEra_Run2_2017_harvesting_trackingOnly(process):
_customise_HPbeamspot(process)
return process

def customisePostEra_Run2_2017_pp_on_XeXe(process):
customisePostEra_Run2_2017(process)
return process

##############################################################################
def customisePPData(process):
Expand Down
8 changes: 4 additions & 4 deletions Configuration/DataProcessing/test/run_CfgTest.sh
Expand Up @@ -13,7 +13,7 @@ function runTest { echo $1 ; python $1 || die "Failure for configuration: $1" $?

runTest "${LOCAL_TEST_DIR}/RunRepack.py --select-events HLT:path1,HLT:path2 --lfn /store/whatever"

declare -a arr=("cosmics" "pp" "cosmicsEra_Run2_25ns" "cosmicsEra_Run2_2016" "ppEra_Run2_50ns" "ppEra_Run2_25ns" "ppEra_Run2_2016" "ppEra_Run2_2016_trackingLowPU" "ppEra_Run2_2016_pA" "cosmicsEra_Run2_2017" "ppEra_Run2_2017" "ppEra_Run2_2017_trackingLowPU")
declare -a arr=("cosmics" "pp" "cosmicsEra_Run2_25ns" "cosmicsEra_Run2_2016" "ppEra_Run2_50ns" "ppEra_Run2_25ns" "ppEra_Run2_2016" "ppEra_Run2_2016_trackingLowPU" "ppEra_Run2_2016_pA" "cosmicsEra_Run2_2017" "ppEra_Run2_2017" "ppEra_Run2_2017_trackingLowPU" "ppEra_Run2_2017_pp_on_XeXe")
for scenario in "${arr[@]}"
do
runTest "${LOCAL_TEST_DIR}/RunExpressProcessing.py --scenario $scenario --global-tag GLOBALTAG --lfn /store/whatever --fevt --dqmio --alcareco TkAlMinBias+SiStripCalMinBias "
Expand All @@ -22,7 +22,7 @@ do
done


declare -a arr=("cosmics" "pp" "cosmicsEra_Run2_25ns" "cosmicsEra_Run2_2016" "AlCaLumiPixels" "AlCaTestEnable" "hcalnzs" "hcalnzsEra_Run2_25ns" "hcalnzsEra_Run2_2016" "ppEra_Run2_2016_trackingLowPU" "ppEra_Run2_2016_pA" "ppEra_Run2_50ns" "ppEra_Run2_25ns" "ppEra_Run2_2016" "cosmicsEra_Run2_2017" "hcalnzsEra_Run2_2017" "ppEra_Run2_2017_trackingLowPU" "ppEra_Run2_2017")
declare -a arr=("cosmics" "pp" "cosmicsEra_Run2_25ns" "cosmicsEra_Run2_2016" "AlCaLumiPixels" "AlCaTestEnable" "hcalnzs" "hcalnzsEra_Run2_25ns" "hcalnzsEra_Run2_2016" "ppEra_Run2_2016_trackingLowPU" "ppEra_Run2_2016_pA" "ppEra_Run2_50ns" "ppEra_Run2_25ns" "ppEra_Run2_2016" "cosmicsEra_Run2_2017" "hcalnzsEra_Run2_2017" "ppEra_Run2_2017_trackingLowPU" "ppEra_Run2_2017_pp_on_XeXe" "ppEra_Run2_2017")
for scenario in "${arr[@]}"
do
runTest "${LOCAL_TEST_DIR}/RunPromptReco.py --scenario $scenario --reco --aod --dqmio --global-tag GLOBALTAG --lfn=/store/whatever --alcareco TkAlMinBias+SiStripCalMinBias"
Expand All @@ -36,14 +36,14 @@ do
done


declare -a arr=("cosmics" "pp" "cosmicsEra_Run2_25ns" "cosmicsEra_Run2_2016" "HeavyIons" "HeavyIonsEra_Run2_HI" "AlCaLumiPixels" "ppEra_Run2_50ns" "ppEra_Run2_25ns" "ppEra_Run2_2016" "ppEra_Run2_2016_trackingLowPU" "ppEra_Run2_2016_pA" "cosmicsEra_Run2_2017" "ppEra_Run2_2017_trackingLowPU" "ppEra_Run2_2017")
declare -a arr=("cosmics" "pp" "cosmicsEra_Run2_25ns" "cosmicsEra_Run2_2016" "HeavyIons" "HeavyIonsEra_Run2_HI" "AlCaLumiPixels" "ppEra_Run2_50ns" "ppEra_Run2_25ns" "ppEra_Run2_2016" "ppEra_Run2_2016_trackingLowPU" "ppEra_Run2_2016_pA" "cosmicsEra_Run2_2017" "ppEra_Run2_2017_trackingLowPU" "ppEra_Run2_2017_pp_on_XeXe" "ppEra_Run2_2017")
for scenario in "${arr[@]}"
do
runTest "${LOCAL_TEST_DIR}/RunAlcaSkimming.py --scenario $scenario --lfn=/store/whatever --global-tag GLOBALTAG --skims SiStripCalZeroBias,SiStripCalMinBias,PromptCalibProd"
runTest "${LOCAL_TEST_DIR}/RunDQMHarvesting.py --scenario $scenario --lfn /store/whatever --run 12345 --dataset /A/B/C --global-tag GLOBALTAG"
done

declare -a arr=("ppEra_Run2_50ns" "ppEra_Run2_25ns" "ppEra_Run2_2016" "ppEra_Run2_2016_trackingLowPU" "ppEra_Run2_2016_pA" "ppEra_Run2_2017" "ppEra_Run2_2017_trackingLowPU")
declare -a arr=("ppEra_Run2_50ns" "ppEra_Run2_25ns" "ppEra_Run2_2016" "ppEra_Run2_2016_trackingLowPU" "ppEra_Run2_2016_pA" "ppEra_Run2_2017" "ppEra_Run2_2017_trackingLowPU" "ppEra_Run2_2017_pp_on_XeXe")
for scenario in "${arr[@]}"
do
runTest "${LOCAL_TEST_DIR}/RunPromptReco.py --scenario $scenario --reco --aod --miniaod --dqmio --global-tag GLOBALTAG --lfn=/store/whatever --alcareco TkAlMinBias+SiStripCalMinBias"
Expand Down
6 changes: 6 additions & 0 deletions Configuration/Eras/python/Era_Run2_2017_pp_on_XeXe_cff.py
@@ -0,0 +1,6 @@
import FWCore.ParameterSet.Config as cms

from Configuration.Eras.Era_Run2_2017_cff import Run2_2017
from Configuration.Eras.Modifier_pp_on_XeXe_2017_cff import pp_on_XeXe_2017

Run2_2017_pp_on_XeXe = cms.ModifierChain(Run2_2017, pp_on_XeXe_2017)
4 changes: 4 additions & 0 deletions Configuration/Eras/python/Modifier_pp_on_XeXe_2017_cff.py
@@ -0,0 +1,4 @@
import FWCore.ParameterSet.Config as cms

pp_on_XeXe_2017 = cms.Modifier()

@@ -0,0 +1,139 @@
import FWCore.ParameterSet.Config as cms

from Configuration.Generator.Pyquen2015Settings_cff import *
from Configuration.Generator.PythiaUESettings_cfi import *

generator = cms.EDFilter("HydjetGeneratorFilter",
aBeamTarget = cms.double(129.0), ## beam/target atomic number
comEnergy = cms.double(5442.0),
qgpInitialTemperature = cms.double(1.), ## initial temperature of QGP; allowed range [0.2,2.0]GeV;
qgpProperTimeFormation = cms.double(0.1), ## proper time of QGP formation; allowed range [0.01,10.0]fm/c;
hadronFreezoutTemperature = cms.double(0.125),
doRadiativeEnLoss = cms.bool(True), ## if true, perform partonic radiative en loss
doCollisionalEnLoss = cms.bool(True),
qgpNumQuarkFlavor = cms.int32(0), ## number of active quark flavors in qgp; allowed values: 0,1,2,3
numQuarkFlavor = cms.int32(0), ## to be removed
sigmaInelNN = cms.double(70),
shadowingSwitch = cms.int32(1),
nMultiplicity = cms.int32(9000),
fracSoftMultiplicity = cms.double(1.),
maxLongitudinalRapidity = cms.double(2.3),
maxTransverseRapidity = cms.double(1.15),
angularSpectrumSelector = cms.int32(1),
rotateEventPlane = cms.bool(True),
allowEmptyEvents = cms.bool(False),
embeddingMode = cms.bool(False),
hydjetMode = cms.string('kHydroQJets'),

PythiaParameters = cms.PSet(
pythiaUESettingsBlock,
hydjetPythiaDefault = cms.vstring('MSEL=0 ! user processes',
'CKIN(3)=8.15',# ! ptMin
'MSTP(81)=1'
),
myParameters = cms.vstring('MDCY(310,1)=0'),
pythiaJets = cms.vstring('MSUB(11)=1', # q+q->q+q
'MSUB(12)=1', # q+qbar->q+qbar
'MSUB(13)=1', # q+qbar->g+g
'MSUB(28)=1', # q+g->q+g
'MSUB(53)=1', # g+g->q+qbar
'MSUB(68)=1' # g+g->g+g
),
pythiaPromptPhotons = cms.vstring('MSUB(14)=1', # q+qbar->g+gamma
'MSUB(18)=1', # q+qbar->gamma+gamma
'MSUB(29)=1', # q+g->q+gamma
'MSUB(114)=1', # g+g->gamma+gamma
'MSUB(115)=1' # g+g->g+gamma
),
pythiaZjets = cms.vstring('MSUB(15)=1',
'MSUB(30)=1'),
pythiaCharmoniumNRQCD = cms.vstring('MSUB(421) = 1',
'MSUB(422) = 1',
'MSUB(423) = 1',
'MSUB(424) = 1',
'MSUB(425) = 1',
'MSUB(426) = 1',
'MSUB(427) = 1',
'MSUB(428) = 1',
'MSUB(429) = 1',
'MSUB(430) = 1',
'MSUB(431) = 1',
'MSUB(432) = 1',
'MSUB(433) = 1',
'MSUB(434) = 1',
'MSUB(435) = 1',
'MSUB(436) = 1',
'MSUB(437) = 1',
'MSUB(438) = 1',
'MSUB(439) = 1'
),

pythiaBottomoniumNRQCD = cms.vstring('MSUB(461) = 1',
'MSUB(462) = 1',
'MSUB(463) = 1',
'MSUB(464) = 1',
'MSUB(465) = 1',
'MSUB(466) = 1',
'MSUB(467) = 1',
'MSUB(468) = 1',
'MSUB(469) = 1',
'MSUB(470) = 1',
'MSUB(471) = 1',
'MSUB(472) = 1',
'MSUB(473) = 1',
'MSUB(474) = 1',
'MSUB(475) = 1',
'MSUB(476) = 1',
'MSUB(477) = 1',
'MSUB(478) = 1',
'MSUB(479) = 1',
),
pythiaWeakBosons = cms.vstring('MSUB(1)=1',
'MSUB(2)=1'),
pythiaQuarkoniaSettings = cms.vstring('PARP(141)=1.16', # Matrix Elements
'PARP(142)=0.0119',
'PARP(143)=0.01',
'PARP(144)=0.01',
'PARP(145)=0.05',
'PARP(146)=9.28',
'PARP(147)=0.15',
'PARP(148)=0.02',
'PARP(149)=0.02',
'PARP(150)=0.085',
# Meson spin
'PARJ(13)=0.60',
'PARJ(14)=0.162',
'PARJ(15)=0.018',
'PARJ(16)=0.054',
# Polarization
'MSTP(145)=0',
'MSTP(146)=0',
'MSTP(147)=0',
'MSTP(148)=1',
'MSTP(149)=1',
# Chi_c branching ratios
'BRAT(861)=0.202',
'BRAT(862)=0.798',
'BRAT(1501)=0.013',
'BRAT(1502)=0.987',
'BRAT(1555)=0.356',
'BRAT(1556)=0.644'
),
# Quarkonia and Weak Bosons added back upon dilepton group's request.
parameterSets = cms.vstring('pythiaUESettings',
'hydjetPythiaDefault',
'myParameters',
'pythiaJets',
'pythiaPromptPhotons',
'pythiaZjets',
'pythiaBottomoniumNRQCD',
'pythiaCharmoniumNRQCD',
'pythiaQuarkoniaSettings',
'pythiaWeakBosons'
)
),
cFlag = cms.int32(1),
bMin = cms.double(0),
bMax = cms.double(30),
bFixed = cms.double(0)
)
2 changes: 2 additions & 0 deletions Configuration/PyReleaseValidation/python/relval_standard.py
Expand Up @@ -423,6 +423,8 @@
workflows[140] = ['',['HydjetQ_B12_5020GeV_2011','DIGIHI2011','RECOHI2011','HARVESTHI2011']]
### Run II cond., 2015
workflows[145] = ['',['HydjetQ_B12_5020GeV_2015','DIGIHI2015','RECOHI2015','HARVESTHI2015']]
### Run II cond., 2017
workflows[148] = ['',['HydjetQ_MinBias_XeXe_5442GeV_2017','DIGIHI2017','RECOHI2017','HARVESTHI2017']]
### Run II cond., 2018
workflows[150] = ['',['HydjetQ_B12_5020GeV_2018','DIGIHI2018','RECOHI2018','HARVESTHI2018']]
workflows[150.1] = ['',['QCD_Pt_80_120_13_HI','DIGIHI2018','RECOHI2018','HARVESTHI2018']]
Expand Down