NtpStRecord.h

Go to the documentation of this file.

// 
// NtpStRecord
//
// NtpStRecord is a standard reconstruction & mc ntuple record.
//
// S. Kasahara 11/04
// Based on Roy Lee's CandEventSR/TTreeSR class

#ifndef NTPSTRECORD_H
#define NTPSTRECORD_H

#include <iosfwd>
#include <vector>

#ifndef RECRECORDIMP_H
#include "Record/RecRecordImp.h" // base class
#endif
#ifndef RECCANDHEADER_H
#include "Record/RecCandHeader.h" 
#endif
#include "MCNtuple/NtpMCSummary.h"
#include "MCNtuple/NtpMCPhotonResult.h"
#include "MCNtuple/NtpMCDetSimResult.h"
#include "CandNtupleSR/NtpSREventSummary.h"
#include "CandNtupleSR/NtpSRShieldSummary.h"
#include "CandNtupleSR/NtpSRCosmicRay.h"
#include "CandNtupleSR/NtpSRDmxStatus.h"
#include "CandNtupleSR/NtpSRDetStatus.h"
#include "CandNtupleSR/NtpSRTimeStatus.h"
#include "CandNtupleSR/NtpSRCalStatus.h"
#include "CandNtupleSR/NtpSRDataQuality.h"
#include "Conventions/ReleaseType.h"

// ROOT forward declaration
class TClonesArray;

// MINOS forward declaration
class NtpSREvent;
class NtpSRShower;
class NtpSRStrip;
class NtpSRTrack;
class NtpMCStdHep;
class NtpMCTruth;
class NtpTHEvent;
class NtpTHTrack;
class NtpTHShower;

class NtpStRecord: public RecRecordImp<RecCandHeader> {

public:
  NtpStRecord();
  NtpStRecord(const RecCandHeader& header);
  virtual ~NtpStRecord(); 

  // State changing methods
  void Clear(Option_t* option = "");
  void ClearStrips(Option_t* option = "");
  
  // State testing methods
  // print methods
  virtual std::ostream& Print(std::ostream& os) const;
  virtual void Print(Option_t* option = "") const;
  ReleaseType::Release_t GetRelease() const;
 
  // standard reconstruction accessors
  std::vector< const NtpSREvent* > GetEvents() const;
  std::vector< const NtpSRTrack* >  GetTracks(int event_index = -1) const;
  std::vector< const NtpSRShower* > GetShowers(int event_index = -1) const;
  std::vector< const NtpSRStrip* >  GetStrips(int event_index = -1) const;

  // mc truth accessors
  std::vector< const NtpMCTruth* > GetMCTruths() const;
  std::vector< const NtpMCStdHep* > GetMCStdHeps(int mctruth_index = -1) const;

  // truth helper accessors
  std::vector< const NtpTHEvent* > GetTHEvents() const;
  std::vector< const NtpTHTrack* > GetTHTracks( int event_index = -1 ) const;
  std::vector< const NtpTHShower* > GetTHShowers(int event_index = -1 ) const;

  // specialized accessors
  const NtpMCTruth* GetEventMCTruth(int event_index) const;
  
 private:
  // private methods
  void Init();

 public:
  // Ntuple is treated like a C-struct with public data members and
  // rule-breaking field data members not prefaced by "f" and all
  // lowercase, by popular demand.

  NtpSREventSummary  evthdr;     // sr summary data
  NtpSRShieldSummary vetohdr;    // veto shield summary data
  NtpSRCosmicRay     crhdr;      // cr data, filled from last recons. trk vtx
  NtpSRDmxStatus     dmxstatus;  // status of demux
  NtpSRDetStatus     detstatus;  // status of detector
  NtpSRTimeStatus    timestatus; // status of timing system, esp ND
  NtpSRCalStatus     calstatus;  // status of calibration
  NtpSRDataQuality   dataquality;// quality of data
  NtpMCSummary       mchdr;      // mc summary data
  NtpMCPhotonResult  photon;     // summary data from photon transport mc
  NtpMCDetSimResult  detsim;     // summary data from det sim mc
  
  // Standard reconstruction arrays
  TClonesArray* vetostp; //-> array of shield strips of type NtpSRShieldStrip
  TClonesArray* vetoexp; //-> array of shield expected hits of type NtpSRShieldExpected 
  TClonesArray* deadchips; //-> array of dead chips of type NtpSRDeadChip

  TClonesArray* stp;  //-> array of strips of type NtpSRStrip 
  TClonesArray* slc;  //-> array of slices of type NtpSRSlice
  TClonesArray* clu;  //-> array of clusters of type NtpSRCluster
  TClonesArray* shw;  //-> array of showers of type NtpSRShower
  TClonesArray* trk;  //-> array of tracks of type NtpSRTrack
  TClonesArray* evt;  //-> array of events of type NtpSREvent

  // MC truth arrays
  TClonesArray* mc;       //-> array of NtpMCTruth 
  TClonesArray* stdhep;   //-> array of NtpMCStdHep 
  TClonesArray* digihit;  //-> array of NtpMCDigiScintHit

  // Truth helper arrays
  TClonesArray* thstp;     //-> array of NtpTHStrip  
  TClonesArray* thslc;     //-> array of NtpTHSlice  
  TClonesArray* thshw;     //-> array of NtpTHShower  
  TClonesArray* thtrk;     //-> array of NtpTHTrack
  TClonesArray* thevt;     //-> array of NtpTHEvent  

  ClassDef(NtpStRecord,10)              
};

#endif // NTPSTRECORD_H

---