1 /**************************************************************************
2 * Copyright(c) 1998-1999, 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 // Trigger class for ZDC
20 // ****************************************************************
24 #include "AliLoader.h"
25 #include "AliRunLoader.h"
26 #include "AliTriggerInput.h"
29 #include "AliZDCDigit.h"
30 #include "AliZDCTrigger.h"
32 //________________________________________________________________
33 ClassImp(AliZDCTrigger)
35 //________________________________________________________________
36 AliZDCTrigger::AliZDCTrigger() :
45 fZDCLeftSemiCentrCut(0),
46 fZDCRightSemiCentrCut(0),
53 SetZDCLeftEMDCuts(0,0);
54 SetZDCRightEMDCuts(0,0);
58 //________________________________________________________________
59 void AliZDCTrigger::CreateInputs()
63 // Do not create inputs again!!
64 if( fInputs.GetEntriesFast() > 0 ) return;
66 fInputs.AddLast(new AliTriggerInput("ZDC_MB_L1", "ZDC", 0));
67 fInputs.AddLast(new AliTriggerInput("ZDC_Semi_L1", "ZDC", 0));
68 fInputs.AddLast(new AliTriggerInput("ZDC_Cent_L1", "ZDC", 0));
69 fInputs.AddLast(new AliTriggerInput("ZDC_EMD_L1", "ZDC", 0));
72 //________________________________________________________________
73 void AliZDCTrigger::Trigger()
78 AliRunLoader *runLoader = gAlice->GetRunLoader();
80 AliLoader *aZDCLoader = runLoader->GetLoader("ZDCLoader");
81 aZDCLoader->LoadDigits("READ");
83 AliZDCDigit* pdigit = &digit;
84 TTree* tD = aZDCLoader->TreeD();
86 cerr<<"AliZDCTrigger: digits tree not found\n";
89 tD->SetBranchAddress("ZDC", &pdigit);
91 Float_t signalZNLeft[]={0,0}, signalZPLeft[]={0,0}, signalZDCLeftSum[]={0,0};
92 Float_t signalZNRight[]={0,0}, signalZPRight[]={0,0}, signalZDCRightSum[]={0,0};
93 Float_t signalZEMSum[]={0,0};
94 for(Int_t iDigit=0; iDigit<tD->GetEntries(); iDigit++){
98 if(digit.GetSector(0)==1)
99 for(Int_t i=0; i<2; i++){ //0=high range; 1=low range
100 signalZNLeft[i] += digit.GetADCValue(i);
101 signalZDCLeftSum[i] += digit.GetADCValue(i);
103 else if(digit.GetSector(0)==2)
104 for(Int_t i=0; i<2; i++){
105 signalZPLeft[i] += digit.GetADCValue(i);
106 signalZDCLeftSum[i] += digit.GetADCValue(i);
108 else if(digit.GetSector(0)==3)
109 for(Int_t i=0; i<2; i++) signalZEMSum[i] += digit.GetADCValue(i);
111 else if(digit.GetSector(0)==4)
112 for(Int_t i=0; i<2; i++){ //0=high range; 1=low range
113 signalZNRight[i] += digit.GetADCValue(i);
114 signalZDCRightSum[i] += digit.GetADCValue(i);
116 else if(digit.GetSector(0)==5)
117 for(Int_t i=0; i<2; i++){
118 signalZPRight[i] += digit.GetADCValue(i);
119 signalZDCRightSum[i] += digit.GetADCValue(i);
122 // *******************************************************************
123 if(signalZDCLeftSum[1]>fZDCLeftMBCut && signalZDCRightSum[1]>fZDCRightMBCut)
124 // *** ZDC minimum bias trigger
125 SetInput("ZDC_1_L1");
126 // *******************************************************************
127 if(signalZDCLeftSum[1]>fZDCLeftCentrCut && signalZDCLeftSum[1]<fZDCLeftSemiCentrCut &&
128 signalZDCRightSum[1]>fZDCRightCentrCut && signalZDCRightSum[1]<fZDCRightSemiCentrCut
129 && signalZEMSum[1]>fZEMCentrCut)
130 // *** ZDC semi-central (10-40%)
131 SetInput("ZDC_2_L1");
132 // *******************************************************************
133 if(signalZDCLeftSum[1]>fZDCLeftMinCut && signalZDCLeftSum[1]<fZDCLeftCentrCut &&
134 signalZDCRightSum[1]>fZDCRightMinCut && signalZDCRightSum[1]<fZDCRightCentrCut &&
135 signalZEMSum[1]>fZEMCentrCut)
136 // *** ZDC central (0-10%)
137 SetInput("ZDC_3_L1");
138 // *******************************************************************
139 if(signalZNLeft[0]>fZDCLeftEMDCuts[0] && signalZNLeft[0]<fZDCLeftEMDCuts[1] &&
140 signalZNRight[0]>fZDCRightEMDCuts[0] && signalZNRight[0]<fZDCRightEMDCuts[1] &&
141 signalZEMSum[1]<fZEMMinCut){ // *** 1n EMD trigger
142 SetInput("ZDC_EMD_L1");
147 //________________________________________________________________
148 void AliZDCTrigger::SetZDCLeftMinCut(Float_t ZDCLeftMinCut)
150 // Set default cut values for ZDC trigger
152 if(ZDCLeftMinCut) fZDCLeftMinCut = ZDCLeftMinCut;
153 else fZDCLeftMinCut = 800.;
155 //________________________________________________________________
156 void AliZDCTrigger::SetZDCRightMinCut(Float_t ZDCRightMinCut)
158 // Set default cut values for ZDC trigger
160 if(ZDCRightMinCut) fZDCRightMinCut = ZDCRightMinCut;
161 else fZDCRightMinCut = 800.;
164 //________________________________________________________________
165 void AliZDCTrigger::SetZEMMinCut(Float_t ZEMMinCut)
167 // Set default cut values for ZDC trigger
169 if(ZEMMinCut) fZEMMinCut = ZEMMinCut;
170 else fZEMMinCut = 80.;
172 //________________________________________________________________
173 void AliZDCTrigger::SetZDCLeftEMDCuts(Float_t* ZDCLeftEMDCuts)
175 // Set default cut values for ZDC trigger
177 if(ZDCLeftEMDCuts) for(int j=0; j<2; j++) fZDCLeftEMDCuts[j] = ZDCLeftEMDCuts[j];
179 fZDCLeftEMDCuts[0] = 600.;
180 fZDCLeftEMDCuts[1] = 1000.;
183 //________________________________________________________________
184 void AliZDCTrigger::SetZDCLeftEMDCuts(Float_t ZDCLeftEMDCutInf,
185 Float_t ZDCLeftEMDCutSup)
187 // Set default cut values for ZDC trigger
189 if(ZDCLeftEMDCutInf && ZDCLeftEMDCutSup){
190 fZDCLeftEMDCuts[0]=ZDCLeftEMDCutInf;
191 fZDCLeftEMDCuts[1]=ZDCLeftEMDCutSup;
194 fZDCLeftEMDCuts[0] = 600.;
195 fZDCLeftEMDCuts[1] = 1000.;
198 //________________________________________________________________
199 void AliZDCTrigger::SetZDCRightEMDCuts(Float_t* ZDCRightEMDCuts)
201 // Set default cut values for ZDC trigger
203 if(ZDCRightEMDCuts) for(int j=0; j<2; j++) fZDCRightEMDCuts[j] = ZDCRightEMDCuts[j];
205 fZDCRightEMDCuts[0] = 600.;
206 fZDCRightEMDCuts[1] = 1000.;
209 //________________________________________________________________
210 void AliZDCTrigger::SetZDCRightEMDCuts(Float_t ZDCRightEMDCutInf,
211 Float_t ZDCRightEMDCutSup)
213 // Set default cut values for ZDC trigger
215 if(ZDCRightEMDCutInf && ZDCRightEMDCutSup){
216 fZDCRightEMDCuts[0]=ZDCRightEMDCutInf;
217 fZDCRightEMDCuts[1]=ZDCRightEMDCutSup;
220 fZDCRightEMDCuts[0] = 600.;
221 fZDCRightEMDCuts[1] = 1000.;
224 //________________________________________________________________
225 void AliZDCTrigger::SetZDCLeftMBCut(Float_t ZDCLeftMBCut)
227 // Set default cut values for ZDC trigger
229 if(ZDCLeftMBCut) fZDCLeftMBCut = ZDCLeftMBCut;
230 else fZDCLeftMBCut = 800.;
232 //________________________________________________________________
233 void AliZDCTrigger::SetZDCRightMBCut(Float_t ZDCRightMBCut)
235 // Set default cut values for ZDC trigger
237 if(ZDCRightMBCut) fZDCRightMBCut = ZDCRightMBCut;
238 else fZDCRightMBCut = 800.;
240 //________________________________________________________________
241 void AliZDCTrigger::SetZDCLeftCentrCut(Float_t ZDCLeftCentrCut)
243 // Set default cut values for ZDC trigger
245 if(ZDCLeftCentrCut) fZDCLeftCentrCut = ZDCLeftCentrCut;
246 else fZDCLeftCentrCut = 10000.;
248 //________________________________________________________________
249 void AliZDCTrigger::SetZDCRightCentrCut(Float_t ZDCRightCentrCut)
251 // Set default cut values for ZDC trigger
253 if(ZDCRightCentrCut) fZDCRightCentrCut = ZDCRightCentrCut;
254 else fZDCRightCentrCut = 10000.;
256 //________________________________________________________________
257 void AliZDCTrigger::SetZDCLeftSemiCentrCut(Float_t ZDCLeftSemiCentrCut)
259 // Set default cut values for ZDC trigger
261 if(ZDCLeftSemiCentrCut) fZDCLeftSemiCentrCut = ZDCLeftSemiCentrCut;
262 else fZDCLeftSemiCentrCut = 18500.;
264 //________________________________________________________________
265 void AliZDCTrigger::SetZDCRightSemiCentrCut(Float_t ZDCRightSemiCentrCut)
267 // Set default cut values for ZDC trigger
269 if(ZDCRightSemiCentrCut) fZDCRightSemiCentrCut = ZDCRightSemiCentrCut;
270 else fZDCRightSemiCentrCut = 18500.;
272 //________________________________________________________________
273 void AliZDCTrigger::SetZEMCentrCut(Float_t ZEMCentrCut)
275 // Set default cut values for ZDC trigger
277 if(ZEMCentrCut) fZEMCentrCut = ZEMCentrCut;
278 else fZEMCentrCut = 210.;