NuFCExperimentGenerator Class Reference

#include <NuFCExperimentGenerator.h>

List of all members.

Public Member Functions
	NuFCExperimentGenerator (const NuFCRunInfo &run, const NuMMParameters &trueparams)
std::auto_ptr< NuMatrix >	GenerateND (const NuSystematicShifts &shifts)
	Generate and return a systematically shifted near detector spectrum.
std::auto_ptr< NuMatrix >	GenerateFD (const NuSystematicShifts &shifts)
	Generate and return a systematically shifted far detector spectrum.
std::auto_ptr< NuMatrix >	GenerateFDSample (const NuFCSample &sample, const NuSystematicShifts &shifts)
	Generate an individual FD sample.
void	SetDelegate (NuFCDelegate *delegate)
	Sets the delegate to use for dumping histogram output.
Private Member Functions
void	GotHistogram (TString name, NuMatrix *matrix)
Private Attributes
const NuFCRunInfo &	fRunInfo
	A reference to the runinfo object we are using.
NuEvent	fNu
	A blank NuEvent to expand into.
const NuMMParameters &	fTrueParams
NuFCDelegate *	fDelegate

---

Detailed Description

Class to generate experiments

Definition at line 20 of file NuFCExperimentGenerator.h.

---

Constructor & Destructor Documentation

NuFCExperimentGenerator::NuFCExperimentGenerator	(	const NuFCRunInfo &	run,
		const NuMMParameters &	trueparams
	)

Constructor. Give the RunInfo for the run you want to generate experiments for, and the true parameters to generate on.

Definition at line 25 of file NuFCExperimentGenerator.cxx.

  : fRunInfo(run),
    fTrueParams(trueparams),
    fDelegate(0)
{

}

---

Member Function Documentation

std::auto_ptr< NuMatrix > NuFCExperimentGenerator::GenerateFD

(

const NuSystematicShifts &

shifts

)

Generate and return a systematically shifted far detector spectrum.

Definition at line 78 of file NuFCExperimentGenerator.cxx.

References NuMatrix::Copy(), NuFCRunInfo::Events(), fRunInfo, GenerateFDSample(), NuFCRunInfo::GetBlankFD(), and NuFCEventManager::Samples().

Referenced by NuFCRunner::GenerateExperiment().

{
  // For each sample:
  //  Poisson Fluctuate
  //  Select a random event, shift and rejection sample
  // Sum the sample histograms
  // Build the Far detector Spectrum
  
  // Create our end result FD spectrum
  auto_ptr<NuMatrix>  fdData(fRunInfo.GetBlankFD()->Copy());
  
  // Get the vector of samples
  vector<NuFCSample> samples = fRunInfo.Events().Samples();
  
  // For every sample,
  vector<NuFCSample>::const_iterator it;
  for (it = samples.begin(); it != samples.end(); ++it)
  {
    NuFCSample smp = *it;
    
    // Generate the histogram for this sample
    auto_ptr<NuMatrix> smpMatrix = GenerateFDSample(smp, shifts);
    
    fdData->Add(*smpMatrix.get());
  }
  
  return fdData;
}

std::auto_ptr< NuMatrix > NuFCExperimentGenerator::GenerateFDSample	(	const NuFCSample &	sample,
		const NuSystematicShifts &	shifts
	)

Generate an individual FD sample.

Definition at line 109 of file NuFCExperimentGenerator.cxx.

References NuMatrix::Copy(), NuFCEvent2::EnergyMC(), NuFCRunInfo::Events(), NuFCEvent2::Expand(), NuFCEventManager::FDEvents(), fRunInfo, NuFCRunInfo::GetBlankFD(), NuFCRunInfo::GetName(), NuFCRunInfo::GetSampleWeight(), NuFCRunInfo::GetTruePrediction(), GotHistogram(), NuMatrixSpectrum::Integral(), Msg::kWarning, MSG, NuEvent::rw, NuFCEventManager::SampleName(), and NuMatrixSpectrum::Spectrum().

Referenced by GenerateFD().

