G3NegVolPars.cc File Reference

#include "globals.hh"
#include "G3VolTable.hh"
#include "G4VSolid.hh"
#include "G3toG4.hh"
#include <cmath>

Go to the source code of this file.

Functions
G4bool	G3CalcParamsFn (G4double *rpar, G4int npar, G4double *rparm, G4String shape, G4String shapem)
G4bool	G3NegVolPars (G4double pars[], G4int *nparpt, G3VolTableEntry *vte, G3VolTableEntry *mvte, const char routine[])

Function Documentation

G4bool G3CalcParamsFn	(	G4double *	rpar,
		G4int	npar,
		G4double *	rparm,
		G4String	shape,
		G4String	shapem
	)

Definition at line 37 of file G3NegVolPars.cc.

References FALSE, G4INCL::Math::min(), and TRUE.

Referenced by G3NegVolPars().

{
  G4bool NegPresent = FALSE;
  // for normal single volume positioning, just substitute for the
  // negative parameters
  // treat only the legal cases
  if (shapem == shape) {
    if (shape == "BOX" || shape == "TRD1" || shape == "TRD2" || 
        shape == "ELTU") {
      for (G4int i=0;i<npar;i++) {
        if (rpar[i] < 0) {
          if (rparm != 0) rpar[i] = rparm[i];
      if (rpar[i] < 0) NegPresent = TRUE;
        }
      }
    }
    if (shape == "TRAP") {
      for (G4int i=0;i<11;i++) {
        if (i != 1 && i != 2 && i != 6 && i != 10) {
          if (rpar[i]<0) {
            if (rparm != 0) rpar[i] = rparm[i];
        if (rpar[i] < 0) NegPresent = TRUE;
      }  
        }
      }
    }
    if (shape == "TUBE" || shape == "TUBS" || shape == "PARA") {
      for (G4int i=0;i<3;i++) {
        if (rpar[i] < 0) {
          if (rparm != 0) rpar[i] = rparm[i];
      if (rpar[i] < 0) NegPresent = TRUE;
        }
      }
    }
    if (shape == "CONE" || shape == "CONS") {
      for (G4int i=0;i<5;i++) {
        if (rpar[i] < 0) {
          if (rparm != 0) rpar[i] = rparm[i];
      if (rpar[i] < 0) NegPresent = TRUE;
        }
      }
    }
    if (shape == "SPHE") {
      for (G4int i=0;i<2;i++) {
        if (rpar[i] < 0) {
          if (rparm != 0) rpar[i] = rparm[i];
      if (rpar[i] < 0) NegPresent = TRUE;
        }
      }
    }
    if (shape == "PGON") {
      G4int nz = G4int(rpar[3]);
      G4int ipl;
      for (G4int i=0;i<nz;i++) {
        ipl = 5 + i*3;
        if (rpar[ipl] < 0) {
          if (rparm != 0) rpar[ipl] = rparm[ipl];
      if (rpar[ipl] < 0) NegPresent = TRUE;
        }
        if (rpar[ipl+1] < 0) {
          if (rparm != 0)  rpar[ipl] = rparm[ipl];
      if (rpar[ipl] < 0) NegPresent = TRUE;
        }
      }
    }
    if (shape == "PCON") {
      G4int nz = G4int(rpar[2]);
      G4int ipl;
      for (G4int i=0;i<nz;i++) {
        ipl = 4 + i*3;
        if (rpar[ipl] < 0) {
          if (rparm != 0) rpar[ipl] = rparm[ipl];
      if (rpar[ipl] < 0) NegPresent = TRUE;
        }
        if (rpar[ipl+1] < 0) {
          // TO DO
      // check - folowing argument might be ipl+1
          if (rparm != 0) rpar[ipl] = rparm[ipl];
      if (rpar[ipl] < 0) NegPresent = TRUE;
        }
      }
    }
  }

  if (shape == "BOX") {
    if (shapem == "TRD1") {
      if (rpar[1] < 0) {
        if (rparm != 0) rpar[1] = rparm[2];
    if (rpar[1] < 0) NegPresent = TRUE;
      }
      if (rpar[2] < 0) {
        if (rparm != 0) rpar[2] = rparm[3];
    if (rpar[2] < 0) NegPresent = TRUE;
      }
      if (rpar[0] < 0) {
        if (rparm != 0) rpar[0] = std::min(rparm[0],rparm[1]) +
                                  std::abs(rparm[0]-rparm[1])*.5*rpar[2]/rparm[3];
    if (rpar[0] < 0) NegPresent = TRUE;
      }
    }
    if (shapem == "TRD2") {
      if (rpar[2] < 0) {
        if (rparm != 0) rpar[2] = rparm[4];
        if (rpar[2]<0) NegPresent = TRUE;
      }
      if (rpar[0] < 0) {
        if (rparm != 0) rpar[0] = std::min(rparm[0],rparm[1]) +
                                  std::abs(rparm[0]-rparm[1])*.5*rpar[2]/rparm[4];
        if (rpar[0]<0) NegPresent = TRUE;
      }
      if (rpar[1] < 0) {
        if (rparm != 0) rpar[1] = std::min(rparm[2],rparm[3]) +
                                  std::abs(rparm[2]-rparm[3])*.5*rpar[2]/rparm[4];
        if (rpar[1]<0) NegPresent = TRUE;
      }
    }
    if (shapem == "TRAP") {
      if (rpar[2] < 0) {
        if (rparm != 0) rpar[2] = rparm[0];
        if (rpar[2] < 0) NegPresent = TRUE;
      }
      if (rpar[0] < 0) {
        if (rparm != 0) {
          G4double xlo = std::min(rparm[4],rparm[8]) +
            std::abs(rparm[4]-rparm[8])*.5*rpar[2]/rparm[0];
          G4double xhi = std::min(rparm[5],rparm[9]) +
            std::abs(rparm[5]-rparm[9])*.5*rpar[2]/rparm[0];
          rpar[0] = std::min(xlo,xhi);
        }
        if (rpar[0] < 0) NegPresent = TRUE;
      }
      if (rpar[1] < 0) {
        if (rparm != 0) rpar[1] = std::min(rparm[3],rparm[7]) +
                        std::abs(rparm[3]-rparm[7])*.5*rpar[2]/rparm[0];
        if (rpar[1] < 0) NegPresent = TRUE;
      }
    }
  }
  return NegPresent;
}

