SimCheck.cxx

Go to the documentation of this file.

#include "SimCheck.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "JobControl/JobCModuleRegistry.h" // For JOBMODULE macro

#include "Digitization/DigiScintHit.h"
#include "Digitization/DigiPE.h"
#include "Digitization/DigiSignal.h"
#include "Record/SimSnarlRecord.h"

#include <set>

using namespace std;

JOBMODULE(SimCheck, "SimCheck",
          "SimCheck");
CVSID("$Id: SimCheck.cxx,v 1.3 2009/02/28 21:46:12 gmieg Exp $");
//......................................................................

SimCheck::SimCheck() 
{
}
//......................................................................

SimCheck::~SimCheck()
{
}

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

JobCResult SimCheck::Ana(const MomNavigator* mom)
{

  SimSnarlRecord* simsnarl = 0;
  TObject* tobj;
  TIter    fragiter = mom->FragmentIter();

  // Get the simsnarl.
  while( ( tobj = fragiter.Next() ) ) {
    simsnarl = dynamic_cast<SimSnarlRecord*>(tobj);
    if(simsnarl) break;
  }
  
  // Verify there IS a simsnarl.
  if(!simsnarl) {
    MSG("Photon",Msg::kError) << "No SimSnarl found. You must run RerootToTruthModule()!" << endl;
    return JobCResult::kFailed;
  }


  const TObjArray* hitlist = 
    dynamic_cast<const TObjArray*>(simsnarl->FindComponent(0,"DigiScintHits"));
  if(!hitlist) return JobCResult::kFailed;

  //const TObjArray* pelist = 
  //  dynamic_cast<const TObjArray*>(simsnarl->FindTemporary(0,"DigiListPe"));

  //if(!pelist) return JobCResult::kFailed;

  const TObjArray* signallist = 
    dynamic_cast<const TObjArray*>(simsnarl->FindComponent(0,"DigiSignal"));
  if(!signallist) return JobCResult::kFailed;

  std::set<const DigiScintHit*> hits;
  std::set<const DigiScintHit*> sighits;
  std::set<const DigiScintHit*> intersection;
  std::set<const DigiScintHit*> difference1;
  std::set<const DigiScintHit*> difference2;
  

  TIter hitarrayIter(hitlist);
  while( (tobj = hitarrayIter.Next()) ) {
    const DigiScintHit* scinthit = dynamic_cast<const DigiScintHit*>(tobj);
    if(scinthit) {
      hits.insert(scinthit);
    }
  } 

  TIter sigarrayIter(signallist);
  while( (tobj = sigarrayIter.Next()) ) {
    const DigiSignal* signal = dynamic_cast<const DigiSignal*>(tobj);
    if(signal) {
      for(UInt_t i=0; i<signal->GetNumberOfHits(); i++) {
        sighits.insert(signal->GetHit(i));
      }
    }
  } 

//gmi These unused declarations are giving an error in GCC 4.4.0
//gmi  set_intersection(hits.begin(),hits.end(),
//gmi              sighits.begin(), sighits.end(),
//gmi              inserter(intersection,intersection.begin()) );
//gmi
//gmi  set_difference(hits.begin(), hits.end(),          
//gmi            sighits.begin(), sighits.end(),
//gmi            inserter(difference1,difference1.begin()) );
//gmi
//gmi  set_difference(sighits.begin(), sighits.end(),
//gmi            hits.begin(), hits.end(),               
//gmi            inserter(difference2,difference2.begin()) );

  cout << "SimCheck-----------------------------------" << endl;
  cout << "In Hit list: " << hits.size() << endl;
  cout << "In DigiSignals: " << sighits.size() << endl;
  cout << "Intersection: " << intersection.size() << endl;
  cout << "Diff1: " << difference1.size() << endl;
  cout << "Diff2: " << difference2.size() << endl;
  

  return JobCResult::kPassed; // kNoDecision, kFailed, etc.
}

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

const Registry& SimCheck::DefaultConfig() const
{
  static Registry r; // Default configuration for module

  // Set name of config
  std::string name = this->GetName();
  name += ".config.default";
  r.SetName(name.c_str());

  // Set values in configuration
  r.UnLockValues();
  r.LockValues();

  return r;
}

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

void SimCheck::Config(const Registry& )
{
}

---