#include <DDG4ProductionCuts.h>

Public Member Functions
	DDG4ProductionCuts (const G4LogicalVolumeToDDLogicalPartMap &, int, const edm::ParameterSet &p)

void	SetVerbosity (int verb)

void	update ()

	~DDG4ProductionCuts ()

Private Member Functions
G4ProductionCuts *	getProductionCuts (G4Region *region)

G4Region *	getRegion (const std::string &region)

void	initialize ()

void	setProdCuts (const DDLogicalPart lpart, G4LogicalVolume *lvolume)

Private Attributes
std::string	m_KeywordRegion

bool	m_protonCut

int	m_Verbosity

G4LogicalVolumeToDDLogicalPartMap	map_

G4LogicalVolumeToDDLogicalPartMap::Vector	vec_

Detailed Description

Definition at line 15 of file DDG4ProductionCuts.h.

Constructor & Destructor Documentation

DDG4ProductionCuts::DDG4ProductionCuts	(	const G4LogicalVolumeToDDLogicalPartMap &	map,
		int	verb,
		const edm::ParameterSet &	p
	)

Definition at line 11 of file DDG4ProductionCuts.cc.

References edm::ParameterSet::getUntrackedParameter(), initialize(), m_KeywordRegion, and m_protonCut.

     : map_(map), m_Verbosity(verb) {
   m_KeywordRegion = "CMSCutsRegion";
   m_protonCut = p.getUntrackedParameter<bool>("CutsOnProton", true);
   initialize();
 }

DDG4ProductionCuts::~DDG4ProductionCuts ( )

Definition at line 20 of file DDG4ProductionCuts.cc.

20 {}

Member Function Documentation

G4ProductionCuts * DDG4ProductionCuts::getProductionCuts ( G4Region * region )

private

Definition at line 152 of file DDG4ProductionCuts.cc.

Referenced by setProdCuts(), and SetVerbosity().

                                                                      {
   G4ProductionCuts* prodCuts = reg->GetProductionCuts();
   if (!prodCuts) {
     prodCuts = new G4ProductionCuts();
     reg->SetProductionCuts(prodCuts);
   }
   return prodCuts;
 }

G4Region * DDG4ProductionCuts::getRegion ( const std::string & region )

private

Definition at line 147 of file DDG4ProductionCuts.cc.

Referenced by initialize(), setProdCuts(), and SetVerbosity().

                                                                 {
   G4Region* reg = G4RegionStore::GetInstance()->FindOrCreateRegion(regName);
   return reg;
 }

void DDG4ProductionCuts::initialize ( void )

private

Definition at line 55 of file DDG4ProductionCuts.cc.

References DDMapper< KeyType, ValueType >::all(), dd_is_greater(), Exception, getRegion(), findQualityFiles::jj, LogDebug, m_KeywordRegion, m_Verbosity, map_, pileupDistInMC::num, AlCaHLTBitMon_QueryRunRegistry::string, DDMapper< KeyType, ValueType >::toString(), and vec_.

Referenced by DDG4ProductionCuts(), and SetVerbosity().

                                     {
   vec_ = map_.all(m_KeywordRegion);
   // sort all root volumes - to get the same sequence at every run of the application.
   // (otherwise, the sequence will depend on the pointer (memory address) of the
   // involved objects, because 'new' does no guarantee that you allways get a
   // higher (or lower) address when allocating an object of the same type ...
   sort(vec_.begin(), vec_.end(), &dd_is_greater);
   if (m_Verbosity > 0) {
     LogDebug("Physics") << " DDG4ProductionCuts (New) : starting\n"
                         << " DDG4ProductionCuts : Got " << vec_.size() << " region roots.\n"
                         << " DDG4ProductionCuts : List of all roots:";
     for (size_t jj = 0; jj < vec_.size(); ++jj)
       LogDebug("Physics") << "   DDG4ProductionCuts : root=" << vec_[jj].second.name();
   }
 
   // Now generate all the regions
   for (G4LogicalVolumeToDDLogicalPartMap::Vector::iterator tit = vec_.begin(); tit != vec_.end(); tit++) {
     std::string regionName;
     unsigned int num = map_.toString(m_KeywordRegion, (*tit).second, regionName);
 
     if (num != 1) {
       throw cms::Exception("SimG4CorePhysics", " DDG4ProductionCuts::initialize: Problem with Region tags.");
     }
     G4Region* region = getRegion(regionName);
     region->AddRootLogicalVolume((*tit).first);
 
     if (m_Verbosity > 0)
       LogDebug("Physics") << " MakeRegions: added " << ((*tit).first)->GetName() << " to region " << region->GetName();
   }
 }

