G4UA::LengthIntegrator Class Reference

A user action used to evaluate thickness of all detectors traversed by outgoing particles. More...

#include <LengthIntegrator.h>

Inheritance diagram for G4UA::LengthIntegrator:

Collaboration diagram for G4UA::LengthIntegrator:

Public Member Functions
	LengthIntegrator (const std::string &histSvcName, bool doHistos)
	Constructor takes the name of the histogram service as argument.
	LengthIntegrator (const LengthIntegrator &)=delete
LengthIntegrator &	operator= (const LengthIntegrator &)=delete
virtual void	BeginOfEventAction (const G4Event *) override
	Called at beginning of G4 event to cache some details about the current primary vertex and particle.
virtual void	EndOfEventAction (const G4Event *) override
	Called at end of G4 event to finalize measurements and fill hists.
virtual void	UserSteppingAction (const G4Step *) override
	Called at every particle step to accumulate thickness.

Private Member Functions
void	fillNtuple ()
std::string	getMaterialClassification (const std::string &name)
std::string	getVolumeType (const std::string &s)
void	addToDetThickMap (const std::string &, double, double)
void	regAndFillHist (const std::string &, const std::pair< double, double > &)
	Setup one set of measurement hists for a detector name.
TProfile2D *	getOrCreateProfile (const std::string &regName, const TString &histoname, const TString &xtitle, int nbinsx, float xmin, float xmax, const TString &ytitle, int nbinsy, float ymin, float ymax, const TString &ztitle)
	this method checks if a histo is on THsvc already and caches a local pointer to it if the histo is not present, it creates and registers it

Private Attributes
G4Pow *	m_g4pow
TTree *	m_tree
int	m_genNPart = 0
float	m_genEta = 0.0F
float	m_genPhi = 0.0F
float	m_genZ = 0.0F
float	m_genR = 0.0F
float	m_total_X0 = 0.0F
float	m_total_L0 = 0.0F
std::vector< double >	m_collected_X0
std::vector< double >	m_collected_L0
std::vector< float >	m_collected_hitr
std::vector< float >	m_collected_hitx
std::vector< float >	m_collected_hity
std::vector< float >	m_collected_hitz
std::vector< float >	m_collected_outhitr
std::vector< float >	m_collected_outhitx
std::vector< float >	m_collected_outhity
std::vector< float >	m_collected_outhitz
std::vector< float >	m_collected_density
std::vector< std::string >	m_collected_material
std::vector< std::string >	m_collected_volume
std::vector< std::string >	m_collected_groupedmaterial
std::vector< std::string >	m_collected_volumetype
bool	m_splitModerator
bool	m_splitPP1
ServiceHandle< ITHistSvc >	m_hSvc
	Handle to the histogram service.
bool	m_doHistos
double	m_etaPrimary
	Cached eta of the current primary.
double	m_phiPrimary
	Cached phi of the current primary.
std::map< std::string, std::pair< double, double > >	m_detThickMap
	Map of detector thickness measurements for current event.
TProfile2D *	m_rzProfRL
	Rad-length profile hist in R-Z.
std::map< std::string, TProfile * >	m_etaMapRL
	Rad-length profile hist in eta.
std::map< std::string, TProfile * >	m_phiMapRL
	Rad-length profile hist in phi.
TProfile2D *	m_rzProfIL
	Int-length profile hist in R-Z.
std::map< std::string, TProfile * >	m_etaMapIL
	Int-length profile hist in eta.
std::map< std::string, TProfile * >	m_phiMapIL
	Int-length profile hist in phi.
std::map< std::string, TProfile2D *, std::less< std::string > >	m_rzMapRL
std::map< std::string, TProfile2D *, std::less< std::string > >	m_xyMapRL
std::map< std::string, TProfile2D *, std::less< std::string > >	m_rzMapIL
std::map< std::string, TProfile2D *, std::less< std::string > >	m_xyMapIL

Detailed Description

A user action used to evaluate thickness of all detectors traversed by outgoing particles.

This user action is currently used only in special runs with geantinos. Thickness is recorded in terms of both rad length and int length.

NOTE: the current design is safe for multi-threading, but not performant due to sharing of the histograms and excessive locking. If this action needs to be used in multi-threaded jobs, we can rewrite it so that each instance has its own copy of the histograms which get merged in finalization of the LengthIntegratorTool.

Definition at line 39 of file LengthIntegrator.h.

Constructor & Destructor Documentation

LengthIntegrator() [1/2]

