[u/mrichter/AliRoot.git] / STEER / AliESDInputHandler.cxx

/**************************************************************************
 * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */

//-------------------------------------------------------------------------
//     Event handler for ESD input 
//     Author: Andreas Morsch, CERN
//-------------------------------------------------------------------------

#include <TTree.h>
#include <TChain.h>
#include <TFile.h>
#include <TArchiveFile.h>
#include <TObjArray.h>
#include <TSystem.h>
#include <TString.h>
#include <TObjString.h>
#include <TProcessID.h>
#include <TMap.h>

#include "AliESDInputHandler.h"
#include "AliESDEvent.h"
#include "AliESD.h"
#include "AliRunTag.h"
#include "AliEventTag.h"
#include "AliLog.h"

ClassImp(AliESDInputHandler)

static Option_t *gESDDataType = "ESD";

//______________________________________________________________________________
AliESDInputHandler::AliESDInputHandler() :
  AliInputEventHandler(),
  fEvent(0x0),
  fAnalysisType(0),
  fNEvents(0),
  fHLTEvent(0x0),
  fHLTTree(0x0),
  fUseHLT(kFALSE),
  fTagCutSumm(0x0),
  fUseTags(kFALSE),
  fChainT(0),
  fTreeT(0),
  fRunTag(0)
{
  // default constructor
}

//______________________________________________________________________________
AliESDInputHandler::~AliESDInputHandler() 
{
  //  destructor
  //  delete fEvent;
}

//______________________________________________________________________________
AliESDInputHandler::AliESDInputHandler(const char* name, const char* title):
    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
}

Bool_t AliESDInputHandler::Init(TTree* tree,  Option_t* opt)
{
    // Initialisation necessary for each new tree 
    fAnalysisType = opt;
    fTree         = tree;
    
    if (!fTree) return kFALSE;
    fTree->GetEntry(0);
    
    // Get pointer to ESD event
    SwitchOffBranches();
    SwitchOnBranches();
    
    if (!fEvent) fEvent = new AliESDEvent();
    fEvent->ReadFromTree(fTree);
    fNEvents = fTree->GetEntries();


    return kTRUE;
}

Bool_t AliESDInputHandler::BeginEvent(Long64_t entry)
{
    // Copy from old to new format if necessary
  AliESD* old = ((AliESDEvent*) fEvent)->GetAliESDOld();
  if (old) {
	((AliESDEvent*)fEvent)->CopyFromOldESD();
	old->Reset();
  }

  if (fHLTTree) {
      fHLTTree->GetEntry(entry);
  }
  
  return kTRUE;
}

Bool_t  AliESDInputHandler::FinishEvent()
{
    // Finish the event 
    if(fEvent)fEvent->Reset();
    return kTRUE;
} 

Bool_t AliESDInputHandler::Notify(const char* path)
{
    // 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")){
	zip = kTRUE;
    } 
    else if (fileName.Contains("AliESDs.root")){
	fileName.ReplaceAll("AliESDs.root", "");
    }
    else if(fileName.Contains("#AliAOD.root")){
	zip = kTRUE;
    }
    else if(fileName.Contains("AliAOD.root")){
	fileName.ReplaceAll("AliAOD.root", "");
    }
    else if(fileName.Contains("#galice.root")){
	// For running with galice and kinematics alone...
	zip = kTRUE;
    }
    else if(fileName.Contains("galice.root")){
	// For running with galice and kinematics alone...
	fileName.ReplaceAll("galice.root", "");
    }

    
    TString pathName("./");
    if (fileName.Length() != 0) {
	pathName = fileName;
    }
    
    printf("AliESDInputHandler::Notify() Path: %s\n", pathName.Data());

    if (fRunTag) {
	fRunTag->Clear();
    } else {
	fRunTag = new AliRunTag();
    }
    
    delete fTreeT; fTreeT = 0;
    
    if (fChainT) {
	delete fChainT;
	fChainT = 0;
    }
    
    if (!fChainT) {
	fChainT = new TChain("T");
    }
    

    const char* tagPattern = "ESD.tag.root";
    const char* name = 0x0;
    TString tagFilename;
    if (zip) {
	TFile* file = fTree->GetCurrentFile();
	TArchiveFile* arch = file->GetArchive();
	TObjArray* arr = arch->GetMembers();
	TIter next(arr);
	
	while ((file = (TFile*) next())) {
	    name = file->GetName();
	    if (strstr(name,tagPattern)) { 
		tagFilename = pathName.Data();
		tagFilename += "#";
		tagFilename += name;
		fChainT->Add(tagFilename);  
		AliInfo(Form("Adding %s to tag chain \n", tagFilename.Data()));
	    }//pattern check
	} // archive file loop
    } else {
	void * dirp = gSystem->OpenDirectory(pathName.Data());
	while((name = gSystem->GetDirEntry(dirp))) {
	    if (strstr(name,tagPattern)) { 
		tagFilename = pathName.Data();
		tagFilename += "/";
		tagFilename += name;
		fChainT->Add(tagFilename);  
		AliInfo(Form("Adding %s to tag chain \n", tagFilename.Data()));
	    }//pattern check
	}//directory loop
    }
    fChainT->SetBranchAddress("AliTAG",&fRunTag);
    fChainT->GetEntry(0);
    return kTRUE;
}


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;
    }
  }

  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;
  
}
Commit	Line	Data
d3dd3d14	1	/**************************************************************************
	2	* Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
	16	/* $Id$ */
	17
	18	//-------------------------------------------------------------------------
	19	// Event handler for ESD input
	20	// Author: Andreas Morsch, CERN
	21	//-------------------------------------------------------------------------
	22
