]> git.uio.no Git - u/mrichter/AliRoot.git/blob - ITS/AliITSTriggerFOProcessor.cxx
load AddTaskEmcalJet.C
[u/mrichter/AliRoot.git] / ITS / AliITSTriggerFOProcessor.cxx
1 /////////////////////////////////////////////////////////////////////
2 // Author: Henrik Tydesjo                                          //
3 //                                                                 //
4 // This class takes care of simulating the output from the pixel   //
5 // trigger system.                                                 //
6 // The fast-or signals are given as input and they are processed   //
7 // to see which algorithm will give a signal to the central        //
8 // trigger system. To avoid doing the same operations many times,  //
9 // there is a method called PreprocessFOSignals, which should      //
10 // always be called for each event before the processing of        //
11 // each individual algorithm is done.                              //
12 //                                                                 //
13 // As soon as a new algorithm has been implemented online, a       //
14 // corresponding method should be implemented here. Which method   //
15 // should be used for a given trigger label is taken care of in    //
16 // ProcessFOSignalsLabel method.                                   //
17 //                                                                 //
18 /////////////////////////////////////////////////////////////////////
19
20 #include "AliITSTriggerFOProcessor.h"
21 #include "AliITSTriggerConditions.h"
22 #include <TError.h>
23
24 /* $Id$ */
25
26 AliITSTriggerFOProcessor::AliITSTriggerFOProcessor() :
27   fFOInner(0), fFOOuter(0), fTriggerCond(NULL)
28 {
29   // default constructor
30 }
31 //______________________________________________________________________
32 AliITSTriggerFOProcessor::AliITSTriggerFOProcessor(AliITSTriggerConditions* ocdbCond) :
33   fFOInner(0), fFOOuter(0), fTriggerCond(ocdbCond)
34 {
35   // constructor
36 }
37 //______________________________________________________________________
38 AliITSTriggerFOProcessor::AliITSTriggerFOProcessor(const AliITSTriggerFOProcessor& handle): 
39   fFOInner(handle.fFOInner), fFOOuter(handle.fFOOuter), fTriggerCond(handle.fTriggerCond)
40 {
41   // copy constructor
42 }
43 //______________________________________________________________________
44 AliITSTriggerFOProcessor::~AliITSTriggerFOProcessor() {
45   // destructor
46 }
47 //______________________________________________________________________
48 AliITSTriggerFOProcessor& AliITSTriggerFOProcessor::operator=(const AliITSTriggerFOProcessor& handle) {
49   // assignment operator
50   if (this!=&handle) {
51     fFOInner = handle.fFOInner;
52     fFOOuter = handle.fFOOuter;
53     fTriggerCond = handle.fTriggerCond;
54   }
55   return *this;
56 }
57 //______________________________________________________________________
58 void AliITSTriggerFOProcessor::SetTriggerConditions(AliITSTriggerConditions* const ocdbCond) {
59   // Method to give pointer to the OCDB conditions entry
60   fTriggerCond = ocdbCond;
61 }
62 //______________________________________________________________________
63 UInt_t AliITSTriggerFOProcessor::GetNumOutputs() const {
64   // return number of outputs (algorithms) in use
65   if (fTriggerCond==NULL) {
66     Error("AliITSTriggerFOProcessor::GetNumOutputs","Trigger conditions entry not yet given.");
67     return 0;
68   }
69   return fTriggerCond->GetNumAlgo();
70 }
71 //______________________________________________________________________
72 const Char_t* AliITSTriggerFOProcessor::GetOutputLabel(Short_t index) const {
73   // return the label of the index'th algorithm in use
74   if (fTriggerCond==NULL) {
75     Error("AliITSTriggerFOProcessor::GetOutputLabel","Trigger conditions entry not yet given.");
76     return "";
77   }
78   return fTriggerCond->GetAlgoLabelI(index);
79 }
80 //______________________________________________________________________
81 void AliITSTriggerFOProcessor::PreprocessFOSignals(AliITSFOSignalsSPD* const signals) {
82   // Pre-process the fast-or signals to retrieve some data needed by most algorithms
83   if (fTriggerCond==NULL) {
84     Error("AliITSTriggerFOProcessor::PreprocessFOSignals","Trigger conditions entry not yet given.");
85     return;
86   }
87   fFOInner = 0;
88   fFOOuter = 0;
89
90   Int_t eq   = -1;
91   Int_t hs   = -1;
92   Int_t chip = -1;
93   while (signals->GetNextSignal(eq,hs,chip)) {
94     if (fTriggerCond->IsChipActive( (UInt_t)eq, (UInt_t)hs, (UInt_t)chip) ) {
95       if (hs<=1) fFOInner++;
96       else       fFOOuter++;
97     }
98   }
99 }
100 //______________________________________________________________________
101 Bool_t AliITSTriggerFOProcessor::ProcessFOSignalsIndex(Short_t index, AliITSFOSignalsSPD* signals) {
102   // NB: For every event - Always call PreprocessFOSignals before calling this method
103   // Process index'th algorithm - returns true if output active
104   return ProcessFOSignalsLabel(GetOutputLabel(index), signals);
105 }
106 //______________________________________________________________________
107 Bool_t AliITSTriggerFOProcessor::ProcessFOSignalsLabel(const Char_t* label, AliITSFOSignalsSPD* signals) {
108   // NB: For every event - Always call PreprocessFOSignals before calling this method
109   // Process algorithm with label ... - returns true if output active 
110   if (fTriggerCond==NULL) {
111     Error("AliITSTriggerFOProcessor::ProcessFOSignalsLabel","Trigger conditions entry not yet given.");
112     return kFALSE;
113   }
114
115   Short_t index = fTriggerCond->GetAlgoIndexL(label);
116   if (index<0) {
117     Error("AliITSTriggerFOProcessor::ProcessFOSignalsLabel", "No conditions for label '%s'.",label);
118     return kFALSE;
119   }
120   
121   if      (strcmp(label, "0SMB") == 0) return ProcessFOSignalsTHRTotalAndTHRInnerAndTHROuter(index, signals);
122   else if (strcmp(label, "0SH1") == 0) return ProcessFOSignalsTHRInnerAndTHROuter(index, signals);
123   else if (strcmp(label, "0SH2") == 0) return ProcessFOSignalsTHRInnerAndTHROuter(index, signals);
124   else if (strcmp(label, "0SH3") == 0) return ProcessFOSignalsTHRInnerAndTHROuter(index, signals);
125   else if (strcmp(label, "0SH4") == 0) return ProcessFOSignalsTHRInnerAndTHROuter(index, signals);
126   else if (strcmp(label, "0SPF") == 0) return ProcessFOSignalsTHRTotalAndTHRInnerAndTHROuter(index, signals);
127   else if (strcmp(label, "0SX1") == 0) return ProcessFOSignalsInnerGTOuterPlusOffset(index, signals);
128   else if (strcmp(label, "0SX2") == 0) return ProcessFOSignalsOuterGTInnerPlusOffset(index, signals);
129   else if (strcmp(label, "0SBK") == 0) return ProcessFOSignalsTHRTotal(index, signals);
130   else if (strcmp(label, "0SCO") == 0) return ProcessFOSignalsCosmic(index, signals);
131
132   else {
133     Error("AliITSTriggerFOProcessor::ProcessFOSignalsLabel", "Algorithm not yet implemented for label '%s'.",label);
134     return kFALSE;
135   }
136 }
137 //______________________________________________________________________
138 Bool_t AliITSTriggerFOProcessor::ProcessFOSignalsTHRTotalAndTHRInnerAndTHROuter(Short_t index, AliITSFOSignalsSPD* /*signals*/) {
139   // NB: For every event - Always call PreprocessFOSignals before calling this method
140   // Process algorithm ' I+O > thr && I > thr && O > thr ' (index is needed to get the correct parameters from the ocdb object)
141   if (fTriggerCond==NULL) {
142     Error("AliITSTriggerFOProcessor::ProcessFOSignalsTHRTotalAndTHRInnerAndTHROuter","Trigger conditions entry not yet given.");
143     return kFALSE;
144   }
145   
146   // Get parameter values:
147   Int_t thIO = fTriggerCond->GetAlgoParamValueIN(index, "total_threshold");
148   if (thIO<0) {
149     Error("AliITSTriggerFOProcessor::ProcessFOSignalsTHRTotalAndTHRInnerAndTHROuter","Parameter 'total_threshold' not defined");
150     return kFALSE;
151   }
152   Int_t thI = fTriggerCond->GetAlgoParamValueIN(index, "inner_threshold");
153   if (thI<0) {
154     Error("AliITSTriggerFOProcessor::ProcessFOSignalsTHRTotalAndTHRInnerAndTHROuter","Parameter 'inner_threshold' not defined");
155     return kFALSE;
156   }
157   Int_t thO = fTriggerCond->GetAlgoParamValueIN(index, "outer_threshold");
158   if (thO<0) {
159     Error("AliITSTriggerFOProcessor::ProcessFOSignalsTHRTotalAndTHRInnerAndTHROuter","Parameter 'outer_threshold' not defined");
160     return kFALSE;
161   }
162
163   // Evaluate:
164   return (fFOInner+fFOOuter >= thIO && fFOInner >= thI && fFOOuter >= thO);
165 }
166 //______________________________________________________________________
167 Bool_t AliITSTriggerFOProcessor::ProcessFOSignalsTHRInnerAndTHROuter(Short_t index, AliITSFOSignalsSPD* /*signals*/) {
168   // NB: For every event - Always call PreprocessFOSignals before calling this method
169   // Process algorithm ' I > thr && O > thr ' (index is needed to get the correct parameters from the ocdb object)
170   if (fTriggerCond==NULL) {
171     Error("AliITSTriggerFOProcessor::ProcessFOSignalsTHRInnerAndTHROuter","Trigger conditions entry not yet given.");
172     return kFALSE;
173   }
174   
175   // Get parameter values:
176   Int_t thI = fTriggerCond->GetAlgoParamValueIN(index, "inner_threshold");
177   if (thI<0) {
178     Error("AliITSTriggerFOProcessor::ProcessFOSignalsTHRInnerAndTHROuter","Parameter 'inner_threshold' not defined");
179     return kFALSE;
180   }
181   Int_t thO = fTriggerCond->GetAlgoParamValueIN(index, "outer_threshold");
182   if (thO<0) {
183     Error("AliITSTriggerFOProcessor::ProcessFOSignalsTHRInnerAndTHROuter","Parameter 'outer_threshold' not defined");
184     return kFALSE;
185   }
186
187   // Evaluate:
188   return (fFOInner >= thI && fFOOuter >= thO);
189 }
190 //______________________________________________________________________
191 Bool_t AliITSTriggerFOProcessor::ProcessFOSignalsTHRTotal(Short_t index, AliITSFOSignalsSPD* /*signals*/) {
192   // NB: For every event - Always call PreprocessFOSignals before calling this method
193   // Process algorithm ' I+O > thr' (index is needed to get the correct parameters from the ocdb object)
194   if (fTriggerCond==NULL) {
195     Error("AliITSTriggerFOProcessor::ProcessFOSignalsTHRTotal","Trigger conditions entry not yet given.");
196     return kFALSE;
197   }
198   
199   // Get parameter values:
200   Int_t thIO = fTriggerCond->GetAlgoParamValueIN(index, "background_threshold_both");
201   if (thIO<0) {
202     Error("AliITSTriggerFOProcessor::ProcessFOSignalsTHRTotal","Parameter 'total_threshold' not defined");
203     return kFALSE;
204   }
205
206   // Evaluate:
207   return (fFOInner + fFOOuter >= thIO);
208 }
209 //______________________________________________________________________
210 Bool_t AliITSTriggerFOProcessor::ProcessFOSignalsInnerGTOuterPlusOffset(Short_t index, AliITSFOSignalsSPD* /*signals*/) {
211   // NB: For every event - Always call PreprocessFOSignals before calling this method
212   // Process algorithm ' I > O+offset ' (index is needed to get the correct parameters from the ocdb object)
213   if (fTriggerCond==NULL) {
214     Error("AliITSTriggerFOProcessor::ProcessFOSignalsGTOuterPlusOffset","Trigger conditions entry not yet given.");
215     return kFALSE;
216   }
217   
218   // Get parameter values:
219   Int_t offset = fTriggerCond->GetAlgoParamValueIN(index, "background_offset_inner");
220   if (offset<0) {
221     Error("AliITSTriggerFOProcessor::ProcessFOSignalsInnerGTOuterPlusOffset","Parameter 'offset' not defined");
222     return kFALSE;
223   }
224
225   // Evaluate:
226   return (fFOInner >=  fFOOuter + offset);
227 }
228 //______________________________________________________________________
229 Bool_t AliITSTriggerFOProcessor::ProcessFOSignalsOuterGTInnerPlusOffset(Short_t index, AliITSFOSignalsSPD* /*signals*/) {
230   // NB: For every event - Always call PreprocessFOSignals before calling this method
231   // Process algorithm ' O > I+offset ' (index is needed to get the correct parameters from the ocdb object)
232   if (fTriggerCond==NULL) {
233     Error("AliITSTriggerFOProcessor::ProcessFOSignalsGTOuterPlusOffset","Trigger conditions entry not yet given.");
234     return kFALSE;
235   }
236   
237   // Get parameter values:
238   Int_t offset = fTriggerCond->GetAlgoParamValueIN(index, "background_offset_outer");
239   if (offset<0) {
240     Error("AliITSTriggerFOProcessor::ProcessFOSignalsOuterGTInnerPlusOffset","Parameter 'offset' not defined");
241     return kFALSE;
242   }
243
244   // Evaluate:
245   return (fFOOuter >=  fFOInner + offset);
246 }
247 //______________________________________________________________________
248 Bool_t AliITSTriggerFOProcessor::ProcessFOSignalsCosmic(Short_t index, AliITSFOSignalsSPD* const signals) {
249   // NB: For every event - Always call PreprocessFOSignals before calling this method
250   // Process algorithm 'cosmic' (index is needed to get the correct parameters from the ocdb object)
251   if (fTriggerCond==NULL) {
252     Error("AliITSTriggerFOProcessor::ProcessFOSignalsCosmic","Trigger conditions entry not yet given.");
253     return kFALSE;
254   }
255
256   // Get parameter values:
257   Int_t cosmicParam = fTriggerCond->GetAlgoParamValueIN(index, "cosmic_mode");
258   if (cosmicParam<0) {
259     Error("AliITSTriggerFOProcessor::ProcessFOSignalsCosmic","Parameter 'cosmic_mode' not defined");
260     return kFALSE;
261   }
262
263   // Evaluate:
264
265   UShort_t topOuter = 0;
266   UShort_t topInner = 0;
267   UShort_t bottomOuter = 0;
268   UShort_t bottomInner = 0;
269
270   Int_t eq   = -1;
271   Int_t hs   = -1;
272   Int_t chip = -1;
273   while (signals->GetNextSignal(eq,hs,chip)) {
274     if (fTriggerCond->IsChipActive( (UInt_t)eq, (UInt_t)hs, (UInt_t)chip) ) {
275       if (hs<=1) {
276         if (eq%10 < 5) topInner++;
277         else        bottomInner++;
278       }
279       else {
280         if (eq%10 < 5) topOuter++;
281         else        bottomOuter++;
282       }
283     }
284   }
285
286   // top outer & bottom outer
287   if (cosmicParam == 0) return (topOuter>0 && bottomOuter>0);
288   // inner & outer
289   if (cosmicParam == 1) return (fFOInner>0 && fFOOuter>0);
290   // double layer ( >=2 of top inner, top outer, bottom inner, bottom outer )
291   if (cosmicParam == 2) {
292     UShort_t nHalfLayers = 0;
293     if (topOuter>0)    nHalfLayers++;
294     if (topInner>0)    nHalfLayers++;
295     if (bottomOuter>0) nHalfLayers++;
296     if (bottomInner>0) nHalfLayers++;
297     return (nHalfLayers>=2);
298   }
299   // top outer & top inner & bottom outer & bottom inner
300   if (cosmicParam == 3) return (topOuter>0 && topInner>0 && bottomOuter>0 && bottomInner>0);
301   // top outer & bottom outer & inner
302   if (cosmicParam == 4) return (topOuter>0 && bottomOuter>0 && fFOInner>0);
303   // global or
304   if (cosmicParam == 5) return (fFOOuter>0 || fFOInner>0);
305
306   Error("AliITSTriggerFOProcessor::ProcessFOSignalsCosmic","'cosmic_algorithm_parameter' = %d not defined",cosmicParam);
307   return kFALSE;
308 }
309