void rec() {
AliReconstruction reco;
reco.SetWriteESDfriend();
reco.SetWriteAlignmentData();
reco.SetRecoParam("ITS",AliITSRecoParam::GetLowFluxParam());
reco.SetRecoParam("TPC",AliTPCRecoParam::GetLowFluxParam());
reco.SetRecoParam("TRD",AliTRDrecoParam::GetLowFluxParam());
reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
reco.SetRecoParam("MUON",AliMUONRecoParam::GetLowFluxParam());
reco.SetRecoParam("EMCAL",AliEMCALRecParam::GetLowFluxParam());
reco.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/?cacheFold=/tmp/CDBCache?operateDisconnected=kFALSE");
reco.SetSpecificStorage("GRP/GRP/Data",
Form("local://%s",gSystem->pwd()));
// add TRD standalone tracks
reco.SetOption("TRD", "sl_tr_1");
reco.SetRunGlobalQA(kTRUE);
TStopwatch timer;
timer.Start();
reco.Run();
timer.Stop();
timer.Print();
}
void rec() {
/* AliLog::SetGlobalLogLevel(AliLog::kInfo);
AliLog::SetClassDebugLevel("AliTOFtracker",2);
AliLog::SetClassDebugLevel("AliTOFtrackerV1",2);
AliLog::SetClassDebugLevel("AliTOFGeometry",2);*/
AliReconstruction reco;
reco.SetWriteESDfriend();
reco.SetWriteAlignmentData();
// reco.SetRecoParam("ITS",AliITSRecoParam::GetHighFluxParam());
AliITSRecoParam * itsRecoParam = AliITSRecoParam::GetHighFluxParam(); //Pb-Pb
// trackleter options ...
itsRecoParam->SetTrackleterOnlyOneTrackletPerC2(kFALSE);
itsRecoParam->SetTrackleterPhiWindow(0.015); //rad
itsRecoParam->SetTrackleterZetaWindow(0.03); //cm
reco.SetRecoParam("ITS",itsRecoParam);
reco.SetRecoParam("TPC",AliTPCRecoParam::GetHighFluxParam());
reco.SetRecoParam("TRD",AliTRDrecoParam::GetHighFluxParam());
reco.SetRecoParam("TOF",AliTOFRecoParam::GetPbPbparam());
reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
reco.SetRecoParam("MUON",AliMUONRecoParam::GetHighFluxParam());
reco.SetRunVertexFinderTracks(kFALSE);
reco.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal");
reco.SetSpecificStorage("GRP/GRP/Data",
Form("local://%s",gSystem->pwd()));
TStopwatch timer;
timer.Start();
reco.Run();
timer.Stop();
timer.Print();
}
void rec(const char *filename="raw.root")
{
/////////////////////////////////////////////////////////////////////////////////////////
//
// Reconstruction script for 2011 RAW data - muon fast reco
//
/////////////////////////////////////////////////////////////////////////////////////////
// Set the CDB storage location
AliCDBManager * man = AliCDBManager::Instance();
man->SetDefaultStorage("alien://folder=/alice/data/2011/OCDB");
AliReconstruction rec;
// Generate or use the local OCDB.root file
// rec.SetFromCDBSnapshot("OCDB.root");
// Set reconstruction flags (skip detectors here if neded with -<detector name>
// rec.SetRunReconstruction("MUON ITS VZERO ZDC T0");
rec.SetRunReconstruction("MUON ITS");
// QA options
rec.SetRunQA("Global MUON:ALL") ;
rec.SetQARefDefaultStorage("local://$ALICE_ROOT/QAref") ;
// MUON only reco - recoparameters
rec.SetRecoParam("ITS",GetSpecialITSRecoParam());
// AliReconstruction settings
rec.SetWriteESDfriend(kTRUE);
rec.SetWriteAlignmentData();
// rec.SetInput(filename);
rec.SetUseTrackingErrorsForAlignment("ITS");
// switch off cleanESD
rec.SetCleanESD(kFALSE);
// Specific reco params for ZDC (why isn't this automatic?)
// rec.SetRecoParam("ZDC",AliZDCRecoParamPbPb::GetHighFluxParam(2760));
//Ignore SetStopOnError
rec.SetStopOnError(kFALSE);
AliLog::Flush();
rec.Run();
}
void rec() {
AliReconstruction MuonRec;
MuonRec.SetInput("raw.root");
MuonRec.SetRunLocalReconstruction("MUON ITS VZERO");
MuonRec.SetRunTracking("MUON ITS VZERO");
MuonRec.SetRunVertexFinder(kTRUE);
MuonRec.SetFillESD("MUON VZERO HLT");
MuonRec.SetRunQA("MUON:ALL");
MuonRec.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
// QA reference
MuonRec.SetQARefDefaultStorage("local://$ALICE_ROOT/QAref") ;
// HLT
MuonRec.SetSpecificStorage("HLT/ConfigMUON/DecisionComponent","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");
MuonRec.SetSpecificStorage("HLT/ConfigMUON/HitReconstructor","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");
MuonRec.SetSpecificStorage("HLT/ConfigMUON/MansoTrackerFSM","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");
MuonRec.SetSpecificStorage("HLT/ConfigMUON/TriggerReconstructor","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");
// CTP
MuonRec.SetSpecificStorage("GRP/CTP/Config","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");
// tracker masks
MuonRec.SetSpecificStorage("MUON/Calib/Gains","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb/rec");
MuonRec.SetSpecificStorage("MUON/Calib/Pedestals","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb/rec");
// reconstruction parameters
AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();
muonRecoParam->SaveFullClusterInESD(kTRUE,100.);
for (Int_t i = 0; i < 10; i++) {
muonRecoParam->SetDefaultNonBendingReso(i, 0.1524);
muonRecoParam->SetDefaultBendingReso(i, 0.05099);
}
//muonRecoParam->Print("FULL");
MuonRec.SetRecoParam("MUON",muonRecoParam);
MuonRec.SetNumberOfEventsPerFile(500);
TStopwatch timer;
timer.Start();
MuonRec.Run();
timer.Stop();
timer.Print();
}
void rec() {
AliReconstruction MuonRec;
MuonRec.SetInput("raw.root");
MuonRec.SetRunLocalReconstruction("MUON ITS VZERO");
MuonRec.SetRunTracking("MUON ITS VZERO");
MuonRec.SetRunVertexFinder(kTRUE);
MuonRec.SetFillESD("MUON VZERO HLT");
MuonRec.SetRunQA("MUON:ALL");
MuonRec.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
// QA reference
MuonRec.SetQARefDefaultStorage("local://$ALICE_ROOT/QAref") ;
// GRP
MuonRec.SetSpecificStorage("GRP/GRP/Data",Form("local://%s",gSystem->pwd()));
// HLT
MuonRec.SetSpecificStorage("HLT/ConfigMUON/DecisionComponent","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");
MuonRec.SetSpecificStorage("HLT/ConfigMUON/HitReconstructor","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");
MuonRec.SetSpecificStorage("HLT/ConfigMUON/MansoTrackerFSM","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");
MuonRec.SetSpecificStorage("HLT/ConfigMUON/TriggerReconstructor","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");
// CTP
MuonRec.SetSpecificStorage("GRP/CTP/Config","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");
// alignment
MuonRec.SetSpecificStorage("MUON/Align/Data","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb/ideal");
// reconstruction parameters
AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();
muonRecoParam->SaveFullClusterInESD(kTRUE,100.);
//muonRecoParam->Print("FULL");
MuonRec.SetRecoParam("MUON",muonRecoParam);
MuonRec.SetNumberOfEventsPerFile(500);
TStopwatch timer;
timer.Start();
MuonRec.Run();
timer.Stop();
timer.Print();
}
void main_recCPass1_OuterDet(const char *filename="raw.root",Int_t nevents=-1, const char *ocdb="raw://", Bool_t useFullITS=kFALSE)
{
// Load some system libs for Grid and monitoring
// Set the CDB storage location
AliCDBManager * man = AliCDBManager::Instance();
man->SetDefaultStorage(ocdb);
// Reconstruction settings
AliReconstruction rec;
// Upload CDB entries from the snapshot (local root file) if snapshot exist
if (gSystem->AccessPathName("OCDB.root", kFileExists)==0) {
man->SetDefaultStorage("local://");
man->SetRaw(kFALSE);
man->SetSnapshotMode("OCDB.root");
}
if (gSystem->AccessPathName("localOCDBaccessConfig.C", kFileExists)==0) {
gInterpreter->ProcessLine("localOCDBaccessConfig();");
}
if(!useFullITS){
// only SPD-trackletting will be done
printf("Special ITS configuration: only SPD-trackletting will be done\n");
rec.SetRecoParam("ITS",GetSpecialITSRecoParam());
}
// All friends
rec.SetFractionFriends(1.0);
// AliReconstruction settings - hardwired MB trigger for calibration
TString newfilename = filename;
newfilename += "?Trigger=kCalibOuter";
rec.SetInput(newfilename.Data());
// Set protection against too many events in a chunk (should not happen)
if (nevents>0) rec.SetEventRange(0,nevents);
// Remove recpoints after each event
rec.SetDeleteRecPoints("ITS MUON EMCAL PHOS VZERO T0");
// Switch off the V0 finder - saves time!
rec.SetRunV0Finder(kFALSE);
//
// QA options - all QA is off
//
rec.SetRunQA(":");
rec.SetRunGlobalQA(kFALSE);
// AliReconstruction settings
rec.SetWriteESDfriend(kFALSE);
if(!useFullITS) rec.SetWriteAlignmentData(kFALSE);
else rec.SetWriteAlignmentData();
rec.SetUseTrackingErrorsForAlignment("ITS");
rec.SetRunReconstruction("ITS MUON EMCAL PHOS VZERO T0");
rec.SetFillESD("ITS MUON EMCAL PHOS VZERO T0");
rec.SetCleanESD(kFALSE);
// Specific reco params for ZDC (why isn't this automatic?)
// rec.SetRecoParam("ZDC",AliZDCRecoParamPbPb::GetHighFluxParam(2760));
//Ignore SetStopOnError
rec.SetStopOnError(kFALSE);
AliLog::Flush();
rec.Run();
}
void rec() {
AliLog::SetGlobalLogLevel(AliLog::kError);
AliReconstruction reco;
reco.SetWriteESDfriend();
reco.SetWriteAlignmentData();
AliITSRecoParam * itspar = AliITSRecoParam::GetLowFluxParam();
itspar->SetStoreLikeSignV0s(kTRUE);
reco.SetRecoParam("ITS",itspar);
reco.SetRecoParam("TPC",AliTPCRecoParam::GetLowFluxParam());
reco.SetRecoParam("TRD",AliTRDrecoParam::GetLowFluxParam());
reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());
reco.SetRecoParam("MUON",AliMUONRecoParam::GetLowFluxParam());
reco.SetRecoParam("EMCAL",AliEMCALRecParam::GetLowFluxParam());
reco.SetRecoParam("GRP",AliGRPRecoParam::GetLowFluxParam());
reco.SetOption("TRD","sl_tr_1"); // Stream Level for the tracker equal to 1
// Only in case of Full misalignment
// AliGRPRecoParam *grpRecoParam = AliGRPRecoParam::GetLowFluxParam();
// grpRecoParam->SetVertexerTracksConstraintITS(kFALSE);
// grpRecoParam->SetVertexerTracksConstraintTPC(kFALSE);
// reco.SetRecoParam("GRP",grpRecoParam);
reco.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
reco.SetSpecificStorage("GRP/GRP/Data/","alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
// No write access to the OCDB => local specific storage
// reco.SetSpecificStorage("GRP/GRP/Data",
// Form("local://%s",gSystem->pwd()));
//-------------------------------------------------------------------------
// Setting the cuts for the V0 and cascade finding
// The values of the cuts below are "reasonable" for pp events
//-------------------------------------------------------------------------
Double_t v0sels[]={33, // max allowed chi2
0.05, // min allowed impact parameter for the 1st daughter
0.05, // min allowed impact parameter for the 2nd daughter
0.5, // max allowed DCA between the daughter tracks
0.99, // max allowed cosine of V0's pointing angle
0.2, // min radius of the fiducial volume
100 // max radius of the fiducial volume
};
AliV0vertexer::SetDefaultCuts(v0sels);
Double_t xisels[]={33., // max allowed chi2 (same as PDC07)
0.025, // min allowed V0 impact parameter (PDC07 was 0.05)
0.010, // "window" around the Lambda mass (PDC07 was 0.008)
0.025, // min allowed bachelor's impact parameter (PDC07 was 0.035)
0.2, // max allowed DCA between the V0 and the bachelor (PDC07 was 0.1)
0.998, // max allowed cosine of the cascade pointing angle (PDC07 was 0.9985)
0.2, // min radius of the fiducial volume (PDC07 was 0.9)
100 // max radius of the fiducial volume (same as PDC07)
};
AliCascadeVertexer::SetDefaultCuts(xisels);
reco.SetRunQA("ALL:ALL");
TStopwatch timer;
timer.Start();
reco.Run();
timer.Stop();
timer.Print();
}