Difference: TrackCalib (1 vs. 37)

Revision 372019-09-03 - FlavioArchilli

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 17 to 17
  %CODE{ lang="bash" num="off" }%
Changed:
<
<
LbLogin -c x86_64-slc6-gcc7-opt
>
>
# (deprecated) LbLogin -c x86_64-slc6-gcc7-opt LbLogin -c x86_64-centos7-gcc8-opt
 lb-dev --nightly lhcb-head Urania/HEAD cd UraniaDev _HEAD git lb-use Urania

Revision 362019-08-14 - JordanDanielRoth

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 22 to 22
 cd UraniaDev _HEAD git lb-use Urania git lb-checkout Urania/master TrackCalib
Added:
>
>
git lb-checkout Urania/master cmake
  make configure make

Revision 352019-06-26 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 270 to 270
 #Hadron track systematics

Tracking systematics for hadrons

Changed:
<
<
TrackCalib is meant for muon tracking efficiency evaluation and doesn't provide any information for hadrons yet! In order to evaluate the systematics related to an hadronic track the fraction of kaons and pions undergo a hadronic interaction before the last T station, around z = 9000mm, has to be evaluated.
>
>
TrackCalib is meant for muon tracking efficiency evaluation and doesn't provide any information for hadrons.

The reason for the difference in efficiency for muons and hadrons are hadronic interactions. Their cross-section depends on the particle type.

A simulation of $B^0 \rightarrow J/&amp;psi; K^{*0}$, where $K^{*0} \rightarrow K^+\pi^-$ shows that about 11% of kaons and 14% of pions cannot be reconstructed due to the hadronic interactions.

The number of tracks that undergo a hadronic interaction before the last T-station depends mostly on its momentum and the amount of material it crosses. Since there is 10% material-budget uncertainty, hadronic efficiencies have an intristic uncertanity of 1.1% for kaons and 1.4% for pions.

In order to evaluate the systematics related to a hadronic track for a specific analysis, the fraction of hadrons that undergo a hadronic interaction before the last T station (around z = 9000mm) has to be evaluated.
Latex rendering error!! dvi file was not created.

  The MCTupleToolKinematic almost does that, but it has to be changed in order to make it to give gives you only hadronic interactions, and not also decays: https://gitlab.cern.ch/lhcb/Analysis/blob/master/Phys/DecayTreeTupleMC/src/MCTupleToolKinematic.cpp#L113

Revision 342019-06-26 - FlavioArchilli

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 267 to 267
 python TrackCalib/scripts/CreateTrackEffTables.py "2015_25ns" --Plot -secondSample "MC,2016_25ns,Sim09b" %ENDCODE%
Added:
>
>
#Hadron track systematics

Tracking systematics for hadrons

TrackCalib is meant for muon tracking efficiency evaluation and doesn't provide any information for hadrons yet! In order to evaluate the systematics related to an hadronic track the fraction of kaons and pions undergo a hadronic interaction before the last T station, around z = 9000mm, has to be evaluated.

The MCTupleToolKinematic almost does that, but it has to be changed in order to make it to give gives you only hadronic interactions, and not also decays: https://gitlab.cern.ch/lhcb/Analysis/blob/master/Phys/DecayTreeTupleMC/src/MCTupleToolKinematic.cpp#L113 All kinematic cuts that are used in the specific analysis has also to be applied, as this can change the final result.

This will then give a fraction, say 14%, and multiplied by the material-budget uncertainty, which is 10%, you'll get the uncertainty per particle type.

 

Meeting contributions

Revision 332019-06-26 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 6 to 6
  TrackCalib was created as tool for creating tracking efficiency correction tables allowing for user-defined track-quality cuts, binning and variable combinations. It also allows to check 1D- and 2D-efficiency dependencies for many variables. The name and idea was shamelessly stolen from PIDCalib.
Changed:
<
<
In case of any questions of comments regarding the TrackCalibTool contact author renata.kopecna@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017 and Tutorial talk at LHCb week, 05.12.2017
>
>
In case of any questions of comments regarding the TrackCalibTool contact author renata.kopecna@cernNOSPAMPLEASE.ch and flavio.archilli@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017 and Tutorial talk at LHCb week, 05.12.2017
 
Line: 30 to 30
  More detailed instructions can be found in the Git4LHCb twiki.
Changed:
<
<
TrackCalib is currently being cleaned up, in case of problems accessing 2016 sample, contact renata.kopecna@cernNOSPAMPLEASE.ch
>
>
TrackCalib is currently being cleaned up, in case of problems accessing 2016 sample, contact renata.kopecna@cernNOSPAMPLEASE.ch

Easy quick fix in the case of a failure is to replace in TrackCalib/python/TrackCalib/Prepare.py, lines 131-143, Data_2016_Mag by Data_2016_25ns_Mag

 

TrackCalib running instructions

Revision 322019-05-07 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 418 to 418
  General naming convension is to have two lines per method: one for μ+ as a probe and one for μ- being a probe. The Jpsi/Upsilon/Z in the line stands for the mother particle of the muons.
Changed:
<
<
Strucutre of the lines is as follows: creation of probe (first, an algorithm creates muon protoParticles, from which a muon Particle is made), tag filter (applies tag cuts on HLT muons), chargeFilter that depending on the line (1 or 2) selects μ+ as a probe and μ- as a tag (or the other way around) and last an algorithm that creates the resonance from tag and probe and applies cats on the resonance.
>
>
Strucutre of the lines is as follows: first, probe μ is created (algorithm starts by creating muon protoParticles, from which a muon Particle is made), then the tag filter selectes tag μ (by applying tag cuts on HLT muons). Then chargeFilter that depending on the line (1 or 2) selects μ+ as a probe and μ- as a tag (or the other way around) is applied. Last, an algorithm that creates the resonance from tag and probe and applies cuts on the resonance is used.
 
Changed:
<
<
The tag, proge and resonance cuts are defined outside of these algorithms to increase transparency of the StrippingLine code. The units of the cuts are always MeV, in the case of p and pT the unit is GeV. All cuts in the stripping lines now agree with Run II Turbo cuts or they are looser.
>
>
The tag, probe and resonance cuts are defined outside of these algorithms to increase transparency of the StrippingLine code. The units of the cuts are always MeV, in the case of p and pT the unit is GeV. All cuts in the stripping lines now agree with Run II Turbo cuts or they are looser.
 
Changed:
<
<
Having StrippingLines in unified structure also allows for using TrackCalib very easily on Run I data.
>
>
Having StrippingLines in an unified structure also allows for using TrackCalib very easily on Run I data.
 

DownMuon (velo method)

StrippingTrackEffDownMuon
Changed:
<
<
The main change is replacing StrippingTrackEffDownMuonNominalLine by DownMuonLine1 and DownMuonLine2. StrippingTrackEffDownMuonNominalLine used two downstream tracks for tag and probe, while DownMuonLine1,2 use downstream track as a probe and long track as a tag. Study was made to check the possible bias on tracks that originate very close to the VELO edge and the effect is of second order and not significant for the efficiencies. The effect on the data/MC ratio is negligible.
>
>
The main change is replacing StrippingTrackEffDownMuonNominalLine by DownMuonLine1 and DownMuonLine2. StrippingTrackEffDownMuonNominalLine used two downstream tracks for tag and probe, while DownMuonLine1,2 use downstream track as a probe and long track as a tag. Study was made to check the possible bias on tracks that might occur to resonances creaty very close to the VELO edge. The effect is of second order and it is not significant for the efficiencies. The effect on the data/MC ratio is negligible.
 

