ExodusValidate.cxx

Go to the documentation of this file.

// $Id: ExodusValidate.cxx,v 1.5 2001/03/13 19:50:24 rhatcher Exp $
//
// ExodusValidate
//
// ExodusValidate is a test harness for RerootExodus package
//
// Author:  R. Hatcher 2000.05.05
//

#include "RerootExodus/ExodusValidate.h"
#include "RerootExodus/RerootExodus.h"

#include "Plex/PlexStripEndId.h"
#include "Plex/PlexSEIdAltL.h"
#include "Plex/PlexHandle.h"

#include "Plex/PlexLoanPool.h"

#include "Validity/VldContext.h"

#include "MessageService/MsgService.h"
CVSID("$Id: ExodusValidate.cxx,v 1.5 2001/03/13 19:50:24 rhatcher Exp $");

ClassImp(ExodusValidate)

//______________________________________________________________________________
ExodusValidate::ExodusValidate()
{
   // Default constructor


   MsgService* msvc = MsgService::Instance();
   MsgStream*  mstr = msvc->GetStream("Exodus");
   
   mstr->SetLogLevel(Msg::kVerbose); // spew out everything for debugging
   mstr->SetFormat(Msg::kVerbose,0xffff);

}

//______________________________________________________________________________
ExodusValidate::~ExodusValidate()
{
   // delete all the owned sub-objects

}

//______________________________________________________________________________
void ExodusValidate::TestPlexHandle(void)
{
   // Test PlexHandle

   VldContext vldc = RerootExodus::BuildVldContext();
   MSG("Exodus",Msg::kVerbose) << 
      "TestPlexHandle dump VldContext " << endl << flush;
   vldc.Print("");

   PlexHandle ph(vldc);

   PlexLoanPool* plp = PlexLoanPool::Instance();
   plp->Print("");

}

//______________________________________________________________________________
void ExodusValidate::TestSteelPlane(void)
{
   // Test REROOT to/from Steel Plane #

   MSG("Exodus",Msg::kVerbose) << "TestSteelPlane " << endl;

   const int n = 21;
   int listpln[n] = {  -1,   0,   1, 
                   118, 119, 120, 121,
                   238, 239, 240, 241, 242, 243,
                   478, 479, 480, 481, 482, 483, 484, 485 };

   for (Int_t i = 0; i < n ; i++) {
      Int_t ipln = listpln[i];
      Int_t spln = RerootExodus::RerootToSteelPlane(ipln);
      Int_t rpln = RerootExodus::SteelToRerootPlane(spln);
      
      MSG("Exodus",Msg::kVerbose) << " reroot " << ipln <<
         " steel " << spln << " reroot " << rpln << endl;
   }

//#define PLANENUMBERINGTEST
#ifdef PLANENUMBERINGTEST
  Int_t nactive  = RerootExodus::NumberOfPlanes(kTRUE);
  Int_t npassive = RerootExodus::NumberOfPlanes(kFALSE);
  Int_t nplanes  = TMath::Max(nactive,npassive);
  Int_t testends = 0;  // set to "1" is you want to test handling bad RPLN #'s
  printf(" this detector has %d active (%d passive) planes\n",nactive,npassive);
  for (Int_t ipln_f77 = 1-testends; ipln_f77<=nplanes+testends; ipln_f77++) {
     if (ipln_f77<10 || (ipln_f77>238 && ipln_f77<250) || ipln_f77>476) {
        Int_t ipln_reroot = ipln_f77-1;
        Int_t steel    = RerootExodus::RerootToSteelPlane(ipln_reroot,kFALSE);
        Int_t active   = RerootExodus::RerootToSteelPlane(ipln_reroot,kTRUE);
        Int_t r_steel  = RerootExodus::SteelToRerootPlane(steel,kFALSE);
        Int_t r_active = RerootExodus::SteelToRerootPlane(active,kTRUE);
        printf(" plane f77=%3d reroot=%3d   steel=%3d  active=%3d  reroot %3d %3d\n",ipln_f77,ipln_reroot,steel,active,r_steel,r_active);
     }
  }
#endif

}

//______________________________________________________________________________
void ExodusValidate::TestStripCount(void)
{
   // Test REROOT strip counting

   for (Int_t rpln = 0; rpln < 11; rpln++) {
      MSG("Exodus",Msg::kInfo) << " REROOT plane " << rpln <<
         " is type '" << RerootExodus::ActvPlaneName(rpln) <<
         "' and has " << RerootExodus::NStripsNonZero(rpln) <<
         " non-zero length strips" << endl;
   }

}
//______________________________________________________________________________

---