27-oct-2005 NvE Memberfunction GetX to access individual vector components also intro...
authornick <nick@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 28 Oct 2005 08:15:05 +0000 (08:15 +0000)
committernick <nick@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 28 Oct 2005 08:15:05 +0000 (08:15 +0000)
                for Ali4Vector.
                Also memberfunction SortJets introduced in AliVertex.
                Misleading comment modified in AliDevice.cxx.

RALICE/Ali4Vector.cxx
RALICE/Ali4Vector.h
RALICE/AliDevice.cxx
RALICE/AliJet.h
RALICE/AliVertex.cxx
RALICE/AliVertex.h
RALICE/history.txt

index ade85d3..8d68a96 100644 (file)
@@ -970,6 +970,41 @@ Double_t Ali4Vector::GetGamma()
  return gamma;
 }
 ///////////////////////////////////////////////////////////////////////////
+Double_t Ali4Vector::GetX(Int_t i,TString f,TString u)
+{
+// Provide i-th vector component according to reference frame f.
+//
+// The string argument "u" allows to choose between different angular units
+// in case e.g. a spherical frame is selected.
+// u = "rad" : angles provided in radians
+//     "deg" : angles provided in degrees
+//
+// The default is u="rad".
+//
+// The vector components are addressed via the generic x0,x1,x2,x3 notation.
+// So, i=0 denotes the scalar component and i=1 denotes the first 3-vector component.
+// The error on the selected component can be obtained via the
+// usual GetResultError() facility.
+ fDresult=0;
+
+ if (i<0 || i>3) return 0;
+
+ Double_t x=0;
+
+ if (i==0)
+ {
+  x=GetScalar();
+ }
+ else
+ {
+  x=fV.GetX(i,f,u);
+  fDresult=fV.GetResultError();
+ }
+
+ return x;
+}
+///////////////////////////////////////////////////////////////////////////
 Double_t Ali4Vector::GetOpeningAngle(Ali4Vector& q,TString u)
 {
 // Provide the opening angle between 3-vector parts with 4-vector q.
index 6956e09..191af7f 100644 (file)
@@ -54,6 +54,7 @@ class Ali4Vector
   Ali3Vector GetBetaVector() const;                 // Provide the beta 3-vector
   Double_t GetBeta();                               // Provide the norm of the beta 3-vector, i.e. v/c
   Double_t GetGamma();                              // Provide the Lorentz gamma factor
+  Double_t GetX(Int_t i,TString f,TString u="rad"); // Provide i-th vector component in frame f in units u
   virtual Double_t GetOpeningAngle(Ali4Vector& q,TString u="rad"); // Opening angle between 3-vector parts in units u
 
  protected:
@@ -67,6 +68,6 @@ class Ali4Vector
   Double_t GetScaTrans(); // Provide "transverse value" of scalar part w.r.t. z-axis
   Double_t GetScaLong();  // Provide "longitudinal value" of scalar part w.r.t. z-axis
 
- ClassDef(Ali4Vector,8) // Handling of Lorentz 4-vectors in various reference frames.
+ ClassDef(Ali4Vector,9) // Handling of Lorentz 4-vectors in various reference frames.
 };
 #endif
index 163b431..9f3715c 100644 (file)
@@ -569,7 +569,7 @@ TObjArray* AliDevice::SortHits(Int_t idx,Int_t mode,TObjArray* hits,Int_t mcal)
     fOrdered->AddAt(fOrdered->At(k-1),k);
    }
    fOrdered->AddAt(s,j); // put hit at empty position
-   break; // go for next matrix module
+   break; // go for next hit
   }
  }
  return fOrdered;
@@ -654,7 +654,7 @@ TObjArray* AliDevice::SortHits(TString name,Int_t mode,TObjArray* hits,Int_t mca
     fOrdered->AddAt(fOrdered->At(k-1),k);
    }
    fOrdered->AddAt(s,j); // put hit at empty position
-   break; // go for next matrix module
+   break; // go for next hit
   }
  }
  return fOrdered;
index 7567538..4c4d80d 100644 (file)
@@ -61,7 +61,7 @@ class AliJet : public TNamed,public Ali4Vector
   TObjArray* fTracks;                    // Array to hold the pointers to the tracks of the jet
   Int_t fTrackCopy;                      // Flag to denote creation of private copies in fTracks
   Int_t fUserId;                         // The user defined identifier
-  TObjArray* fSelected;                  //! Temp. array to hold user selected objects
+  TObjArray* fSelected;                  //! Temp. array to hold user selected or ordered objects
  
  ClassDef(AliJet,12) // Creation and investigation of a jet of particle tracks.
 };
index 2f5df32..fc0e9a1 100644 (file)
@@ -1014,6 +1014,148 @@ void AliVertex::Draw(Int_t secs,Int_t cons,Int_t jets)
  }
 }
 ///////////////////////////////////////////////////////////////////////////