G4bool G3NegVolPars	(	G4double	pars[],
		G4int *	nparpt,
		G3VolTableEntry *	vte,
		G3VolTableEntry *	mvte,
		const char	routine[]
	)

Definition at line 181 of file G3NegVolPars.cc.

References FALSE, G3CalcParamsFn(), G3VolTableEntry::GetNpar(), G3VolTableEntry::GetRpar(), G3VolTableEntry::GetShape(), and TRUE.

Referenced by G4CreateCloneVTE(), and G4ProcessDaughters().

{
  G4bool NegPresent = FALSE;

  // retrieve parameters 

  // the volume
  G4String shape = vte->GetShape();
  G4double* rpar = vte->GetRpar();
  G4int npar = vte->GetNpar();  
  if (npar ==0) {
    // no solid parameters are defined in vte
    npar = *nparpt;
    rpar = pars;
  }
  else {  
    // solid parameters are already defined in vte
    // pars[], nparpt are ignored
    // TO DO: check if g3 ignores them too or resets
    // vte parameters according to this new ones !!
  }  
      
  // mother
  G4String shapem = mvte->GetShape();
  G4double* rparm = mvte->GetRpar();

  if (strcmp(routine,"GSPOS") == 0 || strcmp(routine,"GSVOLU") == 0) {
    NegPresent = G3CalcParamsFn(rpar,npar,rparm,shape,shapem);
  }
  if (strcmp(routine,"GSDVN") == 0) {
    // just set the flag. The parametrization function figures out
    // what to do.
    for (G4int i=0;i<npar;i++) {
      if (rpar[i] < 0) {
        NegPresent = TRUE;
      }
    }
  }
  return NegPresent;
}