From 79aa5519c1b762eeaeac889122837dcb07ba2ce8 Mon Sep 17 00:00:00 2001
From: morsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Wed, 19 Apr 2000 13:42:15 +0000
Subject: [PATCH] Macro to find raw clusters (JB, AM)

---
 RICH/RICHrawclusters.C | 132 +++++++++++++++++++++++++++++++++++++++++
 1 file changed, 132 insertions(+)
 create mode 100644 RICH/RICHrawclusters.C

diff --git a/RICH/RICHrawclusters.C b/RICH/RICHrawclusters.C
new file mode 100644
index 00000000000..348e7f6adff
--- /dev/null
+++ b/RICH/RICHrawclusters.C
@@ -0,0 +1,132 @@
+#include "iostream.h"
+
+void RICHrawclusters (Int_t evNumber1=0,Int_t evNumber2=0) 
+{
+/////////////////////////////////////////////////////////////////////////
+//   This macro is a small example of a ROOT macro
+//   illustrating how to read the output of GALICE
+//   and do some analysis.
+//   
+/////////////////////////////////////////////////////////////////////////
+
+// 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 = new TFile("galice.root","UPDATE");
+
+// 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","Alice test 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");
+    }
+    
+//
+// Set reconstruction models
+//
+// Get pointers to Alice detectors and Digits containers
+    AliRICH *RICH  = (AliRICH*) gAlice->GetModule("RICH");
+
+    RecModel1 = new AliRICHClusterFinder();
+    RecModel1->SetNperMax(90);
+    //RecModel1->SetClusterSize(12);
+    RecModel1->SetClusterSize(100);
+    RecModel1->SetDeclusterFlag(1);
+    RICH->SetReconstructionModel(0,RecModel1);
+    
+    RecModel2 = new AliRICHClusterFinder();
+    RecModel2->SetNperMax(90);
+    //RecModel2->SetClusterSize(12);
+    RecModel2->SetClusterSize(100);
+    RecModel2->SetDeclusterFlag(1);
+    RICH->SetReconstructionModel(1,RecModel2);
+ 
+    RecModel3 = new AliRICHClusterFinder();
+    RecModel3->SetNperMax(90);
+    //RecModel3->SetClusterSize(12);
+    RecModel3->SetClusterSize(100);
+    RecModel3->SetDeclusterFlag(1);
+    RICH->SetReconstructionModel(2,RecModel3);
+    
+    RecModel4 = new AliRICHClusterFinder();
+    RecModel4->SetNperMax(90);
+    //RecModel4->SetClusterSize(12);
+    RecModel4->SetClusterSize(100);
+    RecModel4->SetDeclusterFlag(1);
+    RICH->SetReconstructionModel(3,RecModel4);
+    
+    //RecModel5 = new AliRICHClusterFinderv0();
+    RecModel5 = new AliRICHClusterFinder();
+    RecModel5->SetNperMax(90);
+    //RecModel5->SetClusterSize(15);
+    RecModel5->SetClusterSize(100);
+    RecModel5->SetDeclusterFlag(1);
+    RICH->SetReconstructionModel(4,RecModel5);
+    
+    //RecModel6 = new AliRICHClusterFinderv0();
+    RecModel6 = new AliRICHClusterFinder();
+    RecModel6->SetNperMax(90);
+    //RecModel6->SetClusterSize(15);
+    RecModel6->SetClusterSize(100);
+    RecModel6->SetDeclusterFlag(1);
+    RICH->SetReconstructionModel(5,RecModel6);
+    
+    RecModel7 = new AliRICHClusterFinder();
+    RecModel7->SetNperMax(90);
+    //RecModel7->SetClusterSize(9);
+    RecModel7->SetClusterSize(100);
+    RecModel7->SetDeclusterFlag(1);
+    RICH->SetReconstructionModel(6,RecModel7);
+//
+//   Loop over events              
+//
+    Int_t Nh=0;
+    Int_t Nh1=0;
+    for (int nev=0; nev<= evNumber2; nev++) {
+	Int_t nparticles = gAlice->GetEvent(nev);
+	cout <<endl<< "Processing event:" << nev <<endl;
+	cout << "Particles       :" << nparticles <<endl;
+	if (nev < evNumber1) continue;
+	if (nparticles <= 0) return;
+	
+	TTree *TH = gAlice->TreeH();
+	Int_t ntracks = TH->GetEntries();
+	//cout<<"ntracks "<<ntracks<<endl;
+	
+	Int_t nbytes = 0;
+
+
+	TClonesArray *Particles = gAlice->Particles();
+	TTree *TD = gAlice->TreeD();
+	Int_t nent=TD->GetEntries();
+	//printf("Found %d entries in the tree (must be one per cathode per event!)\n",nent);
+	if (RICH) {
+	    //printf("Finding clusters for event:%d",nev);
+	    RICH->FindClusters(nev,nent-2);
+	}   // end if RICH
+    }   // event loop 
+    //file->ls();
+    file->Close();
+
+
+    //delete gAlice;
+    printf("\nEnd of Macro  *************************************\n");
+}
+
-- 
2.31.1