\r
// vector holding the charges/kinematics of all tracks (charge,y,eta,phi,p0,p1,p2,pt,E)\r
vector<Double_t> *chargeVector[9]; // original charge\r
+ vector<Double_t> *chargeVectorShuffled[9]; // shuffled charge\r
+\r
for(Int_t i = 0; i < 9; i++){\r
- chargeVector[i] = new vector<Double_t>;\r
+ chargeVector[i] = new vector<Double_t>;\r
+ chargeVectorShuffled[i] = new vector<Double_t>;\r
}\r
\r
Double_t v_charge;\r
chargeVector[6]->push_back(v_p[2]);\r
chargeVector[7]->push_back(v_pt);\r
chargeVector[8]->push_back(v_E);\r
+\r
+ if(fRunShuffling) {\r
+ chargeVectorShuffled[0]->push_back(v_charge);\r
+ chargeVectorShuffled[1]->push_back(v_y);\r
+ chargeVectorShuffled[2]->push_back(v_eta);\r
+ chargeVectorShuffled[3]->push_back(v_phi);\r
+ chargeVectorShuffled[4]->push_back(v_p[0]);\r
+ chargeVectorShuffled[5]->push_back(v_p[1]);\r
+ chargeVectorShuffled[6]->push_back(v_p[2]);\r
+ chargeVectorShuffled[7]->push_back(v_pt);\r
+ chargeVectorShuffled[8]->push_back(v_E);\r
+ }\r
\r
} //track loop\r
\r
- // calculate balance function\r
+ // calculate balance function\r
fBalance->FillBalance(fCentrality,chargeVector);\r
+ \r
+ // calculate shuffled balance function\r
+ if(fRunShuffling) {\r
+ random_shuffle(chargeVectorShuffled[0]->begin(), chargeVectorShuffled[0]->end());\r
+ fShuffledBalance->FillBalance(fCentrality,chargeVectorShuffled);\r
+ }\r
+\r
// clean charge vector afterwards\r
for(Int_t i = 0; i < 9; i++){ \r
chargeVector[i]->clear();\r
+ chargeVectorShuffled[i]->clear();\r
}\r
\r
}//Vz cut\r