update in the documentation
authorabercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 23 Jun 2009 18:25:01 +0000 (18:25 +0000)
committerabercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 23 Jun 2009 18:25:01 +0000 (18:25 +0000)
TRD/AliTRDseedV1.cxx
TRD/AliTRDtrackV1.cxx
TRD/AliTRDtrackerV1.cxx

index 876be63..6fd816b 100644 (file)
@@ -558,10 +558,15 @@ Bool_t AliTRDseedV1::CookPID()
 // Output
 //   returns pointer to the probability array and 0x0 if missing DB access 
 //
-// Detailed description
+// Retrieve PID probabilities for e+-, mu+-, K+-, pi+- and p+- from the DB according to tracklet information:
+// - estimated momentum at tracklet reference point 
+// - dE/dx measurements
+// - tracklet length
+// - TRD layer
+// According to the steering settings specified in the reconstruction one of the following methods are used
+// - Neural Network [default] - option "nn"  
+// - 2D Likelihood - option "!nn"  
 
-  
-  // retrive calibration db
   AliTRDcalibDB *calibration = AliTRDcalibDB::Instance();
   if (!calibration) {
     AliError("No access to calibration data");
index ac1edff..8e988a1 100644 (file)
@@ -278,11 +278,20 @@ Bool_t AliTRDtrackV1::CookLabel(Float_t wrong)
 //_______________________________________________________________
 Bool_t AliTRDtrackV1::CookPID()
 {
-  //
-  // Cook the PID information
-  //
+//
+// Cook the PID information for the track by delegating the omonim function of the tracklets. 
+// Computes the number of tracklets used. The tracklet information are considered independent. 
+// For the moment no global track measurement of PID is performed as for example to estimate 
+// bremsstrahlung probability based on global chi2 of the track.
+//
+// The status bit AliESDtrack::kTRDpid is set during the call of AliTRDtrackV1::UpdateESDtrack().The PID performance of the 
+//TRD for tracks with 6 tacklets is displayed below.
+//Begin_Html
+//<img src="TRD/trackPID.gif">
+//End_Html
+//
   
-  // Reset the a priori probabilities
+  /*Reset the a priori probabilities*/
   Double_t pid = 1. / AliPID::kSPECIES;
   for(int ispec=0; ispec<AliPID::kSPECIES; ispec++) fPID[ispec] = pid; 
 
@@ -290,7 +299,7 @@ Bool_t AliTRDtrackV1::CookPID()
   // no tracklet found for PID calculations
   if(!fPIDquality) return kFALSE;
   
-  // slot for PID calculation @ track level
+  // slot for PID calculation @ track level for bremsstrahlung TODO
   
   // normalize probabilities
   Double_t probTotal = 0.0;
index 7a466ce..8f21bb9 100644 (file)
@@ -256,17 +256,27 @@ AliRieman* AliTRDtrackerV1::GetRiemanFitter()
 //_____________________________________________________________________________
 Int_t AliTRDtrackerV1::PropagateBack(AliESDEvent *event) 
 {
-  //
-  // Gets seeds from ESD event. The seeds are AliTPCtrack's found and
-  // backpropagated by the TPC tracker. Each seed is first propagated 
-  // to the TRD, and then its prolongation is searched in the TRD.
-  // If sufficiently long continuation of the track is found in the TRD
-  // the track is updated, otherwise it's stored as originaly defined 
-  // by the TPC tracker.   
-  //  
-
-  // Calibration monitor
-  AliTRDCalibraFillHisto *calibra = AliTRDCalibraFillHisto::Instance();
+// Propagation of ESD tracks from TPC to TOF detectors and building of the TRD track. For building
+// a TRD track an ESD track is used as seed. The informations obtained on the TRD track (measured points,
+// covariance, PID, etc.) are than used to update the corresponding ESD track.
+// Each track seed is first propagated to the geometrical limit of the TRD detector. 
+// Its prolongation is searched in the TRD and if corresponding clusters are found tracklets are 
+// constructed out of them (see AliTRDseedV1::AttachClusters()) and the track is updated. 
+// Otherwise the ESD track is left unchanged.
+// 
+// The following steps are performed:
+// 1. Selection of tracks based on the variance in the y-z plane.
+// 2. Propagation to the geometrical limit of the TRD volume. If track propagation fails the AliESDtrack::kTRDStop is set.
+// 3. Prolongation inside the fiducial volume (see AliTRDtrackerV1::FollowBackProlongation()) and marking
+// the following status bits:
+//   - AliESDtrack::kTRDin - if the tracks enters the TRD fiducial volume
+//   - AliESDtrack::kTRDStop - if the tracks fails propagation
+//   - AliESDtrack::kTRDbackup - if the tracks fulfills chi2 conditions and qualify for refitting
+// 4. Writting to friends, PID, MC label, quality etc. Setting status bit AliESDtrack::kTRDout.
+// 5. Propagation to TOF. If track propagation fails the AliESDtrack::kTRDStop is set.
+//  
+
+  AliTRDCalibraFillHisto *calibra = AliTRDCalibraFillHisto::Instance(); // Calibration monitor
   if (!calibra) AliInfo("Could not get Calibra instance\n");
   
   // Define scalers
@@ -654,7 +664,14 @@ Int_t AliTRDtrackerV1::FollowBackProlongation(AliTRDtrackV1 &t)
 // 5. Register tracklet with the tracker and track; update pulls monitoring.
 //
 // Observation
-//   1. During the propagation a bit map is filled detailing the status of the track in each TRD chamber
+//   1. During the propagation a bit map is filled detailing the status of the track in each TRD chamber. The following errors are being registered for each tracklet:
+// - AliTRDtrackV1::kProlongation : track prolongation failed
+// - AliTRDtrackV1::kPropagation : track prolongation failed
+// - AliTRDtrackV1::kAdjustSector : failed during sector crossing
+// - AliTRDtrackV1::kSnp : too large bending
+// - AliTRDtrackV1::kTrackletInit : fail to initialize tracklet
+// - AliTRDtrackV1::kUpdate : fail to attach clusters or fit the tracklet
+// - AliTRDtrackV1::kUnknown : anything which is not covered before
 //   2. By default the status of the track before first TRD update is saved. 
 // 
 // Debug level 2
@@ -2352,7 +2369,9 @@ Double_t AliTRDtrackerV1::BuildSeedingConfigs(AliTRDtrackingChamber **stack, Int
 Int_t AliTRDtrackerV1::MakeSeeds(AliTRDtrackingChamber **stack, AliTRDseedV1 *sseed, Int_t *ipar)
 {
 //
-// Seed tracklets and build candidate TRD tracks.
+// Seed tracklets and build candidate TRD tracks. The procedure is used during barrel tracking to account for tracks which are 
+// either missed by TPC prolongation or conversions inside the TRD volume. 
+// For stand alone tracking the procedure is used to estimate all tracks measured by TRD. 
 //
 // Parameters :
 //   layers : Array of stack propagation layers containing clusters