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() : AliTriggerDetector()
45 SetZDCLeftEMDCuts(0,0);
46 SetZDCRightEMDCuts(0,0);
49 SetZDCLeftCentrCut(0);
50 SetZDCRightCentrCut(0);
51 SetZDCLeftSemiCentrCut(0);
52 SetZDCRightSemiCentrCut(0);
57 //________________________________________________________________
58 void AliZDCTrigger::CreateInputs()
62 // Do not create inputs again!!
63 if( fInputs.GetEntriesFast() > 0 ) return;
65 fInputs.AddLast(new AliTriggerInput("ZDC_1_L1", "ZDC Minimum Bias", 0x01));
66 fInputs.AddLast(new AliTriggerInput("ZDC_2_L1", "ZDC Central", 0x02));
67 fInputs.AddLast(new AliTriggerInput("ZDC_3_L1", "ZDC Semi-central", 0x04));
68 fInputs.AddLast(new AliTriggerInput("ZDC_EMD_L1", "ZDC EMD events", 0x08));
71 //________________________________________________________________
72 void AliZDCTrigger::Trigger()
77 AliRunLoader *runLoader = gAlice->GetRunLoader();
79 AliLoader *aZDCLoader = runLoader->GetLoader("ZDCLoader");
80 aZDCLoader->LoadDigits("READ");
82 AliZDCDigit* pdigit = &digit;
83 TTree* tD = aZDCLoader->TreeD();
84 if (!tD) cerr<<"AliZDCTrigger: digits tree not found\n";
85 tD->SetBranchAddress("ZDC", &pdigit);
87 Float_t signalZNLeft[2], signalZPLeft[2], signalZDCLeftSum[2];
88 Float_t signalZNRight[2], signalZPRight[2], signalZDCRightSum[2];
89 Float_t signalZEMSum[2];
90 for(Int_t iDigit=0; iDigit<tD->GetEntries(); iDigit++){
94 if(digit.GetSector(0)==1)
95 for(Int_t i=0; i<2; i++){ //0=high range; 1=low range
96 signalZNLeft[i] += digit.GetADCValue(i);
97 signalZDCLeftSum[i] += digit.GetADCValue(i);
99 else if(digit.GetSector(0)==2)
100 for(Int_t i=0; i<2; i++){
101 signalZPLeft[i] += digit.GetADCValue(i);
102 signalZDCLeftSum[i] += digit.GetADCValue(i);
104 else if(digit.GetSector(0)==3)
105 for(Int_t i=0; i<2; i++) signalZEMSum[i] += digit.GetADCValue(i);
107 else if(digit.GetSector(0)==4)
108 for(Int_t i=0; i<2; i++){ //0=high range; 1=low range
109 signalZNRight[i] += digit.GetADCValue(i);
110 signalZDCRightSum[i] += digit.GetADCValue(i);
112 else if(digit.GetSector(0)==5)
113 for(Int_t i=0; i<2; i++){
114 signalZPRight[i] += digit.GetADCValue(i);
115 signalZDCRightSum[i] += digit.GetADCValue(i);
118 // *******************************************************************
119 if(signalZNLeft[0]>fZDCLeftEMDCuts[0] && signalZNLeft[0]<fZDCLeftEMDCuts[1] &&
120 signalZNRight[0]>fZDCRightEMDCuts[0] && signalZNRight[0]<fZDCRightEMDCuts[1] &&
121 signalZEMSum[1]<fZEMMinCut){ // *** 1n EMD trigger
122 SetInput("ZDC_EMD_L1");
124 // *******************************************************************
125 if(signalZDCLeftSum[1]>fZDCLeftMBCut && signalZDCRightSum[1]>fZDCRightMBCut)
126 // *** ZDC minimum bias trigger
127 SetInput("ZDC_1_L1");
128 // *******************************************************************
129 if(signalZDCLeftSum[1]>fZDCLeftMinCut && signalZDCLeftSum[1]<fZDCLeftCentrCut &&
130 signalZDCRightSum[1]>fZDCRightMinCut && signalZDCRightSum[1]<fZDCRightCentrCut &&
131 signalZEMSum[1]>fZEMCentrCut)
132 // *** ZDC central (0-10%)
133 SetInput("ZDC_2_L1");
134 // *******************************************************************
135 if(signalZDCLeftSum[1]>fZDCLeftCentrCut && signalZDCLeftSum[1]<fZDCLeftSemiCentrCut &&
136 signalZDCRightSum[1]>fZDCRightCentrCut && signalZDCRightSum[1]<fZDCRightSemiCentrCut
137 && signalZEMSum[1]>fZEMCentrCut)
138 // *** ZDC semi-central (10-40%)
139 SetInput("ZDC_3_L1");
143 //________________________________________________________________
144 void AliZDCTrigger::SetZDCLeftMinCut(Float_t ZDCLeftMinCut)
146 // Set default cut values for ZDC trigger
148 if(ZDCLeftMinCut) fZDCLeftMinCut = ZDCLeftMinCut;
149 else fZDCLeftMinCut = 800.;
151 //________________________________________________________________
152 void AliZDCTrigger::SetZDCRightMinCut(Float_t ZDCRightMinCut)
154 // Set default cut values for ZDC trigger
156 if(ZDCRightMinCut) fZDCRightMinCut = ZDCRightMinCut;
157 else fZDCRightMinCut = 800.;
160 //________________________________________________________________
161 void AliZDCTrigger::SetZEMMinCut(Float_t ZEMMinCut)
163 // Set default cut values for ZDC trigger
165 if(ZEMMinCut) fZEMMinCut = ZEMMinCut;
166 else fZEMMinCut = 80.;
168 //________________________________________________________________
169 void AliZDCTrigger::SetZDCLeftEMDCuts(Float_t* ZDCLeftEMDCuts)
171 // Set default cut values for ZDC trigger
173 if(ZDCLeftEMDCuts) for(int j=0; j<2; j++) fZDCLeftEMDCuts[j] = ZDCLeftEMDCuts[j];
175 fZDCLeftEMDCuts[0] = 600.;
176 fZDCLeftEMDCuts[1] = 1000.;
179 //________________________________________________________________
180 void AliZDCTrigger::SetZDCLeftEMDCuts(Float_t ZDCLeftEMDCutInf,
181 Float_t ZDCLeftEMDCutSup)
183 // Set default cut values for ZDC trigger
185 if(ZDCLeftEMDCutInf && ZDCLeftEMDCutSup){
186 fZDCLeftEMDCuts[0]=ZDCLeftEMDCutInf;
187 fZDCLeftEMDCuts[1]=ZDCLeftEMDCutSup;
190 fZDCLeftEMDCuts[0] = 600.;
191 fZDCLeftEMDCuts[1] = 1000.;
194 //________________________________________________________________
195 void AliZDCTrigger::SetZDCRightEMDCuts(Float_t* ZDCRightEMDCuts)
197 // Set default cut values for ZDC trigger
199 if(ZDCRightEMDCuts) for(int j=0; j<2; j++) fZDCRightEMDCuts[j] = ZDCRightEMDCuts[j];
201 fZDCRightEMDCuts[0] = 600.;
202 fZDCRightEMDCuts[1] = 1000.;
205 //________________________________________________________________
206 void AliZDCTrigger::SetZDCRightEMDCuts(Float_t ZDCRightEMDCutInf,
207 Float_t ZDCRightEMDCutSup)
209 // Set default cut values for ZDC trigger
211 if(ZDCRightEMDCutInf && ZDCRightEMDCutSup){
212 fZDCRightEMDCuts[0]=ZDCRightEMDCutInf;
213 fZDCRightEMDCuts[1]=ZDCRightEMDCutSup;
216 fZDCRightEMDCuts[0] = 600.;
217 fZDCRightEMDCuts[1] = 1000.;
220 //________________________________________________________________
221 void AliZDCTrigger::SetZDCLeftMBCut(Float_t ZDCLeftMBCut)
223 // Set default cut values for ZDC trigger
225 if(ZDCLeftMBCut) fZDCLeftMBCut = ZDCLeftMBCut;
226 else fZDCLeftMBCut = 800.;
228 //________________________________________________________________
229 void AliZDCTrigger::SetZDCRightMBCut(Float_t ZDCRightMBCut)
231 // Set default cut values for ZDC trigger
233 if(ZDCRightMBCut) fZDCRightMBCut = ZDCRightMBCut;
234 else fZDCRightMBCut = 800.;
236 //________________________________________________________________
237 void AliZDCTrigger::SetZDCLeftCentrCut(Float_t ZDCLeftCentrCut)
239 // Set default cut values for ZDC trigger
241 if(ZDCLeftCentrCut) fZDCLeftCentrCut = ZDCLeftCentrCut;
242 else fZDCLeftCentrCut = 10000.;
244 //________________________________________________________________
245 void AliZDCTrigger::SetZDCRightCentrCut(Float_t ZDCRightCentrCut)
247 // Set default cut values for ZDC trigger
249 if(ZDCRightCentrCut) fZDCRightCentrCut = ZDCRightCentrCut;
250 else fZDCRightCentrCut = 10000.;
252 //________________________________________________________________
253 void AliZDCTrigger::SetZDCLeftSemiCentrCut(Float_t ZDCLeftSemiCentrCut)
255 // Set default cut values for ZDC trigger
257 if(ZDCLeftSemiCentrCut) fZDCLeftSemiCentrCut = ZDCLeftSemiCentrCut;
258 else fZDCLeftSemiCentrCut = 18500.;
260 //________________________________________________________________
261 void AliZDCTrigger::SetZDCRightSemiCentrCut(Float_t ZDCRightSemiCentrCut)
263 // Set default cut values for ZDC trigger
265 if(ZDCRightSemiCentrCut) fZDCRightSemiCentrCut = ZDCRightSemiCentrCut;
266 else fZDCRightSemiCentrCut = 18500.;
268 //________________________________________________________________
269 void AliZDCTrigger::SetZEMCentrCut(Float_t ZEMCentrCut)
271 // Set default cut values for ZDC trigger
273 if(ZEMCentrCut) fZEMCentrCut = ZEMCentrCut;
274 else fZEMCentrCut = 210.;