#include <TString.h>
#include <TObjString.h>
#include <TProcessID.h>
+#include <TMap.h>
#include "AliESDInputHandler.h"
#include "AliESDEvent.h"
+#include "AliVCuts.h"
#include "AliESD.h"
#include "AliRunTag.h"
#include "AliEventTag.h"
ClassImp(AliESDInputHandler)
+static Option_t *gESDDataType = "ESD";
+
//______________________________________________________________________________
AliESDInputHandler::AliESDInputHandler() :
AliInputEventHandler(),
fEvent(0x0),
- fBranches(""),
- fBranchesOn(""),
fAnalysisType(0),
+ fNEvents(0),
+ fHLTEvent(0x0),
+ fHLTTree(0x0),
+ fUseHLT(kFALSE),
+ fTagCutSumm(0x0),
fUseTags(kFALSE),
fChainT(0),
fTreeT(0),
//______________________________________________________________________________
AliESDInputHandler::AliESDInputHandler(const char* name, const char* title):
- AliInputEventHandler(name, title), fEvent(0x0), fBranches(""), fBranchesOn(""), fAnalysisType(0),
- fUseTags(kFALSE), fChainT(0), fTreeT(0), fRunTag(0)
+ AliInputEventHandler(name, title), fEvent(0x0), fAnalysisType(0),
+ fNEvents(0), fHLTEvent(0x0), fHLTTree(0x0), fUseHLT(kFALSE), fTagCutSumm(0x0), fUseTags(kFALSE), fChainT(0), fTreeT(0), fRunTag(0)
{
// Constructor
}
fTree = tree;
if (!fTree) return kFALSE;
+ fTree->GetEntry(0);
+
// Get pointer to ESD event
SwitchOffBranches();
SwitchOnBranches();
- if (fEvent) {
- delete fEvent;
- fEvent = 0;
- }
- fEvent = new AliESDEvent();
-
+ if (!fEvent) fEvent = new AliESDEvent();
fEvent->ReadFromTree(fTree);
+ fNEvents = fTree->GetEntries();
+
+
return kTRUE;
}
-Bool_t AliESDInputHandler::BeginEvent(Long64_t /*entry*/)
+Bool_t AliESDInputHandler::BeginEvent(Long64_t entry)
{
// Copy from old to new format if necessary
AliESD* old = ((AliESDEvent*) fEvent)->GetAliESDOld();
((AliESDEvent*)fEvent)->CopyFromOldESD();
old->Reset();
}
+
+ if (fHLTTree) {
+ fHLTTree->GetEntry(entry);
+ }
+
+ fNewEvent = kTRUE;
+ //
+ // Event selection
+ //
+ if (fEventCuts)
+ fIsSelected = fEventCuts->IsSelected(fEvent);
+
return kTRUE;
}
// Notify a directory change
AliInfo(Form("Directory change %s \n", path));
//
+ if (fUseHLT) {
+ // Get HLTesdTree from current file
+ TTree* cTree = fTree;
+ if (fTree->GetTree()) cTree = fTree->GetTree();
+ TFile* cFile = cTree->GetCurrentFile();
+ cFile->GetObject("HLTesdTree", fHLTTree);
+
+ if (fHLTTree) {
+ if (!fHLTEvent) fHLTEvent = new AliESDEvent();
+ fHLTEvent->ReadFromTree(fHLTTree);
+ }
+ }
+
if (!fUseTags) return (kTRUE);
Bool_t zip = kFALSE;
TString fileName(path);
if(fileName.Contains("#AliESDs.root")){
- fileName.ReplaceAll("#AliESDs.root", "");
zip = kTRUE;
}
else if (fileName.Contains("AliESDs.root")){
fileName.ReplaceAll("AliESDs.root", "");
}
else if(fileName.Contains("#AliAOD.root")){
- fileName.ReplaceAll("#AliAOD.root", "");
zip = kTRUE;
}
else if(fileName.Contains("AliAOD.root")){
}
else if(fileName.Contains("#galice.root")){
// For running with galice and kinematics alone...
- fileName.ReplaceAll("#galice.root", "");
zip = kTRUE;
}
else if(fileName.Contains("galice.root")){
}
printf("AliESDInputHandler::Notify() Path: %s\n", pathName.Data());
-
+
if (fRunTag) {
fRunTag->Clear();
} else {
}
delete fTreeT; fTreeT = 0;
-
+
if (fChainT) {
delete fChainT;
fChainT = 0;
const char* name = 0x0;
TString tagFilename;
if (zip) {
- TFile* file = TFile::Open(fileName.Data());
+ TFile* file = fTree->GetCurrentFile();
TArchiveFile* arch = file->GetArchive();
TObjArray* arr = arch->GetMembers();
TIter next(arr);
}
-void AliESDInputHandler::SwitchOffBranches() const {
- //
- // Switch of branches on user request
- TObjArray * tokens = fBranches.Tokenize(" ");
- Int_t ntok = tokens->GetEntries();
- for (Int_t i = 0; i < ntok; i++) {
- TString str = ((TObjString*) tokens->At(i))->GetString();
- if (str.Length() == 0)
- continue;
- fTree->SetBranchStatus(Form("%s%s%s","*", str.Data(), "*"), 0);
- AliInfo(Form("Branch %s switched off \n", str.Data()));
+
+Option_t *AliESDInputHandler::GetDataType() const
+{
+// Returns handled data type.
+ return gESDDataType;
+}
+
+Int_t AliESDInputHandler::GetNEventAcceptedInFile()
+{
+ // Get number of events in file accepted by the tag cuts
+ // return -1 if no info is available
+ if (!fTagCutSumm) {
+ TList *luo = fTree->GetUserInfo();
+ if (!luo) {
+ AliInfo(Form("No user info in input tree - no tag cut summary\n"));
+ return -1;
+ }
+ for (int iluo=0; iluo<luo->GetEntries(); iluo++) {
+ fTagCutSumm = dynamic_cast<TMap *>(luo->At(iluo));
+ if (fTagCutSumm) break;
}
+ if (!fTagCutSumm) {
+ AliInfo(Form("No tag summary map in input tree\n"));
+ return -1;
+ }
+ }
+
+ TObjString *ostr = 0;
+ if (fTagCutSumm->FindObject(fTree->GetCurrentFile()->GetName()))
+ ostr = (TObjString *) fTagCutSumm->GetValue(fTree->GetCurrentFile()->GetName());
+ else {
+ AliInfo(Form("No tag cut summary for file %s\n", fTree->GetCurrentFile()->GetName()));
+ return -1;
+ }
+ char *iTagInfo;
+ iTagInfo = strdup(ostr->GetString().Data());
+
+ Int_t iAcc = atoi(strtok(iTagInfo, ","));
+
+ AliInfo(Form("Got %i accepted events for file %s", iAcc, fTree->GetCurrentFile()->GetName()));
+
+ free(iTagInfo);
+
+ return iAcc;
}
+Int_t AliESDInputHandler::GetNEventRejectedInFile()
+{
+ // Get number of events in file rejected by the tag cuts
+ // return -1 if no info is available
+ if (!fTagCutSumm) {
+ TList *luo = fTree->GetUserInfo();
+ if (!luo) {
+ AliInfo(Form("No user info in input tree - no tag cut summary\n"));
+ return -1;
+ }
+ for (int iluo=0; iluo<luo->GetEntries(); iluo++) {
+ fTagCutSumm = dynamic_cast<TMap *>(luo->At(iluo));
+ if (fTagCutSumm) break;
+ }
+ if (!fTagCutSumm) {
+ AliInfo(Form("No tag summary map in input tree\n"));
+ return -1;
+ }
+ }
-void AliESDInputHandler::SwitchOnBranches() const {
- //
- // Switch of branches on user request
- TObjArray * tokens = fBranchesOn.Tokenize(" ");
- Int_t ntok = tokens->GetEntries();
-
- for (Int_t i = 0; i < ntok; i++) {
- TString str = ((TObjString*) tokens->At(i))->GetString();
- if (str.Length() == 0)
- continue;
- fTree->SetBranchStatus(Form("%s%s%s","*", str.Data(), "*"), 1);
- AliInfo(Form("Branch %s switched on \n", str.Data()));
+ TObjString *ostr = 0;
+ if (fTagCutSumm->FindObject(fTree->GetCurrentFile()->GetName()))
+ ostr = (TObjString *) fTagCutSumm->GetValue(fTree->GetCurrentFile()->GetName());
+ else {
+ AliInfo(Form("No tag cut summary for file %s\n", fTree->GetCurrentFile()->GetName()));
+ return -1;
+ }
+ char *iTagInfo;
+ iTagInfo = strdup(ostr->GetString().Data());
+
+ strtok(iTagInfo, ",");
+ Int_t iRej = atoi(strtok(NULL, ","));
+
+ AliInfo(Form("Got %i accepted events for file %s", iRej, fTree->GetCurrentFile()->GetName()));
+
+ free(iTagInfo);
+
+ return iRej;
+}
+Bool_t AliESDInputHandler::GetCutSummaryForChain(Int_t *aTotal, Int_t *aAccepted, Int_t *aRejected)
+{
+ // Get number of events in the full chain
+ // Count accepted and rejected events
+ // return kFALSE if no info is available
+ if (!fTagCutSumm) {
+ TList *luo = fTree->GetUserInfo();
+ if (!luo) {
+ AliInfo(Form("No user info in input tree - no tag cut summary\n"));
+ return kFALSE;
+ }
+ for (int iluo=0; iluo<luo->GetEntries(); iluo++) {
+ fTagCutSumm = dynamic_cast<TMap *>(luo->At(iluo));
+ if (fTagCutSumm) break;
+ }
+ if (!fTagCutSumm) {
+ AliInfo(Form("No tag summary map in input tree\n"));
+ return kFALSE;
+ }
+ }
+
+ TMapIter *tIter = new TMapIter(fTagCutSumm);
+
+ Int_t iTotList=0, iAccList=0, iRejList=0;
+
+ TObject *cobj;
+ while ((cobj = tIter->Next())) {
+ TObjString *kstr = (TObjString *) cobj;
+ TObjString *vstr = (TObjString *) fTagCutSumm->GetValue(kstr->GetString().Data());
+ // printf("Got object value %s %s\n", kstr->GetString().Data(), vstr->GetString().Data());
+ char *iTagInfo;
+ iTagInfo = strdup(vstr->GetString().Data());
+
+ Int_t iAcc = atoi(strtok(iTagInfo, ","));
+ Int_t iRej = atoi(strtok(NULL, ","));
+
+ iAccList += iAcc;
+ iRejList += iRej;
+ iTotList += (iAcc+iRej);
}
+
+ *aTotal = iTotList;
+ *aAccepted = iAccList;
+ *aRejected = iRejList;
+
+ return kTRUE;
+}
+
+Int_t AliESDInputHandler::GetNFilesEmpty()
+{
+ // Count number of files in which all events were de-selected
+ // For such files Notify() will NOT be called
+ // return -1 if no info is available
+ if (!fTagCutSumm) {
+ TList *luo = fTree->GetUserInfo();
+ if (!luo) {
+ AliInfo(Form("No user info in input tree - no tag cut summary\n"));
+ return -1;
+ }
+ for (int iluo=0; iluo<luo->GetEntries(); iluo++) {
+ fTagCutSumm = dynamic_cast<TMap *>(luo->At(iluo));
+ if (fTagCutSumm) break;
+ }
+ if (!fTagCutSumm) {
+ AliInfo(Form("No tag summary map in input tree\n"));
+ return -1;
+ }
+ }
+
+ TMapIter *tIter = new TMapIter(fTagCutSumm);
+
+ Int_t iFilesEmpty = 0;
+
+ TObject *cobj;
+ while ((cobj = tIter->Next())) {
+ TObjString *kstr = (TObjString *) cobj;
+ TObjString *vstr = (TObjString *) fTagCutSumm->GetValue(kstr->GetString().Data());
+ // printf("Got object value %s %s\n", kstr->GetString().Data(), vstr->GetString().Data());
+ char *iTagInfo;
+ iTagInfo = strdup(vstr->GetString().Data());
+
+ Int_t iAcc = atoi(strtok(iTagInfo, ","));
+ Int_t iRej = atoi(strtok(NULL, ","));
+
+ if ((iAcc == 0) && ((iRej+iAcc)>0))
+ iFilesEmpty++;
+ }
+
+ return iFilesEmpty;
+
}