Blinder.h

Go to the documentation of this file.

// $Id: Blinder.h,v 1.4 2005/06/20 15:24:32 tagg Exp $
//
// Module to apply blindness for physics analysis.
//
// Blinder::Ana returns kPassed for events that should be in the open dataset
// and kFailed for events that should be only in the blind dataset.
//
// It returns kPassed for all events in the ND, and uses
// the algorithm presented at the Jan 2005 collab meeting to
// select events from the FD.
//
// Nathaniel Tagg
// n.tagg1@physics.ox.ac.uk
#ifndef BLINDER_H
#define BLINDER_H

#include "JobControl/JobCModule.h"

class RawDigitDataBlock;
class TFile;
class TNtuple;
class REROOT_NeuKin;
class TRandom;

class Blinder : public JobCModule
{
public:
  Blinder();
  ~Blinder();

public:
  // Analysis and Reconstruction methods
  virtual JobCResult Ana(const MomNavigator* mom);
  virtual JobCResult Reco(MomNavigator* mom);

  // Module configuration
  const Registry& DefaultConfig() const;
  void            Config(const Registry& r);


private:
  Bool_t GetInfo(const MomNavigator* mom,
                 Int_t &run,
                 Int_t &snarl,
                 Int_t &epoch,
                 Float_t &length,
                 Float_t &energy);

  const REROOT_NeuKin* GetTruth(const MomNavigator* mom);

  // Config data:
  Int_t fSeed;

  // The actual work is done by:

  class BlindingFunction {
  public:
    BlindingFunction();
    void    SetParameters( TRandom* paramChooser, Int_t seed );
    Float_t EventBlindProb( Float_t eventLength, Float_t eventEnergy );
  
    Float_t fLFreq;
    Float_t fEFreq;
    Float_t fLPhase;
    Float_t fEPhase;
  };

  // Private data:
  TFile*   fFile;
  TNtuple* fTree;
  TRandom* fRandom;
  TRandom* fParameterChooser;


  BlindingFunction fBlindingFunction;

};
#endif // BLINDER_H

---