Skip to content

Commit

Permalink
Remove PauliBlocker from the event generation chain for Charm-CCQE
Browse files Browse the repository at this point in the history
  • Loading branch information
@sjgardiner
sjgardiner committed Mar 22, 2019
1 parent fae39cc commit fb7835a
Show file tree
Hide file tree
Showing 5 changed files with 18 additions and 21 deletions.
9 changes: 4 additions & 5 deletions config/EventGenerator.xml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -347,17 +347,16 @@ XSecModel alg Yes Cross section model used at the thread

<param_set name="QEL-CC-CHARM">
<param type="string" name="VldContext"> </param>
<param type="int" name="NModules"> 10 </param>
<param type="int" name="NModules"> 9 </param>
<param type="alg" name="Module-0"> genie::InitialStateAppender/Default </param>
<param type="alg" name="Module-1"> genie::VertexGenerator/Default </param>
<param type="alg" name="Module-2"> genie::FermiMover/Default </param>
<param type="alg" name="Module-3"> genie::QELKinematicsGenerator/CC-Charm-Default </param>
<param type="alg" name="Module-4"> genie::QELPrimaryLeptonGenerator/Default </param>
<param type="alg" name="Module-5"> genie::QELHadronicSystemGenerator/Default </param>
<param type="alg" name="Module-6"> genie::PauliBlocker/Default </param>
<param type="alg" name="Module-7"> genie::NucDeExcitationSim/Default </param>
<param type="alg" name="Module-8"> genie::HadronTransporter/Default </param>
<param type="alg" name="Module-9"> genie::UnstableParticleDecayer/AfterHadronTransport </param>
<param type="alg" name="Module-6"> genie::NucDeExcitationSim/Default </param>
<param type="alg" name="Module-7"> genie::HadronTransporter/Default </param>
<param type="alg" name="Module-8"> genie::UnstableParticleDecayer/AfterHadronTransport </param>
<param type="alg" name="ILstGen"> genie::QELInteractionListGenerator/CC-Charm-Default </param>
</param_set>

Expand Down
9 changes: 4 additions & 5 deletions config/G00_00a/EventGenerator.xml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -254,17 +254,16 @@ XSecModel alg Yes Cross section model used at the thread

<param_set name="QEL-CC-CHARM">
<param type="string" name="VldContext"> </param>
<param type="int" name="NModules"> 10 </param>
<param type="int" name="NModules"> 9 </param>
<param type="alg" name="Module-0"> genie::InitialStateAppender/Default </param>
<param type="alg" name="Module-1"> genie::VertexGenerator/Default </param>
<param type="alg" name="Module-2"> genie::FermiMover/Default </param>
<param type="alg" name="Module-3"> genie::QELKinematicsGenerator/CC-Charm-Default </param>
<param type="alg" name="Module-4"> genie::QELPrimaryLeptonGenerator/Default </param>
<param type="alg" name="Module-5"> genie::QELHadronicSystemGenerator/Default </param>
<param type="alg" name="Module-6"> genie::PauliBlocker/Default </param>
<param type="alg" name="Module-7"> genie::NucDeExcitationSim/Default </param>
<param type="alg" name="Module-8"> genie::HadronTransporter/Default </param>
<param type="alg" name="Module-9"> genie::UnstableParticleDecayer/AfterHadronTransport </param>
<param type="alg" name="Module-6"> genie::NucDeExcitationSim/Default </param>
<param type="alg" name="Module-7"> genie::HadronTransporter/Default </param>
<param type="alg" name="Module-8"> genie::UnstableParticleDecayer/AfterHadronTransport </param>
<param type="alg" name="ILstGen"> genie::QELInteractionListGenerator/CC-Charm-Default </param>
</param_set>

Expand Down
9 changes: 4 additions & 5 deletions config/G00_00b/EventGenerator.xml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -254,17 +254,16 @@ XSecModel alg Yes Cross section model used at the thread

