FillMRCC.cxx

Go to the documentation of this file.

//$Id: FillMRCC.cxx,v 1.7 2010/01/08 23:02:13 jyuko Exp $
// C++
#include <vector>

// ROOT
#include "TClonesArray.h"

// MINOS
#include "Conventions/BeamType.h"
#include "DataUtil/infid.h"
#include "MCNtuple/NtpMCStdHep.h"
#include "MCNtuple/NtpMCTruth.h"
#include "MuonRemoval/NtpMRRecord.h"
#include "MuonRemoval/NtpMREvent.h"
#include "MCReweight/NeugenWeightCalculator.h"
#include "MCReweight/ReweightHelpers.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "StandardNtuple/NtpStRecord.h"

// Local
#include "PhysicsNtuple/Default.h"
#include "PhysicsNtuple/Factory.h"
#include "PhysicsNtuple/Store/FillMRCC.h"

CVSID("$Id: FillMRCC.cxx,v 1.7 2010/01/08 23:02:13 jyuko Exp $");
REGISTER_ANP_OBJECT(AlgStore,FillMRCC)

using namespace std;

//-----------------------------------------------------------------------------------------
Anp::FillMRCC::FillMRCC()
  :fDebug(true)
{    
}

//-----------------------------------------------------------------------------------------
Anp::FillMRCC::~FillMRCC()
{
}

//-----------------------------------------------------------------------------------------
bool Anp::FillMRCC::Run(Record &record, TObject *ptr)
{
   //
   // Get and fill Monte-Carlo truth and StdHep information
   //
  const NtpMRRecord* ntmrcc= dynamic_cast<NtpMRRecord *>(ptr);
  if(!ntmrcc){
    const MomNavigator *mom = dynamic_cast<const MomNavigator *> (ptr);
    if(mom)
      {
        ntmrcc = dynamic_cast<NtpMRRecord *>(mom -> GetFragment("NtpMRRecord")); 
        if(!ntmrcc)
          {
            MSG("FillAlg", Msg::kError) << "Failed to find NtpMRRecord pointer 1" << endl;
            return false;
          }
      }
    else
      {
        MSG("FillAlg", Msg::kError) << "Failed to find MomNavigator pointer" << endl;
        return false;
      }      
  }
  if(!ntmrcc)
    {
      MSG("FillAlg", Msg::kError) << "Failed to get NtpMRRecord pointer 2" << endl;
      return false;
    }

  TClonesArray& evtTca=*(ntmrcc->mrevt);
  //  TClonesArray& truTca=*(ntmrcc->mrtru);

  const Int_t numEvts=evtTca.GetEntriesFast();
  for(int jj=0;jj<numEvts;jj++){
    const NtpMREvent& mrccevt=   *dynamic_cast<NtpMREvent*>(evtTca[jj]);
    MRCCInfo minfo = Fill(mrccevt);
    minfo.infoindex = jj;
    record.Add(minfo);
  }
  if(fDebug)  cout<<"FillMRCC::Run - End "<<endl;
  return true;
}

//----------------------------------------------------------------------
void Anp::FillMRCC::Config(const Registry &reg)
{

  Anp::Read(reg,"FillMRCCDebug",fDebug);
   if(reg.KeyExists("PrintConfig"))
   {
     cout << "FillMRCC::Config" << endl;
     cout << "Debug    "<<fDebug<< endl;
   }
}

//-----------------------------------------------------------------------------------------
const Anp::MRCCInfo Anp::FillMRCC::Fill(const NtpMREvent &mrevent) const
{

  MRCCInfo mrccinfo;

  mrccinfo.mrcc_index = mrevent.best_event;
  mrccinfo.orig_index = mrevent.orig_event;
  mrccinfo.index   =   mrevent.index;
  mrccinfo.ndigit  =   mrevent.ndigit;
  mrccinfo.nstrip  =   mrevent.nstrip;
  //  mrccinfo.stp     =   mrevent.stp;
  //  mrccinfo.mrstp   =   mrevent.mrstp;
  mrccinfo.best_purity =      mrevent.best_purity;
  mrccinfo.best_complete =    mrevent.best_complete;
  mrccinfo.best_purity_phw =  mrevent.best_purity_phw;
  mrccinfo.best_complete_phw= mrevent.best_complete_phw;
  mrccinfo.elec_complete  =   mrevent.elec_complete;
  mrccinfo.elec_complete_phw= mrevent.elec_complete_phw;
  mrccinfo.vtxx  =   mrevent.vtxx;
  mrccinfo.vtxy  =   mrevent.vtxy;
  mrccinfo.vtxz  =   mrevent.vtxz;
  mrccinfo.vtxp  =   mrevent.vtxp;
  mrccinfo.npln  =   mrevent.npln;
  mrccinfo.prng  =   mrevent.prng;
  mrccinfo.pcrv  =   mrevent.pcrv;
  mrccinfo.pvdx  =   mrevent.pvdx;
  mrccinfo.pvdy  =   mrevent.pvdy;
  mrccinfo.pvdz  =   mrevent.pvdz;
  mrccinfo.fitp  =   mrevent.fitp;
  mrccinfo.endc  =   mrevent.endc;
  mrccinfo.pass  =   mrevent.pass;
  mrccinfo.pmux  =   mrevent.pmux;
  mrccinfo.pmuy  =   mrevent.pmuy;
  mrccinfo.pmuz  =   mrevent.pmuz;
  mrccinfo.mxpl  =   mrevent.mxpl;

   return mrccinfo;
}

---