- pVect2[1]=(fEvt+en2)->fTracks[j].fP[0];
- pVect2[2]=(fEvt+en2)->fTracks[j].fP[1];
- pVect2[3]=(fEvt+en2)->fTracks[j].fP[2];
- ch2 = Int_t(((fEvt+en2)->fTracks[j].fCharge + 1)/2.);
-
- for (Int_t k=j+1; k<(fEvt+en3)->fNtracks; k++) {// 3rd particle
- if(en3==0) {
- if(!NormQPairSwitch_E0E0[i][k]) continue;
- if(!NormQPairSwitch_E0E0[j][k]) continue;
- }else if(en3==1){
- if(!NormQPairSwitch_E0E1[i][k]) continue;
- if(!NormQPairSwitch_E0E1[j][k]) continue;
+ pVect3[1]=(fEvt+en3)->fTracks[k].fP[0];
+ pVect3[2]=(fEvt+en3)->fTracks[k].fP[1];
+ pVect3[3]=(fEvt+en3)->fTracks[k].fP[2];
+ ch3 = Int_t(((fEvt+en3)->fTracks[k].fCharge + 1)/2.);
+ Bool_t fill2=kFALSE, fill3=kFALSE, fill4=kFALSE;
+ SetFillBins3(ch1, ch2, ch3, 1, bin1, bin2, bin3, fill2, fill3, fill4);
+
+ Float_t KT3 = sqrt(pow(pVect1[1]+pVect2[1]+pVect3[1],2) + pow(pVect1[2]+pVect2[2]+pVect3[2],2))/3.;
+ if(KT3<=fKT3transition) KT3index=0;
+ else KT3index=1;
+
+ if(en2==0 && en3==0 && en4==0) Charge1[bin1].Charge2[bin2].Charge3[bin3].MB[fMbin].EDB[KT3index].ThreePT[0].fNorm3->Fill(0);
+ if(en2==1 && en3==2 && en4==3) Charge1[bin1].Charge2[bin2].Charge3[bin3].MB[fMbin].EDB[KT3index].ThreePT[4].fNorm3->Fill(0);
+ if(en2==0 && en3==1 && en4==2) {
+ if(fill2) Charge1[bin1].Charge2[bin2].Charge3[bin3].MB[fMbin].EDB[KT3index].ThreePT[1].fNorm3->Fill(0);
+ if(fill3) Charge1[bin1].Charge2[bin2].Charge3[bin3].MB[fMbin].EDB[KT3index].ThreePT[2].fNorm3->Fill(0);
+ if(fill4) Charge1[bin1].Charge2[bin2].Charge3[bin3].MB[fMbin].EDB[KT3index].ThreePT[3].fNorm3->Fill(0);
+ }
+
+
+ for (Int_t l=k+1; l<(fEvt+en4)->fNtracks; l++) {// 4th particle
+ /*if(en4==0){
+ if(!NormQPairSwitch_E0E0[i][l]) continue;
+ if(!NormQPairSwitch_E0E0[j][l]) continue;
+ if(!NormQPairSwitch_E0E0[k][l]) continue;
+ }else if(en4==1){
+ if(!NormQPairSwitch_E0E1[i][l]) continue;
+ if(!NormQPairSwitch_E0E1[j][l]) continue;
+ if(!NormQPairSwitch_E0E1[k][l]) continue;
+ }else if(en4==2){
+ if(!NormQPairSwitch_E0E2[i][l]) continue;
+ if(!NormQPairSwitch_E0E2[j][l]) continue;
+ if(!NormQPairSwitch_E1E2[k][l]) continue;
+ }else{
+ if(!NormQPairSwitch_E0E3[i][l]) continue;
+ if(!NormQPairSwitch_E1E3[j][l]) continue;
+ if(!NormQPairSwitch_E2E3[k][l]) continue;
+ }*/
+ if(en4==0){
+ if(fNormQPairSwitch_E0E0[i]->At(l)=='0') continue;
+ if(fNormQPairSwitch_E0E0[j]->At(l)=='0') continue;
+ if(fNormQPairSwitch_E0E0[k]->At(l)=='0') continue;
+ }else if(en4==1){
+ if(fNormQPairSwitch_E0E1[i]->At(l)=='0') continue;
+ if(fNormQPairSwitch_E0E1[j]->At(l)=='0') continue;
+ if(fNormQPairSwitch_E0E1[k]->At(l)=='0') continue;
+ }else if(en4==2){
+ if(fNormQPairSwitch_E0E2[i]->At(l)=='0') continue;
+ if(fNormQPairSwitch_E0E2[j]->At(l)=='0') continue;
+ if(fNormQPairSwitch_E1E2[k]->At(l)=='0') continue;