// Get the particle masses
// parent
nbins = 1000;
-
+ mass_min = 0.;
+ mass_max = 0.;
pionmass = (TDatabasePDG::Instance()->GetParticle(111))->Mass();
etamass = (TDatabasePDG::Instance()->GetParticle(221))->Mass();
omegamass = (TDatabasePDG::Instance()->GetParticle(223))->Mass();
Double_t pimass = 0.13956995;
+ if(mass>pimass){
corr = vwidth*(vmass/mass)*exp(1.5*log((mass*mass/4.0-pimass*pimass)
/(vmass*vmass/4.0-pimass*pimass)));
+ }
epsilon = (emass/mass)*(emass/mass);