]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TPC/AliTPCtrackerParam.h
Use debug stream only if requested
[u/mrichter/AliRoot.git] / TPC / AliTPCtrackerParam.h
index a9650778f1c65e6269b5f9536f6caab5329eecba..340579a6f6d0b9cdd4cf407fefbef8b7e2d57cc5 100644 (file)
@@ -6,50 +6,41 @@
 /* $Id$ */
 
 
-//-----------------------------------------------------------------------------
-//                    TPC Tracking Parameterization Class
-//
-//   Origin: Andrea Dainese, Padova - e-mail: andrea.dainese@pd.infn.it
-//
-//
-//
-//-----------------------------------------------------------------------------
+/////////////////////////////////////////////////////////////////////////
+//                                                                        
+// This class builds AliTPCtrack objects from generated tracks to feed    
+// ITS tracking (V2). The AliTPCtrack is built from its first hit in      
+// the TPC. The track is assigned a Kalman-like covariance matrix         
+// depending on its pT and pseudorapidity and track parameters are        
+// smeared according to this covariance matrix.                           
+// Output file contains sorted tracks, ready for matching with ITS.        
+//                                                                        
+// See implementation file for more details.  
+//                                  
+//                                                                        
+//  Origin: Andrea Dainese, Padova - e-mail: andrea.dainese@pd.infn.it     
+//  adapted to ESD output by: Marcello Lunardon, Padova
+/////////////////////////////////////////////////////////////////////////
+
 
 //----- Root headers ---------
+class TTree;
 #include <TMatrixD.h>
 //---- AliRoot headers -------
-//#include "alles.h"
-#include "AliRun.h"
-//#include "AliGausCorr.h"
-//#include "AliMagF.h"
+#include "AliConfig.h"
 #include "AliTPCkineGrid.h"
 #include "AliTPCtrack.h"
-//#include "AliTrackReference.h"
+#include "AliESDEvent.h"
 //----------------------------
 
 class AliTPC;
 
 class AliTPCtrackerParam:
   public TObject