G4UA::LengthIntegrator::LengthIntegrator	(	const std::string &	histSvcName,
		bool	doHistos = false )

Constructor takes the name of the histogram service as argument.

Definition at line 101 of file LengthIntegrator.cxx.

    : m_g4pow(0),
      m_hSvc(histSvcName, "LengthIntegrator"),
      m_doHistos(doHistos),
      m_etaPrimary(0), m_phiPrimary(0),
      m_rzProfRL(nullptr), m_rzProfIL(nullptr)
  {
    // Protect concurrent access to the non-thread-safe hist svc
    std::lock_guard<std::mutex> lock(gHistSvcMutex);
 
 
    if(m_doHistos){
      // Register the RZ profiles. The other profiles need to wait until the end
      // of the events as they are created only if used.
      const std::string radName = "/lengths/radLen/RZRadLen";
      if(!getHist(m_hSvc, radName, m_rzProfRL)) {
    m_rzProfRL = new TProfile2D("RZRadLen","RZRadLen",1000,-25000.,25000.,2000,0.,15000.);
    regHist(m_hSvc, radName, m_rzProfRL);
      }
      
      const std::string intName = "/lengths/intLen/RZIntLen";
      if(!getHist(m_hSvc, intName, m_rzProfIL)) {
    m_rzProfIL = new TProfile2D("RZIntLen","RZIntLen",1000,-25000.,25000.,2000,0.,15000.);
    regHist(m_hSvc, intName, m_rzProfIL);
      }
      
    }
    m_tree = new TTree( "TheLarch", "And now, the Larch" );
    //Add Braches to the tree
    //Particle properties
    m_tree->Branch("genNPart",                  &m_genNPart,      "genNPart/I");
    m_tree->Branch("genEta",                    &m_genEta,        "genEta/F");
    m_tree->Branch("genPhi",                    &m_genPhi,        "genPhi/F");
    m_tree->Branch("genZ",                      &m_genZ,          "genZ/F");
    m_tree->Branch("genR",                      &m_genR,          "genR/F");
    m_tree->Branch("total_X0",                  &m_total_X0,      "total_X0/F");
    m_tree->Branch("total_L0",                  &m_total_L0,      "total_L0/F");
 
//    m_tree->Branch("collected_sensitivehit",    &m_collected_sensitivehit); //Vector    
 
    m_tree->Branch("collected_X0",              &m_collected_X0); //Vector
    m_tree->Branch("collected_L0",              &m_collected_L0); //Vector
 
    m_tree->Branch("collected_inhitr",            &m_collected_hitr); //Vector
    m_tree->Branch("collected_inhitz",            &m_collected_hitz); //Vector
    m_tree->Branch("collected_inhitx",            &m_collected_hitx); //Vector
    m_tree->Branch("collected_inhity",            &m_collected_hity); //Vector
 
    m_tree->Branch("collected_outhitr",          &m_collected_outhitr); //Vector
    m_tree->Branch("collected_outhitz",          &m_collected_outhitz); //Vector
    m_tree->Branch("collected_outhitx",          &m_collected_outhitx); //Vector
    m_tree->Branch("collected_outhity",          &m_collected_outhity); //Vector
 
    m_tree->Branch("collected_material",        &m_collected_material); //Vector
    m_tree->Branch("collected_density",         &m_collected_density); //Vector
    m_tree->Branch("collected_volume",          &m_collected_volume); //Vector
    m_tree->Branch("collected_groupedmaterial", &m_collected_groupedmaterial); //Vector
    m_tree->Branch("collected_volumetype",      &m_collected_volumetype); //Vector
 
    regTree(m_hSvc, "/lengths/TheLarch", m_tree );
 
    //find a good way to configure this?
    m_splitModerator = true;
    m_splitPP1 = true;
 
 
    m_g4pow = G4Pow::GetInstance();
 
  }

LengthIntegrator() [2/2]

G4UA::LengthIntegrator::LengthIntegrator ( const LengthIntegrator & )

delete

Member Function Documentation

addToDetThickMap()

void G4UA::LengthIntegrator::addToDetThickMap ( const std::string & name, double thickstepRL, double thickstepIL )

private

Definition at line 511 of file LengthIntegrator.cxx.

  {
    auto it=m_detThickMap.find(name);
    if(it!=m_detThickMap.end()){
      (*it).second.first+=thickstepRL;
      (*it).second.second+=thickstepIL;
    } else {
      m_detThickMap.insert(std::map<std::string,std::pair<double,double>,std::less<std::string> >::value_type( name, std::pair<double,double>( thickstepRL, thickstepIL) ) );
    }
  }

