StackingAction Class Reference

#include <StackingAction.h>

Inheritance diagram for StackingAction:

Public Member Functions
virtual G4ClassificationOfNewTrack	ClassifyNewTrack (const G4Track *aTrack)
void	NewStage ()
void	PrepareNewEvent ()
	StackingAction (EventAction *e, const edm::ParameterSet &ps)
virtual	~StackingAction ()

Private Member Functions
void	initPointer ()
bool	isItLongLived (const G4Track *) const
int	isItPrimaryDecayProductOrConversion (const G4Track *, const G4Track &) const
bool	isThisVolume (const G4VTouchable *, std::vector< G4LogicalVolume * > &) const
bool	rrApplicable (const G4Track *, const G4Track &) const

Private Attributes
std::vector< G4LogicalVolume * >	calo
EventAction *	eventAction_
bool	gRRactive
double	gRusRoCastor
double	gRusRoEcal
double	gRusRoEnerLim
double	gRusRoHcal
double	gRusRoMuonIron
double	gRusRoPreShower
double	gRusRoWorld
bool	killDeltaRay
bool	killHeavy
bool	killInCalo
bool	killInCaloEfH
double	kmaxIon
double	kmaxNeutron
double	kmaxProton
std::vector< std::string >	maxTimeNames
std::vector< G4Region * >	maxTimeRegions
double	maxTrackTime
std::vector< double >	maxTrackTimes
std::vector< G4LogicalVolume * >	muon
NewTrackAction *	newTA
bool	nRRactive
double	nRusRoCastor
double	nRusRoEcal
double	nRusRoEnerLim
double	nRusRoHcal
double	nRusRoMuonIron
double	nRusRoPreShower
double	nRusRoWorld
bool	pRRactive
double	pRusRoCastor
double	pRusRoEcal
double	pRusRoEnerLim
double	pRusRoHcal
double	pRusRoMuonIron
double	pRusRoPreShower
double	pRusRoWorld
G4Region *	regionCastor
G4Region *	regionEcal
G4Region *	regionHcal
G4Region *	regionMuonIron
G4Region *	regionPreShower
G4Region *	regionWorld
bool	saveFirstSecondary
bool	savePDandCinCalo
bool	savePDandCinMuon
bool	savePDandCinTracker
std::vector< G4LogicalVolume * >	tracker
bool	trackNeutrino

Detailed Description

Definition at line 17 of file StackingAction.h.

Constructor & Destructor Documentation

StackingAction::StackingAction	(	EventAction *	e,
		const edm::ParameterSet &	ps
	)

