ZDC/AliZDCRecoParamPbPb.cxx

   1 /**************************************************************************
   2  * Copyright(c) 2007-2009, 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 //                                                                           //
  18 // Class with ZDC reconstruction parameters                                  //
  19 // Origin: Chiara.Oppedisano@to.infn.it                                      //
  20 //                                                                           //
  21 ///////////////////////////////////////////////////////////////////////////////
  22
  23
  24 //#include <TF1.h>
  25 #include "AliZDCRecoParam.h"
  26 #include "AliZDCRecoParamPbPb.h"
  27
  28 ClassImp(AliZDCRecoParamPbPb)
  29
  30 //_____________________________________________________________________________
  31 AliZDCRecoParamPbPb::AliZDCRecoParamPbPb() :
  32   AliZDCRecoParam(),
  33   fZNCen(0),
  34   fZNPer(0),
  35   fZPCen(0),
  36   fZPPer(0),
  37   fZDCCen(0),
  38   fZDCPer(0),
  39   fbCen(0),
  40   fbPer(0),
  41   fZEMn(0),
  42   fZEMp(0),
  43   fZEMsp(0),
  44   fZEMb(0),
  45   fZEMEndValue(0),
  46   fZEMCutFraction(0),
  47   fDZEMSup(0),
  48   fDZEMInf(0),
  49   fEZN1MaxValue(0),
  50   fEZP1MaxValue(0),
  51   fEZDC1MaxValue(0),
  52   fEZN2MaxValue(0),
  53   fEZP2MaxValue(0),
  54   fEZDC2MaxValue(0)
  55 {
  56   //
  57   //Default constructor
  58 }
  59 //_____________________________________________________________________________
  60 AliZDCRecoParamPbPb::~AliZDCRecoParamPbPb()
  61 {
  62 // destructor
  63
  64   delete fZNCen;
  65   delete fZNPer;
  66   delete fZPCen;
  67   delete fZPPer;
  68   delete fZDCCen;
  69   delete fZDCPer;
  70   delete fbCen;
  71   delete fbPer;
  72   delete fZEMn;
  73   delete fZEMp;
  74   delete fZEMsp;
  75   delete fZEMb;
  76
  77 }
  78
  79 //_____________________________________________________________________________
  80 AliZDCRecoParamPbPb *AliZDCRecoParamPbPb::GetPbPbRecoParam()
  81 {
  82   //
  83   // Makes default reconstruction parameters for low flux environment
  84   //
  85   AliZDCRecoParamPbPb *param = new AliZDCRecoParamPbPb();
  86
  87   param->SetfZNCen("(-2.287920+sqrt(2.287920*2.287920-4*(-0.007629)*(11.921710-x)))/(2*(-0.007629))",0.,164.);
  88   param->SetfZNPer("(-37.812280-sqrt(37.812280*37.812280-4*(-0.190932)*(-1709.249672-x)))/(2*(-0.190932))",0.,164.);
  89   param->SetfZPCen("(-1.321353+sqrt(1.321353*1.321353-4*(-0.007283)*(3.550697-x)))/(2*(-0.007283))",0.,60.);
  90   param->SetfZPPer("(-42.643308-sqrt(42.643308*42.643308-4*(-0.310786)*(-1402.945615-x)))/(2*(-0.310786))",0.,60.);
  91   param->SetfZDCCen("(-1.934991+sqrt(1.934991*1.934991-4*(-0.004080)*(15.111124-x)))/(2*(-0.004080))",0.,225.);
  92   param->SetfZDCPer("(-34.380639-sqrt(34.380639*34.380639-4*(-0.104251)*(-2612.189017-x)))/(2*(-0.104251))",0.,225.);
  93   param->SetfbCen("-0.056923+0.079703*x-0.0004301*x*x+0.000001366*x*x*x",0.,220.);
  94   param->SetfbPer("17.943998-0.046846*x+0.000074*x*x",0.,220.);
  95   param->SetfZEMn("121.7-0.1934*x+0.00007565*x*x",0.,1200.);
  96   param->SetfZEMp("80.05-0.1315*x+0.00005327*x*x",0.,1200.);
  97   param->SetfZEMsp("201.7-0.325*x+0.0001292*x*x",0.,1200.);
  98   param->SetfZEMb("13.83-0.02851*x+5.101e-5*x*x-7.305e-8*x*x*x+5.101e-11*x*x*x*x-1.25e-14*x*x*x*x*x",0.,1200.);
  99
 100   param->SetZEMEndValue(1200.);
 101   param->SetZEMCutFraction(0.1);
 102   param->SetDZEMSup(0.04);
 103   param->SetDZEMInf(0.05);
 104   param->SetEZN1MaxValue(161.);
 105   param->SetEZP1MaxValue(59.);
 106   param->SetEZDC1MaxValue(220.);
 107   param->SetEZN2MaxValue(161.);
 108   param->SetEZP2MaxValue(59.);
 109   param->SetEZDC2MaxValue(161.);
 110   //
 111   param->SetBeamEnergy(2.7);
 112
 113   return param;
 114
 115 }
 116
 117 //_____________________________________________________________________________
 118 void AliZDCRecoParamPbPb::PrintParameters() const
 119 {
 120   //
 121   // print reconstruction parameters
 122   //
 123   printf("\n\n\t AliZDCRecoParamPbPb -> parameters set for reconstruction\n");
 124   printf("\t Beam energy (for Pb-Pb): %1.2f TeV\n",fBeamEnergy);
 125   printf("\t Functions for reconstruction of centrality varibles (Pb-Pb):\n");
 126
 127   fZNCen->Print("");
 128   fZNPer->Print("");
 129   fZPCen->Print("");
 130   fZPPer->Print("");
 131   fZDCCen->Print("");
 132   fZDCPer->Print("");
 133   fbCen->Print("");
 134   fbPer->Print("");
 135   fZEMn->Print("");
 136   fZEMp->Print("");
 137   fZEMsp->Print("");
 138   fZEMb->Print("");
 139
 140   printf("\n ####### Parameters from EZDC vs. ZEM correlation #######  \n");
 141   printf("\tZEMEndPoint = %1.2f, ZEMCutFraction = %1.2f \n"
 142     "  DZEMInf = %1.2f, DZEMSup = %1.2f\n",
 143     fZEMEndValue, fZEMCutFraction, fDZEMInf, fDZEMSup);
 144   printf("\n ####### Parameters from EZDC vs. Nspec correlation #######        \n");
 145   printf("\tEZN1MaxValue = %1.2f, EZP1MaxValue = %1.2f, EZDC1MaxValue = %1.2f \n"
 146     "\tEZN2MaxValue = %1.2f, EZP2MaxValue = %1.2f, EZDC2MaxValue = %1.2f \n\n",
 147     fEZN1MaxValue, fEZP1MaxValue, fEZDC1MaxValue,
 148     fEZN2MaxValue, fEZP2MaxValue, fEZDC2MaxValue);
 149 }