fDebugStreamer(0),
fRecoParam(0)
{
+ //
+ // COSNTRUCTOR
+ // param - tpc parameters for given file
+ // recoparam - reconstruction parameters
+ //
fIsOldRCUFormat = kFALSE;
fInput =0;
fOutput=0;
void AliTPCclustererMI::MakeCluster(Int_t k,Int_t max,Float_t *bins, UInt_t /*m*/,
AliTPCclusterMI &c)
{
+ //
+ // k - Make cluster at position k
+ // bins - 2 D array of signals mapped to 1 dimensional array -
+ // max - the number of time bins er one dimension
+ // c - refernce to cluster to be filled
+ //
Int_t i0=k/max; //central pad
Int_t j0=k%max; //central time bin
Int_t cemaxpos= 0;
Double_t ceThreshold=5.*cerms;
Double_t ceSumThreshold=8.*cerms;
- const Int_t cemin=5; // range for the analysis of the ce signal +- channels from the peak
- const Int_t cemax=5;
+ const Int_t kCemin=5; // range for the analysis of the ce signal +- channels from the peak
+ const Int_t kCemax=5;
for (Int_t i=nchannels-2; i>1; i--){
if ( (dsignal[i]-mean06)>ceThreshold && dsignal[i]>=dsignal[i+1] && dsignal[i]>=dsignal[i-1] ){
cemaxpos=i;
}
}
if (cemaxpos!=0){
- for (Int_t i=cemaxpos-cemin; i<cemaxpos+cemax; i++){
+ for (Int_t i=cemaxpos-kCemin; i<cemaxpos+kCemax; i++){
if (i>0 && i<nchannels&&dsignal[i]- cemean>0){
Double_t val=dsignal[i]- cemean;
ceTime+=val*dtime[i];
virtual void SetInput(TTree * tree); // set input tree with digits
virtual void SetOutput(TTree * tree); //set output tree with
private:
- Bool_t IsMaximum(Float_t k, Int_t max, const Float_t *bins);
+ Bool_t IsMaximum(Float_t k, Int_t max, const Float_t *bins) const;
void MakeCluster2(Int_t k,Int_t max,Float_t *bins,UInt_t m,
AliTPCclusterMI &c);
void MakeCluster(Int_t k,Int_t max,Float_t *bins,UInt_t m,
Float_t * fBins; //!digits array
Float_t * fResBins; //!digits array with res. after 1 finder
Int_t fLoop; //loop - cf in 2 loops
- Int_t fMaxBin;
- Int_t fMaxTime;
- Int_t fMaxPad;
+ Int_t fMaxBin; //current ( for current sector) maximal bin
+ Int_t fMaxTime; //current ( for current sector) maximal time
+ Int_t fMaxPad; //current ( for current sector) maximal pad
Int_t fSector; //!current sector
Int_t fRow; //!current row
Float_t fSign; //!current sign
ClassDef(AliTPCclustererMI,1) // Time Projection Chamber digits
};
-inline Bool_t AliTPCclustererMI::IsMaximum(Float_t q,Int_t max,const Float_t *bins){
+inline Bool_t AliTPCclustererMI::IsMaximum(Float_t q,Int_t max,const Float_t *bins) const {
//is this a local maximum ?
if (bins[-max] >= q) return kFALSE;
if (bins[-1 ] >= q) return kFALSE;