-//
-// Generates pad hits (simulated cluster)
-// using the segmentation and the response model
-
- Float_t dx, dy;
- Float_t local[3];
- //Float_t source[3];
- Float_t global[3];
- //
- // Width of the integration area
- //
- dx=(fResponse->SigmaIntegration())*(fResponse->ChargeSpreadX());
- dy=(fResponse->SigmaIntegration())*(fResponse->ChargeSpreadY());
- //
- // Get pulse height from energy loss and generate feedback photons
- Float_t qtot=0;
-
- local[0]=xhit;
- // z-position of the wires relative to the RICH mother volume
- // (2 mmm before CsI) old value: 6.076
- local[1]=1.276 + fGeometry->GetGapThickness()/2 + fGeometry->GetProximityGapThickness() - .39;
- //printf("AliRICHChamber feedback origin:%f",local[1]);
- local[2]=yhit;
-
- LocaltoGlobal(local,global);
-
- Int_t Nfp=0;
-
- if (res==mip) {
- qtot = fResponse->IntPH(eloss);
- Nfp = fResponse->FeedBackPhotons(global,qtot);
- } else if (res==cerenkov) {
- qtot = fResponse->IntPH();
- Nfp = fResponse->FeedBackPhotons(global,qtot);
- }
-
- //printf("Feedbacks:%d\n",Nfp);
-
- //
- // Loop Over Pads
-
- Float_t qcheck=0, qp=0;
-
- nnew=0;
- for (Int_t i=1; i<=fnsec; i++) {
- qcheck=0;
- for (fSegmentation->FirstPad(xhit, yhit, dx, dy);
- fSegmentation->MorePads();
- fSegmentation->NextPad())
- {
- qp= fResponse->IntXY(fSegmentation);
- qp= TMath::Abs(qp);
-
- //printf("Qp:%f\n",qp);
-
- if (qp > 1.e-4) {
- qcheck+=qp;
- //
- // --- store signal information
- newclust[0][nnew]=qtot;
- newclust[1][nnew]=fSegmentation->Ix();
- newclust[2][nnew]=fSegmentation->Iy();
- newclust[3][nnew]=qp * qtot;
- newclust[4][nnew]=fSegmentation->ISector();
- newclust[5][nnew]=(Float_t) i;
- nnew++;
- //printf("Newcluster:%d\n",i);
- }
- } // Pad loop
- } // Cathode plane loop
- //if (fSegmentation->ISector()==2)
- //printf("Nnew:%d\n\n\n\n",nnew);
-}
-
-
-
-
+//debug printout method
+ printf("%s r=%8.3f theta=%5.1f phi=%5.1f x=%8.3f y=%8.3f z=%8.3f %6.2f,%6.2f %6.2f,%6.2f %6.2f,%6.2f\n",fName.Data(),
+ Rho(), ThetaD(),PhiD(), X(), Y(), Z(),
+ ThetaXd(),PhiXd(),ThetaYd(),PhiYd(),ThetaZd(),PhiZd());
+}//Print()
+//__________________________________________________________________________________________________