University of Birmingham Measurement of the B± production asymmetry and the CP asymmetry in B±J/K± decays

• Users may freely distribute the URL that is used to identify this publication. • Users may download and/or print one copy of the publication from the University of Birmingham research portal for the purpose of private study or non-commercial research. • User may use extracts from the document in line with the concept of ‘fair dealing’ under the Copyright, Designs and Patents Act 1988 (?) • Users may not further distribute the material nor use it for the purposes of commercial gain.


I. INTRODUCTION
One of the primary goals of the LHCb experiment is to search for effects of physics beyond the Standard Model through measurements of CP-violating asymmetries in beauty-and charm-hadron decays.A challenge for such measurements in pp collisions is that the heavy flavor production rates differ between particles and antiparticles.These production asymmetries cannot be precisely predicted since they arise in the nonperturbative b or c quark hadronization process [1][2][3].The effects of production asymmetries cancel in measurements of the difference between CP asymmetries of two different decays of the same hadron species.
The CP asymmetries of B þ meson decay rates 1 are often measured relative to that of the decay B þ → J=ψK þ .The leading tree-level diagram for this decay, shown in Fig. 1 (left), is color suppressed, and the total decay amplitude may receive a sizeable contribution from the gluonic loop diagram shown in Fig. 1 (right).Therefore, the B þ → J=ψK þ decay can in principle exhibit a CP asymmetry due to the interference between these amplitudes.The current world average value of the CP asymmetry is [4], and the uncertainty represents a limitation in many B þ meson CP asymmetry measurements that use this channel as a reference.
This analysis exploits the decay B þ → D0 π þ , which is dominated by a Cabibbo-and color-favored tree-level amplitude and is therefore expected to have a CP asymmetry with a smaller value and uncertainty than for the decay mode is used to measure the production asymmetry between the cross sections for B − and B þ mesons, defined as Since the production asymmetry is expected to be a function of the kinematics, the measurement is performed in nine bins of B þ transverse momentum, p T , and rapidity, y, within the fiducial region 2 < p T < 30 GeV=c and 2.1 < y < 4.5.Measurements are performed on two data sets corresponding to integrated luminosities of 1 and 2 fb −1 , recorded at center-of-mass energies of 7 and 8 TeV in 2011 and 2012, respectively.These measurements complement the existing LHCb studies of heavy flavor production asymmetries [5][6][7][8].A combined analysis of decays allows a measurement of the CP asymmetry in the latter mode.The raw charge asymmetry for a flavor-specific decay to the final state f ( f) accessible in decays of B − (B þ ) mesons is defined as For the two decay modes under study, the asymmetries are well approximated by where A det is the detector-induced asymmetry resulting from differences in the detection efficiencies between particles and antiparticles.All contributions to A det are measured on independent control samples from the same data set.The high correlation of A det between the two decay modes implies a partial cancellation in their difference.This cancellation and the low level of CP violation in the

II. LHCB DETECTOR
The LHCb detector [9,10] is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5, designed for the study of particles containing b or c quarks.The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region, a large-area siliconstrip detector located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of siliconstrip detectors and straw drift tubes placed downstream of the magnet.Data samples corresponding to roughly equal integrated luminosities were recorded with configurations in which the magnetic field was pointing vertically upward and downward.This largely canceled any charge asymmetries in the reconstruction efficiency for charged particles.
The tracking system provides a measurement of momentum, p, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200 GeV=c.The minimum distance of a track to a primary vertex (PV), the impact parameter (IP), is measured with a resolution of ð15 þ 29=p T Þ μm, where p T is the component of the momentum transverse to the beam, in GeV=c.Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors.Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers.The online event selection is performed by a trigger [11], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a two-stage software trigger, which applies a full event reconstruction.This analysis makes use of inclusive dimuon and beauty selections at the software trigger stages.