Definition at line 21 of file StackingAction.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), gRRactive, gRusRoCastor, gRusRoEcal, gRusRoEnerLim, gRusRoHcal, gRusRoMuonIron, gRusRoPreShower, gRusRoWorld, i, initPointer(), killDeltaRay, killHeavy, killInCalo, killInCaloEfH, kmaxIon, kmaxNeutron, kmaxProton, maxTimeNames, maxTrackTime, maxTrackTimes, newTA, nRRactive, nRusRoCastor, nRusRoEcal, nRusRoEnerLim, nRusRoHcal, nRusRoMuonIron, nRusRoPreShower, nRusRoWorld, pRRactive, pRusRoCastor, pRusRoEcal, pRusRoEnerLim, pRusRoHcal, pRusRoMuonIron, pRusRoPreShower, pRusRoWorld, regionCastor, regionEcal, regionHcal, regionMuonIron, regionPreShower, regionWorld, saveFirstSecondary, savePDandCinCalo, savePDandCinMuon, savePDandCinTracker, and trackNeutrino.

  : eventAction_(e)
{
  trackNeutrino  = p.getParameter<bool>("TrackNeutrino");
  killHeavy      = p.getParameter<bool>("KillHeavy");
  kmaxIon        = p.getParameter<double>("IonThreshold")*MeV;
  kmaxProton     = p.getParameter<double>("ProtonThreshold")*MeV;
  kmaxNeutron    = p.getParameter<double>("NeutronThreshold")*MeV;
  killDeltaRay   = p.getParameter<bool>("KillDeltaRay");
  maxTrackTime   = p.getParameter<double>("MaxTrackTime")*ns;
  maxTrackTimes  = p.getParameter<std::vector<double> >("MaxTrackTimes");
  for(unsigned int i=0; i<maxTrackTimes.size(); ++i) { maxTrackTimes[i] *= ns; }
  maxTimeNames   = p.getParameter<std::vector<std::string> >("MaxTimeNames");
  savePDandCinTracker = 
    p.getUntrackedParameter<bool>("SavePrimaryDecayProductsAndConversionsInTracker",
                  false);
  savePDandCinCalo = 
    p.getUntrackedParameter<bool>("SavePrimaryDecayProductsAndConversionsInCalo",
                  false);
  savePDandCinMuon = 
    p.getUntrackedParameter<bool>("SavePrimaryDecayProductsAndConversionsInMuon",
                  false);
  saveFirstSecondary = 
    p.getUntrackedParameter<bool>("SaveFirstLevelSecondary",false);
  killInCalo = false;
  killInCaloEfH = false;

  // Russian Roulette
  regionEcal = 0;
  regionHcal = 0;
  regionMuonIron = 0;
  regionPreShower= 0;
  regionCastor = 0;
  regionWorld = 0;

  gRusRoEnerLim = p.getParameter<double>("RusRoGammaEnergyLimit")*MeV;
  nRusRoEnerLim = p.getParameter<double>("RusRoNeutronEnergyLimit")*MeV;
  pRusRoEnerLim = p.getParameter<double>("RusRoProtonEnergyLimit")*MeV;

  gRusRoEcal = p.getParameter<double>("RusRoEcalGamma");
  gRusRoHcal = p.getParameter<double>("RusRoHcalGamma");
  gRusRoMuonIron = p.getParameter<double>("RusRoMuonIronGamma");
  gRusRoPreShower = p.getParameter<double>("RusRoPreShowerGamma");
  gRusRoCastor = p.getParameter<double>("RusRoCastorGamma");
  gRusRoWorld = p.getParameter<double>("RusRoWorldGamma");

  nRusRoEcal = p.getParameter<double>("RusRoEcalNeutron");
  nRusRoHcal = p.getParameter<double>("RusRoHcalNeutron");
  nRusRoMuonIron = p.getParameter<double>("RusRoMuonIronNeutron");
  nRusRoPreShower = p.getParameter<double>("RusRoPreShowerNeutron");
  nRusRoCastor = p.getParameter<double>("RusRoCastorNeutron");
  nRusRoWorld = p.getParameter<double>("RusRoWorldNeutron");

  pRusRoEcal = p.getParameter<double>("RusRoEcalProton");
  pRusRoHcal = p.getParameter<double>("RusRoHcalProton");
  pRusRoMuonIron = p.getParameter<double>("RusRoMuonIronProton");
  pRusRoPreShower = p.getParameter<double>("RusRoPreShowerProton");
  pRusRoCastor = p.getParameter<double>("RusRoCastorProton");
  pRusRoWorld = p.getParameter<double>("RusRoWorldProton");

  gRRactive = false;
  nRRactive = false;
  pRRactive = false;
  if(gRusRoEcal < 1.0 || gRusRoHcal < 1.0 || 
     gRusRoMuonIron < 1.0 || gRusRoPreShower < 1.0 || gRusRoCastor < 1.0 ||
     gRusRoWorld < 1.0) { gRRactive = true; }
  if(nRusRoEcal < 1.0 || nRusRoHcal < 1.0 || 
     nRusRoMuonIron < 1.0 || nRusRoPreShower < 1.0 || nRusRoCastor < 1.0 ||
     nRusRoWorld < 1.0) { nRRactive = true; }
  if(pRusRoEcal < 1.0 || pRusRoHcal < 1.0 || 
     pRusRoMuonIron < 1.0 || pRusRoPreShower < 1.0 || pRusRoCastor < 1.0 ||
     pRusRoWorld < 1.0) { pRRactive = true; }

  if ( p.exists("TestKillingOptions") ) {

    killInCalo = (p.getParameter<edm::ParameterSet>("TestKillingOptions")).getParameter<bool>("KillInCalo");
    killInCaloEfH = (p.getParameter<edm::ParameterSet>("TestKillingOptions")).getParameter<bool>("KillInCaloEfH");
    edm::LogWarning("SimG4CoreApplication") 
      << " *** Activating special test killing options in StackingAction \n"
                                            << " *** Kill secondaries in Calorimetetrs volume = " << killInCalo << "\n"
                                            << " *** Kill electromagnetic secondaries from hadrons in Calorimeters volume = " << killInCaloEfH;

  }

  edm::LogInfo("SimG4CoreApplication") << "StackingAction initiated with"
                       << " flag for saving decay products in "
                       << " Tracker: " << savePDandCinTracker
                                       << " in Calo: " << savePDandCinCalo
                                       << " in Muon: " << savePDandCinMuon
                       << "\n               saveFirstSecondary"
                       << ": " << saveFirstSecondary
                       << " Flag for tracking neutrino: "
                       << trackNeutrino << " Killing Flag "
                       << killHeavy << " protons below " 
                       << kmaxProton <<" MeV, neutrons below "
                       << kmaxNeutron << " MeV and ions"
                       << " below " << kmaxIon << " MeV\n"
                       << "               kill tracks with "
                       << "time larger than " << maxTrackTime
                       << " ns and kill Delta Ray flag set to "
                       << killDeltaRay;
  for (unsigned int i=0; i<maxTrackTimes.size(); i++) {
    maxTrackTimes[i] *= ns;
    edm::LogInfo("SimG4CoreApplication") << "StackingAction::MaxTrackTime for "
                     << maxTimeNames[i] << " is " 
                     << maxTrackTimes[i];
  }
  if(gRRactive) {
    edm::LogInfo("SimG4CoreApplication") 
      << "StackingAction: "
      << "Russian Roulette for gamma Elimit(MeV)= " 
      << nRusRoEnerLim/MeV << "\n"
      << "                 ECAL Prob= " << gRusRoEcal << "\n"
      << "                 HCAL Prob= " << gRusRoHcal << "\n"
      << "             MuonIron Prob= " << gRusRoMuonIron << "\n"
      << "            PreShower Prob= " << gRusRoPreShower << "\n"
      << "               CASTOR Prob= " << gRusRoCastor << "\n"
      << "                World Prob= " << gRusRoWorld;
  }
  if(nRRactive) {
    edm::LogInfo("SimG4CoreApplication") 
      << "StackingAction: "
      << "Russian Roulette for neutron Elimit(MeV)= " 
      << nRusRoEnerLim/MeV << "\n"
      << "                 ECAL Prob= " << nRusRoEcal << "\n"
      << "                 HCAL Prob= " << nRusRoHcal << "\n"
      << "             MuonIron Prob= " << nRusRoMuonIron << "\n"
      << "            PreShower Prob= " << nRusRoPreShower << "\n"
      << "               CASTOR Prob= " << nRusRoCastor << "\n"
      << "                World Prob= " << nRusRoWorld;
  }
  if(pRRactive) {
    edm::LogInfo("SimG4CoreApplication") 
      << "StackingAction: "
      << "Russian Roulette for proton Elimit(MeV)= " 
      << pRusRoEnerLim/MeV << "\n"
      << "                 ECAL Prob= " << pRusRoEcal << "\n"
      << "                 HCAL Prob= " << pRusRoHcal << "\n"
      << "             MuonIron Prob= " << pRusRoMuonIron << "\n"
      << "            PreShower Prob= " << pRusRoPreShower << "\n"
      << "               CASTOR Prob= " << pRusRoCastor << "\n"
      << "                World Prob= " << pRusRoWorld;
  }
  initPointer();
  newTA = new NewTrackAction();
}

