// Get pointers to Alice detectors and Digits containers
AliMUON *MUON = (AliMUON*) gAlice->GetModule("MUON");
for (Int_t i=0; i<10; i++) {
- AliMUONChamber* iChamber= &(MUON->Chamber(i));
- AliMUONResponse* response = iChamber->ResponseModel();
- AliMUONSegmentation* seg1 = iChamber->SegmentationModel(1);
- AliMUONSegmentation* seg2 = iChamber->SegmentationModel(2);
-//
RecModel = new AliMUONClusterFinderVS();
- RecModel->SetNperMax(90);
- RecModel->SetClusterSize(100);
- RecModel->SetDeclusterFlag(0);
- RecModel->SetSegmentation(seg1,seg2);
- RecModel->SetResponse(response);
// RecModel->SetTracks(16,17);
// RecModel->SetTracks(266,267);
+ RecModel->SetGhostChi2Cut(10);
MUON->SetReconstructionModel(i,RecModel);
}
//
//
Int_t Nh=0;
Int_t Nh1=0;
- for (int nev=evNumber1; nev<= evNumber2; nev++) {
- Int_t nparticles = gAlice->GetEvent(nev);
- cout << "nev " << nev <<endl;
- cout << "nparticles " << nparticles <<endl;
- if (nev < evNumber1) continue;
- if (nparticles <= 0) return;
- Int_t nbytes = 0;
- TClonesArray *Particles = gAlice->Particles();
- TTree *TD = gAlice->TreeD();
- Int_t nent=TD->GetEntries();
- if (MUON) {
- MUON->FindClusters(nev,nent-2);
- } // end if MUON
- } // event loop
-
- file->Close();
+ gAlice->RunReco("MUON", evNumber1, evNumber2);
}