ANtpInfoObjectFillerNC.h

Go to the documentation of this file.

// $Id: ANtpInfoObjectFillerNC.h,v 1.18 2009/09/12 14:56:20 tinti Exp $
//
// A class to easily get at the SR ntuple information
//
// B. Rebel 12/2005
#ifndef ANTPINFOOBJECTFILLERNC_H
#define ANTPINFOOBJECTFILLERNC_H

#include "TObjArray.h"
#include "TClonesArray.h"
#include "AnalysisNtuples/Module/ANtpInfoObjectFillerBeam.h"
#include "Conventions/Detector.h"
#include "PhysicsNtuple/Store/Interface.h"

class VldContext;
class VldTimeStamp;
class NtpMCRecord;
class NtpSRRecord;
class NtpStRecord;
class ANtpShowerInfoNC;
class ANtpTrackInfoNC;
class ANtpEventInfoNC;
class ANtpRecoNtpManipulator;

class ANtpInfoObjectFillerNC : public ANtpInfoObjectFillerBeam
{

 public:
  ANtpInfoObjectFillerNC();
  ANtpInfoObjectFillerNC(Detector::Detector_t detector);
  virtual ~ANtpInfoObjectFillerNC();
  
  bool FillInformation(int event, 
                       ANtpRecoNtpManipulator *ntpManipulator, 
                       ANtpEventInfoNC *eventInfo, 
                       ANtpTrackInfoNC *trackInfo,
                       ANtpShowerInfoNC *showerInfo, 
                       ANtpTruthInfoBeam *truthInfo);
  void SetDetector(Detector::Detector_t detector);
  void SetClusterArray(const TClonesArray *clusters);
  void InitializekNN(ANtpRecoNtpManipulator *ntpManipulator);
/*   void InitializeVHSTraining(const int numPlanes=20, */
/*                           const int numStrips=20); */

 private:
  
  void FillEventInformation(ANtpRecoNtpManipulator *ntpManipulator, 
                            NtpSREvent *ntpEvent, 
                            ANtpEventInfoNC *eventInfo);
  void FillTrackInformation(ANtpRecoNtpManipulator *ntpManipulator,
                            NtpSRTrack *ntpTrack, 
                            NtpSREvent *ntpEvent, 
                            ANtpTrackInfoNC *trackInfo,
                            ANtpEventInfoNC *eventInfo);
  void FillShowerInformation(NtpSRShower *ntpShower, 
                             ANtpShowerInfoNC *showerInfo,
                             ANtpEventInfoNC *eventInfo);
  void FillCrossOverInformation(NtpSRTrack *ntpTrack, 
                                NtpSRShower *ntpShower,
                                NtpSREvent *ntpEvent,
                                ANtpRecoNtpManipulator *ntpManipulator,
                                ANtpTrackInfoNC *trackInfo,
                                ANtpShowerInfoNC *showerInfo,
                                ANtpEventInfoNC *eventInfo);
  void FillMCInformation(NtpMCTruth *ntpMCTruth, 
                         TClonesArray *stdHepArray,
                         ANtpTruthInfoBeam *truthInfo);
  void FillStripToPixelMaps();
  void ResetXTalkArrays();
  void FillPlanePixelSignalArrays(NtpSREvent *ntpEvent);
  Float_t FindNearestNeighborPixelSignal(Int_t plane, 
                                         Int_t pixel, 
                                         Int_t pmt,
                                         Float_t planePixelSignal[][64][2]);
  void FindXTalkStrips(NtpSREvent *ntpEvent);
  void FindEarlyActivityWeight(NtpSREvent *ntpEvent, 
                               ANtpEventInfoNC *eventInfo);
  
  void GetEvtVtxWithFixup(ANtpRecoNtpManipulator* ntpManipulator,
                          int& vtxPlane, float& vtxZ) const;

  void GetStripEventTime(const NtpSREvent* ntpEvent,
                         int vtxPlane, float vtxZ,
                         double& evtTime,
                         double& ToFCorrectedEvtTime) const;
  void FillEventTimingAndActivityInformation(ANtpRecoNtpManipulator *ntpManipulator, 
                                             int event,
                                             ANtpEventInfoNC *eventInfo); 
/*   void FillVHSEvtInfo(ANtpRecoNtpManipulator* ntpManipulator, */
/*                    ANtpEventInfoNC*        eventInfo); */

  Int_t fPlaneToPMTMapWest[486];
  Int_t fPlaneToPMTMapEast[486];
  Int_t fStripToPixelMapWest[192];
  Int_t fStripToPixelMapEast[192];
  Int_t fStripToPixelMapNearU1[64];
  Int_t fStripToPixelMapNearV1[64];
  Int_t fStripToPixelMapNearU2[64];
  Int_t fStripToPixelMapNearV2[64];
  Int_t fStripToPixelMapNearU3[64];
  Int_t fStripToPixelMapNearV3[64];
  Float_t fPlanePixelEastSignal[486][64][2];
  Float_t fPlanePixelWestSignal[486][64][2];
  Int_t fStripIsXTalkWest[486][192];
  Int_t fStripIsXTalkEast[486][192];
  Int_t fPlaneCoverage[486];

  Anp::Interface *fAnpInterface;

  const TClonesArray *fClusterArray; // pointer to cluster array in NtpSt
  
  bool fkNNSet;  //flag for whether to set the kNN registry or not

  // Local copy of training vectors etc. for Kregg's image recognition
  int                 fVHSPlanes;    // #planes/window
  int                 fVHSStrips;    // #strips/window
  std::vector<double> fAvgCCe;       // Avg CCe   images
  std::vector<double> fAvgCCmu;      // Avg CCmu  images
  std::vector<double> fAvgCCtau;     // Avg CCtau images
  std::vector<double> fAvgNC;        // Avg NC    images
  std::vector<double> fPCCe;         // Hit probs for CCe   events
  std::vector<double> fPCCmu;        // Hit probs for CCmu  events
  std::vector<double> fPCCtau;       // Hit probs for CCtau events
  std::vector<double> fPNC;          // Hit probs for NC    events
  std::vector<double> fVarCCe;       // Strip variance CCe
  std::vector<double> fVarCCmu;      // Strip variance CCmu
  std::vector<double> fVarCCtau;     // Strip variance CCtau
  std::vector<double> fVarNC;        // Strip variance NC

  ClassDef(ANtpInfoObjectFillerNC, 6) //ANtpInfoObjectFillerNC
};

#endif // ANTPINFOOBECTFILLERNC_H

---