G4ProcessHelper Class Reference

#include <G4ProcessHelper.h>

Public Member Functions
G4bool	ApplicabilityTester (const G4ParticleDefinition &aPart)
	G4ProcessHelper (const edm::ParameterSet &p, CustomParticleFactory *ptr)
ReactionProduct	GetFinalState (const G4Track &aTrack, G4ParticleDefinition *&aTarget)
G4double	GetInclusiveCrossSection (const G4DynamicParticle *aParticle, const G4Element *anElement)
	~G4ProcessHelper ()

Private Member Functions
	G4ProcessHelper (const G4ProcessHelper &)=delete
G4ProcessHelper &	operator= (const G4ProcessHelper &)=delete
G4double	PhaseSpace (const ReactionProduct &aReaction, const G4DynamicParticle *aDynamicParticle)
G4double	Pom (const double boost)
G4bool	ReactionGivesBaryon (const ReactionProduct &aReaction)
G4bool	ReactionIsPossible (const ReactionProduct &aReaction, const G4DynamicParticle *aDynamicParticle)
G4double	ReactionProductMass (const ReactionProduct &aReaction, const G4DynamicParticle *aDynamicParticle)
void	ReadAndParse (const G4String &str, std::vector< G4String > &tokens, const G4String &delimiters=" ")
G4double	Regge (const double boost)

Private Attributes
double	amplitude
G4double	checkfraction
double	ek_0
CustomParticleFactory *	fParticleFactory
double	gamma
TProfile *	h_q_gamma
TProfile *	h_q_p
TH1D *	h_sqrts
TProfile *	h_xsec_cms
TProfile *	h_xsec_lab
double	hadronlifetime
std::map< const G4ParticleDefinition *, G4bool >	known_particles
double	mixing
G4int	n_22
G4int	n_23
ReactionMap	nReactionMap
std::map< G4String, G4double >	parameters
G4ParticleTable *	particleTable
ReactionMap	pReactionMap
bool	reggemodel
bool	resonant
double	suppressionfactor
HistoHelper *	theHistoHelper
G4ParticleDefinition *	theNeutron
G4ParticleDefinition *	theProton
G4ParticleDefinition *	theRbaryoncloud
ReactionMap *	theReactionMap
G4ParticleDefinition *	theRmesoncloud
G4ParticleDefinition *	theTarget

Detailed Description

Definition at line 26 of file G4ProcessHelper.h.

Constructor & Destructor Documentation

G4ProcessHelper::G4ProcessHelper	(	const edm::ParameterSet &	p,
		CustomParticleFactory *	ptr
	)

Definition at line 19 of file G4ProcessHelper.cc.

