PHOS/AliPHOSClusterizerv2.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 //////////////////////////////////////////////////////////////////////////////
  19 //  Clusterization class for IHEP reconstruction.
  20 // Performs clusterization (collects neighbouring active cells)
  21 // It differs from AliPHOSClusterizerv1 in neighbour definition only
  22 //*-- Author: Boris Polichtchouk, IHEP
  23
  24 // --- ROOT system ---
  25 #include "TBenchmark.h"
  26
  27 // --- Standard library ---
  28
  29 // --- AliRoot header files ---
  30 #include "AliLog.h"
  31 #include "AliPHOSClusterizerv2.h"
  32 #include "AliPHOSEvalRecPoint.h"
  33 #include "AliPHOSRecCpvManager.h"
  34 #include "AliPHOSRecEmcManager.h"
  35 #include "AliPHOSGeometry.h"
  36 #include "AliPHOSQualAssDataMaker.h"
  37 #include "AliPHOSDigit.h"
  38
  39 ClassImp(AliPHOSClusterizerv2)
  40
  41 //____________________________________________________________________________
  42 AliPHOSClusterizerv2::AliPHOSClusterizerv2() : AliPHOSClusterizerv1()
  43 {}
  44
  45 //____________________________________________________________________________
  46 AliPHOSClusterizerv2::AliPHOSClusterizerv2(AliPHOSGeometry *geom):
  47 AliPHOSClusterizerv1(geom)
  48 {}
  49
  50 //____________________________________________________________________________
  51 void AliPHOSClusterizerv2::GetNumberOfClustersFound(int* numb) const
  52 {
  53   // Returns the number of found EMC and CPV rec.points
  54
  55   numb[0] = fEMCRecPoints->GetEntries();
  56   numb[1] = fCPVRecPoints->GetEntries();
  57 }
  58
  59 //____________________________________________________________________________
  60 void AliPHOSClusterizerv2::Digits2Clusters(Option_t* option)
  61 {
  62   // Steering method
  63
  64   if(strstr(option,"tim"))
  65     gBenchmark->Start("PHOSClusterizer");
  66
  67   if(strstr(option,"print"))
  68     Print() ;
  69 }
  70 //   AliPHOSGetter * gime = AliPHOSGetter::Instance() ;
  71
  72 //   TFolder* wPoolF =  gime->PhosLoader()->GetDetectorDataFolder();
  73
  74 //   TObjArray* wPool = new TObjArray(400);
  75 //   wPool->SetName("SmartPoints");
  76 //   wPoolF->Add(wPool);
  77 //   wPoolF->Add(this);
  78
  79 //   Int_t nevents = gime->MaxEvent() ;
  80 //   Int_t ievent ;
  81
  82 //   for(ievent = 0; ievent<nevents; ievent++) {
  83
  84 //     gime->Event(ievent,"D") ;
  85
  86 //     AliInfo(Form("MakeClusters invoked..")) ;
  87 //     MakeClusters() ;
  88 //     AliInfo(Form("MakeClusters done.")) ;
  89
  90
  91 //     //SmartRecPoints will communicate with wPool.
  92
  93 //     AliPHOSEvalRecPoint* rp=0;
  94
  95 //     // CPV reconstruction
  96
  97 //     AliPHOSRecCpvManager* recCpv = new AliPHOSRecCpvManager();
  98 //     wPoolF->Add(recCpv);
  99
 100 //     Int_t iPoint; //loop variable
 101
 102 //     for(iPoint=0; iPoint<gime->CpvRecPoints()->GetEntriesFast(); iPoint++) {
 103 //       rp = new AliPHOSEvalRecPoint(iPoint, kTRUE);
 104 //       rp->MakeJob();
 105 //     }
 106
 107 //     AliPHOSEvalRecPoint pt;
 108 //     pt.UpdateWorkingPool();
 109
 110 //     TObjArray * cpvRecPoints = gime->CpvRecPoints() ;
 111 //     Int_t nOldCpv = cpvRecPoints->GetEntries();
 112 //     cpvRecPoints->Delete();
 113 //     cpvRecPoints->Compress();
 114
 115 //     Int_t i; //loop variable
 116
 117 //     for(i=0; i<wPool->GetEntries(); i++)
 118 //       cpvRecPoints->Add(wPool->At(i));
 119
 120 //     wPool->Clear();
 121 //     wPool->Compress();
 122
 123 //     wPoolF->Remove(recCpv);
 124 //     delete recCpv;
 125
 126 //     AliInfo(Form("       %d", gime->CpvRecPoints()->GetEntries() )) ;
 127 //     AliInfo(Form("       %d cpvRecPoints", cpvRecPoints->GetEntries() )) ;
 128
 129
 130 //     // Now Emc reconstruction
 131
 132 //     AliPHOSRecEmcManager* recEmc = new AliPHOSRecEmcManager();
 133 //     wPoolF->Add(recEmc);
 134
 135 //     for(iPoint=0; iPoint<gime->EmcRecPoints()->GetEntriesFast(); iPoint++) {
 136 //       rp = new AliPHOSEvalRecPoint(iPoint, kFALSE);
 137 //       rp->MakeJob();
 138 //     }
 139
 140 //     pt.UpdateWorkingPool();
 141
 142 //     TObjArray * emcRecPoints = gime->EmcRecPoints() ;
 143 //     Int_t nOldEmc = emcRecPoints->GetEntries();
 144 //     emcRecPoints->Delete();
 145 //     emcRecPoints->Compress();
 146
 147 //     for(i=0; i<wPool->GetEntries(); i++)
 148 //       emcRecPoints->Add(wPool->At(i));
 149
 150 //     wPool->Clear();
 151 //     wPool->Compress();
 152
 153 //     wPoolF->Remove(recEmc);
 154 //     delete recEmc;
 155
 156 //     TString message ;
 157 //     message  = "       %d  OLD cpvRecPoints\n" ;
 158 //     message += "       %d\n" ;
 159 //     message += "       %d cpvRecPoints\n" ;
 160
 161 //     message += "       %d OLD emcRecPoints " ;
 162 //     message += "       %d\n" ;
 163 //     message += "       %d emcRecPoints\n" ;
 164
 165 //     AliInfo(Form("%s", message.Data(),
 166 //       nOldCpv,
 167 //       gime->CpvRecPoints()->GetEntries(),cpvRecPoints->GetEntries(),
 168 //       nOldEmc,
 169 //       gime->EmcRecPoints()->GetEntries(), emcRecPoints->GetEntries() ));
 170
 171 //     GetQualAssDataMaker()->Init(AliQualAss::kRECPOINTS) ;
 172 //     GetQualAssDataMaker()->SetData(gime->EmcRecPoints()) ;
 173 //     GetQualAssDataMaker()->Exec(AliQualAss::kRECPOINTS) ;
 174 //     GetQualAssDataMaker()->SetData(gime->CpvRecPoints()) ;
 175 //     GetQualAssDataMaker()->Exec(AliQualAss::kRECPOINTS) ;
 176
 177 //     WriteRecPoints();
 178
 179
 180 //   } // loop over events
 181
 182 //   if(strstr(option,"tim")) {
 183 //     gBenchmark->Stop("PHOSClusterizer");
 184 //     AliInfo(Form("took %f seconds for Clusterizing", gBenchmark->GetCpuTime("PHOSClusterizer") )) ;
 185 //   }
 186 // }
 187
 188 //____________________________________________________________________________
 189 Int_t AliPHOSClusterizerv2::AreNeighbours(AliPHOSDigit* d1, AliPHOSDigit* d2) const
 190 {
 191   // Points are neighbours if they have common edge.
 192   // Points with common vertex are NOT neighbours.
 193   // This treatment of neighbourship is the part of
 194   // IHEP algorithm of clusterization.
 195
 196   // Gives the neighbourness of two digits = 0 are not neighbour but continue searching
 197   //                                       = 1 are neighbour
 198   //                                       = 2 are not neighbour but do not continue searching
 199   // The order of d1 and d2 is important: first (d1) should be a digit already in a cluster
 200   // which is compared to a digit (d2)  not yet in a cluster
 201
 202   Int_t rv = 0 ;
 203
 204   Int_t relid1[4] ;
 205   fGeom->AbsToRelNumbering(d1->GetId(), relid1) ;
 206
 207   Int_t relid2[4] ;
 208   fGeom->AbsToRelNumbering(d2->GetId(), relid2) ;
 209
 210   if ( (relid1[0] == relid2[0]) && (relid1[1]==relid2[1]) ) { // inside the same PHOS module and the same PPSD Module
 211     Int_t rowdiff = TMath::Abs( relid1[2] - relid2[2] ) ;
 212     Int_t coldiff = TMath::Abs( relid1[3] - relid2[3] ) ;
 213
 214     if ( ( (coldiff < 1) && (rowdiff <= 1) ) || ( ( coldiff <= 1 )  && ( rowdiff < 1 ) ) ){
 215       rv = 1 ;
 216     }
 217     else {
 218       if((relid2[2] > relid1[2]) && (relid2[3] > relid1[3]+1))
 219         rv = 2; //  Difference in row numbers is too large to look further
 220     }
 221
 222   }
 223   else {
 224
 225     if( (relid1[0] < relid2[0]) || (relid1[1] < relid2[1]) )
 226       rv=2 ;
 227
 228   }
 229
 230   return rv ;
 231
 232 }