BeginOfEventAction()

void G4UA::LengthIntegrator::BeginOfEventAction ( const G4Event * event )

overridevirtual

Called at beginning of G4 event to cache some details about the current primary vertex and particle.

Also resets some measurements.

Definition at line 174 of file LengthIntegrator.cxx.

  {
    m_detThickMap.clear();
    G4PrimaryVertex* vert = event->GetPrimaryVertex(0);
    G4PrimaryParticle* part = vert->GetPrimary();
    G4ThreeVector mom = part->GetMomentum();
    m_etaPrimary = mom.eta();
    m_phiPrimary = mom.phi();
 
    m_genEta = std::isfinite(mom.eta()) ? mom.eta() : -999;
    m_genPhi = mom.phi();
    m_genZ = vert->GetZ0();
    m_genR = sqrt((vert->GetX0()*vert->GetX0())+(vert->GetY0()*vert->GetY0()));
    m_genNPart++;
 
    m_total_X0 = 0;
    m_total_L0 = 0;
 
 
  }

EndOfEventAction()

void G4UA::LengthIntegrator::EndOfEventAction ( const G4Event * )

overridevirtual

Called at end of G4 event to finalize measurements and fill hists.

Definition at line 198 of file LengthIntegrator.cxx.

  {
    // Lazily protect this whole code from concurrent access
    std::lock_guard<std::mutex> lock(gHistSvcMutex);
    
    if(m_doHistos){
 
      // Loop over volumes
      for (const auto& it : m_detThickMap) {
    
    // If histos already exist, then fill them
    if (m_etaMapRL.find(it.first) != m_etaMapRL.end()) {
      m_etaMapRL[it.first]->Fill(m_etaPrimary, it.second.first, 1.);
      m_phiMapRL[it.first]->Fill(m_phiPrimary, it.second.first, 1.);
      
      m_etaMapIL[it.first]->Fill(m_etaPrimary, it.second.second, 1.);
      m_phiMapIL[it.first]->Fill(m_phiPrimary, it.second.second, 1.);
    }
    // New detector volume; register it
    else {
      regAndFillHist(it.first, it.second);
    }
    
      } // Loop over detectors
      
    }
    
    fillNtuple(); 
    
 }

fillNtuple()

void G4UA::LengthIntegrator::fillNtuple ( )

private

Definition at line 386 of file LengthIntegrator.cxx.

                                   {
      
    m_tree->Fill();
 
      
      //Clean vectors and such
      m_collected_X0.clear();
      m_collected_L0.clear();
      
      m_collected_hitr.clear();
      m_collected_hitx.clear();
      m_collected_hity.clear();
      m_collected_hitz.clear();
 
      m_collected_outhitr.clear();
      m_collected_outhitx.clear();
      m_collected_outhity.clear();
      m_collected_outhitz.clear();
 
      m_collected_material.clear();
      m_collected_density.clear();
      m_collected_volume.clear();
      
      m_collected_groupedmaterial.clear();
      m_collected_volumetype.clear();
 
  }

getMaterialClassification()

std::string G4UA::LengthIntegrator::getMaterialClassification ( const std::string & name )

private

