/home/cern/BDSIM_new/src/BDSDump.cc

//  
//   BDSIM, (C) 2001-2007
//    
//   version 0.3 
//   last modified : 08 May 2007 by agapov@pp.rhul.ac.uk
//  


//
//    beam dumper/reader for online exchange with external codes
//

#include "BDSGlobalConstants.hh" // must be first in include list
#include "BDSDump.hh"
#include "G4Box.hh"
#include "G4Tubs.hh"
#include "G4VisAttributes.hh"
#include "G4LogicalVolume.hh"
#include "G4VPhysicalVolume.hh"
#include "G4PVPlacement.hh"               
#include "G4UserLimits.hh"
#include "BDSOutput.hh"
#include "BDSDumpSD.hh"
#include "G4SDManager.hh"

//#include"MagFieldFunction.hh"
#include <map>



typedef std::map<G4String,int> LogVolCountMap;
extern LogVolCountMap* LogVolCount;

typedef std::map<G4String,G4LogicalVolume*> LogVolMap;
extern LogVolMap* LogVol;
extern BDSOutput bdsOutput;
BDSDumpSD* BDSDumpSensDet;
extern BDSMaterials* theMaterials;
//============================================================

BDSDump::BDSDump (G4String aName,G4double aLength):
  BDSAcceleratorComponent(
                         aName,
                         aLength,0,0,0,
                         SetVisAttributes())
{
  DumpLogicalVolume();
  nDumps++;
  //G4int nDumps=(*LogVolCount)[itsName];

  //BDSRoot->SetDumpNumber(nDumps);

 
}

const int BDSDump::GetNumberOfDumps()
{
  return nDumps;
}

int BDSDump::nDumps=0;

int BDSDump::firstActiveDump=1;

int BDSDump::nUsedDumps=0;

void BDSDump::DumpLogicalVolume()
{
  if(!(*LogVolCount)[itsName])
    {

      G4double SampTransSize;
      SampTransSize=2.*BDSGlobals->GetTunnelRadius();

      itsMarkerLogicalVolume=
        new G4LogicalVolume(
                            new G4Box(itsName+"_solid",
                                      SampTransSize,
                                      SampTransSize,
                                      itsLength/2.0),
                            theMaterials->GetMaterial("Vacuum"),
                            itsName);

      (*LogVolCount)[itsName]=1;
      (*LogVol)[itsName]=itsMarkerLogicalVolume;

      itsOuterUserLimits =new G4UserLimits();
      itsOuterUserLimits->SetMaxAllowedStep(itsLength);
      itsMarkerLogicalVolume->SetUserLimits(itsOuterUserLimits);

      // Sensitive Detector:
      if(true)
        {
          G4SDManager* SDMan = G4SDManager::GetSDMpointer();
          BDSDumpSensDet=new BDSDumpSD(itsName,"plane");
          SDMan->AddNewDetector(BDSDumpSensDet);
          itsMarkerLogicalVolume->SetSensitiveDetector(BDSDumpSensDet);
        }
    }
  else
    {
      (*LogVolCount)[itsName]++;
      itsMarkerLogicalVolume=(*LogVol)[itsName];
      itsMarkerLogicalVolume->SetSensitiveDetector(BDSDumpSensDet);
    }
}

G4VisAttributes* BDSDump::SetVisAttributes()
{
  itsVisAttributes=new G4VisAttributes(G4Colour(1,1,1));
  return itsVisAttributes;
}

BDSDump::~BDSDump()
{
  if(itsVisAttributes) delete itsVisAttributes;
  if(itsUserLimits) delete itsUserLimits;
  nDumps--;
}

---