References CustomPhysics_cfi::amplitude, edm::FileInPath::fullPath(), CustomPhysics_cfi::gamma, edm::ParameterSet::getParameter(), GeV, mps_fire::i, geometryCSVtoXML::line, CustomPhysics_cfi::mixing, parseEventContent::prod, CustomPhysics_cfi::resonant, alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, and edmPickEvents::target.

                                                             {
  fParticleFactory = ptr;

  particleTable = G4ParticleTable::GetParticleTable();

  theProton = particleTable->FindParticle("proton");
  theNeutron = particleTable->FindParticle("neutron");
  
  G4String line;

  edm::FileInPath fp = p.getParameter<edm::FileInPath>("processesDef");
  std::string processDefFilePath = fp.fullPath();
  std::ifstream process_stream (processDefFilePath.c_str());

  resonant = p.getParameter<bool>("resonant");
  ek_0 = p.getParameter<double>("resonanceEnergy")*GeV;
  gamma = p.getParameter<double>("gamma")*GeV;
  amplitude = p.getParameter<double>("amplitude")*millibarn;
  suppressionfactor = p.getParameter<double>("reggeSuppression");
  hadronlifetime = p.getParameter<double>("hadronLifeTime");
  reggemodel = p.getParameter<bool>("reggeModel");
  mixing = p.getParameter<double>("mixing");
  
  edm::LogInfo("SimG4CoreCustomPhysics")
    <<"ProcessHelper: Read in physics parameters:"
    <<"\n Resonant = "<< resonant 
    <<"\n ResonanceEnergy = "<<ek_0/GeV<<" GeV"
    <<"\n Gamma = "<<gamma/GeV<<" GeV"
    <<"\n Amplitude = "<<amplitude/millibarn<<" millibarn"
    <<"ReggeSuppression = "<<100*suppressionfactor<<" %"
    <<"HadronLifeTime = "<<hadronlifetime<<" s"
    <<"ReggeModel = "<< reggemodel 
    <<"Mixing = "<< mixing*100 <<" %";

  checkfraction = 0;
  n_22 = 0;
  n_23 = 0;

  while(getline(process_stream,line)){
    std::vector<G4String> tokens;
    //Getting a line
    ReadAndParse(line,tokens,"#");
    //Important info
    G4String incident = tokens[0];
    
    G4ParticleDefinition* incidentDef = particleTable->FindParticle(incident);
    //particleTable->DumpTable();
    G4int incidentPDG = incidentDef->GetPDGEncoding();
    known_particles[incidentDef]=true;

    G4String target = tokens[1];
    edm::LogInfo("SimG4CoreCustomPhysics")
      <<"ProcessHelper: Incident "<<incident
      <<"; Target "<<target;
    
    // Making a ReactionProduct
    ReactionProduct prod;
    for (size_t i = 2; i != tokens.size();i++){
      G4String part = tokens[i];
      if (particleTable->contains(part))
    {
      prod.push_back(particleTable->FindParticle(part)->GetPDGEncoding());
    } else {
      G4Exception("G4ProcessHelper", "UnkownParticle", FatalException,
              "Initialization: The reaction product list contained an unknown particle");
    }
    }
    if (target == "proton")
      {
    pReactionMap[incidentPDG].push_back(prod);
      } else if (target == "neutron") {
    nReactionMap[incidentPDG].push_back(prod);
      } else {
    G4Exception("G4ProcessHelper", "IllegalTarget", FatalException,
            "Initialization: The reaction product list contained an illegal target particle");
      }
   }

  process_stream.close();

  for(auto part : fParticleFactory->GetCustomParticles()) {
    CustomParticle* particle = dynamic_cast<CustomParticle*>(part);
    if(particle) {
      edm::LogInfo("SimG4CoreCustomPhysics")
    <<"ProcessHelper: Lifetime of "<<part->GetParticleName()<<" set to "
    <<particle->GetPDGLifeTime()/s<<" s;"
    <<" isStable: "<<particle->GetPDGStable();
    }
  }
}

G4ProcessHelper::~G4ProcessHelper

(

)

Definition at line 111 of file G4ProcessHelper.cc.

                                 {
}

G4ProcessHelper::G4ProcessHelper ( const G4ProcessHelper &  )

privatedelete

Member Function Documentation

G4bool G4ProcessHelper::ApplicabilityTester

(

const G4ParticleDefinition &

aPart

)

Definition at line 114 of file G4ProcessHelper.cc.

Referenced by FullModelHadronicProcess::IsApplicable().

                                                                            {
  const G4ParticleDefinition* aP = &aPart; 
  if (known_particles[aP]) return true;
  return false;
}

ReactionProduct G4ProcessHelper::GetFinalState	(	const G4Track &	aTrack,
		G4ParticleDefinition *&	aTarget
	)

Definition at line 193 of file G4ProcessHelper.cc.

References ecalTB2006H4_GenSimDigiReco_cfg::G4cerr, mps_fire::i, createfilelist::int, CustomPhysics_cfi::mixing, CustomPDGParser::s_isMesonino(), CustomPDGParser::s_isRBaryon(), CustomPDGParser::s_isRMeson(), CustomPDGParser::s_isSbaryon(), TauGenJetsDecayModeSelectorAllHadrons_cfi::select, and findQualityFiles::size.