III. SELECTION OF
decays closely follows a recent LHCb analysis involving the same decay channel [12].Events are considered for the analysis if they contain a track with large enough p T and IP to satisfy the requirements of the first stage of the software trigger.An inclusive beauty selection is applied at the second stage of the software trigger.Candidate D0 → decays are constructed from the intersection of two (four) tracks that satisfy appropriate kaon or pion particle identification (PID) criteria and that have a large p T and significant IP with respect to all primary vertices.These candidates must have a mass within AE25 MeV=c 2 of the D0 mass [4].Each D0 candidate is combined with a high-p T track that is identified as a pion to create a displaced vertex that is consistent with the decay of a B þ meson.The B þ candidates are required to have a mass within the range 5079-5899 MeV=c 2 .To reduce to a negligible level the uncertainty related to L0 trigger asymmetries, it is explicitly required that a positive L0 trigger decision was caused by a particle that is distinct from any of the final-state particles that compose the signal candidate.This requirement is independent of whether or not the signal candidate itself also caused a positive L0 trigger decision and is therefore referred to as triggered independently of signal (TIS) [11].
For both the two-and four-body D0 -mode selections, a pair of boosted decision tree (BDT) discriminators [13], implementing the gradient boost algorithm [14], is used to achieve further background suppression.The first of these BDTs is trained to reject candidates with fake D0 decays, and the second is trained to reject backgrounds with real D0 decays.The BDTs are trained using simulated B þ → D0 π þ signal decays and a sample of decays from data with masses in the range 5900-7200 MeV=c 2 to model the combinatorial background in the nominal mass range.For the training of the first BDT, a background sample is provided by candidates with D0 masses that differ by more than AE30 MeV=c 2 from the known D0 mass.The second BDT is trained using a background sample of candidates with D0 masses within AE25 MeV=c 2 of the known D0 mass.A loose cut on the classifier response of the first BDT is FIG. 1. Tree and loop (penguin) applied prior to training the second one.The inputs into the BDTs include properties of each particle (p, p T , and the IP significance) and additional properties of the B and D 0 composite particles (decay time, flight distance, decay vertex quality, radial distance between the decay vertex and the PV, and the angle between the reconstructed momentum vector and the line connecting the production and decay vertex).A further input into the BDTs is an isolation variable, for which the sum is taken over tracks that are not part of the signal candidate but fall within a cone of half-angle ΔR < 1.5 rad, where ðΔRÞ 2 ¼ ðΔθÞ The cut on the second BDT response is optimized by minimizing the expected uncertainty on the asymmetry between the yields of No PID information is used in the BDT training, but the purity of the sample is further improved by requiring all kaon and pion candidates to satisfy PID criteria.Events containing more than one B þ → D0 π þ candidate amount to less than 1%, and in these cases, the candidate with the highestquality B þ decay vertex is selected.
The raw asymmetries between the yields of B − → D 0 π − and B þ → D0 π þ decays are determined by binned maximum likelihood fits to the mass distributions of selected B − and B þ candidates, treating the two-and four-body D0 modes separately.The fit function is built from a signal component and three background components.A sum of two Gaussian functions with asymmetric power-law tails and an additional Gaussian function are combined to model decays have a distribution that is below the signal peak with a tail that extends to lower masses.They are modeled by the sum of two Gaussian functions with asymmetric power-law tail components.Partially reconstructed decays with an additional particle from a D Ã or ρ meson decay form a background at masses lower than that of the signal peak.This component is described by a combination of analytical functions with shapes that depend on the spin parity of the missing particle, following the method described in Ref. [12].A linear function is adequate to describe the combinatorial background distribution.The yield of misidentified B þ → D0 K þ decays is constrained with an independent control sample of these decays, combined with the calibrated particle identification efficiencies and misidentification rates [15].With the exception of the tail parameters, which are fixed to values obtained from simulation, all parameters are allowed to vary in the fit.
Figure 2 shows the fits to the mass distributions in the bin with 4.5 < p T < 9.5 GeV=c and 2.10 < y < 2.85.The subsequent analysis is based on separate fits for the nine kinematic bins and two center-of-mass energies.The signal yields for each of the nine kinematic bins are listed in Table I.The p T and y intervals of each bin are defined in the second and third columns.The yields sum over B AE meson charges and center-of-mass energies.Integrated over the fiducial acceptance, 2 < p T < 30 GeV=c and 2.1 < y < 4.5, the fits return signal yields of around 2.3 × 10 5 decays for the D0 → K þ π − mode and around 1.3 × 10 5 decays for the based on events in which a muon or a generic track, with large p T and IP significance, satisfies the requirements of the first-stage software trigger.Events must be selected based on a dimuon signature by the second-level software trigger.Candidate J=ψ → μ þ μ − decays are reconstructed from high-p T muon candidates with large IPs with respect to all PVs.A mass interval of 3057-3127 MeV=c 2 is imposed on the J=ψ candidates.These candidates are combined with a high-p T identified kaon with a significant IP with respect to all PVs, where the J=ψ candidate invariant mass is constrained to its known value in the combination.The L0 trigger TIS requirement is applied in the same way as for the D0 π þ selection.A single BDT classifier is used to improve the purity of the B þ → J=ψK þ sample.This classifier is trained on a similar set of variables as that for the B þ → D0 π þ selection and exhibits very similar performance in terms of signal efficiency and background rejection.Events containing more than one B þ → J=ψK þ candidate amount to less than 1%, and in these cases, the candidate with the highest-quality B þ decay vertex is selected.
A simultaneous fit of the mass distributions across the kinematic bins is performed, where the same value of is assumed for all bins.The signal peak is described using a Gaussian function with an additional asymmetric power-law tail component.The mean of the Gaussian is constrained to be the same in all kinematic bins, while its width and the tail parameters are allowed to vary between bins.A small background from misidentified B þ → J=ψπ þ decays is described by a similar function, with fixed shape parameters taken from simulation.The yield of this contribution is allowed to vary in each kinematic bin, but a single raw asymmetry is shared between all bins.The contribution from random particle combinations is described by a linear function.The yield of this component and the slope parameter are allowed to vary in each kinematic bin.The yield is also fitted separately for each B AE charge.
Integrated over the full fiducial acceptance, a signal yield of about 2.3 × 10 5 events is measured.Table I lists the yields of each signal decay mode in each of the kinematic bins summing over the two center-of-mass energies.An example of the fit in the bin with 4.5 < p T < 9.5 GeV=c and 2.10 < y < 2.85 is displayed in Fig. 3.

