void SetMinUpDownCutSigma(Float_t th) { fMinUpDownCutSigma=th;} // minimal amplitude up-down - TRF
//
Int_t GetLastSeedRowSec() const { return fLastSeedRowSec;}
+ Int_t GetSeedGapPrim() const { return fSeedGapPrim;}
+ Int_t GetSeedGapSec() const { return fSeedGapSec;}
void SetDoKinks(Bool_t on) { fBKinkFinder=on; }
Bool_t GetDoKinks() const { return fBKinkFinder;}
Double_t GetKinkAngleCutChi2(Int_t index) const {return fKinkAngleCutChi2[index];}
void SetKinkAngleCutChi2(Int_t index,Double_t value) {fKinkAngleCutChi2[index]=value;}
+ void SetSeedGapPrim(Int_t seedGapPrim) { fSeedGapPrim = seedGapPrim;}
+ void SetSeedGapSec(Int_t seedGapSec) { fSeedGapSec = seedGapSec;}
Float_t GetMaxC() const { return fMaxC;}
Bool_t GetSpecialSeeding() const { return fBSpecialSeeding;}
//
Bool_t fBKinkFinder; // do kink finder reconstruction
Double_t fKinkAngleCutChi2[2]; // angular cut for kinks
Int_t fLastSeedRowSec; // Most Inner Row to make seeding for secondaries
+ Int_t fSeedGapPrim; // seeding gap for primary tracks
+ Int_t fSeedGapSec; // seeding gap for secondary tracks
+
//
// Correction switches
//
// to be switched off for pass 0 reconstruction
// Use static function, other option will be to use
// additional specific storage ?
- ClassDef(AliTPCRecoParam, 10)
+ ClassDef(AliTPCRecoParam, 11)
};
TObjArray * seeds = new TObjArray;
TObjArray * arr=0;
+ Int_t fLastSeedRowSec=AliTPCReconstructor::GetRecoParam()->GetLastSeedRowSec();
+ Int_t gapPrim = AliTPCReconstructor::GetRecoParam()->GetSeedGapPrim();
+ Int_t gapSec = AliTPCReconstructor::GetRecoParam()->GetSeedGapSec();
Int_t gap =20;
Float_t cuts[4];
//
//find primaries
cuts[0]=0.0066;
- for (Int_t delta = 0; delta<18; delta+=6){
+ for (Int_t delta = 0; delta<18; delta+=gapPrim){
//
cuts[0]=0.0070;
cuts[1] = 1.5;
//find primaries
cuts[0]=0.0077;
- for (Int_t delta = 20; delta<120; delta+=10){
+ for (Int_t delta = 20; delta<120; delta+=gapPrim){
//
// seed high pt tracks
cuts[0]=0.0060;
SumTracks(seeds,arr);
SignClusters(seeds,fnumber,fdensity);
//
+ arr = Tracking(4,nup-5,nup-5-gap,cuts,-1);
+ SumTracks(seeds,arr);
+ SignClusters(seeds,fnumber,fdensity);
+ //
+ arr = Tracking(4,nup-7,nup-7-gap,cuts,-1);
+ SumTracks(seeds,arr);
+ SignClusters(seeds,fnumber,fdensity);
+ //
+ //
+ arr = Tracking(4,nup-9,nup-9-gap,cuts,-1);
+ SumTracks(seeds,arr);
+ SignClusters(seeds,fnumber,fdensity);
+ //
- for (Int_t delta = 3; delta<30; delta+=5){
+ for (Int_t delta = 9; delta<30; delta+=gapSec){
//
cuts[0] = 0.3;
cuts[1] = 1.5;
fdensity = 2.;
cuts[0]=0.0080;
- Int_t fLastSeedRowSec=AliTPCReconstructor::GetRecoParam()->GetLastSeedRowSec();
// find secondaries
- for (Int_t delta = 30; delta<fLastSeedRowSec; delta+=10){
+ for (Int_t delta = 30; delta<fLastSeedRowSec; delta+=gapSec){
//
cuts[0] = 0.3;
cuts[1] = 3.5;