#include "AliMonitorITS.h"
#include "AliMonitorV0s.h"
#include "AliMonitorHLT.h"
+#include "AliMonitorHLTHough.h"
#include "AliRawReaderRoot.h"
#include "AliLoader.h"
#include "AliRun.h"
#include "AliITStrackerV2.h"
#include "AliITSLoader.h"
#include "AliV0vertexer.h"
-#include <TSystem.h>
+#include <TGrid.h>
#include <TSocket.h>
+#include <TServerSocket.h>
#include <TMessage.h>
#include <TGridResult.h>
#include <TROOT.h>
#ifdef ALI_HLT
+#include <AliLevel3.h>
#include <AliL3Transform.h>
+#include <AliL3Track.h>
+#include <AliL3TrackArray.h>
+#include <AliL3StandardIncludes.h>
+#include <AliL3HoughMaxFinder.h>
+#include <AliL3HoughBaseTransformer.h>
+#include <AliL3Hough.h>
+#include <AliL3ClusterFitter.h>
+#include <AliL3Vertex.h>
+#include <AliL3Fitter.h>
+#include <AliL3DDLDataFileHandler.h>
#endif
ClassImp(AliMonitorProcess)
-const Int_t AliMonitorProcess::kgPort = 9327;
+const Int_t AliMonitorProcess::fgkPort = 9327;
//_____________________________________________________________________________
-AliMonitorProcess::AliMonitorProcess(const char* alienDir,
+AliMonitorProcess::AliMonitorProcess(
+#if ROOT_VERSION_CODE <= 199169 // 3.10/01
+ const char* /*alienHost*/,
+#else
+ const char* alienHost,
+#endif
+ const char* alienDir,
const char* fileNameGalice)
{
// initialize the monitoring process and the monitor histograms
+#if ROOT_VERSION_CODE <= 199169 // 3.10/01
fGrid = TGrid::Connect("alien", gSystem->Getenv("USER"));
+#else
+ fGrid = TGrid::Connect(alienHost, gSystem->Getenv("USER"));
+#endif
if (!fGrid || fGrid->IsZombie() || !fGrid->IsConnected()) {
delete fGrid;
Fatal("AliMonitorProcess", "could not connect to alien");
}
+#if ROOT_VERSION_CODE <= 199169 // 3.10/01
fGrid->cd(alienDir);
+#else
+ fAlienDir = alienDir;
+#endif
fLogicalFileName = "";
fFileName = "";
fRunLoader->LoadgAlice();
gAlice = fRunLoader->GetAliRun();
if (!gAlice) Fatal("AliMonitorProcess", "no gAlice object found");
- AliITS* ITS = (AliITS*) gAlice->GetModule("ITS");
- if (!ITS) Fatal("AliMonitorProcess", "no ITS detector found");
- fITSgeom = ITS->GetITSgeom();
+ AliITS* its = (AliITS*) gAlice->GetModule("ITS");
+ if (!its) Fatal("AliMonitorProcess", "no ITS detector found");
+ fITSgeom = its->GetITSgeom();
if (!fITSgeom) Fatal("AliMonitorProcess", "could not load ITS geometry");
#ifdef ALI_HLT
fRunNumber = 0;
fSubRunNumber = 0;
fNEvents = 0;
- fNEventsMin = 2;
+ fNEventsMin = 1;
fWriteHistoList = kFALSE;
fTopFolder = new TFolder("Monitor", "monitor histograms");
fMonitors.Add(new AliMonitorV0s);
#ifdef ALI_HLT
fMonitors.Add(new AliMonitorHLT(fTPCParam));
+ fMonitors.Add(new AliMonitorHLTHough(fTPCParam));
#endif
for (Int_t iMonitor = 0; iMonitor < fMonitors.GetEntriesFast(); iMonitor++) {
}
gROOT->cd();
- fServerSocket = new TServerSocket(kgPort, kTRUE);
+ fServerSocket = new TServerSocket(fgkPort, kTRUE);
fServerSocket->SetOption(kNoBlock, 1);
fDisplaySocket = NULL;
CheckForConnections();
fStopping = kFALSE;
}
+//_____________________________________________________________________________
+AliMonitorProcess::AliMonitorProcess(const AliMonitorProcess& process) :
+ TObject(process)
+{
+ Fatal("AliMonitorProcess", "copy constructor not implemented");
+}
+
+//_____________________________________________________________________________
+AliMonitorProcess& AliMonitorProcess::operator = (const AliMonitorProcess&
+ /*process*/)
+{
+ Fatal("operator =", "assignment operator not implemented");
+ return *this;
+}
+
//_____________________________________________________________________________
AliMonitorProcess::~AliMonitorProcess()
{
fSockets.Delete();
delete fDisplaySocket;
+#if ROOT_VERSION_CODE <= 199169 // 3.10/01
fGrid->Close();
+#endif
delete fGrid;
fFile->Close();
#ifdef ALI_HLT
delete fHLT;
+ delete fHLTHough;
#endif
}
//_____________________________________________________________________________
void AliMonitorProcess::SetStatus(EStatus status)
{
+// set the current status and process system events
+
fStatus = status;
gSystem->ProcessEvents();
}
{
// check whether a new file was registered in alien
+#if ROOT_VERSION_CODE <= 199169 // 3.10/01
TGridResult* result = fGrid->Ls();
+#else
+ Grid_ResultHandle_t handle = fGrid->OpenDir(fAlienDir);
+ if (!handle) {
+ Error("CheckForNewFile", "could not open alien directory %s",
+ fAlienDir.Data());
+ return kFALSE;
+ }
+ TGridResult* result = fGrid->CreateGridResult(handle);
+#endif
Long_t maxDate = -1;
Long_t maxTime = -1;
TString fileName;
+#if ROOT_VERSION_CODE <= 199169 // 3.10/01
while (const char* entry = result->Next()) {
+#else
+ while (Grid_Result_t* resultEntry = result->Next()) {
+ const char* entry = resultEntry->name.c_str();
+#endif
// entry = host_date_time.root
TString entryCopy(entry);
char* p = const_cast<char*>(entryCopy.Data());
}
}
+ delete result;
if (maxDate < 0) return kFALSE; // no files found
if (fLogicalFileName.CompareTo(fileName) == 0) return kFALSE; // no new file
fLogicalFileName = fileName;
- TGridResult* result2 = fGrid->GetPhysicalFileNames(fLogicalFileName.Data());
- fFileName = result2->Next();
+#if ROOT_VERSION_CODE <= 199169 // 3.10/01
+ result = fGrid->GetPhysicalFileNames(fLogicalFileName.Data());
+ fFileName = result->Next();
+#else
+ fileName = fAlienDir + "/" + fLogicalFileName;
+ handle = fGrid->GetPhysicalFileNames(fileName.Data());
+ if (!handle) {
+ Error("CheckForNewFile", "could not get physical file names for %s",
+ fileName.Data());
+ return kFALSE;
+ }
+ result = fGrid->CreateGridResult(handle);
+ result->Reset();
+ Grid_Result_t* resultEntry = result->Next();
+ if (!resultEntry) {
+ Error("CheckForNewFile", "could not get physical file names for %s",
+ fileName.Data());
+ return kFALSE;
+ }
+ fFileName = resultEntry->name2.c_str();
+ fFileName.ReplaceAll("castor:/", "rfio:/");
+#endif
+ delete result;
return kTRUE;
}
nEvents, fFileName.Data());
#ifdef ALI_HLT
CreateHLT(fFileName);
+ CreateHLTHough(fFileName);
#endif
// loop over the events
if (fStopping) break;
}
+ // monitor only central physics events
+ if (rawReader.GetType() != 7) continue;
+ if ((rawReader.GetAttributes()[0] & 0x02) == 0) continue;
+
if (!ReconstructTPC(&rawReader)) return kFALSE;
if (fStopping) break;
if (!ReconstructITS(&rawReader)) return kFALSE;
if (fStopping) break;
if (!ReconstructHLT(iEvent)) return kFALSE;
if (fStopping) break;
+ if (!ReconstructHLTHough(iEvent)) return kFALSE;
+ if (fStopping) break;
if (fDisplaySocket) fDisplaySocket->Send("new event");
//_____________________________________________________________________________
-UInt_t AliMonitorProcess::GetEventPeriodNumber()
+UInt_t AliMonitorProcess::GetEventPeriodNumber() const
{
// get the period number from the event id
}
//_____________________________________________________________________________
-UInt_t AliMonitorProcess::GetEventOrbitNumber()
+UInt_t AliMonitorProcess::GetEventOrbitNumber() const
{
// get the orbit number from the event id
}
//_____________________________________________________________________________
-UInt_t AliMonitorProcess::GetEventBunchNumber()
+UInt_t AliMonitorProcess::GetEventBunchNumber() const
{
// get the bunch number from the event id
}
//_____________________________________________________________________________
-Int_t AliMonitorProcess::GetNumberOfEvents(const char* fileName)
+Int_t AliMonitorProcess::GetNumberOfEvents(const char* fileName) const
{
// determine the number of events in the given raw data file
fHLT = new AliLevel3(name);
fHLT->Init("./", AliLevel3::kRaw, 1);
- fHLT->SetClusterFinderParam(0, 0, kTRUE);
+ fHLT->SetClusterFinderParam(-1, -1, kTRUE);
- Int_t phi_segments = 50;
- Int_t eta_segments = 100;
+ Int_t phiSegments = 50;
+ Int_t etaSegments = 100;
Int_t trackletlength = 3;
- Int_t tracklength = 40;//40 or 5
+ Int_t tracklength = 20;//40 or 5
Int_t rowscopetracklet = 2;
- Int_t rowscopetrack = 2;
- Double_t min_pt_fit = 0;
- Double_t maxangle = 1.31;
+ Int_t rowscopetrack = 10;
+ Double_t minPtFit = 0;
+ Double_t maxangle = 0.1745;
Double_t goodDist = 5;
- Double_t maxphi = 100;
- Double_t maxeta = 100;
+ Double_t maxphi = 0.1;
+ Double_t maxeta = 0.1;
Double_t hitChi2Cut = 15;//100 or 15
Double_t goodHitChi2 = 5;//20 or 5
Double_t trackChi2Cut = 10;//50 or 10
- fHLT->SetTrackerParam(phi_segments, eta_segments,
+ fHLT->SetTrackerParam(phiSegments, etaSegments,
trackletlength, tracklength,
rowscopetracklet, rowscopetrack,
- min_pt_fit, maxangle, goodDist, hitChi2Cut,
+ minPtFit, maxangle, goodDist, hitChi2Cut,
goodHitChi2, trackChi2Cut, 50, maxphi, maxeta, kTRUE);
fHLT->WriteFiles("./hlt/");
}
+
+//_____________________________________________________________________________
+void AliMonitorProcess::CreateHLTHough(const char* fileName)
+{
+
+// create the HLT Hough transform (L3Hough) object
+
+ if (fHLTHough) delete fHLTHough;
+
+ char name[256];
+ strcpy(name, fileName);
+
+ fHLTHough = new AliL3Hough();
+ fHLTHough->SetThreshold(4);
+ fHLTHough->SetTransformerParams(140,150,0.5,-1);
+ fHLTHough->SetPeakThreshold(9000,-1);// or 6000
+ fHLTHough->Init("./", kFALSE, 50, kFALSE,0,name);
+ fHLTHough->SetAddHistograms();
+ // fHLTHough->GetMaxFinder()->SetThreshold(14000);
+
+}
#endif
//_____________________________________________________________________________
#endif
}
+//_____________________________________________________________________________
+Bool_t AliMonitorProcess::ReconstructHLTHough(
+#ifdef ALI_HLT
+ Int_t iEvent
+#else
+ Int_t /* iEvent */
+#endif
+)
+{
+// run the HLT Hough transformer
+
+ SetStatus(kRecHLT);
+
+#ifndef ALI_HLT
+ Warning("ReconstructHLTHough", "the code was compiled without HLT support");
+ return kTRUE;
+
+#else
+ gSystem->Exec("rm -rf hlt/hough");
+ gSystem->MakeDirectory("hlt/hough");
+ gSystem->Exec("rm -rf hlt/fitter");
+ gSystem->MakeDirectory("hlt/fitter");
+ if (!fHLTHough) return kFALSE;
+
+ // fHLTHough->Process(0, 35);
+ // Loop over TPC sectors and process the data
+ for(Int_t i=0; i<=35; i++)
+ {
+ fHLTHough->ReadData(i,iEvent);
+ fHLTHough->Transform();
+ // if(fHLTHough->fAddHistograms)
+ fHLTHough->AddAllHistograms();
+ fHLTHough->FindTrackCandidates();
+ fHLTHough->AddTracks();
+ }
+ fHLTHough->WriteTracks("./hlt/hough");
+
+ // Run cluster fitter
+ AliL3ClusterFitter *fitter = new AliL3ClusterFitter("./hlt");
+
+ // Set debug flag for the cluster fitter
+ // fitter->Debug();
+
+ // Setting fitter parameters
+ fitter->SetInnerWidthFactor(1,1.5);
+ fitter->SetOuterWidthFactor(1,1.5);
+ fitter->SetNmaxOverlaps(5);
+
+ //fitter->SetChiSqMax(5,kFALSE); //isolated clusters
+ fitter->SetChiSqMax(5,kTRUE); //overlapping clusters
+
+ Int_t rowrange[2] = {0,AliL3Transform::GetNRows()-1};
+
+ // Takes input from global hough tracks produced by HT
+ fitter->LoadSeeds(rowrange,kFALSE,iEvent);
+
+ UInt_t ndigits;
+
+ for(Int_t islice = 0; islice <= 35; islice++)
+ {
+ for(Int_t ipatch = 0; ipatch < AliL3Transform::GetNPatches(); ipatch++)
+ {
+ // Read digits
+ fHLTHough->GetMemHandler(ipatch)->Free();
+ fHLTHough->GetMemHandler(ipatch)->Init(islice,ipatch);
+ AliL3DigitRowData *digits = (AliL3DigitRowData *)fHLTHough->GetMemHandler(ipatch)->AliAltroDigits2Memory(ndigits,iEvent);
+
+ fitter->Init(islice,ipatch);
+ fitter->SetInputData(digits);
+ fitter->FindClusters();
+ fitter->WriteClusters();
+ }
+ }
+
+ // Refit of the clusters
+ AliL3Vertex vertex;
+ //The seeds are the input tracks from circle HT
+ AliL3TrackArray *tracks = fitter->GetSeeds();
+ AliL3Fitter *ft = new AliL3Fitter(&vertex,1);
+
+ ft->LoadClusters("./hlt/fitter/",iEvent,kFALSE);
+ for(Int_t i=0; i<tracks->GetNTracks(); i++)
+ {
+ AliL3Track *track = tracks->GetCheckedTrack(i);
+ if(!track) continue;
+ if(track->GetNHits() < 20) continue;
+ ft->SortTrackClusters(track);
+ ft->FitHelix(track);
+ ft->UpdateTrack(track);
+ }
+ delete ft;
+
+ //Write the final tracks
+ fitter->WriteTracks(20);
+
+ delete fitter;
+
+ // remove the event number from the file names
+ char command[256];
+ sprintf(command, "rename tracks_%d tracks hlt/hough/*.raw", iEvent);
+ gSystem->Exec(command);
+ sprintf(command, "rename tracks_%d tracks hlt/fitter/*.raw", iEvent);
+ gSystem->Exec(command);
+ sprintf(command, "rename points_%d points hlt/fitter/*.raw", iEvent);
+ gSystem->Exec(command);
+ return kTRUE;
+#endif
+}
//_____________________________________________________________________________
Bool_t AliMonitorProcess::WriteHistos()