ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Mods/src/GeneratorMod.cc
Revision: 1.28
Committed: Tue Mar 24 16:13:21 2009 UTC (16 years, 1 month ago) by loizides
Content type: text/plain
Branch: MAIN
CVS Tags: Mit_008
Changes since 1.27: +19 -4 lines
Log Message: 
Added Si Changes for Madgraph.

File Contents

# User Rev Content
1 loizides 1.28 // $Id: GeneratorMod.cc,v 1.27 2009/02/13 14:33:31 loizides Exp $
2 loizides 1.1
3     #include "MitPhysics/Mods/interface/GeneratorMod.h"
4     #include "MitCommon/MathTools/interface/MathUtils.h"
5 loizides 1.8 #include "MitPhysics/Init/interface/ModNames.h"
6 loizides 1.1 #include <TH1D.h>
7     #include <TH2D.h>
8    
9     using namespace mithep;
10    
11     ClassImp(mithep::GeneratorMod)
12    
13     //--------------------------------------------------------------------------------------------------
14 loizides 1.8 GeneratorMod::GeneratorMod(const char *name, const char *title) :
15 loizides 1.1 BaseMod(name,title),
16     fMCPartName(Names::gkMCPartBrn),
17 loizides 1.8 fMCLeptonsName(ModNames::gkMCLeptonsName),
18     fMCAllLeptonsName(ModNames::gkMCAllLeptonsName),
19     fMCTausName(ModNames::gkMCTausName),
20     fMCNeutrinosName(ModNames::gkMCNeutrinosName),
21     fMCQuarksName(ModNames::gkMCQuarksName),
22     fMCqqHsName(ModNames::gkMCqqHsName),
23     fMCBosonsName(ModNames::gkMCBosonsName),
24 ceballos 1.10 fMCPhotonsName(ModNames::gkMCPhotonsName),
25 ceballos 1.14 fPtLeptonMin(0.0),
26     fEtaLeptonMax(5.0),
27     fPtPhotonMin(0.0),
28 loizides 1.15 fEtaPhotonMax(5.0),
29 loizides 1.27 fPdgIdCut(0),
30     fMassMinCut(-FLT_MAX),
31     fMassMaxCut(FLT_MAX),
32 loizides 1.15 fParticles(0)
33 loizides 1.1 {
34     // Constructor.
35     }
36    
37     //--------------------------------------------------------------------------------------------------
38     void GeneratorMod::Process()
39     {
40 loizides 1.5 // Process entries of the tree.
41 loizides 1.1
42 loizides 1.5 // these arrays will be filled in the loop of particles
43 loizides 1.9 MCParticleOArr *GenLeptons = new MCParticleOArr;
44 loizides 1.13 GenLeptons->SetName(fMCLeptonsName);
45 loizides 1.9 MCParticleOArr *GenAllLeptons = new MCParticleOArr;
46 loizides 1.13 GenAllLeptons->SetName(fMCAllLeptonsName);
47 loizides 1.9 MCParticleOArr *GenTaus = new MCParticleOArr;
48 loizides 1.13 GenTaus->SetName(fMCTausName);
49     GenTaus->SetOwner(kTRUE);
50 loizides 1.9 MCParticleOArr *GenNeutrinos = new MCParticleOArr;
51 loizides 1.13 GenNeutrinos->SetName(fMCNeutrinosName);
52 loizides 1.9 MCParticleOArr *GenQuarks = new MCParticleOArr;
53 loizides 1.13 GenQuarks->SetName(fMCQuarksName);
54 loizides 1.9 MCParticleOArr *GenqqHs = new MCParticleOArr;
55 loizides 1.13 GenqqHs->SetName(fMCqqHsName);
56 loizides 1.9 MCParticleOArr *GenBosons = new MCParticleOArr;
57 loizides 1.13 GenBosons->SetName(fMCBosonsName);
58 ceballos 1.10 MCParticleOArr *GenPhotons = new MCParticleOArr;
59 loizides 1.13 GenPhotons->SetName(fMCPhotonsName);
60 loizides 1.1
61 loizides 1.5 // load MCParticle branch
62     LoadBranch(fMCPartName);
63    
64 loizides 1.6 Bool_t isqqH = kFALSE;
65 loizides 1.5 for (UInt_t i=0; i<fParticles->GetEntries(); ++i) {
66 loizides 1.7 const MCParticle *p = fParticles->At(i);
67 loizides 1.5
68     if (!p->IsGenerated()) continue;
69    
70     // muons/electrons from W/Z decays
71 loizides 1.7 if ((p->Is(MCParticle::kEl) || p->Is(MCParticle::kMu)) && p->Status() == 1) {
72 ceballos 1.14 if (p->Pt() > fPtLeptonMin && p->AbsEta() < fEtaLeptonMax) {
73 loizides 1.5 GenAllLeptons->Add(p);
74     }
75 loizides 1.6 Bool_t isGoodLepton = kFALSE;
76 loizides 1.7 const MCParticle *pm = p;
77 loizides 1.5 while (pm->HasMother() && isGoodLepton == kFALSE) {
78 loizides 1.21 if (pm->PdgId() == 92) // string reached, terminate loop
79     break;
80 loizides 1.7 if (pm->Mother()->Is(MCParticle::kZ) || pm->Mother()->Is(MCParticle::kW)) {
81 loizides 1.5 GenLeptons->Add(p);
82     isGoodLepton = kTRUE;
83 loizides 1.7 break;
84 loizides 1.8 } else if (pm->Mother()->Is(MCParticle::kPi0) || pm->Mother()->Is(MCParticle::kEta)) {
85 loizides 1.7 // this is fake, but it is a trick to get rid of these cases and abort the loop
86     break;
87     }
88     pm = pm->Mother();
89 loizides 1.1 }
90 loizides 1.5 }
91    
92 loizides 1.7 // hadronic taus
93     else if (p->Is(MCParticle::kTau) && p->Status() == 2) {
94     if (!p->HasDaughter(MCParticle::kEl) && !p->HasDaughter(MCParticle::kMu)) {
95     const MCParticle *tv = p->FindDaughter(MCParticle::kTauNu);
96     if (tv) {
97     MCParticle *pm_f = new MCParticle(*p);
98     pm_f->SetMom(p->Px()-tv->Px(), p->Py()-tv->Py(),
99     p->Pz()-tv->Pz(), p->E()-tv->E());
100     GenTaus->AddOwned(pm_f);
101     } else {
102 loizides 1.18 SendError(kWarning, "Process", "Could not find a tau neutrino!");
103 loizides 1.5 }
104 loizides 1.1 }
105 loizides 1.5 }
106 loizides 1.1
107 loizides 1.5 // neutrinos
108 loizides 1.7 else if (p->Status() == 1 && p->IsNeutrino()) {
109 loizides 1.5 GenNeutrinos->Add(p);
110     }
111 loizides 1.1
112 loizides 1.5 // quarks from W/Z decays or top particles
113 loizides 1.7 else if (p->IsQuark() && p->HasMother()) {
114     if (p->Mother()->Is(MCParticle::kZ) || p->Mother()->Is(MCParticle::kW) ||
115     p->Is(MCParticle::kTop) || p->Mother()->Is(MCParticle::kTop)) {
116 loizides 1.5 GenQuarks->Add(p);
117 loizides 1.1 }
118 loizides 1.5 }
119 loizides 1.1
120 loizides 1.5 // qqH, information about the forward jets
121 loizides 1.8 else if (isqqH == kFALSE && p->Is(MCParticle::kH)) {
122 loizides 1.5 isqqH = kTRUE;
123 loizides 1.17 const MCParticle *pq1 = fParticles->At(i-1);
124     const MCParticle *pq2 = fParticles->At(i-2);
125 loizides 1.7 if (!pq1 || !pq2) {
126     SendError(kWarning, "Process", "Could not find quark pair!");
127 loizides 1.8 } else if (pq1->IsQuark() && pq2->IsQuark() &&
128     pq1->HasMother() && pq2->HasMother() &&
129 ceballos 1.25 pq1->Mother() == pq2->Mother() &&
130     pq1->Mother() != p->Mother() &&
131     pq1->Mother() != p->Mother()) {
132 loizides 1.5 GenqqHs->Add(pq1);
133     GenqqHs->Add(pq2);
134 loizides 1.1 }
135 ceballos 1.23
136 loizides 1.7 if (p->Status() == 3)
137 loizides 1.24 GenBosons->Add(p); // take higgs boson in account here rather in next else if
138 loizides 1.5 }
139 loizides 1.1
140 loizides 1.5 // information about bosons: W, Z, h, Z', W', H0, A0, H+
141 loizides 1.7 else if (p->Status() == 3 &&
142     (p->Is(MCParticle::kZ) || p->Is(MCParticle::kW) || p->Is(MCParticle::kH) ||
143     p->Is(MCParticle::kZp) || p->Is(MCParticle::kZpp) ||
144     p->Is(MCParticle::kH0) || p->Is(MCParticle::kA0) || p->Is(MCParticle::kHp))) {
145 loizides 1.5 GenBosons->Add(p);
146 loizides 1.1 }
147 ceballos 1.10
148 ceballos 1.14 // photons
149 ceballos 1.10 else if (p->Status() == 1 && p->Is(MCParticle::kGamma) &&
150 ceballos 1.14 p->Pt() > fPtPhotonMin && p->AbsEta() < fEtaPhotonMax) {
151 ceballos 1.10 GenPhotons->Add(p);
152     }
153    
154 loizides 1.18 // W/Z -> lnu for Madgraph
155 loizides 1.28 if (p->IsParton() && p->NDaughters() >= 2) {
156 ceballos 1.16 CompositeParticle *diBoson = new CompositeParticle();
157     diBoson->AddDaughter(p->Daughter(0));
158     diBoson->AddDaughter(p->Daughter(1));
159 loizides 1.18 if (p->HasDaughter(MCParticle::kMu) && p->HasDaughter(MCParticle::kMuNu)) {
160 loizides 1.20 if (GetFillHist())
161 loizides 1.18 hDVMass[0]->Fill(TMath::Min(diBoson->Mass(),199.999));
162 loizides 1.28 const MCParticle *tmp_mu = p->FindDaughter(MCParticle::kMu);
163     while (tmp_mu->HasDaughter(MCParticle::kMu) &&
164     tmp_mu->FindDaughter(MCParticle::kMu)->IsGenerated())
165     tmp_mu = tmp_mu->FindDaughter(MCParticle::kMu);
166     GenLeptons->Add(tmp_mu);
167 ceballos 1.16 }
168 loizides 1.18 else if (p->HasDaughter(MCParticle::kEl) && p->HasDaughter(MCParticle::kElNu)) {
169 loizides 1.20 if (GetFillHist())
170 loizides 1.18 hDVMass[1]->Fill(TMath::Min(diBoson->Mass(),199.999));
171 loizides 1.28 const MCParticle *tmp_e = p->FindDaughter(MCParticle::kEl);
172     while (tmp_e->HasDaughter(MCParticle::kEl) &&
173     tmp_e->FindDaughter(MCParticle::kEl)->IsGenerated())
174     tmp_e = tmp_e->FindDaughter(MCParticle::kEl);
175     GenLeptons->Add(tmp_e);
176 ceballos 1.16 }
177 loizides 1.18 else if (p->HasDaughter(MCParticle::kTau) && p->HasDaughter(MCParticle::kTauNu)) {
178 loizides 1.20 if (GetFillHist())
179 loizides 1.18 hDVMass[2]->Fill(TMath::Min(diBoson->Mass(),199.999));
180     const MCParticle *tau = p->FindDaughter(MCParticle::kTau);
181     if (tau->HasDaughter(MCParticle::kMu))
182     GenLeptons->Add(tau->FindDaughter(MCParticle::kMu));
183     if (tau->HasDaughter(MCParticle::kEl))
184     GenLeptons->Add(tau->FindDaughter(MCParticle::kEl));
185 loizides 1.28 if (tau->HasDaughter(MCParticle::kTau)) {
186     const MCParticle *tau_second = tau->FindDaughter(MCParticle::kTau);
187     if (tau_second->HasDaughter(MCParticle::kMu))
188     GenLeptons->Add(tau_second->FindDaughter(MCParticle::kMu));
189     if (tau_second->HasDaughter(MCParticle::kEl))
190     GenLeptons->Add(tau_second->FindDaughter(MCParticle::kEl));
191     }
192 loizides 1.18 }
193     else if (p->Daughter(0)->Is(MCParticle::kMu) && p->Daughter(1)->Is(MCParticle::kMu)) {
194 loizides 1.20 if (GetFillHist())
195 loizides 1.18 hDVMass[3]->Fill(TMath::Min(diBoson->Mass(),199.999));
196 ceballos 1.16 GenLeptons->Add(p->Daughter(0));
197     GenLeptons->Add(p->Daughter(1));
198     }
199 loizides 1.18 else if (p->Daughter(0)->Is(MCParticle::kEl) && p->Daughter(1)->Is(MCParticle::kEl)) {
200 loizides 1.20 if (GetFillHist())
201 loizides 1.18 hDVMass[4]->Fill(TMath::Min(diBoson->Mass(),199.999));
202 ceballos 1.16 GenLeptons->Add(p->Daughter(0));
203     GenLeptons->Add(p->Daughter(1));
204     }
205 loizides 1.18 else if (p->Daughter(0)->Is(MCParticle::kTau) && p->Daughter(1)->Is(MCParticle::kTau)) {
206 loizides 1.20 if (GetFillHist())
207 loizides 1.18 hDVMass[5]->Fill(TMath::Min(diBoson->Mass(),199.999));
208     const MCParticle *tau0 = p->Daughter(0);
209     if (tau0->HasDaughter(MCParticle::kMu))
210     GenLeptons->Add(tau0->FindDaughter(MCParticle::kMu));
211     if (tau0->HasDaughter(MCParticle::kEl))
212     GenLeptons->Add(tau0->FindDaughter(MCParticle::kEl));
213     const MCParticle *tau1 = p->Daughter(1);
214     if (tau1->HasDaughter(MCParticle::kMu))
215     GenLeptons->Add(tau1->FindDaughter(MCParticle::kMu));
216     if (tau1->HasDaughter(MCParticle::kEl))
217     GenLeptons->Add(tau1->FindDaughter(MCParticle::kEl));
218 ceballos 1.16 }
219     delete diBoson;
220     }
221    
222 loizides 1.18 // t -> lnu for Madgraph
223     if (p->Is(MCParticle::kTop)) {
224 ceballos 1.16 CompositeParticle *diBoson = new CompositeParticle();
225 loizides 1.18 if (p->HasDaughter(MCParticle::kMu) && p->HasDaughter(MCParticle::kMuNu)) {
226 ceballos 1.16 diBoson->AddDaughter(p->FindDaughter(MCParticle::kMu));
227     diBoson->AddDaughter(p->FindDaughter(MCParticle::kMuNu));
228 loizides 1.20 if (GetFillHist())
229 loizides 1.18 hDVMass[6]->Fill(TMath::Min(diBoson->Mass(),199.999));
230 ceballos 1.16 GenLeptons->Add(p->FindDaughter(MCParticle::kMu));
231     }
232 loizides 1.18 else if (p->HasDaughter(MCParticle::kEl) && p->HasDaughter(MCParticle::kElNu)) {
233 ceballos 1.16 diBoson->AddDaughter(p->FindDaughter(MCParticle::kEl));
234     diBoson->AddDaughter(p->FindDaughter(MCParticle::kElNu));
235 loizides 1.20 if (GetFillHist())
236 loizides 1.18 hDVMass[7]->Fill(TMath::Min(diBoson->Mass(),199.999));
237 ceballos 1.16 GenLeptons->Add(p->FindDaughter(MCParticle::kEl));
238     }
239 loizides 1.18 else if (p->HasDaughter(MCParticle::kTau) && p->HasDaughter(MCParticle::kTauNu)) {
240 ceballos 1.16 diBoson->AddDaughter(p->FindDaughter(MCParticle::kTau));
241     diBoson->AddDaughter(p->FindDaughter(MCParticle::kTauNu));
242 loizides 1.20 if (GetFillHist())
243 loizides 1.18 hDVMass[8]->Fill(TMath::Min(diBoson->Mass(),199.999));
244     const MCParticle *tau = p->FindDaughter(MCParticle::kTau);
245     if (tau->HasDaughter(MCParticle::kMu))
246     GenLeptons->Add(tau->FindDaughter(MCParticle::kMu));
247     if (tau->HasDaughter(MCParticle::kEl))
248     GenLeptons->Add(tau->FindDaughter(MCParticle::kEl));
249     }
250     else if (!p->HasDaughter(MCParticle::kW)) {
251     for(UInt_t nd=0; nd<p->NDaughters(); ++nd)
252     if (p->Daughter(nd)->IsNot(MCParticle::kBottom) &&
253     p->Daughter(nd)->IsNot(MCParticle::kGamma))
254     diBoson->AddDaughter(p->Daughter(nd));
255 loizides 1.20 if (GetFillHist())
256 loizides 1.18 hDVMass[9]->Fill(TMath::Min(diBoson->Mass(),199.999));
257 ceballos 1.16 }
258     delete diBoson;
259     }
260 ceballos 1.26
261 loizides 1.27 // mass cut for given pid
262     if(fPdgIdCut && p->Is(fPdgIdCut) &&
263 ceballos 1.26 (p->Mass() < fMassMinCut || p->Mass() > fMassMaxCut)) {
264     SkipEvent();
265     return;
266     }
267    
268 loizides 1.27 } // end loop of particles
269 loizides 1.1
270 loizides 1.19 // sort according to pt
271     GenLeptons->Sort();
272     GenAllLeptons->Sort();
273     GenTaus->Sort();
274     GenNeutrinos->Sort();
275     GenQuarks->Sort();
276     GenqqHs->Sort();
277     GenBosons->Sort();
278     GenPhotons->Sort();
279    
280 loizides 1.8 // add objects to this event for other modules to use
281 loizides 1.13 AddObjThisEvt(GenLeptons);
282     AddObjThisEvt(GenAllLeptons);
283     AddObjThisEvt(GenTaus);
284     AddObjThisEvt(GenNeutrinos);
285     AddObjThisEvt(GenQuarks);
286     AddObjThisEvt(GenqqHs);
287     AddObjThisEvt(GenBosons);
288     AddObjThisEvt(GenPhotons);
289 ceballos 1.16
290 loizides 1.5 // fill histograms if requested
291 loizides 1.20 if (GetFillHist()) {
292 loizides 1.5
293 loizides 1.6 // leptons
294 loizides 1.1 hDGenLeptons[0]->Fill(GenLeptons->GetEntries());
295 loizides 1.5 for(UInt_t i=0; i<GenLeptons->GetEntries(); i++) {
296 loizides 1.1 hDGenLeptons[1]->Fill(GenLeptons->At(i)->Pt());
297 loizides 1.5 hDGenLeptons[2]->Fill(TMath::Abs(GenLeptons->At(i)->Eta()));
298     hDGenLeptons[3]->Fill(GenLeptons->At(i)->PhiDeg());
299     for(UInt_t j=i+1; j<GenLeptons->GetEntries(); j++) {
300 loizides 1.1 CompositeParticle *dilepton = new CompositeParticle();
301     dilepton->AddDaughter(GenLeptons->At(i));
302     dilepton->AddDaughter(GenLeptons->At(j));
303     hDGenLeptons[4]->Fill(dilepton->Mass());
304     delete dilepton;
305     }
306 ceballos 1.22 }
307     // looking at events with two leptons
308     if (GenLeptons->GetEntries() == 2) {
309     hDGenLeptons[5]->Fill(TMath::Min(TMath::Max(TMath::Abs(GenLeptons->At(0)->Eta()),
310     TMath::Abs(GenLeptons->At(1)->Eta())),
311 loizides 1.5 4.999));
312 ceballos 1.22 hDGenLeptons[6]->Fill(TMath::Min(TMath::Min(TMath::Abs(GenLeptons->At(0)->Eta()),
313     TMath::Abs(GenLeptons->At(1)->Eta())),
314 loizides 1.5 4.999));
315 ceballos 1.22 if (TMath::Abs(GenLeptons->At(0)->Eta()) < 2.5 &&
316     TMath::Abs(GenLeptons->At(1)->Eta()) < 2.5) {
317     hDGenLeptons[7]->Fill(TMath::Min(GenLeptons->At(0)->Pt(),199.999));
318     if (GenLeptons->At(0)->Pt() > 20.0) {
319     hDGenLeptons[8]->Fill(TMath::Min(GenLeptons->At(1)->Pt(),199.999));
320     if (GenLeptons->At(1)->Pt() > 10.0) {
321 loizides 1.1 CompositeParticle *dilepton = new CompositeParticle();
322 ceballos 1.22 dilepton->AddDaughter(GenLeptons->At(0));
323     dilepton->AddDaughter(GenLeptons->At(1));
324 loizides 1.1 hDGenLeptons[9]->Fill(TMath::Min(dilepton->Mass(),999.999));
325 ceballos 1.22 if(dilepton->Mass() > 12.0){
326     hDGenLeptons[10]->Fill(MathUtils::DeltaPhi(GenLeptons->At(0)->Phi(),
327     GenLeptons->At(1)->Phi())
328     * 180./ TMath::Pi());
329     hDGenLeptons[11]->Fill(MathUtils::DeltaR(GenLeptons->At(0)->Eta(), GenLeptons->At(0)->Phi(),
330     GenLeptons->At(1)->Eta(), GenLeptons->At(1)->Phi()));
331     }
332 loizides 1.1 delete dilepton;
333     }
334     }
335     }
336     }
337 ceballos 1.22 // looking at events with three leptons
338     if (GenLeptons->GetEntries() == 3) {
339     if (TMath::Abs(GenLeptons->At(0)->Eta()) < 2.5 &&
340     TMath::Abs(GenLeptons->At(1)->Eta()) < 2.5 &&
341     TMath::Abs(GenLeptons->At(2)->Eta()) < 2.5) {
342     hDGenLeptons[12]->Fill(TMath::Min(GenLeptons->At(0)->Pt(),199.999));
343     if (GenLeptons->At(0)->Pt() > 20.0) {
344     hDGenLeptons[13]->Fill(TMath::Min(GenLeptons->At(1)->Pt(),199.999));
345     hDGenLeptons[14]->Fill(TMath::Min(GenLeptons->At(2)->Pt(),199.999));
346     if (GenLeptons->At(1)->Pt() > 10.0 && GenLeptons->At(2)->Pt() > 10.0) {
347     CompositeParticle *dilepton01 = new CompositeParticle();
348     dilepton01->AddDaughter(GenLeptons->At(0));
349     dilepton01->AddDaughter(GenLeptons->At(1));
350     CompositeParticle *dilepton02 = new CompositeParticle();
351     dilepton02->AddDaughter(GenLeptons->At(0));
352     dilepton02->AddDaughter(GenLeptons->At(2));
353     CompositeParticle *dilepton12 = new CompositeParticle();
354     dilepton12->AddDaughter(GenLeptons->At(1));
355     dilepton12->AddDaughter(GenLeptons->At(2));
356     hDGenLeptons[15]->Fill(TMath::Min(dilepton01->Mass(),999.999));
357     hDGenLeptons[15]->Fill(TMath::Min(dilepton02->Mass(),999.999));
358     hDGenLeptons[15]->Fill(TMath::Min(dilepton12->Mass(),999.999));
359     CompositeParticle *trilepton = new CompositeParticle();
360     trilepton->AddDaughter(GenLeptons->At(0));
361     trilepton->AddDaughter(GenLeptons->At(1));
362     trilepton->AddDaughter(GenLeptons->At(2));
363     hDGenLeptons[16]->Fill(TMath::Min(trilepton->Mass(),999.999));
364 loizides 1.27 Double_t deltaR[3] = {MathUtils::DeltaR(GenLeptons->At(0)->Eta(),
365     GenLeptons->At(0)->Phi(),
366     GenLeptons->At(1)->Eta(),
367     GenLeptons->At(1)->Phi()),
368     MathUtils::DeltaR(GenLeptons->At(0)->Eta(),
369     GenLeptons->At(0)->Phi(),
370     GenLeptons->At(2)->Eta(),
371     GenLeptons->At(2)->Phi()),
372     MathUtils::DeltaR(GenLeptons->At(1)->Eta(),
373     GenLeptons->At(1)->Phi(),
374     GenLeptons->At(2)->Eta(),
375     GenLeptons->At(2)->Phi())};
376     Double_t deltaRMin = deltaR[0];
377     for(Int_t i=1; i<3; i++) if(deltaRMin > deltaR[i]) deltaRMin = deltaR[i];
378 ceballos 1.22 hDGenLeptons[17]->Fill(deltaRMin);
379    
380     delete dilepton01;
381     delete dilepton02;
382     delete dilepton12;
383     delete trilepton;
384     }
385     }
386     }
387     }
388     // looking at events with four leptons
389     if (GenLeptons->GetEntries() == 4) {
390     if (TMath::Abs(GenLeptons->At(0)->Eta()) < 2.5 &&
391     TMath::Abs(GenLeptons->At(1)->Eta()) < 2.5 &&
392     TMath::Abs(GenLeptons->At(2)->Eta()) < 2.5 &&
393     TMath::Abs(GenLeptons->At(3)->Eta()) < 2.5) {
394     hDGenLeptons[18]->Fill(TMath::Min(GenLeptons->At(0)->Pt(),199.999));
395     if (GenLeptons->At(0)->Pt() > 20.0) {
396     hDGenLeptons[19]->Fill(TMath::Min(GenLeptons->At(1)->Pt(),199.999));
397     hDGenLeptons[20]->Fill(TMath::Min(GenLeptons->At(2)->Pt(),199.999));
398     hDGenLeptons[21]->Fill(TMath::Min(GenLeptons->At(3)->Pt(),199.999));
399     if (GenLeptons->At(1)->Pt() > 10.0 && GenLeptons->At(2)->Pt() > 10.0 &&
400     GenLeptons->At(3)->Pt() > 10.0) {
401     CompositeParticle *dilepton01 = new CompositeParticle();
402     dilepton01->AddDaughter(GenLeptons->At(0));
403     dilepton01->AddDaughter(GenLeptons->At(1));
404     CompositeParticle *dilepton02 = new CompositeParticle();
405     dilepton02->AddDaughter(GenLeptons->At(0));
406     dilepton02->AddDaughter(GenLeptons->At(2));
407     CompositeParticle *dilepton03 = new CompositeParticle();
408     dilepton03->AddDaughter(GenLeptons->At(0));
409     dilepton03->AddDaughter(GenLeptons->At(3));
410     CompositeParticle *dilepton12 = new CompositeParticle();
411     dilepton12->AddDaughter(GenLeptons->At(1));
412     dilepton12->AddDaughter(GenLeptons->At(2));
413     CompositeParticle *dilepton13 = new CompositeParticle();
414     dilepton13->AddDaughter(GenLeptons->At(1));
415     dilepton13->AddDaughter(GenLeptons->At(3));
416     CompositeParticle *dilepton23 = new CompositeParticle();
417     dilepton23->AddDaughter(GenLeptons->At(2));
418     dilepton23->AddDaughter(GenLeptons->At(3));
419     hDGenLeptons[22]->Fill(TMath::Min(dilepton01->Mass(),999.999));
420     hDGenLeptons[22]->Fill(TMath::Min(dilepton02->Mass(),999.999));
421     hDGenLeptons[22]->Fill(TMath::Min(dilepton03->Mass(),999.999));
422     hDGenLeptons[22]->Fill(TMath::Min(dilepton12->Mass(),999.999));
423     hDGenLeptons[22]->Fill(TMath::Min(dilepton13->Mass(),999.999));
424     hDGenLeptons[22]->Fill(TMath::Min(dilepton23->Mass(),999.999));
425     CompositeParticle *fourlepton = new CompositeParticle();
426     fourlepton->AddDaughter(GenLeptons->At(0));
427     fourlepton->AddDaughter(GenLeptons->At(1));
428     fourlepton->AddDaughter(GenLeptons->At(2));
429     fourlepton->AddDaughter(GenLeptons->At(3));
430     hDGenLeptons[23]->Fill(TMath::Min(fourlepton->Mass(),999.999));
431 loizides 1.27 Double_t deltaR[6] = {MathUtils::DeltaR(GenLeptons->At(0)->Eta(),
432     GenLeptons->At(0)->Phi(),
433     GenLeptons->At(1)->Eta(),
434     GenLeptons->At(1)->Phi()),
435     MathUtils::DeltaR(GenLeptons->At(0)->Eta(),
436     GenLeptons->At(0)->Phi(),
437     GenLeptons->At(2)->Eta(),
438     GenLeptons->At(2)->Phi()),
439     MathUtils::DeltaR(GenLeptons->At(0)->Eta(),
440     GenLeptons->At(0)->Phi(),
441     GenLeptons->At(3)->Eta(),
442     GenLeptons->At(3)->Phi()),
443     MathUtils::DeltaR(GenLeptons->At(1)->Eta(),
444     GenLeptons->At(1)->Phi(),
445     GenLeptons->At(2)->Eta(),
446     GenLeptons->At(2)->Phi()),
447     MathUtils::DeltaR(GenLeptons->At(1)->Eta(),
448     GenLeptons->At(1)->Phi(),
449     GenLeptons->At(3)->Eta(),
450     GenLeptons->At(3)->Phi()),
451     MathUtils::DeltaR(GenLeptons->At(2)->Eta(),
452     GenLeptons->At(2)->Phi(),
453     GenLeptons->At(3)->Eta(),
454     GenLeptons->At(3)->Phi())};
455     Double_t deltaRMin = deltaR[0];
456     for(Int_t i=1; i<6; i++) if(deltaRMin > deltaR[i]) deltaRMin = deltaR[i];
457 ceballos 1.22 hDGenLeptons[24]->Fill(deltaRMin);
458    
459     delete dilepton01;
460     delete dilepton02;
461     delete dilepton03;
462     delete dilepton12;
463     delete dilepton13;
464     delete dilepton23;
465     delete fourlepton;
466     }
467     }
468     }
469     }
470 loizides 1.1
471 loizides 1.6 // all leptons
472 ceballos 1.3 hDGenAllLeptons[0]->Fill(GenAllLeptons->GetEntries());
473 loizides 1.5 for(UInt_t i=0; i<GenAllLeptons->GetEntries(); i++) {
474 ceballos 1.3 hDGenAllLeptons[1]->Fill(GenAllLeptons->At(i)->Pt());
475     hDGenAllLeptons[2]->Fill(GenAllLeptons->At(i)->Eta());
476 loizides 1.5 hDGenAllLeptons[3]->Fill(GenAllLeptons->At(i)->PhiDeg());
477 ceballos 1.3 }
478    
479 loizides 1.6 // taus
480 loizides 1.1 hDGenTaus[0]->Fill(GenTaus->GetEntries());
481 loizides 1.5 for(UInt_t i=0; i<GenTaus->GetEntries(); i++) {
482 loizides 1.1 hDGenTaus[1]->Fill(GenTaus->At(i)->Pt());
483     hDGenTaus[2]->Fill(GenTaus->At(i)->Eta());
484 loizides 1.5 hDGenTaus[3]->Fill(GenTaus->At(i)->PhiDeg());
485 loizides 1.1 }
486    
487 loizides 1.6 // neutrinos
488 loizides 1.1 hDGenNeutrinos[0]->Fill(GenNeutrinos->GetEntries());
489     CompositeParticle *neutrinoTotal = new CompositeParticle();
490 loizides 1.5 for(UInt_t i=0; i<GenNeutrinos->GetEntries(); i++) {
491     if (GenNeutrinos->At(i)->HasMother())
492 loizides 1.1 neutrinoTotal->AddDaughter(GenNeutrinos->At(i));
493     }
494 loizides 1.5 if (GenNeutrinos->GetEntries() > 0) {
495 loizides 1.1 hDGenNeutrinos[1]->Fill(neutrinoTotal->Pt());
496     hDGenNeutrinos[2]->Fill(neutrinoTotal->Eta());
497 loizides 1.5 hDGenNeutrinos[3]->Fill(neutrinoTotal->PhiDeg());
498 loizides 1.1 }
499     delete neutrinoTotal;
500    
501 loizides 1.6 // quarks
502 loizides 1.1 hDGenQuarks[0]->Fill(GenQuarks->GetEntries());
503 loizides 1.5 for(UInt_t i=0; i<GenQuarks->GetEntries(); i++) {
504     for(UInt_t j=i+1; j<GenQuarks->GetEntries(); j++) {
505 loizides 1.1 CompositeParticle *dijet = new CompositeParticle();
506     dijet->AddDaughter(GenQuarks->At(i));
507     dijet->AddDaughter(GenQuarks->At(j));
508     hDGenQuarks[1]->Fill(dijet->Pt());
509     hDGenQuarks[2]->Fill(dijet->Mass());
510 loizides 1.5 if (TMath::Abs(GenQuarks->At(i)->Eta()) < 2.5 &&
511     TMath::Abs(GenQuarks->At(j)->Eta()) < 2.5) {
512 loizides 1.1 hDGenQuarks[3]->Fill(dijet->Pt());
513     hDGenQuarks[4]->Fill(dijet->Mass());
514     }
515     delete dijet;
516     }
517 ceballos 1.2 // b quark info
518 loizides 1.5 if (GenQuarks->At(i)->AbsPdgId() == 5) {
519 ceballos 1.2 hDGenQuarks[5]->Fill(GenQuarks->At(i)->Pt());
520     hDGenQuarks[6]->Fill(GenQuarks->At(i)->Eta());
521     hDGenQuarks[7]->Fill(GenQuarks->At(i)->Phi());
522 ceballos 1.22 if (GenLeptons->GetEntries() >= 2 &&
523     GenLeptons->At(0)->Pt() > 20 &&
524     GenLeptons->At(1)->Pt() > 15) {
525     if (TMath::Abs(GenLeptons->At(0)->Eta()) < 2.5 &&
526     TMath::Abs(GenLeptons->At(1)->Eta()) < 2.5) {
527 ceballos 1.2 hDGenQuarks[8]->Fill(GenQuarks->At(i)->Pt());
528     hDGenQuarks[9]->Fill(GenQuarks->At(i)->Eta());
529     hDGenQuarks[10]->Fill(GenQuarks->At(i)->Phi());
530     }
531     }
532     }
533     // t quark info
534 loizides 1.5 else if (GenQuarks->At(i)->AbsPdgId() == 6) {
535 ceballos 1.2 hDGenQuarks[11]->Fill(GenQuarks->At(i)->Pt());
536     hDGenQuarks[12]->Fill(GenQuarks->At(i)->Eta());
537     hDGenQuarks[13]->Fill(GenQuarks->At(i)->Phi());
538     }
539     // light quark info
540     else {
541     hDGenQuarks[14]->Fill(GenQuarks->At(i)->Pt());
542     hDGenQuarks[15]->Fill(GenQuarks->At(i)->Eta());
543     hDGenQuarks[16]->Fill(GenQuarks->At(i)->Phi());
544     }
545 loizides 1.1 }
546    
547 loizides 1.6 // wbf
548 loizides 1.5 if (GenqqHs->GetEntries() == 2) {
549 loizides 1.1 hDGenWBF[0]->Fill(MathUtils::DeltaPhi(GenqqHs->At(0)->Phi(),
550     GenqqHs->At(1)->Phi()) * 180./ TMath::Pi());
551 loizides 1.5 hDGenWBF[1]->Fill(TMath::Abs(GenqqHs->At(0)->Eta()-GenqqHs->At(1)->Eta()));
552 loizides 1.1 hDGenWBF[2]->Fill(TMath::Max(GenqqHs->At(0)->Pt(),GenqqHs->At(1)->Pt()));
553     hDGenWBF[3]->Fill(TMath::Min(GenqqHs->At(0)->Pt(),GenqqHs->At(1)->Pt()));
554     CompositeParticle *diqq = new CompositeParticle();
555     diqq->AddDaughter(GenqqHs->At(0));
556     diqq->AddDaughter(GenqqHs->At(1));
557     hDGenWBF[4]->Fill(diqq->Mass());
558     delete diqq;
559     }
560    
561 loizides 1.6 // bosons
562 loizides 1.1 hDGenBosons[0]->Fill(GenBosons->GetEntries());
563 loizides 1.5 for(UInt_t i=0; i<GenBosons->GetEntries(); i++) {
564 loizides 1.1 hDGenBosons[1]->Fill(GenBosons->At(i)->Pt());
565     hDGenBosons[2]->Fill(GenBosons->At(i)->Eta());
566     hDGenBosons[3]->Fill(GenBosons->At(i)->Mass());
567     }
568 ceballos 1.10
569     // photons
570     hDGenPhotons[0]->Fill(GenPhotons->GetEntries());
571     for(UInt_t i=0; i<GenPhotons->GetEntries(); i++) {
572     hDGenPhotons[1]->Fill(GenPhotons->At(i)->Pt());
573     hDGenPhotons[2]->Fill(GenPhotons->At(i)->Eta());
574     }
575 loizides 1.1 }
576     }
577    
578     //--------------------------------------------------------------------------------------------------
579     void GeneratorMod::SlaveBegin()
580     {
581 loizides 1.5 // Book branch and histograms if wanted.
582    
583 loizides 1.1 ReqBranch(fMCPartName, fParticles);
584    
585 loizides 1.5 // fill histograms
586 loizides 1.20 if (GetFillHist()) {
587 loizides 1.5 char sb[1024];
588 loizides 1.6 // leptons from W
589 loizides 1.1 sprintf(sb,"hDGenLeptons_%d", 0); hDGenLeptons[0] = new TH1D(sb,sb,10,-0.5,9.5);
590     sprintf(sb,"hDGenLeptons_%d", 1); hDGenLeptons[1] = new TH1D(sb,sb,100,0.0,200.0);
591     sprintf(sb,"hDGenLeptons_%d", 2); hDGenLeptons[2] = new TH1D(sb,sb,50,0.0,5.0);
592     sprintf(sb,"hDGenLeptons_%d", 3); hDGenLeptons[3] = new TH1D(sb,sb,90,0.0,180.0);
593     sprintf(sb,"hDGenLeptons_%d", 4); hDGenLeptons[4] = new TH1D(sb,sb,1000,0.0,1000.0);
594     sprintf(sb,"hDGenLeptons_%d", 5); hDGenLeptons[5] = new TH1D(sb,sb,50,0.0,5.0);
595     sprintf(sb,"hDGenLeptons_%d", 6); hDGenLeptons[6] = new TH1D(sb,sb,50,0.0,5.0);
596     sprintf(sb,"hDGenLeptons_%d", 7); hDGenLeptons[7] = new TH1D(sb,sb,100,0.0,200.0);
597     sprintf(sb,"hDGenLeptons_%d", 8); hDGenLeptons[8] = new TH1D(sb,sb,100,0.0,200.0);
598     sprintf(sb,"hDGenLeptons_%d", 9); hDGenLeptons[9] = new TH1D(sb,sb,1000,0.0,1000.0);
599     sprintf(sb,"hDGenLeptons_%d",10); hDGenLeptons[10] = new TH1D(sb,sb,90,0.0,180.0);
600 ceballos 1.22 sprintf(sb,"hDGenLeptons_%d",11); hDGenLeptons[11] = new TH1D(sb,sb,100,0.0,5.0);
601     sprintf(sb,"hDGenLeptons_%d",12); hDGenLeptons[12] = new TH1D(sb,sb,100,0.0,200.0);
602     sprintf(sb,"hDGenLeptons_%d",13); hDGenLeptons[13] = new TH1D(sb,sb,100,0.0,200.0);
603     sprintf(sb,"hDGenLeptons_%d",14); hDGenLeptons[14] = new TH1D(sb,sb,100,0.0,200.0);
604     sprintf(sb,"hDGenLeptons_%d",15); hDGenLeptons[15] = new TH1D(sb,sb,1000,0.0,1000.0);
605     sprintf(sb,"hDGenLeptons_%d",16); hDGenLeptons[16] = new TH1D(sb,sb,1000,0.0,1000.0);
606     sprintf(sb,"hDGenLeptons_%d",17); hDGenLeptons[17] = new TH1D(sb,sb,100,0.0,5.0);
607     sprintf(sb,"hDGenLeptons_%d",18); hDGenLeptons[18] = new TH1D(sb,sb,100,0.0,200.0);
608     sprintf(sb,"hDGenLeptons_%d",19); hDGenLeptons[19] = new TH1D(sb,sb,100,0.0,200.0);
609     sprintf(sb,"hDGenLeptons_%d",20); hDGenLeptons[20] = new TH1D(sb,sb,100,0.0,200.0);
610     sprintf(sb,"hDGenLeptons_%d",21); hDGenLeptons[21] = new TH1D(sb,sb,100,0.0,200.0);
611     sprintf(sb,"hDGenLeptons_%d",22); hDGenLeptons[22] = new TH1D(sb,sb,1000,0.0,1000.0);
612     sprintf(sb,"hDGenLeptons_%d",23); hDGenLeptons[23] = new TH1D(sb,sb,1000,0.0,1000.0);
613     sprintf(sb,"hDGenLeptons_%d",24); hDGenLeptons[24] = new TH1D(sb,sb,100,0.0,5.0);
614     for(Int_t i=0; i<25; i++) AddOutput(hDGenLeptons[i]);
615 loizides 1.1
616 loizides 1.6 // all leptons
617 ceballos 1.3 sprintf(sb,"hDGenAllLeptons_%d", 0); hDGenAllLeptons[0] = new TH1D(sb,sb,10,-0.5,9.5);
618     sprintf(sb,"hDGenAllLeptons_%d", 1); hDGenAllLeptons[1] = new TH1D(sb,sb,100,0.0,200.0);
619 ceballos 1.14 sprintf(sb,"hDGenAllLeptons_%d", 2); hDGenAllLeptons[2] = new TH1D(sb,sb,100,-5.0,5.0);
620 ceballos 1.3 sprintf(sb,"hDGenAllLeptons_%d", 3); hDGenAllLeptons[3] = new TH1D(sb,sb,90,0.0,180.0);
621 loizides 1.5 for(Int_t i=0; i<4; i++) AddOutput(hDGenAllLeptons[i]);
622 ceballos 1.3
623 loizides 1.6 // taus
624 loizides 1.1 sprintf(sb,"hDGenTaus_%d", 0); hDGenTaus[0] = new TH1D(sb,sb,10,-0.5,9.5);
625     sprintf(sb,"hDGenTaus_%d", 1); hDGenTaus[1] = new TH1D(sb,sb,100,0.0,200.0);
626     sprintf(sb,"hDGenTaus_%d", 2); hDGenTaus[2] = new TH1D(sb,sb,100,-5.0,5.0);
627     sprintf(sb,"hDGenTaus_%d", 3); hDGenTaus[3] = new TH1D(sb,sb,90,0.0,180.0);
628 loizides 1.5 for(Int_t i=0; i<4; i++) AddOutput(hDGenTaus[i]);
629 loizides 1.1
630 loizides 1.6 // neutrinos
631 loizides 1.1 sprintf(sb,"hDGenNeutrinos_%d", 0); hDGenNeutrinos[0] = new TH1D(sb,sb,10,-0.5,9.5);
632     sprintf(sb,"hDGenNeutrinos_%d", 1); hDGenNeutrinos[1] = new TH1D(sb,sb,100,0.0,200.0);
633     sprintf(sb,"hDGenNeutrinos_%d", 2); hDGenNeutrinos[2] = new TH1D(sb,sb,100,-5.0,5.0);
634     sprintf(sb,"hDGenNeutrinos_%d", 3); hDGenNeutrinos[3] = new TH1D(sb,sb,90,0.0,180.0);
635 loizides 1.5 for(Int_t i=0; i<4; i++) AddOutput(hDGenNeutrinos[i]);
636 loizides 1.1
637 loizides 1.6 // quarks
638 loizides 1.1 sprintf(sb,"hDGenQuarks_%d", 0); hDGenQuarks[0] = new TH1D(sb,sb,10,-0.5,9.5);
639     sprintf(sb,"hDGenQuarks_%d", 1); hDGenQuarks[1] = new TH1D(sb,sb,200,0.0,400.);
640     sprintf(sb,"hDGenQuarks_%d", 2); hDGenQuarks[2] = new TH1D(sb,sb,2000,0.0,2000.);
641     sprintf(sb,"hDGenQuarks_%d", 3); hDGenQuarks[3] = new TH1D(sb,sb,200,0.0,400.);
642     sprintf(sb,"hDGenQuarks_%d", 4); hDGenQuarks[4] = new TH1D(sb,sb,2000,0.0,2000.);
643 ceballos 1.2 sprintf(sb,"hDGenQuarks_%d", 5); hDGenQuarks[5] = new TH1D(sb,sb,200,0.0,400.);
644     sprintf(sb,"hDGenQuarks_%d", 6); hDGenQuarks[6] = new TH1D(sb,sb,200,-10.0,10.);
645 loizides 1.5 sprintf(sb,"hDGenQuarks_%d", 7); hDGenQuarks[7] = new TH1D(sb,sb,200,-TMath::Pi(),
646     TMath::Pi());
647 ceballos 1.2 sprintf(sb,"hDGenQuarks_%d", 8); hDGenQuarks[8] = new TH1D(sb,sb,200,0.0,400.);
648     sprintf(sb,"hDGenQuarks_%d", 9); hDGenQuarks[9] = new TH1D(sb,sb,200,-10.0,10.);
649 loizides 1.5 sprintf(sb,"hDGenQuarks_%d",10); hDGenQuarks[10] = new TH1D(sb,sb,200,-TMath::Pi(),
650     TMath::Pi());
651 ceballos 1.2 sprintf(sb,"hDGenQuarks_%d",11); hDGenQuarks[11] = new TH1D(sb,sb,200,0.0,400.);
652     sprintf(sb,"hDGenQuarks_%d",12); hDGenQuarks[12] = new TH1D(sb,sb,200,-10.0,10.);
653 loizides 1.5 sprintf(sb,"hDGenQuarks_%d",13); hDGenQuarks[13] = new TH1D(sb,sb,200,-TMath::Pi(),
654     TMath::Pi());
655 ceballos 1.2 sprintf(sb,"hDGenQuarks_%d",14); hDGenQuarks[14] = new TH1D(sb,sb,200,0.0,400.);
656     sprintf(sb,"hDGenQuarks_%d",15); hDGenQuarks[15] = new TH1D(sb,sb,200,-10.0,10.);
657 loizides 1.5 sprintf(sb,"hDGenQuarks_%d",16); hDGenQuarks[16] = new TH1D(sb,sb,200,-TMath::Pi(),
658     TMath::Pi());
659     for(Int_t i=0; i<17; i++) AddOutput(hDGenQuarks[i]);
660 loizides 1.1
661     // qqH
662     sprintf(sb,"hDGenWBF_%d", 0); hDGenWBF[0] = new TH1D(sb,sb,90,0.0,180.);
663     sprintf(sb,"hDGenWBF_%d", 1); hDGenWBF[1] = new TH1D(sb,sb,100,0.0,10.);
664     sprintf(sb,"hDGenWBF_%d", 2); hDGenWBF[2] = new TH1D(sb,sb,200,0.0,400.);
665     sprintf(sb,"hDGenWBF_%d", 3); hDGenWBF[3] = new TH1D(sb,sb,200,0.0,400.);
666     sprintf(sb,"hDGenWBF_%d", 4); hDGenWBF[4] = new TH1D(sb,sb,200,0.0,4000.);
667 loizides 1.5 for(Int_t i=0; i<5; i++) AddOutput(hDGenWBF[i]);
668 loizides 1.1
669 loizides 1.6 // bosons
670 loizides 1.1 sprintf(sb,"hDGenBosons_%d", 0); hDGenBosons[0] = new TH1D(sb,sb,10,-0.5,9.5);
671     sprintf(sb,"hDGenBosons_%d", 1); hDGenBosons[1] = new TH1D(sb,sb,200,0.0,400.0);
672     sprintf(sb,"hDGenBosons_%d", 2); hDGenBosons[2] = new TH1D(sb,sb,100,-5.0,5.0);
673     sprintf(sb,"hDGenBosons_%d", 3); hDGenBosons[3] = new TH1D(sb,sb,2000,0.0,2000.0);
674 loizides 1.5 for(Int_t i=0; i<4; i++) AddOutput(hDGenBosons[i]);
675 ceballos 1.10
676     // photons
677     sprintf(sb,"hDGenPhotons_%d", 0); hDGenPhotons[0] = new TH1D(sb,sb,10,-0.5,9.5);
678     sprintf(sb,"hDGenPhotons_%d", 1); hDGenPhotons[1] = new TH1D(sb,sb,200,0.0,400.0);
679 ceballos 1.14 sprintf(sb,"hDGenPhotons_%d", 2); hDGenPhotons[2] = new TH1D(sb,sb,100,-5.0,5.0);
680 ceballos 1.10 for(Int_t i=0; i<3; i++) AddOutput(hDGenPhotons[i]);
681 ceballos 1.16
682     // auxiliar
683     sprintf(sb,"hDVMass_%d", 0); hDVMass[0] = new TH1D(sb,sb,200,0.,200.);
684     sprintf(sb,"hDVMass_%d", 1); hDVMass[1] = new TH1D(sb,sb,200,0.,200.);
685     sprintf(sb,"hDVMass_%d", 2); hDVMass[2] = new TH1D(sb,sb,200,0.,200.);
686     sprintf(sb,"hDVMass_%d", 3); hDVMass[3] = new TH1D(sb,sb,200,0.,200.);
687     sprintf(sb,"hDVMass_%d", 4); hDVMass[4] = new TH1D(sb,sb,200,0.,200.);
688     sprintf(sb,"hDVMass_%d", 5); hDVMass[5] = new TH1D(sb,sb,200,0.,200.);
689     sprintf(sb,"hDVMass_%d", 6); hDVMass[6] = new TH1D(sb,sb,200,0.,200.);
690     sprintf(sb,"hDVMass_%d", 7); hDVMass[7] = new TH1D(sb,sb,200,0.,200.);
691     sprintf(sb,"hDVMass_%d", 8); hDVMass[8] = new TH1D(sb,sb,200,0.,200.);
692     sprintf(sb,"hDVMass_%d", 9); hDVMass[9] = new TH1D(sb,sb,200,0.,200.);
693     for(Int_t i=0; i<10; i++) AddOutput(hDVMass[i]);
694 loizides 1.1 }
695     }