GeoMedium.cxx

Go to the documentation of this file.

//
// GeoMedium
//
// GeoMedium IsA TGeoMedium with additional data members to store medium
// specific cut and process flags.
//
// Author:  S. Kasahara 01/2008

#include "MessageService/MsgService.h"
#include "Util/UtilMCFlag.h"
using namespace UtilMCFlag;
#include "GeoGeometry/GeoMedium.h"

ClassImp(GeoMedium)

CVSID("$Id: GeoMedium.cxx,v 1.2 2008/01/09 07:00:05 schubert Exp $");

//_____________________________________________________________________________
GeoMedium::GeoMedium(const char* medname, Int_t medId, Int_t matId,
                     Int_t isvol, Geo::ESwimMethod ifield, Double_t fieldm, 
                     Double_t tmaxfd, Double_t stemax, Double_t deemax,
                     Double_t epsil, Double_t stmin) : 
  TGeoMedium(medname,medId,matId,isvol,ifield,fieldm,tmaxfd,stemax,
             deemax,epsil,stmin) {
  // Normal constructor. 

  MSG("Geo",Msg::kVerbose) << "GeoMedium normal ctor @ " << this
                           << " for " << medname << "." << endl;

}

//_____________________________________________________________________________
GeoMedium::GeoMedium(const char* medname, Int_t medId, const TGeoMaterial* mat,
                     Double_t *params) : TGeoMedium(medname,medId,mat,params) {
  // Normal constructor. 

  MSG("Geo",Msg::kVerbose) << "GeoMedium normal ctor @ " << this
                           << " for " << medname << endl;

}

//_____________________________________________________________________________
GeoMedium::GeoMedium() {
  // Default constructor. 

  MSG("Geo",Msg::kVerbose) << "GeoMedium default ctor @ " << this << endl;

}

//_____________________________________________________________________________
GeoMedium::~GeoMedium() {
  // Destructor. 

  MSG("Geo",Msg::kVerbose) << "GeoMedium destructor @ " << this << endl;

}

//_____________________________________________________________________________
void GeoMedium::PrintHeader(Option_t* /* option */) {
  // Static method to Print media table header.

    cout << setiosflags(ios::left) << setfill(' ') 
         << "Medium specific physics process & cut flags are listed after "
         << "medium definition if appropriate." << endl;  
    cout << setiosflags(ios::left) << setw(27) << "Medium" 
         << setw(19)  << "Material" 
         << setw(6)   << "isvol"
         << setw(5)   << "ifld"
         << setw(5)   << "fldm"
         << setw(7)   << "tmaxfd"
         << setw(7)   << "stemax"
         << setw(7)   << "deemax"
         << setw(6)   << "epsil"
         << setw(6)   << "stmin"
         << endl;
    cout << setiosflags(ios::left) 
         << setfill('-') <<  setw(6) << "-" << setfill(' ') << setw(21) << " "
         << setfill('-') <<  setw(8) << "-" << setfill(' ') << setw(11) << " "
         << setfill('-') <<  setw(5) << "-" << setfill(' ') << setw(1) << " "
         << setfill('-') <<  setw(4) << "-" << setfill(' ') << setw(1) << " "
         << setfill('-') <<  setw(4) << "-" << setfill(' ') << setw(1) << " "
         << setfill('-') <<  setw(6) << "-" << setfill(' ') << setw(1) << " "
         << setfill('-') <<  setw(6) << "-" << setfill(' ') << setw(1) << " "
         << setfill('-') <<  setw(6) << "-" << setfill(' ') << setw(1) << " "
         << setfill('-') <<  setw(5) << "-" << setfill(' ') << setw(1) << " "
         << setfill('-') <<  setw(5) << "-" << setfill(' ') << setw(1) << " "
         << endl;
  
}

  
//_____________________________________________________________________________
void GeoMedium::Print(Option_t* /* option */) const {
  // Print medium.

  cout << setiosflags(ios::left) << setw(27) << GetName() 
       << setw(19) << GetMaterial()->GetName() 
       << setw(6)  << GetParam(0) // isvol
       << setw(5)  << GetParam(1) // ifield
       << setw(5)  << GetParam(2) // fieldm
       << setw(7)  << GetParam(3) // tmaxfd
       << setw(7)  << GetParam(4) // stemax
       << setw(7)  << GetParam(5) // deemax
       << setw(6)  << GetParam(6) // epsil
       << setw(6)  << GetParam(7); // stmin
  
  Int_t nflag = 0;
  const ProcessMap& processmap = GetProcessMap();
  for ( ProcessMapConstItr citr  = processmap.begin(); 
                           citr != processmap.end(); citr++ ) {
    if ( nflag%6 == 0 && nflag != 0 ) cout << "\n" << setw(95) << " ";
    cout << setw(5) << UtilMCFlag::AsString(citr->first) 
         << setw(2) << citr->second;
    nflag++;
  }

  const CutMap& cutmap = GetCutMap();
  for ( CutMapConstItr citr = cutmap.begin(); citr != cutmap.end(); citr++ ) {
    if ( nflag%3 == 0 && nflag !=0 ) cout << "\n" << setw(95) << " ";
    cout << setw(7) << UtilMCFlag::AsString(citr->first) 
         << setw(7) << citr->second;
    nflag++;
  }
  cout << endl;
  
}

//_____________________________________________________________________________
Double_t GeoMedium::GetCut(UtilMCFlag::ECut cutflg) const {
  // Get medium specific cut value corresponding to cutflg.  
  // Return -1 if not in map.
  
  CutMapConstItr citr = fCutMap.find(cutflg);
  if ( citr != fCutMap.end() ) return citr -> second;
  else return -1;
  
}

//_____________________________________________________________________________
Int_t GeoMedium::GetProcess(UtilMCFlag::EProcess processflg) const {
  // Get medium specific process value corresponding to processflg.
  // Return -1 if not in map. 
  
  ProcessMapConstItr citr = fProcessMap.find(processflg);
  if ( citr != fProcessMap.end() ) return citr -> second;
  else return -1;

}

//_____________________________________________________________________________
Double_t GeoMedium::GetTracking(UtilMCFlag::ETracking trackingflg) const {
  // Get medium specific tracking value corresponding to trackingflg.
  // Return -1 if unknown.

  if ( trackingflg == UtilMCFlag::kUnknownTracking ) {
    MSG("Geo",Msg::kWarning) 
       << "GeoMedium::GetTracking called with trackingflg "
       << UtilMCFlag::AsString(trackingflg) << "." << endl;
    return -1;
  }
  
  return GetParam(trackingflg);

}

//_____________________________________________________________________________
void GeoMedium::SetTracking(UtilMCFlag::ETracking trackingflg, 
                            Double_t tracking) {
  // Set tracking value associated with trackingflg

  if ( trackingflg == UtilMCFlag::kUnknownTracking ) {
    MSG("Geo",Msg::kWarning)
      << "GeoMedium::SetTracking called with trackingflg "
      << UtilMCFlag::AsString(trackingflg) << "." << endl;
    return;
  }
  if ( GetParam(trackingflg) != tracking ) {
    MSG("Geo",Msg::kDebug)
      << "Medium tracking parameter " << UtilMCFlag::AsString(trackingflg) 
      << " set to " << tracking << " for medium " << GetName()
      << " overriding existing " << GetParam(trackingflg) << "." << endl;    
  }
  
  SetParam(trackingflg,tracking);
  
}

---