-void rec(const char *filename="raw.root", const Int_t mfield=1)\r
+void rec(const char *filename="raw.root")\r
{\r
/////////////////////////////////////////////////////////////////////////////////////////\r
//\r
- // Second version of the reconstruction\r
- // script for the 2008 cosmic data (LHC08b) \r
+ // First version of the reconstruction\r
+ // script for the 2009 data (LHC09b) \r
//\r
/////////////////////////////////////////////////////////////////////////////////////////\r
- //AliLog::SetGlobalLogLevel(AliLog::kWarning);\r
- AliLog::SetGlobalLogLevel(AliLog::kError);\r
-\r
gSystem->Load("libRAliEn.so");\r
gSystem->Load("libNet.so");\r
gSystem->Load("libMonaLisa.so");\r
\r
// Set the CDB storage location\r
AliCDBManager * man = AliCDBManager::Instance();\r
- // man->SetDefaultStorage("local://$ALICE_ROOT");\r
- man->SetDefaultStorage("alien://folder=/alice/data/2008/LHC08b/OCDB/");\r
- man->SetSpecificStorage("ITS/Calib/*","local://$ALICE_ROOT");\r
+ man->SetDefaultStorage("raw://");\r
\r
- // Example in case a specific CDB storage is needed\r
- // man->SetSpecificStorage("ITS/Calib/MapsAnodeSDD","local://$ALICE_ROOT");\r
-\r
// Reconstruction settings\r
AliReconstruction rec;\r
\r
- // ITS settings\r
- AliITSRecoParam * itsRecoParam = AliITSRecoParam::GetCosmicTestParam();\r
- itsRecoParam->SetFactorSAWindowSizes(20);\r
- itsRecoParam->SetClusterErrorsParam(2);\r
- itsRecoParam->SetFindV0s(kFALSE);\r
- itsRecoParam->SetAddVirtualClustersInDeadZone(kFALSE);\r
- itsRecoParam->SetUseAmplitudeInfo(kFALSE);\r
- // In case we want to switch off a layer\r
- // itsRecoParam->SetLayerToSkip(<N>);\r
- // itsRecoParam->SetLayerToSkip(4);\r
- // itsRecoParam->SetLayerToSkip(5);\r
- itsRecoParam->SetLayerToSkip(2);\r
- itsRecoParam->SetLayerToSkip(3);\r
- //itsRecoParam->SetSAOnePointTracks();\r
- itsRecoParam->SetClusterMisalError(0.1); // [cm]\r
- itsRecoParam->SetSAUseAllClusters();\r
- rec.SetRecoParam("ITS",itsRecoParam);\r
-\r
- // TPC settings\r
- //AliLog::SetClassDebugLevel("AliTPCclustererMI",2);\r
- AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kFALSE);\r
- tpcRecoParam->SetTimeInterval(60,940);\r
- Double_t sysError[5]={0.3,1, 0.3/150., 1./150.,0.3/(150*150.)};\r
- tpcRecoParam->SetSystematicError(sysError);\r
- tpcRecoParam->SetMinMaxCutAbs(4.);\r
- tpcRecoParam->SetMinLeftRightCutAbs(6.);\r
- tpcRecoParam->SetMinUpDownCutAbs(6.);\r
- // tpcRecoParam->Dump();\r
- rec.SetRecoParam("TPC",tpcRecoParam);\r
- AliTPCReconstructor::SetStreamLevel(1);\r
-\r
- // TRD setting\r
- // Settings for the TRD Raw Reader\r
- AliTRDrawStreamBase::SetRawStreamVersion("TB");\r
- AliTRDrawStreamTB::SetNoErrorWarning();\r
- AliTRDrawStreamTB::AllowCorruptedData();\r
- AliTRDrawStreamTB::DisableStackNumberChecker();\r
- AliTRDrawStreamTB::DisableStackLinkNumberChecker();\r
- AliTRDrawStreamTB::SetSubtractBaseline(10);\r
- \r
- // TRD reconstruction params\r
- AliTRDrecoParam *fTRDrecoParam = AliTRDrecoParam::GetCosmicTestParam();\r
- rec.SetRecoParam("TRD",fTRDrecoParam);\r
- AliTRDtrackerV1::SetNTimeBins(30);\r
-\r
- // PHOS settings\r
- AliPHOSRecoParam* recPHOS = new AliPHOSRecoParam();\r
- recPHOS->SetEMCSubtractPedestals(kTRUE);\r
- recPHOS->SetEMCMinE(0.05);\r
- recPHOS->SetEMCClusteringThreshold(0.10);\r
- rec.SetRecoParam("PHOS",recPHOS);\r
-\r
- // T0 settings\r
- //AliLog::SetModuleDebugLevel("T0", 10);\r
+ // QA options\r
+ rec.SetRunQA("ALL:ALL") ;\r
+ rec.SetQARefDefaultStorage("local://$ALICE_ROOT/QAref") ;\r
\r
- // MUON settings\r
- //AliLog::SetClassDebugLevel("AliMUONRawStreamTracker",3);\r
- AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();\r
- muonRecoParam->CombineClusterTrackReco(kTRUE);\r
- muonRecoParam->SetCalibrationMode("NOGAIN");\r
- //muonRecoParam->SetClusteringMode("PEAKFIT");\r
- //muonRecoParam->SetClusteringMode("PEAKCOG");\r
- muonRecoParam->Print("FULL");\r
- rec.SetRecoParam("MUON",muonRecoParam);\r
- \r
- // Tracking settings\r
- AliMagFMaps* field;\r
- if (mfield)\r
- field = new AliMagFMaps("Maps","Maps", 2, 1., 10., AliMagFMaps::k5kG);\r
- else\r
- field = new AliMagFMaps("Maps","Maps", 2, 0., 10., 2);\r
- AliTracker::SetFieldMap(field,1);\r
- Double_t mostProbPt=0.35;\r
- AliExternalTrackParam::SetMostProbablePt(mostProbPt);\r
+ AliGRPRecoParam *grpRecoParam = AliGRPRecoParam::GetLowFluxParam();\r
+ grpRecoParam->SetVertexerTracksConstraintITS(kFALSE);\r
+ grpRecoParam->SetVertexerTracksConstraintTPC(kFALSE);\r
+ grpRecoParam->SetMostProbablePt(3.0);\r
+ rec.SetRecoParam("GRP",grpRecoParam);\r
\r
// AliReconstruction settings\r
- rec.SetUniformFieldTracking(kFALSE);\r
rec.SetWriteESDfriend(kTRUE);\r
rec.SetWriteAlignmentData();\r
rec.SetInput(filename);\r
- // rec.SetRunReconstruction("ALL");\r
+ rec.SetRunReconstruction("ALL");\r
rec.SetUseTrackingErrorsForAlignment("ITS");\r
\r
- // In case some detectors have to be switched off...\r
- rec.SetRunReconstruction("ITS TPC TRD TOF HMPID PHOS MUON FMD PMD T0 VZERO ZDC ACORDE");\r
-\r
- // Enable vertex finder - it is needed for cosmic track reco\r
- rec.SetRunVertexFinder(kTRUE);\r
-\r
- // To be enabled if some equipment IDs are not set correctly by DAQ\r
- // rec.SetEquipmentIdMap("EquipmentIdMap.data");\r
-\r
- // Detector options if any\r
- rec.SetOption("ITS","cosmics");\r
- rec.SetOption("MUON","SAVEDIGITS");\r
- rec.SetOption("T0","cosmic");\r
-\r
- // Enabled when CTP readout starts\r
- rec.SetFillTriggerESD(kTRUE);\r
-\r
- // all events in one single file\r
- rec.SetNumberOfEventsPerFile(-1);\r
-\r
// switch off cleanESD\r
rec.SetCleanESD(kFALSE);\r
\r
- //rec.SetEventRange(0,15);\r
-\r
- rec.SetRunQA("ITS TPC:ESD RECPOINT");\r
- rec.SetRunGlobalQA(kTRUE);\r
AliLog::Flush();\r
rec.Run();\r
\r
- //cout << "-----------------------------------------------------------------" << endl;\r
- //cout << "-----------------------------------------------------------------" << endl;\r
- //cout << "--------- Reconstruction Completed. Start merging QAs -----------" << endl;\r
- //cout << "-----------------------------------------------------------------" << endl;\r
- //cout << "-----------------------------------------------------------------" << endl;\r
- //AliQADataMakerSteer qas;\r
- //qas.Merge();\r
}\r