V. MEASUREMENT OF THE B þ PRODUCTION ASYMMETRY
The B þ production asymmetry is determined in the nine bins of p T and y according to where A Dπ raw and A Dπ CP are the raw charge asymmetry and CP asymmetry in the B þ → D0 π þ decay, respectively.The four A det terms correct for detector-induced asymmetries and will be described in the following.All terms other than A Dπ CP are evaluated separately for the four disjoint data sets corresponding to the two center-of-mass energies and the two magnet polarities.An average of the A Dπ raw values for the two D0 decay modes is computed with weights that are chosen to minimize the uncertainty.The same weights are used to compute averages over the two D0 decay modes for TABLE I.The p T and y intervals for each kinematic bin and the corresponding signal yields in each of the B þ decay modes, summing over the two center-of-mass energies.all other terms in Eq. ( 5) apart from A TIS det , which is independent of the D0 decay.Tables II and III list the values of the first five terms in Eq. ( 5) for the 7 and 8 TeV data sets, respectively.The overall detection asymmetry has two main contributions.The first arises because K − mesons have a larger nuclear interaction cross section than K þ mesons.This means that more K − mesons than K þ mesons interact inelastically with the detector material before they leave enough hits to be reconstructed in the tracking stations.The resulting K − -K þ detection asymmetry is around 10 −2 .The second cause of asymmetry is the different trajectories of positively and negatively charged particles, which therefore have different sensitivities to misalignments and inhomogeneities of the detector.This source contributes to all detection asymmetry terms.It is partially cancelled when averaging measurements over data recorded with the dipole magnet in the two polarities.
The D0 detection asymmetry, A Kπ det , is measured using samples of D − mesons that are produced in the primary pp interactions and decay to the K þ π − π − and K 0 S π − final TABLE II.A summary of the terms that enter the production asymmetry determination [Eq.( 5)] in the 7 TeV data set.The p T and y intervals of each bin are provided in Table I.The L0 trigger asymmetry A TIS det is omitted from this table since it is assumed to be independent of the B þ kinematics.All uncertainties are statistical.