Referenced by FullModelHadronicProcess::PostStepDoIt().

                                                                                                   {

  const G4DynamicParticle* aDynamicParticle = aTrack.GetDynamicParticle();

  //-----------------------------------------------
  // Choose n / p as target
  // and get ReactionProductList pointer
  //-----------------------------------------------

  G4Material* aMaterial = aTrack.GetMaterial();
  const G4ElementVector* theElementVector = aMaterial->GetElementVector() ;
  const G4double* NbOfAtomsPerVolume = aMaterial->GetVecNbOfAtomsPerVolume();

  G4double NumberOfProtons=0;
  G4double NumberOfNucleons=0;

  for ( size_t elm=0 ; elm < aMaterial->GetNumberOfElements() ; elm++ )
    {
      //Summing number of protons per unit volume
      NumberOfProtons += NbOfAtomsPerVolume[elm]*(*theElementVector)[elm]->GetZ();
      //Summing nucleons (not neutrons)
      NumberOfNucleons += NbOfAtomsPerVolume[elm]*(*theElementVector)[elm]->GetN();
    }

  if(G4UniformRand()<NumberOfProtons/NumberOfNucleons)
    {
      theReactionMap = &pReactionMap;
      theTarget = theProton;
    } else {
      theReactionMap = &nReactionMap;
      theTarget = theNeutron;
    }
  aTarget = theTarget;

  G4int theIncidentPDG = aDynamicParticle->GetDefinition()->GetPDGEncoding();

  if(reggemodel
     &&CustomPDGParser::s_isMesonino(theIncidentPDG)
     && G4UniformRand()*mixing>0.5
     &&aDynamicParticle->GetDefinition()->GetPDGCharge()==0.
     ) 
    {
      //      G4cout<<"Oscillating..."<<G4endl;
      theIncidentPDG *= -1;
    }


  bool baryonise=false;
  
  if(reggemodel
     && G4UniformRand()>0.9
     &&(
    (CustomPDGParser::s_isMesonino(theIncidentPDG)&&theIncidentPDG>0)
    ||
    CustomPDGParser::s_isRMeson(theIncidentPDG)
    )
     )  
    baryonise=true;

  //Making a pointer directly to the ReactionProductList we are looking at. Makes life easier :-)
  ReactionProductList*  aReactionProductList = &((*theReactionMap)[theIncidentPDG]);

  //-----------------------------------------------
  // Count processes
  // kinematic check
  // compute number of 2 -> 2 and 2 -> 3 processes
  //-----------------------------------------------

  G4int N22 = 0; //Number of 2 -> 2 processes
  G4int N23 = 0; //Number of 2 -> 3 processes. Couldn't think of more informative names
  
  //This is the list to be populated
  ReactionProductList theReactionProductList;
  std::vector<bool> theChargeChangeList;

  for (ReactionProductList::iterator prod_it = aReactionProductList->begin();
       prod_it != aReactionProductList->end();
       prod_it++){
    G4int secondaries = prod_it->size();
    // If the reaction is not possible we will not consider it
    if(ReactionIsPossible(*prod_it,aDynamicParticle)
       && (!reggemodel ||
      (baryonise&&ReactionGivesBaryon(*prod_it)) 
      ||
       (!baryonise&&!ReactionGivesBaryon(*prod_it))
      ||
      (CustomPDGParser::s_isSbaryon(theIncidentPDG))
      ||
      (CustomPDGParser::s_isRBaryon(theIncidentPDG))
      )
       )
      {
      // The reaction is possible. Let's store and count it
      theReactionProductList.push_back(*prod_it);
      if (secondaries == 2){
    N22++;
      } else if (secondaries ==3) {
    N23++;
      } else {
    G4cerr << "ReactionProduct has unsupported number of secondaries: "<<secondaries<<G4endl;
      }
    } /*else {
      edm::LogInfo("SimG4CoreCustomPhysics")<<"There was an impossible process"<<G4endl;
      }*/
  }
  //  edm::LogInfo("SimG4CoreCustomPhysics")<<"The size of the ReactionProductList is: "<<theReactionProductList.size()<<G4endl;

  if (theReactionProductList.empty()) G4Exception("G4ProcessHelper", "NoProcessPossible", FatalException,
                            "GetFinalState: No process could be selected from the given list.");

  // For the Regge model no phase space considerations. We pick a process at random
  if(reggemodel)
    {
      int n_rps = theReactionProductList.size();
      int select = (int)( G4UniformRand()*n_rps);
      //      G4cout<<"Possible: "<<n_rps<<", chosen: "<<select<<G4endl;
      return theReactionProductList[select];
    }

  // Fill a probability map. Remember total probability
  // 2->2 is 0.15*1/n_22 2->3 uses phase space
  G4double p22 = 0.15;
  G4double p23 = 1-p22; // :-)

  std::vector<G4double> Probabilities;
  std::vector<G4bool> TwotoThreeFlag;
  
  G4double CumulatedProbability = 0;

  // To each ReactionProduct we assign a cumulated probability and a flag
  // discerning between 2 -> 2 and 2 -> 3
  for (unsigned int i = 0; i != theReactionProductList.size(); i++){
    if (theReactionProductList[i].size() == 2) {
      CumulatedProbability += p22/N22;
      TwotoThreeFlag.push_back(false);
    } else {
      CumulatedProbability += p23/N23;
      TwotoThreeFlag.push_back(true);
    }
    Probabilities.push_back(CumulatedProbability);
    //    edm::LogInfo("SimG4CoreCustomPhysics")<<"Pushing back cumulated probability: "<<CumulatedProbability<<G4endl;
  }

  //Renormalising probabilities
  //  edm::LogInfo("SimG4CoreCustomPhysics")<<"Probs: ";
  for (std::vector<G4double>::iterator it = Probabilities.begin();
       it != Probabilities.end();
       it++)
    {
      *it /= CumulatedProbability;
      //      edm::LogInfo("SimG4CoreCustomPhysics")<<*it<<" ";
    }
  //  edm::LogInfo("SimG4CoreCustomPhysics")<<G4endl;

  // Choosing ReactionProduct

  G4bool selected = false;
  G4int tries = 0;
  //  ReactionProductList::iterator prod_it;

  //Keep looping over the list until we have a choice, or until we have tried 100 times  
  unsigned int i;
  while(!selected && tries < 100){
    i=0;
    G4double dice = G4UniformRand();
    // edm::LogInfo("SimG4CoreCustomPhysics")<<"What's the dice?"<<dice<<G4endl;
    while(dice>Probabilities[i] && i<theReactionProductList.size()){
      //      edm::LogInfo("SimG4CoreCustomPhysics")<<"i: "<<i<<G4endl;
      i++;
    }

    //    edm::LogInfo("SimG4CoreCustomPhysics")<<"Chosen i: "<<i<<G4endl;

    if(!TwotoThreeFlag[i]) {
      // 2 -> 2 processes are chosen immediately
      selected = true;
    } else {
      // 2 -> 3 processes require a phase space lookup
      if (PhaseSpace(theReactionProductList[i],aDynamicParticle)>G4UniformRand()) selected = true;
      //selected = true;
    }
    //    double suppressionfactor=0.5;
    if(selected&&particleTable->FindParticle(theReactionProductList[i][0])->GetPDGCharge()!=aDynamicParticle->GetDefinition()->GetPDGCharge())
      {
    /*
    edm::LogInfo("SimG4CoreCustomPhysics")<<"Incoming particle "<<aDynamicParticle->GetDefinition()->GetParticleName()
          <<" has charge "<<aDynamicParticle->GetDefinition()->GetPDGCharge()<<G4endl;
    edm::LogInfo("SimG4CoreCustomPhysics")<<"Suggested particle "<<particleTable->FindParticle(theReactionProductList[i][0])->GetParticleName()
          <<" has charge "<<particleTable->FindParticle(theReactionProductList[i][0])->GetPDGCharge()<<G4endl;
    */
    if(G4UniformRand()<suppressionfactor) selected = false;
      }
    tries++;
    //    edm::LogInfo("SimG4CoreCustomPhysics")<<"Tries: "<<tries<<G4endl;
  }
  if(tries>=100) G4cerr<<"Could not select process!!!!"<<G4endl;

  //  edm::LogInfo("SimG4CoreCustomPhysics")<<"So far so good"<<G4endl;
  //  edm::LogInfo("SimG4CoreCustomPhysics")<<"Sec's: "<<theReactionProductList[i].size()<<G4endl;
  
  //Updating checkfraction:
  if (theReactionProductList[i].size()==2) {
    n_22++;
  } else {
    n_23++;
  }

  checkfraction = (1.0*n_22)/(n_22+n_23);
  //  edm::LogInfo("SimG4CoreCustomPhysics")<<"n_22: "<<n_22<<" n_23: "<<n_23<<" Checkfraction: "<<checkfraction<<G4endl;
  //  edm::LogInfo("SimG4CoreCustomPhysics") <<"Biig number: "<<n_22+n_23<<G4endl;
  //Return the chosen ReactionProduct
  return theReactionProductList[i];
}

