PWG2/FLOW/AliFlowCommon/AliFlowLYZHist1.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 //#include "Riostream.h"  //in case one wants to make print statements
  18 #include "TProfile.h"
  19 #include "TString.h"
  20 #include "TComplex.h"
  21 #include "TList.h"
  22 #include "AliFlowLYZConstants.h"
  23 #include "AliFlowLYZHist1.h"
  24 #include "TBrowser.h"
  25
  26 class TH1D;
  27
  28 // Class to organize the histograms in the first run
  29 // in the Lee Yang Zeros Flow analysis.
  30 // Also contains methods to find the first minimum R0
  31 // of the generating function.
  32 // author: N. van der Kolk (kolk@nikhef.nl)
  33
  34 ClassImp(AliFlowLYZHist1)
  35
  36 //-----------------------------------------------------------------------
  37
  38   AliFlowLYZHist1::AliFlowLYZHist1():
  39     TNamed(),
  40     fHistGtheta(0),
  41     fHistProReGtheta(0),
  42     fHistProImGtheta(0),
  43     fHistList(NULL)
  44 {
  45   //default constructor
  46 }
  47
  48 //-----------------------------------------------------------------------
  49
  50   AliFlowLYZHist1::AliFlowLYZHist1(Int_t theta, const char *anInput,Bool_t useSum):
  51     TNamed(anInput,anInput),
  52     fHistGtheta(0),
  53     fHistProReGtheta(0),
  54     fHistProImGtheta(0),
  55     fHistList(NULL)
  56 {
  57
  58   //constructor creating histograms
  59   Int_t iNbins      = AliFlowLYZConstants::GetMaster()->GetNbins();
  60   Double_t dMaxSUM  = AliFlowLYZConstants::GetMaster()->GetMaxSUM();
  61   Double_t dMaxPROD = AliFlowLYZConstants::GetMaster()->GetMaxPROD();
  62   Double_t dMin     = 0.;
  63
  64   TString name, addlast;
  65
  66   if (useSum) { addlast = "LYZSUM"; }
  67   else { addlast = "LYZPROD"; }
  68
  69   //fHistGtheta
  70   name = "First_Flow_Gtheta";
  71   name +=theta;
  72   name +=addlast;
  73   if (useSum) { fHistGtheta = new TH1D(name.Data(),name.Data(),iNbins,dMin,dMaxSUM); }
  74   else { fHistGtheta = new TH1D(name.Data(),name.Data(),iNbins,dMin,dMaxPROD); }
  75   fHistGtheta->SetXTitle("r");
  76   fHistGtheta->SetYTitle("|G^{#theta}(ir)|^{2}");
  77
  78   //fHistProReGtheta
  79   name = "First_FlowPro_ReGtheta";
  80   name +=theta;
  81   name +=addlast;
  82   if (useSum) { fHistProReGtheta = new TProfile(name.Data(),name.Data(),iNbins,dMin,dMaxSUM); }
  83   else { fHistProReGtheta = new TProfile(name.Data(),name.Data(),iNbins,dMin,dMaxPROD); }
  84   fHistProReGtheta->SetXTitle("r");
  85   fHistProReGtheta->SetYTitle("Re G^{#theta}(ir)");
  86
  87   //fHistProImGtheta
  88   name = "First_FlowPro_ImGtheta";
  89   name +=theta;
  90   name +=addlast;
  91   if (useSum) { fHistProImGtheta = new TProfile(name.Data(),name.Data(),iNbins,dMin,dMaxSUM); }
  92   else { fHistProImGtheta = new TProfile(name.Data(),name.Data(),iNbins,dMin,dMaxPROD); }
  93   fHistProImGtheta->SetXTitle("r");
  94   fHistProImGtheta->SetYTitle("Im G^{#theta}(ir)");
  95
  96   //list of histograms
  97   fHistList = new TList();
  98   fHistList-> Add(fHistGtheta);
  99   fHistList-> Add(fHistProReGtheta);
 100   fHistList-> Add(fHistProImGtheta);
 101
 102 }
 103
 104 //-----------------------------------------------------------------------
 105
 106 AliFlowLYZHist1::~AliFlowLYZHist1()
 107 {
 108   //deletes histograms
 109   delete fHistGtheta;
 110   delete fHistProReGtheta;
 111   delete fHistProImGtheta;
 112   delete fHistList;
 113 }
 114
 115 //-----------------------------------------------------------------------
 116
 117 void AliFlowLYZHist1::Fill(Double_t f, TComplex c)
 118 {
 119   //fill the histograms
 120
 121   fHistProReGtheta->Fill(f, c.Re());
 122   fHistProImGtheta->Fill(f, c.Im());
 123 }
 124
 125 //-----------------------------------------------------------------------
 126
 127 TH1D* AliFlowLYZHist1::FillGtheta()
 128 {
 129   //fill the fHistGtheta histograms
 130   Int_t iNbins = fHistGtheta->GetNbinsX();
 131   for (Int_t bin=1;bin<=iNbins;bin++)
 132         {
 133           //get bincentre of bins in histogram
 134           Double_t dBin = fHistGtheta->GetBinCenter(bin);
 135           Double_t dRe = fHistProReGtheta->GetBinContent(bin);
 136           Double_t dIm = fHistProImGtheta->GetBinContent(bin);
 137           TComplex cGtheta(dRe,dIm);
 138           //fill fHistGtheta with the modulus squared of cGtheta
 139           fHistGtheta->Fill(dBin,cGtheta.Rho2());
 140         }
 141
 142   return fHistGtheta;
 143 }
 144
 145 //-----------------------------------------------------------------------
 146
 147 Double_t AliFlowLYZHist1::GetR0()
 148 {
 149   //find the first minimum of the square of the modulus of Gtheta
 150
 151   Int_t iNbins = fHistGtheta->GetNbinsX();
 152   Double_t dR0 = 0.;
 153
 154   for (Int_t b=2;b<iNbins;b++)
 155     {
 156       Double_t dG0 = fHistGtheta->GetBinContent(b);
 157       Double_t dGnext = fHistGtheta->GetBinContent(b+1);
 158       Double_t dGnextnext = fHistGtheta->GetBinContent(b+2);
 159
 160       if (dGnext > dG0 && dGnextnext > dG0)
 161         {
 162           Double_t dGlast = fHistGtheta->GetBinContent(b-1);
 163           Double_t dXlast = fHistGtheta->GetBinCenter(b-1);
 164           Double_t dX0 = fHistGtheta->GetBinCenter(b);
 165           Double_t dXnext = fHistGtheta->GetBinCenter(b+1);
 166
 167           dR0 = dX0 - ((dX0-dXlast)*(dX0-dXlast)*(dG0-dGnext) - (dX0-dXnext)*(dX0-dXnext)*(dG0-dGlast))/
 168             (2.*((dX0-dXlast)*(dG0-dGnext) - (dX0-dXnext)*(dG0-dGlast))); //interpolated minimum
 169
 170           break; //stop loop if minimum is found
 171         } //if
 172
 173     }//b
 174
 175
 176   return dR0;
 177 }
 178
 179 //-----------------------------------------------------------------------
 180 Double_t AliFlowLYZHist1::GetBinCenter(Int_t i)
 181 {
 182   //gets bincenter of histogram
 183   Double_t dR = fHistGtheta->GetBinCenter(i);
 184   return dR;
 185 }
 186
 187 //-----------------------------------------------------------------------
 188
 189 Int_t AliFlowLYZHist1::GetNBins()
 190 {
 191   //gets iNbins
 192   Int_t iNbins = fHistGtheta->GetNbinsX();
 193   return iNbins;
 194 }
 195
 196 //-----------------------------------------------------------------------
 197  Double_t AliFlowLYZHist1::Merge(TCollection *aList)
 198 {
 199   //merge fuction
 200   if (!aList) return 0;
 201   if (aList->IsEmpty()) return 0; //no merging is needed
 202
 203   Int_t iCount = 0;
 204   TIter next(aList); // list is supposed to contain only objects of the same type as this
 205   AliFlowLYZHist1 *toMerge;
 206   // make a temporary list
 207   TList *pTemp = new TList();
 208   while ((toMerge=(AliFlowLYZHist1*)next())) {
 209     pTemp->Add(toMerge->GetHistList());
 210     iCount++;
 211   }
 212   // Now call merge for fHistList providing temp list
 213   fHistList->Merge(pTemp);
 214   // Cleanup
 215   delete pTemp;
 216
 217   return (double)iCount;
 218
 219 }
 220
 221 //-----------------------------------------------------------------------
 222 void AliFlowLYZHist1::Print(Option_t *option) const
 223 {
 224   //   -*-*-*-*-*Print some global quantities for this histogram collection class *-*-*-*-*-*-*-*
 225   //             ===============================================
 226   //   printf( "TH1.Print Name  = %s, Entries= %d, Total sum= %g\n",GetName(),Int_t(fEntries),GetSumOfWeights());
 227   printf( "Class.Print Name = %s, Histogram list:\n",GetName());
 228
 229   if (fHistList) {
 230     fHistList->Print(option);
 231   }
 232   else
 233     {
 234       printf( "Empty histogram list \n");
 235     }
 236 }
 237
 238 //-----------------------------------------------------------------------
 239  void AliFlowLYZHist1::Browse(TBrowser *b)
 240 {
 241
 242   if (!b) return;
 243   if (fHistList) b->Add(fHistList,"AliFlowLYZHist1List");
 244 }
 245
 246
 247