StackingAction::~StackingAction ( )

virtual

Definition at line 168 of file StackingAction.cc.

References newTA.

{
  delete newTA;
}

Member Function Documentation

G4ClassificationOfNewTrack StackingAction::ClassifyNewTrack ( const G4Track *  aTrack )

virtual

Definition at line 173 of file StackingAction.cc.

References funct::abs(), calo, gRRactive, gRusRoCastor, gRusRoEcal, gRusRoEnerLim, gRusRoHcal, gRusRoMuonIron, gRusRoPreShower, gRusRoWorld, isItLongLived(), isItPrimaryDecayProductOrConversion(), isThisVolume(), ke, killDeltaRay, killHeavy, killInCalo, killInCaloEfH, kmaxIon, kmaxNeutron, kmaxProton, LogDebug, muon, newTA, nRRactive, nRusRoCastor, nRusRoEcal, nRusRoEnerLim, nRusRoHcal, nRusRoMuonIron, nRusRoPreShower, nRusRoWorld, NewTrackAction::primary(), mix_2012_Summer_inTimeOnly_cff::prob, pRRactive, pRusRoCastor, pRusRoEcal, pRusRoEnerLim, pRusRoHcal, pRusRoMuonIron, pRusRoPreShower, pRusRoWorld, regionCastor, regionEcal, regionHcal, regionMuonIron, regionPreShower, regionWorld, rrApplicable(), savePDandCinCalo, savePDandCinMuon, savePDandCinTracker, NewTrackAction::secondary(), CurrentG4Track::track(), tracker, and trackNeutrino.