VeloMuon (T-station method)

StrippingTrackEffVeloMuon

Changed:
<
<
Variables used for cuts in this stripping line were redefined in a way that there is a clear distinction between prob, tag and resonance (Jpsi/Z/Upsilon) cuts.
>
>
Variables used for cuts in this stripping lines were redefined in a way that there is a clear distinction between probe, tag and resonance (Jpsi/Z/Upsilon) cuts.
 

MuonTT (Long method)

StrippingTrackEffMuonTT

Revision 312019-05-06 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 409 to 409
 
    • Tuple tool designed for creating TrackCalib input from Turbo. Replaced by the TupleTool in Castelao!
Added:
>
>

Stripping lines

In the restripping campaign in 2019 (restripping of both Run I and Run II data), significant changes in the stripping lines were made. Main changes are summarized in this talk, however, the list in the talk is not exhaustive.

General naming convension is to have two lines per method: one for μ+ as a probe and one for μ- being a probe. The Jpsi/Upsilon/Z in the line stands for the mother particle of the muons.

Strucutre of the lines is as follows: creation of probe (first, an algorithm creates muon protoParticles, from which a muon Particle is made), tag filter (applies tag cuts on HLT muons), chargeFilter that depending on the line (1 or 2) selects μ+ as a probe and μ- as a tag (or the other way around) and last an algorithm that creates the resonance from tag and probe and applies cats on the resonance.

The tag, proge and resonance cuts are defined outside of these algorithms to increase transparency of the StrippingLine code. The units of the cuts are always MeV, in the case of p and pT the unit is GeV. All cuts in the stripping lines now agree with Run II Turbo cuts or they are looser.

Having StrippingLines in unified structure also allows for using TrackCalib very easily on Run I data.

DownMuon (velo method)

StrippingTrackEffDownMuon

The main change is replacing StrippingTrackEffDownMuonNominalLine by DownMuonLine1 and DownMuonLine2. StrippingTrackEffDownMuonNominalLine used two downstream tracks for tag and probe, while DownMuonLine1,2 use downstream track as a probe and long track as a tag. Study was made to check the possible bias on tracks that originate very close to the VELO edge and the effect is of second order and not significant for the efficiencies. The effect on the data/MC ratio is negligible.

VeloMuon (T-station method)

StrippingTrackEffVeloMuon

Variables used for cuts in this stripping line were redefined in a way that there is a clear distinction between prob, tag and resonance (Jpsi/Z/Upsilon) cuts.

MuonTT (Long method)

StrippingTrackEffMuonTT

All the previously hard-coded cuts were now moved to the config.

LowMultJpsiLines were added. These lines are designed to select Jpsis from events with very low occupancy, making them very useful for very-low multiplicity production studies.

Revision 302019-04-02 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 228 to 228
  Release date: TBA
Changed:
<
<
Will include 2017 DATA and MC. Will be released once the WG production of 2017 data and 2017 Sim09D MC is finished.
>
>
Will include 2017 DATA and MC. Will be released once the WG production of 2017 and 2018 data and 2017 Sim09D, 2018 SIM09E MC is finished.

Release date: TBA

Fixes:

2018 and 2017 available!

Added division-by-zero check into Fit.py

Organized the WG directory in a more natural way.

Prepared for including Run I efficiencies.

New empty default of -ignoreBins: default: "".

New options:

-plotMassFits Plot also mass fits in all bins for given sample and variables. default: False

-plotOnlyOne Plots only firstSample (see below) efficiency. default: False

New features:

Disabled the annoying pop-up canvases in Fit.py

-firstSample -secondSample Allows for plotting a comparison of two arbitrary samples, including the ratio of these efficiencies. Has to be in the format of mode, year, simVer (if mode is MC). default: -firstSample: "data,year," -secondSample: "MC,year,simVer".

As and example, let's say we want to compare 2015 Data (25ns) (which is our default, therefore -firstSample doesn't have to be specified) and 2016 Sim09b sample:

<!-- SyntaxHighlightingPlugin -->
python TrackCalib/scripts/CreateTrackEffTables.py "2015_25ns" --Plot  -secondSample "MC,2016_25ns,Sim09b"
<!-- end SyntaxHighlightingPlugin -->
 

Meeting contributions

Revision 292019-02-20 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 30 to 30
  More detailed instructions can be found in the Git4LHCb twiki.
Added:
>
>
TrackCalib is currently being cleaned up, in case of problems accessing 2016 sample, contact renata.kopecna@cernNOSPAMPLEASE.ch
 

TrackCalib running instructions

Line: 242 to 244
 
Changed:
<
<
>
>

 

Bugs

Revision 282018-11-23 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 242 to 242
 
Added:
>
>
 

Bugs

Revision 272018-10-18 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 278 to 278
 
Changed:
<
<
      • MC is MC? Default: false=
      • =Turbo use Turbo or stripping output? Default: true
>
>
      • MC is MC? Default: false
      • Turbo use Turbo or stripping output? Default: true
 
    • Turbo-related options:
      • MatchedLocation1 location of muons in turbo data. Default LowStat trigger line
      • MatchedLocation2 location of muons in turbo data. Default HighStat trigger line
    • Stripping-related options:
Changed:
<
<
      • Method Which method is used? Options are "MuonTT" (long), "VeloMuon" (T), "DownMuon" (velo) or "Downstream" (velo). Default: MuonTT=
      • VeloAssocFrac= Minimal fraction of common VELO hits when a track is considered associate. Default: =0.5
      • MuonAssocFrac= Minimal fraction of common MUON hits when a track is considered associate. Default: =0.4
      • TTAssocFrac= Minimal fraction of common TT hits when a track is considered associate. Default: =0.4
      • ITAssocFrac= Minimal fraction of common IT hits when a track is considered associate. Default: =0.4
      • =OTAssocFrac Minimal fraction of common OT hits when a track is considered associate. Default: 0.4
