message: Updated initialization of D* fits (Grazia)
diff: @@ -232,6 +232,10 @@ Int_t DmesonsFlowEvShapeAnalysis(Int_t cutmeth, Int_t analysismeth) {
  Int_t nMassBins=histtofit->GetNbinsX();
  Double_t hmin=histtofit->GetBinLowEdge(2); // need wide range for <pt>
  Double_t hmax=histtofit->GetBinLowEdge(nMassBins-2); // need wide range for <pt>
  + if (partname.Contains("Dstar")) {
  + if (hmin < 0.140) hmin=0.140;
  + if (hmax > 0.175) hmax=0.175;
  + }
  AliHFMassFitterVAR* fitter=new AliHFMassFitterVAR(histtofit,hmin,hmax,1,typeb,types);
  if(useTemplD0Refl){
  Printf("USE TEMPLATE FOR AVERAGE Pt");
  @@ -1224,12 +1228,12 @@ void FillSignalGraph(TList *masslist,TGraphAsymmErrors **gSignal,TGraphAsymmErro
  if(partname.Contains("Dplus")){
  massD=TDatabasePDG::Instance()->GetParticle(411)->Mass();
  }
  - if(partname.Contains("Dstar")) {
  - massD=(TDatabasePDG::Instance()->GetParticle(413)->Mass() - TDatabasePDG::Instance()->GetParticle(421)->Mass());
  - }
  if(partname.Contains("Ds")) {
  massD=(TDatabasePDG::Instance()->GetParticle(431)->Mass());
  }
  + if(partname.Contains("Dstar")) {
  + massD=(TDatabasePDG::Instance()->GetParticle(413)->Mass() - TDatabasePDG::Instance()->GetParticle(421)->Mass());
  + }
  }
  else {cerr << "Error: Wrong TDirectoryFile name " << dirname << ". Exit." << endl;}
  @@ -1280,6 +1284,9 @@ void FillSignalGraph(TList *masslist,TGraphAsymmErrors **gSignal,TGraphAsymmErro
  }
  fitter->SetInitialGaussianMean(massD);
  fitter->SetInitialGaussianSigma(0.012);
  + if (partname.Contains("Dstar")) {
  + fitter->SetInitialGaussianSigma(0.0004);
  + }
  fitter->SetUseLikelihoodFit();
  Bool_t ok=fitter->MassFitter(kFALSE);
  Double_t sigmaforcounting=0;
  @@ -1328,6 +1335,7 @@ void FillSignalGraph(TList *masslist,TGraphAsymmErrors **gSignal,TGraphAsymmErro
  fitter->SetFixReflOverS(sOverRef,kTRUE);
  }
  fitter->SetInitialGaussianMean(massD);
  + if (partname.Contains("Dstar")) fitter->SetInitialGaussianSigma(0.0004);
  fitter->SetUseLikelihoodFit();
  Bool_t ok=fitter->MassFitter(kFALSE);
  if(ok){
  @@ -1352,6 +1360,7 @@ void FillSignalGraph(TList *masslist,TGraphAsymmErrors **gSignal,TGraphAsymmErro
  fitter2->SetFixReflOverS(sOverRef,kTRUE);
  }
  fitter2->SetInitialGaussianMean(massD);
  + if (partname.Contains("Dstar")) fitter2->SetInitialGaussianSigma(0.0004);
  fitter2->SetFixGaussianSigma(sigma);
  if(fixAlsoMass) fitter2->SetFixGaussianMean(massFromFit);
  fitter2->SetUseLikelihoodFit();
  