{
  // G4 interface part
  G4ClassificationOfNewTrack classification = fUrgent;
  if (aTrack->GetCreatorProcess()==0 || aTrack->GetParentID()==0) {
    /*
    std::cout << "StackingAction: primary weight= " 
          << aTrack->GetWeight() << " "
          << aTrack->GetDefinition()->GetParticleName()
          << " " << aTrack->GetKineticEnergy()
          << " Id=" << aTrack->GetTrackID()
          << "  trackInfo " << aTrack->GetUserInformation() 
          << std::endl;
    */
    newTA->primary(aTrack);
    /*
    if (!trackNeutrino) {
      int pdg = std::abs(aTrack->GetDefinition()->GetPDGEncoding());
      if (pdg == 12 || pdg == 14 || pdg == 16 || pdg == 18) {
    classification = fKill;
      }
    }
    */
  } else if (aTrack->GetTouchable() == 0) {
    edm::LogError("SimG4CoreApplication")
      << "StackingAction: no touchable for track " << aTrack->GetTrackID()
      << " from " << aTrack->GetParentID()
      << " with PDG code " << aTrack->GetDefinition()->GetParticleName();
    classification = fKill;
  } else {
    int pdg = aTrack->GetDefinition()->GetPDGEncoding();

    if (aTrack->GetTrackStatus() == fStopAndKill) { classification = fKill; }
    if (killHeavy && classification != fKill) {
      double ke  = aTrack->GetKineticEnergy()/MeV;
      if (((pdg/1000000000 == 1) && (((pdg/10000)%100) > 0) && 
       (((pdg/10)%100) > 0) && (ke<kmaxIon)) || 
      ((pdg == 2212) && (ke < kmaxProton)) ||
      ((pdg == 2112) && (ke < kmaxNeutron))) { classification = fKill; }
    }
    if (!trackNeutrino  && classification != fKill) {
      if (pdg == 12 || pdg == 14 || pdg == 16 || pdg == 18) 
    classification = fKill;
    }
    if (classification != fKill && isItLongLived(aTrack)) 
      { classification = fKill; }
    if (killDeltaRay && classification != fKill) {
      if (aTrack->GetCreatorProcess()->GetProcessType() == fElectromagnetic &&
      aTrack->GetCreatorProcess()->GetProcessSubType() == fIonisation)
    classification = fKill;
    }
    if (killInCalo && classification != fKill) {
      if ( isThisVolume(aTrack->GetTouchable(),calo)) { 
        classification = fKill; 
      }
    }
    // Russian roulette && MC truth
    if(classification != fKill) {
      const G4Track * mother = CurrentG4Track::track();
      int flag = 0;
      if (savePDandCinTracker && isThisVolume(aTrack->GetTouchable(),tracker)) {
    flag = isItPrimaryDecayProductOrConversion(aTrack, *mother);
      }
      if (savePDandCinCalo && 0 == flag 
      && isThisVolume(aTrack->GetTouchable(),calo)) {
    flag = isItPrimaryDecayProductOrConversion(aTrack, *mother);
      }
      if (savePDandCinMuon && 0 == flag 
      && isThisVolume(aTrack->GetTouchable(),muon)) {
    flag = isItPrimaryDecayProductOrConversion(aTrack, *mother);
      }

      // Russian roulette
      if(2112 == pdg || 22 == pdg || 2212 == pdg) {
    double currentWeight = aTrack->GetWeight();

    if(1.0 >= currentWeight) {
      double prob = 1.0;
      double elim = 0.0;

      // neutron
      if(nRRactive && pdg == 2112) {
        elim = nRusRoEnerLim;
        G4Region* reg = aTrack->GetVolume()->GetLogicalVolume()->GetRegion();
        if(reg == regionEcal)             { prob = nRusRoEcal; }
        else if(reg == regionHcal)        { prob = nRusRoHcal; }
        else if(reg == regionMuonIron)    { prob = nRusRoMuonIron; }
        else if(reg == regionPreShower)   { prob = nRusRoPreShower; }
        else if(reg == regionCastor)      { prob = nRusRoCastor; }
        else if(reg == regionWorld)       { prob = nRusRoWorld; }

        // gamma
      } else if(gRRactive && pdg == 22) {
        elim = gRusRoEnerLim;

        G4Region* reg = aTrack->GetVolume()->GetLogicalVolume()->GetRegion();
        if(reg == regionEcal || reg == regionPreShower) {
              if(rrApplicable(aTrack, *mother)) {
        if(reg == regionEcal)         { prob = gRusRoEcal; }
                else                          { prob = gRusRoPreShower; }
          }
        } else {
          if(reg == regionHcal)           { prob = gRusRoHcal; }
          else if(reg == regionMuonIron)  { prob = gRusRoMuonIron; }
          else if(reg == regionCastor)    { prob = gRusRoCastor; }
          else if(reg == regionWorld)     { prob = gRusRoWorld; }
        }

        // proton
      } else if(pRRactive && pdg == 2212) {
        elim = pRusRoEnerLim;
        G4Region* reg = aTrack->GetVolume()->GetLogicalVolume()->GetRegion();
        if(reg == regionEcal)             { prob = pRusRoEcal; }
        else if(reg == regionHcal)        { prob = pRusRoHcal; }
        else if(reg == regionMuonIron)    { prob = pRusRoMuonIron; }
        else if(reg == regionPreShower)   { prob = pRusRoPreShower; }
        else if(reg == regionCastor)      { prob = pRusRoCastor; }
        else if(reg == regionWorld)       { prob = pRusRoWorld; }
      }
      if(prob < 1.0 && aTrack->GetKineticEnergy() < elim) {
        if(G4UniformRand() < prob) {
          const_cast<G4Track*>(aTrack)->SetWeight(currentWeight/prob);
        } else {
          classification = fKill;
        }
      }
    }
      }
    
      if(classification != fKill && killInCaloEfH) {
    int pdgMother = mother->GetDefinition()->GetPDGEncoding();
    if ( (pdg == 22 || std::abs(pdg) == 11) && 
         (std::abs(pdgMother) < 11 || std::abs(pdgMother) > 17) && 
         pdgMother != 22  ) {
      if ( isThisVolume(aTrack->GetTouchable(),calo)) { 
        classification = fKill; 
      }
    }
      }
      
      if (classification != fKill) {
    newTA->secondary(aTrack, *mother, flag);
      }
    }
  /*
  double wt2 = aTrack->GetWeight();
  if(wt2 != 1.0) { 
    G4Region* reg = aTrack->GetVolume()->GetLogicalVolume()->GetRegion();
    std::cout << "StackingAction: weight= " << wt2 << " "
          << aTrack->GetDefinition()->GetParticleName()
          << " " << aTrack->GetKineticEnergy()
          << " Id=" << aTrack->GetTrackID()
          << " IdP=" << aTrack->GetParentID();
    const G4VProcess* pr = aTrack->GetCreatorProcess();
    if(pr) std::cout << " from  " << pr->GetProcessName();
    if(reg) std::cout << " in  " << reg->GetName()
              << "  trackInfo " << aTrack->GetUserInformation(); 
    std::cout << std::endl;
  }
  */
  }
#ifdef DebugLog
  LogDebug("SimG4CoreApplication") << "StackingAction:Classify Track "
                   << aTrack->GetTrackID() << " Parent " 
                   << aTrack->GetParentID() << " Type "
                   << aTrack->GetDefinition()->GetParticleName() 
                   << " K.E. " << aTrack->GetKineticEnergy()/MeV
                   << " MeV from process/subprocess " 
                   << aTrack->GetCreatorProcess()->GetProcessType() 
                   << "|"
                   <<aTrack->GetCreatorProcess()->GetProcessSubType()
                   << " as " << classification << " Flag " << flag;
#endif
  
  return classification;
}

