Add class to access simulation parameters, AliEMCALSimParam, to be set in configurati...
[u/mrichter/AliRoot.git] / ZDC / AliZDCTrigger.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, 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 //      Trigger class for ZDC
19 //
20 // ****************************************************************
21
22 #include "AliLog.h"
23 #include "AliRun.h"
24 #include "AliLoader.h"
25 #include "AliRunLoader.h"
26 #include "AliTriggerInput.h"
27
28 #include "AliZDC.h"
29 #include "AliZDCDigit.h"
30 #include "AliZDCTrigger.h"
31
32 //________________________________________________________________
33 ClassImp(AliZDCTrigger)
34
35 //________________________________________________________________
36 AliZDCTrigger::AliZDCTrigger() : 
37    AliTriggerDetector(), 
38    fZDCLeftMinCut(0),
39    fZDCRightMinCut(0),
40    fZEMMinCut(0),
41    fZDCLeftMBCut(0),
42    fZDCRightMBCut(0),
43    fZDCLeftCentrCut(0),
44    fZDCRightCentrCut(0),
45    fZDCLeftSemiCentrCut(0),
46    fZDCRightSemiCentrCut(0),
47    fZEMCentrCut(0)
48 {  
49    // Constructor
50    SetName("ZDC");
51    CreateInputs();
52    //
53    SetZDCLeftEMDCuts(0,0);
54    SetZDCRightEMDCuts(0,0);
55
56 }
57
58 //________________________________________________________________
59 void AliZDCTrigger::CreateInputs()
60 {
61    // Trigger inputs
62    
63    // Do not create inputs again!!
64    if( fInputs.GetEntriesFast() > 0 ) return;
65    
66    fInputs.AddLast(new AliTriggerInput("ZDC_1_L1",   "ZDC", 1));
67    fInputs.AddLast(new AliTriggerInput("ZDC_2_L1", "ZDC", 1));
68    fInputs.AddLast(new AliTriggerInput("ZDC_3_L1", "ZDC", 1));
69    fInputs.AddLast(new AliTriggerInput("ZDC_EMD_L1",  "ZDC", 1));
70 }
71
72 //________________________________________________________________
73 void AliZDCTrigger::Trigger()
74 {
75
76    // Trigger selection
77    //
78    AliRunLoader *runLoader = AliRunLoader::Instance();
79
80    AliLoader *aZDCLoader = runLoader->GetLoader("ZDCLoader");
81    
82    aZDCLoader->LoadDigits("READ");
83    AliZDCDigit digit;
84    AliZDCDigit* pdigit = &digit;
85    TTree* tD = aZDCLoader->TreeD();
86    if (!tD) {
87      cerr<<"AliZDCTrigger: digits tree not found\n";
88      return;
89    }
90    tD->SetBranchAddress("ZDC", &pdigit);
91    //
92    Float_t signalZNLeft[]={0,0}, signalZPLeft[]={0,0}, signalZDCLeftSum[]={0,0};
93    Float_t signalZNRight[]={0,0}, signalZPRight[]={0,0}, signalZDCRightSum[]={0,0};
94    Float_t signalZEMSum[]={0,0};
95    for(Int_t iDigit=0; iDigit<tD->GetEntries(); iDigit++){
96       tD->GetEntry(iDigit);
97       //
98       // *** ZDC LEFT
99       if(digit.GetSector(0)==1)
100          for(Int_t i=0; i<2; i++){ //0=high range; 1=low range
101             signalZNLeft[i] += digit.GetADCValue(i);
102             signalZDCLeftSum[i] += digit.GetADCValue(i);
103           }
104       else if(digit.GetSector(0)==2)
105          for(Int_t i=0; i<2; i++){
106             signalZPLeft[i] += digit.GetADCValue(i);
107             signalZDCLeftSum[i] += digit.GetADCValue(i);
108           }
109       else if(digit.GetSector(0)==3)
110          for(Int_t i=0; i<2; i++) signalZEMSum[i] += digit.GetADCValue(i);
111       // *** ZDC RIGHT
112       else if(digit.GetSector(0)==4)
113          for(Int_t i=0; i<2; i++){ //0=high range; 1=low range
114             signalZNRight[i] += digit.GetADCValue(i);
115             signalZDCRightSum[i] += digit.GetADCValue(i);
116           }
117       else if(digit.GetSector(0)==5)
118          for(Int_t i=0; i<2; i++){
119             signalZPRight[i] += digit.GetADCValue(i);
120             signalZDCRightSum[i] += digit.GetADCValue(i);
121           }
122    }
123    // *******************************************************************
124    if(signalZDCLeftSum[1]>fZDCLeftMBCut && signalZDCRightSum[1]>fZDCRightMBCut) 
125        // *** ZDC minimum bias trigger
126        SetInput("ZDC_1_L1");
127    // *******************************************************************
128    if(signalZDCLeftSum[1]>fZDCLeftCentrCut && signalZDCLeftSum[1]<fZDCLeftSemiCentrCut &&
129       signalZDCRightSum[1]>fZDCRightCentrCut && signalZDCRightSum[1]<fZDCRightSemiCentrCut
130       && signalZEMSum[1]>fZEMCentrCut) 
131        // *** ZDC semi-central (10-40%)
132        SetInput("ZDC_2_L1");
133    // *******************************************************************
134    if(signalZDCLeftSum[1]>fZDCLeftMinCut && signalZDCLeftSum[1]<fZDCLeftCentrCut &&
135       signalZDCRightSum[1]>fZDCRightMinCut && signalZDCRightSum[1]<fZDCRightCentrCut &&
136       signalZEMSum[1]>fZEMCentrCut) 
137        // *** ZDC central (0-10%)
138        SetInput("ZDC_3_L1");
139    // *******************************************************************
140    if(signalZNLeft[0]>fZDCLeftEMDCuts[0] && signalZNLeft[0]<fZDCLeftEMDCuts[1] && 
141       signalZNRight[0]>fZDCRightEMDCuts[0] && signalZNRight[0]<fZDCRightEMDCuts[1] &&
142       signalZEMSum[1]<fZEMMinCut){ // *** 1n EMD trigger
143         SetInput("ZDC_EMD_L1");
144    }
145    
146 }
147
148 //________________________________________________________________
149 void AliZDCTrigger::SetZDCLeftMinCut(Float_t ZDCLeftMinCut) 
150 {
151   // Set default cut values for ZDC trigger
152   //
153   if(ZDCLeftMinCut)  fZDCLeftMinCut = ZDCLeftMinCut;
154   else  fZDCLeftMinCut = 800.;
155 }
156 //________________________________________________________________
157 void AliZDCTrigger::SetZDCRightMinCut(Float_t ZDCRightMinCut) 
158 {
159   // Set default cut values for ZDC trigger
160   //
161   if(ZDCRightMinCut)  fZDCRightMinCut = ZDCRightMinCut;
162   else  fZDCRightMinCut = 800.;
163 }
164
165 //________________________________________________________________
166 void AliZDCTrigger::SetZEMMinCut(Float_t ZEMMinCut) 
167 {
168   // Set default cut values for ZDC trigger
169   //
170   if(ZEMMinCut)  fZEMMinCut = ZEMMinCut;
171   else  fZEMMinCut = 80.;
172 }
173 //________________________________________________________________
174 void AliZDCTrigger::SetZDCLeftEMDCuts(Float_t* ZDCLeftEMDCuts) 
175 {
176   // Set default cut values for ZDC trigger
177   //
178   if(ZDCLeftEMDCuts) for(int j=0; j<2; j++) fZDCLeftEMDCuts[j] = ZDCLeftEMDCuts[j];
179   else{
180     fZDCLeftEMDCuts[0] = 600.;
181     fZDCLeftEMDCuts[1] = 1000.;
182   }
183 }
184 //________________________________________________________________
185 void AliZDCTrigger::SetZDCLeftEMDCuts(Float_t ZDCLeftEMDCutInf, 
186         Float_t ZDCLeftEMDCutSup) 
187 {
188   // Set default cut values for ZDC trigger
189   //
190   if(ZDCLeftEMDCutInf && ZDCLeftEMDCutSup){
191     fZDCLeftEMDCuts[0]=ZDCLeftEMDCutInf; 
192     fZDCLeftEMDCuts[1]=ZDCLeftEMDCutSup;
193   }     
194   else{
195     fZDCLeftEMDCuts[0] = 600.;
196     fZDCLeftEMDCuts[1] = 1000.;
197   }
198 }
199 //________________________________________________________________
200 void AliZDCTrigger::SetZDCRightEMDCuts(Float_t* ZDCRightEMDCuts) 
201 {
202   // Set default cut values for ZDC trigger
203   //
204   if(ZDCRightEMDCuts) for(int j=0; j<2; j++) fZDCRightEMDCuts[j] = ZDCRightEMDCuts[j];
205   else{
206     fZDCRightEMDCuts[0] = 600.;
207     fZDCRightEMDCuts[1] = 1000.;
208   }
209 }
210 //________________________________________________________________
211 void AliZDCTrigger::SetZDCRightEMDCuts(Float_t ZDCRightEMDCutInf, 
212         Float_t ZDCRightEMDCutSup) 
213 {
214   // Set default cut values for ZDC trigger
215   //
216   if(ZDCRightEMDCutInf && ZDCRightEMDCutSup){
217     fZDCRightEMDCuts[0]=ZDCRightEMDCutInf; 
218     fZDCRightEMDCuts[1]=ZDCRightEMDCutSup;
219   }     
220   else{
221     fZDCRightEMDCuts[0] = 600.;
222     fZDCRightEMDCuts[1] = 1000.;
223   }
224 }
225 //________________________________________________________________
226 void AliZDCTrigger::SetZDCLeftMBCut(Float_t ZDCLeftMBCut) 
227 {
228   // Set default cut values for ZDC trigger
229   //
230   if(ZDCLeftMBCut) fZDCLeftMBCut = ZDCLeftMBCut;
231   else fZDCLeftMBCut = 800.;
232 }
233 //________________________________________________________________
234 void AliZDCTrigger::SetZDCRightMBCut(Float_t ZDCRightMBCut) 
235 {
236   // Set default cut values for ZDC trigger
237   //
238   if(ZDCRightMBCut) fZDCRightMBCut = ZDCRightMBCut;
239   else fZDCRightMBCut = 800.;
240 }
241 //________________________________________________________________
242 void AliZDCTrigger::SetZDCLeftCentrCut(Float_t ZDCLeftCentrCut) 
243 {
244   // Set default cut values for ZDC trigger
245   //
246   if(ZDCLeftCentrCut) fZDCLeftCentrCut = ZDCLeftCentrCut;
247   else fZDCLeftCentrCut = 10000.;
248 }
249 //________________________________________________________________
250 void AliZDCTrigger::SetZDCRightCentrCut(Float_t ZDCRightCentrCut) 
251 {
252   // Set default cut values for ZDC trigger
253   //
254   if(ZDCRightCentrCut) fZDCRightCentrCut = ZDCRightCentrCut;
255   else fZDCRightCentrCut = 10000.;
256 }
257 //________________________________________________________________
258 void AliZDCTrigger::SetZDCLeftSemiCentrCut(Float_t ZDCLeftSemiCentrCut) 
259 {
260   // Set default cut values for ZDC trigger
261   //
262   if(ZDCLeftSemiCentrCut) fZDCLeftSemiCentrCut = ZDCLeftSemiCentrCut;
263   else fZDCLeftSemiCentrCut = 18500.;
264 }
265 //________________________________________________________________
266 void AliZDCTrigger::SetZDCRightSemiCentrCut(Float_t ZDCRightSemiCentrCut) 
267 {
268   // Set default cut values for ZDC trigger
269   //
270   if(ZDCRightSemiCentrCut) fZDCRightSemiCentrCut = ZDCRightSemiCentrCut;
271   else fZDCRightSemiCentrCut = 18500.;
272 }
273 //________________________________________________________________
274 void AliZDCTrigger::SetZEMCentrCut(Float_t ZEMCentrCut) 
275 {
276   // Set default cut values for ZDC trigger
277   //
278   if(ZEMCentrCut) fZEMCentrCut = ZEMCentrCut;
279   else fZEMCentrCut = 210.;
280 }