optimize parameters for the upsilon->mm fit

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

Hide file tree

Showing 3 changed files with 9 additions and 9 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')

diff --git a/DQMOffline/Alignment/test/DiMuonTkAlDQMValidator_cfg.py b/DQMOffline/Alignment/test/DiMuonTkAlDQMValidator_cfg.py
@@ -58,7 +58,7 @@
 elif (options.resonance == 'Upsilon'):
     # upsilon
     print('',30*"#",'\n # will process Upsilon->mm data\n',30*"#")
-    readFiles = ['/store/relval/CMSSW_12_5_0_pre2/RelValUpsilon1SToMuMu_14/ALCARECO/TkAlUpsilonMuMu-124X_mcRun3_2022_realistic_v3-v1/2580000/fca73916-5076-4c9f-a460-2481588825ab.root']
+    readFiles = ['/store/relval/CMSSW_12_5_0_pre2/RelValUpsilon1SToMuMu_14/ALCARECO/TkAlUpsilonMuMu-124X_mcRun3_2022_realistic_v3-v1/2580000/fca73916-5076-4c9f-a460-2481588825ab.root']    
     process.seqALCARECOTkAlDQM = cms.Sequence(process.ALCARECOTkAlUpsilonMuMuVtxDQM + process.ALCARECOTkAlUpsilonMassBiasDQM)
 else:
     print('unrecongnized %s resonance',options.resonance)