diff --git a/src/snippets/Diffraction.in b/src/snippets/Diffraction.in --- a/src/snippets/Diffraction.in +++ b/src/snippets/Diffraction.in @@ -1,114 +1,120 @@ ################################################## # Create separate SubProcessHandler for Diffraction ################################################## cd /Herwig/MatrixElements create Herwig::MEDiffraction MEDiffractionLeft set MEDiffractionLeft:DiffDirection Left create Herwig::MEDiffraction MEDiffractionRight set MEDiffractionRight:DiffDirection Right create Herwig::MEDiffraction MEDiffractionDouble set MEDiffractionDouble:DiffDirection Both create Herwig::MEDiffraction MEDiffractionDeltaLeft set MEDiffractionDeltaLeft:DiffDirection Left create Herwig::MEDiffraction MEDiffractionDeltaRight set MEDiffractionDeltaRight:DiffDirection Right create Herwig::MEDiffraction MEDiffractionDeltaDouble set MEDiffractionDeltaDouble:DiffDirection Both # Make a parton extractor for diffraction cd /Herwig/Partons cp PPExtractor DiffPPExtractor set DiffPPExtractor:FirstPDF /Herwig/Partons/NoPDF set DiffPPExtractor:SecondPDF /Herwig/Partons/NoPDF cd /Herwig/MatrixElements/ # Create Diffraction SubProcessHandler cp SubProcess QCDDiffraction # Assign the PartonExtractor to the SubProcessHandler set QCDDiffraction:PartonExtractor /Herwig/Partons/DiffPPExtractor # Use only Delta as final excited state (Yes/No) set MEDiffractionLeft:DeltaOnly No set MEDiffractionRight:DeltaOnly No set MEDiffractionDouble:DeltaOnly No set MEDiffractionDeltaLeft:DeltaOnly Yes set MEDiffractionDeltaRight:DeltaOnly Yes set MEDiffractionDeltaDouble:DeltaOnly Yes # Set weight for Diffraction set MEDiffractionLeft:DiffractionAmplitude 12 set MEDiffractionRight:DiffractionAmplitude 12 set MEDiffractionDouble:DiffractionAmplitude 8 set MEDiffractionDeltaLeft:DiffractionAmplitude 4 set MEDiffractionDeltaRight:DiffractionAmplitude 4 set MEDiffractionDeltaDouble:DiffractionAmplitude 2 # Set soft diffraction paramters # Parameter values from arxiv/0709.0395 set MEDiffractionLeft:ProtonPomeronSlope 10.1 set MEDiffractionLeft:SoftPomeronIntercept 1.08 set MEDiffractionLeft:SoftPomeronSlope 0.25 set MEDiffractionRight:ProtonPomeronSlope 10.1 set MEDiffractionRight:SoftPomeronIntercept 1.08 set MEDiffractionRight:SoftPomeronSlope 0.25 set MEDiffractionDouble:ProtonPomeronSlope 10.1 set MEDiffractionDouble:SoftPomeronIntercept 1.08 set MEDiffractionDouble:SoftPomeronSlope 0.25 set MEDiffractionDeltaLeft:ProtonPomeronSlope 10.1 set MEDiffractionDeltaLeft:SoftPomeronIntercept 1.08 set MEDiffractionDeltaLeft:SoftPomeronSlope 0.25 set MEDiffractionDeltaRight:ProtonPomeronSlope 10.1 set MEDiffractionDeltaRight:SoftPomeronIntercept 1.08 set MEDiffractionDeltaRight:SoftPomeronSlope 0.25 set MEDiffractionDeltaDouble:ProtonPomeronSlope 10.1 set MEDiffractionDeltaDouble:SoftPomeronIntercept 1.08 set MEDiffractionDeltaDouble:SoftPomeronSlope 0.25 # Set number of clusters for dissociation set MEDiffractionLeft:DissociationDecay One set MEDiffractionRight:DissociationDecay One set MEDiffractionDouble:DissociationDecay One set MEDiffractionDeltaLeft:DissociationDecay One set MEDiffractionDeltaRight:DissociationDecay One set MEDiffractionDeltaDouble:DissociationDecay One # Insert matrix elements insert QCDDiffraction:MatrixElements[0] MEDiffractionLeft insert QCDDiffraction:MatrixElements[0] MEDiffractionRight insert QCDDiffraction:MatrixElements[0] MEDiffractionDouble insert QCDDiffraction:MatrixElements[0] MEDiffractionDeltaLeft insert QCDDiffraction:MatrixElements[0] MEDiffractionDeltaRight #insert QCDDiffraction:MatrixElements[0] MEDiffractionDeltaDouble # No cluster fission (set to E/2, if needed. Here E=7TeV.) #set /Herwig/Hadronization/ClusterFissioner:ClMaxLight 3500 # The DiffractionMEs need to know the MPIHandler for the cross section. set MEDiffractionLeft:MPIHandler /Herwig/UnderlyingEvent/MPIHandler set MEDiffractionRight:MPIHandler /Herwig/UnderlyingEvent/MPIHandler set MEDiffractionDouble:MPIHandler /Herwig/UnderlyingEvent/MPIHandler set MEDiffractionDeltaLeft:MPIHandler /Herwig/UnderlyingEvent/MPIHandler set MEDiffractionDeltaRight:MPIHandler /Herwig/UnderlyingEvent/MPIHandler set MEDiffractionDeltaDouble:MPIHandler /Herwig/UnderlyingEvent/MPIHandler set QCDDiffraction:CascadeHandler NULL set /Herwig/Generators/EventGenerator:EventHandler:CascadeHandler NULL insert /Herwig/Generators/EventGenerator:EventHandler:SubProcessHandlers[0] QCDDiffraction + + +## These cuts are also set in MB.in +# If only diffraction is needed, we want the same cuts: +set EventGenerator:EventHandler:Cuts /Herwig/Cuts/MinBiasCuts + ########################################################################## ##########################################################################