Bin
A Dπ raw (×10 −2 )  III.A summary of the terms that enter the production asymmetry determination [Eq.( 5)] in the 8 TeV data set.The L0 trigger asymmetry A TIS det is omitted from this table since it is assumed to be independent of the B þ kinematics.All uncertainties are statistical.candidates in the bin with 4.5 < p T < 9.5 GeV=c and 2.10 < y < 2.85.These distributions sum over the two center-of-mass energies.B − candidates are displayed on the left, and B þ candidates are on the right.The signal components are displayed as red dotted lines, while the background from combinatorial events is shown by the green dashed lines.
The fit in other kinematic bins is similar, aside from the specific signal and background component yields.
states.The K 0 S mesons are reconstructed in their decay to π þ π − .Within a small phase-space region in terms of the D − decay products, it is assumed that the detection asymmetry for a K þ π − pair can be determined using with a small correction for the effects of CP violation in K 0 − K0 mixing and the different material interactions of K 0 and K0 .For each of the such that the kinematic distributions of the K þ and the matched π − agree with those from the signal D0 decays.For the sample, the procedure is repeated for each of the two possible pions with opposite charge to the kaon, averaging over the two.Each D − → K 0 S π − candidate is assigned a weight, such that the π − kinematic distributions agree with those of the unmatched π − in the weighted D − → K þ π − π − sample, and the D − kinematic distributions are equalized between the two D − decay modes.This ensures cancellation of the D − production asymmetry and means that any detection asymmetry associated with the unmatched π − is cancelled with a corresponding asymmetry affecting the D − → K 0 S π − sample.This weighting procedure is performed for each of the nine B þ kinematic bins.The raw asymmetries that enter Eq. ( 6) are determined by fitting the weighted mass spectra for the four combinations of D AE decay modes and charges.
Using a detailed description of the LHCb detector and cross section measurements from fixed target experiments [4], the nuclear interaction contribution to the pion asymmetry is estimated to be negligibly small.The tracking asymmetry can therefore be assumed to be the same for pions and muons.The π þ tracking asymmetry, A π det , is therefore inferred from that of muons measured using a sample of J=ψ → μ þ μ − decays in which one of the muons is reconstructed without requiring hits in all tracking stations [16].Weights are assigned to the J=ψ candidates such that the kinematic distributions of this muon match those of the The PID requirements on the B þ → D0 π þ decays can introduce asymmetries.Corrections are determined using a control sample of D Ãþ → D 0 π þ decays, with D 0 → K − π þ , in which no PID requirements are imposed on the K − or π þ from the D 0 decay.The asymmetry associated with PID requirements on the D0 decays is partially accounted for in the A Kπ det correction, since PID requirements are imposed on the final-state kaons and pions in the D − control samples.The requirements are tighter in these control samples, and so a residual correction must still be applied.The sum of this correction, and a corresponding correction for the PID requirement on the π þ from the D Ãþ → D 0 π þ decays, is denoted A PID det .The asymmetry associated with the TIS trigger efficiency, A TIS det , is determined using a sample of b-hadron decays to the final state D0 μ þ ν μ X with D0 → K þ π − .An unbiased probe of the TIS trigger efficiency is provided by the subset of these in which the muon prompted a positive decision by the L0 muon trigger.The corresponding asymmetries do not exhibit any kinematic dependence, and so a single correction is determined for each center-ofmass energy and is applied to all kinematic bins.The measured A TIS det values are ðþ0.16AE 0.16Þ × 10 −2 and ðþ0.02AE 0.10Þ × 10 −2 for the 7 and 8 TeV data sets, respectively.
The CP asymmetry, A Dπ CP , is estimated from measurements of the angle γ of the unitarity triangle of the Cabibbo-Kobayashi-Maskawa matrix [17,18] and the hadronic parameters of B þ → D0 π þ decays [19].Different values are obtained for the D0 → decay modes due to the smaller coherence factor from the competing hadronic resonances in the four-body mode.The asymmetries are with a 55% correlation between the uncertainties on these two quantities.The A Dπ CP values reported in Tables II and III are averaged over the two-and four-body modes.These values vary between the kinematic bins due to the different weights of the two-and four-body modes.
Several sources of systematic uncertainty arise in the determination of the production asymmetries.Their contributions are listed in Table IV.Variations in the weighting procedures that are used to determine A Kπ det and A π det yield uncertainties of 0.07 × 10 −2 and 0.04 × 10 −2 , respectively.An uncertainty of 0.04 × 10 −2 is assigned to a possible pion nuclear interaction asymmetry that is not accounted for in the tracking efficiency measurements with muons from J=ψ decays.Finally, the A Dπ CP uncertainties are included in the total systematic uncertainty, which is taken to be correlated between the kinematic bins.