<param_set name="QEL-CC-CHARM">
<param type="string" name="VldContext"> </param>
<param type="int" name="NModules"> 10 </param>
<param type="int" name="NModules"> 9 </param>
<param type="alg" name="Module-0"> genie::InitialStateAppender/Default </param>
<param type="alg" name="Module-1"> genie::VertexGenerator/Default </param>
<param type="alg" name="Module-2"> genie::FermiMover/Default </param>
<param type="alg" name="Module-3"> genie::QELKinematicsGenerator/CC-Charm-Default </param>
<param type="alg" name="Module-4"> genie::QELPrimaryLeptonGenerator/Default </param>
<param type="alg" name="Module-5"> genie::QELHadronicSystemGenerator/Default </param>
<param type="alg" name="Module-6"> genie::PauliBlocker/Default </param>
<param type="alg" name="Module-7"> genie::NucDeExcitationSim/Default </param>
<param type="alg" name="Module-8"> genie::HadronTransporter/Default </param>
<param type="alg" name="Module-9"> genie::UnstableParticleDecayer/AfterHadronTransport </param>
<param type="alg" name="Module-6"> genie::NucDeExcitationSim/Default </param>
<param type="alg" name="Module-7"> genie::HadronTransporter/Default </param>
<param type="alg" name="Module-8"> genie::UnstableParticleDecayer/AfterHadronTransport </param>
<param type="alg" name="ILstGen"> genie::QELInteractionListGenerator/CC-Charm-Default </param>
</param_set>

Expand Down
9 changes: 4 additions & 5 deletions config/G00_00z/EventGenerator.xml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -254,17 +254,16 @@ XSecModel alg Yes Cross section model used at the thread

<param_set name="QEL-CC-CHARM">
<param type="string" name="VldContext"> </param>
<param type="int" name="NModules"> 10 </param>
<param type="int" name="NModules"> 9 </param>
<param type="alg" name="Module-0"> genie::InitialStateAppender/Default </param>
<param type="alg" name="Module-1"> genie::VertexGenerator/Default </param>
<param type="alg" name="Module-2"> genie::FermiMover/Default </param>
<param type="alg" name="Module-3"> genie::QELKinematicsGenerator/CC-Charm-Default </param>
<param type="alg" name="Module-4"> genie::QELPrimaryLeptonGenerator/Default </param>
<param type="alg" name="Module-5"> genie::QELHadronicSystemGenerator/Default </param>
<param type="alg" name="Module-6"> genie::PauliBlocker/Default </param>
<param type="alg" name="Module-7"> genie::NucDeExcitationSim/Default </param>
<param type="alg" name="Module-8"> genie::HadronTransporter/Default </param>
<param type="alg" name="Module-9"> genie::UnstableParticleDecayer/AfterHadronTransport </param>
<param type="alg" name="Module-6"> genie::NucDeExcitationSim/Default </param>
<param type="alg" name="Module-7"> genie::HadronTransporter/Default </param>
<param type="alg" name="Module-8"> genie::UnstableParticleDecayer/AfterHadronTransport </param>
<param type="alg" name="ILstGen"> genie::QELInteractionListGenerator/CC-Charm-Default </param>
</param_set>

Expand Down
3 changes: 2 additions & 1 deletion src/Physics/NuclearState/PauliBlocker.cxx
View file
Original file line number Diff line number Diff line change
Expand Up @@ -168,7 +168,8 @@ void PauliBlocker::LoadModelType(void){
double PauliBlocker::GetFermiMomentum(const Target& tgt, int pdg_Nf,
double radius) const
{
assert(pdg::IsProton(pdg_Nf) || pdg::IsNeutron(pdg_Nf));
// Pauli blocking should only be applied for nucleons
assert( pdg::IsProton(pdg_Nf) || pdg::IsNeutron(pdg_Nf) );
double kF = 0.;
if ( fLFG ) {
int A = tgt.A();
Expand Down

0 comments on commit fb7835a

Please sign in to comment.