>
>
      • Method Which method is used? Options are "MuonTT" (long), "VeloMuon" (T), "DownMuon" (velo) or "Downstream" (velo). Default: MuonTT
      • VeloAssocFrac Minimal fraction of common VELO hits when a track is considered associate. Default: 0.5=
      • MuonAssocFrac Minimal fraction of common MUON hits when a track is considered associate. Default: 0.4=
      • TTAssocFrac Minimal fraction of common TT hits when a track is considered associate. Default: 0.4=
      • ITAssocFrac Minimal fraction of common IT hits when a track is considered associate. Default: 0.4=
      • OTAssocFrac Minimal fraction of common OT hits when a track is considered associate. Default: 0.4
    • Output:
      • -head_Assoc: Does the probe track have an associated track?
      • -head_hasTTHits: Does the probe track have an associated track?
      • -head_Overlap_TTFraction: Maximum overlap fraction of the probe track in TT
      • -head_Overlap_TFraction: Maximum overlap fraction of the probe track in T (meaning OT+IT)
      • -head_Overlap_ITFraction: Maximum overlap fraction of the probe track in IT
      • -head_Overlap_OTFraction: Maximum overlap fraction of the probe track in OT
      • -head_Overlap_VeloFraction: Maximum overlap fraction of the probe track in VELO
      • -head_Overlap_MuonFraction: Maximum overlap fraction of the probe track in MUONs
      • -head_maxOverlap: Maximum sum of all previous overlap fractions.
      • -head_Matched_StateCTB_tx: state closest to beam Tx
      • -head_Matched_StateCTB_ty: state closest to beam Ty
      • -head_Matched_StateCTB_x: state closest to beam X
      • -head_Matched_StateCTB_y: state closest to beam Y
      • -head_Matched_StateCTB_z: state closest to beam Z
      • -head_Matched_Px: x component of momentum of track matched to the probe track.
      • -head_Matched_Py: y component of momentum of track matched to the probe track.
      • -head_Matched_Pz: z component of momentum of track matched to the probe track.
      • -head_Matched_dP: Momentum error of associated track
      • -head_Matched_probChi2: ProbChi2 of track matched to the probe track.
      • -head_Matched_Chi2NDoF: Chi2 / ndof of track matched to the probe track.
      • -head_Matched_VeloChi2: ProbChi2 of track matched to the probe track.
      • -head_Matched_VeloNDoF: Chi2 / ndof of velo track matched to the probe track.
      • -head_Matched_GhostProb: Chi2 / ndof of velo track matched to the probe track.
      • -head_DeltaR: Delta R between the probe track and track matched to the probe track.
      • -head_DeltaP: Momentum difference between the probe track and track matched to the probe track.
      • -head_DeltaPhi: Phi difference between the probe track and track matched to the probe track.
      • -head_DeltaEta: Pseudorapidity difference between the probe track and track matched to the probe track.
      • -head_Matched_SameMCPart: probe MC particle has an associated MC particle
      • -head_MCPart: Does the probe track have an MC particle?
      • -head_Matched_MCPart: Does the associated track have an MC particle?
      • -head_MatchedMCPx: Associated MC particle x component of momentum
      • -head_MatchedMC_Py:Associated MC particle y component of momentum
      • -head_MatchedMC_Pz: Associated MC particle z component of momentum
      • -head_MatchedMC_P: Associated MC particle momentum
      • -head_MatchedMC_PID: Associated MC particle PID
 
Changed:
<
<
* ToBeFilled
>
>
    • TupleTool for creating long method (MuonTT tracks) tuples suitable for efficiency studies from Stripping.
    • Parameters of the tool:
      • TTAssocFraction: Minimal fraction of common TT hits when a track is considered associate
      • MuonAssocFraction: Minimal fraction of common muon hits when a track is considered associate
      • InverseMatch: Start tag-and-probe from J/psi
      • InverseMatchLocation: Location of the J/psi from stripping
      • MotherLocation: Location of the J/psi from stripping (for purityAssoc)
      • MC: Use MC output?
    • Basic output (the MuonTT suffix is there in order to distinguish output from this tool from the output of TupleToolTrackEffMatch)
      • -head_Overlap_MuonFraction_MuonTT: The maximum fraction of muon hits matched to the long track
      • -head_TFraction_MuonTT: The maximum fraction of muon hits matched to the long track
      • -head_hasTTHits_MuonTT: Has the probe track TT hits?
      • -head_TTHitAssoc_MuonT: Is MuonTT track matched to a long track with TT hits?
      • -head_noTTHitAssoc_MuonTT: Is MuonTT track matched to a long track without TT hits?
    • if Verbose is true:
      • -head_upgradedMass_MuonTT: Invariant mass of tag track and track associated to probe track
      • -head_nTTHits_MuonTT: Number of TT hits of MuonTT track.
      • -head_Matched_nTTHits_MuonTT: Number of TT hits on track matched to MuonTT track.
      • -head_Matched_Px_MuonTT: x component of momentum of track matched to MuonTT track.
      • -head_Matched_Py_MuonTT: y component of momentum of track matched to MuonTT track.
      • -head_Matched_Pz_MuonTT: z component of momentum of track matched to MuonTT track.
      • -head_Matched_P_MuonTT: absolute value of momentum of track matched to MuonTT track.
      • -head_Matched_dP_MuonTT: error on first track state of track matched to MuonTT track.
      • -head_Matched_probChi2_MuonTT: ProbChi2 of track matched to MuonTT track.
      • -head_Matched_Chi2NDoF_MuonTT: Chi2 / ndof of track matched to MuonTT track.
      • -head_DeltaR_MuonTT: Delta R between MuonTT track and track matched to MuonTT track.
      • -head_DeltaP_MuonTT: Momentum difference between MuonTT track and track matched to MuonTT track.
      • -head_wrongDeltaR_MuonTT: Delta R between MuonTT track and track closest to MuonTT track, not being the matched one. A track was matched in this event.
      • -head_wrongDeltaP_MuonTT: Momentum difference between MuonTT track and track closest to MuonTT track, not being the matched one. A track was matched in this event.
      • -head_wrongP_MuonTT: Momentum of track closest to MuonTT track, not being the matched one. A track was matched in this event.
      • -head_wrongDeltaR2_MuonTT: Delta R between MuonTT track and track closest to MuonTT track. No track was matched in this event.
      • -head_wrongDeltaP2_MuonTT: Momentum difference between MuonTT track and track closest to MuonTT track. No track was matched in this event.
      • -head_wrongP2_MuonTT: Momentum of track closest to MuonTT track. No track was matched in this event.
      • -head_Matched_MCPart_T: Has the probe track matched MC particle?
      • -head_MatchedMC_Px_MuonTT: x component of momentum of MC particle mached to the track matched to MuonTT track.
      • -head_MatchedMC_Py_MuonTT: y component of momentum of MC particle mached to the track matched to MuonTT track.
      • -head_MatchedMC_Pz_MuonTT: z component of momentum of MC particle mached to the track matched to MuonTT track.
      • -head_MatchedMC_P_MuonTT: Absolute value of momentum of MC particle mached to the track matched to MuonTT track.
      • -head_MatchedMC_PID_MuonTT: PID of MC particle mached to the track matched to MuonTT track.
 

Revision 262018-10-18 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 278 to 277
 

Tracking efficiency TupleTools

Changed:
<
<
* Modified version ToBeFilled
>
>
    • TupleTool for creating TrackCalib input from both Turbo and Stripping. Input parameters are:
      • MC is MC? Default: false=
      • =Turbo use Turbo or stripping output? Default: true
    • Turbo-related options:
      • MatchedLocation1 location of muons in turbo data. Default LowStat trigger line
      • MatchedLocation2 location of muons in turbo data. Default HighStat trigger line
    • Stripping-related options:
      • Method Which method is used? Options are "MuonTT" (long), "VeloMuon" (T), "DownMuon" (velo) or "Downstream" (velo). Default: MuonTT=
      • VeloAssocFrac= Minimal fraction of common VELO hits when a track is considered associate. Default: =0.5
      • MuonAssocFrac= Minimal fraction of common MUON hits when a track is considered associate. Default: =0.4
      • TTAssocFrac= Minimal fraction of common TT hits when a track is considered associate. Default: =0.4
      • ITAssocFrac= Minimal fraction of common IT hits when a track is considered associate. Default: =0.4
      • =OTAssocFrac Minimal fraction of common OT hits when a track is considered associate. Default: 0.4
 
Deleted:
<
<
  * ToBeFilled
Deleted:
<
<
 
Changed:
<
<
* Tuple tool designed for Turbo output. Replaced by the TupleTool in Castelao!
>
>
    • Tuple tool designed for creating TrackCalib input from Turbo. Replaced by the TupleTool in Castelao!

Revision 252018-10-18 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 228 to 228
  Will include 2017 DATA and MC. Will be released once the WG production of 2017 data and 2017 Sim09D MC is finished.
Changed:
<
<
#Talks
>
>
 

Meeting contributions

Added:
>
>
  * Tracking and Alignment meeting, 21.12.2017 Bug in the tracking efficiency HLT lines report
Added:
>
>
 
Changed:
<
<
#Bugs
>
>
 