G4double G4ProcessHelper::GetInclusiveCrossSection	(	const G4DynamicParticle *	aParticle,
		const G4Element *	anElement
	)

Definition at line 120 of file G4ProcessHelper.cc.

References CustomPhysics_cfi::amplitude, CustomPhysics_cfi::gamma, dataset::name, funct::pow(), dttmaxenums::R, CustomPhysics_cfi::resonant, CustomPDGParser::s_containedQuarks(), CustomPDGParser::s_isMesonino(), CustomPDGParser::s_isRBaryon(), CustomPDGParser::s_isRGlueball(), CustomPDGParser::s_isRMeson(), CustomPDGParser::s_isSbaryon(), and mathSSE::sqrt().

Referenced by FullModelHadronicProcess::GetMicroscopicCrossSection().

                                                      {

  //We really do need a dedicated class to handle the cross sections. They might not always be constant


  //Disassemble the PDG-code

  G4int thePDGCode = aParticle->GetDefinition()->GetPDGEncoding();
  double boost = (aParticle->GetKineticEnergy()+aParticle->GetMass())/aParticle->GetMass();
  //  G4cout<<"thePDGCode: "<<thePDGCode<<G4endl;
  G4double theXsec = 0;
  G4String name = aParticle->GetDefinition()->GetParticleName();
  if(!reggemodel)
  {
  //Flat cross section
  if(CustomPDGParser::s_isRGlueball(thePDGCode)) {
    theXsec = 24 * millibarn;
  } else {
    std::vector<G4int> nq=CustomPDGParser::s_containedQuarks(thePDGCode);
    //    edm::LogInfo("SimG4CoreCustomPhysics")<<"Number of quarks: "<<nq.size()<<G4endl;
    for (std::vector<G4int>::iterator it = nq.begin();
     it != nq.end();
     it++)
      {
    //    edm::LogInfo("SimG4CoreCustomPhysics")<<"Quarkvector: "<<*it<<G4endl;
    if (*it == 1 || *it == 2) theXsec += 12 * millibarn;
    if (*it == 3) theXsec += 6 * millibarn;
      }
   }
  } 
  else 
  {  //reggemodel
    double R = Regge(boost);
    double P = Pom(boost);
    if(thePDGCode>0)
      {
    if(CustomPDGParser::s_isMesonino(thePDGCode)) theXsec=(P+R)*millibarn;
    if(CustomPDGParser::s_isSbaryon(thePDGCode)) theXsec=2*P*millibarn;
    if(CustomPDGParser::s_isRMeson(thePDGCode)||CustomPDGParser::s_isRGlueball(thePDGCode)) theXsec=(R+2*P)*millibarn;
    if(CustomPDGParser::s_isRBaryon(thePDGCode)) theXsec=3*P*millibarn;
      }
    else
      {
    if(CustomPDGParser::s_isMesonino(thePDGCode)) theXsec=P*millibarn;
    if(CustomPDGParser::s_isSbaryon(thePDGCode)) theXsec=(2*(P+R)+30/sqrt(boost))*millibarn;
    if(CustomPDGParser::s_isRMeson(thePDGCode)||CustomPDGParser::s_isRGlueball(thePDGCode)) theXsec=(R+2*P)*millibarn;
    if(CustomPDGParser::s_isRBaryon(thePDGCode)) theXsec=3*P*millibarn;
      }
  }
  //Adding resonance
  if(resonant)
    {
      double e_0 = ek_0 + aParticle->GetDefinition()->GetPDGMass(); //Now total energy

      e_0 = sqrt(aParticle->GetDefinition()->GetPDGMass()*aParticle->GetDefinition()->GetPDGMass()
         + theProton->GetPDGMass()*theProton->GetPDGMass()
         + 2.*e_0*theProton->GetPDGMass());
      //      edm::LogInfo("SimG4CoreCustomPhysics")<<e_0/GeV<<G4endl;
      //      edm::LogInfo("SimG4CoreCustomPhysics")<<ek_0/GeV<<" "<<aParticle->GetDefinition()->GetPDGMass()/GeV<<" "<<theProton->GetPDGMass()/GeV<<G4endl;
      double sqrts=sqrt(aParticle->GetDefinition()->GetPDGMass()*aParticle->GetDefinition()->GetPDGMass()
            + theProton->GetPDGMass()*theProton->GetPDGMass() + 2*aParticle->GetTotalEnergy()*theProton->GetPDGMass());

      double res_result = amplitude*(gamma*gamma/4.)/((sqrts-e_0)*(sqrts-e_0)+(gamma*gamma/4.));//Non-relativistic Breit Wigner

      theXsec += res_result;
      //      if(fabs(aParticle->GetKineticEnergy()/GeV-200)<10)  std::cout<<sqrts/GeV<<" "<<theXsec/millibarn<<std::endl;
    }

  //  std::cout<<"Xsec/nucleon: "<<theXsec/millibarn<<"millibarn, Total Xsec: "<<theXsec * anElement->GetN()/millibarn<<" millibarn"<<std::endl;
  return theXsec * pow(anElement->GetN(),0.7)*1.25;// * 0.523598775598299;
}

