-// Double_t efit ;
-
-// AliEMCALDigit * digit ;
-// Int_t iDigit ;
-
-// for( iDigit = 0 ; iDigit < nOdigits ; iDigit++) {
-
-// digit = (AliEMCALDigit*) digits->At( emcDigits[iDigit] ) ;
-
-// Int_t relid[4] ;
-// Float_t xDigit ;
-// Float_t zDigit ;
-
-// geom->AbsToRelNumbering(digit->GetId(), relid) ;
-
-// geom->RelPosInModule(relid, xDigit, zDigit) ;
-
-// if(iflag == 2){ // calculate gradient
-// Int_t iParam = 0 ;
-// efit = 0 ;
-// while(iParam < nPar ){
-// Double_t distance = (xDigit - x[iParam]) * (xDigit - x[iParam]) ;
-// iParam++ ;
-// distance += (zDigit - x[iParam]) * (zDigit - x[iParam]) ;
-// distance = TMath::Sqrt( distance ) ;
-// iParam++ ;
-// efit += x[iParam] * ShowerShape(distance) ;
-// iParam++ ;
-// }
-// Double_t sum = 2. * (efit - emcEnergies[iDigit]) / emcEnergies[iDigit] ; // Here we assume, that sigma = sqrt(E)
-// iParam = 0 ;
-// while(iParam < nPar ){
-// Double_t xpar = x[iParam] ;
-// Double_t zpar = x[iParam+1] ;
-// Double_t epar = x[iParam+2] ;
-// Double_t dr = TMath::Sqrt( (xDigit - xpar) * (xDigit - xpar) + (zDigit - zpar) * (zDigit - zpar) );
-// Double_t shape = sum * ShowerShape(dr) ;
-// Double_t r4 = dr*dr*dr*dr ;
-// Double_t r295 = TMath::Power(dr,2.95) ;
-// Double_t deriv =-4. * dr*dr * ( 2.32 / ( (2.32 + 0.26 * r4) * (2.32 + 0.26 * r4) ) +
-// 0.0316 * (1. + 0.0171 * r295) / ( ( 1. + 0.0652 * r295) * (1. + 0.0652 * r295) ) ) ;
-
-// Grad[iParam] += epar * shape * deriv * (xpar - xDigit) ; // Derivative over x
-// iParam++ ;
-// Grad[iParam] += epar * shape * deriv * (zpar - zDigit) ; // Derivative over z
-// iParam++ ;
-// Grad[iParam] += shape ; // Derivative over energy
-// iParam++ ;
-// }
-// }
-// efit = 0;
-// iparam = 0 ;
-
-// while(iparam < nPar ){
-// Double_t xpar = x[iparam] ;
-// Double_t zpar = x[iparam+1] ;
-// Double_t epar = x[iparam+2] ;
-// iparam += 3 ;
-// Double_t distance = (xDigit - xpar) * (xDigit - xpar) + (zDigit - zpar) * (zDigit - zpar) ;
-// distance = TMath::Sqrt(distance) ;
-// efit += epar * ShowerShape(distance) ;
-// }
-
-// fret += (efit-emcEnergies[iDigit])*(efit-emcEnergies[iDigit])/emcEnergies[iDigit] ;
-// // Here we assume, that sigma = sqrt(E)
-// }
-
-}
-
-//____________________________________________________________________________
-void AliEMCALClusterizerv1::Print(Option_t * option)const
-{
- // Print clusterizer parameters
-
- if( strcmp(GetName(), "") !=0 ){
-
- // Print parameters
-
- TString taskName(GetName()) ;
- taskName.ReplaceAll(Version(), "") ;
-
- cout << "---------------"<< taskName.Data() << " " << GetTitle()<< "-----------" << endl
- << "Clusterizing digits from the file: " << fHeaderFileName.Data() << endl
- << " Branch: " << fDigitsBranchTitle.Data() << endl
- << endl
- << " EMC Clustering threshold = " << fTowerClusteringThreshold << endl
- << " EMC Local Maximum cut = " << fTowerLocMaxCut << endl
- << " EMC Logarothmic weight = " << fW0 << endl
- << endl
- << " CPV Clustering threshold = " << fPreShoClusteringThreshold << endl
- << " CPV Local Maximum cut = " << fPreShoLocMaxCut << endl
- << " CPV Logarothmic weight = " << fW0CPV << endl
- << endl ;
- if(fToUnfold)
- cout << " Unfolding on " << endl ;
- else
- cout << " Unfolding off " << endl ;
-
- cout << "------------------------------------------------------------------" <<endl ;
- }
- else
- cout << " AliEMCALClusterizerv1 not initialized " << endl ;