Bugs

2016 TurCal bug

Line: 259 to 269
 

TupleToolTrackEffMatch Bug

Does not affect the efficiency results. This bug effects 2015 and 2016 results. Due to a wrong if, Overlap variables in the tuples can only have values 0 or 1. This means custom cuts on the Overlap are not possible. The cuts used directly in the trigger are OverlapTT, OverlapIT, OverlapMuon > 0.4, OverlapVelo > 0.5.

Added:
>
>

For experts

WG production

Tracking efficiency TupleTools

* Modified version ToBeFilled

* ToBeFilled

* Tuple tool designed for Turbo output. Replaced by the TupleTool in Castelao!

Revision 242018-10-18 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 228 to 228
  Will include 2017 DATA and MC. Will be released once the WG production of 2017 data and 2017 Sim09D MC is finished.
Added:
>
>
#Talks

Meeting contributions

* Tracking and Alignment meeting, 21.12.2017 Bug in the tracking efficiency HLT lines report

 #Bugs

Bugs

Revision 232018-09-24 - CarlosVazquezSierra

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 17 to 17
  %CODE{ lang="bash" num="off" }%
Changed:
<
<
lb-dev Urania/v7r0 cd UraniaDev_v7r0
>
>
LbLogin -c x86_64-slc6-gcc7-opt lb-dev --nightly lhcb-head Urania/HEAD cd UraniaDev_HEAD
 git lb-use Urania git lb-checkout Urania/master TrackCalib

Revision 222018-09-24 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 109 to 109
  -matchCrit "Probe:<criterion1,criterion2>;Global:<criterion3>": Define tighter matching criterion (when is an event considered efficient), e.g. Matched_GhostProb < 0.3, separation as above
Changed:
<
<
-maxEntries 1234567: Sets the maximum number of events used per method and charge. Needed to limit memory consumption of fit script; default: 4.e7
>
>
-maxEntries 1234567: Sets the maximum number of events used per method and charge. Needed to limit memory consumption of fit script; default: 2.5e8
  As an example, running Prepare script over MC only, Sim08h with 100000 entries using only T-station method without multiplicity reweighting requires the following command:
Line: 196 to 196
  Release date: 14.09.2018
Changed:
<
<
Features fixes:
>
>
Fixes:
 
Changed:
<
<
-MaxEntries the new default is -1. This is especially handy when using special cuts. Plotting either Velo of T-station now doesn't crash (found a way around an if needed by TGraphAsymErrors
>
>
-MaxEntries the new default is 2.5e8. This allows to access the full data sample in 2016. This will be especially handy when using special cuts.
 
Changed:
<
<
Fix of the naming convention when plotting the efficiencies in order to produce ready-to-use .C plots.
>
>
Plotting either Velo of T-station now doesn't crash
 
Changed:
<
<
Fix using two methods at the same time.
>
>
Fixed of the naming convention when plotting the efficiencies in order to produce ready-to-use .C plots.

Fixed using two methods at the same time.

 
New options:
Line: 218 to 219
 
New features:

When fitting, a two log files are produced. First, called trackEff_mode_method_status.out returns the ROOT fit status (should be 0), χ^2 /NDOF of the fit, all parameter values with errors, its the relative errors and final efficiency values.

Changed:
<
<
Second, file called WARNINGS_trackEff_mode_method.out is created. If the fit fails, the χ^2 /NDOF>chi2Limit or when the relative error is bigger than relErrLimit, an entry is created in the WARNINGS file. Always check at least the warnings file for problematic fits.
>
>
Second, file called WARNINGS_trackEff_mode_method.out is created. If the fit fails, the χ^2 /NDOF>chi2Limit or when the relative error is bigger than relErrLimit, an entry is created in the WARNINGS file. The limits at the moment are chi2Limit = 5.0 and relErrLimit = 0.5. The values can be changed here Always check at least the warnings file for problematic fits.
 

Version 2.0

Line: 249 to 250
  Do not cut on/bin in anything that is produced by TupleToolTrackInfo (such as TRACK CHI2NDOF). Overlap information or matching information is not safe neither.
Changed:
<
<

TupleToolTrackEffMatch

>
>

TupleToolTrackEffMatch Bug

  Does not affect the efficiency results. This bug effects 2015 and 2016 results. Due to a wrong if, Overlap variables in the tuples can only have values 0 or 1. This means custom cuts on the Overlap are not possible. The cuts used directly in the trigger are OverlapTT, OverlapIT, OverlapMuon > 0.4, OverlapVelo > 0.5.
Deleted:
<
<
Bin borders separated by comma, variable+binning sets by semicolon. Format follows the usual ROOT definition (Int_t nbinsx, Double_t xlow, Double_t xup)

Revision 212018-09-11 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 88 to 88
  -polarity "MagUp/MagDown": specify to run only over MagUp or MagDown. Does not affect the Prepare script
Changed:
<
<
Hence, for example to run over 2016 collision data only, use Long and Velo methods for PT and ETA with GhostProb < 0.5, run:
>
>
Hence, for example to run over 2016 collision data only, use Long and Velo methods for PT and ETA with GhostProb < 0.5 of the probe, run:
  %CODE{ lang="bash" num="off" }%
Changed:
<
<
python TrackCalib/scripts/CreateTrackEffTables.py "2016_25ns" -mode "Data" -method "Long, Velo" -variables "PT,ETA" -binning "PT:0,1000,5000" -variables2D "PT-ETA" -matchCrit "Probe: GhostProb <0.5;Global:GhostProb<0.5"
>
>
python TrackCalib/scripts/CreateTrackEffTables.py "2016_25ns" -mode "Data" -method "Long,Velo" -variables "PT,ETA" -binning "PT:0,1000,5000" -variables2D "PT-ETA" -matchCrit "Probe: GhostProb<0.5"
 %ENDCODE%

Note that the -matchCrit is listed in the Prepare options. This example might be useful to check the effect of looser cut on GhostProb on PT and ETA on data. Note that since PT is not a default variable (P and ETA are), binning in PT needs to be specified. For more information about the binning option see Fit options.

Line: 111 to 111
  -maxEntries 1234567: Sets the maximum number of events used per method and charge. Needed to limit memory consumption of fit script; default: 4.e7
Changed:
<
<
As an example, running Prepare script over MC only, Sim08h with 100.000 entries using only T-station method without multiplicity reweighting requires the following command:
>
>
As an example, running Prepare script over MC only, Sim08h with 100000 entries using only T-station method without multiplicity reweighting requires the following command:
  %CODE{ lang="bash" num="off" }% python TrackCalib/scripts/CreateTrackEffTables.py "2015_50ns" -mode "MC" -simVer "Sim08h" -noWeight -method "T" -maxEntries 100 000 --Prepare
Line: 203 to 203
  Fix of the naming convention when plotting the efficiencies in order to produce ready-to-use .C plots.
Added:
>
>
Fix using two methods at the same time.
 
New options:

-fineBinning Creates default tables in the new fine binning

Deleted:
<
<
-autoBinning "<var1>:<NoOfBins>,<binBorderLow2>,<binBorderUp>;<"=: binning scheme used for the customised variables defined as above. Bin borders separated by comma, variable+binning sets by semicolon. Format follows the usual ROOT definition (Int_t nbinsx, Double_t xlow, Double_t xup)
 
