NuMadAnalysis.cxx

Go to the documentation of this file.

// Coded by ... Dec/2005       
//                                                                    
// Contact: 

#include "TClonesArray.h"

#include "MessageService/MsgService.h"
#include "Mad/MadQuantities.h"

#include "NuMuBar/NuMadAnalysis.h"

using std::endl;
using std::cout;
using std::map;
using std::vector;

CVSID("$Id: NuMadAnalysis.cxx,v 1.1 2007/08/20 21:45:40 hartnell Exp $");

//......................................................................

NuMadAnalysis::NuMadAnalysis() : MadQuantities()
{
  MSG("NuMadAnalysis",Msg::kInfo)
    <<"Running NuMadAnalysis Constructor..."<<endl;
  
  //I don't want to give it a chain...
  Bool_t chainSR=false;

  //mad base does this:

  fCurRun = -1;
  lastEntry = -1;
  if(!chainSR) {
    record = 0;
    emrecord = 0;
    mcrecord = 0;
    threcord = 0;
    strecord = 0;
    Clear();
    whichSource = -1;
    isMC = true;
    isTH = true;
    isEM = true;
    isST = false;
    Nentries = -1;
    fCurRun = -1;
    return;
  }

}

//......................................................................

NuMadAnalysis::~NuMadAnalysis()
{
  MSG("NuMadAnalysis",Msg::kDebug)
    <<"Running NuMadAnalysis Destructor..."<<endl;
  

  MSG("NuMadAnalysis",Msg::kDebug)
    <<"Finished NuMadAnalysis Destructor"<<endl;
}

//......................................................................

void NuMadAnalysis::SetEntry(const NtpStRecord *ntp)
{
  MSG("NuMadAnalysis",Msg::kVerbose)
    <<"SetEntry..."<<endl;


  //have to mess about const_casting everything because the 
  //base class members aren't const for some reason

  //set the data member
  strecord=const_cast<NtpStRecord*>(ntp);
  isST=true;

  //strecord=const_cast<NtpStRecord*>(strecord);

  //set other things:
  if(isST){
    ntpHeader=&(strecord->GetHeader());
    mcHeader=const_cast<NtpMCSummary*>(&(strecord->mchdr));
    shieldSummary=const_cast<NtpSRShieldSummary*>(&(strecord->vetohdr));
    eventSummary=const_cast<NtpSREventSummary*>( &(strecord->evthdr));
    dmxStatus=const_cast<NtpSRDmxStatus*>(&(strecord->dmxstatus));
    detStatus=const_cast<NtpSRDetStatus*>(&(strecord->detstatus));
    if(mcHeader->nmc==0) {
      isMC = false; isTH = false;
    }
    else {
      isMC = true;
      TClonesArray &array = *(strecord->thstp);
      if(array.GetEntries()>0) isTH = true;      
    }
  }
  fCurRun = ntpHeader->GetRun();

  LoadSlice(0);  
  LoadEvent(0);
  LoadTrack(0);
  LoadCluster(0);
  LoadShower(0);
  LoadShieldStrip(0);
  
  MSG("NuMadAnalysis",Msg::kVerbose)
    <<"SetEntry... done!"<<endl;
}

//......................................................................

---