MUON/AliMUONClusterFinderCOG.cxx

   1 /**************************************************************************
   2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3 *                                                                        *
   4 * Author: The ALICE Off-line Project.                                    *
   5 * Contributors are mentioned in the code where appropriate.              *
   6 *                                                                        *
   7 * Permission to use, copy, modify and distribute this software and its   *
   8 * documentation strictly for non-commercial purposes is hereby granted   *
   9 * without fee, provided that the above copyright notice appears in all   *
  10 * copies and that both the copyright notice and this permission notice   *
  11 * appear in the supporting documentation. The authors make no claims     *
  12 * about the suitability of this software for any purpose. It is          *
  13 * provided "as is" without express or implied warranty.                  *
  14 **************************************************************************/
  15
  16 // $Id$
  17
  18 #include "AliMUONClusterFinderCOG.h"
  19
  20 #include "AliLog.h"
  21 #include "AliMUONCluster.h"
  22 #include "AliMUONDigit.h"
  23 #include "AliMUONPad.h"
  24 #include "AliMUONPreClusterFinder.h"
  25 #include "AliMpArea.h"
  26 #include "TClonesArray.h"
  27 #include "TVector2.h"
  28
  29 /// \class AliMUONClusterFinderCOG
  30 ///
  31 /// A very basic (and mostly useless, probably) cluster finder.
  32 ///
  33 /// We use AliMUONPreClusterFinder to actually build the cluster,
  34 /// and then we simply use center-of-gravity to get the coordinates
  35 /// of the cluster.
  36 /// Only point to note is that we compute separately both
  37 /// cathodes when doing, in order to take the positions from the
  38 /// direction with the better resolution.
  39 ///
  40 /// \author Laurent Aphecetche
  41
  42 /// \cond CLASSIMP
  43 ClassImp(AliMUONClusterFinderCOG)
  44 /// \endcond
  45
  46 //_____________________________________________________________________________
  47 AliMUONClusterFinderCOG::AliMUONClusterFinderCOG()
  48 : AliMUONVClusterFinder(),
  49 fPreClusterFinder(0x0)
  50 {
  51   /// ctor
  52 }
  53
  54 //_____________________________________________________________________________
  55 AliMUONClusterFinderCOG::~AliMUONClusterFinderCOG()
  56 {
  57   /// dtor
  58   delete fPreClusterFinder;
  59 }
  60
  61 //_____________________________________________________________________________
  62 Bool_t
  63 AliMUONClusterFinderCOG::Prepare(const AliMpVSegmentation* segmentations[2],
  64                                        TClonesArray* digits[2])
  65 {
  66   /// Prepare for clustering
  67
  68   // Find out the DetElemId
  69   Int_t detElemId(-1);
  70
  71   for ( Int_t i = 0; i < 2; ++i )
  72   {
  73     AliMUONDigit* d = static_cast<AliMUONDigit*>(digits[i]->First());
  74     if (d)
  75     {
  76       detElemId = d->DetElemId();
  77       break;
  78     }
  79   }
  80
  81   if ( detElemId < 0 )
  82   {
  83     AliWarning("Could not find DE. Probably no digits at all ?");
  84     return kFALSE;
  85   }
  86
  87   delete fPreClusterFinder;
  88   fPreClusterFinder = new AliMUONPreClusterFinder;
  89   return fPreClusterFinder->Prepare(segmentations,digits);
  90 }
  91
  92 //_____________________________________________________________________________
  93 AliMUONCluster*
  94 AliMUONClusterFinderCOG::NextCluster()
  95 {
  96   /// Get next cluster
  97
  98   if ( !fPreClusterFinder ) return 0x0;
  99   AliMUONCluster* cluster = fPreClusterFinder->NextCluster();
 100   if ( cluster )
 101   {
 102     ComputePosition(*cluster);
 103
 104     if ( cluster->Charge() < 7 )
 105     {
 106       // skip that one
 107       return NextCluster();
 108     }
 109   }
 110   return cluster;
 111 }
 112
 113 //_____________________________________________________________________________
 114 void
 115 AliMUONClusterFinderCOG::ComputePosition(AliMUONCluster& cluster)
 116 {
 117   /// Compute a first estimate of cluster position by a basic center-of-gravity
 118
 119   Double_t xmin = 1E9;
 120   Double_t ymin = 1E9;
 121   Double_t xmax = -1E9;
 122   Double_t ymax = -1E9;
 123
 124   Double_t x[] = { 0.0, 0.0 };
 125   Double_t y[] = { 0.0, 0.0 };
 126
 127   Double_t xsize[] = { 0.0, 0.0 } ;
 128   Double_t ysize[] = { 0.0, 0.0 } ;
 129
 130   for ( Int_t cathode = 0; cathode < 2; ++cathode )
 131   {
 132     for ( Int_t i = 0; i < cluster.Multiplicity(); ++i )
 133     {
 134       AliMUONPad* pad = cluster.Pad(i);
 135       TVector2 padPosition = pad->Position();
 136       AliMpArea area(pad->Position(),pad->Dimensions());
 137       xmin = TMath::Min(area.LeftBorder(),xmin);
 138       xmax = TMath::Max(area.RightBorder(),xmax);
 139       ymin = TMath::Min(area.DownBorder(),ymin);
 140       ymax = TMath::Max(area.UpBorder(),ymax);
 141       if ( cathode == pad->Cathode() )
 142       {
 143         x[cathode] += padPosition.X()*pad->Charge();
 144         y[cathode] += padPosition.Y()*pad->Charge();
 145         xsize[cathode] += pad->Dimensions().X();
 146         ysize[cathode] += pad->Dimensions().Y();
 147       }
 148     }
 149     if ( cluster.Charge(cathode) )
 150     {
 151       x[cathode] /= cluster.Charge(cathode);
 152       y[cathode] /= cluster.Charge(cathode);
 153     }
 154     if ( cluster.Multiplicity(cathode) )
 155     {
 156       xsize[cathode] /= cluster.Multiplicity(cathode);
 157       ysize[cathode] /= cluster.Multiplicity(cathode);
 158     }
 159   }
 160
 161   Double_t xCOG = 0;
 162   Double_t yCOG = 0;
 163
 164   // take the positions from the direction with the better resolution
 165   xCOG = ( xsize[0] < xsize[1] ) ? x[0] : x[1];
 166   yCOG = ( ysize[0] < ysize[1] ) ? y[0] : y[1];
 167
 168   AliDebug(1,Form("Cluster mult %d (x,y)=(%e,%e) boundaries=(xmin,ymin,xmax,ymax)=(%e,%e,%e,%e)"
 169                   " (x0,y0,x1,y1)=(%e,%e,%e,%e) ",
 170                   cluster.Multiplicity(),xCOG,yCOG,xmin,ymin,xmax,ymax,
 171                   x[0],y[0],x[1],y[1]));
 172
 173   cluster.SetPosition(TVector2(xCOG,yCOG),cluster.Pad(0)->Dimensions()); // FIXME: what to put as an error here ?
 174 }
 175
 176
 177