void StackingAction::initPointer ( )

private

Definition at line 354 of file StackingAction.cc.

References calo, gRusRoCastor, gRusRoEcal, gRusRoHcal, gRusRoMuonIron, gRusRoPreShower, i, killInCalo, killInCaloEfH, maxTimeNames, maxTimeRegions, maxTrackTimes, muon, mergeVDriftHistosByStation::name, nRusRoCastor, nRusRoEcal, nRusRoHcal, nRusRoMuonIron, nRusRoPreShower, nRusRoWorld, pileupDistInMC::num, pRusRoCastor, pRusRoEcal, pRusRoHcal, pRusRoMuonIron, pRusRoPreShower, pRusRoWorld, regionCastor, regionEcal, regionHcal, regionMuonIron, regionPreShower, regionWorld, savePDandCinCalo, savePDandCinMuon, savePDandCinTracker, and tracker.

Referenced by StackingAction().

                                 {

  const G4LogicalVolumeStore * lvs = G4LogicalVolumeStore::GetInstance();
  if (lvs) {
    std::vector<G4LogicalVolume*>::const_iterator lvcite;
    for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) {
      if (savePDandCinTracker) {
        if (strcmp("Tracker",(*lvcite)->GetName().c_str()) == 0) {
      tracker.push_back(*lvcite);
    }
        if (strcmp("BEAM",(*lvcite)->GetName().substr(0,4).c_str()) == 0) {
      tracker.push_back(*lvcite);
    }
      }
      if (savePDandCinCalo || killInCalo || killInCaloEfH ) {
        if (strcmp("CALO",(*lvcite)->GetName().c_str()) == 0) {
      calo.push_back(*lvcite);
    }
        if (strcmp("VCAL",(*lvcite)->GetName().c_str()) == 0) {
      calo.push_back(*lvcite);
    }
      }
      if (savePDandCinMuon) {
        if (strcmp("MUON",(*lvcite)->GetName().c_str()) == 0) {
      muon.push_back(*lvcite);
    }
      }
    }
    edm::LogInfo("SimG4CoreApplication") << "# of LV for Tracker " 
                     << tracker.size() << " for Calo "
                                         << calo.size() << " for Muon "
                     << muon.size();
    for (unsigned int i=0; i<tracker.size(); ++i)
      edm::LogInfo("SimG4CoreApplication") << "Tracker vol " << i << " name "
                       << tracker[i]->GetName();
    for (unsigned int i=0; i<calo.size(); ++i)
      edm::LogInfo("SimG4CoreApplication") << "Calorimeter vol " <<i <<" name "
                       << calo[i]->GetName();
    for (unsigned int i=0; i<muon.size(); ++i)
      edm::LogInfo("SimG4CoreApplication") << "Muon vol " << i << " name "
                       << muon[i]->GetName();
  }

  std::vector<double> tofs;

  // Russian roulette
  const G4RegionStore * rs = G4RegionStore::GetInstance();
  std::vector<G4Region*>::const_iterator rcite;
  for (rcite = rs->begin(); rcite != rs->end(); rcite++) {
    if ((gRusRoEcal < 1.0 || nRusRoEcal < 1.0 || pRusRoEcal < 1.0) && 
    (*rcite)->GetName() == "EcalRegion") {  
      regionEcal = (*rcite); 
    }
    if ((gRusRoHcal < 1.0 || nRusRoHcal < 1.0 || pRusRoHcal < 1.0) && 
    (*rcite)->GetName() == "HcalRegion") {  
      regionHcal = (*rcite); 
    }
    if ((gRusRoMuonIron < 1.0 || nRusRoMuonIron < 1.0 || pRusRoMuonIron < 1.0) && 
    (*rcite)->GetName() == "MuonIron") {  
      regionMuonIron = (*rcite); 
    }
    if ((gRusRoPreShower < 1.0 || nRusRoPreShower < 1.0 || pRusRoPreShower < 1.0) 
    && (*rcite)->GetName() == "PreshowerRegion") {  
      regionPreShower = (*rcite); 
    }
    if ((gRusRoCastor < 1.0 || nRusRoCastor < 1.0 || pRusRoCastor < 1.0) && 
    (*rcite)->GetName() == "CastorRegion") {  
      regionCastor = (*rcite); 
    }
    if ((nRusRoWorld < 1.0 || nRusRoWorld < 1.0 || pRusRoWorld < 1.0) && 
    (*rcite)->GetName() == "DefaultRegionForTheWorld") {  
      regionWorld = (*rcite); 
    }

    // time limits
    unsigned int num = maxTimeNames.size();
    if (num > 0) {
      for (unsigned int i=0; i<num; i++) {
    if ((*rcite)->GetName() == (G4String)(maxTimeNames[i])) {
      maxTimeRegions.push_back(*rcite);
      tofs.push_back(maxTrackTimes[i]);
      break;
    }
      }
    }
  }
  if(0 < tofs.size()) {
    for (unsigned int i=0; i<tofs.size(); i++) {
      maxTrackTimes[i] = tofs[i];
      G4String name = "Unknown";
      if (maxTimeRegions[i]) { name = maxTimeRegions[i]->GetName(); }
      edm::LogInfo("SimG4CoreApplication") << name << " with pointer " 
                       << maxTimeRegions[i]<<" KE cut off "
                       << maxTrackTimes[i];
    }
  }
}

bool StackingAction::isItLongLived ( const G4Track *  aTrack ) const

private

Definition at line 514 of file StackingAction.cc.

References i, maxTimeRegions, maxTrackTime, maxTrackTimes, and cond::rpcobgas::time.

Referenced by ClassifyNewTrack().

{
  bool   flag = false;
  double time = (aTrack->GetGlobalTime())/nanosecond;
  double tofM = maxTrackTime;
  if (maxTimeRegions.size() > 0) {
    G4Region* reg = 
      aTrack->GetTouchable()->GetVolume(0)->GetLogicalVolume()->GetRegion();
    for (unsigned int i=0; i<maxTimeRegions.size(); i++) {
      if (reg == maxTimeRegions[i]) {
    tofM = maxTrackTimes[i];
    break;
      }
    }
  }
  if (time > tofM) { flag = true; }
  return flag;
}

int StackingAction::isItPrimaryDecayProductOrConversion ( const G4Track *  aTrack, const G4Track &  mother  ) const

private

Definition at line 467 of file StackingAction.cc.

References TrackInformation::isPrimary().

Referenced by ClassifyNewTrack().

{
  int flag = 0;
  TrackInformationExtractor extractor;
  const TrackInformation & motherInfo(extractor(mother));
  // Check whether mother is a primary
  if (motherInfo.isPrimary()) {
    if (aTrack->GetCreatorProcess()->GetProcessType() == fDecay) {
      flag = 1;
    } else if (aTrack->GetCreatorProcess()->GetProcessSubType()==fGammaConversion) {
      flag = 2;
    } else {
      flag = 3;
    } 
  }   
  return flag;
}

bool StackingAction::isThisVolume ( const G4VTouchable *  touch, std::vector< G4LogicalVolume * > &  lvs  ) const

private

Definition at line 452 of file StackingAction.cc.

References prof2calltree::count, cuy::ii, and testEve_cfg::level.

Referenced by ClassifyNewTrack().

{
  bool flag = false;
  if (lvs.size() > 0 && touch !=0) {
    int level = ((touch->GetHistoryDepth())+1);
    if (level >= 3) {
      unsigned int  ii = (unsigned int)(level - 3);
      flag = (std::count(lvs.begin(),lvs.end(),
             (touch->GetVolume(ii)->GetLogicalVolume())) != 0);
    }
  }
  return flag;
}

void StackingAction::NewStage

(

)

Definition at line 350 of file StackingAction.cc.

{}

void StackingAction::PrepareNewEvent

(

)

Definition at line 352 of file StackingAction.cc.

{}

bool StackingAction::rrApplicable ( const G4Track *  aTrack, const G4Track &  mother  ) const

private

Definition at line 486 of file StackingAction.cc.

References TrackInformation::genParticlePID().

Referenced by ClassifyNewTrack().

{
  bool flag = true;
  TrackInformationExtractor extractor;
  const TrackInformation & motherInfo(extractor(mother));
  // Check whether mother is a primary
  //if (motherInfo.isPrimary()) {
  // }   
  int genID = motherInfo.genParticlePID();
  if(22 == genID || 11 == genID || -11 == genID) { flag = false; }
    
  /*
  //check if the primary was g, e+, e-
  int genID = motherInfo.genParticlePID();
  double genp = motherInfo.genParticleP();
  std::cout << "Track# " << aTrack->GetTrackID() << "  " 
        << aTrack->GetDefinition()->GetParticleName()  
        << "  E(MeV)= " << aTrack->GetKineticEnergy()/MeV 
        << " mother: " << mother.GetTrackID()
        << "  " << mother.GetDefinition()->GetParticleName()
        << " E(GeV)= " <<  mother.GetKineticEnergy()/GeV
        << " flag: " << flag << " genID= " << genID 
        << " p(GeV)= " << genp/GeV << std::endl; 
    */
  return flag;
}

Member Data Documentation

std::vector<G4LogicalVolume*> StackingAction::calo

private

Definition at line 46 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and initPointer().

EventAction* StackingAction::eventAction_

private

Definition at line 36 of file StackingAction.h.

bool StackingAction::gRRactive

private

Definition at line 83 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

double StackingAction::gRusRoCastor

private

Definition at line 76 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::gRusRoEcal

private

Definition at line 64 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::gRusRoEnerLim

private

Definition at line 59 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

double StackingAction::gRusRoHcal

private

Definition at line 67 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::gRusRoMuonIron

private

Definition at line 70 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::gRusRoPreShower

private

Definition at line 73 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::gRusRoWorld

private

Definition at line 79 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

bool StackingAction::killDeltaRay

private

Definition at line 40 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

bool StackingAction::killHeavy

private

Definition at line 40 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

bool StackingAction::killInCalo

private

Definition at line 39 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

bool StackingAction::killInCaloEfH

private

Definition at line 39 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::kmaxIon

private

Definition at line 41 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

double StackingAction::kmaxNeutron

private

Definition at line 41 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

double StackingAction::kmaxProton

private

Definition at line 41 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

std::vector<std::string> StackingAction::maxTimeNames

private

Definition at line 44 of file StackingAction.h.

Referenced by initPointer(), and StackingAction().

std::vector<G4Region*> StackingAction::maxTimeRegions

private

Definition at line 45 of file StackingAction.h.

Referenced by initPointer(), and isItLongLived().

double StackingAction::maxTrackTime

private

Definition at line 42 of file StackingAction.h.

Referenced by isItLongLived(), and StackingAction().

std::vector<double> StackingAction::maxTrackTimes

private

Definition at line 43 of file StackingAction.h.

Referenced by initPointer(), isItLongLived(), and StackingAction().

std::vector<G4LogicalVolume*> StackingAction::muon

private

Definition at line 46 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and initPointer().

NewTrackAction* StackingAction::newTA

private

Definition at line 48 of file StackingAction.h.

Referenced by ClassifyNewTrack(), StackingAction(), and ~StackingAction().

bool StackingAction::nRRactive

private

Definition at line 84 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

double StackingAction::nRusRoCastor

private

Definition at line 77 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::nRusRoEcal

private

Definition at line 65 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::nRusRoEnerLim

private

Definition at line 60 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

double StackingAction::nRusRoHcal

private

Definition at line 68 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::nRusRoMuonIron

private

Definition at line 71 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::nRusRoPreShower

private

Definition at line 74 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::nRusRoWorld

private

Definition at line 80 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

bool StackingAction::pRRactive

private

Definition at line 85 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

double StackingAction::pRusRoCastor

private

Definition at line 78 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::pRusRoEcal

private

Definition at line 66 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::pRusRoEnerLim

private

Definition at line 61 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().

double StackingAction::pRusRoHcal

private

Definition at line 69 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::pRusRoMuonIron

private

Definition at line 72 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::pRusRoPreShower

private

Definition at line 75 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

double StackingAction::pRusRoWorld

private

Definition at line 81 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

G4Region* StackingAction::regionCastor

private

Definition at line 55 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

G4Region* StackingAction::regionEcal

private

Definition at line 51 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

G4Region* StackingAction::regionHcal

private

Definition at line 52 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

G4Region* StackingAction::regionMuonIron

private

Definition at line 53 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

G4Region* StackingAction::regionPreShower

private

Definition at line 54 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

G4Region* StackingAction::regionWorld

private

Definition at line 56 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

bool StackingAction::saveFirstSecondary

private

Definition at line 38 of file StackingAction.h.

Referenced by StackingAction().

bool StackingAction::savePDandCinCalo

private

Definition at line 37 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

bool StackingAction::savePDandCinMuon

private

Definition at line 38 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

bool StackingAction::savePDandCinTracker

private

Definition at line 37 of file StackingAction.h.

Referenced by ClassifyNewTrack(), initPointer(), and StackingAction().

std::vector<G4LogicalVolume*> StackingAction::tracker

private

Definition at line 46 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and initPointer().

bool StackingAction::trackNeutrino

private

Definition at line 40 of file StackingAction.h.

Referenced by ClassifyNewTrack(), and StackingAction().