numuCCBwdSelection.hxx

#ifndef numuCCBwdSelection_h
#define numuCCBwdSelection_h

#include "SelectionBase.hxx"
#include "Parameters.hxx"
#include "CutUtils.hxx"
#include "ToyBoxTracker.hxx"

class numuCCBwdSelection: public SelectionBase{

public:
  numuCCBwdSelection(bool forceBreak=true, SubDetId::SubDetEnum aFV = SubDetId::kFGD1);
  virtual ~numuCCBwdSelection(){}

  void DefineSteps();
  void DefineDetectorFV();  
  void InitializeEvent(AnaEventC& event);
  inline ToyBoxB* MakeToyBox() { return new ToyBoxTracker();  }
  bool FillEventSummary(AnaEventC& event, Int_t allCutsPassed[]);
  
  SampleId::SampleEnum GetSampleEnum(){
    if( _detectorFV[0] == SubDetId::kFGD1 ){
     return SampleId::kFGD1NuMuCCBwd;
    }else if( _detectorFV[0] == SubDetId::kFGD2 ){
     return SampleId::kFGD2NuMuCCBwd;
    }
    return SampleId::kUnassigned;   
  }
  
  bool CheckRedoSelection(const AnaEventC& event, const ToyBoxB& PreviousToyBox, Int_t& redoFromStep);
  bool IsRelevantRecObjectForSystematic(const AnaEventC& event, AnaRecObjectC* recObj, SystId_h systId, Int_t branch) const;
  bool IsRelevantSystematic(const AnaEventC& event, const ToyBoxB& box, SystId_h systId, Int_t branch) const;

protected:

  SubDetId::SubDetEnum _aFV;

};


namespace numuCCBwd{

  /// reduced FGD1 FV for Bwd
  const   Float_t FVdefminFGD1[3] = {57.66, 57.66, 0.0};
  const   Float_t FVdefmaxFGD1[3] = {57.66, 57.66, 0.0};

  /// reduced FGD2 FV for Bwd
  const   Float_t FVdefminFGD2[3] = {57.66, 57.66, 0.0};
  const   Float_t FVdefmaxFGD2[3] = {57.66, 57.66, 0.0};
  
  //Modified track quality cut
  bool TrackQualityCut(const AnaTrackB& track, SubDetId::SubDetEnum tpc);
  
  // In FGD FV
  SubDetId::SubDetEnum InFGDFV( const Float_t pos[4] );
  
  // Ends in FGD FV
  SubDetId::SubDetEnum EndsInFGDFV(const AnaTrackB& track);
  
  /// Event failed FGD1 selection when running FGD2 selection in a FGD1 + FGD2 analysis
  //********************************************************************
  class Exclusive: public StepBase{
  //********************************************************************
  public:
    using StepBase::Apply;
    bool Apply(AnaEventC& event, ToyBoxB& box) const;
    StepBase* MakeClone(){return new Exclusive();}
  };
  
  /// Muon candidate has a positive time of flight
  //********************************************************************
  class ToF: public StepBase{
  //********************************************************************
  public:
    ToF() {
      _FGD1_FGD2  = (Float_t) ND::params().GetParameterD("numuCCBwdAnalysis.Cuts.ToF.FGD1_FGD2"); 
      _P0D_FGD1  = (Float_t) ND::params().GetParameterD("numuCCBwdAnalysis.Cuts.ToF.P0D_FGD1"); 
      _ECal_FGD1  = (Float_t) ND::params().GetParameterD("numuCCBwdAnalysis.Cuts.ToF.ECal_FGD1"); 
      _ECal_FGD2  = (Float_t) ND::params().GetParameterD("numuCCBwdAnalysis.Cuts.ToF.ECal_FGD2"); 
    }
        
    
    using StepBase::Apply;
    bool Apply(AnaEventC& event, ToyBoxB& box) const;
    StepBase* MakeClone(){return new ToF();}
    
  private:
    Float_t _FGD1_FGD2, _P0D_FGD1, _ECal_FGD1, _ECal_FGD2;
  };
  
  /// Track has more than one node in ecal or p0d
  //********************************************************************
  class NoShower: public StepBase{
  //********************************************************************
  public:
    using StepBase::Apply;
    bool Apply(AnaEventC& event, ToyBoxB& box) const;
    StepBase* MakeClone(){return new NoShower();}
  };
  
  /// Timing available
  //********************************************************************
  class AvailableToF: public StepBase{
  //********************************************************************
  public:
    using StepBase::Apply;
    bool Apply(AnaEventC& event, ToyBoxB& box) const;
    StepBase* MakeClone(){return new AvailableToF();}
  };
  
  /// Muon PID cut
  //********************************************************************
  class MuonPIDCut: public StepBase{
  //********************************************************************
  public:
    MuonPIDCut(){
      _muon_pid = (Float_t)ND::params().GetParameterD("numuCCBwdAnalysis.Cuts.MuonPID"); 
    }
    
    using StepBase::Apply;
    bool Apply(AnaEventC& event, ToyBoxB& box) const;
    StepBase* MakeClone(){return new MuonPIDCut();}
    
  private:
    Float_t _muon_pid;
  };
  
  /// Define vertex
  //********************************************************************
  class FindVertexAction: public StepBase{
  //********************************************************************
  public:
    using StepBase::Apply;
    bool Apply(AnaEventC& event, ToyBoxB& box) const;
    StepBase* MakeClone(){return new FindVertexAction();}
  };
  
  /// Switch back TPC tracks sense
  //********************************************************************
  class SwitchBackAction: public StepBase{
  //********************************************************************
  public:
    using StepBase::Apply;
    bool Apply(AnaEventC& event, ToyBoxB& box) const;
    StepBase* MakeClone(){return new SwitchBackAction();}
  };
  
  /// Find Muon candidate
  //********************************************************************
  class FindMuonCandidateAction: public StepBase{
  //********************************************************************
  public:
    using StepBase::Apply;
    bool Apply(AnaEventC& event, ToyBoxB& box) const;
    StepBase* MakeClone(){return new FindMuonCandidateAction();}

  };

}
  
#endif