// position range. Fail otherwise
int mult = event->NumberOfTracks();
double vertexZPos = event->PrimVertPos().z();
- cout << "AliFemtoBasicEventCut:: mult: " << fEventMult[0] << " < " << mult << " < " << fEventMult[1] << endl;
- cout << "AliFemtoBasicEventCut:: VertexZPos: " << fVertZPos[0] << " < " << vertexZPos << " < " << fVertZPos[1] << endl;
- cout << "AliFemtoBasicEventCut:: VertexZErr: " << event->PrimVertCov()[4] << endl;
+// cout << "AliFemtoBasicEventCut:: mult: " << fEventMult[0] << " < " << mult << " < " << fEventMult[1] << endl;
+// cout << "AliFemtoBasicEventCut:: VertexZPos: " << fVertZPos[0] << " < " << vertexZPos << " < " << fVertZPos[1] << endl;
+// cout << "AliFemtoBasicEventCut:: VertexZErr: " << event->PrimVertCov()[4] << endl;
bool goodEvent =
((mult >= fEventMult[0]) &&
(mult <= fEventMult[1]) &&
(vertexZPos < fVertZPos[1]) &&
(fAcceptBadVertex || (event->PrimVertCov()[4] > -1000.0)));
goodEvent ? fNEventsPassed++ : fNEventsFailed++ ;
- cout << "AliFemtoBasicEventCut:: return : " << goodEvent << endl;
+// cout << "AliFemtoBasicEventCut:: return : " << goodEvent << endl;
return (goodEvent);
}
//------------------------------
if ( iy<0 || iy >= fBinsy) return 0;
if ( iz<0 || iz >= fBinsz) return 0;
int bin = ix + iy*fBinsx + iz*fBinsy*fBinsx;
- cout << " AliFemtoPicoEventCollectionVectorHideAway::PicoEventCollection(...) - bin(ix,iy,iz): ";
- cout << bin << "(" << ix <<"," << iy << "," << iz <<")" << endl;
+// cout << " AliFemtoPicoEventCollectionVectorHideAway::PicoEventCollection(...) - bin(ix,iy,iz): ";
+// cout << bin << "(" << ix <<"," << iy << "," << iz <<")" << endl;
return fCollectionVector[bin];
}
// -----------------------------------
return *this;
}
+unsigned int AliFemtoPicoEventCollectionVectorHideAway::GetBinXNumber(double x) { return (int)floor( (x-fMinx)/fStepx ); }
+unsigned int AliFemtoPicoEventCollectionVectorHideAway::GetBinYNumber(double y) { return (int)floor( (y-fMiny)/fStepy ); }
+unsigned int AliFemtoPicoEventCollectionVectorHideAway::GetBinZNumber(double z) { return (int)floor( (z-fMinz)/fStepz ); }
AliFemtoPicoEventCollection* PicoEventCollection(int bx, int by, int bz);
AliFemtoPicoEventCollection* PicoEventCollection(double x, double y=0, double z=0);
+ unsigned int GetBinXNumber(double x);
+ unsigned int GetBinYNumber(double y);
+ unsigned int GetBinZNumber(double z);
private:
int fBinsTot; // Total number of bins
int fBinsx,fBinsy,fBinsz; // Number of bins on x, y, z axis
bool tmpPassEvent = fEventCut->Pass(hbtEvent);
fEventCut->FillCutMonitor(hbtEvent, tmpPassEvent);
if (tmpPassEvent) {
- cout << "AliFemtoSimpleAnalysis::ProcessEvent() - Event has passed cut - build picoEvent from " <<
- hbtEvent->TrackCollection()->size() << " tracks in TrackCollection" << endl;
+// cout << "AliFemtoSimpleAnalysis::ProcessEvent() - Event has passed cut - build picoEvent from " <<
+// hbtEvent->TrackCollection()->size() << " tracks in TrackCollection" << endl;
+// cout << "Event has passed cut with " << hbtEvent->TrackCollection()->size() << " tracks" << endl;
// OK, analysis likes the event-- build a pico event from it, using tracks the analysis likes...
fPicoEvent = new AliFemtoPicoEvent; // this is what we will make pairs from and put in Mixing Buffer
// no memory leak. we will delete picoevents when they come out of the mixing buffer
FillHbtParticleCollection(fFirstParticleCut,(AliFemtoEvent*)hbtEvent,fPicoEvent->FirstParticleCollection());
if ( !(AnalyzeIdenticalParticles()) )
FillHbtParticleCollection(fSecondParticleCut,(AliFemtoEvent*)hbtEvent,fPicoEvent->SecondParticleCollection());
- cout <<"AliFemtoSimpleAnalysis::ProcessEvent - #particles in First, Second Collections: " <<
+// cout <<"AliFemtoSimpleAnalysis::ProcessEvent - #particles in First, Second Collections: " <<
+// fPicoEvent->FirstParticleCollection()->size() << " " <<
+// fPicoEvent->SecondParticleCollection()->size() << endl;
+
+ cout << "#particles in Collection 1, 2: " <<
fPicoEvent->FirstParticleCollection()->size() << " " <<
fPicoEvent->SecondParticleCollection()->size() << endl;