933fd60f	23	#include <TTree.h>
1aa09fd1	24	#include <TChain.h>
	25	#include <TFile.h>
	26	#include <TArchiveFile.h>
	27	#include <TObjArray.h>
	28	#include <TSystem.h>
300d5701	29	#include <TString.h>
300d5701	30	#include <TObjString.h>
80e1b60a	31	#include <TProcessID.h>
f0f29480	32	#include <TMap.h>
d3dd3d14	33
	34	#include "AliESDInputHandler.h"
	35	#include "AliESDEvent.h"
	36	#include "AliESD.h"
1aa09fd1	37	#include "AliRunTag.h"
1aa09fd1	38	#include "AliEventTag.h"
300d5701	39	#include "AliLog.h"
d3dd3d14	40
	41	ClassImp(AliESDInputHandler)
	42
c2b6979d	43	static Option_t *gESDDataType = "ESD";
c2b6979d	44
d3dd3d14	45	//______________________________________________________________________________
6989bff3	46	AliESDInputHandler::AliESDInputHandler() :
6989bff3	47	AliInputEventHandler(),
0cd61c1d	48	fEvent(0x0),
1aa09fd1	49	fAnalysisType(0),
4d2a7370	50	fNEvents(0),
6ccd727b	51	fHLTEvent(0x0),
	52	fHLTTree(0x0),
	53	fUseHLT(kFALSE),
f0f29480	54	fTagCutSumm(0x0),
1aa09fd1	55	fUseTags(kFALSE),
	56	fChainT(0),
	57	fTreeT(0),
	58	fRunTag(0)
d3dd3d14	59	{
	60	// default constructor
	61	}
	62
	63	//______________________________________________________________________________
	64	AliESDInputHandler::~AliESDInputHandler()
	65	{
1aa09fd1	66	// destructor
6989bff3	67	// delete fEvent;
d3dd3d14	68	}
	69
	70	//______________________________________________________________________________
	71	AliESDInputHandler::AliESDInputHandler(const char* name, const char* title):
f16ba708	72	AliInputEventHandler(name, title), fEvent(0x0), fAnalysisType(0),
f0f29480	73	fNEvents(0), fHLTEvent(0x0), fHLTTree(0x0), fUseHLT(kFALSE), fTagCutSumm(0x0), fUseTags(kFALSE), fChainT(0), fTreeT(0), fRunTag(0)
d3dd3d14	74	{
e12e402c	75	// Constructor
d3dd3d14	76	}
d3dd3d14	77
1aa09fd1	78	Bool_t AliESDInputHandler::Init(TTree* tree, Option_t* opt)
d3dd3d14	79	{
300d5701	80	// Initialisation necessary for each new tree
1aa09fd1	81	fAnalysisType = opt;
1aa09fd1	82	fTree = tree;
300d5701	83
6073f8c9	84	if (!fTree) return kFALSE;
6ccd727b	85	fTree->GetEntry(0);
6ccd727b	86
d3dd3d14	87	// Get pointer to ESD event
300d5701	88	SwitchOffBranches();
3858dd64	89	SwitchOnBranches();
300d5701	90
efea3f54	91	if (!fEvent) fEvent = new AliESDEvent();
60996693	92	fEvent->ReadFromTree(fTree);
4d2a7370	93	fNEvents = fTree->GetEntries();
6ccd727b	94
50aaaf00	95
d3dd3d14	96	return kTRUE;
	97	}
	98