Changed:
<
<
-ignoreBins1D "<var1><bin1>-<bin2>,..."=: Sets which variable bins should be set to 1+-5%; default: "". Given bins are set to 1 +- 5%
>
>
-autoBinning "<var1>:<NoOfBins>,<binBorderLow2>,<binBorderUp>;<var2>:<NoOfBins>,<binBorderLow2>,<binBorderUp>": binning scheme used for the customised variables defined as above. Bin borders separated by comma, variable+binning sets by semicolon. Format follows the usual ROOT definition (Int_t nbinsx, Double_t xlow, Double_t xup)

-ignoreBins1D "<var1><bin1>;<var2><bin3>,...": Sets which variable bins should be set to 1+-5%; default: "". Given bins are set to 1 +- 5%

  =-MaxError1D=Sets the maximal acceptable error for the efficiencies; default: 0.03. If the error is >maxError, the efficiency is set to 1 +- 5%
Line: 248 to 251
 

TupleToolTrackEffMatch

Changed:
<
<
Does not affect the efficiency results. This bug effects 2015 and 2016 results. Due to a wrong if, Overlap variables in the tuples can only have values 0 or 1. This means custom cuts on the Overlap are not possible. The cuts used directly in the trigger are OverlapTT, OverlapIT, OverlapMuon > 0.4, OverlapVelo > 0.5.
>
>
Does not affect the efficiency results. This bug effects 2015 and 2016 results. Due to a wrong if, Overlap variables in the tuples can only have values 0 or 1. This means custom cuts on the Overlap are not possible. The cuts used directly in the trigger are OverlapTT, OverlapIT, OverlapMuon > 0.4, OverlapVelo > 0.5.

Bin borders separated by comma, variable+binning sets by semicolon. Format follows the usual ROOT definition (Int_t nbinsx, Double_t xlow, Double_t xup)

Revision 202018-09-11 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 149 to 149
  -vartitles "<var1>:<vartitle1>;<var2>:<vartitle2>": defines the axis title for the custom variables as defined above, standard TLatex synthax can be used
Changed:
<
<
-maxError: Sets the maximal acceptable error for the efficiencies; default: 0.03. If the error is >maxError, the efficiency is set to 1 +- 5%
>
>
-maxError: Sets the maximal acceptable error for the efficiency ratios; default: 0.03. If the error is >maxError, the efficiency ratio is set to 1 +- 5%
  -maxDeviation: Sets the maximal deviation of the efficiency from 1; default: 0.3. If the deviation is >maxDeviation, the efficiency is set to 1 +- 5%
Line: 185 to 185
  ”2016_25ns”: [”Sim09b”] Data (Reco16/Turbo02a) and MC (Sim09b)
Added:
>
>
#Versions

Versions

Version 1.0

First release version.

Version 1.1

Release date: 14.09.2018

Features fixes:

-MaxEntries the new default is -1. This is especially handy when using special cuts. Plotting either Velo of T-station now doesn't crash (found a way around an if needed by TGraphAsymErrors

Fix of the naming convention when plotting the efficiencies in order to produce ready-to-use .C plots.

New options:

-fineBinning Creates default tables in the new fine binning -autoBinning "<var1>:<NoOfBins>,<binBorderLow2>,<binBorderUp>;<"=: binning scheme used for the customised variables defined as above. Bin borders separated by comma, variable+binning sets by semicolon. Format follows the usual ROOT definition (Int_t nbinsx, Double_t xlow, Double_t xup)

-ignoreBins1D "<var1><bin1>-<bin2>,..."=: Sets which variable bins should be set to 1+-5%; default: "". Given bins are set to 1 +- 5%

=-MaxError1D=Sets the maximal acceptable error for the efficiencies; default: 0.03. If the error is >maxError, the efficiency is set to 1 +- 5%

New features:

When fitting, a two log files are produced. First, called trackEff_mode_method_status.out returns the ROOT fit status (should be 0), χ^2 /NDOF of the fit, all parameter values with errors, its the relative errors and final efficiency values. Second, file called WARNINGS_trackEff_mode_method.out is created. If the fit fails, the χ^2 /NDOF>chi2Limit or when the relative error is bigger than relErrLimit, an entry is created in the WARNINGS file. Always check at least the warnings file for problematic fits.

Version 2.0

Release date: TBA

Will include 2017 DATA and MC. Will be released once the WG production of 2017 data and 2017 Sim09D MC is finished.

 #Bugs

Bugs

Added:
>
>
 

2016 TurCal bug

Only affects VELO and T-station methods in 2016. All public plots and correction tables are OK. However, there is a possibility for users to be affected by this bug in special cases as explained below. At the moment the recommendation is to use Long method for 2016.

Line: 206 to 245
 Dangerous variables

Do not cut on/bin in anything that is produced by TupleToolTrackInfo (such as TRACK CHI2NDOF). Overlap information or matching information is not safe neither.

Added:
>
>

TupleToolTrackEffMatch

Does not affect the efficiency results. This bug effects 2015 and 2016 results. Due to a wrong if, Overlap variables in the tuples can only have values 0 or 1. This means custom cuts on the Overlap are not possible. The cuts used directly in the trigger are OverlapTT, OverlapIT, OverlapMuon > 0.4, OverlapVelo > 0.5.

Revision 192018-09-11 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 105 to 105
  -weightVar: reweighting MC to fit the Data multiplicity in the weighting variable; default = "nSPDHits"
Changed:
<
<
-cuts "Probe:<cut1,cut2,cut3>;Tag:<cut4>;Global:<cut6,cut7>": Define additional cuts applied to the probe/tag muon as well as global cuts (e.g. Jpsi cuts, runNumber, etc.). Separate individual cuts by comma, categories by semicolon
>
>
-cuts "Probe:<cut1,cut2,cut3>;Tag:<cut4>;Global:<cut6,cut7>": Define additional cuts applied to the probe/tag muon as well as global cuts (e.g. Jpsi cuts, runNumber, etc.). Separate individual cuts by comma, categories by semicolon. Name convention for particles is muplus, muminus, J_psi_1S, a cut example for all particles having η>2 would be -cuts "Probe:ETA>2;Tag:ETA>2;Global:J_psi_1S_ETA > 2"
  -matchCrit "Probe:<criterion1,criterion2>;Global:<criterion3>": Define tighter matching criterion (when is an event considered efficient), e.g. Matched_GhostProb < 0.3, separation as above

Revision 182018-07-04 - AlbertBursche

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 18 to 18
 %CODE{ lang="bash" num="off" }%

lb-dev Urania/v7r0

Changed:
<
<
cd UraniaDeV _v7r0
>
>
cd UraniaDev_v7r0
 git lb-use Urania git lb-checkout Urania/master TrackCalib

Revision 172018-07-04 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 161 to 161
 python TrackCalib /scripts/CreateTrackEffTables.py "2016_25ns" -variables2D "PT-ETA" -vartitles "PT:P_{T} [MeV]" -ignoreBins "PT-ETA:1-1" -maxError 0.02 --Plot %ENDCODE%
Added:
>
>
There are known issues when trying to plot only Velo or T methods. The fix will be included in the next release.
 

Adding global variables

Line: 184 to 186
 ”2016_25ns”: [”Sim09b”] Data (Reco16/Turbo02a) and MC (Sim09b)

#Bugs

Changed:
<
<

2016 TurCal bug

>
>

Bugs

2016 TurCal bug

  Only affects VELO and T-station methods in 2016. All public plots and correction tables are OK. However, there is a possibility for users to be affected by this bug in special cases as explained below. At the moment the recommendation is to use Long method for 2016.

Revision 162018-06-22 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 140 to 140
 python TrackCalib /scripts/CreateTrackEffTables.py "2016_25ns" -variables "PT,ETA" -binning "PT:0,1000,5000" -binnedFit --Fit %ENDCODE%
