1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
|
import FWCore.ParameterSet.Config as cms
from DQMOffline.Alignment.DiMuonMassBiasClient_cfi import DiMuonMassBiasClient as diMuonMassBiasClient
# Z-> mm
__selectionName = 'TkAlDiMuonAndVertex'
ALCARECOTkAlZMuMuMassBiasClient = diMuonMassBiasClient.clone(
FolderName = "AlCaReco/"+__selectionName,
fitBackground = True,
useRooCMSShape = True,
fit_par = dict(mean_par = [90., 89., 100.], # parameters of the mass pole
width_par = [2.49, 2.48, 2.50], # natural width of the Z
sigma_par = [1.5, 0.0, 10.0]) # detector smearing
)
alcaTkAlZMuMuBiasClients = cms.Sequence(ALCARECOTkAlZMuMuMassBiasClient)
# J/psi -> mm
__selectionName = 'TkAlJpsiMuMu'
ALCARECOTkAlJpsiMuMuMassBiasClient = diMuonMassBiasClient.clone(
FolderName = "AlCaReco/"+__selectionName,
fitBackground = True,
useRooCBShape = True, # crystal-ball fit
useRooCMSShape = True,
fit_par = dict(mean_par = [3.09, 2.7, 3.4], # parameters of the mass pole
width_par = [1.0, 0.0, 5.0], # not actually used
sigma_par = [0.01, 0.0, 5.0]) # width of the CB
)
alcaTkAlJpsiMuMuBiasClients = cms.Sequence(ALCARECOTkAlJpsiMuMuMassBiasClient)
# Upsilon -> mm
__selectionName = 'TkAlUpsilonMuMu'
ALCARECOTkAlUpsilonMuMuMassBiasClient = diMuonMassBiasClient.clone(
FolderName = "AlCaReco/"+__selectionName,
fitBackground = True,
useRooCBShape = True, # crystal-ball fit
useRooCMSShape = False, # using the exponential is useful to model the onset of Y(2S) peak
fit_par = dict(mean_par = [9.46, 8.9, 9.9], # parameters of the mass pole
width_par = [1.0, 0.0, 5.0], # not actually used
sigma_par = [1.0, 0.0, 5.0]) # width of the CB
)
alcaTkAlUpsilonMuMuBiasClients = cms.Sequence(ALCARECOTkAlUpsilonMuMuMassBiasClient)
|