211 |
|
|
212 |
|
if (GetEventHeader()->IsMC()) { |
213 |
|
LoadBranch(fPileupInfoName); |
214 |
< |
Int_t npu[4] = {0,0,0,0}; |
214 |
> |
Double_t npu[4] = {0.,0.,0.,0.}; |
215 |
|
for (UInt_t i=0; i<fPileupInfo->GetEntries(); ++i) { |
216 |
|
const PileupInfo *puinfo = fPileupInfo->At(i); |
217 |
|
if (puinfo->GetBunchCrossing()==0) npu[0]= puinfo->GetPU_NumInteractions(); |
219 |
|
else if (puinfo->GetBunchCrossing()==1) npu[2] = puinfo->GetPU_NumInteractions(); |
220 |
|
if (puinfo->GetBunchCrossing()==0) npu[3]= puinfo->GetPU_NumMean(); |
221 |
|
} |
222 |
– |
|
222 |
|
hNPU->Fill(npu[0]); |
223 |
|
hNPU50ns->Fill(npu[0],npu[1],npu[2]); |
224 |
|
hNPUTrue->Fill(npu[3]); |
260 |
|
hNPU50ns = new TH3D("hNPU50ns", "hNPU50ns", 201, -0.5, 200.5, 201, -0.5, 200.5, 201, -0.5, 200.5); |
261 |
|
AddOutput(hNPU50ns); |
262 |
|
|
263 |
< |
hNPUTrue = new TH1D("hNPUTrue", "hNPUTrue", 201, -0.5, 200.5); |
263 |
> |
hNPUTrue = new TH1D("hNPUTrue", "hNPUTrue", 2000, 0.0, 200.0); |
264 |
|
AddOutput(hNPUTrue); |
265 |
|
|
266 |
|
} |