Enable dipole subtraction for processes with decays into coloured particles

[marek]

Amegic currently cannot build subtraction terms for processes with decays into coloured particles. As far as I can see, the problem is that when determining the colour indices of the insertion operators, the decays are not accounted for. Attached is a relevant run-card for ttbar production. Steffen, could you please have a look?

(run){
  EVENTS 100;
  BATCH_MODE 0; OUTPUT 2;
  SHOWER_GENERATOR None;
  FRAGMENTATION Off;
}(run);

(processes){
  Process 93 93 -> 6[a] -6[b];
    Decay 6[a] -> 5 -11 12;  
    Decay -6[b] -> -5 13 -14;
  NLO_QCD_Part BIRS;
  Order_EW 0;
  Print_Graphs 1;
  End process;
}(processes);

(model){
  ACTIVE[25] 0; 
  MASSIVE[5] 1; 
}(model)

(me){
  ME_SIGNAL_GENERATOR Amegic;      
  EVENT_GENERATION_MODE Weighted;
  NLO_Mode 1;
  SCALES VAR{Abs2(p[0]+p[1])};
}(me);

(beam){
  BEAM_1  2212; BEAM_ENERGY_1 3500;
  BEAM_2  2212; BEAM_ENERGY_2 3500;
  K_PERP_MEAN_1 0.; K_PERP_SIGMA_1 0.;
  K_PERP_MEAN_2 0.; K_PERP_SIGMA_2 0.;
}(beam);

[Steffen Schumann]

I will look into that asap. Steffen