-// Author: Benjamin Hess 23/09/2008
+// Author: Benjamin Hess 29/01/2010
/*************************************************************************
- * Copyright (C) 2008, Alexandru Bercuci, Benjamin Hess. *
+ * Copyright (C) 2009-2010, Alexandru Bercuci, Benjamin Hess. *
* All rights reserved. *
*************************************************************************/
+
//////////////////////////////////////////////////////////////////////////
// //
// AliEveTRDTrackList //
// the list in the same way as for the TEveElementList). In general, //
// please use AddMacro(...) for this purpose. //
// Macros that are no longer needed can be removed from the list via //
-// RemoveSelectionMacros(...) or RemoveProcessMacros(...) respectively. //
-// This function takes an iterator of the list of entries that are to //
+// RemoveSelectedMacros(...).This function takes an iterator of the //
+// list of macros that are to be removed. //
// be removed. An entry looks like: //
-// "MacroName.C (Path: MacroPath)". This is the way, the information //
-// about a macro is stored in the AliEveTRDTrackList. If you have path //
-// and name of a macro, use MakeMacroEntry(...) to get the corresponding//
-// entry. The type of the macros is stored in a map. You can get the //
-// macro type via GetMacroType(...). //
+// The data for each macro consists of path, name, type and the command //
+// that will be used to apply the macro. This stuff is stored in a map //
+// which takes the macro name for the key and the above mentioned data //
+// in a TMacroData-object for the value. //
+// You can get the macro type via GetMacroType(...). //
// With ApplySTSelectionMacros(...) or ApplyProcessMacros(...) //
// respectively you can apply the macros to the track list via //
-// iterators (same style like for RemoveProcessMacros(...) (cf. above)).//
+// iterators (same style like for RemoveSelectedMacros(...)(cf. above)).//
// Selection macros (de-)select macros according to a selection rule //
// by setting the rnr-state of the tracks. //
// If multiple selection macros are applied, a track is selected, if //
// TH1* YourMacro(const AliTRDtrackV1*); //
// TH1* YourMacro(const AliTRDtrackV1*, const AliTRDtrackV1*); //
// //
-// The macros which take 2 tracks are applied to all pairs //
-// fullfilling the selection criteria. //
+// The macros which take 2 tracks are applied to all track pairs //
+// (whereby BOTH tracks of the pair have to be selected by the single //
+// track selection macros and have to be unequal, otherwise they will //
+// be skipped) that have been selected by ALL correlated tracks //
+// selection macros. The selection macros with 2 tracks do NOT affect //
+// process macros that process only a single track! //
//////////////////////////////////////////////////////////////////////////
// Uncomment to display debugging infos
//#define ALIEVETRDTRACKLIST_DEBUG
-#include "AliEveTRDTrackList.h"
-
-#include <AliTRDReconstructor.h>
#include <TFile.h>
#include <TFunction.h>
#include <TH1.h>
#include <TSystem.h>
#include <TTree.h>
#include <TTreeStream.h>
+#include <TMethodCall.h>
+
+#include <AliTRDReconstructor.h>
+
+#include <EveDet/AliEveTRDTrackList.h>
+#include <EveDet/AliEveTRDTrackListEditor.h>
+
+#include <../PWG1/TRD/AliTRDrecoTask.h>
+#include <../PWG1/TRD/AliTRDpwg1Helper.h>
ClassImp(AliEveTRDTrackList)
///////////////////////////////////////////////////////////
AliEveTRDTrackList::AliEveTRDTrackList(const Text_t* n, const Text_t* t, Bool_t doColor):
TEveElementList(n, t, doColor),
- fMacroList(0),
- fMacroSelList(0),
- fDataFromMacroList(0),
- fMacroTypes(0),
- fDataTree(0),
+ fEditor(0x0),
+ fDataFromMacroList(0x0),
+ fMacroList(0x0),
+ fDataTree(0x0),
fHistoDataSelected(0),
fMacroListSelected(0),
- fMacroSelListSelected(0),
fSelectedTab(1), // Standard tab: "Apply macros" (index 1)
fSelectedStyle(0)
{
SetChildClass(AliEveTRDTrack::Class());
// Allocate memory for the lists and declare them as owners of their contents
- fMacroList = new TList();
- fMacroList->TCollection::SetOwner(kTRUE);
- fMacroSelList = new TList();
- fMacroSelList->TCollection::SetOwner(kTRUE);
fDataFromMacroList = new TList();
fDataFromMacroList->TCollection::SetOwner(kTRUE);
- fMacroTypes = new TMap();
+ fMacroList = new TMap();
// Set map to owner of it's objects to delete them, if they are removed from the map
- fMacroTypes->SetOwnerKeyValue(kTRUE, kTRUE);
+ fMacroList->SetOwnerKeyValue(kTRUE, kTRUE);
// Set the build directory for AClic
- gSystem->SetBuildDir("$HOME/.trdQArec");
+ if(gSystem->AccessPathName(Form("%s/.trdQArec" , gSystem->Getenv("HOME")))) gSystem->Exec("mkdir $HOME/.trdQArec");
+ gSystem->SetBuildDir(Form("%s/.trdQArec", gSystem->Getenv("HOME")));
- AddStandardMacros();
+ AddStandardContent();
}
//______________________________________________________
{
// Frees allocated memory (lists etc.).
- if (fMacroList != 0)
+ // Let the editor know that the list will be destroyed -> The editor will save the data
+ if (fEditor != 0)
{
- fMacroList->Delete();
- delete fMacroList;
- fMacroList = 0;
+ fEditor->SaveMacroList(fMacroList);
+ fEditor = 0;
}
- if (fMacroSelList != 0)
- {
- fMacroSelList->Delete();
- delete fMacroSelList;
- fMacroSelList = 0;
- }
+
if (fDataFromMacroList != 0)
{
fDataFromMacroList->Delete();
delete fDataTree;
fDataTree = 0;
}
- if (fMacroTypes != 0)
+ if (fMacroList != 0)
{
- fMacroTypes->DeleteAll();
- delete fMacroTypes;
- fMacroTypes = 0;
+ fMacroList->DeleteAll();
+ delete fMacroList;
+ fMacroList = 0;
}
// Note: gSystem->AccessPathName(...) returns kTRUE, if the access FAILED!
if(!gSystem->AccessPathName(Form("/tmp/TRD.TrackListMacroData_%s.root", gSystem->Getenv("USER"))))
// void YourMacro(const AliTRDtrackV1*, const AliTRDtrackV1*,
// Double_t*&, Int_t&)
// TH1* YourMacro(const AliTRDtrackV1*)
- // TH1* YourMacro(const AliTRDtrackV1*, const AliTRDtrackV1*)
- //
- // The macros which take 2 tracks are applied to all pairs
- // fullfilling the selection criteria.
-
- Char_t* entryName = MakeMacroEntry(path, nameC);
+ // TH1* YourMacro(const AliTRDtrackV1*, const AliTRDtrackV1*)
Char_t pathname[fkMaxMacroPathNameLength];
memset(pathname, '\0', sizeof(Char_t) * fkMaxMacroPathNameLength);
// Expand the path and create the pathname
Char_t* systemPath = gSystem->ExpandPathName(path);
- sprintf(pathname, "%s/%s", systemPath, nameC);
+ snprintf(pathname, fkMaxMacroPathNameLength, "%s/%s", systemPath, nameC);
delete systemPath;
systemPath = 0;
for (UInt_t ind = 0; ind < fkMaxMacroNameLength && ind < strlen(nameC) - 2; ind++) name[ind] = nameC[ind];
// Check, if files exists
- FILE* fp = 0;
-
- fp = fopen(pathname, "rb");
- if (fp != 0)
- {
+ FILE* fp = 0x0;
+ if((fp = fopen(pathname, "rb"))){
fclose(fp);
- fp = 0;
- }
- else
- {
- if (entryName != 0) delete entryName;
- entryName = 0;
-
- return NOT_EXIST_ERROR;
- }
-
+ fp = 0x0;
+ } else return NOT_EXIST_ERROR;
+
// Clean up root, load the desired macro and then check the type of the macro
- gROOT->Reset();
+ // A.B. gROOT->Reset();
- if (forceReload) gROOT->ProcessLineSync(Form(".L %s++", pathname));
- else gROOT->ProcessLineSync(Form(".L %s+", pathname));
+ gROOT->ProcessLineSync(Form(".L %s+%c", pathname, forceReload ? '+' : ' '));
- AliEveTRDTrackListMacroType type = GetMacroType(entryName, kFALSE);
+ // We need this line... otherwise, in some cases, there will be problems concerning ACLIC
+ gROOT->ProcessLineSync(Form(".L %s", pathname));
+
+ AliEveTRDTrackListMacroType type = GetMacroType(name, kFALSE);
// Clean up again
- gROOT->Reset();
+ // A.B. gROOT->Reset();
// Has not the correct signature!
- if (type == kUnknown)
- {
- if (entryName != 0) delete entryName;
- entryName = 0;
- return SIGNATURE_ERROR;
- }
-
- Int_t returnValue = WARNING;
+ if (type == kUnknown) return SIGNATURE_ERROR;
// Only add macro, if it is not already in the list
- if ((type == kSingleTrackAnalyse || type == kSingleTrackHisto
- || type == kCorrelTrackAnalyse || type == kCorrelTrackHisto) && fMacroList->FindObject(entryName) == 0)
- {
- fMacroList->Add(new TObjString(entryName));
- fMacroList->Sort();
-
- fMacroTypes->Add(new TObjString(entryName), new TObjString(Form("%d", type)));
-
- // We do not know, where the element has been inserted - deselect this list
- fMacroListSelected = 0;
-
- returnValue = SUCCESS;
- }
- else if ((type == kSingleTrackSelect || type == kCorrelTrackSelect) && fMacroSelList->FindObject(entryName) == 0)
- {
- fMacroSelList->Add(new TObjString(entryName));
- fMacroSelList->Sort();
-
- fMacroTypes->Add(new TObjString(entryName), new TObjString(Form("%d", type)));
-
- // We do not know, where the element has been inserted - deselect this list
- fMacroSelListSelected = 0;
-
- returnValue = SUCCESS;
+ Int_t returnValue = WARNING;
+ if(fMacroList->GetValue(name) == 0) {
+ returnValue = AddMacroFast(path, name, type) ? SUCCESS : ERROR;
}
- else returnValue = WARNING;
-
- if (entryName != 0) delete entryName;
- entryName = 0;
-
return returnValue;
}
//______________________________________________________
-void AliEveTRDTrackList::AddMacroFast(const Char_t* entry, AliEveTRDTrackListMacroType type)
+Bool_t AliEveTRDTrackList::AddMacroFast(const Char_t* path, const Char_t* name, AliEveTRDTrackListMacroType type)
{
- // Adds an entry to the corresponding list (cf. overloaded function).
+ // Adds a macro (path/name) to the corresponding list. No checks are performed (file exist,
+ // macro already in list/map, signature correct), no libraries are created!
+ // You can use this function only, if the macro has been added successfully before
+ // (and then maybe was removed). The function is very fast. On success kTRUE is returned, otherwise: kFALSE;
+
+ Bool_t success = kFALSE;
switch (type)
{
case kSingleTrackSelect:
case kCorrelTrackSelect:
- fMacroSelList->Add(new TObjString(entry));
- fMacroSelList->Sort();
-
- fMacroTypes->Add(new TObjString(entry), new TObjString(Form("%d", type)));
-
- // We do not know, where the element has been inserted - deselect this list
- fMacroSelListSelected = 0;
-
- break;
case kSingleTrackAnalyse:
case kSingleTrackHisto:
case kCorrelTrackAnalyse:
case kCorrelTrackHisto:
- fMacroList->Add(new TObjString(entry));
- fMacroList->Sort();
-
- fMacroTypes->Add(new TObjString(entry), new TObjString(Form("%d", type)));
+ fMacroList->Add(new TObjString(name), new TMacroData(name, path, type));
// We do not know, where the element has been inserted - deselect this list
fMacroListSelected = 0;
+
+ success = kTRUE;
+
+#ifdef ALIEVETRDTRACKLIST_DEBUG
+ // Successfull add will only be displayed in debug mode
+ printf("AliEveTRDTrackList::AddMacroFast: Added macro \"%s/%s\" to the corresponding list\n", path, name);
+#endif
+
break;
- default:
- Error("AliEveTRDTrackList::AddMacroFast", Form("Unknown macro type for entry \"%s\"!", entry));
- break;
- }
-}
-//______________________________________________________
-void AliEveTRDTrackList::AddMacroFast(const Char_t* path, const Char_t* name, AliEveTRDTrackListMacroType type)
-{
- // Adds a macro (path/name) to the list associated with the "type" parameter.
- // No checks are performed (fast) and no libraries are loaded.
- // Do use only, if library already exists!
+ default:
+ // Error will always be displayed
+ printf("AliEveTRDTrackList::AddMacroFast: ERROR: Could not add macro \"%s/%s\" to the corresponding list\n",
+ path, name);
- Char_t* entry = MakeMacroEntry(path, name);
- if (entry != 0)
- {
- AddMacroFast(entry, type);
+ success = kFALSE;
-#ifdef ALIEVETRDTRACKLIST_DEBUG
- // Successfull add will only be displayed in debug mode
- printf("#AliEveTRDTrackList::AddMacroFast: Added macro \"%s/%s\" to the corresponding list\n", path, name);
-#endif
-
- delete entry;
- entry = 0;
- }
- else
- {
- // Error will always be displayed
- printf("#AliEveTRDTrackList::AddMacroFast: ERROR: Could not add macro \"%s/%s\" to the corresponding list\n",
- path, name);
+ break;
}
+
+ return success;
}
//______________________________________________________
-void AliEveTRDTrackList::AddStandardMacros()
+void AliEveTRDTrackList::AddStandardContent()
{
- // Adds standard macros to the lists.
+ // Adds standard macros to the macro list.
- // Add your standard macros here, e.g.:
+ // Add your standard macros here, e.g.:
// To add a macro use:
// AddMacro("$(ALICE_ROOT)/myFolder", "myMacroName.C");
// -> If the file does not exist, nothing happens. So if you want to handle this,
// use the return value of AddMacro (NOT_EXIST_ERROR is returned, if file does not exist)
// (-> You can also check for other return values (see AddMacro(...)))
- AddMacro("$(ALICE_ROOT)/TRD/qaRec/macros", "clusterSelection.C");
- AddMacro("$(ALICE_ROOT)/TRD/qaRec/macros", "chargeDistr.C");
- AddMacro("$(ALICE_ROOT)/TRD/qaRec/macros", "clusterResiduals.C");
- AddMacro("$(ALICE_ROOT)/TRD/qaRec/macros", "PH.C");
+
+ const Char_t *libs[] = {"libANALYSIS.so", "libANALYSISalice.so", "libTENDER.so", "libPWG1.so"};
+ Int_t nlibs = static_cast<Int_t>(sizeof(libs)/sizeof(Char_t *));
+ for(Int_t ilib=0; ilib<nlibs; ilib++){
+ if(gSystem->Load(libs[ilib]) >= 0) continue;
+ AliError(Form("Fail loading %s.", libs[ilib]));
+ return;
+ }
+
+ const Char_t *fgkTRDPWG1taskClassName[AliTRDpwg1Helper::kNTRDQATASKS] = {
+ "AliTRDcheckESD"
+ ,"AliTRDinfoGen"
+ ,"AliTRDcheckDET"
+ ,"AliTRDefficiency"
+ ,"AliTRDresolution"
+ ,"AliTRDcheckPID"
+ ,"AliTRDv0Monitor"
+ };
+ AliTRDrecoTask *task(NULL);
+ TList *fPlots(NULL);
+ for(Int_t it=2; it<AliTRDpwg1Helper::kNTRDQATASKS; it++){
+ TClass c(fgkTRDPWG1taskClassName[it]);
+ task = (AliTRDrecoTask*)c.New();
+ task->SetMCdata(kFALSE);
+ if(!(fPlots = task->GetPlotFunctors())){
+ //AliWarning(Form("No Plot functors defined for task \"%s\"", fgkTRDtaskClassName[it]));
+ delete task;
+ continue;
+ }
+ if(!(task->Histos())){
+ //AliWarning(Form("No Ref Histograms defined for task \"%s\"", fgkTRDtaskClassName[it]));
+ delete task;
+ continue;
+ }
+
+ // export task to CINT and add functions
+ gROOT->ProcessLine(Form("%s* %s = (%s*)%p;", fgkTRDPWG1taskClassName[it], task->GetName(), fgkTRDPWG1taskClassName[it], (void*)task));
+ TIter iter(fPlots); TMethodCall *m = 0x0;
+ while((m = dynamic_cast<TMethodCall*>(iter()))){
+ AddMacroFast("", Form("%s->%s", task->GetName(), m->GetMethodName()), kSingleTrackHisto);
+ }
+ }
}
+
//______________________________________________________
Bool_t AliEveTRDTrackList::ApplyProcessMacros(const TList* selIterator, const TList* procIterator)
{
if (procIterator->GetEntries() <= 0) return kTRUE;
// Clear root
- gROOT->Reset();
+ // A.B. gROOT->Reset();
// Clear old data and re-allocate
- if (fDataTree == 0) fDataTree = new TTreeSRedirector(Form("/tmp/TRD.TrackListMacroData_%s.root",
- gSystem->Getenv("USER")));
- if (!fDataTree)
- {
- Error("Apply process macros", Form("File \"/tmp/TRD.TrackListMacroData_%s.root\" could not be accessed properly!",
- gSystem->Getenv("USER")));
+ if (!fDataTree){
+ TDirectory *cwd = gDirectory;
+ fDataTree = new TTreeSRedirector(Form("/tmp/TRD.TrackListMacroData_%s.root", gSystem->Getenv("USER")));
+ cwd->cd();
+ }
+ if (!fDataTree){
+ Error("Apply process macros", Form("File \"/tmp/TRD.TrackListMacroData_%s.root\" could not be accessed properly!", gSystem->Getenv("USER")));
return kFALSE;
}
- if (fDataFromMacroList != 0)
- {
+ if (fDataFromMacroList != 0) {
fDataFromMacroList->Delete();
delete fDataFromMacroList;
}
fHistoDataSelected = 0;
- Char_t name[fkMaxMacroNameLength];
- Char_t** procCmds = new Char_t*[procIterator->GetEntries()];
- Char_t** selCmds = new Char_t*[selIterator->GetEntries()];
+ TMacroData* macro = 0;
+
+ TString* procCmds = new TString[procIterator->GetEntries()];
AliEveTRDTrackListMacroType* mProcType = new AliEveTRDTrackListMacroType[procIterator->GetEntries()];
- AliEveTRDTrackListMacroType* mSelType = new AliEveTRDTrackListMacroType[selIterator->GetEntries()];
+ TString* selCmds(NULL);
+ AliEveTRDTrackListMacroType* mSelType(NULL);
+ if (selIterator->GetEntries() > 0) {
+ selCmds = new TString[selIterator->GetEntries()];
+ mSelType = new AliEveTRDTrackListMacroType[selIterator->GetEntries()];
+ }
+
Bool_t selectedByCorrSelMacro = kFALSE;
AliEveTRDTrackListMacroType macroType = kUnknown;
Int_t numHistoMacros = 0;
- TH1** histos = 0;
+ TH1** histos(NULL);
- AliEveTRDTrack* track1 = 0;
- AliEveTRDTrack* track2 = 0;
- TH1* returnedHist = 0x0;
+ AliEveTRDTrack* track1(NULL);
+ AliEveTRDTrack* track2(NULL);
// Collect the commands for each process macro and add them to "data-from-list"
- for (Int_t i = 0; i < procIterator->GetEntries(); i++)
- {
- memset(name, '\0', sizeof(Char_t) * fkMaxMacroNameLength);
-
- procCmds[i] = new Char_t[(fkMaxMacroPathNameLength + fkMaxApplyCommandLength)];
- memset(procCmds[i], '\0', sizeof(Char_t) * (fkMaxMacroNameLength + fkMaxApplyCommandLength));
+ for (Int_t i = 0; i < procIterator->GetEntries(); i++){
+ macro = (TMacroData*)fMacroList->GetValue(procIterator->At(i)->GetTitle());
-#ifdef ALIEVETRDTRACKLIST_DEBUG
- printf("AliEveTRDTrackList: Applying process macro: %s\n", procIterator->At(i)->GetTitle());
-#endif
-
- // Extract the name
- sscanf(procIterator->At(i)->GetTitle(), "%s (Path: %*s)", name);
-
- // Delete ".C" at the end
- // -> Note: Physical address pointer, do NOT delete. / Changes "name" as well!
- Char_t* dotC = (Char_t*)strrchr(name, '.');
- if (dotC != 0)
- {
- *dotC = '\0';
- dotC++;
- *dotC = '\0';
- }
-
- // Find the type of the process macro
- macroType = GetMacroType(procIterator->At(i)->GetTitle(), kTRUE);
- if (macroType == kSingleTrackHisto)
- {
- mProcType[i] = macroType;
- numHistoMacros++;
- // Create the command
- sprintf(procCmds[i], "%s(automaticTrackV1_1);", name);
-
- // Add to "data-from-list" -> Mark as a histo macro with the substring "(histo macro)"
- fDataFromMacroList->Add(new TObjString(Form("%s (histo macro)", name)));
+ if (!macro){
+ Error("Apply process macros",
+ Form("Macro list is corrupted: Macro \"%s\" is not registered!",
+ procIterator->At(i)->GetTitle()));
+ continue;
}
- else if (macroType == kSingleTrackAnalyse)
- {
- mProcType[i] = macroType;
- // Create the command
- sprintf(procCmds[i], "%s(automaticTrackV1_1, results, n);", name);
- // Add to "data-from-list"
- fDataFromMacroList->Add(new TObjString(name));
- }
- else if (macroType == kCorrelTrackHisto)
- {
+#ifdef ALIEVETRDTRACKLIST_DEBUG
+ printf("AliEveTRDTrackList: Checking process macro: %s\n", macro->GetName());
+#endif
+
+ // Find the type of the process macro
+ macroType = macro->GetType();
+ if (macroType == kSingleTrackHisto || macroType == kCorrelTrackHisto){
mProcType[i] = macroType;
numHistoMacros++;
// Create the command
- sprintf(procCmds[i], "%s(automaticTrackV1_1, automaticTrackV1_2);", name);
+ procCmds[i] = macro->GetCmd();
// Add to "data-from-list" -> Mark as a histo macro with the substring "(histo macro)"
- fDataFromMacroList->Add(new TObjString(Form("%s (histo macro)", name)));
- }
- else if (macroType == kCorrelTrackAnalyse)
- {
+ fDataFromMacroList->Add(new TObjString(Form("%s (histo macro)", macro->GetName())));
+ } else if (macroType == kSingleTrackAnalyse || macroType == kCorrelTrackAnalyse) {
mProcType[i] = macroType;
// Create the command
- sprintf(procCmds[i], "%s(automaticTrackV1_1, automaticTrackV1_2, results, n);", name);
+ procCmds[i] = macro->GetCmd();
// Add to "data-from-list"
- fDataFromMacroList->Add(new TObjString(name));
- }
- else
- {
+ fDataFromMacroList->Add(new TObjString(macro->GetName()));
+ } else {
Error("Apply process macros",
- Form("Process macro list corrupted: Macro \"%s\" is not registered as a process macro!", name));
+ Form("Macro list corrupted: Macro \"%s/%s.C\" is not registered as a process macro!",
+ macro->GetPath(), macro->GetName()));
mProcType[i] = kUnknown;
}
}
// Collect the commands for each selection macro and add them to "data-from-list"
- for (Int_t i = 0; i < selIterator->GetEntries(); i++)
- {
- memset(name, '\0', sizeof(Char_t) * fkMaxMacroNameLength);
-
- selCmds[i] = new Char_t[(fkMaxMacroPathNameLength + fkMaxApplyCommandLength)];
- memset(selCmds[i], '\0', sizeof(Char_t) * (fkMaxMacroNameLength + fkMaxApplyCommandLength));
+ for (Int_t i = 0; i < selIterator->GetEntries(); i++){
+ macro = (TMacroData*)fMacroList->GetValue(selIterator->At(i)->GetTitle());
+
+ if (!macro){
+ Error("Apply process macros",
+ Form("Macro list is corrupted: Macro \"%s\" is not registered!",
+ selIterator->At(i)->GetTitle()));
+ continue;
+ }
#ifdef ALIEVETRDTRACKLIST_DEBUG
- printf("AliEveTRDTrackList: Applying selection macro (correlated tracks): %s\n", selIterator->At(i)->GetTitle());
+ printf("AliEveTRDTrackList: Checking selection macro: %s\n", macro->GetName());
#endif
-
- // Extract the name
- sscanf(selIterator->At(i)->GetTitle(), "%s (Path: %*s)", name);
-
- // Delete ".C" at the end
- // -> Note: Physical address pointer, do NOT delete. / Changes "name" as well!
- Char_t* dotC = (Char_t*)strrchr(name, '.');
- if (dotC != 0)
- {
- *dotC = '\0';
- dotC++;
- *dotC = '\0';
- }
// Find the type of the process macro
- macroType = GetMacroType(selIterator->At(i)->GetTitle(), kTRUE);
+ macroType = macro->GetType();
+
// Single track select macro
- if (macroType == kSingleTrackSelect)
- {
+ if (macroType == kSingleTrackSelect) {
// Has already been processed by ApplySTSelectionMacros(...)
- mSelType[i] = macroType;
+ if(mSelType) mSelType[i] = macroType;
}
// Correlated tracks select macro
- else if (macroType == kCorrelTrackSelect)
- {
- mSelType[i] = macroType;
+ else if (macroType == kCorrelTrackSelect) {
+ if(mSelType) mSelType[i] = macroType;
// Create the command
- sprintf(selCmds[i], "%s(automaticTrackV1_1, automaticTrackV1_2);", name);
- }
- else
- {
+ if(selCmds) selCmds[i] = macro->GetCmd();
+ } else {
Error("Apply process macros",
- Form("Selection macro list corrupted: Macro \"%s\" is not registered as a selection macro!", name));
- mProcType[i] = kUnknown;
+ Form("Macro list corrupted: Macro \"%s/%s.C\" is not registered as a selection macro!",
+ macro->GetPath(), macro->GetName()));
+ if(mSelType) mSelType[i] = kUnknown;
}
}
// Allocate memory for the histograms
- if (numHistoMacros > 0) histos = new TH1*[numHistoMacros];
- for (Int_t i = 0; i < numHistoMacros; i++) histos[i] = 0;
-
- // Walk through the list of tracks
- for (TEveElement::List_i iter = this->BeginChildren(); iter != this->EndChildren(); ++iter)
- {
- track1 = dynamic_cast<AliEveTRDTrack*>(*iter);
+ if (numHistoMacros > 0){
+ histos = new TH1*[numHistoMacros];
+ memset(histos, 0, numHistoMacros*sizeof(TH1*));
+ }
+
+ //////////////////////////////////
+ // WALK THROUGH THE LIST OF TRACKS
+ //////////////////////////////////
+ for (TEveElement::List_i iter = this->BeginChildren(); iter != this->EndChildren(); ++iter){
+ if(!(track1 = dynamic_cast<AliEveTRDTrack*>(*iter))) continue;
- if (!track1) continue;
-
// Skip tracks that have not been selected
if (!track1->GetRnrState()) continue;
-
- track1->ExportToCINT((Text_t*)"automaticTrack");
+
// Cast to AliTRDtrackV1
+ gROOT->ProcessLineSync(Form("AliEveTRDTrack *automaticTrack = (AliEveTRDTrack*)%p;", (void*)track1));
gROOT->ProcessLineSync("AliTRDtrackV1* automaticTrackV1_1 = (AliTRDtrackV1*)automaticTrack->GetUserData();");
// Collect data for each macro
- for (Int_t i = 0, histoIndex = 0; i < procIterator->GetEntries(); i++)
- {
+ for (Int_t i = 0, histoIndex = 0; i < procIterator->GetEntries(); i++){
// Single track histo
- if (mProcType[i] == kSingleTrackHisto)
- {
- returnedHist = (TH1*)gROOT->ProcessLineSync(procCmds[i]);
- if (returnedHist != 0x0)
- {
- if (histos[histoIndex] == 0) histos[histoIndex] = returnedHist;
- else
- {
- histos[histoIndex]->Add((const TH1*)returnedHist);
- delete returnedHist;
- returnedHist = 0;
- }
- }
- histoIndex++;
- }
- // Correlated tracks histo
- else if (mProcType[i] == kCorrelTrackHisto)
- {
+ if (mProcType[i] == kSingleTrackHisto){
+ if(histos) histos[histoIndex++] = (TH1*)gROOT->ProcessLineSync(procCmds[i]);
+ // Correlated tracks histo
+ } else if (mProcType[i] == kCorrelTrackHisto) {
// Loop over all pairs behind the current one - together with the other loop this will be a loop
// over all pairs. We have a pair of tracks, if and only if both tracks of the pair are selected (Rnr-state)
// and are not equal.
- // The correlated tracks process macro will applied to all pairs that will be additionally selected by
+ // The correlated tracks process macro will be applied to all pairs that will be additionally selected by
// all correlated tracks selection macros.
TEveElement::List_i iter2 = iter;
iter2++;
for ( ; iter2 != this->EndChildren(); ++iter2)
{
- track2 = dynamic_cast<AliEveTRDTrack*>(*iter2);
+ if(!(track2 = dynamic_cast<AliEveTRDTrack*>(*iter2))) continue;
- if (!track2) continue;
-
// Skip tracks that have not been selected
if (!track2->GetRnrState()) continue;
- track2->ExportToCINT((Text_t*)"automaticTrack");
// Cast to AliTRDtrackV1
+ gROOT->ProcessLineSync(Form("AliEveTRDTrack *automaticTrack = (AliEveTRDTrack*)%p;", (void*)track2));
gROOT->ProcessLineSync("AliTRDtrackV1* automaticTrackV1_2 = (AliTRDtrackV1*)automaticTrack->GetUserData();");
// Select track by default (so it will be processed, if there are no correlated tracks selection macros!)
selectedByCorrSelMacro = kTRUE;
- for (Int_t j = 0; j < selIterator->GetEntries(); j++)
- {
- if (mSelType[j] == kCorrelTrackSelect)
- {
+ for (Int_t j = 0; j < selIterator->GetEntries(); j++){
+ if (mSelType && mSelType[j] == kCorrelTrackSelect){
selectedByCorrSelMacro = (Bool_t)gROOT->ProcessLineSync(selCmds[j]);
if (!selectedByCorrSelMacro) break;
}
// If the pair has not been selected by the correlated tracks selection macros, skip it!
if (!selectedByCorrSelMacro) continue;
- returnedHist = (TH1*)gROOT->ProcessLineSync(procCmds[i]);
- if (returnedHist != 0x0)
- {
- if (histos[histoIndex] == 0) histos[histoIndex] = returnedHist;
- else
- {
- histos[histoIndex]->Add((const TH1*)returnedHist);
-
- delete returnedHist;
- returnedHist = 0;
- }
- }
+ if(histos) histos[histoIndex] = (TH1*)gROOT->ProcessLineSync(procCmds[i]);
}
histoIndex++;
}
// Single track analyse
- else if (mProcType[i] == kSingleTrackAnalyse)
- {
+ else if (mProcType[i] == kSingleTrackAnalyse) {
// Create data pointers in CINT, execute the macro and get the data
gROOT->ProcessLineSync("Double_t* results = 0;");
gROOT->ProcessLineSync("Int_t n = 0;");
Double_t* results = (Double_t*)gROOT->ProcessLineSync("results;");
Int_t nResults = (Int_t)gROOT->ProcessLineSync("n;");
- if (results == 0)
- {
+ if (results == 0) {
Error("Apply macros", Form("Error reading data from macro \"%s\"", procIterator->At(i)->GetTitle()));
continue;
}
- for (Int_t resInd = 0; resInd < nResults; resInd++)
- {
+ for (Int_t resInd = 0; resInd < nResults; resInd++){
(*fDataTree) << Form("TrackData%d", i) << Form("Macro%d=", i) << results[resInd] << (Char_t*)"\n";
}
results = 0;
}
// Correlated tracks analyse
- else if (mProcType[i] == kCorrelTrackAnalyse)
- {
+ else if (mProcType[i] == kCorrelTrackAnalyse){
// Loop over all pairs behind the current one - together with the other loop this will be a loop
// over all pairs. We have a pair of tracks, if and only if both tracks of the pair are selected (Rnr-state)
// and are not equal.
- // The correlated tracks process macro will applied to all pairs that will be additionally selected by
+ // The correlated tracks process macro will be applied to all pairs that will be additionally selected by
// all correlated tracks selection macros.
TEveElement::List_i iter2 = iter;
iter2++;
- for ( ; iter2 != this->EndChildren(); ++iter2)
- {
- track2 = dynamic_cast<AliEveTRDTrack*>(*iter2);
- if (!track2) continue;
-
+ for ( ; iter2 != this->EndChildren(); ++iter2) {
+ if(!(track2 = dynamic_cast<AliEveTRDTrack*>(*iter2))) continue;
+
// Skip tracks that have not been selected
if (!track2->GetRnrState()) continue;
-
- track2->ExportToCINT((Text_t*)"automaticTrack");
+
// Cast to AliTRDtrackV1
+ gROOT->ProcessLineSync(Form("AliEveTRDTrack *automaticTrack = (AliEveTRDTrack*)%p;", (void*)track2));
gROOT->ProcessLineSync("AliTRDtrackV1* automaticTrackV1_2 = (AliTRDtrackV1*)automaticTrack->GetUserData();");
// Select track by default (so it will be processed, if there are no correlated tracks selection macros!)
selectedByCorrSelMacro = kTRUE;
- for (Int_t j = 0; j < selIterator->GetEntries(); j++)
- {
- if (mSelType[j] == kCorrelTrackSelect)
- {
+ for (Int_t j = 0; j < selIterator->GetEntries(); j++) {
+ if (mSelType && mSelType[j] == kCorrelTrackSelect) {
selectedByCorrSelMacro = (Bool_t)gROOT->ProcessLineSync(selCmds[j]);
if (!selectedByCorrSelMacro) break;
}
gROOT->ProcessLineSync(procCmds[i]);
Double_t* results = (Double_t*)gROOT->ProcessLineSync("results;");
Int_t nResults = (Int_t)gROOT->ProcessLineSync("n;");
-
- if (results == 0)
- {
+
+ if (results == 0) {
Error("Apply macros", Form("Error reading data from macro \"%s\"", procIterator->At(i)->GetTitle()));
continue;
}
- for (Int_t resInd = 0; resInd < nResults; resInd++)
- {
+ for (Int_t resInd = 0; resInd < nResults; resInd++) {
(*fDataTree) << Form("TrackData%d", i) << Form("Macro%d=", i) << results[resInd] << (Char_t*)"\n";
}
}
}
- for (Int_t i = 0, histoIndex = 0; i < procIterator->GetEntries() && histoIndex < numHistoMacros; i++)
- {
- if (mProcType[i] == kSingleTrackHisto || mProcType[i] == kCorrelTrackHisto)
- {
+ for (Int_t i = 0, histoIndex = 0; i < procIterator->GetEntries() && histoIndex < numHistoMacros; i++) {
+ if (mProcType[i] == kSingleTrackHisto || mProcType[i] == kCorrelTrackHisto) {
// Might be empty (e.g. no tracks have been selected)!
- if (histos[histoIndex] != 0)
- {
+ if (histos[histoIndex]) {
(*fDataTree) << Form("TrackData%d", i) << Form("Macro%d=", i) << histos[histoIndex] << (Char_t*)"\n";
}
histoIndex++;
}
}
- if (fDataTree != 0) delete fDataTree;
- fDataTree = 0;
+ if (fDataTree) delete fDataTree;
+ fDataTree = NULL;
- if (procCmds != 0) delete [] procCmds;
- if (mProcType != 0) delete mProcType;
- mProcType = 0;
+ if (procCmds) delete [] procCmds;
+ procCmds = NULL;
+ if (mProcType) delete [] mProcType;
+ mProcType = NULL;
- if (selCmds != 0) delete [] selCmds;
- if (mSelType != 0) delete mSelType;
- mSelType = 0;
+ if (selCmds) delete [] selCmds;
+ selCmds = NULL;
+ if (mSelType) delete [] mSelType;
+ mSelType = NULL;
- if (histos != 0) delete [] histos;
- histos = 0;
+ if (histos) delete [] histos;
+ histos = NULL;
// Clear root
- gROOT->Reset();
+ // A.B. gROOT->Reset();
// If there is data, select the first data set
if (procIterator->GetEntries() > 0) SETBIT(fHistoDataSelected, 0);
// "ST" stands for "single track". This means that only single track selection macros are applied.
// Correlated tracks selection macros will be used inside the call of ApplyProcessMacros(...)!
- Char_t name[fkMaxMacroNameLength];
- Char_t cmd[(fkMaxMacroNameLength + fkMaxApplyCommandLength)];
-
+ TMacroData* macro = 0;
AliEveTRDTrackListMacroType macroType = kUnknown;
AliEveTRDTrack* track1 = 0;
Bool_t selectedByMacro = kFALSE;
// Clear root
- gROOT->Reset();
+ // A.B. gROOT->Reset();
- // Select all tracks at first. A track is then deselect, if at least one selection macro
- // returns kFALSE for this track
- // Enable all tracks (Note: EnableListElements(..) will call "ElementChanged", which will cause unforeseen behavior!)
- for (TEveElement::List_i iter = this->BeginChildren(); iter != this->EndChildren(); ++iter)
- {
- ((TEveElement*)(*iter))->SetRnrState(kTRUE);
- }
+ // Select all tracks at first. A track is then deselected, if at least one selection macro
+ // returns kFALSE for this track.
+ // Enable all tracks (Note: EnableListElements(..) will call "ElementChanged", which will cause unforeseen behaviour!)
+ for (TEveElement::List_i iter = this->BeginChildren(); iter != this->EndChildren(); ++iter) ((TEveElement*)(*iter))->SetRnrState(kTRUE);
SetRnrState(kTRUE);
- for (Int_t i = 0; i < iterator->GetEntries(); i++)
- {
+ for (Int_t i = 0; i < iterator->GetEntries(); i++){
+ macro = (TMacroData*)fMacroList->GetValue(iterator->At(i)->GetTitle());
- memset(name, '\0', sizeof(Char_t) * fkMaxMacroNameLength);
- memset(cmd, '\0', sizeof(Char_t) * (fkMaxMacroNameLength + fkMaxApplyCommandLength));
+ if (!macro){
+ Error("Apply selection macros",
+ Form("Macro list is corrupted: Macro \"%s\" is not registered!", iterator->At(i)->GetTitle()));
+ continue;
+ }
#ifdef ALIEVETRDTRACKLIST_DEBUG
- printf("AliEveTRDTrackList: Applying selection macro: %s\n", iterator->At(i)->GetTitle());
+ printf("AliEveTRDTrackList: Applying selection macro: %s\n", macro->GetName());
#endif
- // Extract the name
- sscanf(iterator->At(i)->GetTitle(), "%s (Path: %*s)", name);
- // Delete ".C" at the end
- // -> Note: Physical address pointer, do NOT delete. / Changes "name" as well!
- Char_t* dotC = (Char_t*)strrchr(name, '.');
- if (dotC != 0)
- {
- *dotC = '\0';
- dotC++;
- *dotC = '\0';
- }
-
// Determine macro type
- macroType = GetMacroType(iterator->At(i)->GetTitle(), kTRUE);
+ macroType = macro->GetType();
// Single track select macro
- if (macroType == kSingleTrackSelect)
- {
- // Create the command
- sprintf(cmd, "%s(automaticTrackV1);", name);
-
+ if (macroType == kSingleTrackSelect){
// Walk through the list of tracks
for (TEveElement::List_i iter = this->BeginChildren(); iter != this->EndChildren(); ++iter)
{
// If the track has already been deselected, nothing is to do here
if (!track1->GetRnrState()) continue;
- track1->ExportToCINT((Text_t*)"automaticTrack");
// Cast to AliTRDtrackV1
- gROOT->ProcessLineSync("AliTRDtrackV1* automaticTrackV1 = (AliTRDtrackV1*)automaticTrack->GetUserData();");
- selectedByMacro = (Bool_t)gROOT->ProcessLineSync(cmd);
+ gROOT->ProcessLineSync(Form("AliEveTRDTrack *automaticTrack = (AliEveTRDTrack*)%p;", (void*)track1));
+ gROOT->ProcessLineSync("AliTRDtrackV1* automaticTrackV1_1 = (AliTRDtrackV1*)automaticTrack->GetUserData();");
+ selectedByMacro = (Bool_t)gROOT->ProcessLineSync(macro->GetCmd());
track1->SetRnrState(selectedByMacro && track1->GetRnrState());
}
}
// Correlated tracks select macro
- else if (macroType == kCorrelTrackSelect)
- {
+ else if (macroType == kCorrelTrackSelect){
// Will be processed in ApplyProcessMacros(...)
continue;
- }
- else
- {
+ } else {
Error("Apply selection macros",
- Form("Selection macro list corrupted: Macro \"%s\" is not registered as a selection macro!", name));
+ Form("Macro list corrupted: Macro \"%s/%s.C\" is not registered as a selection macro!",
+ macro->GetPath(), macro->GetName()));
}
}
// Clear root
- gROOT->Reset();
+ // A.B. gROOT->Reset();
}
//______________________________________________________
-AliEveTRDTrackList::AliEveTRDTrackListMacroType AliEveTRDTrackList::GetMacroType(const Char_t* entry, Bool_t UseList) const
+AliEveTRDTrackList::AliEveTRDTrackListMacroType AliEveTRDTrackList::GetMacroType(const Char_t* name, Bool_t UseList) const
{
- // Returns the type of the macro of the corresponding entry (i.e. "macro.C (Path: path)").
- // If you have only the name and the path, you can simply use MakeMacroEntry.
+ // Returns the type of the corresponding macro.
// If "UseList" is kTRUE, the type will be looked up in the internal list (very fast). But if this list
// does not exist, you have to use kFALSE for this parameter. Then the type will be determined by the
// prototype! NOTE: It is assumed that the macro has been compiled! If not, the return value is not
// predictable, but normally will be kUnknown.
- // Note: AddMacro(Fast) will update the internal list and RemoveProcess(/Selection)Macros respectively.
+ // Note: AddMacro(Fast) will update the internal list and RemoveMacros respectively.
AliEveTRDTrackListMacroType type = kUnknown;
- // Re do the check of the macro type
- if (!UseList)
- {
- Char_t name[fkMaxMacroNameLength];
-
- memset(name, '\0', sizeof(Char_t) * fkMaxMacroNameLength);
-
- // Extract the name
- sscanf(entry, "%s (Path: %*s)", name);
-
- // Delete ".C" at the end
- // -> Note: Physical address pointer, do NOT delete. / Changes "name" as well!
- Char_t* dotC = (Char_t*)strrchr(name, '.');
- if (dotC != 0)
- {
- *dotC = '\0';
- dotC++;
- *dotC = '\0';
- }
-
+ // Re-do the check of the macro type
+ if (!UseList){
// Single track select macro or single track histo macro?
TFunction* f = gROOT->GetGlobalFunctionWithPrototype(name, "const AliTRDtrackV1*", kTRUE);
if (f != 0x0)
if (strstr(f->GetMangledName(), "oPconstsPAliTRDtrackV1mUsP") != 0x0 &&
strstr(f->GetMangledName(), "cOconstsPAliTRDtrackV1mUsP") != 0x0)
{
- // Single track select macro?
+ // Correlated track select macro?
if (!strcmp(f->GetReturnTypeName(), "Bool_t"))
{
type = kCorrelTrackSelect;
}
- // single track histo macro?
+ // Correlated track histo macro?
else if (!strcmp(f->GetReturnTypeName(), "TH1*"))
{
type = kCorrelTrackHisto;
// Use list to look up the macro type
else
{
- TObjString* objEntry = 0;
- objEntry = (TObjString*)fMacroTypes->GetValue(entry);
- if (objEntry == 0) return kUnknown;
+ TMacroData* macro = 0;
+ macro = (TMacroData*)fMacroList->GetValue(name);
+ if (macro == 0) return kUnknown;
- type = (AliEveTRDTrackListMacroType)objEntry->GetString().Atoi();
+ type = macro->GetType();
switch (type)
{
case kSingleTrackSelect:
}
//______________________________________________________
-Char_t* AliEveTRDTrackList::MakeMacroEntry(const Char_t* path, const Char_t* name) const
-{
- // Constructs an entry for the macro lists with path and name.
-
- Char_t* entry = new Char_t[(fkMaxMacroPathNameLength + 30)];
- memset(entry, '\0', sizeof(Char_t) * (fkMaxMacroPathNameLength + 30));
-
- Char_t* systemPath = gSystem->ExpandPathName(path);
- sprintf(entry, "%s (Path: %s)", name, systemPath);
- delete systemPath;
- systemPath = 0;
-
- return entry;
-}
-
-//______________________________________________________
-void AliEveTRDTrackList::RemoveProcessMacros(const TList* iterator)
+void AliEveTRDTrackList::RemoveSelectedMacros(const TList* iterator)
{
- // Uses the iterator (for the selected process macros) to remove the process macros from
+ // Uses the iterator (for the selected macros) to remove the selected macros from
// the corresponding list.
- TObjString* obj = 0;
+ TObject* key = 0;
+ TPair* entry = 0;
for (Int_t i = 0; i < iterator->GetEntries(); i++)
{
- // Key and value will be deleted, too, since fMacroTypes is the owner of them
- fMacroTypes->DeleteEntry(fMacroTypes->FindObject(iterator->At(i)->GetTitle()));
+ entry = (TPair*)fMacroList->FindObject(iterator->At(i)->GetTitle());
- obj = (TObjString*)fMacroList->Remove(fMacroList->FindObject(iterator->At(i)->GetTitle()));
- if (obj != 0) delete obj;
- }
- obj = 0;
-}
+ if (entry == 0)
+ {
+ Error("AliEveTRDTrackList::RemoveSelectedMacros", Form("Macro \"%s\" not found in list!",
+ iterator->At(i)->GetTitle()));
+ continue;
+ }
+ key = entry->Key();
-//______________________________________________________
-void AliEveTRDTrackList::RemoveSelectionMacros(const TList* iterator)
-{
- // Uses the iterator (for the selected selection macros) to remove the selection macros from
- // the corresponding list.
-
- TObjString* obj = 0;
- for (Int_t i = 0; i < iterator->GetEntries(); i++)
- {
- // Key and value will be deleted, too, since fMacroTypes is the owner of them
- fMacroTypes->DeleteEntry(fMacroTypes->FindObject(iterator->At(i)->GetTitle()));
+ if (key == 0)
+ {
+ Error("AliEveTRDTrackList::RemoveSelectedMacros", Form("Key for macro \"%s\" not found in list!",
+ iterator->At(i)->GetTitle()));
+ continue;
+ }
+
+ // Key and value will be deleted, too, since fMacroList is the owner of them
+ Bool_t rem = fMacroList->DeleteEntry(key);
- obj = (TObjString*)fMacroSelList->Remove(fMacroSelList->FindObject(iterator->At(i)->GetTitle()));
- if (obj != 0) delete obj;
+ if (rem)
+ {
+#ifdef ALIEVETRDTRACKLIST_DEBUG
+ printf("AliEveTRDTrackList::RemoveSelectedMacros(): Removed macro: %s\n", iterator->At(i)->GetTitle());
+#endif
+ }
+ else
+ {
+ Error("AliEveTRDTrackList::RemoveSelectedMacros", Form("Macro \"%s\" could not be removed from the list!",
+ iterator->At(i)->GetTitle()));
+ }
}
- obj = 0;
}
//______________________________________________________
CLRBIT(fSelectedStyle, AliEveTRDTrack::kSource);
switch(ss)
{
- case AliTRDReconstructor::kLQPID:
+ case AliTRDpidUtil::kLQ:
CLRBIT(fSelectedStyle, AliEveTRDTrack::kPID);
break;
- case AliTRDReconstructor::kNNPID:
+ case AliTRDpidUtil::kNN:
SETBIT(fSelectedStyle, AliEveTRDTrack::kPID);
break;
}
CLRBIT(fSelectedStyle, AliEveTRDTrack::kTrackModel);
break;
case AliEveTRDTrack::kKalman:
- AliWarning("Kalman fit under testing for the moment.");
- //SETBIT(fSelectedStyle, AliEveTRDTrack::kTrackModel);
+ //AliWarning("Kalman fit under testing for the moment.");
+ SETBIT(fSelectedStyle, AliEveTRDTrack::kTrackModel);
break;
}
break;