TABLE IV. Systematic uncertainties on the
CP uncertainty varies between the kinematic bins, and the range is indicated.All systematic uncertainties are considered to be correlated between kinematic bins.

Source
Size (×10 −2 ) A Kπ det method AE0.07A π det method AE0.04 Pion nuclear interactions AE0.04A Dπ CP AEð0.04 − 0.05Þ The measured A prod ðB þ Þ values for each kinematic bin are listed in Table V for both center-of-mass energies.They are shown as a function of rapidity for the three p T ranges in Fig. 4. Samples of simulated B AE decays are produced using PYTHIA 8 [20,21] with a specific LHCb configuration [22] and are used to determine the weights that are assigned to each of the nine bins, such that the sum corresponds to the asymmetry integrated over the full fiducial region covering 2 < p T < 30 GeV=c and 2.1 < y < 4.5.These weights are listed in Table V.The integrated asymmetries, which are also reported in Table V, are ffiffi ffi s p ¼ 7 TeVÞ ¼ ð−0.41 AE 0.49 AE 0.10Þ × 10 −2 ; where the first uncertainty is statistical and includes contributions from A Dπ raw and the detection asymmetry corrections which are inherently statistical in nature.The second uncertainty is systematic.Several cross-checks are performed.The measured value of A prod ðB þ Þ is found to have no statistically significant dependence on the B þ decay time or kaon momentum.Statistically compatible results are obtained for the two magnet polarities.

VI. MEASUREMENT OF
where A ψK raw is the raw asymmetry of B AE → J=ψK AE decays and δA Kπ det corrects for the different detection asymmetries of the two decay modes.The two final states differ by the transformation of a π þ π − pair to a μ þ μ − pair, where the only significant contribution to the difference between the overall detection asymmetries arises from the charged kaon asymmetry.The method used to determine A Kπ det , as described in the previous section, is applied to the J=ψK þ final state by considering the muon with opposite charge to the kaon as a pion.The difference between this and the corresponding asymmetry for the B þ → D0 π þ mode is defined as δA Kπ det ¼ A Kπ det ðB → J=ψKÞ − A Kπ det ðB → D0 πÞ.The uncertainties are cancelled to a large degree in this difference.TABLE V.The measured A prod values for each kinematic bin and integrated over the full kinematic acceptance, 2 < p T < 30 GeV=c and 2.1 < y < 4.5.The integrated values sum over the asymmetries in each bin, weighted by the values, w, in the second and fourth columns for the two center-of-mass energies.The first uncertainty is the statistical uncertainty on A Dπ raw and is uncorrelated between the kinematic bins.The second uncertainty is the statistical uncertainty on the detection asymmetry corrections and is taken to be correlated between the kinematic bins.The third uncertainty is purely systematic and is assumed to be correlated between bins.

