2 //get trigger decision and write it in TreeR of MUON.RecPoints.root
4 void MUONtrigger (char* filename="galice.root",
5 Int_t evNumber1=0, Int_t evNumber2=9999)
7 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
8 // Creating Run Loader and openning file containing Hits
9 AliRunLoader * RunLoader = AliRunLoader::Open(filename,"MUONFolder","UPDATE");
10 if (RunLoader ==0x0) {
11 printf(">>> Error : Error Opening %s file \n",filename);
15 // Loading AliRun master
16 RunLoader->UnloadgAlice();
17 RunLoader->LoadgAlice();
18 gAlice = RunLoader->GetAliRun();
20 // Loading MUON subsystem
21 AliMUON * MUON = (AliMUON *) gAlice->GetDetector("MUON");
22 AliLoader * MUONLoader = RunLoader->GetLoader("MUONLoader");
23 AliMUONData * muondata = MUON->GetMUONData();
24 muondata->SetLoader(MUONLoader);
26 Int_t ievent, nevents;
27 nevents = RunLoader->GetNumberOfEvents();
29 MUONLoader->LoadDigits("READ");
30 MUONLoader->LoadRecPoints("UPDATE");
31 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
32 if (evNumber2>nevents) evNumber2=nevents;
33 for (Int_t ievent=evNumber1; ievent<evNumber2; ievent++) { // event loop
34 printf("event %d\n",ievent);
35 RunLoader->GetEvent(ievent);
36 if (MUONLoader->TreeR() == 0x0) {
37 MUONLoader->MakeTree("R");
39 muondata->MakeBranch("GLT");
40 muondata->SetTreeAddress("D,GLT");
41 MUON->Trigger(ievent);
42 muondata->ResetDigits();
43 muondata->ResetTrigger();
45 MUONLoader->UnloadDigits();
46 MUONLoader->UnloadRecPoints();