G4ProcessHelper& G4ProcessHelper::operator= ( const G4ProcessHelper &  )

privatedelete

G4double G4ProcessHelper::PhaseSpace ( const ReactionProduct &  aReaction, const G4DynamicParticle *  aDynamicParticle  )

private

Definition at line 438 of file G4ProcessHelper.cc.

References GeV, funct::pow(), and mathSSE::sqrt().

                                                                                                              {
  G4double qValue = ReactionProductMass(aReaction,aDynamicParticle);
  G4double phi = sqrt(1+qValue/(2*0.139*GeV))*pow(qValue/(1.1*GeV),3./2.);
  return (phi/(1+phi));
}

double G4ProcessHelper::Pom ( const double  boost )

private

Definition at line 477 of file G4ProcessHelper.cc.

References a, b, EnergyCorrector::c, edmIntegrityCheck::d, funct::pow(), and mathSSE::sqrt().

{
  double a=4.138224000651535;
  double b=1.50377557581421;
  double c=-0.05449742257808247;
  double d=0.0008221235048211401;
  return a + b*sqrt(boost) + c*boost + d*pow(boost,1.5);
}

G4bool G4ProcessHelper::ReactionGivesBaryon ( const ReactionProduct &  aReaction )

private

Definition at line 432 of file G4ProcessHelper.cc.

