virtual Long64_t Merge(TCollection* list);
virtual void Copy(TObject &c) const;
+ AliESDtrackCuts(const AliESDtrackCuts& pd); // Copy Constructor
+ AliESDtrackCuts &operator=(const AliESDtrackCuts &c);
//######################################################
// track quality cut setters
TH1F* fhCutStatistics; //-> statistics of what cuts the tracks did not survive
TH2F* fhCutCorrelation; //-> 2d statistics plot
- private:
-
- AliESDtrackCuts(const AliESDtrackCuts& pd); // Copy Constructor
- AliESDtrackCuts &operator=(const AliESDtrackCuts &c);
-
-
ClassDef(AliESDtrackCuts, 1)
};
// read the user objects, called from slavebegin and begin
if (!fEsdTrackCutsAll && fInput)
- fEsdTrackCutsAll = dynamic_cast<AliESDtrackCuts*> (fInput->FindObject("esdTrackCutsAll"));
+ fEsdTrackCutsAll = dynamic_cast<AliESDtrackCuts*> (fInput->FindObject("esdTrackCutsAll")->Clone());
if (!fEsdTrackCutsPri && fInput)
- fEsdTrackCutsPri = dynamic_cast<AliESDtrackCuts*> (fInput->FindObject("esdTrackCutsPri"));
+ fEsdTrackCutsPri = dynamic_cast<AliESDtrackCuts*> (fInput->FindObject("esdTrackCutsPri")->Clone());
if (!fEsdTrackCutsSec && fInput)
- fEsdTrackCutsSec = dynamic_cast<AliESDtrackCuts*> (fInput->FindObject("esdTrackCutsSec"));
+ fEsdTrackCutsSec = dynamic_cast<AliESDtrackCuts*> (fInput->FindObject("esdTrackCutsSec")->Clone());
if (!fEsdTrackCutsAll && tree)
fEsdTrackCutsAll = dynamic_cast<AliESDtrackCuts*> (tree->GetUserInfo()->FindObject("esdTrackCutsAll"));
AliDebug(AliLog::kError, Form("UNEXPECTED: part with label %d not found in stack (track loop).", label));
continue;
}
- if (AliPWG0Helper::IsPrimaryCharged(particle, nPrim) == kFALSE)
+
+ if (label < nPrim)
fEsdTrackCutsPri->AcceptTrack(esdTrack);
else
fEsdTrackCutsSec->AcceptTrack(esdTrack);
AliSelectorRL::Terminate();
+ if (fOutput)
+ fOutput->Print();
+
+ fEsdTrackCutsAll = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("esdTrackCutsAll"));
+ fEsdTrackCutsPri = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("esdTrackCutsPri"));
+ fEsdTrackCutsSec = dynamic_cast<AliESDtrackCuts*> (fOutput->FindObject("esdTrackCutsSec"));
+
+ // check if the esd track cut objects are there
+ if (!fEsdTrackCutsAll || !fEsdTrackCutsPri || !fEsdTrackCutsSec) {
+ AliDebug(AliLog::kError, Form("fEsdTrackCutsXXX not available %p %p %p", fEsdTrackCutsAll, fEsdTrackCutsPri, fEsdTrackCutsSec));
+ return;
+ }
+
TFile* file = TFile::Open("trackCuts.root", "RECREATE");
fEsdTrackCutsAll->SaveHistograms("esdTrackCutsAll");
fEsdTrackCutsPri->SaveHistograms("esdTrackCutsPri");