+TObjArray* AliVertex::SortJets(Int_t mode,TObjArray* jets)
+{
+// Order the references to an array of jets by looping over the input array "jets"
+// and checking the value of a certain observable.
+// The ordered array is returned as a TObjArray.
+// In case jets=0 (default), the registered jets of the current vertex are used. 
+// Note that the original jet array is not modified.
+// Via the "mode" argument the user can specify the observable to be checked upon
+// and specify whether sorting should be performed in decreasing order (mode<0)
+// or in increasing order (mode>0).
+//
+// The convention for the observable selection is the following :
+// mode : 1 ==> Number of tracks in the jet
+//        2 ==> Jet energy
+//        3 ==> Jet momentum
+//        4 ==> Invariant mass of the jet
+//        5 ==> Transverse momentum of the jet
+//        6 ==> Longitudinal momentum of the jet
+//        7 ==> Transverse energy of the jet
+//        8 ==> Longitudinal energy of the jet
+//        9 ==> Transverse mass of the jet
+//       10 ==> Jet rapidity
+//       11 ==> Pseudo-rapidity of the jet
+//
+// The default is mode=-1.
+//
+// Note : This sorting routine uses a common area in memory, which is used
+//        by various other sorting facilities as well.
+//        This means that the resulting sorted TObjArray may be overwritten
+//        when another sorting is invoked.
+//        To retain the sorted list of pointers, the user is advised to copy
+//        the pointers contained in the returned TObjArray into a private
+//        TObjArray instance.
+
+ if (fSelected)
+ {
+  delete fSelected;
+  fSelected=0;
+ }
+
+ if (!jets) jets=fJets;
+ if (abs(mode)>11 || !jets) return fSelected;
+
+ Int_t njets=jets->GetEntries();
+ if (!njets)
+ {
+  return fSelected;
+ }
+ else
+ {
+  fSelected=new TObjArray(njets);
+ }
+
+ Double_t val1,val2; // Values of the observable to be tested upon
+ Int_t nord=0;
+ for (Int_t i=0; i<njets; i++) // Loop over all jets of the array
+ {
+  AliJet* jx=(AliJet*)jets->At(i);
+
+  if (!jx) continue;
+  if (nord == 0) // store the first jet at the first ordered position
+  {
+   nord++;
+   fSelected->AddAt(jx,nord-1);
+   continue;
+  }
+  for (Int_t j=0; j<=nord; j++) // put jet in the right ordered position
+  {
+   if (j == nord) // jet has smallest (mode<0) or largest (mode>0) observable value seen so far
+   {
+    nord++;
+    fSelected->AddAt(jx,j); // add jet at the end
+    break; // go for next jet
+   }
+   
+   switch (abs(mode))
+   {
+    case 1:
+     val1=jx->GetNtracks();
+     val2=((AliJet*)fSelected->At(j))->GetNtracks();
+     break;
+    case 2:
+     val1=jx->GetEnergy();
+     val2=((AliJet*)fSelected->At(j))->GetEnergy();
+     break;
+    case 3:
+     val1=jx->GetMomentum();
+     val2=((AliJet*)fSelected->At(j))->GetMomentum();
+     break;
+    case 4:
+     val1=jx->GetInvmass();
+     val2=((AliJet*)fSelected->At(j))->GetInvmass();
+     break;
+    case 5:
+     val1=jx->GetPt();
+     val2=((AliJet*)fSelected->At(j))->GetPt();
+     break;
+    case 6:
+     val1=jx->GetPl();
+     val2=((AliJet*)fSelected->At(j))->GetPl();
+     break;
+    case 7:
+     val1=jx->GetEt();
+     val2=((AliJet*)fSelected->At(j))->GetEt();
+     break;
+    case 8:
+     val1=jx->GetEl();
+     val2=((AliJet*)fSelected->At(j))->GetEl();
+     break;
+    case 9:
+     val1=jx->GetMt();
+     val2=((AliJet*)fSelected->At(j))->GetMt();
+     break;
+    case 10:
+     val1=jx->GetRapidity();
+     val2=((AliJet*)fSelected->At(j))->GetRapidity();
+     break;
+    case 11:
+     val1=jx->GetPseudoRapidity();
+     val2=((AliJet*)fSelected->At(j))->GetPseudoRapidity();
+     break;
+   }
+
+   if (mode<0 && val1 < val2) continue;
+   if (mode>0 && val1 > val2) continue;
+   nord++;
+   for (Int_t k=nord-1; k>j; k--) // create empty position
+   {
+    fSelected->AddAt(fSelected->At(k-1),k);
+   }
+   fSelected->AddAt(jx,j); // put jet at empty position
+   break; // go for next jet
+  }
+ }
+ return fSelected;
+}
+///////////////////////////////////////////////////////////////////////////
 TObject* AliVertex::Clone(const char* name) const
 {
 // Make a deep copy of the current object and provide the pointer to the copy.
index f4fcb22..4fa3345 100644 (file)
@@ -48,6 +48,7 @@ class AliVertex : public AliJet,public AliPosition
   Int_t IsJetTrack(AliTrack* t) const;     // Indicate if track is resulting from jet addition
   virtual void Draw(Option_t*) { Draw(1,1,0); } // Override TObject::Draw for default event display
   virtual void Draw(Int_t secs,Int_t cons=1,Int_t jets=0); // Draw the vertex in an event display
+  TObjArray* SortJets(Int_t mode=-1,TObjArray* jets=0); // Sort jets by the number of tracks
 
  protected:
   void Init();          // Initialisation of pointers etc... 
@@ -66,6 +67,6 @@ class AliVertex : public AliJet,public AliPosition
  private:
   void Dumps(AliVertex* v,Int_t n,TString f,TString u); // Recursively print all sec. vertices
  
- ClassDef(AliVertex,10) // Creation and investigation of an AliVertex.
+ ClassDef(AliVertex,11) // Creation and investigation of an AliVertex.
 };
 #endif
index 7c22f69..aed03f4 100644 (file)
                 Memberfunction GetOpeningAngle introduced in Ali3Vector and Ali4Vector.
 25-oct-2005 NvE Support introduced for user selectable angular units (rad or deg) in Ali3Vector,
                 Ali4Vector and all derived classes.
+27-oct-2005 NvE Memberfunction GetX to access individual vector components also introduced
+                for Ali4Vector.
+                Also memberfunction SortJets introduced in AliVertex.
+                Misleading comment modified in AliDevice.cxx.