References CustomPDGParser::s_isRBaryon(), and CustomPDGParser::s_isSbaryon().

                                                                           {
  for (ReactionProduct::const_iterator it = aReaction.begin();it!=aReaction.end();it++)
    if(CustomPDGParser::s_isSbaryon(*it)||CustomPDGParser::s_isRBaryon(*it)) return true;
  return false;
}

G4bool G4ProcessHelper::ReactionIsPossible ( const ReactionProduct &  aReaction, const G4DynamicParticle *  aDynamicParticle  )

private

Definition at line 427 of file G4ProcessHelper.cc.

                                                                                                                    {
  if (ReactionProductMass(aReaction,aDynamicParticle)>0) return true;
  return false;
}

G4double G4ProcessHelper::ReactionProductMass ( const ReactionProduct &  aReaction, const G4DynamicParticle *  aDynamicParticle  )

private

Definition at line 407 of file G4ProcessHelper.cc.

References mathSSE::sqrt().

                                                                                                                       {
  // Incident energy:
  G4double E_incident = aDynamicParticle->GetTotalEnergy();
  //edm::LogInfo("SimG4CoreCustomPhysics")<<"Total energy: "<<E_incident<<" Kinetic: "<<aDynamicParticle->GetKineticEnergy()<<G4endl;
  // sqrt(s)= sqrt(m_1^2 + m_2^2 + 2 E_1 m_2)
  G4double m_1 = aDynamicParticle->GetDefinition()->GetPDGMass();
  G4double m_2 = theTarget->GetPDGMass();
  //edm::LogInfo("SimG4CoreCustomPhysics")<<"M_R: "<<m_1/GeV<<" GeV, M_np: "<<m_2/GeV<<" GeV"<<G4endl;
  G4double sqrts = sqrt(m_1*m_1 + m_2*(m_2 + 2 * E_incident));
  //edm::LogInfo("SimG4CoreCustomPhysics")<<"sqrt(s) = "<<sqrts/GeV<<" GeV"<<G4endl;
  // Sum of rest masses after reaction:
  G4double M_after = 0;
  for (ReactionProduct::const_iterator r_it = aReaction.begin(); r_it !=aReaction.end(); r_it++){
    //edm::LogInfo("SimG4CoreCustomPhysics")<<"Mass contrib: "<<(particleTable->FindParticle(*r_it)->GetPDGMass())/MeV<<" MeV"<<G4endl;
    M_after += particleTable->FindParticle(*r_it)->GetPDGMass();
  }
  //edm::LogInfo("SimG4CoreCustomPhysics")<<"Intending to return this ReactionProductMass: "<<(sqrts - M_after)/MeV<<" MeV"<<G4endl;
  return sqrts - M_after;
}

