optimize parameters for the upsilon->mm fit

cms-sw · Nov 7, 2022 · 993d99c · 993d99c
1 parent 262dc89
commit 993d99c
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 8 deletions.
diff --git a/DQMOffline/Alignment/python/DiMuonMassBiasHarvesting_cff.py b/DQMOffline/Alignment/python/DiMuonMassBiasHarvesting_cff.py
@@ -23,8 +23,8 @@
     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],
-                   sigma_par = [0.01, 0.0, 5.0])
+                   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)
@@ -34,11 +34,11 @@
 ALCARECOTkAlUpsilonMuMuMassBiasClient = diMuonMassBiasClient.clone(
     FolderName = "AlCaReco/"+__selectionName,
     fitBackground = True,
-    useRooCBShape = True, # crystal-ball fit
-    useRooCMSShape = 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],
-                   sigma_par = [1.0, 0.0, 5.0])
+                   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)
diff --git a/DQMOffline/Alignment/test/DiMuonTkAlDQMHarvester_cfg.py b/DQMOffline/Alignment/test/DiMuonTkAlDQMHarvester_cfg.py
@@ -114,8 +114,8 @@
     _useRooCBShape = cms.bool(True)
     _useRooCMSShape = cms.bool(True)
     _fitPar =  cms.PSet(mean_par = cms.vdouble(3.09, 2.7, 3.4),
-                        width_par = cms.vdouble(0.001, 0.0, 5e-3),
-                        sigma_par = cms.vdouble(0.1, 0.0, 5.0))
+                        width_par = cms.vdouble(1.0, 0.0, 5.0),
+                        sigma_par = cms.vdouble(0.01, 0.0, 5.0))
 elif (options.resonance == 'Upsilon'):
     print('',30*"#",'\n # will harvest Upsilon file \n',30*"#")
     _folderName =  cms.string('AlCaReco/TkAlUpsilonMuMu')