ScintHitAna.cxx

Go to the documentation of this file.

// $Id: ScintHitAna.cxx,v 1.3 2006/12/01 20:12:11 rhatcher Exp $
//
// FILL_IN: [Document your code!!]
//
// n.tagg1@physics.ox.ac.uk
#include "ScintHitAna.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "JobControl/JobCModuleRegistry.h" // For JOBMODULE macro

#include "Record/SimSnarlRecord.h"
#include "Record/SimSnarlHeader.h"
#include "UgliGeometry/UgliGeomHandle.h"

JOBMODULE(ScintHitAna, "ScintHitAna",
          "ScintHitAna");


ClassImp(MyGeom)

CVSID("$Id: ScintHitAna.cxx,v 1.3 2006/12/01 20:12:11 rhatcher Exp $");
//......................................................................

ScintHitAna::ScintHitAna()
{
  //======================================================================
  // Make a tree.
  //======================================================================
  fTreeFile = new TFile("scinthits.root","RECREATE");
  fTree = new TTree("shtree","shtree");
  fTree->Branch("scinthit","DigiScintHit",&fHit,32000,1);
  fTree->Branch("geom","MyGeom",&fGeom,32000,1);
  fTree->Branch("event",&fEvent,"event/I");
  fEvent = 0;
  fGeom = new MyGeom;
}

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

ScintHitAna::~ScintHitAna()
{
//======================================================================
// Close the file
//======================================================================
  fTreeFile->cd();
  fTree->Write();
  fTreeFile->Close();
  delete fTreeFile;
  delete fTree;
}

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

JobCResult ScintHitAna::Ana(const MomNavigator* mom)
{
//======================================================================
// Fill the tree.
//======================================================================
  
  SimSnarlRecord* simsnarl = 0;
  TObject* tobj;
  TIter    fragiter = mom->FragmentIter();
  
  // Get the simsnarl.
  while( ( tobj = fragiter.Next() ) ) {
    simsnarl = dynamic_cast<SimSnarlRecord*>(tobj);
    if(simsnarl) break;
  }
  const SimSnarlHeader* simHeader = simsnarl->GetSimSnarlHeader();
  if(simHeader ==0){
    MSG("Photon",Msg::kError) << "Cannot find SimSnarlHeader in SimSnarl." << endl;
    return JobCResult::kFailed;
  }
  
  // We need a context to do table lookups
  VldContext simContext = simHeader->GetVldContext();
  
  
  // Get the DigiScintHits.
  // Get the scint hit array.
  const TObjArray* hitArray = 
    dynamic_cast<const TObjArray*>(simsnarl->FindComponent(0,"DigiScintHits"));
  if(hitArray==0) {
    cout << "Can't find scint hit array.\n"; 
    return JobCResult::kPassed;  
  };
  
  
  // Fill the tree.
  UgliGeomHandle ugli(simContext);
  
  TIter hitarrayIter(hitArray);
  while( (tobj = hitarrayIter.Next()) ) {
    const DigiScintHit* scinthit = dynamic_cast<DigiScintHit*>(tobj);
    if(scinthit) {
      fHit = const_cast<DigiScintHit*>(scinthit);

      UgliStripHandle strip = ugli.GetStripHandle(scinthit->StripEndId());
      fGeom->length = strip.GetHalfLength();
      fGeom->width  = strip.GetHalfWidth();
      fGeom->thickness = strip.GetHalfThickness();

      fTree->Fill();

    }
  }
  

  fEvent++;

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

---