From fb41bb4974fde9c1f90a28690adc8a69efc616ab Mon Sep 17 00:00:00 2001
From: auras <antonio.uras@cern.ch>
Date: Wed, 8 Oct 2014 01:04:14 +0200
Subject: [PATCH] Updated MFT geometry and cluster information

---
 MFT/AliMFT.cxx              |   6 +++---
 MFT/AliMFTCluster.cxx       |   7 +++++--
 MFT/AliMFTCluster.h         |   5 ++++-
 MFT/AliMFTClusterFinder.cxx |   7 +++++--
 MFT/AliMFTGeometry.root     | Bin 7399 -> 7399 bytes
 MFT/SetMFTGeometry.C        |   2 +-
 6 files changed, 18 insertions(+), 9 deletions(-)

diff --git a/MFT/AliMFT.cxx b/MFT/AliMFT.cxx
index 9445b0d1eab..a88417ec83c 100644
--- a/MFT/AliMFT.cxx
+++ b/MFT/AliMFT.cxx
@@ -63,7 +63,7 @@ AliMFT::AliMFT():
   fChargeDispersion(25.e-4),
   fSingleStepForChargeDispersion(0),
   fNStepForChargeDispersion(4),
-  fDensitySupportOverSi(0.05),
+  fDensitySupportOverSi(0.036),
   fFileNameForUnderyingEvent(0), 
   fFileNameForPileUpEvents(0),
   fNPileUpEvents(0), 
@@ -93,7 +93,7 @@ AliMFT::AliMFT(const Char_t *name, const Char_t *title):
   fChargeDispersion(25.e-4),
   fSingleStepForChargeDispersion(0),
   fNStepForChargeDispersion(4),
-  fDensitySupportOverSi(0.05),
+  fDensitySupportOverSi(0.036),
   fFileNameForUnderyingEvent(0), 
   fFileNameForPileUpEvents(0),
   fNPileUpEvents(0), 
@@ -126,7 +126,7 @@ AliMFT::AliMFT(const Char_t *name, const Char_t *title, Char_t *nameGeomFile):
   fChargeDispersion(25.e-4),
   fSingleStepForChargeDispersion(0),
   fNStepForChargeDispersion(4),
-  fDensitySupportOverSi(0.05),
+  fDensitySupportOverSi(0.036),
   fFileNameForUnderyingEvent(0), 
   fFileNameForPileUpEvents(0),
   fNPileUpEvents(0), 
diff --git a/MFT/AliMFTCluster.cxx b/MFT/AliMFTCluster.cxx
index a764cd35945..3f18ebf73f9 100644
--- a/MFT/AliMFTCluster.cxx
+++ b/MFT/AliMFTCluster.cxx
@@ -48,7 +48,8 @@ AliMFTCluster::AliMFTCluster():
   fTrackChi2(0),
   fLocalChi2(0),
   fDigitsInCluster(0),
