if(trackpair == 0x0) return;
Double_t q = trackpair->GetQInv();
- Bool_t fill = kFALSE;
Double_t weight = partpair->GetWeight();
fNumerator->Fill(q,weight);
- if ( (q < 0.15) && (fNPairsFitArea < 2.e+5))
- {
- fNPairsFitArea++;
- fill = kTRUE;
- }
- if ( (q > 0.15) && (q < 0.3) && (fNPairsFitArea < 1.e+5))
- {
- fNPairsNormArea++;
- fill = kTRUE;
- }
-
- if (fill)
- {
- const AliVAODParticle& p1 = *(trackpair->Particle1());
- const AliVAODParticle& p2 = *(trackpair->Particle2());
- fNtuple->Fill(p1.Px(),p1.Py(),p1.Pz(),p1.E(),
- p2.Px(),p2.Py(),p2.Pz(),p2.E());
- }
}
/****************************************************************/
trackpair = CheckPair(trackpair);
// partpair = CheckPair(partpair);
if ( trackpair && partpair)
- {
- fDenominator->Fill(trackpair->GetQInv());
- }
+ {
+ Double_t q = trackpair->GetQInv();
+
+ Bool_t fill = kFALSE;
+
+ if ( (q < 0.15) && (fNPairsFitArea < 2.e+5))
+ {
+ fNPairsFitArea++;
+ fill = kTRUE;
+ }
+
+ if ( (q > 0.15) && (q < 0.3) && (fNPairsFitArea < 1.e+5))
+ {
+ fNPairsNormArea++;
+ fill = kTRUE;
+ }
+
+ if (fill)
+ {
+ const AliVAODParticle& p1 = *(trackpair->Particle1());
+ const AliVAODParticle& p2 = *(trackpair->Particle2());
+ fNtuple->Fill(p1.Px(),p1.Py(),p1.Pz(),p1.E(),
+ p2.Px(),p2.Py(),p2.Pz(),p2.E());
+ }
+
+ fDenominator->Fill(q);
+ }
}
/*****************************************************************/