1 Int_t AliITSFindClusters() {
3 printf("FindClusters\n");
5 TFile *in=TFile::Open("galice.root","UPDATE");
6 if (!in->IsOpen()) {cerr<<"Can't open galice.root !\n"; return 2;}
10 if (!(gAlice=(AliRun*)in->Get("gAlice"))) {
11 cerr<<"gAlice have not been found on galice.root !\n";
18 AliITS *ITS = (AliITS*)gAlice->GetDetector("ITS");
20 cerr<<"ITSFindClusters.C : AliITS object not found on file\n";
23 Int_t ver = ITS->IsVersion();
24 cerr<<"ITS version "<<ver<<" has been found !\n";
27 // Set the models for cluster finding
28 AliITSgeom *geom = ITS->GetITSgeom();
31 AliITSDetType *iDetType=ITS->DetType(0);
32 AliITSsegmentationSPD *seg0=(AliITSsegmentationSPD*)iDetType->GetSegmentationModel();
33 TClonesArray *dig0 = ITS->DigitsAddress(0);
34 TClonesArray *recp0 = ITS->ClustersAddress(0);
35 AliITSClusterFinderSPD *rec0=new AliITSClusterFinderSPD(seg0,dig0,recp0);
36 ITS->SetReconstructionModel(0,rec0);
38 printf("SPD dimensions %f %f \n",seg0->Dx(),seg0->Dz());
39 printf("SPD npixels %d %d \n",seg0->Npz(),seg0->Npx());
43 Float_t baseline = 10.;
45 Float_t thres = baseline+3.*noise;
47 AliITSDetType *iDetType=ITS->DetType(1);
48 AliITSsegmentationSDD *seg1=(AliITSsegmentationSDD*)iDetType->GetSegmentationModel();
49 if (!seg1) seg1 = new AliITSsegmentationSDD(geom);
50 AliITSresponseSDD *res1 = (AliITSresponseSDD*)iDetType->GetResponseModel();
51 if (!res1) res1=new AliITSresponseSDD();
52 res1->SetMagicValue(900.);
53 Float_t magic = res1->MagicValue();
54 Float_t top = res1->MaxAdc();
56 res1->SetNoiseParam(noise,baseline);
57 TClonesArray *dig1 = ITS->DigitsAddress(1);
58 TClonesArray *recp1 = ITS->ClustersAddress(1);
59 AliITSClusterFinderSDD *rec1=new AliITSClusterFinderSDD(seg1,res1,dig1,recp1);
60 rec1->SetMinNCells(6);
61 rec1->SetTimeCorr(70.);
62 rec1->SetCutAmplitude((int)thres);
63 ITS->SetReconstructionModel(1,rec1);
67 AliITSDetType *iDetType=ITS->DetType(2);
68 AliITSsegmentationSSD *seg2=(AliITSsegmentationSSD*)iDetType->GetSegmentationModel();
69 TClonesArray *dig2 = ITS->DigitsAddress(2);
70 //TClonesArray *recp2 = ITS->ClustersAddress(2);
71 AliITSClusterFinderSSD *rec2=new AliITSClusterFinderSSD(seg2,dig2);
72 ITS->SetReconstructionModel(2,rec2);
74 printf("SSD dimensions %f %f \n",seg2->Dx(),seg2->Dz());
75 printf("SSD nstrips %d %d \n",seg2->Npz(),seg2->Npx());
79 if(!gAlice->TreeR()) gAlice->MakeTree("R");
87 cerr<<"Looking for clusters...\n";
90 ITS->DigitsToRecPoints(0,1,"All");
94 cerr<<"Invalid ITS version !\n";
98 timer.Stop(); timer.Print();
105 delete gAlice; gAlice=0;