Definition at line 230 of file LengthIntegrator.cxx.

  {
 
    if((name.find("DM_Atlas_Air") != std::string::npos) || (name.find("DM_Atlas") != std::string::npos) ||
       (name.find("pix::HEX") != std::string::npos)){
      return "NONE";
    }
 
    if((name.find("Silicon") != std::string::npos) || (name.find("SiMetal") != std::string::npos)){
      return "ActiveSensors";
    }
 
 
    //These material categories and groupings are based on the ITk
    //for other purposes these need to be different
 
    if(name.find("SCT_TiMetal_heat") != std::string::npos) return "HeatExchangers";
    if(name.find("pix::HEX") != std::string::npos) return "HeatExchangers";
    if(name.find("HeatExchangers") != std::string::npos) return "PP1"; // comment out this line to see the Heat Exchangers separately in the final plot
    
    //----HYBRID----:
    if(name.find("matEC_Hybrid") != std::string::npos) return "StripChips"; //"Hybrid";
    else if(name.find("matB_HybridPCB") != std::string::npos) return "StripChips"; //"Hybrid";
    //----SERVICES----:
    else if(name.find("matSV_Endcap") != std::string::npos) return "Services";
    else if(name.find("matSV_Barrel") != std::string::npos) return "Services";
    //----ADHESIVE----
      else if(name.find("BoronNitrideEpoxy") != std::string::npos) return "SupportStructure"; //"Adhesive";
      else if(name.find("SE4445") != std::string::npos) return "SupportStructure"; //"Adhesive";
      //----SUPPORT STRUCTURE----
      else if(name.find("Peek") != std::string::npos) return "SupportStructure";
      else if(name.find("CFRP") != std::string::npos) return "SupportStructure";
      else if(name.find("CFoam") != std::string::npos) return "SupportStructure";
      else if(name.find("K13C2U") != std::string::npos) return "SupportStructure";
      else if(name.find("K13D2U") != std::string::npos) return "SupportStructure";
      else if(name.find("Rohacell110A") != std::string::npos) return "SupportStructure";
      //----SUPPORT HONEYCOMB----
      else if(name.find("Honeycomb10pcf") != std::string::npos) return "SupportStructure"; //"HoneyComb";
      else if(name.find("Honeycomb2pcf") != std::string::npos) return "SupportStructure"; //"HoneyComb";
      else if(name.find("Honeycomb3pcf") != std::string::npos) return "SupportStructure"; //"HoneyComb";
      //----Copper----
      else if(name.find("CuMetal") != std::string::npos) return "Services"; //"Copper";
      else if(name.find("Copper") != std::string::npos) return "Services"; //"Copper";
      //----Active Sensors---- (includes sensors, ABC and HCC chips)
      else if(name.find("SiMetal") != std::string::npos) return "ActiveSensors";
      //----Air----
      else if(name.find("Air") != std::string::npos) return "Air";
      else if(name.find("N2") != std::string::npos) return "Air";
      //----Cooling----
      else if(name.find("CO2Liquid") != std::string::npos) return "Cooling";
      else if(name.find("k9Allcomp") != std::string::npos) return "Cooling";
      else if(name.find("TiMetal") != std::string::npos) return "Cooling";
      //----BusTape----
      else if(name.find("Kapton") != std::string::npos) return "Services"; //"BusTape";
      else if(name.find("matPetalBusKapton") != std::string::npos) return "Services"; //"BusTape";
      //----closeouts and connectors----
      else if(name.find("T300CF") != std::string::npos) return "Cooling"; //"CloseoutsAndConnectors";
      else if(name.find("Torlon") != std::string::npos) return "Cooling"; //"CloseoutsAndConnectors";
      //----strip chips----
      else if(name.find("matDCDC") != std::string::npos) return "StripChips";
      else if(name.find("matEOS") != std::string::npos) return "StripChips";
      //else return "OtherSCT";
 
      // Older stuff
      if(name.find("CO2") != std::string::npos) return "Cooling";
      //if(name.find("BoronNitrideEpoxy") != std::string::npos) return "Cooling";
      if(name.find("BoronNitrideEpoxy") != std::string::npos) return "SupportStructure";
      if(name.find("FwdSpineOutBase") != std::string::npos) return "SupportStructure";
      if(name.find("Rohacell") != std::string::npos) return "SupportStructure";
      if(name.find("Honeycomb") != std::string::npos) return "SupportStructure";
      if(name.find("matSV") != std::string::npos) return "Services";
      if(name.find("matEOS") != std::string::npos) return "Services";
      if(name.find("matDCDC") != std::string::npos) return "Services";
      if(name.find("PCB") != std::string::npos) return "Services";
 
      if(name.find("N2") != std::string::npos) return "DryNitrogen";
      if(name.find("TiMetal") != std::string::npos) return "Services";
      if(name.find("CuMetal") != std::string::npos) return "Services";
      if(name.find("Aluminium") != std::string::npos) return "BeamPipe";
 
    // Do you want to split cooling and services where possible?
    if(name.find("Cooling") != std::string::npos) return "Services";
 
    if(name.find("Inconel") != std::string::npos) return "BeamPipe";
    if(name.find("Aerogel") != std::string::npos) return "BeamPipe";
    if(name.find("Beryllium") != std::string::npos) return "BeamPipe";
    if(name.find("Getter") != std::string::npos) return "BeamPipe";
    if(name.find("Vacuum") != std::string::npos) return "BeamPipe";
    
    
    if(name.find("Iron") != std::string::npos) return "SupportStructure";
    if(name.find("Peek") != std::string::npos) return "SupportStructure";
    if(name.find("CFRP") != std::string::npos) return "SupportStructure";
    if(name.find("CFoam") != std::string::npos) return "SupportStructure";
    if(name.find("K13D2U") != std::string::npos) return "SupportStructure";
    if(name.find("BoratedPolyethylene") != std::string::npos) return "Moderator";
      
    if(name.find("Alpine") != std::string::npos) return "SupportStructure";
    
    if(name.find("Support") != std::string::npos) return "SupportStructure";
    if(name.find("Carbon") != std::string::npos) return "SupportStructure";
    if(name.find("Default") != std::string::npos) return "SupportStructure";
    if(name.find("Moderator") != std::string::npos) return "Moderator";
    if(name.find("Steel") != std::string::npos) return "Steel";
    if(name.find("BarrelStrip") != std::string::npos) return "Services";
    if(name.find("Brl") != std::string::npos) return "Services";
    if(name.find("Svc") != std::string::npos) return "Services";
    if(name.find("InnerIST") != std::string::npos) return "Services";
    if(name.find("InnerPST") != std::string::npos) return "Services";
    if(name.find("BarrelPixel") != std::string::npos) return "Services";
    if(name.find("EndcapPixel") != std::string::npos) return "Services";
    if(name.find("InnerPixel") != std::string::npos) return "Services";
    if(name.find("OuterPixel") != std::string::npos) return "Services";
    if(name.find("pix::Chip") != std::string::npos) return "PixelChips";
    if(name.find("pix::Hybrid") != std::string::npos) return "PixelChips";
    if(name.find("PP0") != std::string::npos) return "PP0";
    if(name.find("PP1") != std::string::npos) return "PP1";
    
    if(name.find("PST") != std::string::npos) return "SupportStructure";
    if(name.find("IST") != std::string::npos) return "SupportStructure";
    if(name.find("Silicon") != std::string::npos) return "ActiveSensors";
    if(name.find("Straw") != std::string::npos) return "ActiveSensors";
    if(name.find("Diamond") != std::string::npos) return "ActiveSensors";
    if(name.find("SiMetal") != std::string::npos) return "ActiveSensors";
    if(name.find("Air") != std::string::npos) return "Air";
    
    //FIXME hack while we fix the air->nitrogen issue
    if(name.find("DryNitrogen") != std::string::npos) return "Air";
 
    if(name.find("CurlyPipeMountain") != std::string::npos) return "SupportStructure";
    if(name.find("Flange") != std::string::npos) return "SupportStructure";
 
    if(name.find("Pigtail") != std::string::npos) return "Services"; //Overwritten by some other statement...
   
    return "NONE";
 
  }