void G4ProcessHelper::ReadAndParse ( const G4String &  str, std::vector< G4String > &  tokens, const G4String &  delimiters = " "  )

private

Definition at line 444 of file G4ProcessHelper.cc.

References groupFilesInBlocks::temp.

{
  // Skip delimiters at beginning.
  G4String::size_type lastPos = str.find_first_not_of(delimiters, 0);
  // Find first "non-delimiter".
  G4String::size_type pos     = str.find_first_of(delimiters, lastPos);
  
  while (G4String::npos != pos || G4String::npos != lastPos)
    {
      //Skipping leading / trailing whitespaces
      G4String temp = str.substr(lastPos, pos - lastPos);
      while(temp.c_str()[0] == ' ') temp.erase(0,1);
      while(temp[temp.size()-1] == ' ') temp.erase(temp.size()-1,1);
      // Found a token, add it to the vector.
      tokens.push_back(temp);
      // Skip delimiters.  Note the "not_of"
      lastPos = str.find_first_not_of(delimiters, pos);
      // Find next "non-delimiter"
      pos = str.find_first_of(delimiters, lastPos);
    }
}

double G4ProcessHelper::Regge ( const double  boost )

private

Definition at line 468 of file G4ProcessHelper.cc.

References a, b, EnergyCorrector::c, JetChargeProducer_cfi::exp, and cmsBatch::log.

