FMD/AliFMDAnaESD.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 //
  17 // Utility class for analysing ESD data.
  18 // This class does sharing and background correction
  19 //
  20 #include <AliESDFMD.h>
  21 #include "AliFMDAnaESD.h"
  22 #include "AliFMDAnaRing.h"
  23 #include <TMath.h>
  24 #include <TBrowser.h>
  25 // clude <AliLog.h>
  26
  27 //____________________________________________________________________
  28 AliFMDAnaESD::AliFMDAnaESD()
  29   : fNEvents(0)
  30 {
  31   // Constructor
  32   // Parameters:
  33   //     lower cut.
  34   //
  35   AddLoad(kESD);
  36   for (size_t i = 0; i < 5; i++) fRing[i] = 0;
  37 }
  38 //____________________________________________________________________
  39 Bool_t
  40 AliFMDAnaESD::Init()
  41 {
  42   // Called at beginning of run
  43   // Parameters:
  44   // Return
  45   //     @c false on error
  46   //
  47   for (int i = 0; i < 5; i++) if (fRing[i]) fRing[i]->Init();
  48   return AliFMDInput::Init();
  49 }
  50 //____________________________________________________________________
  51 Bool_t
  52 AliFMDAnaESD::Begin(Int_t ev)
  53 {
  54   // Begining of event
  55   // Parameters:
  56   //     ev Event number
  57   // Return
  58   //     @c false on error
  59   //
  60   fNEvents++;
  61   for (int i = 0; i < 5; i++) if (fRing[i]) fRing[i]->Begin();
  62   Bool_t ret = AliFMDInput::Begin(ev);
  63   return ret;
  64 }
  65 //____________________________________________________________________
  66 Bool_t
  67 AliFMDAnaESD::ProcessESDs()
  68 {
  69   // Loop over all ESD data, and call ProcessESD for each entry.
  70   // Parameters:
  71   // Return
  72   //      @c false on error
  73   //
  74
  75   // Process event summary data
  76   if (!fESD) return kFALSE;
  77   for (UShort_t det = 1; det <= 3; det++) {
  78     Char_t rings[] = { 'I', (det == 1 ? '\0' : 'O'), '\0' };
  79     for (Char_t* rng = rings; *rng != '\0'; rng++) {
  80       UShort_t   nsec = (*rng == 'I' ?  20 :  40);
  81       UShort_t   nstr = (*rng == 'I' ? 512 : 256);
  82       Int_t      ridx = FindRing(det, *rng);
  83       if (ridx < 0 || !fRing[ridx]) continue;
  84       for (UShort_t sec = 0; sec < nsec; sec++) {
  85         for (UShort_t str = 0; str < nstr; str++) {
  86           Float_t eta    = fESD->Eta(det,*rng,sec,str);
  87           Float_t mult   = fESD->Multiplicity(det,*rng,sec,str);
  88           Float_t multp1 = (str == nstr-1 ? 0 :
  89                             fESD->Multiplicity(det,*rng,sec,str+1));
  90           if (!fESD->IsAngleCorrected()) {
  91             double c=TMath::Abs(TMath::Cos(2.*TMath::ATan(TMath::Exp(-eta))));
  92             mult    *= c;
  93             multp1  *= c;
  94           }
  95
  96           // Check signal isn't pedestal
  97           if (mult < 0.001) continue;
  98
  99           // Phi
 100           Double_t dPhi = 2*TMath::Pi() / fRing[ridx]->NSeq();
 101           Double_t phi;
 102           if (det == 3) phi =  TMath::Pi() - (sec + .5) * dPhi;
 103           else          phi =  (sec + .5) * dPhi;
 104           if (phi < 0)  phi += 2 * TMath::Pi();
 105
 106           if (fRing[ridx]->ProcessESD(phi, eta, mult, multp1)) str++;
 107           // "mult" possibly updated
 108           Fill(phi, eta, mult);
 109         }
 110       }
 111     }
 112   }
 113   return kTRUE;
 114 }
 115 //____________________________________________________________________
 116 Bool_t
 117 AliFMDAnaESD::End()
 118 {
 119   for (int i = 0; i < 5; i++) if (fRing[i]) fRing[i]->End();
 120   return AliFMDInput::End();
 121 }
 122 //____________________________________________________________________
 123 Bool_t
 124 AliFMDAnaESD::Finish()
 125 {
 126   // Called at the end of run
 127   // Parameters:
 128   // Return
 129   //     @c false in case of errors
 130   //
 131   for (int i = 0; i < 5; i++) if (fRing[i]) fRing[i]->Finish();
 132   return AliFMDInput::Finish();
 133 }
 134 //__________________________________________________________________
 135 void
 136 AliFMDAnaESD::AddRing(AliFMDAnaRing* ring)
 137 {
 138   // Add a ring
 139   // Parameters:
 140   //     ring Ring object
 141   //
 142   fRing[FindRing(ring->Detector(),ring->Ring())] = ring;
 143 }
 144 //____________________________________________________________________
 145 Int_t
 146 AliFMDAnaESD::FindRing(UShort_t det, Char_t ring) const
 147 {
 148   // Find a ring index
 149   // Parameters:
 150   //     det  Detector number
 151   //     ring Ring id
 152   // Return
 153   //     Index of ring object
 154   //
 155   Int_t idx = -1;
 156   switch (det) {
 157   case 1: idx = 0; break;
 158   case 2: idx = 1 + (ring == 'O' || ring == 'o' ? 1 : 0); break;
 159   case 3: idx = 3 + (ring == 'O' || ring == 'o' ? 1 : 0); break;
 160   }
 161   return idx;
 162 }
 163
 164 //____________________________________________________________________
 165 void
 166 AliFMDAnaESD::Browse(TBrowser* b)
 167 {
 168   for (size_t i = 0; i < 5; i++)
 169     if (fRing[i]) b->Add(fRing[i], fRing[i]->Name());
 170 }
 171
 172 //____________________________________________________________________
 173 //
 174 // EOF
 175 //