FMD/FMDMerger.C

   1 void FMDMerger (Int_t evNumber1=0, Int_t evNumber2=0, Int_t ibg=0, Int_t bgr=10)
   2 {
   3 // Dynamically link some shared libs
   4
   5    if (gClassTable->GetID("AliRun") < 0) {
   6       gROOT->LoadMacro("loadlibs.C");
   7       loadlibs();
   8    }
   9
  10 // Connect the Root Galice file containing Geometry, Kine and Hits
  11
  12    TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("mgalice.root");
  13    if (file) file->Close();
  14    file = new TFile("mgalice.root","UPDATE");
  15    cout<<" Just starting... file "<<file<<endl;
  16
  17 // Get AliRun object from file or create it if not on file
  18
  19    if (!gAlice) {
  20        gAlice = (AliRun*)file->Get("gAlice");
  21        if (gAlice) printf("AliRun object found on file\n");
  22        if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
  23    }
  24    AliFMD *FMD  = (AliFMD*) gAlice->GetModule("FMD");
  25   cout<<" FMDMerger:: FMD "<<FMD<<" gAlice "<<gAlice<<endl;
  26
  27    if (FMD) {
  28 // creation
  29        AliFMDMerger* merger = new AliFMDMerger();
  30 // configuration
  31        if (ibg) {
  32          merger->SetMode(ibg);
  33          merger->SetBackgroundFileName("bg.root");
  34       }
  35        // pass
  36        //  merger->SetSignalFileName("mgalice.root");
  37        FMD->SetMerger(merger);
  38    }
  39 // Action !
  40 //
  41 //   Loop over events
  42 //
  43     for (int nev=evNumber1; nev<= evNumber2; nev++) {
  44         Int_t nparticles = gAlice->GetEvent(nev);
  45         cout << "nev         " << nev <<endl;
  46         cout << "nparticles  " << nparticles <<endl;
  47         if (nev < evNumber1) continue;
  48         if (nparticles <= 0) return;
  49         Int_t nbgr_ev = Int_t(nev*bgr/(evNumber2+1));
  50         cout<<" nbgr_ev "<<nbgr_ev<<endl;
  51
  52         if (ibg) {
  53             merger->SetBackgroundEventNumber(nbgr_ev);
  54         }
  55
  56         gAlice->SDigits2Digits("FMD");
  57
  58         char hname[30];
  59         sprintf(hname,"TreeD%d",nev);
  60         //      gAlice->TreeD()->Write(hname);
  61         //      cout<<hname<<" was written in file"<<file<<endl;
  62         //      gAlice->TreeD()->Print();
  63         //reset tree
  64         gAlice->TreeD()->Reset();
  65
  66     }   // event loop
  67 }
  68
  69
  70
  71
  72
  73
  74
  75
  76
  77
  78
  79
  80
  81