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 **************************************************************************/
19 /// \file MUONTrigger.C
20 /// \brief This macro is to be used to check the trigger algorithm w/o having to
21 /// (re-)perform simulation and digitalization.
23 /// See full description on the \ref README_trigger page.
25 /// \author P.Crochet (LPC)
27 #if !defined(__CINT__) || defined(__MAKECINT__)
30 #include "AliMUONDigit.h"
31 #include "AliMUONTriggerElectronics.h"
32 #include "AliMUONCalibrationData.h"
33 #include "AliCDBManager.h"
34 #include "AliMUONDataInterface.h"
35 #include "AliMUONMCDataInterface.h"
36 #include "AliMUONVTriggerStore.h"
37 #include "AliMUONDigitStoreV1.h"
38 #include <TClonesArray.h>
43 void MUONTrigger(const char* filename)
45 // Creating Run Loader and openning file containing Digits
46 AliRunLoader * RunLoader = AliRunLoader::Open(filename,"MUONLoader","UPDATE");
47 if (RunLoader ==0x0) {
48 printf(">>> Error : Error Opening %s file \n",filename);
51 // Loading AliRun master
52 if (RunLoader->GetAliRun() == 0x0) RunLoader->LoadgAlice();
53 gAlice = RunLoader->GetAliRun();
55 // Loading MUON subsystem
56 AliLoader* MUONLoader = RunLoader->GetDetectorLoader("MUON");
57 MUONLoader->LoadDigits("READ");
58 MUONLoader->LoadRecPoints("UPDATE"); // absolutely essential !!!
60 // Creating MUONTriggerDecision
61 AliCDBManager* cdbManager = AliCDBManager::Instance();
62 cdbManager->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
65 cdbManager->SetRun(runnumber);
66 AliMpCDB::LoadDDLStore();
68 AliMUONCalibrationData *CalibrationData = new AliMUONCalibrationData(runnumber);
69 AliMUONTriggerElectronics *TriggerProcessor = new AliMUONTriggerElectronics(CalibrationData);
71 Int_t nevents = RunLoader->GetNumberOfEvents();
72 AliMUONVDigitStore* digitStore=0x0;
73 AliMUONVTriggerStore* triggerStore=0x0;
75 for(Int_t ievent = 0; ievent < nevents; ievent++) {
76 printf(">>> Event %i out of %i \n",ievent,nevents);
77 RunLoader->GetEvent(ievent);
79 MUONLoader->LoadRecPoints("update");
80 MUONLoader->CleanRecPoints();
81 MUONLoader->MakeRecPointsContainer();
82 TTree* clustersTree = MUONLoader->TreeR();
83 TFile* cfile = clustersTree->GetCurrentFile();
86 cout << " could not find Cluster file " << endl;
90 MUONLoader->LoadDigits("read");
91 TTree* digitsTree = MUONLoader->TreeD();
92 TFile* dfile = digitsTree->GetCurrentFile();
95 cout << " could not find Digit file " << endl;
99 // here start reconstruction
100 if (!digitStore) digitStore = AliMUONVDigitStore::Create(*digitsTree);
101 if (!triggerStore) triggerStore = AliMUONVTriggerStore::Create(*digitsTree);
102 // insure we start with empty stores
106 Bool_t alone = ( triggerStore ? kFALSE : kTRUE );
107 Bool_t ok = digitStore->Connect(*digitsTree,alone);
110 cerr << "Could not connect digitStore to digitsTree \n";
114 cerr << "digitStore does not exist " << "\n";
118 digitsTree->GetEvent(0);
120 // process trigger response
121 TriggerProcessor->Digits2Trigger(*digitStore,*triggerStore);
123 //triggerStore->Print();
126 if ( triggerStore ) {
127 ok = triggerStore->Connect(*clustersTree,kTRUE);
130 cerr << "Could not create triggerStore branches in TreeR " << "\n";
134 cerr << "triggerStore does not exist " << "\n";
139 clustersTree->Fill();
140 MUONLoader->UnloadDigits();
141 MUONLoader->WriteRecPoints("OVERWRITE");
142 MUONLoader->UnloadRecPoints();