Fixes for bug #52499: Field polarities inconsistiency
[u/mrichter/AliRoot.git] / STEER / AliCentralTrigger.cxx
CommitLineData
a5a091ce 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/* $Id$ */
17
18///////////////////////////////////////////////////////////////////////////////
19// //
20// This class for running the Central Trigger Processor //
21// //
22// //
51f6d619 23// Load Configuration //
24// Make a list the trigger detectors involved (from the classes) //
a5a091ce 25// For the each event //
26// Run the Trigger for the each detector //
27// Get the inputs //
51f6d619 28// Check the trigger classes //
a5a091ce 29// Create the class mask //
30// Save result //
31// //
32// //
33///////////////////////////////////////////////////////////////////////////////
34
35#include <TString.h>
36#include <TObjString.h>
37#include <TObjArray.h>
38#include <TStopwatch.h>
bacbe0fd 39#include <TFile.h>
40#include <TTree.h>
a5a091ce 41
42#include "AliLog.h"
43#include "AliRun.h"
44#include "AliRunLoader.h"
45#include "AliModule.h"
46
47#include "AliTriggerInput.h"
48#include "AliTriggerDetector.h"
51f6d619 49#include "AliTriggerConfiguration.h"
50#include "AliTriggerClass.h"
51#include "AliTriggerCluster.h"
a5a091ce 52#include "AliCentralTrigger.h"
bacbe0fd 53#include "AliDetectorEventHeader.h"
54#include "AliHeader.h"
55
51f6d619 56#include "AliCDBManager.h"
57#include "AliCDBPath.h"
58#include "AliCDBEntry.h"
a5a091ce 59
60ClassImp( AliCentralTrigger )
61
62//_____________________________________________________________________________
63AliCentralTrigger::AliCentralTrigger() :
64 TObject(),
75e3794b 65 fClassMask(0),
51f6d619 66 fClusterMask(0),
d3f1f232 67 fL0TriggerInputs(0),
68 fL1TriggerInputs(0),
69 fL2TriggerInputs(0),
51f6d619 70 fConfiguration(NULL)
a5a091ce 71{
72 // Default constructor
596735df 73 SetOwner();
a5a091ce 74}
75
76//_____________________________________________________________________________
51f6d619 77AliCentralTrigger::AliCentralTrigger( TString & config ) :
a5a091ce 78 TObject(),
75e3794b 79 fClassMask(0),
51f6d619 80 fClusterMask(0),
d3f1f232 81 fL0TriggerInputs(0),
82 fL1TriggerInputs(0),
83 fL2TriggerInputs(0),
51f6d619 84 fConfiguration(NULL)
a5a091ce 85{
86 // Default constructor
51f6d619 87 LoadConfiguration( config );
bacbe0fd 88}
89
90//_____________________________________________________________________________
a5a091ce 91AliCentralTrigger::~AliCentralTrigger()
92{
51f6d619 93 // Destructor
94 DeleteConfiguration();
bacbe0fd 95}
96
97//_____________________________________________________________________________
51f6d619 98void AliCentralTrigger::DeleteConfiguration()
bacbe0fd 99{
51f6d619 100 // Delete the active configuration
101 fClassMask = 0;
102 fClusterMask = 0;
d3f1f232 103 fL0TriggerInputs = 0;
104 fL1TriggerInputs = 0;
105 fL2TriggerInputs = 0;
596735df 106 if (fConfiguration) {
107 if (IsOwner()) delete fConfiguration;
108 fConfiguration = 0x0;
109 }
a5a091ce 110}
111
112//_____________________________________________________________________________
bacbe0fd 113void AliCentralTrigger::Reset()
114{
51f6d619 115 // Reset Class Mask and classes
bacbe0fd 116 fClassMask = 0;
51f6d619 117 fClusterMask = 0;
d3f1f232 118 fL0TriggerInputs = 0;
119 fL1TriggerInputs = 0;
120 fL2TriggerInputs = 0;
51f6d619 121
122 if (fConfiguration) {
123 const TObjArray& classesArray = fConfiguration->GetClasses();
124 Int_t nclasses = classesArray.GetEntriesFast();
125 for( Int_t j=0; j<nclasses; j++ ) {
126 AliTriggerClass* trclass = (AliTriggerClass*)classesArray.At( j );
127 trclass->Reset();
128 }
bacbe0fd 129 }
130}
131
132//_____________________________________________________________________________
133void AliCentralTrigger::MakeBranch( TString name, TTree * tree )
134{
135 // Make a branch to store only trigger class mask event by event
136
137 if( tree ) {
138 AliDebug( 1, "Got Tree from folder." );
139 TBranch* branch = tree->GetBranch( name );
140 if( branch == 0x0 ) {
141 AliDebug( 1, "Creating new branch" );
d3f1f232 142 branch = tree->Branch( name, &(this->fClassMask), "fClassMask/l:fClusterMask/i:fL0TriggerInputs/i:fL1TriggerInputs/i:fL2TriggerInputs/s" );
bacbe0fd 143 branch->SetAutoDelete( kFALSE );
144 }
145 else {
146 AliDebug( 1, "Got Branch from Tree" );
147 branch->SetAddress( &(this->fClassMask) );
148 }
149 }
150}
151
152//_____________________________________________________________________________
51f6d619 153Bool_t AliCentralTrigger::LoadConfiguration( TString & config )
a5a091ce 154{
51f6d619 155 // Load one and only one pre-created COnfiguration from database/file that match
156 // with the input string 'config'
157 // Ej: "p-p", "Pb-Pb" or "p-p-DIMUON CALIBRATION-CENTRAL-BARREL"
158
159 // Delete the active configuration, if any
160 DeleteConfiguration();
161
162 // Load the selected configuration
163 if (!config.IsNull()) {
164 fConfiguration = AliTriggerConfiguration::LoadConfiguration( config );
596735df 165 SetOwner();
51f6d619 166 if(fConfiguration)
167 return kTRUE;
168 else {
169 AliError( Form( "Valid TriggerConfiguration (%s) is not found ! Disabling the trigger simulation !", config.Data() ) );
170 return kFALSE;
171 }
172 }
173 else {
174 // Load one and only one trigger descriptor from CDB
7e88424f 175 AliInfo( "Getting trigger configuration from OCDB!" );
51f6d619 176
177 AliCDBPath path( "GRP", "CTP", "Config" );
178
179 AliCDBEntry *entry=AliCDBManager::Instance()->Get(path.GetPath());
596735df 180 SetOwner(kFALSE);
51f6d619 181 if( !entry ) AliFatal( "Couldn't load trigger description data from CDB!" );
182
183 fConfiguration = (AliTriggerConfiguration *)entry->GetObject();
184 if(fConfiguration)
185 return kTRUE;
186 else {
187 AliError( "No valid configuration is found in the CDB ! Disabling the trigger simulation !" );
188 return kFALSE;
189 }
a5a091ce 190 }
a5a091ce 191}
192
193//_____________________________________________________________________________
194TString AliCentralTrigger::GetDetectors()
195{
51f6d619 196 // return TString with the detectors (modules) to be used for triggering
a5a091ce 197
198 TString result;
199
51f6d619 200 if (fConfiguration)
201 result = fConfiguration->GetTriggeringModules();
a5a091ce 202
203 return result;
204}
205
206//_____________________________________________________________________________
8480396b 207Bool_t AliCentralTrigger::RunTrigger( AliRunLoader* runLoader, const char *detectors )
a5a091ce 208{
209 // run the trigger
bacbe0fd 210
51f6d619 211 if( !fConfiguration ) {
212 AliError( "No trigger configuration loaded, skipping trigger" );
bacbe0fd 213 return kFALSE;
214 }
215
216 TTree *tree = runLoader->TreeCT();
217 if( !tree ) {
51f6d619 218 AliError( "No folder with trigger loaded, skipping trigger" );
a5a091ce 219 return kFALSE;
220 }
221
222 TStopwatch stopwatch;
223 stopwatch.Start();
bacbe0fd 224
01eed601 225 AliInfo( Form(" Triggering Detectors: %s \n", GetDetectors().Data() ) );
226 AliInfo( Form(" Detectors with digits: %s \n", detectors ) );
227
a5a091ce 228 // Process each event
229 for( Int_t iEvent = 0; iEvent < runLoader->GetNumberOfEvents(); iEvent++ ) {
01eed601 230 AliInfo( Form("Processing event %d", iEvent) );
bacbe0fd 231 runLoader->GetEvent( iEvent );
a5a091ce 232 // Get detectors involve
233 TString detStr = GetDetectors();
8480396b 234 TString detWithDigits = detectors;
a5a091ce 235 TObjArray* detArray = runLoader->GetAliRun()->Detectors();
236 // Reset Mask
237 fClassMask = 0;
895affe8 238 fClusterMask = 0;
239 // Reset configuration object (inputs and classes)
240 fConfiguration->Reset();
bacbe0fd 241 TObjArray trgdetArray; // use as garbage collector
a5a091ce 242 for( Int_t iDet = 0; iDet < detArray->GetEntriesFast(); iDet++ ) {
243 AliModule* det = (AliModule*) detArray->At( iDet );
244 if( !det || !det->IsActive() ) continue;
8480396b 245 if( IsSelected(det->GetName(), detStr) &&
246 IsSelected(det->GetName(), detWithDigits) ) {
bacbe0fd 247
01eed601 248 AliDebug(1,Form("Triggering from digits for %s", det->GetName() ) );
a5a091ce 249 AliTriggerDetector* trgdet = det->CreateTriggerDetector();
a4d25339 250 trgdet->AssignInputs(fConfiguration->GetInputs());
a5a091ce 251 TStopwatch stopwatchDet;
252 stopwatchDet.Start();
253 trgdet->Trigger();
01eed601 254 AliDebug(1, Form("Execution time for %s: R:%.2fs C:%.2fs",
a5a091ce 255 det->GetName(), stopwatchDet.RealTime(), stopwatchDet.CpuTime() ) );
bacbe0fd 256
bacbe0fd 257 trgdetArray.AddLast( trgdet );
258
259 // Write trigger detector in Event folder in Digits file
260 TString loadername = det->GetName();
261 loadername.Append( "Loader" );
262 AliLoader * loader = runLoader->GetLoader( loadername );
263 if( loader ) {
264 AliDataLoader * dataLoader = loader->GetDigitsDataLoader();
265 if( !dataLoader->IsFileOpen() ) {
266 if( dataLoader->OpenFile( "UPDATE" ) ) {
267 AliWarning( Form( "\n\nCan't write trigger for %s\n", det->GetName() ) );
268 }
269 }
270 dataLoader->Cd();
271 if( gFile && !gFile->IsWritable() ) {
272 gFile->ReOpen( "UPDATE" );
273 dataLoader->Cd();
274 }
275 trgdet->Write( "Trigger", TObject::kOverwrite );
276 dataLoader->CloseFile();
277 }
278 else AliWarning( Form( "Not loader found for %s", det->GetName() ) );
a5a091ce 279 }
280 }
281
282 // Check trigger conditions and create the trigger class mask
51f6d619 283 TriggerClasses();
d3f1f232 284 // Calculate trigger Input pattern
285 TriggerInputs();
bacbe0fd 286
a5a091ce 287 // Clear trigger detectors
288 trgdetArray.SetOwner();
289 trgdetArray.Delete();
290
291 if( (detStr.CompareTo( "ALL" ) != 0) && !detStr.IsNull() ) {
292 AliError( Form("the following detectors were not found: %s",
293 detStr.Data()));
e229562b 294 //JF return kFALSE;
a5a091ce 295 }
296
bacbe0fd 297 // Save trigger mask
298 tree->Fill();
01eed601 299 AliInfo( Form("Trigger Class Mask:0x%X", fClassMask ) );
a5a091ce 300 } // end event loop
a5a091ce 301
bacbe0fd 302 Reset();
303// cout << endl << " Print " << endl;
304// Print();
305
306 // Write
307 runLoader->WriteTrigger( "OVERWRITE" );
a5a091ce 308
bacbe0fd 309 return kTRUE;
310}
d3f1f232 311//----------------------------------------------------------------------------
312void AliCentralTrigger::TriggerInputs()
313{
314 // Find which inputs are in configuration
315 // and calculate input pattern
316 fL0TriggerInputs=0;
317 fL1TriggerInputs=0;
318 fL2TriggerInputs=0;
319 if(fConfiguration){
320 const TObjArray& inputsArray = fConfiguration->GetInputs();
321 Int_t ninputs = inputsArray.GetEntriesFast();
322 for( Int_t j=0; j<ninputs; j++ ) {
323 AliTriggerInput* input = (AliTriggerInput*)inputsArray.At( j );
324 if(input->GetValue()){
325 UChar_t level=input->GetLevel();
326 if(level == 0) fL0TriggerInputs |= (input->GetMask());
327 else if(level == 1) fL1TriggerInputs |= (input->GetMask());
328 else if(level == 2) fL2TriggerInputs |= (input->GetMask());
329 else{
330 AliError(Form("Unknown input level:%c:",level));
331 }
332 }
333 }
334 }
335}
a5a091ce 336//_____________________________________________________________________________
51f6d619 337ULong64_t AliCentralTrigger::TriggerClasses()
a5a091ce 338{
51f6d619 339 // Check trigger conditions and create the trigger class
340 // and trigger cluster masks
7e88424f 341 fClassMask = 0;
342 fClusterMask = 0;
51f6d619 343 if (fConfiguration) {
344 const TObjArray& classesArray = fConfiguration->GetClasses();
345 Int_t nclasses = classesArray.GetEntriesFast();
346 for( Int_t j=0; j<nclasses; j++ ) {
347 AliTriggerClass* trclass = (AliTriggerClass*)classesArray.At( j );
348 trclass->Trigger( fConfiguration->GetInputs(), fConfiguration->GetFunctions() );
349 fClassMask |= trclass->GetValue();
7e88424f 350 if (trclass->GetStatus()) {
351 AliTriggerCluster *trclust = trclass->GetCluster();
352 fClusterMask |= AliDAQ::DetectorPattern(trclust->GetDetectorsInCluster());
353 }
51f6d619 354 }
355 }
356 return fClassMask;
a5a091ce 357}
358//_____________________________________________________________________________
51f6d619 359TObjArray* AliCentralTrigger::GetFiredClasses() const
a5a091ce 360{
361 // return only the true conditions
362
363 TObjArray* result = new TObjArray();
364
51f6d619 365 if (fConfiguration) {
366 const TObjArray& classesArray = fConfiguration->GetClasses();
367 Int_t nclasses = classesArray.GetEntriesFast();
368 for( Int_t j=0; j<nclasses; j++ ) {
369 AliTriggerClass* trclass = (AliTriggerClass*)classesArray.At( j );
370 if( trclass->GetStatus() ) result->AddLast( trclass );
371 }
a5a091ce 372 }
373
374 return result;
375}
376
a5a091ce 377//_____________________________________________________________________________
378void AliCentralTrigger::Print( const Option_t* ) const
379{
380 // Print
381 cout << "Central Trigger: " << endl;
382 cout << " Trigger Class Mask: 0x" << hex << fClassMask << dec << endl;
51f6d619 383 if (fConfiguration) fConfiguration->Print();
a5a091ce 384 cout << endl;
385}
386
387
388//////////////////////////////////////////////////////////////////////////////
389// Helper method
390
391//_____________________________________________________________________________
392Bool_t AliCentralTrigger::IsSelected( TString detName, TString& detectors ) const
393{
394 // check whether detName is contained in detectors
395 // if yes, it is removed from detectors
396
397 // check if all detectors are selected
398 if( (detectors.CompareTo("ALL") == 0 ) ||
399 detectors.BeginsWith("ALL ") ||
400 detectors.EndsWith(" ALL") ||
401 detectors.Contains(" ALL ") ) {
402 detectors = "ALL";
403 return kTRUE;
404 }
405
406 // search for the given detector
407 Bool_t result = kFALSE;
408 if( (detectors.CompareTo( detName ) == 0) ||
409 detectors.BeginsWith( detName+" " ) ||
410 detectors.EndsWith( " "+detName ) ||
411 detectors.Contains( " "+detName+" " ) ) {
412 detectors.ReplaceAll( detName, "" );
413 result = kTRUE;
414 }
415
416 // clean up the detectors string
417 while( detectors.Contains(" ") ) detectors.ReplaceAll( " ", " " );
418 while( detectors.BeginsWith(" ") ) detectors.Remove( 0, 1 );
419 while( detectors.EndsWith(" ") ) detectors.Remove( detectors.Length()-1, 1 );
420
421 return result;
422}
67327b72 423
424//_____________________________________________________________________________
7e88424f 425Bool_t AliCentralTrigger::CheckTriggeredDetectors() const
67327b72 426{
427 // Check the trigger mask, finds which trigger classes
428 // have been fired, load the corresponding trigger clusters and
429 // finally makes a list of the detectors that have been readout
7e88424f 430 // for each particular event. This list is then compared to the
431 // one stored in fClusterMask. Return value:
432 // true = two lists are equal
433 // false = two lists are not equal meaning wrong trigger config
434 // is loaded.
67327b72 435
436 if (!fConfiguration) {
7e88424f 437 AliError("The trigger confiration has not yet been loaded! Cross-check is not possible!");
438 return kFALSE;
67327b72 439 }
7e88424f 440 else {
67327b72 441
7e88424f 442 // Make a cross-check so that to exclude wrong trigger configuration used
443 // Loop over the trigger classes
444 UInt_t clusterMask = 0;
445 const TObjArray& classesArray = fConfiguration->GetClasses();
446 Int_t nclasses = classesArray.GetEntriesFast();
447 for( Int_t j=0; j<nclasses; j++ ) {
448 AliTriggerClass* trclass = (AliTriggerClass*)classesArray.At( j );
449 if (trclass->GetMask() & fClassMask) { // class was fired
450 AliTriggerCluster *trclust = trclass->GetCluster();
451 clusterMask |= AliDAQ::DetectorPattern(trclust->GetDetectorsInCluster());
67327b72 452 }
453 }
7e88424f 454 // Compare the stored cluster mask with the one
455 // that we get from trigger classes
7e88424f 456 if (clusterMask != fClusterMask) {
454d6738 457 if ((clusterMask & fClusterMask) == clusterMask) {
458 AliInfo(Form("Cluster mask from trigger classes (%x) and from data (%x) differ. Concurrent DAQ run(s) could be the reason.",
459 (UInt_t)clusterMask,(UInt_t)fClusterMask));
460 return kTRUE;
461 }
462 else {
463 AliError(Form("Wrong cluster mask from trigger classes (%x), expecting (%x)! Loaded trigger configuration is possibly wrong!",
464 (UInt_t)clusterMask,(UInt_t)fClusterMask));
465 return kFALSE;
466 }
7e88424f 467 }
67327b72 468 }
469
7e88424f 470 return kTRUE;
67327b72 471}