+Int_t GetClusCoor(TObject *ge, TClonesArray *ITSrec, Int_t mod, TH2F *h2, TH1F *h1, Bool_t verb){
+
+ AliITSgeom *geom = (AliITSgeom*)ge;
+ Int_t nrecp = ITSrec->GetEntries();
+ if(nrecp>0){
+ Float_t lc[3]; for(Int_t i=0; i<3; i++) lc[i]=0.;
+ Float_t gc[3]; for(Int_t i=0; i<3; i++) gc[i]=0.;
+ if(verb){
+ cout<<"-------------------------------------------------------"<<endl;
+ cout<<"Number of CLUSTERS for module "<<mod<<": "<<nrecp<<endl;
+ }
+ for(Int_t irec=0;irec<nrecp;irec++) {
+ AliITSRecPoint *recp = (AliITSRecPoint*)ITSrec->At(irec);
+ Double_t rot[9];
+ geom->GetRotMatrix(mod,rot);
+ Int_t lay,lad,det;
+ geom->GetModuleId(mod,lay,lad,det);
+ Float_t tx,ty,tz;
+ geom->GetTrans(lay,lad,det,tx,ty,tz);
+
+ Double_t alpha=TMath::ATan2(rot[1],rot[0])+TMath::Pi();
+ Double_t phi1=TMath::Pi()/2+alpha;
+ if(lay==1) phi1+=TMath::Pi();
+
+ Float_t cp=TMath::Cos(phi1), sp=TMath::Sin(phi1);
+ Float_t r=tx*cp+ty*sp;
+ gc[0]= r*cp - recp->GetY()*sp;
+ gc[1]= r*sp + recp->GetY()*cp;
+ gc[2]=recp->GetZ();
+
+ if(verb){
+ Float_t r=sqrt(gc[0]*gc[0]+gc[1]*gc[1]);
+ cout<<"Global coor. + radius "<<gc[0]<<" "<<gc[1]<<" "<<gc[2]<<" ";
+ cout<<r<<endl;
+ cout<<"Associated track "<<recp->GetLabel(0)<<endl;
+ }
+ h2->Fill(gc[0],gc[1]);
+ h1->Fill(gc[2]);
+
+ }
+ }
+ return nrecp;
+}