-  fIsClusterEditable(kTRUE)
+  fIsClusterEditable(kTRUE),
+  fIsClusterFront(kTRUE)
 {
 
   // default constructor
@@ -77,7 +78,8 @@ AliMFTCluster::AliMFTCluster(const AliMFTCluster& cluster):
   fTrackChi2(cluster.fTrackChi2),
   fLocalChi2(cluster.fLocalChi2),
   fDigitsInCluster(NULL),
-  fIsClusterEditable(cluster.fIsClusterEditable)
+  fIsClusterEditable(cluster.fIsClusterEditable),
+  fIsClusterFront(cluster.fIsClusterFront)
 {
 
   // copy constructor
@@ -122,6 +124,7 @@ AliMFTCluster& AliMFTCluster::operator=(const AliMFTCluster& cluster) {
   fTrackChi2         = cluster.fTrackChi2;
   fLocalChi2         = cluster.fLocalChi2;
   fIsClusterEditable = cluster.fIsClusterEditable;
+  fIsClusterFront    = cluster.fIsClusterFront;
 
   for (Int_t iTrack=0; iTrack<fNMaxMCTracks; iTrack++) fMCLabel[iTrack] = (cluster.fMCLabel)[iTrack];
   fDigitsInCluster      = new TClonesArray(*(cluster.fDigitsInCluster));
diff --git a/MFT/AliMFTCluster.h b/MFT/AliMFTCluster.h
index 7015349ad59..f01d3488f77 100644
--- a/MFT/AliMFTCluster.h
+++ b/MFT/AliMFTCluster.h
@@ -86,6 +86,9 @@ public:
 
   Double_t GetDistanceFromPixel(AliMFTDigit *pixel);
 
+  void SetClusterFront(Bool_t clusterFront) { if(fIsClusterEditable) fIsClusterFront = clusterFront; }
+  Bool_t IsClusterFront() { return fIsClusterFront; }
+
   AliMUONRawCluster* CreateMUONCluster();
   
 private:
@@ -108,7 +111,7 @@ private:
   
   TClonesArray *fDigitsInCluster;   //! (Temporary) Array of the digits composing the cluster
 
-  Bool_t fIsClusterEditable;
+  Bool_t fIsClusterEditable, fIsClusterFront;
 
   ClassDef(AliMFTCluster, 1)
 
diff --git a/MFT/AliMFTClusterFinder.cxx b/MFT/AliMFTClusterFinder.cxx
index f09aa007bbc..c56361ecdf0 100644
--- a/MFT/AliMFTClusterFinder.cxx
+++ b/MFT/AliMFTClusterFinder.cxx
@@ -230,12 +230,15 @@ void AliMFTClusterFinder::DigitsToClusters(const TObjArray *pDigitList) {
       for (Int_t iCluster=0; iCluster<clustersPerDetElem[iDetElem]->GetEntries(); iCluster++) {
 	newCluster = (AliMFTCluster*) (clustersPerDetElem[iDetElem]->At(iCluster));
 	newCluster -> TerminateCluster();
+	newCluster -> SetClusterEditable(kTRUE);
+	if (TMath::Abs(newCluster->GetZ())<TMath::Abs(fSegmentation->GetPlane(iPlane)->GetZCenter())) newCluster->SetClusterFront(kTRUE);
+	else                                                                                          newCluster->SetClusterFront(kFALSE);
 	if (fApplyMisalignment) {
-	  newCluster -> SetClusterEditable(kTRUE);
 	  newCluster -> SetX(newCluster->GetX()+misalignmentX);
 	  newCluster -> SetY(newCluster->GetY()+misalignmentY);
-	  newCluster -> SetClusterEditable(kFALSE);
 	}
+	newCluster -> SetClusterEditable(kFALSE);
+
 	new ((*fClustersPerPlane[iPlane])[fClustersPerPlane[iPlane]->GetEntries()]) AliMFTCluster(*newCluster);
       }
     }
diff --git a/MFT/AliMFTGeometry.root b/MFT/AliMFTGeometry.root
index 8b0007dc9469fb6abe42dc5463084761eaba16f0..9b93593725b06b3d4d4d3d6ebd2f6a243e9a58ff 100644
GIT binary patch
delta 271
zcmaEE`P_1XzJPSZ)K`9af={N4Gcq#N?|VPdQH_^@p^brsfq@~=Zvtc4#3FlTAO&Iq
zg=&G=3y4#I7zkvv85kiZOuVAV3{w3g6+}%=U`ztjN1?P9QwEqn6+%y9&INM|m;)3y
z8ynlBVW7!{%x{3C2*}#W51I48s)Jba!SoT9DlqNHS_!82vetp<&30^yS-^}ZTtR|h
l8A};6Fm1<XCI&Rl9_%eHuwP6W7=i9zG1*Z@Ve&m0P5={DOpyQp

delta 271
zcmaEE`P_1XzCg_B`>QPV1fNV7XJllk-}io^qZ%&*LmLAN0|P^#UufaViADCxKnlbJ
z3e^I!7Z9faF%ZaTGcZC-n0Q5z8KnA0Du|k#z?cN4k3wlJrVKEDDukZIoD1d_Fb60$
zFfiDoVW7!{%x{3C2*}#W51I48s)Jba!SoT9DlqNHS_!82vetp<&30^yS-^}ZTtR|h
l8A};6Fm1<XCI&Rl9_%eHuwP6W7=i9zG1*Z@Ve&m0P5?z>QicEk

diff --git a/MFT/SetMFTGeometry.C b/MFT/SetMFTGeometry.C
index 119f600be10..d5cd48b4afc 100644
--- a/MFT/SetMFTGeometry.C
+++ b/MFT/SetMFTGeometry.C
@@ -19,7 +19,7 @@ void SetMFTGeometry() {
   const Float_t pixelSizeY[nPlanes]       = { 28.e-4,  28.e-4,  28.e-4,  28.e-4,  28.e-4 };   // expressed in cm
 
   const Float_t thicknessActive[nPlanes]  = {  50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4 };   // expressed in cm
-  const Float_t thicknessSupport[nPlanes] = {   1.   ,    1.   ,    1.   ,    1.   ,    1.    };   // expressed in cm
+  const Float_t thicknessSupport[nPlanes] = {   1.4  ,    1.4  ,    1.4  ,    1.4  ,    1.4   };   // expressed in cm
   const Float_t thicknessReadout[nPlanes] = {  50.e-4,   50.e-4,   50.e-4,   50.e-4,   50.e-4 };   // expressed in cm
 
   const Float_t equivalentSilicon[nPlanes]            = { 600.e-4, 600.e-4, 600.e-4, 600.e-4, 600.e-4 };    // expressed in cm
-- 
2.43.0