6ccd727b	99	Bool_t AliESDInputHandler::BeginEvent(Long64_t entry)
d3dd3d14	100	{
d3dd3d14	101	// Copy from old to new format if necessary
d7749dec	102	AliESD* old = ((AliESDEvent*) fEvent)->GetAliESDOld();
d7749dec	103	if (old) {
d3dd3d14	104	((AliESDEvent*)fEvent)->CopyFromOldESD();
d3dd3d14	105	old->Reset();
d7749dec	106	}
6ccd727b	107
	108	if (fHLTTree) {
	109	fHLTTree->GetEntry(entry);
	110	}
	111
d7749dec	112	return kTRUE;
d3dd3d14	113	}
d3dd3d14	114
e12e402c	115	Bool_t AliESDInputHandler::FinishEvent()
	116	{
	117	// Finish the event
	118	if(fEvent)fEvent->Reset();
	119	return kTRUE;
6989bff3	120	}
6989bff3	121
1aa09fd1	122	Bool_t AliESDInputHandler::Notify(const char* path)
	123	{
	124	// Notify a directory change
a5112f7d	125	AliInfo(Form("Directory change %s \n", path));
a5112f7d	126	//
50aaaf00	127	if (fUseHLT) {
	128	// Get HLTesdTree from current file
	129	TTree* cTree = fTree;
	130	if (fTree->GetTree()) cTree = fTree->GetTree();
	131	TFile* cFile = cTree->GetCurrentFile();
	132	cFile->GetObject("HLTesdTree", fHLTTree);
	133
	134	if (fHLTTree) {
	135	if (!fHLTEvent) fHLTEvent = new AliESDEvent();
	136	fHLTEvent->ReadFromTree(fHLTTree);
	137	}
	138	}
	139
1aa09fd1	140	if (!fUseTags) return (kTRUE);
	141
	142	Bool_t zip = kFALSE;
	143
	144	TString fileName(path);
	145	if(fileName.Contains("#AliESDs.root")){
1aa09fd1	146	zip = kTRUE;
	147	}
	148	else if (fileName.Contains("AliESDs.root")){
	149	fileName.ReplaceAll("AliESDs.root", "");
	150	}
	151	else if(fileName.Contains("#AliAOD.root")){
1aa09fd1	152	zip = kTRUE;
	153	}
	154	else if(fileName.Contains("AliAOD.root")){
	155	fileName.ReplaceAll("AliAOD.root", "");
	156	}
	157	else if(fileName.Contains("#galice.root")){
	158	// For running with galice and kinematics alone...
1aa09fd1	159	zip = kTRUE;
	160	}
	161	else if(fileName.Contains("galice.root")){
	162	// For running with galice and kinematics alone...
	163	fileName.ReplaceAll("galice.root", "");
	164	}
	165
1aa09fd1	166
be4763e2	167	TString pathName("./");
	168	if (fileName.Length() != 0) {
	169	pathName = fileName;
	170	}
	171
	172	printf("AliESDInputHandler::Notify() Path: %s\n", pathName.Data());
6ccd727b	173
1aa09fd1	174	if (fRunTag) {
	175	fRunTag->Clear();
	176	} else {
	177	fRunTag = new AliRunTag();
	178	}
	179
	180	delete fTreeT; fTreeT = 0;
6ccd727b	181
1aa09fd1	182	if (fChainT) {
	183	delete fChainT;
	184	fChainT = 0;
	185	}
	186
	187	if (!fChainT) {
	188	fChainT = new TChain("T");
	189	}
	190
	191
	192
	193	const char* tagPattern = "ESD.tag.root";
	194	const char* name = 0x0;
	195	TString tagFilename;
	196	if (zip) {
8ac5c1d7	197	TFile* file = fTree->GetCurrentFile();
1aa09fd1	198	TArchiveFile* arch = file->GetArchive();
	199	TObjArray* arr = arch->GetMembers();
	200	TIter next(arr);
	201
7cec0871	202	while ((file = (TFile*) next())) {
1aa09fd1	203	name = file->GetName();
1aa09fd1	204	if (strstr(name,tagPattern)) {
be4763e2	205	tagFilename = pathName.Data();
1aa09fd1	206	tagFilename += "#";
	207	tagFilename += name;
	208	fChainT->Add(tagFilename);
	209	AliInfo(Form("Adding %s to tag chain \n", tagFilename.Data()));
	210	}//pattern check
	211	} // archive file loop
	212	} else {
be4763e2	213	void * dirp = gSystem->OpenDirectory(pathName.Data());
1aa09fd1	214	while((name = gSystem->GetDirEntry(dirp))) {
1aa09fd1	215	if (strstr(name,tagPattern)) {
be4763e2	216	tagFilename = pathName.Data();
1aa09fd1	217	tagFilename += "/";
	218	tagFilename += name;
	219	fChainT->Add(tagFilename);
	220	AliInfo(Form("Adding %s to tag chain \n", tagFilename.Data()));
	221	}//pattern check
	222	}//directory loop
	223	}
	224	fChainT->SetBranchAddress("AliTAG",&fRunTag);
	225	fChainT->GetEntry(0);
	226	return kTRUE;
	227	}
	228
	229
