RICH/RICHMerger.C

   1 // 0 = all
   2 // 1 = not pion
   3 // 2 = not kaon
   4 // 3 = not proton
   5 // 4 = not muon
   6 // 5 = not electron
   7 // 6 = not neutron
   8
   9
  10 Int_t particle_type=0;
  11
  12 #include "iostream.h"
  13
  14 void RICHMerger (Int_t evNumber1=0,Int_t evNumber2=0)
  15 {
  16 /////////////////////////////////////////////////////////////////////////
  17 //   This macro is a small example of a ROOT macro
  18 //   illustrating how to read the output of GALICE
  19 //   and do some analysis.
  20 //
  21 /////////////////////////////////////////////////////////////////////////
  22
  23 // Dynamically link some shared libs
  24
  25
  26   if (gClassTable->GetID("AliRun") < 0) {
  27       gROOT->LoadMacro("loadlibs.C");
  28       loadlibs();
  29   }else {
  30     delete gAlice;
  31     gAlice = 0;
  32   }
  33
  34   galice=0;
  35
  36 // Connect the Root Galice file containing Geometry, Kine and Hits
  37
  38    TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
  39    if (file) file->Close();
  40    file = new TFile("galice.root","UPDATE");
  41 //   file->ls();
  42 // Get AliRun object from file or create it if not on file
  43
  44
  45
  46    if (gClassTable->GetID("AliRun") < 0) {
  47         gROOT->LoadMacro("loadlibs.C");
  48         loadlibs();
  49     }
  50
  51
  52    if (!gAlice) {
  53        gAlice = (AliRun*)file->Get("gAlice");
  54        if (gAlice) printf("AliRun object found on file\n");
  55        if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
  56    } else {
  57       delete gAlice;
  58       gAlice = (AliRun*)file->Get("gAlice");
  59         if (gAlice) printf("AliRun object found on file\n");
  60         if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
  61     }
  62
  63    AliRICH *RICH  = (AliRICH*) gAlice->GetDetector("RICH");
  64
  65    AliRICHChamber*       iChamber;
  66
  67    printf("Generating tresholds...\n");
  68
  69    for(Int_t i=0;i<7;i++)
  70      {
  71        iChamber = &(RICH->Chamber(i));
  72        iChamber->GenerateTresholds();
  73      }
  74
  75    // creation
  76    AliRICHMerger* merger = new AliRICHMerger();
  77    // configuration
  78    merger->SetMode(1);
  79    merger->SetSignalEventNumber(0);
  80    merger->SetBackgroundEventNumber(0);
  81    merger->SetBackgroundFileName("bg.root");
  82
  83 // pass
  84    RICH->SetMerger(merger);
  85
  86 //
  87 // Event Loop
  88 //
  89    for (int nev=0; nev<= evNumber2; nev++) {
  90        Int_t nparticles = gAlice->GetEvent(nev);
  91        cout <<endl<< "Processing event:" <<nev<<endl;
  92        cout << "Particles       :" <<nparticles<<endl;
  93        if (nev < evNumber1) continue;
  94        if (nparticles <= 0) return;
  95        if (RICH)
  96          {
  97            //gAlice->MakeTree("D");
  98            //RICH->MakeBranch("D");
  99            merger->Digitise(nev, particle_type);
 100          }
 101        //if (RICH) gAlice->SDigits2Digits("RICH");
 102        //char hname[30];
 103        //sprintf(hname,"TreeD%d",nev);
 104        //gAlice->TreeD()->Write(hname);
 105        //gAlice->TreeD()->Reset();
 106    } // event loop
 107    file->Close();
 108
 109    //delete gAlice;
 110    printf("\nEnd of Macro  *************************************\n");
 111 }
 112
 113
 114
 115
 116
 117
 118
 119
 120
 121
 122
 123
 124
 125