FillTruth.h

Go to the documentation of this file.

#ifndef ANP_FILLTRUTH_H
#define ANP_FILLTRUTH_H

// $Id: FillTruth.h,v 1.11 2008/01/31 22:18:17 rustem Exp $
//
// This algorithm fills all Monte-Carlo truth information. It is a 
// most complicated AlgStore algorithm. The following information
// is processed by this algorithm:
//
// 1) Neutrino truth information from beam and detector simulation,
//
// 2) Stdhep truth information for recorded stdhep blocks,
// 
// 3) Mapping between true and reconstructed objects (events and tracks),
// 
// 4) Input parameters for neugen cross section weights
//

// Local
#include "PhysicsNtuple/AlgStore.h"
#include "PhysicsNtuple/Truth.h"

class NtpMCTruth;
class NtpStRecord;
class NtpTHEvent;
class NtpTHTrack;
class NtpTHShower;
class TClonesArray;

namespace Anp
{
   class FillTruth : public AlgStore
   {
   public:
      
      FillTruth();
      virtual ~FillTruth();     
      
      bool Run(Record &record, TObject *ptr = 0);

      void Config(const Registry &reg);

   private:
      
      const Truth Fill(const NtpMCTruth &ntptruth) const;

      void FillStdHep(Record &record, const TClonesArray &array) const;

      void Fill(Truth &truth, const vector<const NtpTHEvent *>  &vec);
      void Fill(Truth &truth, const vector<const NtpTHTrack *>  &vec);
      void Fill(Truth &truth, const vector<const NtpTHShower *> &vec);
      
      bool IsConsistent(const Record &record) const;

   private:

      bool fStdHep;
      bool fCheck;

      int fHepKey;

      mutable unsigned int fNTruth;
      mutable unsigned int fNStHep;
      mutable unsigned int fNRecoEvent;
      mutable unsigned int fNRecoTrack;
      mutable unsigned int fNRecoShower;
   };

   template<class T> 
      std::vector<const T *> FillRecoTruth(const TClonesArray *array);
}

#endif

---