The BsToJPsiPhi analyzer

NOTE of the Authors: If you don't find you answer here, look at our working group twiki page. CMS.BsJpsiPhi_AWG That's certainly up-to-date.

Motivation (Purpose of analyzer)

. Measurement of the width difference of the weak eigenstates (untagged analysis) . Bs lifetime

Code package in CVS and tags

on web:
http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/HeavyFlavorAnalysis/BsToJpsiPhi/ [called H= home smile ]
important files:
.. python file: H/test/BsToJpsiPhiAnalysis_cfg.py
.. main analyzer: H/src/BsToJpsiPhiAnalysis.cc, and
H/src/BsToJpsiPhiRootTree.cc
.. python configuration files for GEN-SIM-HLT-RECO (tag:V00-226-02):
H/test/python/BsToJpsiPhi_uptoHLT_1.py
H/test/python/PYTHIA6_Bs2JpsiPhi_10TeV_cff_py_RAW2DIGI_RECO_IDEAL_noHLTcut.py
H/test/python/PYTHIA6_Bs2JpsiPhi_10TeV_cff_py_RAW2DIGI_RECO_IDEAL_withHLTcut.py

Working release

CMSSW_2_2_X

last tags and bugs found (in inverse order)

- 30/06/2009:
A small bug was found which affects all the 2XY releases. Two packages have therefore to be checked-out from CVS into your area and compiled. The packages with the correct tags are:

* RecoVertex/KinematicFitPrimitives: V01-00-03
* RecoVertex/KinematicFit: V01-00-03
(see https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideKinematicVertexFit for more details)

- 23/06/2009:
New offline reconstruction code committed with tag V00-226-07. This new versions contains two flag to separate the Bd-->Jpsi K*(Kpi) and Bd-->Jpsi Ks(pipi) components in the MC Summer '08 Jpsi X sample and a check for the MC match with a further flag.

- 16/06/2009:
New offline reconstruction code committed with tag V00-226-06. This new versions contains two flag to evaluate mumuKpi and mumupipi non resonant components in the MC Summer '08 Jpsi X sample.

Configuration files

H/test/BsToJpsiPhiAnalysis_cfg.py

How to run the 3_1_X analyzer ----- NEW

- Create a work area and build the application (we are working nowadays in CMSSW_3_1_2)

- cvs co HeavyFlavorAnalysis/BsToJpsiPhi
- cvs co GeneratorInterface/GenFilters
(of the GeneratorInterfacest pkg these the important files:
in GeneratorInterface/GenFilters/src:
-- BdecayFilter.cc
-- SealModule.cc
in GeneratorInterface/GenFilters/interface:
-- BdecayFilter.h,

(Note
We have a private version of these files adapted for our code.
1] create the dir structure doing the cvs co of the package.
2] substitute the above mentioned files in interface and src directories, copying their modified version from ~azzolini/public/GenFilters_3_1_2 )

- run (from CMSSW_3_1_2 ) and solve the problems if there are any:
scramv1 b -j 4

- cd HeavyFlavorAnalysis/BsToJpsiPhi/test
- choose a dataset and put it in input in the python file in "cms.Source":
-- from DBS ( https://cmsweb.cern.ch/dbs_discovery/ ) using the windows or typing,e.g. " find dataset where dataset like *Jpsimu*Summer09*"
, NOTE: 2009 pre-production /BtoJPsiMuMu/Summer09-MC_31X_V3_preproduction_312_AODSIM-v2/AODSIM
-- from a more working group local repository, not yet available )
- try to run:
cmsRun BsToJpsiPhiAnalysis_cfg.py

- you should obtain 2 output root files.

How to run the 2_2_X analyzer

- Create a work area and build the application (we are working nowadays in CMSSW_2_2_6) ( how to: https://twiki.cern.ch/twiki/bin/view/CMS/WorkBookSetComputerNode#CreateWork )

- cvs co HeavyFlavorAnalysis/BsToJpsiPhi
- cvs co GeneratorInterface/GenFilters
(of the last one the important files are:
in GeneratorInterface/GenFilters/src:
-- BdecayFilter.cc
-- PythiaFilter.cc
-- SealModule.cc
in GeneratorInterface/GenFilters/interface:
-- BdecayFilter.h,
-- PythiaFilter.h)

(Note
create the dir structure doing the cvs co of the packages.
Then cancel from interface and src directory all the files and substitute them
copying their modified version from ~azzolini/public/GenFilters )

- run (from CMSSW_2_2_X ) and solve the problems if there are any:
scramv1 b

- cd HeavyFlavorAnalysis/BsToJpsiPhi/test
- choose a dataset and put it in input in the python file in "cms.Source":
-- from DBS ( https://cmsweb.cern.ch/dbs_discovery/ ) using the windows or typing,e.g. " find dataset where dataset like *Jpsimu*Summer08*"
-- from a more working group local repository, e.g., rfio:/castor/cern.ch/user/a/azzolini/BsToJpsiPhi_uptoRECO_1_CMSSW223.root )
- try to run:
cmsRun BsToJpsiPhiAnalysis_cfg.py

- you should obtain 2 output root files.

Parameter settings (see python file)

- cut on "candJpsi": [3.0 < mass < 3.22 & pt > 4.0] --> [3.0 < mass < 3.22]
- cut on "candphi ": [1.0 < mass < 1.04 & pt > 1.] --> [1.0 < mass < 1.04]
- cut on "kTracks" : [pt > 0.8 & abs(eta) < 2.5] --> [pt > 0. & abs(eta) < 2.5] because obsolete
- cut on "candBs" : [pt > 5.0] --> [pt > 0. && 5.3 < mass < 5.5]

Output (format)

2 output root files

How to run the analyzer on the grid

add crab file

Repository of the root files:

- azzolini: /castor/cern.ch/user/a/azzolini/BsToJpsiPhi
- cerizza: /castor/cern.ch/user/c/cerizza/BsToJpsiPhi

scripts

The background sample "InclusiveMuX" is done, if you DON'T WANT to run accessing CASTOR, download directly from Vir's account with the following script. This script allows to download files from CASTOR just changing the source directory.

http://hep.phys.utk.edu/~cerizza/copyFromCastor.py

How to run the analyzer with real data

Quoting a notable answer, listened during a meeting : "... (silence)... mmm but we don't have data".
We will let you know asap.

J/psi phi analysis with real data

Documentation:
CADI BPH-08-005
Authors: V. Azzolini (University of Helsinki)
G. Cerizza (University of Tennessee)
B. Millan Mejias (University of Zurich)
L. Wilke (University of Zurich)

J/psi phi MC analysis published in 2006

Documentation:
CADI BPH-08-005  
Study of the decay Bs-> J/psi phi -> mu+mu- K+K- CMS NOTE-2006/121 2006/11/22
Study of the decay Bs-> J/psi phi -> mu+mu- K+K- CMS AN-2006/090 2006/04/10
On Delta M (Bs) Measurement by Bs => J/Psi Phi Decay Channel CMS NOTE-2000/021 (broken link) ?

Authors: K. Prokofiev, T. Speer (University of Zurich)
V. Ciulli, N. Magini (Firenze)
S.Shulga, T. Ilicheva (JINR)
P. Galumian (Strasbourg)

Talks

. B PAG 19/05/09
http://indico.cern.ch/conferenceDisplay.py?confId=59179

Competition and other links

Contact

Any feedback to the authors is welcomed!

-- VirginiaAzzolini - 12-Mar-2010

Edit | Attach | Watch | Print version | History: r7 < r6 < r5 < r4 < r3 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r7 - 2010-03-12 - VirginiaAzzolini
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    CMSPublic All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback