added more efficient initialization in AbsDecay::fillWignerDs

PwaUtils/AbsDecay.cc

@@ -304,10 +304,14 @@ void AbsDecay::fillWignerDs(std::map<std::string, Vector4<double> >& fsMap, Vect
     for (Spin lam12=-lam12Max; lam12<=lam12Max; ++lam12){
       double thePhi=0.;
       if(_hasMotherPart) thePhi=daughter2HelMother.Phi();
-      Id3StringType IdSpinMotherLamMotherLam12=FunctionUtils::spin3Index(spinMother, lamMother, lam12); 
+      Id3StringType IdSpinMotherLamMotherLam12=FunctionUtils::spin3Index(spinMother, lamMother, lam12);
+
+      std::map<Id3StringType, complex<double> >::iterator found = evtData->WignerDStringId[_wignerDKey].find(IdSpinMotherLamMotherLam12);
+      if(found != evtData->WignerDStringId[_wignerDKey].end()){
+        continue;
+      } 
       evtData->WignerDStringId[_wignerDKey][IdSpinMotherLamMotherLam12]=Wigner_D(thePhi,daughter2HelMother.Theta(),0,spinMother,lamMother,lam12);
-      //      evtData->WignerDsString[_wignerDKey][spinMother][lamMother][lam12]=Wigner_D(thePhi,daughter2HelMother.Theta(),0,spinMother,lamMother,lam12);
-      //      if(evtData->WignerDsString[_wignerDKey][spinMother][lamMother][lam12] != evtData->WignerDsString[_wignerDKey][spinMother][lamMother][lam12]){
+
       if(evtData->WignerDStringId[_wignerDKey][IdSpinMotherLamMotherLam12] != evtData->WignerDStringId[_wignerDKey][IdSpinMotherLamMotherLam12]){
 	DebugMsg << "WignerDsString nan:\t" << evtData->WignerDsString[_wignerDKey][spinMother][lamMother][lam12] << endmsg;
     	DebugMsg << "daughter2HelMother.Theta():\t" << daughter2HelMother.Theta() << endmsg;
@@ -328,7 +332,6 @@ void AbsDecay::fillWignerDs(std::map<std::string, Vector4<double> >& fsMap, Vect
      }
      if (_useProdBarrier){
        double qVal=daughter2HelMother.P();
-       //       double qValNorm=breakupMomQ(mother4Vec.M(), daughter1_fsMassSum, daughter2_fsMassSum).real();
        double qValNorm=breakupMomQ(mother4Vec.M(), massSumFsParticlesDec1(), massSumFsParticlesDec2()).real();
        evtData->DoubleString[_wignerDKey]=qVal;
        evtData->DoubleString[_wignerDKey+"qNorm"] = qValNorm;