-{
-   
-  /////////////////////////////////////////////////////////////////////////
-  //                                                                        
-  // This class builds AliTPCtrack objects from generated tracks to feed    
-  // ITS tracking (V2). The AliTPCtrack is built from its first hit in      
-  // the TPC. The track is assigned a Kalman-like covariance matrix         
-  // depending on its pT and pseudorapidity and track parameters are        
-  // smeared according to this covariance matrix.                           
-  // Output file contains sorted tracks, ready for matching with ITS.        
-  //                                                                        
-  // See implementation file for more details.  
-  //                                  
-  //                                                                        
-  //  Origin: Andrea Dainese, Padova - e-mail: andrea.dainese@pd.infn.it     
-  //                                                                        
-  /////////////////////////////////////////////////////////////////////////
+{   
  public:
-  AliTPCtrackerParam(Int_t coll=0, Double_t Bz=0.4, Int_t n=1,
+  
+  AliTPCtrackerParam(Int_t coll=0, Double_t Bz=0.4,
                     const char* evfoldname = AliConfig::GetDefaultEventFolderName());
   virtual ~AliTPCtrackerParam();
 
@@ -57,14 +48,16 @@ class AliTPCtrackerParam:
   //
   AliTPCtrackerParam(const AliTPCtrackerParam& p);
   //
-  Int_t BuildTPCtracks(const TFile *inp, TFile *out);
+  Int_t Init();
+  Int_t BuildTPCtracks(AliESDEvent* event);
   // these functions are used to create a DB of cov. matrices,
   // including regularization, efficiencies and dE/dx
   void  AllGeantTracks() { fSelAndSmear=kFALSE; return; }
   void  AnalyzedEdx(const Char_t *outName,Int_t pdg);
   void  AnalyzePulls(const Char_t *outName);
   void  AnalyzeResolutions(Int_t pdg);
-  void  CompareTPCtracks(const Char_t *galiceName="galice.root",
+  void  CompareTPCtracks(Int_t nEvents=1,
+                        const Char_t *galiceName="galice.root",
                         const Char_t *trkGeaName="AliTPCtracksGeant.root",
                         const Char_t *trkKalName="AliTPCtracksSorted.root",
                         const Char_t *covmatName="CovMatrix.root",
@@ -86,18 +79,31 @@ class AliTPCtrackerParam:
     void AssignMass(Double_t mass) {SetMass(mass); return;}
     
   private:
-
   };
   //********* end of internal class ***********
 
   //********* Internal class definition *******
   class AliTPCseedGeant : public TObject {
   public:
+    AliTPCseedGeant():TObject(),
+      fXg(0.),
+      fYg(0.),
+      fZg(0.),
+      fPx(0.),
+      fPy(0.),
+      fPz(0.),
+      fAlpha(0.),
+      fLabel(0),
+      fSector(0){}
+
     AliTPCseedGeant(Double_t x=0.,Double_t y=0.,Double_t z=0.,
                    Double_t px=0.,Double_t py=0.,Double_t pz=0.,
                    Int_t lab=0);
     Int_t    GetLabel() const { return fLabel; }
     Double_t GetAlpha() const { return fAlpha; }      
+    Double_t GetXG() const { return fXg; }
+    Double_t GetYG() const { return fYg; }
+    Double_t GetZG() const { return fZg; }
     Double_t GetXL() const 
       { return fXg*TMath::Cos(fAlpha)+fYg*TMath::Sin(fAlpha); }
     Double_t GetYL() const 
@@ -129,9 +135,10 @@ class AliTPCtrackerParam:
   //******* end of internal class ****************
   
  private:
+  AliTPCtrackerParam & operator=(const AliTPCtrackerParam & );
+
   TString fEvFolderName;//! name of data folder
 
-  Int_t           fNevents;     // number of events in the file to be processed
   Double_t        fBz;          // value of the z component of L3 field (Tesla)
   Int_t           fColl;        // collision code (0: PbPb6000; 1: pp)
   Bool_t          fSelAndSmear; // if kFALSE returns GEANT tracks 
@@ -142,6 +149,11 @@ class AliTPCtrackerParam:
 
   TTree          *fCovTree;  // tree with regularized cov matrices 
                              // for the current track
+  TTree           *fCovTreePi[30];
+  TTree           *fCovTreeKa[30];
+  TTree           *fCovTreePr[30];
+  TTree           *fCovTreeEl[30];
+  TTree           *fCovTreeMu[30];
 
   AliTPCkineGrid *fDBgrid;   // grid for the cov matrix look-up table  
   AliTPCkineGrid  fDBgridPi; //               "                  for pions  
@@ -187,8 +199,6 @@ class AliTPCtrackerParam:
 
 
   void     BuildTrack(AliTPCseedGeant *s,Int_t ch);
-  Int_t    CheckLabel(AliTPCseedGeant *s,Int_t nPart,
-                     Double_t *ptkine,Double_t *pzkine) const;
   void     CookdEdx(Double_t pt,Double_t eta);
   void     CookTrack(Double_t pt,Double_t eta);
   TMatrixD GetSmearingMatrix(Double_t *cc, Double_t pt,Double_t eta) const;
@@ -202,6 +212,7 @@ class AliTPCtrackerParam:
   Int_t    ReadEffs(const Char_t *inName);
   Int_t    ReadPulls(const Char_t *inName);
   Int_t    ReadRegParams(const Char_t *inName,Int_t pdg);
+  Int_t    ReadCovTrees(const Char_t* inName); 
   Bool_t   SelectedTrack(Double_t pt, Double_t eta) const;
   void     SetParticle(Int_t pdg);  
   void     SmearTrack(Double_t *xx,Double_t *xxsm,TMatrixD cov) const;