c2b6979d	230
	231	Option_t *AliESDInputHandler::GetDataType() const
	232	{
	233	// Returns handled data type.
	234	return gESDDataType;
	235	}
f0f29480	236
	237	Int_t AliESDInputHandler::GetNEventAcceptedInFile()
	238	{
	239	// Get number of events in file accepted by the tag cuts
	240	// return -1 if no info is available
	241	if (!fTagCutSumm) {
	242	TList *luo = fTree->GetUserInfo();
	243	if (!luo) {
	244	AliInfo(Form("No user info in input tree - no tag cut summary\n"));
	245	return -1;
	246	}
	247	for (int iluo=0; iluo<luo->GetEntries(); iluo++) {
	248	fTagCutSumm = dynamic_cast<TMap *>(luo->At(iluo));
	249	if (fTagCutSumm) break;
	250	}
	251	if (!fTagCutSumm) {
	252	AliInfo(Form("No tag summary map in input tree\n"));
	253	return -1;
	254	}
	255	}
	256
	257	TObjString *ostr = 0;
	258	if (fTagCutSumm->FindObject(fTree->GetCurrentFile()->GetName()))
	259	ostr = (TObjString *) fTagCutSumm->GetValue(fTree->GetCurrentFile()->GetName());
	260	else {
	261	AliInfo(Form("No tag cut summary for file %s\n", fTree->GetCurrentFile()->GetName()));
	262	return -1;
	263	}
	264	char *iTagInfo;
	265	iTagInfo = strdup(ostr->GetString().Data());
	266
	267	Int_t iAcc = atoi(strtok(iTagInfo, ","));
	268
	269	AliInfo(Form("Got %i accepted events for file %s", iAcc, fTree->GetCurrentFile()->GetName()));
	270
	271	free(iTagInfo);
	272
	273	return iAcc;
	274	}
	275	Int_t AliESDInputHandler::GetNEventRejectedInFile()
	276	{
	277	// Get number of events in file rejected by the tag cuts
	278	// return -1 if no info is available
	279	if (!fTagCutSumm) {
	280	TList *luo = fTree->GetUserInfo();
	281	if (!luo) {
	282	AliInfo(Form("No user info in input tree - no tag cut summary\n"));
	283	return -1;
	284	}
	285	for (int iluo=0; iluo<luo->GetEntries(); iluo++) {
	286	fTagCutSumm = dynamic_cast<TMap *>(luo->At(iluo));
	287	if (fTagCutSumm) break;
	288	}
	289	if (!fTagCutSumm) {
	290	AliInfo(Form("No tag summary map in input tree\n"));
	291	return -1;
	292	}
	293	}
	294
	295	TObjString *ostr = 0;
	296	if (fTagCutSumm->FindObject(fTree->GetCurrentFile()->GetName()))
	297	ostr = (TObjString *) fTagCutSumm->GetValue(fTree->GetCurrentFile()->GetName());
	298	else {
	299	AliInfo(Form("No tag cut summary for file %s\n", fTree->GetCurrentFile()->GetName()));
300	return -1;
301	}
302	char *iTagInfo;
303	iTagInfo = strdup(ostr->GetString().Data());
304
305	strtok(iTagInfo, ",");
306	Int_t iRej = atoi(strtok(NULL, ","));
307
308	AliInfo(Form("Got %i accepted events for file %s", iRej, fTree->GetCurrentFile()->GetName()));
309
310	free(iTagInfo);
311
312	return iRej;
313	}
c334d0d3	314	Bool_t AliESDInputHandler::GetCutSummaryForChain(Int_t aTotal, Int_t aAccepted, Int_t *aRejected)
	315	{
	316	// Get number of events in the full chain
	317	// Count accepted and rejected events
	318	// return kFALSE if no info is available
	319	if (!fTagCutSumm) {
	320	TList *luo = fTree->GetUserInfo();
	321	if (!luo) {
	322	AliInfo(Form("No user info in input tree - no tag cut summary\n"));
	323	return kFALSE;
	324	}
	325	for (int iluo=0; iluo<luo->GetEntries(); iluo++) {
	326	fTagCutSumm = dynamic_cast<TMap *>(luo->At(iluo));
	327	if (fTagCutSumm) break;
	328	}
	329	if (!fTagCutSumm) {
	330	AliInfo(Form("No tag summary map in input tree\n"));
	331	return kFALSE;
	332	}
	333	}
	334
	335	TMapIter *tIter = new TMapIter(fTagCutSumm);
	336
	337	Int_t iTotList=0, iAccList=0, iRejList=0;
	338
	339	TObject *cobj;
9a587f77	340	while ((cobj = tIter->Next())) {
c334d0d3	341	TObjString kstr = (TObjString ) cobj;
	342	TObjString vstr = (TObjString ) fTagCutSumm->GetValue(kstr->GetString().Data());
	343	// printf("Got object value %s %s\n", kstr->GetString().Data(), vstr->GetString().Data());
	344	char *iTagInfo;
	345	iTagInfo = strdup(vstr->GetString().Data());
	346
	347	Int_t iAcc = atoi(strtok(iTagInfo, ","));
	348	Int_t iRej = atoi(strtok(NULL, ","));
	349
	350	iAccList += iAcc;
	351	iRejList += iRej;
	352	iTotList += (iAcc+iRej);
	353	}
	354
	355	*aTotal = iTotList;
	356	*aAccepted = iAccList;
	357	*aRejected = iRejList;
	358
	359	return kTRUE;
	360	}
	361
	362	Int_t AliESDInputHandler::GetNFilesEmpty()
	363	{
	364	// Count number of files in which all events were de-selected
	365	// For such files Notify() will NOT be called
	366	// return -1 if no info is available
	367	if (!fTagCutSumm) {
	368	TList *luo = fTree->GetUserInfo();
	369	if (!luo) {
	370	AliInfo(Form("No user info in input tree - no tag cut summary\n"));
	371	return -1;
	372	}
	373	for (int iluo=0; iluo<luo->GetEntries(); iluo++) {
	374	fTagCutSumm = dynamic_cast<TMap *>(luo->At(iluo));
	375	if (fTagCutSumm) break;
	376	}
	377	if (!fTagCutSumm) {
	378	AliInfo(Form("No tag summary map in input tree\n"));
	379	return -1;
	380	}
	381	}
	382
	383	TMapIter *tIter = new TMapIter(fTagCutSumm);
	384
	385	Int_t iFilesEmpty = 0;
	386
	387	TObject *cobj;
9a587f77	388	while ((cobj = tIter->Next())) {
c334d0d3	389	TObjString kstr = (TObjString ) cobj;
	390	TObjString vstr = (TObjString ) fTagCutSumm->GetValue(kstr->GetString().Data());
	391	// printf("Got object value %s %s\n", kstr->GetString().Data(), vstr->GetString().Data());
	392	char *iTagInfo;
	393	iTagInfo = strdup(vstr->GetString().Data());
	394
	395	Int_t iAcc = atoi(strtok(iTagInfo, ","));
	396	Int_t iRej = atoi(strtok(NULL, ","));
	397
	398	if ((iAcc == 0) && ((iRej+iAcc)>0))
	399	iFilesEmpty++;
	400	}
	401
	402	return iFilesEmpty;
	403
	404	}