MUON/MUONtrigger.C

   1 #include "iostream.h"
   2 //get trigger decision and write it in TreeR of MUON.RecPoints.root
   3
   4 void MUONtrigger (char* filename="galice.root",
   5                   Int_t evNumber1=0, Int_t evNumber2=0)
   6 {
   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);
  12     return;
  13   }
  14
  15   // Loading AliRun master
  16   RunLoader->UnloadgAlice();
  17   RunLoader->LoadgAlice();
  18   gAlice = RunLoader->GetAliRun();
  19
  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);
  25
  26   Int_t ievent, nevents;
  27   nevents = RunLoader->GetNumberOfEvents();
  28
  29   MUONLoader->LoadDigits("READ");
  30 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  31
  32    for (Int_t ievent=evNumber1; ievent<evNumber2; ievent++) { // event loop
  33        printf("event %d\n",ievent);
  34        RunLoader->GetEvent(ievent);
  35        if (MUONLoader->TreeR() == 0x0) MUONLoader->MakeTree("R");
  36        muondata->MakeBranch("GLT");
  37        muondata->SetTreeAddress("D,GLT");
  38        MUON->Trigger(ievent);
  39    } // event loop
  40    MUONLoader->UnloadDigits();
  41 }
  42
  43
  44
  45
  46
  47
  48
  49
  50
  51
  52
  53
  54
  55