HMdcDeDx

http://www-hades.gsi.de/docs/hydra/classDocumentation/dev/HMdcDeDx.html

http://www-hades.gsi.de/docs/hydra/classDocumentation/dev/HMdcDeDxMaker.html

http://www-hades.gsi.de/docs/hydra/classDocumentation/dev/HMdcDeDxHists.html




 

   

HMdcDeDx

This transformation class calculates the "pseudo dEdx" from t2-t1 (time above threshold) of all fired drift cells included in one HMdcSeg. The transformation is performed doing a normalization of the measured t2-t1 with impact angle and minimum distance to wire (MDCCAL2 parameterization) and the algorithm to normalize the single measurements and average over all cells included in the segment.

Calculation of dEdx:
Float_t calcDeDx(*HMdcSeg ,Float_t,Float_t*,UChar_t*,Float_t*,UChar_t*,Int_t,Int_t)
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 increasing order, the mean and sigma is calculated. Afterwards the truncated 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 (default): Input filling from HMdcSeg (wires, thrown out by the fitter are ignored)
vers==1 : Input filling from HMdcClusFit (wires, thrown out by the fitter are taken into account)
If the fit has failed (chi2=-1), vers 0 is used
mod==0 : calc dEdx from 1st module in segment
mod==1 : calc dEdx from 2nd module in segment
mod==2 (default) : calc dEdx from whole segment

Settings of the truncated Mean method:
The truncated mean is calculated according to the set window (in sigma units)
(setWindow( Float_t window)) around the mean. For the calculation of the truncated mean only drift cells with a dEdx inside the window around the mean (calculated from all drift cells of the segment) will be taken into account.
With setMinimumWires(Int_t minwires) the algorithm can be forced to keep a minimum number of wires. The method will stop removing drift cells from the active sample if the minimum number is reached.
 

 



 

    HMdcDeDxMaker

This transformation class calculates the pseudo dEdx from t2-t1 (time above threshold) of all fired drift cells included in one HMdcSeg. The transformation is performed using the parameter container HMdcDeDx, which holds the parameters for normalization of the measured t2-t1 with impact angle and minimum distance to wire (CAL2 parameterization) and the algorithm to normalize the single measurements and average over all cells included in the segment. Inside the execute() function a loop over the category of HMdcTrkCand is performed and the dEdx calculation is done for each inner and outer HMdcSeg separately. The result is filled to HMdcTrkCand.
OPTIONS:
setDeDxMethod(Int_t vers) :
vers=0 (default) calculate from HMdcSeg
vers=1 calculate from HMdcClusFit
Explanations can be looked up in HMdcDeDx.

 
 



 

    HMdcDeDxHists

Fills momentum vs MDC dEdx ( HSplineTrack or HKickTrack ) for full tracks (seg1+seg2), seg1, seg2 , system 0/1 and sector wise.
OPTIONS:
Momentum version:-----------------------------------------------------------
The histograms can be filled by HKicktrack or HSplineTrack:
The selection can be done with
setUseSpline(Bool_t) :
kTRUE (default) - Fills from HSplineTrack
: kFALSE - Fills from HKickTrack
dEdx calculation:-----------------------------------------------------------
The mode of HMdcDeDx::calcDeDx(......., mode,module) can be switched
via setDeDxMode(Int_t m) : 0 (default) - Fills from HMdcSeg (only cells kept by the fitter)
1 - Fills from HMdcClusFit (cells with weight 0 AND 1 are used)
via setDeDxModule(Int_t m): 0 - Fills for 1st module in HMdcSeg
1 - Fills for 2nd module in HMdcSeg
2 (default) - Fills for both modules in HMdcSeg

via setDeDxWindow(Float_t w) the window of the truncated mean can be set
in units of sigma (default 5.0);
via setDeDxMinimumWires(Int_t n) the minimal number of wires after the truncated
mean procedure can be set.
via setDeDxUseShifts(Bool_t shift) the normalization procedure of the dEdx calculation
can be switched on or off (default=ON).
CPR mode:-------------------------------------------------------------------
Single wire distribution or ntuple filled for single or doubles of Jaro's CPR condition
Runs on special prepared input file (Stripped single or doubles with HKickTrack containing
1 for singles and 2 for doubles )
setCPR(Bool_t cpr) : kTRUE run in the special mode
setCPR(Bool_t cpr) : kFALSE (default) run in the normal mode
DEBUG mode:-----------------------------------------------------------------
If Debug mode is switched on (default OFF),
HMdcTrkCand.getFirstCandInd() =p*charge;
HMdcTrkCand.getNCandForSeg1() =single_double (only if cpr, 1=single,2=double);
The out put dEdx of HMdcTrkCand will be refilled with the calculated values.
setDebug(Bool_t d)
 
 




-- JochenMarkert - 18 Feb 2004

I Attachment ActionSorted ascending Size Date Who Comment
hmdcdedx.htmlhtml hmdcdedx.html manage 2 K 2004-02-26 - 15:15 JochenMarkert hmdcdedx
hmdcdedxhists.htmlhtml hmdcdedxhists.html manage 3 K 2004-02-26 - 15:16 JochenMarkert hmdcdedxhists
hmdcdedxmaker.htmlhtml hmdcdedxmaker.html manage 1 K 2004-02-26 - 15:16 JochenMarkert hmdcdedxmaker
Topic revision: r9 - 2004-03-03, JochenMarkert
Copyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding Foswiki Send feedback | Imprint | Privacy Policy (in German)