--- /dev/null
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * *
+ * Author: The ALICE Off-line Project. *
+ * Contributors are mentioned in the code where appropriate. *
+ * *
+ * Permission to use, copy, modify and distribute this software and its *
+ * documentation strictly for non-commercial purposes is hereby granted *
+ * without fee, provided that the above copyright notice appears in all *
+ * copies and that both the copyright notice and this permission notice *
+ * appear in the supporting documentation. The authors make no claims *
+ * about the suitability of this software for any purpose. It is *
+ * provided "as is" without express or implied warranty. *
+ **************************************************************************/
+
+
+///////////////////////////////////////////////////////////////////////////////
+// //
+// Time Projection Chamber //
+// Comparison macro for reconstructed tracks - ESDs V0s //
+// responsible:
+// marian.ivanov@cern.ch //
+//
+//
+
+
+
+
+
+#include <stdio.h>
+#include <string.h>
+//ROOT includes
+#include "Rtypes.h"
+//
+//ALIROOT includes
+//
+#include "AliESDtrack.h"
+#include "AliTPCParam.h"
+#include "AliTrackReference.h"
+#include "AliTPCParamSR.h"
+#include "AliESD.h"
+#include "AliESDfriend.h"
+#include "AliESDtrack.h"
+#include "AliTPCseed.h"
+#include "AliITStrackMI.h"
+#include "AliTRDtrack.h"
+#include "AliHelix.h"
+#include "AliESDVertex.h"
+#include "AliExternalTrackParam.h"
+#include "AliESDkink.h"
+#include "AliESDv0.h"
+#include "AliV0.h"
+//
+#include "AliTreeDraw.h"
+#include "AliGenInfo.h"
+#include "AliESDRecInfo.h"
+
+
+
+ClassImp(AliESDRecInfo)
+
+
+
+
+AliTPCParam * GetTPCParam(){
+ AliTPCParamSR * par = new AliTPCParamSR;
+ par->Update();
+ return par;
+}
+
+
+
+
+AliESDRecInfo::AliESDRecInfo():
+ fITSOn(0), // ITS refitted inward
+ fTRDOn(0), // ITS refitted inward
+ fDeltaP(0), //delta of momenta
+ fSign(0), // sign
+ fReconstructed(0), //flag if track was reconstructed
+ fFake(0), // fake track
+ fMultiple(0), // number of reconstructions
+ fTPCOn(0), // TPC refitted inward
+ fBestTOFmatch(0), //best matching between times
+ fESDtrack(0), // esd track
+ fTrackF(0), // friend track
+ fTPCtrack(0), // tpc track
+ fITStrack(0), // its track
+ fTRDtrack(0) // trd track
+{
+ //
+ // default constructor
+ //
+}
+
+
+AliESDRecInfo::AliESDRecInfo(const AliESDRecInfo& recinfo):
+ TObject()
+{
+ //
+ //
+ //
+ memcpy(this,&recinfo, sizeof(recinfo));
+ fESDtrack=0; fTrackF=0; fTPCtrack=0;fITStrack=0;fTRDtrack=0;
+ SetESDtrack(recinfo.GetESDtrack());
+}
+
+
+AliESDRecInfo::~AliESDRecInfo()
+
+{
+ //
+ // destructor
+ //
+ if (fESDtrack) { delete fESDtrack; fESDtrack=0;}
+ if (fTrackF) { delete fTrackF; fTrackF=0;}
+ if (fTPCtrack) { delete fTPCtrack; fTPCtrack=0;}
+ if (fITStrack) { delete fITStrack; fITStrack=0;}
+ if (fTRDtrack) { delete fTRDtrack; fTRDtrack=0;}
+
+}
+
+
+
+void AliESDRecInfo::Reset()
+{
+ //
+ // reset info
+ //
+ fMultiple =0;
+ fFake =0;
+ fReconstructed=0;
+ if (fESDtrack) { delete fESDtrack; fESDtrack=0;}
+ if (fTrackF) { delete fTrackF; fTrackF=0;}
+ if (fTPCtrack) { delete fTPCtrack; fTPCtrack=0;}
+ if (fITStrack) { delete fITStrack; fITStrack=0;}
+ if (fTRDtrack) { delete fTRDtrack; fTRDtrack=0;}
+}
+
+void AliESDRecInfo::SetESDtrack(const AliESDtrack *track){
+ //
+ //
+ //
+ if (fESDtrack) delete fESDtrack;
+ fESDtrack = (AliESDtrack*)track->Clone();
+ if (0 &&track->GetFriendTrack()){
+ if (fTrackF) delete fTrackF;
+ fTrackF = (AliESDfriendTrack*)track->GetFriendTrack()->Clone();
+ if (fTrackF->GetCalibObject(0)){
+ if (fTPCtrack) delete fTPCtrack;
+ fTPCtrack = (AliTPCseed*)fTrackF->GetCalibObject(0)->Clone();
+ }
+ }
+
+}
+
+void AliESDRecInfo::UpdatePoints(AliESDtrack*track)
+{
+ //
+ //
+ Int_t iclusters[200];
+ Float_t density[160];
+ for (Int_t i=0;i<160;i++) density[i]=-1.;
+ fTPCPoints[0]= 160;
+ fTPCPoints[1] = -1;
+ //
+ if (fTPCPoints[0]<fTPCPoints[1]) return;
+ // Int_t nclusters=track->GetTPCclusters(iclusters);
+
+ Int_t ngood=0;
+ Int_t undeff=0;
+ Int_t nall =0;
+ Int_t range=20;
+ for (Int_t i=0;i<160;i++){
+ Int_t last = i-range;
+ if (nall<range) nall++;
+ if (last>=0){
+ if (iclusters[last]>0&& (iclusters[last]&0x8000)==0) ngood--;
+ if (iclusters[last]==-1) undeff--;
+ }
+ if (iclusters[i]>0&& (iclusters[i]&0x8000)==0) ngood++;
+ if (iclusters[i]==-1) undeff++;
+ if (nall==range &&undeff<range/2) density[i-range/2] = Float_t(ngood)/Float_t(nall-undeff);
+ }
+ Float_t maxdens=0;
+ Int_t indexmax =0;
+ for (Int_t i=0;i<160;i++){
+ if (density[i]<0) continue;
+ if (density[i]>maxdens){
+ maxdens=density[i];
+ indexmax=i;
+ }
+ }
+ //
+ //max dens point
+ fTPCPoints[3] = maxdens;
+ fTPCPoints[1] = indexmax;
+ //
+ // last point
+ for (Int_t i=indexmax;i<160;i++){
+ if (density[i]<0) continue;
+ if (density[i]<maxdens/2.) {
+ break;
+ }
+ fTPCPoints[2]=i;
+ }
+ //
+ // first point
+ for (Int_t i=indexmax;i>0;i--){
+ if (density[i]<0) continue;
+ if (density[i]<maxdens/2.) {
+ break;
+ }
+ fTPCPoints[0]=i;
+ }
+ //
+ // Density at the last 30 padrows
+ //
+ //
+ nall = 0;
+ ngood = 0;
+ for (Int_t i=159;i>0;i--){
+ if (iclusters[i]==-1) continue; //dead zone
+ nall++;
+ if (iclusters[i]>0) ngood++;
+ if (nall>20) break;
+ }
+ fTPCPoints[4] = Float_t(ngood)/Float_t(nall);
+ //
+ if ((track->GetStatus()&AliESDtrack::kITSrefit)>0) fTPCPoints[0]=-1;
+
+
+}
+
+//
+//
+void AliESDRecInfo::Update(AliMCInfo* info,AliTPCParam * /*par*/, Bool_t reconstructed)
+{
+ //
+ //
+ //calculates derived variables
+ //
+ //
+ UpdatePoints(fESDtrack);
+ fBestTOFmatch=1000;
+ AliTrackReference * ref = &(info->fTrackRef);
+ fTPCinR0[0] = info->fTrackRef.X();
+ fTPCinR0[1] = info->fTrackRef.Y();
+ fTPCinR0[2] = info->fTrackRef.Z();
+ fTPCinR0[3] = TMath::Sqrt(fTPCinR0[0]*fTPCinR0[0]+fTPCinR0[1]*fTPCinR0[1]);
+ fTPCinR0[4] = TMath::ATan2(fTPCinR0[1],fTPCinR0[0]);
+ //
+ fTPCinP0[0] = ref->Px();
+ fTPCinP0[1] = ref->Py();
+ fTPCinP0[2] = ref->Pz();
+ fTPCinP0[3] = ref->Pt();
+ fTPCinP0[4] = ref->P();
+ fDeltaP = (ref->P()-info->fParticle.P())/info->fParticle.P();
+ //
+ //
+ if (fTPCinP0[3]>0.0000001){
+ //
+ fTPCAngle0[0] = TMath::ATan2(fTPCinP0[1],fTPCinP0[0]);
+ fTPCAngle0[1] = TMath::ATan(fTPCinP0[2]/fTPCinP0[3]);
+ }
+ //
+ //
+ fITSinP0[0]=info->fParticle.Px();
+ fITSinP0[1]=info->fParticle.Py();
+ fITSinP0[2]=info->fParticle.Pz();
+ fITSinP0[3]=info->fParticle.Pt();
+ //
+ fITSinR0[0]=info->fParticle.Vx();
+ fITSinR0[1]=info->fParticle.Vy();
+ fITSinR0[2]=info->fParticle.Vz();
+ fITSinR0[3] = TMath::Sqrt(fITSinR0[0]*fITSinR0[0]+fITSinR0[1]*fITSinR0[1]);
+ fITSinR0[4] = TMath::ATan2(fITSinR0[1],fITSinR0[0]);
+ //
+ //
+ if (fITSinP0[3]>0.0000001){
+ fITSAngle0[0] = TMath::ATan2(fITSinP0[1],fITSinP0[0]);
+ fITSAngle0[1] = TMath::ATan(fITSinP0[2]/fITSinP0[3]);
+ }
+ //
+ for (Int_t i=0;i<4;i++) fStatus[i] =0;
+ fReconstructed = kFALSE;
+ fTPCOn = kFALSE;
+ fITSOn = kFALSE;
+ fTRDOn = kFALSE;
+ if (reconstructed==kFALSE) return;
+
+ fLabels[0] = info->fLabel;
+ fLabels[1] = info->fPrimPart;
+ fReconstructed = kTRUE;
+ fTPCOn = ((fESDtrack->GetStatus()&AliESDtrack::kTPCrefit)>0) ? kTRUE : kFALSE;
+ fITSOn = ((fESDtrack->GetStatus()&AliESDtrack::kITSrefit)>0) ? kTRUE : kFALSE;
+ fTRDOn = ((fESDtrack->GetStatus()&AliESDtrack::kTRDrefit)>0) ? kTRUE : kFALSE;
+ //
+ //
+ if ((fESDtrack->GetStatus()&AliESDtrack::kTPCrefit)>0){
+ fStatus[1] =3;
+ }
+ else{
+ if ((fESDtrack->GetStatus()&AliESDtrack::kTPCout)>0){
+ fStatus[1] =2;
+ }
+ else{
+ if ((fESDtrack->GetStatus()&AliESDtrack::kTPCin)>0)
+ fStatus[1]=1;
+ }
+ }
+ //
+ if ((fESDtrack->GetStatus()&AliESDtrack::kITSout)>0){
+ fStatus[0] =2;
+ }
+ else{
+ if ((fESDtrack->GetStatus()&AliESDtrack::kITSrefit)>0){
+ fStatus[0] =1;
+ }
+ else{
+ fStatus[0]=0;
+ }
+ }
+
+ //
+ //
+ if ((fESDtrack->GetStatus()&AliESDtrack::kTRDrefit)>0){
+ fStatus[2] =2;
+ }
+ else{
+ if ((fESDtrack->GetStatus()&AliESDtrack::kTRDout)>0){
+ fStatus[2] =1;
+ }
+ }
+ if ((fESDtrack->GetStatus()&AliESDtrack::kTRDStop)>0){
+ fStatus[2] =10;
+ }
+
+ //
+ //TOF
+ //
+ if (((fESDtrack->GetStatus()&AliESDtrack::kTOFout)>0)){
+ //
+ // best tof match
+ Double_t times[5];
+ fESDtrack->GetIntegratedTimes(times);
+ for (Int_t i=0;i<5;i++){
+ if ( TMath::Abs(fESDtrack->GetTOFsignal()-times[i]) <TMath::Abs(fBestTOFmatch) ){
+ fBestTOFmatch = fESDtrack->GetTOFsignal()-times[i];
+ }
+ }
+ Int_t toflabel[3];
+ fESDtrack->GetTOFLabel(toflabel);
+ Bool_t toffake=kTRUE;
+ Bool_t tofdaughter=kFALSE;
+ for (Int_t i=0;i<3;i++){
+ if (toflabel[i]<0) continue;
+ if (toflabel[i]== TMath::Abs(fESDtrack->GetLabel())) toffake=kFALSE;
+ if (toflabel[i]==info->fParticle.GetDaughter(0) || (toflabel[i]==info->fParticle.GetDaughter(1))) tofdaughter=kTRUE; // decay product of original particle
+ fStatus[3]=1;
+ }
+ if (toffake) fStatus[3] =3; //total fake
+ if (tofdaughter) fStatus[3]=2; //fake because of decay
+ }else{
+ fStatus[3]=0;
+ }
+
+
+ if (fStatus[1]>0 &&info->fNTPCRef>0&&TMath::Abs(fTPCinP0[3])>0.0001){
+ //TPC
+ fESDtrack->GetInnerXYZ(fTPCinR1);
+ fTPCinR1[3] = TMath::Sqrt(fTPCinR1[0]*fTPCinR1[0]+fTPCinR1[1]*fTPCinR1[1]);
+ fTPCinR1[4] = TMath::ATan2(fTPCinR1[1],fTPCinR1[0]);
+ fESDtrack->GetInnerPxPyPz(fTPCinP1);
+ fTPCinP1[3] = TMath::Sqrt(fTPCinP1[0]*fTPCinP1[0]+fTPCinP1[1]*fTPCinP1[1]);
+ fTPCinP1[4] = TMath::Sqrt(fTPCinP1[3]*fTPCinP1[3]+fTPCinP1[2]*fTPCinP1[2]);
+ //
+ //
+ if (fTPCinP1[3]>0.000000000000001){
+ fTPCAngle1[0] = TMath::ATan2(fTPCinP1[1],fTPCinP1[0]);
+ fTPCAngle1[1] = TMath::ATan(fTPCinP1[2]/fTPCinP1[3]);
+ }
+ Double_t cov[15], param[5],x, alpha;
+ fESDtrack->GetInnerExternalCovariance(cov);
+ fESDtrack->GetInnerExternalParameters(alpha, x,param);
+ if (x<50) return ;
+ //
+ fTPCDelta[0] = (fTPCinR0[4]-fTPCinR1[4])*fTPCinR1[3]; //delta rfi
+ fTPCPools[0] = fTPCDelta[0]/TMath::Sqrt(cov[0]);
+ fTPCDelta[1] = (fTPCinR0[2]-fTPCinR1[2]); //delta z
+ fTPCPools[1] = fTPCDelta[1]/TMath::Sqrt(cov[2]);
+ fTPCDelta[2] = (fTPCAngle0[0]-fTPCAngle1[0]);
+ fTPCPools[2] = fTPCDelta[2]/TMath::Sqrt(cov[5]);
+ fTPCDelta[3] = (TMath::Tan(fTPCAngle0[1])-TMath::Tan(fTPCAngle1[1]));
+ fTPCPools[3] = fTPCDelta[3]/TMath::Sqrt(cov[9]);
+ fTPCDelta[4] = (fTPCinP0[3]-fTPCinP1[3]);
+ Double_t sign = (param[4]>0)? 1.:-1;
+ fSign =sign;
+ fTPCPools[4] = sign*(1./fTPCinP0[3]-1./fTPCinP1[3])/TMath::Sqrt(TMath::Abs(cov[14]));
+ }
+ if (fITSOn){
+ // ITS
+ Double_t param[5],x;
+ fESDtrack->GetExternalParameters(x,param);
+ // fESDtrack->GetConstrainedExternalParameters(x,param);
+ Double_t cov[15];
+ fESDtrack->GetExternalCovariance(cov);
+ //fESDtrack->GetConstrainedExternalCovariance(cov);
+ if (TMath::Abs(param[4])<0.0000000001) return;
+
+ fESDtrack->GetXYZ(fITSinR1);
+ fESDtrack->GetPxPyPz(fITSinP1);
+ fITSinP1[3] = TMath::Sqrt(fITSinP1[0]*fITSinP1[0]+fITSinP1[1]*fITSinP1[1]);
+ //
+ fITSinR1[3] = TMath::Sqrt(fITSinR1[0]*fITSinR1[0]+fITSinR1[1]*fITSinR1[1]);
+ fITSinR1[4] = TMath::ATan2(fITSinR1[1],fITSinR1[0]);
+ //
+ //
+ if (fITSinP1[3]>0.0000001){
+ fITSAngle1[0] = TMath::ATan2(fITSinP1[1],fITSinP1[0]);
+ fITSAngle1[1] = TMath::ATan(fITSinP1[2]/fITSinP1[3]);
+ }
+ //
+ //
+ fITSDelta[0] = (fITSinR0[4]-fITSinR1[4])*fITSinR1[3]; //delta rfi
+ fITSPools[0] = fITSDelta[0]/TMath::Sqrt(cov[0]);
+ fITSDelta[1] = (fITSinR0[2]-fITSinR1[2]); //delta z
+ fITSPools[1] = fITSDelta[1]/TMath::Sqrt(cov[2]);
+ fITSDelta[2] = (fITSAngle0[0]-fITSAngle1[0]);
+ fITSPools[2] = fITSDelta[2]/TMath::Sqrt(cov[5]);
+ fITSDelta[3] = (TMath::Tan(fITSAngle0[1])-TMath::Tan(fITSAngle1[1]));
+ fITSPools[3] = fITSDelta[3]/TMath::Sqrt(cov[9]);
+ fITSDelta[4] = (fITSinP0[3]-fITSinP1[3]);
+ Double_t sign = (param[4]>0) ? 1:-1;
+ fSign = sign;
+ fITSPools[4] = sign*(1./fITSinP0[3]-1./fITSinP1[3])/TMath::Sqrt(cov[14]);
+ }
+
+}
+
+
--- /dev/null
+#ifndef ALIESDRECINFO_H
+#define ALIESDRECINFO_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice */
+
+
+
+//////////////////////////////////////////////////////////////////////////////
+// Class AliRecInfo //
+// collect together MC info and Rec info for comparison purposes
+// - effieciency studies and so on// //
+// marian.ivanov@cern.ch //
+//////////////////////////////////////////////////////////////////////////////
+
+
+#include "TObject.h"
+#include "AliESDkink.h"
+class AliESDEvent;
+class AliESDtrack;
+class AliV0;
+class AliESDkink;
+class AliESDfriendTrack;
+class AliITStrackMI;
+class AliTRDtrack;
+
+class AliTPCseed;
+
+/////////////////////////////////////////////////////////////////////////
+class AliESDRecInfo: public TObject {
+ friend class AliRecInfoMaker;
+ friend class AliESDRecV0Info;
+ friend class AliESDRecKinkInfo;
+
+public:
+ AliESDRecInfo();
+ AliESDRecInfo(const AliESDRecInfo& recinfo);
+ ~AliESDRecInfo();
+ void UpdatePoints(AliESDtrack* track);
+ void Update(AliMCInfo* info,AliTPCParam * par, Bool_t reconstructed);
+ void Reset();
+ //
+ void SetESDtrack(const AliESDtrack *track);
+ AliESDtrack *GetESDtrack() const { return fESDtrack;}
+ AliESDfriendTrack *GetTrackF() const { return fTrackF;}
+ AliTPCseed *GetTPCtrack() const { return fTPCtrack;}
+ AliITStrackMI *GetITStrack() const { return fITStrack;}
+ AliTRDtrack *GetTRDtrack() const { return fTRDtrack;}
+ Int_t GetStatus(Int_t i) { return fStatus[i];}
+protected:
+ //
+ Float_t fTPCPoints[10]; //start , biggest end points,max density .. density at the last 30 pad-rows
+ Double_t fTPCinR0[5]; //generated position of the track at inner tpc - radius [3] and fi [4]
+ Double_t fTPCinR1[5]; //reconstructed postion of the track - radius [3] and fi [
+ Double_t fTPCinP0[5]; //generated position of the track at inner tpc
+ Double_t fTPCinP1[5]; //reconstructed postion of the track
+ Double_t fTPCAngle0[2]; // generated angle
+ Double_t fTPCAngle1[2]; //refconstructed angle
+ Double_t fTPCDelta[5]; // deltas
+ Double_t fTPCPools[5]; // pools
+ Double_t fITSinR0[5]; //generated position of the track at inner tpc
+ Double_t fITSinR1[5]; //reconstructed postion of the track
+ Double_t fITSinP0[5]; //generated position of the track at inner tpc
+ Double_t fITSinP1[5]; //reconstructed postion of the track
+ Double_t fITSAngle0[2]; // generated angle
+ Double_t fITSAngle1[2]; //refconstructed angle
+ Double_t fITSDelta[5]; // deltas
+ Double_t fITSPools[5]; // pools
+ Float_t fTRLocalCoord[3]; //local coordinates of the track ref.
+ Int_t fStatus[4]; // status -0 not found - 1 -only in - 2 -in-out -3 -in -out-refit
+ Int_t fLabels[2]; // labels
+
+ Bool_t fITSOn; // ITS refitted inward
+ Bool_t fTRDOn; // ITS refitted inward
+ Float_t fDeltaP; //delta of momenta
+ Double_t fSign; // sign
+ Int_t fReconstructed; //flag if track was reconstructed
+ Int_t fFake; // fake track
+ Int_t fMultiple; // number of reconstructions
+ Bool_t fTPCOn; // TPC refitted inward
+ Float_t fBestTOFmatch; //best matching between times
+
+private:
+ AliESDtrack *fESDtrack; // esd track
+ AliESDfriendTrack *fTrackF; // friend track
+ AliTPCseed *fTPCtrack; // tpc track
+ AliITStrackMI *fITStrack; // its track
+ AliTRDtrack *fTRDtrack; // trd track
+
+ ClassDef(AliESDRecInfo,2) // container for
+};
+
+
+#endif
--- /dev/null
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * *
+ * Author: The ALICE Off-line Project. *
+ * Contributors are mentioned in the code where appropriate. *
+ * *
+ * Permission to use, copy, modify and distribute this software and its *
+ * documentation strictly for non-commercial purposes is hereby granted *
+ * without fee, provided that the above copyright notice appears in all *
+ * copies and that both the copyright notice and this permission notice *
+ * appear in the supporting documentation. The authors make no claims *
+ * about the suitability of this software for any purpose. It is *
+ * provided "as is" without express or implied warranty. *
+ **************************************************************************/
+
+
+///////////////////////////////////////////////////////////////////////////////
+// //
+// Time Projection Chamber //
+// Comparison macro for reconstructed tracks - ESDs V0s //
+// responsible:
+// marian.ivanov@cern.ch //
+//
+//
+
+
+
+
+
+#include <stdio.h>
+#include <string.h>
+//ROOT includes
+#include "Rtypes.h"
+//
+//ALIROOT includes
+//
+#include "AliESDtrack.h"
+#include "AliTPCParam.h"
+#include "AliTrackReference.h"
+#include "AliTPCParamSR.h"
+#include "AliESD.h"
+#include "AliESDfriend.h"
+#include "AliESDtrack.h"
+#include "AliTPCseed.h"
+#include "AliITStrackMI.h"
+#include "AliTRDtrack.h"
+#include "AliHelix.h"
+#include "AliESDVertex.h"
+#include "AliExternalTrackParam.h"
+#include "AliESDkink.h"
+#include "AliESDv0.h"
+#include "AliV0.h"
+//
+#include "AliTreeDraw.h"
+#include "AliGenInfo.h"
+#include "AliESDRecKinkInfo.h"
+
+
+
+ClassImp(AliESDRecKinkInfo)
+
+
+
+
+
+////
+void AliESDRecKinkInfo::Update()
+{
+
+ if ( (fT1.fTPCOn)&& (fT2.fTPCOn)){
+ //
+ // IF BOTH RECONSTRUCTED
+ Float_t distance1,distance2;
+ Double_t xx[3],pp[3];
+ //
+ Double_t xd[3],pd[3],signd;
+ Double_t xm[3],pm[3],signm;
+ for (Int_t i=0;i<3;i++){
+ xd[i] = fT2.fTPCinR1[i];
+ pd[i] = fT2.fTPCinP1[i];
+ xm[i] = fT1.fTPCinR1[i];
+ pm[i] = fT1.fTPCinP1[i];
+ }
+ signd = fT2.fSign<0 ? -1:1;
+ signm = fT1.fSign<0 ? -1:1;
+
+ AliHelix dhelix1(xd,pd,signd);
+ dhelix1.GetMomentum(0,pp,0);
+ dhelix1.Evaluate(0,xx);
+ //
+ // Double_t x2[3],p2[3];
+ //
+ AliHelix mhelix(xm,pm,signm);
+ //
+ //find intersection linear
+ //
+ Double_t phase[2][2],radius[2];
+ Int_t points = dhelix1.GetRPHIintersections(mhelix, phase, radius,200);
+ Double_t delta1=10000,delta2=10000;
+
+ if (points==1){
+ fMinR = TMath::Sqrt(radius[0]);
+ }
+ if (points==2){
+ fMinR =TMath::Min(TMath::Sqrt(radius[0]),TMath::Sqrt(radius[1]));
+ }
+
+ if (points>0){
+ dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
+ dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
+ dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
+ }
+ if (points==2){
+ dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
+ dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
+ dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
+ }
+ if (points==1){
+ fMinR = TMath::Sqrt(radius[0]);
+ fDistMinR = delta1;
+ }
+ if (points==2){
+ if (radius[0]<radius[1]){
+ fMinR = TMath::Sqrt(radius[0]);
+ fDistMinR = delta1;
+ }
+ else{
+ fMinR = TMath::Sqrt(radius[1]);
+ fDistMinR = delta2;
+ }
+ }
+ //
+ //
+ distance1 = TMath::Min(delta1,delta2);
+ //
+ //find intersection parabolic
+ //
+ points = dhelix1.GetRPHIintersections(mhelix, phase, radius);
+ delta1=10000,delta2=10000;
+
+ if (points>0){
+ dhelix1.ParabolicDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
+ }
+ if (points==2){
+ dhelix1.ParabolicDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
+ }
+
+ distance2 = TMath::Min(delta1,delta2);
+ if (delta1<delta2){
+ //get V0 info
+ dhelix1.Evaluate(phase[0][0],fXr);
+ dhelix1.GetMomentum(phase[0][0],fPdr);
+ mhelix.GetMomentum(phase[0][1],fPm);
+ dhelix1.GetAngle(phase[0][0],mhelix,phase[0][1],fAngle);
+ fRr = TMath::Sqrt(radius[0]);
+ }
+ else{
+ dhelix1.Evaluate(phase[1][0],fXr);
+ dhelix1.GetMomentum(phase[1][0], fPdr);
+ mhelix.GetMomentum(phase[1][1], fPm);
+ dhelix1.GetAngle(phase[1][0],mhelix,phase[1][1],fAngle);
+ fRr = TMath::Sqrt(radius[1]);
+ }
+ fDist1 = TMath::Sqrt(distance1);
+ fDist2 = TMath::Sqrt(distance2);
+
+ if (fDist2<10.5){
+ Double_t x,alpha,param[5],cov[15];
+ //
+ fT1.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
+ fT1.GetESDtrack()->GetInnerExternalCovariance(cov);
+ AliExternalTrackParam paramm(x,alpha,param,cov);
+ //
+ fT2.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
+ fT2.GetESDtrack()->GetInnerExternalCovariance(cov);
+ AliExternalTrackParam paramd(x,alpha,param,cov);
+ /*
+ AliESDkink kink;
+ kink.Update(¶mm,¶md);
+ // kink.Dump();
+ Double_t diff = kink.fRr-fRr;
+ Double_t diff2 = kink.fDist2-fDist2;
+ printf("Diff\t%f\t%f\n",diff,diff2);
+ */
+ }
+
+ //
+ //
+ }
+
+}
+
--- /dev/null
+#ifndef ALIESDRECKINKINFO_H
+#define ALIESDRECKINKINFO_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice */
+
+
+
+//////////////////////////////////////////////////////////////////////////////
+// Class AliRecInfo //
+// collect together MC info and Rec info for comparison purposes
+// - effieciency studies and so on// //
+// marian.ivanov@cern.ch //
+//////////////////////////////////////////////////////////////////////////////
+
+
+#include "TObject.h"
+#include "AliESDRecInfo.h"
+
+
+
+class AliESDRecKinkInfo: public TObject {
+friend class AliRecInfoMaker;
+public:
+ void Update();
+protected:
+ AliESDRecInfo fT1; //track1
+ AliESDRecInfo fT2; //track2
+ AliESDkink fKink; //kink
+ Double_t fDist1; //info about closest distance according closest MC - linear DCA
+ Double_t fDist2; //info about closest distance parabolic DCA
+ Double_t fInvMass; //reconstructed invariant mass -
+ //
+ Double_t fPdr[3]; //momentum at vertex daughter - according approx at DCA
+ Double_t fXr[3]; //rec. position according helix
+ //
+ Double_t fPm[3]; //momentum at the vertex mother
+ Double_t fAngle[3]; //three angles
+ Double_t fRr; // rec position of the vertex
+ Double_t fMinR; // minimum radius in rphi intersection
+ Double_t fDistMinR; // distance at minimal radius
+ Int_t fLab[2]; //MC label of the partecle
+ Float_t fPointAngleFi; //point angle fi
+ Float_t fPointAngleTh; //point angle theta
+ Float_t fPointAngle; //point angle full
+ Int_t fStatus; //status -tracks
+ Int_t fRecStatus; //kink -status- 0 - not found 1-good - fake
+ Int_t fMultiple; // how many times was kink reconstructed
+ Int_t fKinkMultiple; // how many times was kink reconstructed
+ ClassDef(AliESDRecKinkInfo,1) // container for
+};
+
+#endif
--- /dev/null
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * *
+ * Author: The ALICE Off-line Project. *
+ * Contributors are mentioned in the code where appropriate. *
+ * *
+ * Permission to use, copy, modify and distribute this software and its *
+ * documentation strictly for non-commercial purposes is hereby granted *
+ * without fee, provided that the above copyright notice appears in all *
+ * copies and that both the copyright notice and this permission notice *
+ * appear in the supporting documentation. The authors make no claims *
+ * about the suitability of this software for any purpose. It is *
+ * provided "as is" without express or implied warranty. *
+ **************************************************************************/
+
+
+///////////////////////////////////////////////////////////////////////////////
+// //
+// //
+// Comparison class for V0 information //
+// responsible:
+// marian.ivanov@cern.ch //
+//
+//
+
+
+
+
+
+#include <stdio.h>
+#include <string.h>
+//ROOT includes
+#include "Rtypes.h"
+//
+//ALIROOT includes
+//
+#include "AliESDtrack.h"
+#include "AliTPCParam.h"
+#include "AliTrackReference.h"
+#include "AliTPCParamSR.h"
+#include "AliESD.h"
+#include "AliESDfriend.h"
+#include "AliESDtrack.h"
+#include "AliTPCseed.h"
+#include "AliITStrackMI.h"
+#include "AliTRDtrack.h"
+#include "AliHelix.h"
+#include "AliESDVertex.h"
+#include "AliExternalTrackParam.h"
+#include "AliESDkink.h"
+#include "AliESDv0.h"
+#include "AliV0.h"
+//
+#include "AliTreeDraw.h"
+#include "AliGenInfo.h"
+
+#include "AliESDRecV0Info.h"
+
+
+
+ClassImp(AliESDRecV0Info)
+
+
+
+
+
+void AliESDRecV0Info::Update(Float_t vertex[3])
+{
+
+ if ( (fT1.fStatus[1]>0)&& (fT2.fStatus[1]>0)){
+ Float_t distance1,distance2;
+ Double_t xx[3],pp[3];
+ //
+ Double_t xd[3],pd[3],signd;
+ Double_t xm[3],pm[3],signm;
+ //
+ //
+ if (fT1.fITSOn&&fT2.fITSOn){
+ for (Int_t i=0;i<3;i++){
+ xd[i] = fT2.fITSinR1[i];
+ pd[i] = fT2.fITSinP1[i];
+ xm[i] = fT1.fITSinR1[i];
+ pm[i] = fT1.fITSinP1[i];
+ }
+ }
+ else{
+
+ for (Int_t i=0;i<3;i++){
+ xd[i] = fT2.fTPCinR1[i];
+ pd[i] = fT2.fTPCinP1[i];
+ xm[i] = fT1.fTPCinR1[i];
+ pm[i] = fT1.fTPCinP1[i];
+ }
+ }
+ //
+ //
+ signd = fT2.fSign<0 ? -1:1;
+ signm = fT1.fSign<0 ? -1:1;
+
+ AliHelix dhelix1(xd,pd,signd);
+ dhelix1.GetMomentum(0,pp,0);
+ dhelix1.Evaluate(0,xx);
+ //
+ // Double_t x2[3],p2[3];
+ //
+ AliHelix mhelix(xm,pm,signm);
+ //
+ //find intersection linear
+ //
+ Double_t phase[2][2],radius[2];
+ Int_t points = dhelix1.GetRPHIintersections(mhelix, phase, radius,200);
+ Double_t delta1=10000,delta2=10000;
+
+ if (points==1){
+ fRs[0] = TMath::Sqrt(radius[0]);
+ fRs[1] = TMath::Sqrt(radius[0]);
+ }
+ if (points==2){
+ fRs[0] =TMath::Min(TMath::Sqrt(radius[0]),TMath::Sqrt(radius[1]));
+ fRs[1] =TMath::Max(TMath::Sqrt(radius[0]),TMath::Sqrt(radius[1]));
+ }
+
+ if (points>0){
+ dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
+ dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
+ dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
+ }
+ if (points==2){
+ dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
+ dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
+ dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
+ }
+ if (points==1){
+ fRs[0] = TMath::Sqrt(radius[0]);
+ fRs[1] = TMath::Sqrt(radius[0]);
+ fDistMinR = delta1;
+ }
+ if (points==2){
+ if (radius[0]<radius[1]){
+ fRs[0] = TMath::Sqrt(radius[0]);
+ fRs[1] = TMath::Sqrt(radius[1]);
+ fDistMinR = delta1;
+ }
+ else{
+ fRs[0] = TMath::Sqrt(radius[1]);
+ fRs[1] = TMath::Sqrt(radius[0]);
+ fDistMinR = delta2;
+ }
+ }
+ //
+ //
+ distance1 = TMath::Min(delta1,delta2);
+ //
+ //find intersection parabolic
+ //
+ points = dhelix1.GetRPHIintersections(mhelix, phase, radius);
+ delta1=10000,delta2=10000;
+
+ if (points>0){
+ dhelix1.ParabolicDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
+ }
+ if (points==2){
+ dhelix1.ParabolicDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
+ }
+
+ distance2 = TMath::Min(delta1,delta2);
+ if (distance2>100) fDist2 =100;
+ return;
+ if (delta1<delta2){
+ //get V0 info
+ dhelix1.Evaluate(phase[0][0],fXr);
+ dhelix1.GetMomentum(phase[0][0],fPdr);
+ mhelix.GetMomentum(phase[0][1],fPm);
+ dhelix1.GetAngle(phase[0][0],mhelix,phase[0][1],fAngle);
+ fRr = TMath::Sqrt(radius[0]);
+ }
+ else{
+ dhelix1.Evaluate(phase[1][0],fXr);
+ dhelix1.GetMomentum(phase[1][0], fPdr);
+ mhelix.GetMomentum(phase[1][1], fPm);
+ dhelix1.GetAngle(phase[1][0],mhelix,phase[1][1],fAngle);
+ fRr = TMath::Sqrt(radius[1]);
+ }
+ fDist1 = TMath::Sqrt(distance1);
+ fDist2 = TMath::Sqrt(distance2);
+
+ if (fDist2<10.5){
+ Double_t x,alpha,param[5],cov[15];
+ //
+ fT1.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
+ fT1.GetESDtrack()->GetInnerExternalCovariance(cov);
+ AliExternalTrackParam paramm(x,alpha,param,cov);
+ //
+ fT2.GetESDtrack()->GetInnerExternalParameters(alpha,x,param);
+ fT2.GetESDtrack()->GetInnerExternalCovariance(cov);
+ AliExternalTrackParam paramd(x,alpha,param,cov);
+ }
+ //
+ //
+
+ Float_t v[3] = {fXr[0]-vertex[0],fXr[1]-vertex[1],fXr[2]-vertex[2]};
+ Float_t p[3] = {fPdr[0]+fPm[0], fPdr[1]+fPm[1],fPdr[2]+fPm[2]};
+
+ Float_t vnorm2 = v[0]*v[0]+v[1]*v[1];
+ Float_t vnorm3 = TMath::Sqrt(v[2]*v[2]+vnorm2);
+ vnorm2 = TMath::Sqrt(vnorm2);
+ Float_t pnorm2 = p[0]*p[0]+p[1]*p[1];
+ Float_t pnorm3 = TMath::Sqrt(p[2]*p[2]+pnorm2);
+ pnorm2 = TMath::Sqrt(pnorm2);
+
+ fPointAngleFi = (v[0]*p[0]+v[1]*p[1])/(vnorm2*pnorm2);
+ fPointAngleTh = (v[2]*p[2]+vnorm2*pnorm2)/(vnorm3*pnorm3);
+ fPointAngle = (v[0]*p[0]+v[1]*p[1]+v[2]*p[2])/(vnorm3*pnorm3);
+ }
+}
+
--- /dev/null
+#ifndef ALIESDRECV0INFO_H
+#define ALIESDRECV0INFO_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice */
+
+
+
+#include "TObject.h"
+#include "AliESDRecInfo.h"
+
+
+
+
+/////////////////////////////////////////////////////////
+/////////////////////////////////////////////////////////
+/////////////////////////////////////////////////////////
+
+
+class AliESDRecV0Info: public TObject {
+ friend class AliRecInfoMaker;
+public:
+ void Update(Float_t vertex[3]);
+protected:
+ AliESDRecInfo fT1; //track1
+ AliESDRecInfo fT2; //track2
+ Double_t fDist1; //info about closest distance according closest MC - linear DCA
+ Double_t fDist2; //info about closest distance parabolic DCA
+ Double_t fInvMass; //reconstructed invariant mass -
+ //
+ Double_t fPdr[3]; //momentum at vertex daughter - according approx at DCA
+ Double_t fXr[3]; //rec. position according helix
+ //
+ Double_t fRs[2]; // minimum radius in rphi intersection
+ Double_t fDistMinR; // distance at minimal radius
+ Double_t fPm[3]; //momentum at the vertex mother
+ Double_t fAngle[3]; //three angles
+ Double_t fRr; // rec position of the vertex
+ Int_t fLab[2]; //MC label of the partecle
+ Float_t fPointAngleFi; //point angle fi
+ Float_t fPointAngleTh; //point angle theta
+ Float_t fPointAngle; //point angle full
+ Int_t fV0Status; // status of the kink
+ AliV0* fV0tpc; // Vo information from reconsturction according TPC
+ AliV0* fV0its; // Vo information from reconsturction according ITS
+ AliV0* fV0rec; // V0 information form the reconstruction
+ Int_t fMultiple; // how man times V0 was recostructed
+ Int_t fV0Multiple; // how man times was V0 reconstucted
+ Int_t fRecStatus; // status form the reconstuction
+ ClassDef(AliESDRecV0Info,2) // container for
+};
+
+
+
+#endif