getOrCreateProfile()

TProfile2D * G4UA::LengthIntegrator::getOrCreateProfile ( const std::string & regName, const TString & histoname, const TString & xtitle, int nbinsx, float xmin, float xmax, const TString & ytitle, int nbinsy, float ymin, float ymax, const TString & ztitle )

private

this method checks if a histo is on THsvc already and caches a local pointer to it if the histo is not present, it creates and registers it

note that this should be called from a section protected by a mutex, since it talks to the THitSvc

Definition at line 484 of file LengthIntegrator.cxx.

                                                                                                                                                                                                                                              {
 
    if(m_hSvc->exists(regName)){
      TH2* histo=nullptr;
      if(m_hSvc->getHist(regName,histo).isSuccess())
    return dynamic_cast<TProfile2D*>(histo);
    } else {
      TProfile2D* result= new TProfile2D(histoname,histoname,nbinsx,xmin,xmax,nbinsy,ymin,ymax);
      result->GetXaxis()->SetTitle(xtitle);
      result->GetYaxis()->SetTitle(ytitle);
      result->GetZaxis()->SetTitle(ztitle);
 
      if (m_hSvc && m_hSvc->regHist(regName,result).isFailure()){
    //ATH_MSG_FATAL( "Registration of histogram " << rznameReg << " failed" );
    throw GaudiException("Registration of histogram " + regName + " failed", "RegHistErr", StatusCode::FAILURE);
      }
      return result;
    }
 
    // should never be here
    G4cout<<"ERROR something went wrong in handling of THistSvc "<<regName <<" "<<histoname<<G4endl;
    return nullptr;
  }

getVolumeType()

std::string G4UA::LengthIntegrator::getVolumeType ( const std::string & s )

private

Definition at line 368 of file LengthIntegrator.cxx.

                                                             {
 
    std::string type = "";
    if(m_splitModerator && ( s.find("Moderator") != std::string::npos || s.find("BoratedPolyethylene") != std::string::npos)) type = "Moderator";
    else if(m_splitPP1 && ( s.find("PP0") != std::string::npos || s.find("PP1") != std::string::npos)) type = "PP1";
    else{
      auto colsPos = s.find("::");
      if (colsPos != std::string::npos)
    type = s.substr(0, colsPos);
      else
    type=s;
    }
 
    return type;
  }