void DDG4ProductionCuts::setProdCuts	(	const DDLogicalPart	lpart,
		G4LogicalVolume *	lvolume
	)

private

Definition at line 86 of file DDG4ProductionCuts.cc.

References Exception, getProductionCuts(), getRegion(), LogDebug, m_KeywordRegion, m_protonCut, m_Verbosity, map_, pileupDistInMC::num, AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, DDMapper< KeyType, ValueType >::toDouble(), and DDMapper< KeyType, ValueType >::toString().

Referenced by SetVerbosity(), and update().

                                                                                      {
   if (m_Verbosity > 0)
     LogDebug("Physics") << " DDG4ProductionCuts: inside setProdCuts";
 
   G4Region* region = nullptr;
 
   std::string regionName;
   unsigned int num = map_.toString(m_KeywordRegion, lpart, regionName);
 
   if (num != 1) {
     throw cms::Exception("SimG4CorePhysics", " DDG4ProductionCuts::setProdCuts: Problem with Region tags.");
   }
   if (m_Verbosity > 0)
     LogDebug("Physics") << "Using region " << regionName;
 
   region = getRegion(regionName);
 
   //
   // search for production cuts
   // you must have four of them: e+ e- gamma proton
   //
   double gammacut;
   double electroncut;
   double positroncut;
   double protoncut = 0.0;
   int temp = map_.toDouble("ProdCutsForGamma", lpart, gammacut);
   if (temp != 1) {
     throw cms::Exception(
         "SimG4CorePhysics",
         " DDG4ProductionCuts::setProdCuts: Problem with Region tags - no/more than one ProdCutsForGamma.");
   }
   temp = map_.toDouble("ProdCutsForElectrons", lpart, electroncut);
   if (temp != 1) {
     throw cms::Exception(
         "SimG4CorePhysics",
         " DDG4ProductionCuts::setProdCuts: Problem with Region tags - no/more than one ProdCutsForElectrons.");
   }
   temp = map_.toDouble("ProdCutsForPositrons", lpart, positroncut);
   if (temp != 1) {
     throw cms::Exception(
         "SimG4CorePhysics",
         " DDG4ProductionCuts::setProdCuts: Problem with Region tags - no/more than one ProdCutsForPositrons.");
   }
   //
   // For the moment I assume all of the three are set
   //
   G4ProductionCuts* prodCuts = getProductionCuts(region);
   prodCuts->SetProductionCut(gammacut, idxG4GammaCut);
   prodCuts->SetProductionCut(electroncut, idxG4ElectronCut);
   prodCuts->SetProductionCut(positroncut, idxG4PositronCut);
   // For recoil use the same cut as for e-
   if (m_protonCut) {
     protoncut = electroncut;
   }
   prodCuts->SetProductionCut(protoncut, idxG4ProtonCut);
   if (m_Verbosity > 0) {
     LogDebug("Physics") << "DDG4ProductionCuts : Setting cuts for " << regionName << "\n    Electrons: " << electroncut
                         << "\n    Positrons: " << positroncut << "\n    Gamma    : " << gammacut;
   }
 }

void DDG4ProductionCuts::SetVerbosity ( int verb )

inline

Definition at line 20 of file DDG4ProductionCuts.h.

References getProductionCuts(), getRegion(), initialize(), m_Verbosity, setProdCuts(), and AlCaHLTBitMon_QueryRunRegistry::string.

                               {
     m_Verbosity = verb;
     return;
   }

void DDG4ProductionCuts::update ( )

Definition at line 46 of file DDG4ProductionCuts.cc.

References setProdCuts(), and vec_.

                                 {
   //
   // Loop over all DDLP and provide the cuts for each region
   //
   for (G4LogicalVolumeToDDLogicalPartMap::Vector::iterator tit = vec_.begin(); tit != vec_.end(); tit++) {
     setProdCuts((*tit).second, (*tit).first);
   }
 }