class HMdcDeDx2: public HParCond

 destructor

 prints the parameters of HMdcDeDx2 to the screen.

 intitializes the container from an input

 create the maximum allowed value of ToT to keep
 the numerical calulations inside the range of double.
 The procedure loops over all s,m,abin,dbin to set
 the boundary automatically. This function needs
 the parameter array for the fit functions to be initialized
 before. The boundary value are stored inside the container.

 writes the container to an output

 Puts all params of HMdcDeDx2 to the parameter list of
 HParamList (which ist used by the io);

 Puts the measurement values into decreasing order.

 calculates dedx of a MDC segment by doing normalization for
 impact angle and distance from wire for the fired cells of
 the segment. The measurements are sorted in decreasing order,
 the mean and sigma is calculated. Afterwards the truncated mean
 according to the set window (in sigma units) arround the mean
 is calculated and returned as dedx. Mean,sigma,number of wires before
 truncating,truncated mean,truncated mean sigma and number of wires
 after truncating can be obtained by the returned pointers.
 Different methods can be selected:
 vers==0 : Input filling from inner HMdcSeg
 vers==1 : Input filling from outer HMdcSeg
 vers==2 : Input filling from inner+outer HMdcSeg (default)
 mod==0  : calc dedx from 1st module in segment
 mod==1  : calc dedx from 2nd module in segment
 mod==2  : calc dedx from whole segment (default)
 useTruncMean: kTRUE (default) apply truncated Mean
 truncMeanWindow (unit: SIGMA RMS of mean TOT): -99 (default) use standard window
 inputselect==0 : fill from segment (default, wires rejected by fit are missing)
 inputselect==1 : fill from cluster
 returns -99 if nothing was calculated

 Fills array of measurements from HMdcSeg and returns the number of fired wires in Segment
 The single measurements are normalized by impact angle and distance from wire. Only the first
 wires < MAX_WIRES are accepted.
 inputselect==0 : fill from segment (default,wires rejected by fit are missing)
 inputselect==1 : fill from cluster

 loops over the measurement array with nWire values and puts values to
 -99 if the measurements are outside the wanted window arround the mean value mean
 (set with setWindow() (sigma units)).
 The procedure keeps at least a number of minimum wires set with setMinimumWires().

 set the values for the dedx functions by s, m, angle bin, dist bin. Needs pointer
 to parameter array[N_PARAM]

 set all values for the dedx functions. Needs pointer
 to parameter array[6][4][N_ANGLE][N_DIST][N_PARAM]

 set all values for the dedx functions. Needs pointer
 to parameter array[6][4][N_ANGLE][N_DIST][N_PARAM]

 set the values for the dedx functions by s, m, angle bin, dist bin.

 set all values for the dedx functions. Needs pointer
 to parameter array[6][4][N_ANGLE][N_DIST]

 set all values for the dedx functions. Needs pointer
 to parameter array[6][4][N_ANGLE][N_DIST]

 set the values fpr the dedx width functions by s, m, angle bin, dist bin. Needs pointer
 to parameter array[N_SHIFT_PARAM]

 set all values for the dedx width functions. Needs pointer
 to parameter array[6][4][N_ANGLE][N_DIST][N_SHIFT_PARAM]

 set all values for the dedx width functions. Needs pointer
 to parameter array[6][4][N_ANGLE][N_DIST][N_SHIFT_PARAM]

 calculates 2 coordinate points from segment

 returns asymmetric width of ToT for single drift cell
 shift == 0 (default) : returns 0
       ==-1           : low   bound width of ToT distribution (sigma)
       == 1           : upper bound width of ToT distribution (sigma)

 calculation of ToT -> dEdX for single drift cell

 calculation of dEdX -> ToT for single drift cell

 calibrate ToT :
 dEdX=TotTodEdX(t2-t1) ----> dEdXToToT(dEdX) for reference module,angle,distbin

 Calculates the dEdx (MeV 1/g cm2) of an particle with GEANT ID id
 and momentum p (MeV) for He/i-butan gas mixture with He fraction hefr
 (He (hefr) / i-butan (1-hefr)). The fomular is taken from PDG and doesn't
 include the density correction term. The values for the mean excitation
 energy are taken from Sauli.

 creates a TGraph for particle with GEANT ID id and
 and He/i-butan gas mixture with He fraction hefr
 (He (hefr) / i-butan (1-hefr) dEdx vs p,beta,1/beta2 or betagamma
 depending on opt (p,beta,1/beta2,betagamma). exchange=kTRUE
 will exchange x-axis and y-axis.

 calculated the t2 of drift cell measurent
 from dedx of a given particle with momentum p.
 The assumption is that all drift cell measurements
 belonging to one track can be treated independent.
 return t2 and smeared error to pointer t2err for an
 single measurement. If the result t2-t1 would
 be < 0 a random value for t2 0-20 ns larger than
 t1(already including error) will be created.
 The smeared error will be returned to the pointer t2err.

set user values

energy loss