97 |
|
) |
98 |
|
) |
99 |
|
) { |
100 |
< |
if (debug_ > 0) { |
100 |
> |
if (debug_ > 1) { |
101 |
|
std::cout << "** initial parton #" << imc ; |
102 |
|
printParticle(lPart); |
103 |
|
} |
106 |
|
else n_antiquarks++; |
107 |
|
p_partonFlavour->Fill(lPart.pdgId()); |
108 |
|
} |
109 |
< |
else if ((debug_ > 0 && imc < 15) || debug_ > 1) { |
109 |
> |
else if ((debug_ > 1 && imc < 15) || debug_ > 2) { |
110 |
|
std::cout << "** not selected: #" << imc; |
111 |
|
printParticle(lPart); |
112 |
|
} |
136 |
|
nTotal_ += nPartons; |
137 |
|
assert (nPartons == n_quarks+n_antiquarks); |
138 |
|
|
139 |
< |
if (flavour_ > 3) { |
140 |
< |
assert(n_quarks >= 1 && n_antiquarks >= 1); |
141 |
< |
assert (nPartons == 2); |
142 |
< |
} |
139 |
> |
// if (flavour_ > 3) { |
140 |
> |
// assert(n_quarks >= 1 && n_antiquarks >= 1); |
141 |
> |
// assert (nPartons == 2); |
142 |
> |
// } |
143 |
|
|
144 |
|
p_nQuarks->Fill(n_quarks); |
145 |
|
p_nAntiQuarks->Fill(n_antiquarks); |