Added:
>
>
Sometimes, the unbinned fit requires too much resources and gets killed. In such a case, using -binnedFit: is recommended.
 

'Plot' script

Revision 152018-04-26 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 6 to 6
  TrackCalib was created as tool for creating tracking efficiency correction tables allowing for user-defined track-quality cuts, binning and variable combinations. It also allows to check 1D- and 2D-efficiency dependencies for many variables. The name and idea was shamelessly stolen from PIDCalib.
Changed:
<
<
In case of any questions of comments regarding the TrackCalibTool contact authors michael.kolpin@cernNOSPAMPLEASE.ch and renata.kopecna@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017 and Tutorial talk at LHCb week, 05.12.2017
>
>
In case of any questions of comments regarding the TrackCalibTool contact author renata.kopecna@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017 and Tutorial talk at LHCb week, 05.12.2017
 
Line: 180 to 180
 ”2015_5TeV”: [”Sim09b”] 5TeV Data (Reco15a/Turbo01aEM) and MC (Sim09b)

”2016_25ns”: [”Sim09b”] Data (Reco16/Turbo02a) and MC (Sim09b)

Added:
>
>
#Bugs

2016 TurCal bug

Only affects VELO and T-station methods in 2016. All public plots and correction tables are OK. However, there is a possibility for users to be affected by this bug in special cases as explained below. At the moment the recommendation is to use Long method for 2016.

A bug in Tesla was discovered, that rewrites candidate tracks by long tracks. Details to the bug can be found in this JIRA task: https://its.cern.ch/jira/browse/LHCBPS-1774

This means probe downstream and VELOmuon tracks were replaced by long tracks. Therefore applying any offline cut on variables produced by DaVinci using direct track information is potentially dangerous, since long tracks behave differently when vertex fitting.

Default TrackCalib setting is safe to be used. It uses information directly from the trigger without reprocessing anything. However, several TupleTools access the track information directly. List of variables that are safe and that are potentially dangerous to cut on or bin is can be found below.

In case of doubt, do not hesitate to contact the TrackCalib authors.

Safe variables

It is safe to cut on global event variables, tag mu and Jpsi. It is also safe to use any kinematic variable (produced by TupleToolKinematics, TupleToolGeometry, ETA) or any PID variable.

Dangerous variables

Do not cut on/bin in anything that is produced by TupleToolTrackInfo (such as TRACK CHI2NDOF). Overlap information or matching information is not safe neither.

Revision 142018-04-06 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 126 to 126
  -binning "<var1>:<binborder1>,<binborder2>,<binborder3>;<var2>:<binborder4>,...": binning scheme used for the customised variables defined as above. Bin borders separated by comma, variable+binning sets by semicolon.
Changed:
<
<
-simFit "True/False": controls the use of a simultaneous fit to matched and failed events to evaluate efficiency. Should be set to False for binnings/datasets with low statistics; default: True
>
>
-simFit "True/False": controls the use of a simultaneous fit to matched and failed events to evaluate efficiency. Should be set to False for binnings/datasets with low statistic. default: True

-simpleFit "True/False":Controls the fit shape.Should be set to False for binnings/datasets with low statistic. If True fit shape is 2 CBs (sig) + exp (bkg), ifFalse fit shape is 2 Gaussians (sig) + exp (bkg); default: False

  -binnedFit: Sets the use of a binned fit to evaluate efficiency. Should be set to True for datasets with large statistics; default: False

Revision 132018-02-06 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 6 to 6
  TrackCalib was created as tool for creating tracking efficiency correction tables allowing for user-defined track-quality cuts, binning and variable combinations. It also allows to check 1D- and 2D-efficiency dependencies for many variables. The name and idea was shamelessly stolen from PIDCalib.
Changed:
<
<
In case of any questions of comments regarding the TrackCalibTool contact authors michael.kolpin@cernNOSPAMPLEASE.ch and renata.kopecnan@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017
>
>
In case of any questions of comments regarding the TrackCalibTool contact authors michael.kolpin@cernNOSPAMPLEASE.ch and renata.kopecna@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017 and Tutorial talk at LHCb week, 05.12.2017
 
Line: 51 to 51
  This loops over both MC and Data, all three methods (Long, Velo and T), creates a reduced tuple (with multiplicity weights), fits in the default binning, and create plots of the results. The year_conditions define the dataset used. It is possible to also run this step-by-step or to add other options.
Changed:
<
<
Valid year/conditions strings: are "2015_50ns" (EM), "2015_25ns", "2015_5TeV" and "2016_25ns"
>
>
Valid year/conditions strings: are "2015_50ns" (EM), "2015_25ns", "2015_5TeV" and "2016_25ns". For more information see the detailed list.
 

Running step-by-step

Line: 167 to 167
  totCandidates, nTracks, nLongTracks, nTTTracks, nVeloTracks, nPVs, nSPDHits, nVeloClusters, nITClusters, nOTClusters, nTTClusters, nRich1Hits, nRich2Hits, nMuonTracks, matched, Polarity, runNumber

In case of any doubt contact the authors. \ No newline at end of file

Added:
>
>

Available datasets

”2015_50ns”: [”Sim09b”] EM Data (Reco15em/Turbo01em) and MC (Sim09b)

”2015_25ns”: [”Sim09b”] Data (Reco15a/Turbo01a) and MC (Sim09b)

”2015_5TeV”: [”Sim09b”] 5TeV Data (Reco15a/Turbo01aEM) and MC (Sim09b)

”2016_25ns”: [”Sim09b”] Data (Reco16/Turbo02a) and MC (Sim09b)

Revision 122018-01-15 - LuciaGrillo

Line: 1 to 1
Changed:
<
<
META TOPICPARENT name="LHCbTrackingEfficiencies"
>
>
META TOPICPARENT name="LHCbInternal.LHCbTrackingEfficiencies"
 
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->

TrackCalib Package

TrackCalib was created as tool for creating tracking efficiency correction tables allowing for user-defined track-quality cuts, binning and variable combinations. It also allows to check 1D- and 2D-efficiency dependencies for many variables. The name and idea was shamelessly stolen from PIDCalib.

Changed:
<
<
In case of any questions of comments regarding the TrackCalibTool contact authors michael.kolpin@cernNOSPAMPLEASE.ch and renata.kopecnan@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017
>
>
In case of any questions of comments regarding the TrackCalibTool contact authors michael.kolpin@cernNOSPAMPLEASE.ch and renata.kopecnan@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017
 

Revision 112017-12-04 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 13 to 13
 

TrackCalib setup instructions

Changed:
<
<
TrackCalib is now integrated in Urania. To get Urania, standart setup is used:
>
>
TrackCalib is now integrated in Urania. To get Urania, standard setup is used:
  %CODE{ lang="bash" num="off" }%

lb-dev Urania/v7r0

Changed:
<
<
cd UraniaDev _v7r0
>
>
cd UraniaDeV _v7r0
 git lb-use Urania git lb-checkout Urania/master TrackCalib

Revision 102017-11-20 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->

Revision 92017-11-17 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 114 to 114
 As an example, running Prepare script over MC only, Sim08h with 100.000 entries using only T-station method without multiplicity reweighting requires the following command:

%CODE{ lang="bash" num="off" }%

Changed:
<
<
python TrackCalib/scripts/CreateTrackEffTables.py "2015_50ns" -mode "MC" -simVer "Sim08h" -noWeight -method "T" -maxEntries "100.000 --Prepare
>
>
python TrackCalib /scripts/CreateTrackEffTables.py "2015_50ns" -mode "MC" -simVer "Sim08h" -noWeight -method "T" -maxEntries 100 000 --Prepare
 %ENDCODE%
