|   | 1 |  | 
          
          
            |   | 2 |  | 
          
          
            |   | 3 | === How do I implement my own renormalization/factorization scales in !MadEvent? === | 
          
          
            |   | 4 |  | 
          
          
            |   | 5 | Since v.4.2.17, there is a new default dynamical scale choice in [http://cp3wks05.fynu.ucl.ac.be/cgi-bin/ViewVC/viewvc.cgi/MG_ME/Template/Cards/run_card.dat?view=markup run_card.dat] file:  | 
          
          
            |   | 6 |  | 
          
          
            |   | 7 | {{{ | 
          
          
            |   | 8 | #********************************************************************* | 
          
          
            |   | 9 | # Renormalization and factorization scales                           * | 
          
          
            |   | 10 | #********************************************************************* | 
          
          
            |   | 11 |  F        = fixed_ren_scale  ! if .true. use fixed ren scale | 
          
          
            |   | 12 |  F        = fixed_fac_scale  ! if .true. use fixed fac scale | 
          
          
            |   | 13 |  91.1880  = scale            ! fixed ren scale | 
          
          
            |   | 14 |  91.1880  = dsqrt_q2fact1    ! fixed fact scale for pdf1 | 
          
          
            |   | 15 |  91.1880  = dsqrt_q2fact2    ! fixed fact scale for pdf2 | 
          
          
            |   | 16 |  1        = scalefact        ! scale factor for event-by-event scales | 
          
          
            |   | 17 | }}} | 
          
          
            |   | 18 |  | 
          
          
            |   | 19 |  | 
          
          
            |   | 20 | i.e., variable factorization and renormalization scales,set to the central mT^2 scale after kT-clustering of the event. | 
          
          
            |   | 21 | For a single heavy particle, this corresponds to the M^2+pT^2,  | 
          
          
            |   | 22 | for pair of heavy particles it is the geometric mean of M^2+pT^2 for each particle,  | 
          
          
            |   | 23 | for massless particles it is the pT^2 of the last pair after clustering.  | 
          
          
            |   | 24 | If the last clustering is a t-channel colorless exchange (as, e.g., in WBF),  the two fact. scales are set to the last mT^2 values on either side.  | 
          
          
            |   | 25 | These scales are overridden by scales defined in [http://cp3wks05.fynu.ucl.ac.be/cgi-bin/ViewVC/viewvc.cgi/MG_ME/Template/SubProcesses/setscales.f?view=markup setscales.f], if they are returned as non-0. | 
          
          
            |   | 26 |  | 
          
          
            |   | 27 | Fixed values for the scale (not event-by-event)  can be set in the same file, by setting the two flags {{{ fixed_ren_scale, fixed_fac_scale}}}  to True and choosing the desired value for them : {{{ scale, dsqrt_q2fact1, dsqrt_q2fact2 }}} | 
          
          
            |   | 28 |  | 
          
          
            |   | 29 |  | 
          
          
            |   | 30 | If the user wants to have a different choice for the dynamical scale, one needs | 
          
          
            |   | 31 |  | 
          
          
            |   | 32 |    1. to know a bit of fortran | 
          
          
            |   | 33 |    2. to download the code (either after process generation or from our Download page) and run on your local machine.  | 
          
          
            |   | 34 |  | 
          
          
            |   | 35 | The file to be edited is  [http://cp3wks05.fynu.ucl.ac.be/cgi-bin/ViewVC/viewvc.cgi/MG_ME/Template/SubProcesses/setscales.f?view=markup SubProcesses/setscales.f]. | 
          
          
            |   | 36 |  | 
          
          
            |   | 37 |  | 
          
          
            |   | 38 |  | 
          
          
            |   | 39 | -- Main.FabioMaltoni - 02 Mar 2009 |