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