Changed:
<
<
This allowes for a quicker check of the efficiencies in different MC samples.
>
>
This allows for a quicker check of the efficiencies in different MC samples.
 

'Fit' script

Revision 82017-11-16 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->

Revision 72017-11-13 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->

Revision 62017-11-08 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
Line: 164 to 164
  Standard global variables are:
Changed:
<
<
totCandidates, nTracks, nLongTracks, nTTTracks, nVeloTracks, nPVs, nSPDHits, nVeloClusters, nITClusters, nOTClusters, nTTClusters, nRich1Hits, nRich2Hits, nMuonTracks, matched, Polarity
>
>
totCandidates, nTracks, nLongTracks, nTTTracks, nVeloTracks, nPVs, nSPDHits, nVeloClusters, nITClusters, nOTClusters, nTTClusters, nRich1Hits, nRich2Hits, nMuonTracks, matched, Polarity, runNumber
  In case of any doubt contact the authors.

Revision 52017-11-08 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbTrackingEfficiencies"
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->

TrackCalib Package

Changed:
<
<
TrackCalib wass created as tool for creating tracking efficiency correction tables allowing for user-defined track-quality cuts, binning and variable combinations. It also allows to check 1D- and 2D-efficiency dependencies for many variables. The name and idea was shamelessly stolen from PIDCalib.
>
>
TrackCalib was created as tool for creating tracking efficiency correction tables allowing for user-defined track-quality cuts, binning and variable combinations. It also allows to check 1D- and 2D-efficiency dependencies for many variables. The name and idea was shamelessly stolen from PIDCalib.
 
Changed:
<
<
In case of any questions of comments regarding the TrackCalibTool contact authors michael.kolpin@cernNOSPAMPLEASE.ch and renata.kopecnan@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017
>
>
In case of any questions of comments regarding the TrackCalibTool contact authors michael.kolpin@cernNOSPAMPLEASE.ch and renata.kopecnan@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017
 
Line: 27 to 27
  %ENDCODE%
Changed:
<
<
More detailed instructions can be found in the Git4LHCb twiki.
>
>
More detailed instructions can be found in the Git4LHCb twiki.
 

TrackCalib running instructions

Line: 88 to 88
  -polarity "MagUp/MagDown": specify to run only over MagUp or MagDown. Does not affect the Prepare script
Changed:
<
<
Hence, for example to run over 2016 collision data only, use Long and Velo methods for PT and ETA with GhostProb '&lt' 0.5, run:
>
>
Hence, for example to run over 2016 collision data only, use Long and Velo methods for PT and ETA with GhostProb < 0.5, run:
  %CODE{ lang="bash" num="off" }%
Changed:
<
<
python TrackCalib/scripts/CreateTrackEffTables.py "2016_25ns" -mode "Data" -method "Long, Velo" -variables "PT,ETA" -binning "PT:0,1000,5000" -variables2D "PT-ETA" -matchCrit "Probe: GhostProb<0.5;Global:GhostProb<0.5"
>
>
python TrackCalib /scripts/CreateTrackEffTables.py "2016_25ns" -mode "Data" -method "Long, Velo" -variables "PT,ETA" -binning "PT:0,1000,5000" -variables2D "PT-ETA" -matchCrit "Probe: GhostProb <0.5;Global:GhostProb<0.5"
 %ENDCODE%

Note that the -matchCrit is listed in the Prepare options. This example might be useful to check the effect of looser cut on GhostProb on PT and ETA on data. Note that since PT is not a default variable (P and ETA are), binning in PT needs to be specified. For more information about the binning option see Fit options.

Line: 120 to 121
 

'Fit' script

Changed:
<
<
Executes a simultaneous fit to matched and failed events in each variable and bin to evaluate efficiency. It also creates a separate ROOT file for each variable, containing fitted distributions and efficiency histogram. Output should be checked (at least) for convergence and reasonable data description. Output files are created at results/conditions/(simVer)/trackEffData/MC_variable_method.root
>
>
Executes a simultaneous fit to matched and failed events in each variable and bin to evaluate efficiency. It also creates a separate ROOT file for each variable, containing fitted distributions and efficiency histogram. Output should be checked (at least) for convergence and reasonable data description. Output files are created at results/conditions/(simVer)/trackEffData/MC_variable_method.root
  -binning "<var1>:<binborder1>,<binborder2>,<binborder3>;<var2>:<binborder4>,...": binning scheme used for the customised variables defined as above. Bin borders separated by comma, variable+binning sets by semicolon.

Revision 42017-11-08 - RenataKopecna1

Line: 1 to 1
 
META TOPICPARENT name="LHCbTrackingEfficiencies"
Changed:
<
<
Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples.
>
>
<!--- Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples. Tutorial talk at A&S week, 06.04.2017 -- MichaelKolpin - 2017-04-05--->
 
Changed:
<
<
Tutorial talk at A&S week, 06.04.2017
>
>

TrackCalib Package

 
Deleted:
<
<
-- MichaelKolpin - 2017-04-05
 \ No newline at end of file
Added:
>
>
TrackCalib wass created as tool for creating tracking efficiency correction tables allowing for user-defined track-quality cuts, binning and variable combinations. It also allows to check 1D- and 2D-efficiency dependencies for many variables. The name and idea was shamelessly stolen from PIDCalib.

In case of any questions of comments regarding the TrackCalibTool contact authors michael.kolpin@cernNOSPAMPLEASE.ch and renata.kopecnan@cernNOSPAMPLEASE.ch. More general information regarding tracking efficiencies can be found in the Tracking efficiencies twiki. Also see the Tutorial talk at A&S week, 06.04.2017

TrackCalib setup instructions

TrackCalib is now integrated in Urania. To get Urania, standart setup is used:

<!-- SyntaxHighlightingPlugin -->
lb-dev Urania/v7r0
cd UraniaDev_v7r0
git lb-use Urania
git lb-checkout Urania/master TrackCalib

make configure
make
<!-- end SyntaxHighlightingPlugin -->

More detailed instructions can be found in the Git4LHCb twiki.

TrackCalib running instructions

In order to run the TrackCalib package, a bash shell needs to be loaded

<!-- SyntaxHighlightingPlugin -->
./run bash --norc
<!-- end SyntaxHighlightingPlugin -->

The package is run from the Urania folder.

Minimal preparation script execution:

<!-- SyntaxHighlightingPlugin -->
python TrackCalib/scripts/CreateTrackEffTables.py "2016_25ns"
<!-- end SyntaxHighlightingPlugin -->

This loops over both MC and Data, all three methods (Long, Velo and T), creates a reduced tuple (with multiplicity weights), fits in the default binning, and create plots of the results. The year_conditions define the dataset used. It is possible to also run this step-by-step or to add other options.

Valid year/conditions strings: are "2015_50ns"=(EM), ="2015_25ns", "2015_5TeV" and "2016_25ns"

Running step-by-step

The three sub-scripts can be run separately by specifying:

--Prepare: only runs the Prepare script, reducing the dataset to specified variables, weighting MC in multiplicity, and applying additional matching criteria

--Fit: creates the efficiency fits in a given set of variables, mode, and method with an earlier created tuple

--Plot: Plots the efficiencies from the fit, and creates data/MC correction tables in all 2D variables

Additional options

There are several additional options. Either they are global or they affect the sub-scripts only, as listed below:

Global options

