Possibility to stop tracking at a given layer (default is 1). Added new method AliITS...

[u/mrichter/AliRoot.git] / RICH / RICHdetect.C

RICHdetect (Int_t evNumber1=0,Int_t evNumber2=0, Int_t type, Int_t version) {
// Dynamically link some shared libs
 
    if (gClassTable->GetID("AliRun") < 0) {
	gROOT->LoadMacro("loadlibs.C");
	loadlibs();
    }
    else {
      //delete gAlice;
      gAlice = 0;
    }
// Connect the Root Galice file containing Geometry, Kine and Hits
    
    TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
    if (file) file->Close(); 
    file = new TFile("galice.root","UPDATE");
    file->ls();
    
    //printf ("I'm after Map \n");
    
// Get AliRun object from file or create it if not on file
    
    if (!gAlice) {
	gAlice = (AliRun*)file->Get("gAlice");
	if (gAlice) printf("AliRun object found on file\n");
	if (!gAlice) gAlice = new AliRun("gAlice","RICH reconstruction program");
    } else {
      delete gAlice;
      gAlice = (AliRun*)file->Get("gAlice");
      	if (gAlice) printf("AliRun object found on file\n");
	if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
    }
    
    //printf ("I'm after gAlice \n");
    
    AliRICH *RICH  = (AliRICH*) gAlice->GetDetector("RICH"); 

    // Create Recontruction algorithm object
    
    //V0
    if(version == 0)
      AliRICHDetect *detect = new AliRICHDetect("RICH reconstruction algorithm","Reconstruction");

    if(version == 1)
      AliRICHDetect *detect = new AliRICHDetectV1("RICH reconstruction algorithm","Reconstruction");
    

// Reconstruct
    // Event Loop
    //
    for (int nev=0; nev<= evNumber2; nev++) {
      Int_t nparticles = gAlice->GetEvent(nev);
      printf("\nProcessing event: %d\n",nev);
      printf("\nParticles       : %d\n",nparticles);
      if (nev < evNumber1) continue;
      if (nparticles <= 0) return;
      if (RICH) detect->Detect(nev, type);
      char hname[30];
      sprintf(hname,"TreeR%d",nev);
      gAlice->TreeR()->Write(hname);
      gAlice->TreeR()->Reset();
   } // event loop  
   
    //delete gAlice;
    printf("\nEnd of Macro  *************************************\n");
}