THcNPSApparatus.cxx

Go to the documentation of this file.

 
#include "THcNPSApparatus.h"
#include "THaNonTrackingDetector.h"
 
#include "TDatime.h"
 
#include "THcGlobals.h"
#include "THcParmList.h"
 
#ifdef WITH_DEBUG
#include <iostream>
#endif
using namespace std;
 
THcNPSApparatus::THcNPSApparatus(
  const char* name, const char* description
) :
  THaApparatus(name, description),
  fKwPrefix(""),
  fStagesDone(0), fListInit(kFALSE)
{
  // Constructorn
 
  fTracks = 0;
  fNonTrackingDetectors = new TList;
  // Construct the kwPrefix here. Better place would be in `Setup` method.
  string kwPrefix = name;
  std::transform(kwPrefix.begin(), kwPrefix.end(), kwPrefix.begin(), ::tolower);
  fKwPrefix = kwPrefix;
}
 
 
THcNPSApparatus::~THcNPSApparatus() {}
 
//_____________________________________________________________________________
void THcNPSApparatus::ListInit()
{
  // Initialize lists of specialized detectors. 
  // Private function called by Reconstruct().
 
  fNonTrackingDetectors->Clear();
 
  TIter next(fDetectors);
  while( THaDetector* theDetector = 
         static_cast<THaDetector*>( next() )) {
 
    //    if( theDetector->IsTracking() )
    //      fTrackingDetectors->Add( theDetector );
    //    else
    fNonTrackingDetectors->Add( theDetector );
 
      //    if( theDetector->IsPid() )
      //      fPidDetectors->Add( theDetector );
  }
 
  // Set up PIDinfo and vertex objects that can be associated with tracks
 
  //  UInt_t ndet  = GetNpidDetectors();
  //  UInt_t npart = GetNpidParticles();
  //  TClonesArray& pid  = *fTrackPID;
 
  //  for( int i = 0; i < kInitTrackMultiplicity; i++ ) {
  //    new( pid[i] )  THaPIDinfo( ndet, npart );
  //  }
  
  fListInit = kTRUE;
}
//_____________________________________________________________________________
Int_t THcNPSApparatus::CoarseReconstruct()
{
  // 2nd step: Coarse processing.  Pass the coarse tracks to the remaining
  // detectors for any processing that can be done at this stage.
  // This may include clustering and preliminary PID.
  // PID information is tacked onto the tracks as a THaPIDinfo object.
 
  if( !fListInit )
    ListInit();
  
  // No Tracking in NPS apparatus
  //  if( !IsDone(kCoarseTrack))a
  //    CoarseTrack();
 
  TIter next( fNonTrackingDetectors );
  while( THaNonTrackingDetector* theNonTrackDetector =
         static_cast<THaNonTrackingDetector*>( next() )) {
#ifdef WITH_DEBUG
    if( fDebug>1 ) cout << "Call CoarseProcess() for " 
                        << theNonTrackDetector->GetName() << "... ";
#endif
    theNonTrackDetector->CoarseProcess( *fTracks );
#ifdef WITH_DEBUG
    if( fDebug>1 ) cout << "done.\n";
#endif
  }
 
  fStagesDone |= kCoarseRecon;
  return 0;
}
 
//_____________________________________________________________________________
Int_t THcNPSApparatus::Reconstruct()
{
  // This method implements a fairly generic algorithm for spectrometer
  // reconstruction which should be useful for most situations.
  // Additional, equipment-specific processing can be done in
  // a derived class that calls this method first.
  //
  // The algorithm is as follows:
  //
  // For all tracking detectors:
  //   CoarseTrack()
  // For all non-tracking detectors:
  //   CoarseProcess()
  // For all tracking detectors:
  //   FineTrack()
  // Reconstruct tracks to target
  // For all non-tracking detectors:
  //   FineProcess()
  // Compute additional attributes of tracks (e.g. momentum, beta)
  // Find "Golden Track"
  // Combine all PID detectors to get overall PID for each track
  //
 
  // Do prior analysis stages if not done yet
  //  if( !IsDone(kTracking))
  //    Track();
 
  // Fine processing.  Pass the precise tracks to the
  // remaining detectors for any precision processing.
  // PID likelihoods should be calculated here.
 
  TIter next( fNonTrackingDetectors );
  while( THaNonTrackingDetector* theNonTrackDetector =
         static_cast<THaNonTrackingDetector*>( next() )) {
#ifdef WITH_DEBUG
    if( fDebug>1 ) cout << "Call FineProcess() for " 
                        << theNonTrackDetector->GetName() << "... ";
#endif
    theNonTrackDetector->FineProcess( *fTracks );
#ifdef WITH_DEBUG
    if( fDebug>1 ) cout << "done.\n";
#endif
  }
 
  // Compute additional track properties (e.g. beta)
  // Find "Golden Track" if appropriate.
  
  //TrackCalc();
 
  // Compute combined PID
 
  //  if( fPID ) CalcPID();
 
  fStagesDone |= kReconstruct;
  return 0;
}
 
 
TString THcNPSApparatus::GetKwPrefix() {
  return fKwPrefix;
}
 
Int_t THcNPSApparatus::ReadRunDatabase(const TDatime& date) {
  
  return kOK;
}
 
 
Int_t THcNPSApparatus::ReadDatabase(const TDatime& date) {
 
  Double_t shmsThetaLab;
  Double_t npsThetaOffset;
  DBRequest list[]={
    {"stheta_lab", &shmsThetaLab, kDouble},
    {0}
  };
  gHcParms->LoadParmValues((DBRequest*)&list);
 
  DBRequest listprefixed[]={
    {"_theta_offset", &npsThetaOffset, kDouble},
    {0}
  };
  gHcParms->LoadParmValues((DBRequest*)&listprefixed,fKwPrefix.c_str());
 
  fTheta_lab = shmsThetaLab + npsThetaOffset;
  
  return kOK;
}
 
 
Int_t THcNPSApparatus::DefineVariables(THaAnalysisObject::EMode mode) {
  if (mode == kDefine && fIsSetup) return kOK;
  fIsSetup = (mode == kDefine);
 
  return kOK;
  //  std::vector<RVarDef> vars;
  //  RVarDef end {0};
  //  vars.push_back(end);
 
  //  return DefineVarsFromList(vars.data(), mode);
}
 
 
THcNPSApparatus::THcNPSApparatus() {}
 
 
ClassImp(THcNPSApparatus)