STEER/ESD/AliESDMuonCluster.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 /// \class AliESDMuonCluster
  20 ///
  21 /// Class to describe the MUON clusters in the Event Summary Data
  22 ///
  23 /// \author Philippe Pillot, Subatech
  24 //-----------------------------------------------------------------------------
  25
  26 #include "AliESDEvent.h"
  27 #include "AliESDMuonCluster.h"
  28 #include "AliESDMuonPad.h"
  29
  30 #include "AliLog.h"
  31
  32 #include <TClonesArray.h>
  33 #include <Riostream.h>
  34
  35 using std::endl;
  36 using std::cout;
  37 /// \cond CLASSIMP
  38 ClassImp(AliESDMuonCluster)
  39 /// \endcond
  40
  41 //_____________________________________________________________________________
  42 AliESDMuonCluster::AliESDMuonCluster()
  43 : TObject(),
  44   fCharge(0.),
  45   fChi2(0.),
  46   fPads(0x0),
  47   fNPads(0),
  48   fPadsId(0x0),
  49   fLabel(-1)
  50 {
  51   /// default constructor
  52   fXYZ[0] = fXYZ[1] = fXYZ[2] = 0.;
  53   fErrXY[0] = fErrXY[1] = 0.;
  54 }
  55
  56 //_____________________________________________________________________________
  57 AliESDMuonCluster::AliESDMuonCluster (const AliESDMuonCluster& cluster)
  58 : TObject(cluster),
  59   fCharge(cluster.fCharge),
  60   fChi2(cluster.fChi2),
  61   fPads(0x0),
  62   fNPads(cluster.fNPads),
  63   fPadsId(0x0),
  64   fLabel(cluster.fLabel)
  65 {
  66   /// Copy constructor
  67   fXYZ[0] = cluster.fXYZ[0];
  68   fXYZ[1] = cluster.fXYZ[1];
  69   fXYZ[2] = cluster.fXYZ[2];
  70   fErrXY[0] = cluster.fErrXY[0];
  71   fErrXY[1] = cluster.fErrXY[1];
  72
  73   if (cluster.fPads) {
  74     fPads = new TClonesArray("AliESDMuonPad",cluster.fPads->GetEntriesFast());
  75     AliESDMuonPad *pad = (AliESDMuonPad*) cluster.fPads->First();
  76     while (pad) {
  77       new ((*fPads)[fPads->GetEntriesFast()]) AliESDMuonPad(*pad);
  78       pad = (AliESDMuonPad*) cluster.fPads->After(pad);
  79     }
  80   }
  81
  82   if (cluster.fPadsId) fPadsId = new TArrayI(*(cluster.fPadsId));
  83 }
  84
  85 //_____________________________________________________________________________
  86 AliESDMuonCluster& AliESDMuonCluster::operator=(const AliESDMuonCluster& cluster)
  87 {
  88   /// Equal operator
  89   if (this == &cluster) return *this;
  90
  91   TObject::operator=(cluster); // don't forget to invoke the base class' assignment operator
  92
  93   fXYZ[0] = cluster.fXYZ[0];
  94   fXYZ[1] = cluster.fXYZ[1];
  95   fXYZ[2] = cluster.fXYZ[2];
  96   fErrXY[0] = cluster.fErrXY[0];
  97   fErrXY[1] = cluster.fErrXY[1];
  98
  99   fCharge = cluster.fCharge;
 100   fChi2 = cluster.fChi2;
 101   fLabel = cluster.fLabel;
 102
 103   delete fPads;
 104   if (cluster.fPads) {
 105     fPads = new TClonesArray("AliESDMuonPad",cluster.fPads->GetEntriesFast());
 106     AliESDMuonPad *pad = (AliESDMuonPad*) cluster.fPads->First();
 107     while (pad) {
 108       new ((*fPads)[fPads->GetEntriesFast()]) AliESDMuonPad(*pad);
 109       pad = (AliESDMuonPad*) cluster.fPads->After(pad);
 110     }
 111   } else fPads = 0x0;
 112
 113   SetPadsId(cluster.fNPads, cluster.GetPadsId());
 114
 115   return *this;
 116 }
 117
 118 //_____________________________________________________________________________
 119 void AliESDMuonCluster::Copy(TObject &obj) const {
 120
 121   /// This overwrites the virtual TOBject::Copy()
 122   /// to allow run time copying without casting
 123   /// in AliESDEvent
 124
 125   if(this==&obj)return;
 126   AliESDMuonCluster *robj = dynamic_cast<AliESDMuonCluster*>(&obj);
 127   if(!robj)return; // not an AliESDMuonCluster
 128   *robj = *this;
 129
 130 }
 131
 132 //__________________________________________________________________________
 133 AliESDMuonCluster::~AliESDMuonCluster()
 134 {
 135   /// Destructor
 136   delete fPads;
 137   delete fPadsId;
 138 }
 139
 140 //__________________________________________________________________________
 141 void AliESDMuonCluster::Clear(Option_t* opt)
 142 {
 143   /// Clear arrays
 144   if (opt && opt[0] == 'C') {
 145     if (fPads) fPads->Clear("C");
 146   } else {
 147     delete fPads; fPads = 0x0;
 148   }
 149   delete fPadsId; fPadsId = 0x0;
 150   fNPads = 0;
 151 }
 152
 153 //_____________________________________________________________________________
 154 void AliESDMuonCluster::AddPadId(UInt_t padId)
 155 {
 156   /// Add the given pad Id to the list associated to the cluster
 157   if (!fPadsId) fPadsId = new TArrayI(10);
 158   if (fPadsId->GetSize() <= fNPads) fPadsId->Set(fNPads+10);
 159   fPadsId->AddAt(static_cast<Int_t>(padId), fNPads++);
 160 }
 161
 162 //_____________________________________________________________________________
 163 void AliESDMuonCluster::SetPadsId(Int_t nPads, const UInt_t *padsId)
 164 {
 165   /// Fill the list pads'Id associated to the cluster with the given list
 166
 167   if (nPads <= 0 || !padsId) {
 168     delete fPadsId;
 169     fPadsId = 0x0;
 170     fNPads = 0;
 171     return;
 172   }
 173
 174   if (!fPadsId) fPadsId = new TArrayI(nPads, reinterpret_cast<const Int_t*>(padsId));
 175   else fPadsId->Set(nPads, reinterpret_cast<const Int_t*>(padsId));
 176   fNPads = nPads;
 177
 178 }
 179
 180 //_____________________________________________________________________________
 181 void AliESDMuonCluster::MovePadsToESD(AliESDEvent &esd)
 182 {
 183   /// move the pads to the new ESD structure
 184   if (!fPads) return;
 185   for (Int_t i = 0; i < fPads->GetEntriesFast(); i++) {
 186     AliESDMuonPad *pad = static_cast<AliESDMuonPad*>(fPads->UncheckedAt(i));
 187     AliESDMuonPad *newPad = esd.NewMuonPad();
 188     *newPad = *pad;
 189     AddPadId(newPad->GetUniqueID());
 190   }
 191   delete fPads;
 192   fPads = 0x0;
 193 }
 194
 195 //_____________________________________________________________________________
 196 void AliESDMuonCluster::Print(Option_t */*option*/) const
 197 {
 198   /// print cluster content
 199   UInt_t cId = GetUniqueID();
 200
 201   cout<<Form("clusterID=%u (ch=%d, det=%d, index=%d)",
 202              cId,GetChamberId(),GetDetElemId(),GetClusterIndex())<<endl;
 203
 204   cout<<Form("  position=(%5.2f, %5.2f, %5.2f), sigma=(%5.2f, %5.2f, 0.0)",
 205              GetX(),GetY(),GetZ(),GetErrX(),GetErrY())<<endl;
 206
 207   cout<<Form("  charge=%5.2f, chi2=%5.2f, MClabel=%d", GetCharge(), GetChi2(), GetLabel())<<endl;
 208
 209   if (PadsStored()) {
 210     cout<<"  pad infos:"<<endl;
 211     for (Int_t iPad=0; iPad<GetNPads(); iPad++) cout<<"  "<<GetPadId(iPad)<<endl;
 212   }
 213 }
 214