// Remark: cumulants calculated here are NOT corrected for non-uniform acceptance. This correction is applied in the method ...
// to be improved (description)
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
{
// calculate final results for integrated flow of RPs and POIs
- Int_t typeFlag = -1;
+ Int_t typeFlag = 0;
if(type == "RP")
{
{
// fill common result histograms for differential flow
- Int_t typeFlag = -1;
- //Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ //Int_t ptEtaFlag = 0;
if(type == "RP")
{
// 2: <<6'>>
// 3: <<8'>>
- Int_t t = -1; // type flag
- Int_t pe = -1; // ptEta flag
+ Int_t t = 0; // type flag
+ Int_t pe = 0; // ptEta flag
if(type == "RP")
{
// Calculate sums of various event weights for reduced correlations.
// (These quantitites are needed in expressions for unbiased estimators relevant for the statistical errors.)
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
// Table:
// [0=<2>,1=<2'>,2=<4>,3=<4'>,4=<6>,5=<6'>,6=<8>,7=<8'>] x [0=<2>,1=<2'>,2=<4>,3=<4'>,4=<6>,5=<6'>,6=<8>,7=<8'>]
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
{
// Transfer profiles into histograms and calculate statistical errors correctly.
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
// <4'><8>, <4'><8'>, <6><6'>, <6><8'>, <6'><8>,
// <6'><8'>, <8><8'>.
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
// 5th bin: Cov(<2'>,<4'>) * (sum_{i=1}^{N} w_{<2'>}_i w_{<4'>}_i )/[(sum_{i=1}^{N} w_{<2'>}_i) * (sum_{j=1}^{N} w_{<4'>}_j)]
// ...
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
{
// calculate differential flow from differential cumulants and previously obtained integrated flow: (to be improved: description)
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
{
// Calculate all correlations needed for differential flow using particle weights.
- Int_t t = -1; // type flag
- Int_t pe = -1; // ptEta flag
+ Int_t t = 0; // type flag
+ Int_t pe = 0; // ptEta flag
if(type == "RP")
{
//Double_t dImQ3n = (*fImQ)(2,0);
//Double_t dImQ4n = (*fImQ)(3,0);
- Int_t t = -1; // type flag
- Int_t pe = -1; // ptEta flag
+ Int_t t = 0; // type flag
+ Int_t pe = 0; // ptEta flag
if(type == "RP")
{
//Double_t dImQ3n = (*fImQ)(2,0);
//Double_t dImQ4n = (*fImQ)(3,0);
- Int_t t = -1; // type flag
- Int_t pe = -1; // ptEta flag
+ Int_t t = 0; // type flag
+ Int_t pe = 0; // ptEta flag
if(type == "RP")
{
// to be improved: debugged - I do not correctly transfer all profiles into histos (bug appears only after merging)
- Int_t t = -1; // type flag
- Int_t pe = -1; // ptEta flag
+ Int_t t = 0; // type flag
+ Int_t pe = 0; // ptEta flag
if(type == "RP")
{
{
// Calculate generalized differential flow Q-cumulants (corrected for non-uniform acceptance)
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
// to be improved (rewritten completely)
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
// <4'> = <cos(n*(psi1+phi2-phi3-phi4))>
// ...
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
typeFlag = 0;
{
// Compare correlations needed for diff. flow calculated with nested loops and those calculated from Q-vectors
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
typeFlag = 0;
// <4'> = <w2 w3 w4 cos(n*(psi1+phi2-phi3-phi4))>
// ...
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
typeFlag = 0;
// 5:
// 6:
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
typeFlag = 0;
{
// Compare corrections temrs for non-uniform acceptance needed for diff. flow calculated with nested loops and those calculated from Q-vectors
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
typeFlag = 0;
Double_t dSM1p2k = (*fSMpk)(0,2);
Double_t dSM2p1k = (*fSMpk)(1,1);
- Int_t t = -1; // type flag
- Int_t pe = -1; // ptEta flag
+ Int_t t = 0; // type flag
+ Int_t pe = 0; // ptEta flag
if(type == "RP")
{
Double_t dSM1p2k = (*fSMpk)(0,2);
Double_t dSM2p1k = (*fSMpk)(1,1);
- Int_t t = -1; // type flag
- Int_t pe = -1; // ptEta flag
+ Int_t t = 0; // type flag
+ Int_t pe = 0; // ptEta flag
if(type == "RP")
{
// 5:
// 6:
- Int_t typeFlag = -1;
- Int_t ptEtaFlag = -1;
+ Int_t typeFlag = 0;
+ Int_t ptEtaFlag = 0;
if(type == "RP")
{
typeFlag = 0;