Options for the efficiency algorithm and how to set them
they refer to DaVinciEff v2r8 and higher ( DaVinci v9r0, DaVinci v8r0 compatible)
1. Options relative to the MC truth decays
1.1 You must set the decay you are interested in for the MC
truth part. If you do not do it the programs will stop
at initialization.
EffXXX.MCDecay = "description of the decay according to
MCDecayFinder tool";
( type std::string, def="" )
A description of the syntax can be found in the doxygen
detailed description of the MCDecayFinderTool
A list of decay descriptors for the signal files generated
can be found out of the physics web page
1.2 You must "tag" (putting the character ^ before the
specified particles) the end decay particles in order
to have the rec'tible/rec'ted numbers (see definition).
The executable will stop after the first event if this
is not done.
For example for Bs -> J/Psi (mumu) phi (KK) with possible
radiative gamma where the stable particles are the muons and
kaons you would write:
EffXXX.MCDecay = "B_s0 -> (J/psi(1S) -> ^mu+ ^mu- {, gamma}) (phi(1020) -> ^K+ ^K-)"
(note that the "tag" in not done in the decay descriptors
provided in the list)
1.3 You can specify one parent particle type for which the decays
has to originate, meaning that for the decay described in
the EffXXX.MCDecay only those that will come from the particle
specified here will be considered. This cannot be the head of the
tree as specified in MCDecay because a particle cannot
have itself as a parent.
The default value = "Any", implies that the "parent" of the
decay of interest is not checked.
EffXXX.MCDoIFromParent = "particle name from which the decay
must come from";
( type std::string, def="Any" )
This can be used to look at partial decays of a B, for
example for Ks that come from a B0 in which case one would
write:
EffXXX.MCDecay = "Ks -> ^pi+ ^pi-";
EffXXX.MCDoIFromParent = "B0";
2. Options relative to the selected decay
2.1 You must set the decay you are interested in for the selection
part. At the moment no protection are implemented and it is
directly delegated to the DecayFinder tool. This will be set
in a future release similarly to the MC part.
EffXXX.DecayFinder.Decay = "description of the decay according to
DecayFinder tool";
( type std::string, def="B0 -> pi+ pi-" )
2.2 You must set where in the Transient Store the head of the selected
decay is. Note that at the moment only a single location is used.
EffXXX.SelLocation = "TES path where head particle for selected
decay is located";
( type std::string, def="" )
2.3 You can set a mass window for which separate numbers of selected and
associated numbers will be reported.
EffXXX.MotherMassLow = value of lower limit of mass window;
( type double, def = 4.0 GeV )
EffXXX.MotherMassUpp = value of upper limit of mass window;
( type double, def = 6.0 GeV )
3. Options relative to the associated decay
To do the association the efficiency uses two associators at the
same time (the Chi2 and the Composite Link). Note that at the
moment for some decays with intermediate optional resonances only
the Chi2 gives sensible results.
Property of the associators need to be set.
Both the Chi2 and the Composite Link are forwarding where to find
the input particles to the necessary algorithms.
Hence only their "InputData" has to be set.
You must set the InputData for both the composite link associator
and for the chi2 associators.
As an example for Bd->J/Psi(mumu)Ks the necessary options for the
associators are:
// Associator with Links
EffBd2Jpsi2MuMu_Ks2PiPi.CompP2MCAsct.Location = "Phys/Relations/CompPart2MCfromLinks";
EffBd2Jpsi2MuMu_Ks2PiPi.CompP2MCAsct.AlgorithmType = "CompositeParticle2MCLinks";
EffBd2Jpsi2MuMu_Ks2PiPi.CompP2MCAsct.AlgorithmName = "CompositeParticle2MCLinks";
CompositeParticle2MCLinks.InputData = { "Phys/JPsi/Particles",
"Phys/Kshort/Particles",
"Phys/B0/Particles" };
// Associator with Chi2
EffBd2Jpsi2MuMu_Ks2PiPi.Chi2P2MCAsct.Location = "Phys/Relations/Particle2MC";
EffBd2Jpsi2MuMu_Ks2PiPi.Chi2P2MCAsct.AlgorithmType = "Particle2MCChi2";
EffBd2Jpsi2MuMu_Ks2PiPi.Chi2P2MCAsct.AlgorithmName = "Particle2MCChi2";
Particle2MCChi2.InputData = { "Phys/B0/Particles" };
Note that:
- for the Composite link associator you need to specify all the TES
locations where particles contributing to your decay are located
- for the Chi2 associator it is sufficient to specify the TES location
where the head of the decay tree is.