- fXEta=new Double_t[fNin]; fXPhi=new Double_t[fNin];
- fPx = new TVectorD(fNin);
- for (Int_t iIn=0; iIn<fNin; iIn++){
- TLorentzVector *lv=(TLorentzVector*)lvArray->At(iIn);
- fXEta[iIn] = lv->Eta();
- fXPhi[iIn] = lv->Phi()<0 ? lv->Phi() + 2*TMath::Pi() : lv->Phi();
- (*fPx)(iIn)=lv->Pt();
- dEtSum+=(*fPx)(iIn);
+ TClonesArray *lvArray = fReader->GetMomentumArray();
+ Int_t nEntr = lvArray->GetEntries();
+ fNin=0;
+ for (Int_t iEn=0; iEn<nEntr; iEn++) if (fReader->GetCutFlag(iEn)==1) fNin++;
+ Double_t *xEta = new Double_t[fNin];
+ Double_t *xPhi = new Double_t[fNin];
+ xData[0]=xEta; xData[1]=xPhi;
+ vPx->ResizeTo(fNin);
+ Int_t iIn=0;
+ for (Int_t iEn=0; iEn<nEntr; iEn++){
+ if (fReader->GetCutFlag(iEn)==0) continue;
+ TLorentzVector *lv=(TLorentzVector*)lvArray->At(iEn);
+ xEta[iIn] = lv->Eta();
+ xPhi[iIn] = lv->Phi()<0 ? lv->Phi() + 2*TMath::Pi() : lv->Phi();
+ (*vPx)(iIn)=lv->Pt();
+ dEtSum+=(*vPx)(iIn);
+ iIn++;