DDG4ProductionCuts Class Reference

#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 10 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 17 of file DDG4ProductionCuts.cc.

                                       {
}

Member Function Documentation

G4ProductionCuts * DDG4ProductionCuts::getProductionCuts ( G4Region *  region )

private

Definition at line 154 of file DDG4ProductionCuts.cc.

Referenced by setProdCuts().

                                                                         {

  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 149 of file DDG4ProductionCuts.cc.

Referenced by initialize(), and setProdCuts().

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

void DDG4ProductionCuts::initialize ( )

private

Definition at line 55 of file DDG4ProductionCuts.cc.

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

Referenced by DDG4ProductionCuts().

                                    {

  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 SimG4Exception("DDG4ProductionCuts: 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 93 of file DDG4ProductionCuts.cc.

References 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 update().

                                             {  
  
  if ( m_Verbosity > 0 ) 
    LogDebug("Physics") <<" DDG4ProductionCuts: inside setProdCuts";
  
  G4Region * region = 0;
  
  std::string  regionName;
  unsigned int num= map_.toString(m_KeywordRegion,lpart,regionName);
  
  if (num != 1) 
    throw SimG4Exception("DDG4ProductionCuts: 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 SimG4Exception("DDG4ProductionCuts: Problem with Region tags: no/more than one ProdCutsForGamma.");
  }
  temp =  map_.toDouble("ProdCutsForElectrons",lpart,electroncut);
  if (temp != 1){
    throw SimG4Exception("DDG4ProductionCuts: Problem with Region tags: no/more than one ProdCutsForElectrons.");
  }
  temp =  map_.toDouble("ProdCutsForPositrons",lpart,positroncut);
  if (temp != 1) {
    throw SimG4Exception("DDG4ProductionCuts: 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 22 of file DDG4ProductionCuts.h.

References m_Verbosity, and reco::return().

{ m_Verbosity = verb; return ; }

void DDG4ProductionCuts::update

(

void

)

Definition at line 44 of file DDG4ProductionCuts.cc.

References setProdCuts(), and vec_.

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), PhysicsList::SetCuts(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                {
  //
  // 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);
  }
}

Member Data Documentation

std::string DDG4ProductionCuts::m_KeywordRegion

private

Definition at line 32 of file DDG4ProductionCuts.h.

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

bool DDG4ProductionCuts::m_protonCut

private

Definition at line 34 of file DDG4ProductionCuts.h.

Referenced by DDG4ProductionCuts(), and setProdCuts().

int DDG4ProductionCuts::m_Verbosity

private

Definition at line 33 of file DDG4ProductionCuts.h.

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

G4LogicalVolumeToDDLogicalPartMap DDG4ProductionCuts::map_

private

Definition at line 31 of file DDG4ProductionCuts.h.

Referenced by initialize(), and setProdCuts().

G4LogicalVolumeToDDLogicalPartMap::Vector DDG4ProductionCuts::vec_

private

Definition at line 35 of file DDG4ProductionCuts.h.

Referenced by initialize(), and update().