Bin
wð7 TeVÞ A prod ðB þ ; 7 TeVÞ (×10  det are positive, since the kaons in the J=ψK þ decays tend to have higher momenta than those in the B þ → D0 π þ decays.A further asymmetry could result from differences between the kinematic distributions of the pion in the B þ → D0 π þ decay compared to the μ þ in the J=ψK þ decay, but this is estimated to be negligibly small. The values of A ψK raw in each bin are corrected according to Eq. ( 7) using measurements of A Dπ raw , δA Kπ det , and A Dπ CP in order to extract Gaussian constraints are applied to the values of A Dπ raw and δA Kπ det , such that the statistical uncertainty on these parameters is included in the overall statistical uncertainty for A systematic uncertainty of 0.02 × 10 −2 is assigned for the use of fixed parameters in the mass fits, while a systematic uncertainty of 0.05 × 10 −2 is assigned for the method used to measure δA Kπ det .The A Dπ CP values contribute a systematic uncertainty of 0.04 × 10 −2 .The final result is where the first uncertainty is statistical and the second is systematic.By fixing all Gaussian constrained parameters to have zero uncertainty, the contribution from the finite B þ → J=ψK þ statistics is found to be AE0.20 × 10 −2 .This result is consistent with, and improves upon, the current world average value of

VII. SUMMARY AND CONCLUSIONS
The B þ meson production asymmetry is a crucial input in the measurement of CP asymmetries in B þ decays.A sample of B þ → D0 π þ decays is used to measure the production asymmetry.The analyzed data set corresponds to integrated luminosities of 1 and 2 fb −1 recorded during 2011 and 2012 at proton-proton center-of-mass energies of 7 and 8 TeV, respectively.The production asymmetries are measured in nine bins of transverse momenta and rapidity, covering the region 2 < p T < 30 GeV=c and 2.1 < y < 4.5, and separately for the two center-of-mass energies.The measurements are generally consistent with zero asymmetry within typical uncertainties of roughly 10 −2 , which is in agreement with b -quark hadronization models [1][2][3].Integrated over the full p T and y ranges, the production asymmetries are measured to be where the first uncertainty accounts for all statistical sources and the second accounts for all systematic sources.A simultaneous study of the B þ → J=ψK þ decay allows a measurement of its CP asymmetry,

FIG. 2 .
FIG. 2. Mass distributions of selected (top) B AE → ½K AE π ∓ D π AE and (bottom) B AE → ½K AE π AE π ∓ π ∓ D π AE candidates in the bin with 4.5 < p T < 9.5 GeV=c and 2.10 < y < 2.85.These distributions sum over the two center-of-mass energies.B − candidates are displayed on the left, and B þ candidates are on the right.The red dotted lines indicate the contribution from B AE → Dπ AE decays.The purple dashed-dotted lines indicate the contribution from misidentified B AE → DK AE decays.The gray shaded regions at low values of reconstructed mass indicate the contribution from various partially reconstructed B decays, and the green dashed lines indicate the combinatorial background.The total fit function is shown by the blue solid lines.The fit in other kinematic bins is similar, aside from the specific signal and background component yields.
2þ ðΔϕÞ 2 , and Δθ and Δϕ are the differences in the polar and azimuthal angles of each track with respect to the B þ candidate direction.Tracks are only considered in the isolation cone if they are associated, by the smallest IP, to the same primary vertex as the signal candidate.Signal decays are expected to have larger values of I p T than background.
The measured A prod ðB þ Þ as a function of rapidity of the B meson in three bins of p T .The ranges of p T are indicated in the legends.The left-and right-hand figures correspond to 7 and 8 TeV center-of-mass energies, respectively.

TABLE VI .
Residual differences δA Kπ det , measured in each bin of B kinematics.These are the effective values after summing over center-of-mass energies and averaging over the two D0 decay modes.Also at P.N.Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia.c Also at Università di Milano Bicocca, Milano, Italy.d Also at Università di Modena e Reggio Emilia, Modena, Italy.e Also at Novosibirsk State University, Novosibirsk, Russia.f Also at LIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain.g Also at Università di Bologna, Bologna, Italy.h Also at Università di Roma Tor Vergata, Roma, Italy.i Also at Università di Genova, Genova, Italy.j Also at Scuola Normale Superiore, Pisa, Italy.k Also at Università di Cagliari, Cagliari, Italy.l Also at Università di Padova, Padova, Italy.m Also at Università di Bari, Bari, Italy.n Also at Laboratoire Leprince-Ringuet, Palaiseau, France.o Also at Università degli Studi di Milano, Milano, Italy.p Also at Universidade Federal do Triângulo Mineiro (UFTM), Uberaba-MG, Brazil.q Also at AGH -University of Science and Technology, Faculty of Computer Science, Electronics and Telecommunications, Kraków, Poland.r Also at Iligan Institute of Technology (IIT), Iligan, Philippines.s Also at Hanoi University of Science, Hanoi, Viet Nam.t Also at Università di Pisa, Pisa, Italy.u Also at Università di Roma La Sapienza, Roma, Italy.v Also at Università della Basilicata, Potenza, Italy.w Also at Università di Urbino, Urbino, Italy.