operator=()

LengthIntegrator & G4UA::LengthIntegrator::operator= ( const LengthIntegrator & )

delete

regAndFillHist()

void G4UA::LengthIntegrator::regAndFillHist ( const std::string & detName, const std::pair< double, double > & thicks )

private

Setup one set of measurement hists for a detector name.

Definition at line 525 of file LengthIntegrator.cxx.

  {
    TProfile* profEtaRL = nullptr;
    TProfile* profEtaIL = nullptr;
    TProfile* profPhiRL = nullptr;
    TProfile* profPhiIL = nullptr;
 
    auto pathEtaRL = "/lengths/radLen/" + detName + "_RL";
    auto pathEtaIL = "/lengths/intLen/" + detName + "_IL";
    auto pathPhiRL = "/lengths/radLen/" + detName + "Phi_RL";
    auto pathPhiIL = "/lengths/intLen/" + detName + "Phi_IL";
 
    // Eta rad profile
    if(!getHist(m_hSvc, pathEtaRL, profEtaRL)) {
      const std::string name(detName+"_RL");
      profEtaRL = new TProfile(name.c_str(), name.c_str(), 1000, -6., 6.);
      profEtaRL->GetXaxis()->SetTitle("#eta");
      profEtaRL->GetYaxis()->SetTitle("%X0");
      regHist(m_hSvc, pathEtaRL, profEtaRL);
    }
    // Eta int profile
    if(!getHist(m_hSvc, pathEtaIL, profEtaIL)) {
      const std::string name(detName+"_IL");
      profEtaIL = new TProfile(name.c_str(), name.c_str(), 1000, -6., 6.);
      profEtaIL->GetXaxis()->SetTitle("#eta");
      profEtaIL->GetYaxis()->SetTitle("#lambda");
      regHist(m_hSvc, pathEtaIL, profEtaIL);
    }
    // Phi rad profile
    if(!getHist(m_hSvc, pathPhiRL, profPhiRL)) {
      const std::string name(detName+"Phi_RL");
      profPhiRL = new TProfile(name.c_str(), name.c_str(), 500, -M_PI, M_PI);
      profPhiRL->GetXaxis()->SetTitle("#phi");
      profPhiRL->GetYaxis()->SetTitle("%X0");
      regHist(m_hSvc, pathPhiRL, profPhiRL);
    }
    // Phi int profile
    if(!getHist(m_hSvc, pathPhiIL, profPhiIL)) {
      const std::string name(detName+"Phi_IL");
      profPhiIL = new TProfile(name.c_str(), name.c_str(), 500, -M_PI, M_PI);
      profPhiIL->GetXaxis()->SetTitle("#phi");
      profPhiIL->GetYaxis()->SetTitle("#lambda");
      regHist(m_hSvc, pathPhiIL, profPhiIL);
    }
 
    m_etaMapRL[detName] = profEtaRL;
    m_etaMapIL[detName] = profEtaIL;
    m_phiMapRL[detName] = profPhiRL;
    m_phiMapIL[detName] = profPhiIL;
 
    profEtaRL->Fill(m_etaPrimary, thicks.first, 1.);
    profEtaIL->Fill(m_etaPrimary, thicks.second, 1.);
    profPhiRL->Fill(m_phiPrimary, thicks.first, 1.);
    profPhiIL->Fill(m_phiPrimary, thicks.second, 1.);
  }

UserSteppingAction()

void G4UA::LengthIntegrator::UserSteppingAction ( const G4Step * aStep )

overridevirtual

Called at every particle step to accumulate thickness.