-v: verbose printout

-vROOT: sets ROOT printout level to kInfo, default: kFatal

-mode "MC/Data": specify to only run over Data or MC; default: both

-simVer "SimXYZ": defines the Sim version of the MC; default: Sim09b

-method "<method1>,<method2>": specify method to use, separate by comma; default: Long, T, and Velo

-variables "<var1>,<var2>,<var3>": additional variables for which the efficiency dependency is evaluated; default: P, ETA, nPVs, nSPDHits

-variables2D "<var1>-<var2>,<var3>-<var4>": additional sets of 2D variables for which the efficiency dependency is evaluated. Minus connects the variables within a set, comma separates sets; default: P-ETA

-polarity "MagUp/MagDown": specify to run only over MagUp or MagDown. Does not affect the Prepare script

Hence, for example to run over 2016 collision data only, use Long and Velo methods for PT and ETA with GhostProb '&lt' 0.5, run:

<!-- SyntaxHighlightingPlugin -->
python TrackCalib/scripts/CreateTrackEffTables.py "2016_25ns" -mode "Data" -method "Long, Velo" -variables "PT,ETA" -binning "PT:0,1000,5000" -variables2D "PT-ETA" -matchCrit "Probe: GhostProb<0.5;Global:GhostProb<0.5"
<!-- end SyntaxHighlightingPlugin -->

Note that the -matchCrit is listed in the Prepare options. This example might be useful to check the effect of looser cut on GhostProb on PT and ETA on data. Note that since PT is not a default variable (P and ETA are), binning in PT needs to be specified. For more information about the binning option see Fit options.

'Prepare' script options

Prepare.py runs over the master tuple for the given data taking conditions/Sim version, applies additional selection cuts and reweights MC in nSPDHits. The outputfile is local at tuples/trackEffTuple_Data/MC_conditions(_simVer).root

-noWeight: set to disable reweighting MC to fit the Data multiplicity

-weightVar: reweighting MC to fit the Data multiplicity in the weighting variable; default = "nSPDHits"

-cuts "Probe:<cut1,cut2,cut3>;Tag:<cut4>;Global:<cut6,cut7>": Define additional cuts applied to the probe/tag muon as well as global cuts (e.g. Jpsi cuts, runNumber, etc.). Separate individual cuts by comma, categories by semicolon

-matchCrit "Probe:<criterion1,criterion2>;Global:<criterion3>": Define tighter matching criterion (when is an event considered efficient), e.g. Matched_GhostProb < 0.3, separation as above

-maxEntries 1234567: Sets the maximum number of events used per method and charge. Needed to limit memory consumption of fit script; default: 4.e7

As an example, running Prepare script over MC only, Sim08h with 100.000 entries using only T-station method without multiplicity reweighting requires the following command:

<!-- SyntaxHighlightingPlugin -->
python TrackCalib/scripts/CreateTrackEffTables.py "2015_50ns" -mode "MC" -simVer "Sim08h" -noWeight -method "T" -maxEntries "100.000 --Prepare
<!-- end SyntaxHighlightingPlugin -->

This allowes for a quicker check of the efficiencies in different MC samples.

'Fit' script

Executes a simultaneous fit to matched and failed events in each variable and bin to evaluate efficiency. It also creates a separate ROOT file for each variable, containing fitted distributions and efficiency histogram. Output should be checked (at least) for convergence and reasonable data description. Output files are created at results/conditions/(simVer)/trackEffData/MC_variable_method.root

-binning "<var1>:<binborder1>,<binborder2>,<binborder3>;<var2>:<binborder4>,...": binning scheme used for the customised variables defined as above. Bin borders separated by comma, variable+binning sets by semicolon.

-simFit "True/False": controls the use of a simultaneous fit to matched and failed events to evaluate efficiency. Should be set to False for binnings/datasets with low statistics; default: True

-binnedFit: Sets the use of a binned fit to evaluate efficiency. Should be set to True for datasets with large statistics; default: False

-binsInFit 123: Sets the number of bins in the fit; default: 100

As an example, the following command runs the Fit script for both MC and 2016 data, using Long method for PT and ETA with bins of PT 0MeV-1000MeV-5000MeV and standard ETA bins using binned fits:

<!-- SyntaxHighlightingPlugin -->
python TrackCalib/scripts/CreateTrackEffTables.py "2016_25ns" -variables "PT,ETA" -binning "PT:0,1000,5000" -binnedFit --Fit
<!-- end SyntaxHighlightingPlugin -->

'Plot' script

This script uses the results of previous Fit script, plots all 1D efficiencies and compares data with MC. It also creates the data/MC ratio for 2D variables and saves the ratio as a ROOT file (meaning it creates the correction tables). It also combines the results from Velo and T station methods (called combined). As the final step it combines the combined method with Long method (called final). For 2016 results, the T-station method is not fully functional for now , so use only Long method as an input for your analysis! The plots and tables are stored in plots/conditions_simVer/.

-vartitles "<var1>:<vartitle1>;<var2>:<vartitle2>": defines the axis title for the custom variables as defined above, standard TLatex synthax can be used

-maxError: Sets the maximal acceptable error for the efficiencies; default: 0.03. If the error is >maxError, the efficiency is set to 1 +- 5%

-maxDeviation: Sets the maximal deviation of the efficiency from 1; default: 0.3. If the deviation is >maxDeviation, the efficiency is set to 1 +- 5%

-ignoreBins "<var1>-<var2>:<bin1var1>-<bin1var2>,<bin2var1>-<bin2var2>,...": Sets which variable bins should be set to 1+-5%; default: P-ETA:1-2,5-1. Given bins are set to 1 +- 5%

In order to set the efficiency to 1+-5% in the first bin of PT-ETA as well as in any bin that has the efficiency error larger than 2%, run:

<!-- SyntaxHighlightingPlugin -->
python TrackCalib/scripts/CreateTrackEffTables.py "2016_25ns" -variables2D "PT-ETA" -vartitles "PT:P_{T} [MeV]" -ignoreBins "PT-ETA:1-1" -maxError 0.02 --Plot
<!-- end SyntaxHighlightingPlugin -->

Adding global variables

Note that adding non-standard global variables requires manual changes in Prepare.py.

Standard global variables are:

totCandidates, nTracks, nLongTracks, nTTTracks, nVeloTracks, nPVs, nSPDHits, nVeloClusters, nITClusters, nOTClusters, nTTClusters, nRich1Hits, nRich2Hits, nMuonTracks, matched, Polarity

In case of any doubt contact the authors.

 \ No newline at end of file

Revision 32017-04-07 - MichaelKolpin

Line: 1 to 1
 
META TOPICPARENT name="LHCbTrackingEfficiencies"
Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples.
Changed:
<
<
Tutorial talk at A&S week, 06.03.2017
>
>
Tutorial talk at A&S week, 06.04.2017
  -- MichaelKolpin - 2017-04-05 \ No newline at end of file

Revision 22017-04-06 - MichaelKolpin

Line: 1 to 1
 
META TOPICPARENT name="LHCbTrackingEfficiencies"
Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples.
Added:
>
>
Tutorial talk at A&S week, 06.03.2017
 -- MichaelKolpin - 2017-04-05 \ No newline at end of file

Revision 12017-04-05 - MichaelKolpin

Line: 1 to 1
Added:
>
>
META TOPICPARENT name="LHCbTrackingEfficiencies"
Once created, this page will give instructions on how to use the TrackCalib tool, which options are available, and a list of variables available in the "Master" tuples.

-- MichaelKolpin - 2017-04-05

 
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