GeoCheckOverlaps.cc

Go to the documentation of this file.

//
// GeoCheckOverlaps
//
// Package: (Geo) GeoGeometry
//
// S. Kasahara 10/2005
//
// Purpose: Class for use in validating the GeoGeometry package.

#include <iostream>
using namespace std;

#include "TGeoManager.h"
#include "GeoGeometry/test/GeoCheckOverlaps.h"
#include "UgliGeometry/UgliLoanPool.h"
#include "UgliGeometry/UgliGeomHandle.h"

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

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

//_____________________________________________________________________________
GeoCheckOverlaps::GeoCheckOverlaps(float precision, Int_t detmask, 
                                                    Int_t simmask) : 
                       fVldContext(),fGeoManager(0),fPrecision(precision),
                       fDetMask(detmask),fSimMask(simmask) {
  // Default constructor

}

//_____________________________________________________________________________
GeoCheckOverlaps::~GeoCheckOverlaps() {
  // Destructor
  // fGeoManager is not owned  
}

//_____________________________________________________________________________
bool GeoCheckOverlaps::TestCheckOverlaps() {
  //  Run CheckOverlaps on all volumes in detector
  bool pass = true;
  
  if ( !fGeoManager ) {
    cerr << "GeoCheckOverlaps::TestCheckOverlaps Error! fGeoManager is Null!" 
         << endl;
    pass = false;
    return pass;
  }
    
  fGeoManager->CheckOverlaps(fPrecision);
  
  TObjArray* overlaps = fGeoManager->GetListOfOverlaps();
  if ( overlaps ) {
    if ( overlaps -> GetEntriesFast() > 0 ) {
      cout << "Found " << overlaps->GetEntriesFast() 
           << " overlaps/extrusions." << endl;
      fGeoManager->PrintOverlaps();
      pass = false;
    }
    fGeoManager->ClearOverlaps(); // reset before next call
  }
   
  return pass;

}

//_____________________________________________________________________________
bool GeoCheckOverlaps::RunAllTests() {
  //
  //  Purpose:  Run sequence of GeoGeometry package validity tests.
  //
  //  Arguments: none.
  //
  //  Return:    pass (if all tests successful) or fail
  //
  //  Contact:   S. Kasahara
  // 

  // Test 3 detector types in sequence, with kData DB values
  const int ndet = 3;
  Detector::Detector_t dettype[ndet] = {Detector::kFar,
                                        Detector::kNear,
                                        Detector::kCalDet};

  const int nsim = 2;
  SimFlag::SimFlag_t simflag[nsim] = { SimFlag::kMC, SimFlag::kData };
      
  VldTimeStamp currentTime;
  
  cout << "*** GeoValidate::RunAllTests ***" << endl;
  cout << "*** CheckOverlaps w/Precision  = "
       << fPrecision << "(cm), Detector: "
       << Detector::MaskToString(fDetMask) 
       << ", SimFlag: " << SimFlag::MaskToString(fSimMask)
       << " ***" << endl;
  
 
  Bool_t passed[ndet][nsim];
  VldContext vldc[ndet][nsim];

  for ( int isim = 0; isim < nsim; isim++ ) {
    if ( !(simflag[isim] & fSimMask ) ) continue;
    for ( int idet = 0; idet < ndet; idet++ ) {
      if ( !(dettype[idet] & fDetMask) ) continue;
    
      passed[idet][isim] = true;
      
      if (dettype[idet] != Detector::kCalDet) {
        // Build vldc using current time
        fVldContext = VldContext(dettype[idet],simflag[isim],currentTime);
      }
      else {
        // Current time isn't valid for caldet
        fVldContext = VldContext(dettype[idet],simflag[isim],
                                 VldTimeStamp(2003,9,15,0,0,0));
      }
      vldc[idet][isim] = fVldContext; 
    
      cout << "\n*** Testing geometry with vld " 
           << fVldContext << " ***" << endl;
    
      UgliLoanPool::Instance()->SetUseGeo(true);
      UgliGeomHandle geo(fVldContext);
      fGeoManager = gGeoManager;
    
      // Run test methods
      cout << "*** TestCheckOverlaps ***" << endl;
      if ( TestCheckOverlaps() ) {
        cout << "*** All CheckOverlaps tests Passed ***" << endl;
      }
      else {
        cout << "*** At least one CheckOverlaps test FAILED ***" << endl;
        passed[idet][isim] = false;
      }
    
      std::string detstatus = " Passed all Tests ***";
      if ( !passed[idet][isim] ) detstatus = " FAILED at least one Test ***";
    
      cout << "\n*** Detector w/vld " << fVldContext << detstatus.c_str() 
           << endl;
    }
  }

  cout << "\n*** Finished all tests:" << endl;
  Bool_t allPassed = kTRUE;
  for ( int isim = 0; isim < nsim; isim++ ) {
    if ( !(simflag[isim] & fSimMask ) ) continue;
    for ( int idet = 0; idet < ndet; idet++ ) {
      if ( !(dettype[idet] & fDetMask) ) continue;

      std::string detstatus = " Passed all Tests ***";
      if ( !passed[idet][isim] ) {
        detstatus = " FAILED at least one Test ***";
        allPassed = kFALSE;
      }
      cout << "*** Detector w/vld " << vldc[idet][isim] 
           << detstatus.c_str() << endl;
    }
  }
  
  return allPassed;

}

---