{
  // Generate the prediction and weighting histograms
  NuMatrixSpectrum prediction = fRunInfo.GetTruePrediction(sample);
  const NuMatrixSpectrum &weights = fRunInfo.GetSampleWeight(sample);
  
  // Now work out the fluctuated event count
  Double_t averageEvents = prediction.Integral();
  Int_t poissonEvents = gRandom->Poisson(averageEvents);
  // MSG("NuFCExperimentGenerator",Msg::kInfo)
  //   << "Prediction for " <<  fRunInfo.Events().SampleName(sample)
  //   << " has " << averageEvents << " events. "
  //   << "Fluctuated is " << poissonEvents << " events" << endl;
  
  // Prepare the variables for drawing and filling
  const NuEventVector<NuFCEvent2> &evts = fRunInfo.Events().FDEvents(sample);
  NuEvent nu;
  // Get a new, blank FD spectrum
  auto_ptr<NuMatrix>  fdSample(fRunInfo.GetBlankFD()->Copy());
  fdSample->SetFillRecoEnergy();
  Int_t evtCount = 0, evtDraws = 0;
  // Int_t rejectExtrap = 0, rejectRw = 0;
  
  // Pull out all the events from the distribution
  // for (Int_t evtN = 0; evtN < poissonEvents; ++evtN)
  // Keep pulling until we get the events we need
  while (evtCount < poissonEvents)
  {
    Double_t rand = gRandom->Uniform();
    
    // Pick a random event from the sample
    Int_t evtIndex = gRandom->Integer(evts.size());
    const NuFCEvent2 &evt = evts[evtIndex];
    ++evtDraws;
    // Get the weight from the truth energy ratio
    Int_t bin = weights.Spectrum()->GetXaxis()->FindFixBin(evt.EnergyMC());
    Double_t evtWeight = weights.Spectrum()->GetBinContent(bin);
  
    // Add in the events intrinsic reweight to this
    // Is this necessary? In principle, It should not affect the final energy
    // spectrum because shape is determined by the rejection sampling..
    // What it should do, however, is ensure that the distribution of points within
    // a single true energy bin is correct. This might not have any effect, as in
    // priciple all events in a single bin should have the same rw value (they don't)
    // and this in turn should probably not have any effect on the final answer (as that
    // woudl probably indicate some realiance on pt/pz for final reconstructed values, 
    // which doesn't make any sense).
    // UPDATE:
    // This is disabled for now because whilst the rw->reco distributions are independant,
    // the shape of the rw distribution in an individual bin compared to another bin can
    // affect the probability of rejection for a single bin and therefore change the shape
    // of the end distribution by changing a bin-dependent normalisation. The way around
    // this is to fluctuate each bin on it's own and control the normalisation that way,
    // but that is a bit of a severe change for an effect which we expect to make little,
    // if any difference.
    // evtWeight *= evt.Rw();
  
    // Store this pre-systematic weight
    // Double_t preSystematicsWeight = evtWeight;
  
    // TODO: Systematics code
    // apply the systematic weight here
    // Adjust the event count number if the systematics have affected the event count
    
    // Warn if we are getting unrejectable probabilities
    if (evtWeight > 1.0) {
      MSG("NuFCExperimentGenerator",Msg::kWarning)
        << "Got rejection probability of " << 1-evtWeight
        << "; impossible to reject this event" << endl;
    }
    // Does this event fail the rejection test?
    if (rand > evtWeight) continue;
    
    // We have passed the event - add to our histogram!
    evt.Expand(nu);
    nu.rw = 1.0;
    fdSample->Fill(nu);
    ++evtCount;
  }
  GotHistogram("FD_"+fRunInfo.GetName()+"_"+fRunInfo.Events().SampleName(sample), fdSample.get());
  // MSG("NuFCExperimentGenerator",Msg::kInfo)
  //   << "Generated " << poissonEvents << " events" << endl;
  // MSG("NuFCExperimentGenerator",Msg::kInfo)
  //   << "Drew a total of " << evtDraws << " events." << endl;
  
  return fdSample;
}

std::auto_ptr< NuMatrix > NuFCExperimentGenerator::GenerateND

(

const NuSystematicShifts &

shifts

)

Generate and return a systematically shifted near detector spectrum.

Definition at line 45 of file NuFCExperimentGenerator.cxx.

References NuFCRunInfo::Events(), fNu, fRunInfo, NuMMRun::GetNDData(), NuFCRunInfo::GetRun(), NuFCEventManager::NDEvents(), and NuEventVector< T >::pot().

Referenced by NuFCRunner::GenerateExperiment().

{
  // Get a blank copy of the near detector spectrum
  assert(fRunInfo.GetRun().GetNDData().get());
  auto_ptr<NuMatrix> nd(fRunInfo.GetRun().GetNDData()->Copy());
  nd->SetFillRecoEnergy();
  nd->Scale(0);
  
  // Now, get the list of events in the run
  const NuEventVector<NuFCEvent2> &events = fRunInfo.Events().NDEvents();
  
  // Set the new ND POT (this should not change, but do it here anwyay)
  nd->ResetPOT(events.pot());
  
  // Copy these events into the new ND spectrum
  NuEventVector<NuFCEvent2>::const_iterator it;
  for (it = events.begin(); it != events.end(); ++it)
  {
    (*it).Expand(fNu);
    
    // Each experiment should be systematically shifted here
    // NuSystematic syst;
    // syst.SetShiftsAsValues(shifts);
    // sys.Shift(fNu);
    
    nd->Fill(fNu);
  }
  return nd;
}

void NuFCExperimentGenerator::GotHistogram	(	TString	name,
		NuMatrix *	matrix
	)			[private]

Definition at line 36 of file NuFCExperimentGenerator.cxx.

References fDelegate, and NuFCDelegate::GotHistogram().

Referenced by GenerateFDSample().

{
  if (!fDelegate) return;
  
  fDelegate->GotHistogram(name, matrix);
}

void NuFCExperimentGenerator::SetDelegate

(

NuFCDelegate *

delegate

)

[inline]

Sets the delegate to use for dumping histogram output.

Definition at line 58 of file NuFCExperimentGenerator.h.

References fDelegate.

Referenced by NuFCRunner::GenerateExperiment().

{ fDelegate = delegate; }

---

Member Data Documentation

NuFCDelegate* NuFCExperimentGenerator::fDelegate [private]

Definition at line 34 of file NuFCExperimentGenerator.h.

Referenced by GotHistogram(), and SetDelegate().

NuEvent NuFCExperimentGenerator::fNu [private]

A blank NuEvent to expand into.

Definition at line 26 of file NuFCExperimentGenerator.h.

Referenced by GenerateND().

const NuFCRunInfo& NuFCExperimentGenerator::fRunInfo [private]

A reference to the runinfo object we are using.

Definition at line 23 of file NuFCExperimentGenerator.h.

Referenced by GenerateFD(), GenerateFDSample(), and GenerateND().

const NuMMParameters& NuFCExperimentGenerator::fTrueParams [private]

Definition at line 29 of file NuFCExperimentGenerator.h.

---

The documentation for this class was generated from the following files:

• NuFCExperimentGenerator.h
• NuFCExperimentGenerator.cxx

---