DemoNtupleModule.cxx

Go to the documentation of this file.

//
// DemoRecordModule.cxx
//
// A JobControl Module for demonstrating record construction.
//
#include <iostream>
#include <cassert>
using std::endl;

#include "TClonesArray.h"
#include "TVector3.h"

#include "JobControl/JobCModuleRegistry.h"
#include "JobControl/JobCommand.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "RawData/RawRecord.h"
#include "RawData/RawDaqSnarlHeader.h"

#include "Record/test/DemoNtupleModule.h"
#include "Record/demo/DemoNtupleRecord.h"

ClassImp(DemoNtupleModule)

//   Definition of static data members
//   *********************************

CVSID("$Id: DemoNtupleModule.cxx,v 1.4 2005/02/01 18:51:22 tagg Exp $");
JOBMODULE(DemoNtupleModule, "DemoNtupleModule",
         "A module for demonstrating ntuple construction.");


// Definition of methods (alphabetical order)
// ***************************************************

JobCResult DemoNtupleModule::Reco(MomNavigator *mom) {
  //
  //  Purpose:  Create ntuple record.
  //
  //  Arguments: none.
  //  
  //  Return: none.
  //
  //  Contact:   S. Kasahara
  // 


  JobCResult result(JobCResult::kPassed);  
  MSG("DemoNtuple",Msg::kVerbose) << "DemoNtupleModule::Reco" << endl;

  // Check that mom exists.
  assert(mom);

  // Create a DemoNtupleRecord from a RawRecord (DaqSnarl) 
  // Should iterate over rawrecords in mom, but this is just a demo
  RawRecord *rawRec = dynamic_cast<RawRecord*>(mom->GetFragment("RawRecord"));
  if ( !rawRec ) {
    MSG("RawRecord",Msg::kWarning) << "No RawRecord in Mom " <<endl;
    result.SetWarning().SetFailed();
    return result;
  }

  const RawDaqSnarlHeader* oldHdr = dynamic_cast<const RawDaqSnarlHeader*>
                                         (rawRec -> GetRawHeader());
  // Use validity context to create header for ntuple record
  RecHeader newHdr(oldHdr->GetVldContext());

  // Create new record
  DemoNtupleRecord* ntp = new DemoNtupleRecord(newHdr);
  ntp -> SetMyData(oldHdr->GetSnarl());

  // For each new vertex, a TVector3 object is created and added to
  // the TClonesArray.
  int nvectors = 0;
  for ( int i = 0; i < 2; i++ ) {
    // This uses "new with placement" to reuse previously allocated memory
    if ( oldHdr->GetSnarl()%2 == 0 ) {
      // just for fun, every other entry has different values
      new((ntp->GetVectors())[nvectors++]) TVector3(1,2,3);
    }
    else {
      new((ntp->GetVectors())[nvectors++]) TVector3(4,5,6);
    }
  }

  mom -> AdoptFragment(ntp);

  return result;

}

---