Merging fixes from v4-04-Release
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 24 May 2006 10:28:15 +0000 (10:28 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 24 May 2006 10:28:15 +0000 (10:28 +0000)
RICH/AliRICHRecon.cxx
TPC/AliTPCseed.cxx

index e8652088bb85001e2d8ca7bd115eb55805dafa06..b7f906dc1136989f53646a056427bb0efd9ee39d 100644 (file)
 #define NPointsOfRing 201
 
 //__________________________________________________________________________________________________
-AliRICHRecon::AliRICHRecon()
-             :TTask       ("RichRec","RichPat")
+AliRICHRecon::AliRICHRecon():
+  TTask       ("RichRec","RichPat"),
+  fClusters(0x0),
+  fParam(AliRICHParam::Instance()),
+  fPhotonsNumber(0),
+  fPhotonIndex(0), // should we use -1?
+  fFittedHoughPhotons(0),
+  fMinNumPhots(3),
+  fTrackTheta(0),
+  fTrackPhi(0),
+  fMinDist(999),
+  fTrackBeta(0),
+  fXtoentr(999),
+  fYtoentr(999),
+  fThetaPhotonInTRS(0),
+  fPhiPhotonInTRS(0),
+  fThetaPhotonInDRS(0),
+  fPhiPhotonInDRS(0),
+  fThetaAtQuartz(0),
+  fXEmiss(-999),
+  fYEmiss(-999),
+  fXInner(-999),
+  fYInner(-999),
+  fXOuter(-999),
+  fYOuter(-999),
+  fInnerRadius(-999),
+  fOuterRadius(-999),
+  fEphot(0),
+  fLengthEmissionPoint(0),
+  fPhotonLimitX(999),
+  fPhotonLimitY(999),
+  fDistanceFromCluster(999),
+  fCerenkovAnglePad(999),
+  fThetaPhotonCerenkov(0),
+  fShiftX(0),
+  fShiftY(0),
+  fXcoord(999),
+  fYcoord(999),
+  fIntersectionX(999),
+  fIntersectionY(999),
+  fMassHypotesis(0.139567),
+  fThetaOfRing(0),
+  fAreaOfRing(0),
+  fPortionOfRing(0),
+  fHoughArea(0),
+  fHoughRMS(999),
+  fCandidatePhotonX(0x0),
+  fCandidatePhotonY(0x0),
+  fFittedTrackTheta(0),
+  fFittedTrackPhi(0),
+  fFittedThetaCerenkov(0),
+  fThetaMin(0.0),
+  fThetaMax(0.75),
+  fIsWEIGHT(kFALSE),
+  fIsBACKGROUND(kFALSE),
+  fDTheta(0.001),
+  fWindowWidth(0.045),
+  fNumEtaPhotons(0),
+  fnPhotBKG(0),
+  fThetaCerenkov(0),
+  fWeightThetaCerenkov(0),
+  fThetaPeakPos(0)
 {
 // main ctor
-  fThetaMin = 0.0; fThetaMax = 0.75; 
-  fDTheta       = 0.001;   fWindowWidth  = 0.045;
-  fMinNumPhots = 3;
-  fParam=AliRICHParam::Instance(); //get the pointer to AliRICHParam
+  for (Int_t i=0; i<3000; i++) {
+    fPhotonFlag[i] = 0;
+    fPhiPoint[i] = -999;
+    fPhotonEta[i] = -999;
+    fPhotonWeight[i] = 0.0;
+    fEtaFlag[i] = 0;
+    fEtaPhotons[i] = 0;
+    fWeightPhotons[i] = 0;
+  }
 }
 //__________________________________________________________________________________________________
 Double_t AliRICHRecon::ThetaCerenkov(AliRICHHelix *pHelix,TClonesArray *pClusters,Int_t &iMipId)
index 6c13f61021f39c169b44175a066e3694daaf961a..46dbc370effa27100474f82927174f5eb0eebb6e 100644 (file)
@@ -56,7 +56,7 @@ AliTPCseed::AliTPCseed():
   fCircular(0)
 {
   //
-  for (Int_t i=0;i<200;i++) SetClusterIndex2(i,-3);
+  for (Int_t i=0;i<160;i++) SetClusterIndex2(i,-3);
   for (Int_t i=0;i<160;i++) fClusterPointer[i]=0;
   for (Int_t i=0;i<3;i++)   fKinkIndexes[i]=0;
   for (Int_t i=0;i<AliPID::kSPECIES;i++)   fTPCr[i]=0.2;
@@ -108,7 +108,16 @@ AliTPCseed::AliTPCseed(const AliTPCseed &s, Bool_t clusterOwner):
     fTrackPoints[i] = s.fTrackPoints[i];
   }
   for (Int_t i=0;i<160;i++) fIndex[i] = s.fIndex[i];
+  for (Int_t i=0;i<AliPID::kSPECIES;i++)   fTPCr[i]=s.fTPCr[i];
+  for (Int_t i=0;i<4;i++) {
+    fDEDX[i] = s.fDEDX[i];
+    fSDEDX[i] = s.fSDEDX[i];
+    fNCDEDX[i] = s.fNCDEDX[i];
+  }
+  for (Int_t i=0;i<12;i++) fOverlapLabels[i] = s.fOverlapLabels[i];
 }
+
+
 AliTPCseed::AliTPCseed(const AliTPCtrack &t):
   AliTPCtrack(t),
   fEsd(0x0),
@@ -139,7 +148,6 @@ AliTPCseed::AliTPCseed(const AliTPCtrack &t):
   // Constructor from AliTPCtrack
   //
   fFirstPoint =0;
-  for (Int_t i=0;i<3;i++)   fKinkIndexes[i]=t.GetKinkIndex(i);
   for (Int_t i=0;i<5;i++)   fTPCr[i]=0.2;
   for (Int_t i=0;i<160;i++) {
     fClusterPointer[i] = 0;
@@ -191,9 +199,8 @@ AliTPCseed::AliTPCseed(UInt_t index,  const Double_t xx[5],
   // Constructor
   //
   fFirstPoint =0;
-  for (Int_t i=0;i<200;i++) SetClusterIndex2(i,-3);
+  for (Int_t i=0;i<160;i++) SetClusterIndex2(i,-3);
   for (Int_t i=0;i<160;i++) fClusterPointer[i]=0;
-  for (Int_t i=0;i<3;i++)   fKinkIndexes[i]=0;
   for (Int_t i=0;i<5;i++)   fTPCr[i]=0.2;
   for (Int_t i=0;i<4;i++) {
     fDEDX[i] = 0.;