#include "AliTracker.h"
#include "AliESDtrack.h"
#include "AliESDVertex.h"
+#include "TTreeStream.h"
ClassImp(AliTPCtrack)
}
//_____________________________________________________________________________
-AliTPCtrack::AliTPCtrack(const AliESDtrack& t) :
+AliTPCtrack::AliTPCtrack(const AliESDtrack& t, TTreeSRedirector *pcstream) :
AliKalmanTrack(),
fdEdx(t.GetTPCsignal()),
fSdEdx(1e10),
//-----------------------------------------------------------------
// Conversion AliESDtrack -> AliTPCtrack.
//-----------------------------------------------------------------
+ const Double_t kmaxC[4]={10,10,0.1,0.1}; // cuts on the rms /fP0,fP1,fP2,fP3
SetNumberOfClusters(t.GetTPCclusters(fIndex));
SetLabel(t.GetLabel());
SetMass(t.GetMass());
for (Int_t i=0; i<12;i++) fKinkPoint[i]=0.;
for (Int_t i=0; i<3;i++) fKinkIndexes[i]=0;
for (Int_t i=0; i<3;i++) fV0Indexes[i]=0;
+ //
+ // choose parameters to start
+ //
+ Int_t reject=0;
+ AliExternalTrackParam param(t);
+ const AliExternalTrackParam *tpcout=(t.GetFriendTrack())? ((AliESDfriendTrack*)(t.GetFriendTrack()))->GetTPCOut():0;
+ const AliExternalTrackParam *tpcin = t.GetInnerParam();
+ const AliExternalTrackParam *tpc=(tpcout)?tpcout:tpcin;
+ if (!tpc) tpc=¶m;
+ Bool_t isOK=kTRUE;
+ if (param.GetCovariance()[0]>kmaxC[0]*kmaxC[0]) isOK=kFALSE;
+ if (param.GetCovariance()[2]>kmaxC[1]*kmaxC[1]) isOK=kFALSE;
+ if (param.GetCovariance()[5]>kmaxC[2]*kmaxC[2]) isOK=kFALSE;
+ if (param.GetCovariance()[9]>kmaxC[3]*kmaxC[3]) isOK=kFALSE;
+ if (!isOK){
+ param=*tpc;
+ isOK=kTRUE;
+ reject=1;
+ }
+ param.Rotate(tpc->GetAlpha());
+ isOK=AliTracker::PropagateTrackToBxByBz(¶m,tpc->GetX(),t.GetMass(),2.,kFALSE);
+ if (param.GetCovariance()[0]>kmaxC[0]*kmaxC[0]) isOK=kFALSE;
+ if (param.GetCovariance()[2]>kmaxC[1]*kmaxC[1]) isOK=kFALSE;
+ if (param.GetCovariance()[5]>kmaxC[2]*kmaxC[2]) isOK=kFALSE;
+ if (param.GetCovariance()[9]>kmaxC[3]*kmaxC[3]) isOK=kFALSE;
+ if (!isOK){
+ param=*tpc;
+ isOK=kTRUE;
+ reject=2;
+ }
+ if (reject>0){
+ param.ResetCovariance(4.); // reset covariance if start from backup param
+ }
+ //
+ //
+ if (pcstream){
+ AliExternalTrackParam dummy;
+ AliExternalTrackParam *ptpc=(AliExternalTrackParam *)tpc;
+ if (!ptpc) ptpc=&dummy;
+ AliESDtrack *esd= (AliESDtrack *)&t;
+ (*pcstream)<<"trackP"<<
+ "reject="<<reject<< // flag - rejection of current esd track parameters
+ "esd.="<<esd<< // original esd track
+ "tr.="<<¶m<< // starting track parameters
+ "out.="<<ptpc<< // backup tpc parameters
+ "\n";
+ }
- Set(t.GetX(),t.GetAlpha(),t.GetParameter(),t.GetCovariance());
+ Set(param.GetX(),param.GetAlpha(),param.GetParameter(),param.GetCovariance());
if ((t.GetStatus()&AliESDtrack::kTIME) == 0) return;
StartTimeIntegral();
#include "AliExternalTrackParam.h"
class AliESDtrack;
class AliESDVertex;
-
+class TTreeSRedirector;
//_____________________________________________________________________________
class AliTPCtrack : public AliKalmanTrack {
public:
AliTPCtrack();
AliTPCtrack(Double_t x, Double_t alpha, const Double_t p[5],
const Double_t cov[15], Int_t index);
- AliTPCtrack(const AliESDtrack& t);
+ AliTPCtrack(const AliESDtrack& t, TTreeSRedirector *pcstream=0);
AliTPCtrack(const AliTPCtrack& t);
AliTPCtrack& operator=(const AliTPCtrack& o);
virtual ~AliTPCtrack() {}