Definition at line 418 of file LengthIntegrator.cxx.

  {
 
    const G4TouchableHistory* touchHist =
      static_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
    G4LogicalVolume* lv = touchHist->GetVolume()->GetLogicalVolume();
    G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
    G4ThreeVector endPoint = aStep->GetPostStepPoint()->GetPosition();
    G4Material* mat = lv->GetMaterial();
 
    const std::string& volName = lv->GetName();
    const std::string& matName = mat->GetName();
 
    double radl = mat->GetRadlen();
    double intl = mat->GetNuclearInterLength();
    double stepl = aStep->GetStepLength();
    double density = mat->GetDensity()*CLHEP::cm3/CLHEP::gram;
 
    // FIXME: using DBL_MAX just ensures double overflow below.
    double thickstepRL = radl != 0 ? stepl/radl : DBL_MAX;
    double thickstepIL = intl != 0 ? stepl/intl : DBL_MAX;
 
    //Fill information for the step
    m_collected_X0.push_back(thickstepRL);
    m_collected_L0.push_back(thickstepIL);
    
    m_collected_material.push_back(matName);
    m_collected_density.push_back(density);
    m_collected_volume.push_back(volName);
 
    m_total_X0+=thickstepRL;
    m_total_L0+=thickstepIL;
 
    m_collected_hitr.push_back(hitPoint.perp());
    m_collected_hitz.push_back(hitPoint.z());
    m_collected_hitx.push_back(hitPoint.x());
    m_collected_hity.push_back(hitPoint.y());
 
    m_collected_outhitr.push_back(endPoint.perp());
    m_collected_outhitz.push_back(endPoint.z());
    m_collected_outhitx.push_back(endPoint.x());
    m_collected_outhity.push_back(endPoint.y());
 
    std::string groupmaterial = getMaterialClassification(matName);
 
    std::string volumetype = getVolumeType(matName);
 
    m_collected_groupedmaterial.push_back(std::move(groupmaterial));
    m_collected_volumetype.push_back(std::move(volumetype));
 
    if(m_doHistos){
      // Protect concurrent histo filling
      {
    static std::mutex mutex_instance;
    std::lock_guard<std::mutex> lock(mutex_instance);
    m_rzProfRL->Fill( hitPoint.z() , hitPoint.perp() , thickstepRL , 1. );
    m_rzProfIL->Fill( hitPoint.z() , hitPoint.perp() , thickstepIL , 1. );
    m_rzProfRL->Fill( endPoint.z() , endPoint.perp() , thickstepRL , 1. );
    m_rzProfIL->Fill( endPoint.z() , endPoint.perp() , thickstepIL , 1. );
      }
    }
    
  }

Member Data Documentation

m_collected_density

std::vector<float> G4UA::LengthIntegrator::m_collected_density

private

Definition at line 93 of file LengthIntegrator.h.

m_collected_groupedmaterial

std::vector<std::string> G4UA::LengthIntegrator::m_collected_groupedmaterial

private

Definition at line 97 of file LengthIntegrator.h.

m_collected_hitr

std::vector<float> G4UA::LengthIntegrator::m_collected_hitr

private

Definition at line 83 of file LengthIntegrator.h.

m_collected_hitx

std::vector<float> G4UA::LengthIntegrator::m_collected_hitx

private

Definition at line 84 of file LengthIntegrator.h.

m_collected_hity

std::vector<float> G4UA::LengthIntegrator::m_collected_hity

private

Definition at line 85 of file LengthIntegrator.h.

m_collected_hitz

std::vector<float> G4UA::LengthIntegrator::m_collected_hitz

private

Definition at line 86 of file LengthIntegrator.h.

m_collected_L0

std::vector<double> G4UA::LengthIntegrator::m_collected_L0

private

Definition at line 81 of file LengthIntegrator.h.

m_collected_material

std::vector<std::string> G4UA::LengthIntegrator::m_collected_material

private

Definition at line 94 of file LengthIntegrator.h.

m_collected_outhitr

std::vector<float> G4UA::LengthIntegrator::m_collected_outhitr

private

Definition at line 88 of file LengthIntegrator.h.

m_collected_outhitx

std::vector<float> G4UA::LengthIntegrator::m_collected_outhitx

private

Definition at line 89 of file LengthIntegrator.h.

m_collected_outhity

std::vector<float> G4UA::LengthIntegrator::m_collected_outhity

private

Definition at line 90 of file LengthIntegrator.h.

m_collected_outhitz

std::vector<float> G4UA::LengthIntegrator::m_collected_outhitz

private

Definition at line 91 of file LengthIntegrator.h.

m_collected_volume

std::vector<std::string> G4UA::LengthIntegrator::m_collected_volume

private

Definition at line 95 of file LengthIntegrator.h.

m_collected_volumetype

std::vector<std::string> G4UA::LengthIntegrator::m_collected_volumetype

private

Definition at line 98 of file LengthIntegrator.h.

m_collected_X0

std::vector<double> G4UA::LengthIntegrator::m_collected_X0

private

Definition at line 80 of file LengthIntegrator.h.

m_detThickMap

std::map<std::string, std::pair<double, double> > G4UA::LengthIntegrator::m_detThickMap

private

