From aad7688e8f3eb86abd23511298ec9958d8f6d875 Mon Sep 17 00:00:00 2001 From: marian Date: Tue, 16 Oct 2007 12:58:04 +0000 Subject: [PATCH] removing files - files splitted to 3 separate files --- PWG1/AliRecInfo.cxx | 719 -------------------------------------------- PWG1/AliRecInfo.h | 163 ---------- 2 files changed, 882 deletions(-) delete mode 100644 PWG1/AliRecInfo.cxx delete mode 100644 PWG1/AliRecInfo.h diff --git a/PWG1/AliRecInfo.cxx b/PWG1/AliRecInfo.cxx deleted file mode 100644 index f58eb1842ae..00000000000 --- a/PWG1/AliRecInfo.cxx +++ /dev/null @@ -1,719 +0,0 @@ -/************************************************************************** - * 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 -#include -//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 "AliRecInfo.h" - - - -ClassImp(AliESDRecInfo) -ClassImp(AliESDRecV0Info) -ClassImp(AliESDRecKinkInfo) - - - - -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]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=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 &&undeffmaxdens){ - 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]0;i--){ - if (density[i]<0) continue; - if (density[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]) 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]); - } - -} - - -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]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 (delta1GetInnerExternalParameters(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); - } -} - -//// -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]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 (delta1GetInnerExternalParameters(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); - */ - } - - // - // - } - -} - diff --git a/PWG1/AliRecInfo.h b/PWG1/AliRecInfo.h deleted file mode 100644 index aaa655fe635..00000000000 --- a/PWG1/AliRecInfo.h +++ /dev/null @@ -1,163 +0,0 @@ -#ifndef ALIRECINFO_H -#define ALIRECINFO_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 "AliESDEvent.h" -#include "AliESDtrack.h" -#include "AliV0.h" -#include "AliESDkink.h" -#include "AliESDfriendTrack.h" -#include "AliITStrackMI.h" -#include "AliTRDtrack.h" -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 -}; - - - -///////////////////////////////////////////////////////// -///////////////////////////////////////////////////////// -///////////////////////////////////////////////////////// - - -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 -}; - - - -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 -- 2.43.0