1 /**************************************************************************
2 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
16 ///////////////////////////////////////////////////////////////////////////////
18 // Class with ZDC reconstruction parameters //
19 // Origin: Chiara.Oppedisano@to.infn.it //
21 ///////////////////////////////////////////////////////////////////////////////
25 #include "AliZDCRecoParam.h"
26 #include "AliZDCRecoParamPbPb.h"
28 ClassImp(AliZDCRecoParamPbPb)
30 //_____________________________________________________________________________
31 AliZDCRecoParamPbPb::AliZDCRecoParamPbPb() :
59 //_____________________________________________________________________________
60 AliZDCRecoParamPbPb::~AliZDCRecoParamPbPb()
79 //_____________________________________________________________________________
80 AliZDCRecoParamPbPb *AliZDCRecoParamPbPb::GetPbPbRecoParam()
83 // Makes default reconstruction parameters for low flux environment
85 AliZDCRecoParamPbPb *param = new AliZDCRecoParamPbPb();
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.);
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.);
111 param->SetBeamEnergy(2.7);
117 //_____________________________________________________________________________
118 void AliZDCRecoParamPbPb::PrintParameters() const
121 // print reconstruction parameters
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");
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);