@@ -333,33 +333,30 @@ void GFWWeights::createNUA(bool IntegrateOverCentAndPt)
333333 TH1D * h1;
334334 if (!fW_data || fW_data ->GetEntries () < 1 )
335335 return ;
336- if (IntegrateOverCentAndPt) {
337- if (fAccInt )
338- delete fAccInt ;
339- fAccInt = reinterpret_cast <TH3D *>(fW_data ->At (0 )->Clone (" IntegratedAcceptance" ));
340- fAccInt ->SetDirectory (0 );
341- fAccInt ->Sumw2 ();
342- for (int etai = 1 ; etai <= fAccInt ->GetNbinsY (); etai++) {
343- fAccInt ->GetYaxis ()->SetRange (etai, etai);
336+ if (fAccInt )
337+ delete fAccInt ;
338+ fAccInt = reinterpret_cast <TH3D *>(fW_data ->At (0 )->Clone (" IntegratedAcceptance" ));
339+ fAccInt ->SetDirectory (0 );
340+ fAccInt ->Sumw2 ();
341+ for (int etai = 1 ; etai <= fAccInt ->GetNbinsY (); etai++) {
342+ fAccInt ->GetYaxis ()->SetRange (etai, etai);
343+ if (fAccInt ->Integral () < 1 )
344+ continue ;
345+ for (int vzi = 1 ; vzi <= fAccInt ->GetNbinsZ (); vzi++) {
346+ fAccInt ->GetZaxis ()->SetRange (vzi, vzi);
344347 if (fAccInt ->Integral () < 1 )
345348 continue ;
346- for (int vzi = 1 ; vzi <= fAccInt ->GetNbinsZ (); vzi++) {
347- fAccInt ->GetZaxis ()->SetRange (vzi, vzi);
348- if (fAccInt ->Integral () < 1 )
349- continue ;
350- h1 = reinterpret_cast <TH1D *>(fAccInt ->Project3D (" x" ));
351- double maxv = h1->GetMaximum ();
352- for (int phii = 1 ; phii <= h1->GetNbinsX (); phii++) {
353- fAccInt ->SetBinContent (phii, etai, vzi, fAccInt ->GetBinContent (phii, etai, vzi) / maxv);
354- fAccInt ->SetBinError (phii, etai, vzi, fAccInt ->GetBinError (phii, etai, vzi) / maxv);
355- }
356- delete h1;
349+ h1 = reinterpret_cast <TH1D *>(fAccInt ->Project3D (" x" ));
350+ double maxv = h1->GetMaximum ();
351+ for (int phii = 1 ; phii <= h1->GetNbinsX (); phii++) {
352+ fAccInt ->SetBinContent (phii, etai, vzi, fAccInt ->GetBinContent (phii, etai, vzi) / maxv);
353+ fAccInt ->SetBinError (phii, etai, vzi, fAccInt ->GetBinError (phii, etai, vzi) / maxv);
357354 }
358- fAccInt -> GetZaxis ()-> SetRange ( 1 , fAccInt -> GetNbinsZ ()) ;
355+ delete h1 ;
359356 }
360- fAccInt ->GetYaxis ()->SetRange (1 , fAccInt ->GetNbinsY ());
361- return ;
357+ fAccInt ->GetZaxis ()->SetRange (1 , fAccInt ->GetNbinsZ ());
362358 }
359+ fAccInt ->GetYaxis ()->SetRange (1 , fAccInt ->GetNbinsY ());
363360};
364361TH1D * GFWWeights::getdNdPhi ()
365362{
@@ -393,21 +390,18 @@ void GFWWeights::createNUE(bool IntegrateOverCentrality)
393390 TH3D * den = 0 ;
394391 if (!fW_mcrec || !fW_mcgen || fW_mcrec ->GetEntries () < 1 || fW_mcgen ->GetEntries () < 1 )
395392 return ;
396- if (IntegrateOverCentrality) {
397- num = reinterpret_cast <TH3D *>(fW_mcrec ->At (0 ));
398- den = reinterpret_cast <TH3D *>(fW_mcgen ->At (0 ));
399- num->Sumw2 ();
400- den->Sumw2 ();
401- num->RebinY (2 );
402- den->RebinY (2 );
403- num->RebinZ (5 );
404- den->RebinZ (5 );
405- delete fEffInt ;
406- fEffInt = reinterpret_cast <TH3D *>(num->Clone (" Efficiency_Integrated" ));
407- fEffInt ->SetDirectory (0 );
408- fEffInt ->Divide (den);
409- return ;
410- }
393+ num = reinterpret_cast <TH3D *>(fW_mcrec ->At (0 ));
394+ den = reinterpret_cast <TH3D *>(fW_mcgen ->At (0 ));
395+ num->Sumw2 ();
396+ den->Sumw2 ();
397+ num->RebinY (2 );
398+ den->RebinY (2 );
399+ num->RebinZ (5 );
400+ den->RebinZ (5 );
401+ delete fEffInt ;
402+ fEffInt = reinterpret_cast <TH3D *>(num->Clone (" Efficiency_Integrated" ));
403+ fEffInt ->SetDirectory (0 );
404+ fEffInt ->Divide (den);
411405};
412406void GFWWeights::readAndMerge (const TString& filelinks, const TString& listName, bool addData, bool addRec, bool addGen)
413407{
0 commit comments