Map of detector thickness measurements for current event.

Definition at line 129 of file LengthIntegrator.h.

m_doHistos

bool G4UA::LengthIntegrator::m_doHistos

private

Definition at line 121 of file LengthIntegrator.h.

m_etaMapIL

std::map<std::string, TProfile*> G4UA::LengthIntegrator::m_etaMapIL

private

Int-length profile hist in eta.

Definition at line 143 of file LengthIntegrator.h.

m_etaMapRL

std::map<std::string, TProfile*> G4UA::LengthIntegrator::m_etaMapRL

private

Rad-length profile hist in eta.

Definition at line 134 of file LengthIntegrator.h.

m_etaPrimary

double G4UA::LengthIntegrator::m_etaPrimary

private

Cached eta of the current primary.

Definition at line 124 of file LengthIntegrator.h.

m_g4pow

G4Pow* G4UA::LengthIntegrator::m_g4pow

private

Definition at line 64 of file LengthIntegrator.h.

m_genEta

float G4UA::LengthIntegrator::m_genEta = 0.0F

private

Definition at line 71 of file LengthIntegrator.h.

m_genNPart

int G4UA::LengthIntegrator::m_genNPart = 0

private

Definition at line 70 of file LengthIntegrator.h.

m_genPhi

float G4UA::LengthIntegrator::m_genPhi = 0.0F

private

Definition at line 72 of file LengthIntegrator.h.

m_genR

float G4UA::LengthIntegrator::m_genR = 0.0F

private

Definition at line 74 of file LengthIntegrator.h.

m_genZ

float G4UA::LengthIntegrator::m_genZ = 0.0F

private

Definition at line 73 of file LengthIntegrator.h.

m_hSvc

ServiceHandle<ITHistSvc> G4UA::LengthIntegrator::m_hSvc

private

Handle to the histogram service.

Definition at line 118 of file LengthIntegrator.h.

m_phiMapIL

std::map<std::string, TProfile*> G4UA::LengthIntegrator::m_phiMapIL

private

Int-length profile hist in phi.

Definition at line 145 of file LengthIntegrator.h.

m_phiMapRL

std::map<std::string, TProfile*> G4UA::LengthIntegrator::m_phiMapRL

private

Rad-length profile hist in phi.

Definition at line 136 of file LengthIntegrator.h.

m_phiPrimary

double G4UA::LengthIntegrator::m_phiPrimary

private

Cached phi of the current primary.

Definition at line 126 of file LengthIntegrator.h.

m_rzMapIL

std::map<std::string,TProfile2D*,std::less<std::string> > G4UA::LengthIntegrator::m_rzMapIL

private

Definition at line 151 of file LengthIntegrator.h.

m_rzMapRL

std::map<std::string,TProfile2D*,std::less<std::string> > G4UA::LengthIntegrator::m_rzMapRL

private

Definition at line 148 of file LengthIntegrator.h.

m_rzProfIL

TProfile2D* G4UA::LengthIntegrator::m_rzProfIL

private

Int-length profile hist in R-Z.

Definition at line 141 of file LengthIntegrator.h.

m_rzProfRL

TProfile2D* G4UA::LengthIntegrator::m_rzProfRL

private

Rad-length profile hist in R-Z.

Definition at line 132 of file LengthIntegrator.h.

m_splitModerator

bool G4UA::LengthIntegrator::m_splitModerator

private

Definition at line 100 of file LengthIntegrator.h.

m_splitPP1

bool G4UA::LengthIntegrator::m_splitPP1

private

Definition at line 101 of file LengthIntegrator.h.

m_total_L0

float G4UA::LengthIntegrator::m_total_L0 = 0.0F

private

Definition at line 78 of file LengthIntegrator.h.

m_total_X0

float G4UA::LengthIntegrator::m_total_X0 = 0.0F

private

Definition at line 77 of file LengthIntegrator.h.

m_tree

TTree* G4UA::LengthIntegrator::m_tree

private

Definition at line 67 of file LengthIntegrator.h.

m_xyMapIL

std::map<std::string,TProfile2D*,std::less<std::string> > G4UA::LengthIntegrator::m_xyMapIL

private

Definition at line 152 of file LengthIntegrator.h.

m_xyMapRL

std::map<std::string,TProfile2D*,std::less<std::string> > G4UA::LengthIntegrator::m_xyMapRL

private

Definition at line 149 of file LengthIntegrator.h.

---

The documentation for this class was generated from the following files:

• LengthIntegrator.h
• LengthIntegrator.cxx