Refitting tracks with TwoBodyDecay constraints

Introduction

Description of the TwoBodyDecay algorithm

Related documents

Author of TwoBodyDecay package : Edmund Widl, R. Fruhwirth. Reference document : CMS NOTE-2007/032

Implementation and validation of the refitting using TwoBodyDecay constraint (ongoing) : Jeremy Andrea, Eric Conte, Christophe Goetzmann. Main talks :

• DPG Tracker - Alignment : Thursday, the 8th December 2011
• DPG Tracker - Alignment : Thursday, the 1st December 2011
• DPG Tracker - Alignment : Thursday, the 14th July 2011
• Alignement workshop @ DESY : Monday, the 30th May 2011
• DPG Tracker - Alignment : Thursday, the 5th May 2011

Installation of the package

The current package has been developed and validated on CMSSW release 4.4.X.

scramv1 project CMSSW CMSSW_4_4_2
cd CMSSW_4_4_2/src
cmsenv
cvs co UserCode/IPHCalignment/installPackage.sh
source UserCode/IPHCalignment/installPackage.sh

Directories called 'Alignement/TwoBodyDecay' and 'RecoTracker/TrackProducer' will be created in 'CMSSW_4_4_2/src'. To compile, use the classical CMSSW command 'scramv1 b -jX' with X the number of cores.

Details about the installator script 'installPackage.sh' can be found here : http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/IPHCalignment/installPackage.sh?view=markup

The available recipes

There is a python 'cff' file dedicated to each kind of constraint in the directory 'RecoTrack/TrackProducer/python'.

• 'TwoBodyDecayVertexConstraintProducer_cff.py' for vertex constraint (3 parameters)
• 'TwoBodyDecayMomentumConstraintProducer_cff.py' for momentum magnitude constraint (1 parameter)
• 'TwoBodyDecayConstraintProducer_cff.py' for trajectory state constraint (5 parameters)

The module sequences of these files are based on the same sketch :

• selecting the tracks with a TrackSelector
• refitting the tracks according to the alignment and calibration configuration
• constructing a constraint with TwoBodyDecay from the selected tracks. The constraint is calculated from the estimation of the parameters of a model describing the Z to mu mu kinematics. A threshold on the quality of the fit is applied.
• refitting the selected track with the constraint produced by TwoBodyDecay If TwoBodyDecay does not manage to produce the constraint (or if the constraint is rejected by the chi2 cut), the refit is done without constraint.

In order to use these files, you have just to include in your main python script the following line (in this example case : the momentum constraint) :

process.load("refitter.TwoBodyDecayMomentumConstraintProducer_cff.py")

A complete example of main python script can be found here : http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/IPHCalignment/RecoTracker/TrackProducer/test/RefitWithTBDConstraint_cfg.py?view=markup

Nota Bene :

• The collection of tracks to give in input of the sequence is by default 'generalTracks'. This statement can be changed by adding the simple line in your main script :

AlignmentTrackSelector.src = cms.InputTag( "generalTracks" )
 

• The sequence will return a collection of 'reco::Track' with the label 'TrackRefitterTBDConstraint'.

Advanced use : change recipe settings

Track selection

TwoBodyDecay

Several parts in the configuration.

• input for TwoBodyDecay fit

Parameters	Description
'beamSpot'	InputTag of the beamSpot
'secondaryMass'	mean value of the muon mass
'primaryMass'	mean value of the Z mass
'primaryWidth'	width of the Z boson = weight to the Z mass in the chi2 estimation

• internal parameter to the TwoBodyDecay fit

• cut

• error propagation to the constraint

Note : In the case of vertex position constraint, the error on X, Y and Z (coordinates of the primary vertex) are propagated as the covariance matrix of those parameters, calculated by TwoBodyDecay.

-- EricConte - 14-Feb-2012 -- ChristopheGoetzmann - 15-Feb-2012