{
  double a=2.165635078566177;
  double b=0.1467453738547229;
  double c=-0.9607903711871166;
  return 1.5*exp(a+b/boost+c*log(boost));
}

Member Data Documentation

double G4ProcessHelper::amplitude

private

Definition at line 84 of file G4ProcessHelper.h.

G4double G4ProcessHelper::checkfraction

private

Definition at line 50 of file G4ProcessHelper.h.

double G4ProcessHelper::ek_0

private

Definition at line 82 of file G4ProcessHelper.h.

CustomParticleFactory* G4ProcessHelper::fParticleFactory

private

Definition at line 97 of file G4ProcessHelper.h.

double G4ProcessHelper::gamma

private

Definition at line 83 of file G4ProcessHelper.h.

TProfile* G4ProcessHelper::h_q_gamma

private

Definition at line 104 of file G4ProcessHelper.h.

TProfile* G4ProcessHelper::h_q_p

private

Definition at line 103 of file G4ProcessHelper.h.

TH1D* G4ProcessHelper::h_sqrts

private

Definition at line 102 of file G4ProcessHelper.h.

TProfile* G4ProcessHelper::h_xsec_cms

private

Definition at line 101 of file G4ProcessHelper.h.

TProfile* G4ProcessHelper::h_xsec_lab

private

Definition at line 100 of file G4ProcessHelper.h.

double G4ProcessHelper::hadronlifetime

private

Definition at line 87 of file G4ProcessHelper.h.

std::map<const G4ParticleDefinition*,G4bool> G4ProcessHelper::known_particles

private

Definition at line 75 of file G4ProcessHelper.h.

double G4ProcessHelper::mixing

private

Definition at line 88 of file G4ProcessHelper.h.

G4int G4ProcessHelper::n_22

private

Definition at line 51 of file G4ProcessHelper.h.

G4int G4ProcessHelper::n_23

private

Definition at line 52 of file G4ProcessHelper.h.

ReactionMap G4ProcessHelper::nReactionMap

private

Definition at line 95 of file G4ProcessHelper.h.

std::map<G4String,G4double> G4ProcessHelper::parameters

private

Definition at line 78 of file G4ProcessHelper.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

G4ParticleTable* G4ProcessHelper::particleTable

private

Definition at line 98 of file G4ProcessHelper.h.

ReactionMap G4ProcessHelper::pReactionMap

private

Definition at line 92 of file G4ProcessHelper.h.

bool G4ProcessHelper::reggemodel

private

Definition at line 86 of file G4ProcessHelper.h.

bool G4ProcessHelper::resonant

private

Definition at line 81 of file G4ProcessHelper.h.

double G4ProcessHelper::suppressionfactor

private

Definition at line 85 of file G4ProcessHelper.h.

HistoHelper* G4ProcessHelper::theHistoHelper

private

Definition at line 99 of file G4ProcessHelper.h.

G4ParticleDefinition* G4ProcessHelper::theNeutron

private

Definition at line 56 of file G4ProcessHelper.h.

G4ParticleDefinition* G4ProcessHelper::theProton

private

Definition at line 55 of file G4ProcessHelper.h.

G4ParticleDefinition* G4ProcessHelper::theRbaryoncloud

private

Definition at line 58 of file G4ProcessHelper.h.

ReactionMap* G4ProcessHelper::theReactionMap

private

Definition at line 60 of file G4ProcessHelper.h.

G4ParticleDefinition* G4ProcessHelper::theRmesoncloud

private

Definition at line 57 of file G4ProcessHelper.h.

G4ParticleDefinition* G4ProcessHelper::theTarget

private

Definition at line 54 of file G4ProcessHelper.h.