Classes for trigger processing in simulation
[u/mrichter/AliRoot.git] / ZDC / AliZDCRecoParamPbPb.cxx
CommitLineData
7bff3766 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
28ClassImp(AliZDCRecoParamPbPb)
29
30//_____________________________________________________________________________
31AliZDCRecoParamPbPb::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//_____________________________________________________________________________
60AliZDCRecoParamPbPb::~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//_____________________________________________________________________________
80AliZDCRecoParamPbPb *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.);
7bff3766 110
111 return param;
112
113}
114
115//_____________________________________________________________________________
116void AliZDCRecoParamPbPb::PrintParameters() const
117{
118 //
119 // print reconstruction parameters
120 //
121 printf("\n\n\t AliZDCRecoParamPbPb -> parameters set for reconstruction\n");
7bff3766 122 printf("\t Functions for reconstruction of centrality varibles (Pb-Pb):\n");
123
124 fZNCen->Print("");
125 fZNPer->Print("");
126 fZPCen->Print("");
127 fZPPer->Print("");
128 fZDCCen->Print("");
129 fZDCPer->Print("");
130 fbCen->Print("");
131 fbPer->Print("");
132 fZEMn->Print("");
133 fZEMp->Print("");
134 fZEMsp->Print("");
135 fZEMb->Print("");
136
137 printf("\n ####### Parameters from EZDC vs. ZEM correlation ####### \n");
138 printf("\tZEMEndPoint = %1.2f, ZEMCutFraction = %1.2f \n"
139 " DZEMInf = %1.2f, DZEMSup = %1.2f\n",
140 fZEMEndValue, fZEMCutFraction, fDZEMInf, fDZEMSup);
141 printf("\n ####### Parameters from EZDC vs. Nspec correlation ####### \n");
142 printf("\tEZN1MaxValue = %1.2f, EZP1MaxValue = %1.2f, EZDC1MaxValue = %1.2f \n"
143 "\tEZN2MaxValue = %1.2f, EZP2MaxValue = %1.2f, EZDC2MaxValue = %1.2f \n\n",
144 fEZN1MaxValue, fEZP1MaxValue, fEZDC1MaxValue,
145 fEZN2MaxValue, fEZP2MaxValue, fEZDC2MaxValue);
146}