// Calculate final results for integrated flow of RPs and POIs.
// (to be improved - this method can be implemented much better)
- Int_t rp = -1;
+ Int_t rp = 0;
if(rpPoi == "RP")
{
// Fill common result histograms for differential flow.
// (to be improved - this method can be implemented much better)
- Int_t rp = -1;
+ Int_t rp = 0;
if(rpPoi == "RP")
{
{
// Calculate differential flow for RPs/POIs vs pt/eta from cumulants.
- Int_t rp = -1; // RP or POI
- Int_t pe = -1; // pt or eta
+ Int_t rp = 0; // RP or POI
+ Int_t pe = 0; // pt or eta
if(rpPoi == "RP")
{
{
// Calculate errors of differential flow.
- Int_t rp = -1; // RP or POI
- Int_t pe = -1; // pt or eta
+ Int_t rp = 0; // RP or POI
+ Int_t pe = 0; // pt or eta
if(rpPoi == "RP")
{
{
// Calculate cumulants for differential flow.
- Int_t rp = -1; // RP or POI
- Int_t pe = -1; // pt or eta
+ Int_t rp = 0; // RP or POI
+ Int_t pe = 0; // pt or eta
if(rpPoi == "RP")
{
if(fWeightsList->FindObject("phi_weights"))
{
fPhiWeights = dynamic_cast<TH1F*>(fWeightsList->FindObject("phi_weights"));
+ if(!fPhiWeights){printf("\n WARNING (GFC): !fPhiWeights !!!!\n");exit(0);}
if(TMath::Abs(fPhiWeights->GetBinWidth(1)-fPhiBinWidth)>pow(10.,-6.))
{
cout<<endl;
if(fWeightsList->FindObject("pt_weights"))
{
fPtWeights = dynamic_cast<TH1D*>(fWeightsList->FindObject("pt_weights"));
+ if(!fPtWeights){printf("\n WARNING (GFC): !fPtWeights !!!!\n");exit(0);}
if(TMath::Abs(fPtWeights->GetBinWidth(1)-fPtBinWidth)>pow(10.,-6.))
{
cout<<endl;
if(fWeightsList->FindObject("eta_weights"))
{
fEtaWeights = dynamic_cast<TH1D*>(fWeightsList->FindObject("eta_weights"));
+ if(!fEtaWeights){printf("\n WARNING (GFC): !fEtaWeights !!!!\n");exit(0);}
if(TMath::Abs(fEtaWeights->GetBinWidth(1)-fEtaBinWidth)>pow(10.,-6.))
{
cout<<endl;
if(fWeightsList->FindObject("phi_weights"))
{
fPhiWeights = dynamic_cast<TH1F*>(fWeightsList->FindObject("phi_weights"));
+ if(!fPhiWeights){printf("\n WARNING (FQD): !fPhiWeights !!!!\n");exit(0);}
if(TMath::Abs(fPhiWeights->GetBinWidth(1)-fPhiBinWidth)>pow(10.,-6.))
{
cout<<endl;
if(fWeightsList->FindObject("pt_weights"))
{
fPtWeights = dynamic_cast<TH1D*>(fWeightsList->FindObject("pt_weights"));
+ if(!fPtWeights){printf("\n WARNING (FQD): !fPtWeights !!!!\n");exit(0);}
if(TMath::Abs(fPtWeights->GetBinWidth(1)-fPtBinWidth)>pow(10.,-6.))
{
cout<<endl;
if(fWeightsList->FindObject("eta_weights"))
{
fEtaWeights = dynamic_cast<TH1D*>(fWeightsList->FindObject("eta_weights"));
+ if(!fEtaWeights){printf("\n WARNING (FQD): !fEtaWeights !!!!\n");exit(0);}
if(TMath::Abs(fEtaWeights->GetBinWidth(1)-fEtaBinWidth)>pow(10.,-6.))
{
cout<<endl;
//RP:
TH1F* fHistPtRP = NULL;
- if(fCommonHists)
+ if(fCommonHists && fCommonHists->GetHistPtRP())
{
fHistPtRP = fCommonHists->GetHistPtRP();
}
}
//POI:
- TH1F* fHistPtPOI = fCommonHists->GetHistPtPOI();
+ TH1F* fHistPtPOI = NULL;
+ if(fCommonHists && fCommonHists->GetHistPtPOI())
+ {
+ fHistPtPOI = fCommonHists->GetHistPtPOI();
+ }
Double_t dYieldPtPOI = 0.;
Double_t dVPOI = 0.;
Double_t dErrVPOI = 0.;
// Get pointers for histograms with particle weights.
TList *weightsList = dynamic_cast<TList*>(fHistList->FindObject("Weights"));
- if(weightsList) this->SetWeightsList(weightsList);
+ if(!weightsList){printf("\n WARNING (QC): weightsList is NULL in AFAWQC::GPFPWH() !!!!\n");exit(0);}
+ this->SetWeightsList(weightsList);
TString fUseParticleWeightsName = "fUseParticleWeightsQC"; // to be improved (hirdwired label QC)
fUseParticleWeightsName += fAnalysisLabel->Data();
TProfile *useParticleWeights = dynamic_cast<TProfile*>(weightsList->FindObject(fUseParticleWeightsName.Data()));