private void OnCLEMAnimalBreeding(object sender, EventArgs e)
{
this.Status = ActivityStatus.NotNeeded;
NumberConceived = 0;
// get list of all pregnant females
List <RuminantFemale> pregnantherd = CurrentHerd(true).OfType <RuminantFemale>().Where(a => a.IsPregnant).ToList();
// determine all fetus and newborn mortality of all pregnant females.
foreach (RuminantFemale female in pregnantherd)
{
// calculate fetus and newborn mortality
// total mortality / (gestation months + 1) to get monthly mortality
// done here before births to account for post birth motality as well..
// IsPregnant status does not change until births occur in next section so will include mortality in month of birth
// needs to be calculated for each offspring carried.
for (int i = 0; i < female.CarryingCount; i++)
{
var rnd = RandomNumberGenerator.Generator.NextDouble();
if (rnd < (female.BreedParams.PrenatalMortality / (female.BreedParams.GestationLength + 1)))
{
female.OneOffspringDies();
if (female.NumberOfOffspring == 0)
{
// report conception status changed when last multiple birth dies.
female.BreedParams.OnConceptionStatusChanged(new Reporting.ConceptionStatusChangedEventArgs(Reporting.ConceptionStatus.Failed, female, clock.Today));
}
}
}
if (female.BirthDue)
{
int numberOfNewborn = female.CarryingCount;
for (int i = 0; i < numberOfNewborn; i++)
{
bool isMale = RandomNumberGenerator.Generator.NextDouble() <= female.BreedParams.ProportionOffspringMale;
Sex sex = isMale ? Sex.Male : Sex.Female;
double weight = female.BreedParams.SRWBirth * female.StandardReferenceWeight * (1 - 0.33 * (1 - female.Weight / female.StandardReferenceWeight));
Ruminant newSucklingRuminant = Ruminant.Create(sex, female.BreedParams, 0, weight);
newSucklingRuminant.HerdName = female.HerdName;
newSucklingRuminant.Breed = female.BreedParams.Breed;
newSucklingRuminant.ID = HerdResource.NextUniqueID;
newSucklingRuminant.Location = female.Location;
newSucklingRuminant.Mother = female;
newSucklingRuminant.Number = 1;
newSucklingRuminant.SetUnweaned();
// suckling/calf weight from Freer
newSucklingRuminant.PreviousWeight = newSucklingRuminant.Weight;
newSucklingRuminant.SaleFlag = HerdChangeReason.Born;
// add attributes inherited from mother
foreach (var attribute in female.Attributes.Items)
{
if (attribute.Value != null)
{
newSucklingRuminant.Attributes.Add(attribute.Key, attribute.Value.GetInheritedAttribute() as IIndividualAttribute);
}
}
HerdResource.AddRuminant(newSucklingRuminant, this);
// add to sucklings
female.SucklingOffspringList.Add(newSucklingRuminant);
// this now reports for each individual born not a birth event as individual wean events are reported
female.BreedParams.OnConceptionStatusChanged(new Reporting.ConceptionStatusChangedEventArgs(Reporting.ConceptionStatus.Birth, female, clock.Today));
}
female.UpdateBirthDetails();
this.Status = ActivityStatus.Success;
}
}
// Perform breeding
IEnumerable <Ruminant> herd = null;
if (useControlledMating && controlledMating.TimingOK)
{
// determined by controlled mating and subsequent timer (e.g. smart milking)
herd = controlledMating.BreedersToMate();
}
else if (!useControlledMating && TimingOK)
{
// whole herd for activity including males
herd = CurrentHerd(true);
}
if (herd != null && herd.Any())
{
// group by location
var breeders = from ind in herd
where ind.IsAbleToBreed
group ind by ind.Location into grp
select grp;
// identify not ready for reporting and tracking
var notReadyBreeders = herd.Where(a => a.Sex == Sex.Female).Cast <RuminantFemale>().Where(a => a.IsBreeder && !a.IsAbleToBreed && !a.IsPregnant);
foreach (RuminantFemale female in notReadyBreeders)
{
female.BreedParams.OnConceptionStatusChanged(new Reporting.ConceptionStatusChangedEventArgs(Reporting.ConceptionStatus.NotReady, female, clock.Today));
}
int numberPossible = breeders.Sum(a => a.Count());
int numberServiced = 1;
List <Ruminant> maleBreeders = new List <Ruminant>();
// for each location where parts of this herd are located
foreach (var location in breeders)
{
numberPossible = -1;
if (useControlledMating)
{
numberPossible = Convert.ToInt32(location.OfType <RuminantFemale>().Count(), CultureInfo.InvariantCulture);
}
else
{
numberPossible = 0;
// uncontrolled conception
if (location.GroupBy(a => a.Sex).Count() == 2)
{
int maleCount = location.OfType <RuminantMale>().Count();
// get a list of males to provide attributes when incontrolled mating.
if (maleCount > 0 && location.FirstOrDefault().BreedParams.IncludedAttributeInheritanceWhenMating)
{
maleBreeders = location.Where(a => a.Sex == Sex.Male).ToList();
}
int femaleCount = location.Where(a => a.Sex == Sex.Female).Count();
numberPossible = Convert.ToInt32(Math.Ceiling(maleCount * location.FirstOrDefault().BreedParams.MaximumMaleMatingsPerDay * 30), CultureInfo.InvariantCulture);
}
}
numberServiced = 0;
lastJoinIndex = -1;
int cnt = 0;
// shuffle the not pregnant females when obtained to avoid any inherant order by creation of individuals affecting which individuals are available first
var notPregnantFemales = location.OfType <RuminantFemale>().Where(a => !a.IsPregnant).OrderBy(a => RandomNumberGenerator.Generator.Next()).ToList();
int totalToBreed = notPregnantFemales.Count;
while (cnt < totalToBreed)
{
RuminantFemale female = notPregnantFemales.ElementAt(cnt);
Reporting.ConceptionStatus status = Reporting.ConceptionStatus.NotMated;
if (numberServiced < numberPossible)
{
// calculate conception
double conceptionRate = ConceptionRate(female, out status);
// if mandatory attributes are present in the herd, save male value with female details.
// update male for both successful and failed matings (next if statement
if (female.BreedParams.IncludedAttributeInheritanceWhenMating)
{
object male = null;
if (useControlledMating)
{
bool newJoining = needsNewJoiningMale(controlledMating.JoiningsPerMale, numberServiced);
// save all male attributes
AddMalesAttributeDetails(female, controlledMating.SireAttributes, newJoining);
}
else
{
male = maleBreeders[RandomNumberGenerator.Generator.Next(0, maleBreeders.Count() - 1)];
female.LastMatingStyle = ((male as RuminantMale).IsWildBreeder ? MatingStyle.WildBreeder : MatingStyle.Natural);
// randomly select male
AddMalesAttributeDetails(female, male as Ruminant);
}
}
if (conceptionRate > 0)
{
if (RandomNumberGenerator.Generator.NextDouble() <= conceptionRate)
{
female.UpdateConceptionDetails(female.CalulateNumberOfOffspringThisPregnancy(), conceptionRate, 0);
if (useControlledMating)
{
female.LastMatingStyle = MatingStyle.Controlled;
}
status = Reporting.ConceptionStatus.Conceived;
NumberConceived++;
}
else
{
status = Reporting.ConceptionStatus.Unsuccessful;
}
}
numberServiced++;
this.Status = ActivityStatus.Success;
}
// report change in breeding status
// do not report for -1 (controlled mating outside timing)
if (numberPossible >= 0 && status != Reporting.ConceptionStatus.NotAvailable)
{
female.BreedParams.OnConceptionStatusChanged(new Reporting.ConceptionStatusChangedEventArgs(status, female, clock.Today));
}
cnt++;
}
// report a natural mating locations for transparency via a message
if (numberServiced > 0 & !useControlledMating)
{
string warning = $"Natural (uncontrolled) mating ocurred in [r={(location.Key ?? "Not specified - general yards")}]";
Warnings.CheckAndWrite(warning, Summary, this, MessageType.Information);
}
}
}
// report that this activity was performed as it does not use base GetResourcesRequired
this.TriggerOnActivityPerformed();
}
private void OnCLEMAnimalBreeding(object sender, EventArgs e)
{
this.Status = ActivityStatus.NotNeeded;
NumberConceived = 0;
// get list of all pregnant females
List <RuminantFemale> pregnantherd = CurrentHerd(true).Where(a => a.Gender == Sex.Female).Cast <RuminantFemale>().Where(a => a.IsPregnant).ToList();
// determine all fetus and newborn mortality of all pregnant females.
foreach (RuminantFemale female in pregnantherd)
{
// calculate fetus and newborn mortality
// total mortality / (gestation months + 1) to get monthly mortality
// done here before births to account for post birth motality as well..
// IsPregnant status does not change until births occur in next section so will include mortality in month of birth
// needs to be calculated for each offspring carried.
for (int i = 0; i < female.CarryingCount; i++)
{
if (RandomNumberGenerator.Generator.NextDouble() < (female.BreedParams.PrenatalMortality / (female.BreedParams.GestationLength + 1)))
{
female.OneOffspringDies();
if (female.NumberOfOffspring == 0)
{
// report conception status changed when last multiple birth dies.
female.BreedParams.OnConceptionStatusChanged(new Reporting.ConceptionStatusChangedEventArgs(Reporting.ConceptionStatus.Failed, female, Clock.Today));
}
}
}
if (female.BirthDue)
{
int numberOfNewborn = female.CarryingCount;
for (int i = 0; i < numberOfNewborn; i++)
{
object newCalf = null;
bool isMale = (RandomNumberGenerator.Generator.NextDouble() <= female.BreedParams.ProportionOffspringMale);
double weight = female.BreedParams.SRWBirth * female.StandardReferenceWeight * (1 - 0.33 * (1 - female.Weight / female.StandardReferenceWeight));
if (isMale)
{
newCalf = new RuminantMale(0, Sex.Male, weight, female.BreedParams);
}
else
{
newCalf = new RuminantFemale(0, Sex.Female, weight, female.BreedParams);
}
Ruminant newCalfRuminant = newCalf as Ruminant;
newCalfRuminant.HerdName = female.HerdName;
newCalfRuminant.Breed = female.BreedParams.Breed;
newCalfRuminant.ID = Resources.RuminantHerd().NextUniqueID;
newCalfRuminant.Location = female.Location;
newCalfRuminant.Mother = female;
newCalfRuminant.Number = 1;
newCalfRuminant.SetUnweaned();
// calf weight from Freer
newCalfRuminant.PreviousWeight = newCalfRuminant.Weight;
newCalfRuminant.SaleFlag = HerdChangeReason.Born;
// add attributes inherited from mother
foreach (var attribute in female.Attributes)
{
newCalfRuminant.AddAttribute(attribute.Key, attribute.Value.GetInheritedAttribute() as ICLEMAttribute);
}
Resources.RuminantHerd().AddRuminant(newCalfRuminant, this);
// add to sucklings
female.SucklingOffspringList.Add(newCalfRuminant);
// this now reports for each individual born not a birth event as individual wean events are reported
female.BreedParams.OnConceptionStatusChanged(new Reporting.ConceptionStatusChangedEventArgs(Reporting.ConceptionStatus.Birth, female, Clock.Today));
}
female.UpdateBirthDetails();
this.Status = ActivityStatus.Success;
}
}
// Perform breeding
IEnumerable <Ruminant> herd = null;
if (useControlledMating && controlledMating.TimingOK)
{
// determined by controlled mating and subsequent timer (e.g. smart milking)
herd = controlledMating.BreedersToMate();
this.TriggerOnActivityPerformed();
}
else if (!useControlledMating && TimingOK)
{
// whole herd for activity
herd = CurrentHerd(true);
// report that this activity was performed as it does not use base GetResourcesRequired
this.TriggerOnActivityPerformed();
}
if (herd != null && herd.Count() > 0)
{
// group by location
var breeders = from ind in herd
where ind.IsAbleToBreed
group ind by ind.Location into grp
select grp;
int breedersCount = breeders.Count();
int numberPossible = breedersCount;
int numberServiced = 1;
List <Ruminant> maleBreeders = new List <Ruminant>();
// for each location where parts of this herd are located
foreach (var location in breeders)
{
numberPossible = -1;
if (useControlledMating)
{
numberPossible = Convert.ToInt32(location.Where(a => a.Gender == Sex.Female).Count(), CultureInfo.InvariantCulture);
}
else
{
numberPossible = 0;
// uncontrolled conception
if (location.GroupBy(a => a.Gender).Count() == 2)
{
int maleCount = location.Where(a => a.Gender == Sex.Male).Count();
// get a list of males to provide attributes when incontrolled mating.
if (maleCount > 0 && location.FirstOrDefault().BreedParams.IncludedAttributeInheritanceWhenMating)
{
maleBreeders = location.Where(a => a.Gender == Sex.Male).ToList();
}
int femaleCount = location.Where(a => a.Gender == Sex.Female).Count();
numberPossible = Convert.ToInt32(Math.Ceiling(maleCount * location.FirstOrDefault().BreedParams.MaximumMaleMatingsPerDay * 30), CultureInfo.InvariantCulture);
}
}
numberServiced = 1;
foreach (RuminantFemale female in location.Where(a => a.Gender == Sex.Female).Cast <RuminantFemale>().Where(a => !a.IsPregnant & a.Age <= a.BreedParams.MaximumAgeMating).ToList())
{
Reporting.ConceptionStatus status = Reporting.ConceptionStatus.NotMated;
if (numberServiced <= numberPossible)
{
// calculate conception
double conceptionRate = ConceptionRate(female, out status);
if (conceptionRate > 0)
{
if (RandomNumberGenerator.Generator.NextDouble() <= conceptionRate)
{
female.UpdateConceptionDetails(female.CalulateNumberOfOffspringThisPregnancy(), conceptionRate, 0);
// if mandatory attributes are present in the herd, save male value with female details.
if (female.BreedParams.IncludedAttributeInheritanceWhenMating)
{
if (useControlledMating)
{
// save all male attributes
AddMalesAttributeDetails(female, controlledMating.SireAttributes);
}
else
{
// randomly select male
AddMalesAttributeDetails(female, maleBreeders[RandomNumberGenerator.Generator.Next(0, maleBreeders.Count() - 1)]);
}
}
status = Reporting.ConceptionStatus.Conceived;
NumberConceived++;
}
}
numberServiced++;
this.Status = ActivityStatus.Success;
}
// report change in breeding status
// do not report for -1 (controlled mating outside timing)
if (numberPossible >= 0 && status != Reporting.ConceptionStatus.NotAvailable)
{
female.BreedParams.OnConceptionStatusChanged(new Reporting.ConceptionStatusChangedEventArgs(status, female, Clock.Today));
}
}
// report a natural mating locations for transparency via a message
if (numberServiced > 0 & !useControlledMating)
{
string warning = "Natural (uncontrolled) mating ocurred in [r=" + location.Key + "]";
if (!Warnings.Exists(warning))
{
Warnings.Add(warning);
Summary.WriteMessage(this, warning);
}
}
}
}
}