From 1749398fc352949972990021ba16e7bc6f3679bf Mon Sep 17 00:00:00 2001 From: schutz Date: Mon, 21 Jul 2008 14:28:12 +0000 Subject: [PATCH] updated with changes from Peter --- test/QA/recraw.C | 138 +++++++++++++++++++++++++++++++++++++---------- 1 file changed, 110 insertions(+), 28 deletions(-) diff --git a/test/QA/recraw.C b/test/QA/recraw.C index 916e09c8507..33c1d32c3c1 100644 --- a/test/QA/recraw.C +++ b/test/QA/recraw.C @@ -1,41 +1,123 @@ void recraw() { const char * kYear = "08" ; - AliCDBManager * man = AliCDBManager::Instance(); - //man->SetDefaultStorage("alien://Folder=/alice/simulation/2007/PDC07_v4-09-Rev-00/Ideal/CDB/"); + + AliLog::SetGlobalLogLevel(AliLog::kError); + + gSystem->Load("libRAliEn.so"); + gSystem->Load("libNet.so"); + gSystem->Load("libMonaLisa.so"); + new TMonaLisaWriter(0, "GridAliRoot-rec.C", 0, 0, "global"); + gSystem->Setenv("APMON_INTERVAL", "120"); + + // Set the CDB storage location + AliCDBManager * man = AliCDBManager::Instance(); man->SetDefaultStorage("local://$ALICE_ROOT"); - man->SetSpecificStorage("EMCAL/*","local://DB"); - - AliReconstruction reco; + //man->SetDefaultStorage("alien://Folder=/alice/data/2008/LHC08b/OCDB/"); - reco.SetWriteESDfriend(); - reco.SetWriteAlignmentData(); + // Example in case a specific CDB storage is needed + man->SetSpecificStorage("EMCAL/*","local://DB"); - AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetLowFluxParam(); + // ITS settings + AliITSRecoParam * itsRecoParam = AliITSRecoParam::GetCosmicTestParam(); + itsRecoParam->SetFactorSAWindowSizes(20); + itsRecoParam->SetClusterErrorsParam(2); + itsRecoParam->SetFindV0s(kFALSE); + itsRecoParam->SetAddVirtualClustersInDeadZone(kFALSE); + itsRecoParam->SetUseAmplitudeInfo(kFALSE); + // In case we want to switch off a layer + // itsRecoParam->SetLayerToSkip(); + // itsRecoParam->SetLayerToSkip(4); + // itsRecoParam->SetLayerToSkip(5); + itsRecoParam->SetLayerToSkip(2); + itsRecoParam->SetLayerToSkip(3); + //itsRecoParam->SetSAOnePointTracks(); + itsRecoParam->SetClusterMisalError(1.0); // [cm] + itsRecoParam->SetSAUseAllClusters(); + AliITSReconstructor::SetRecoParam(itsRecoParam); + + // TPC settings + //AliLog::SetClassDebugLevel("AliTPCclustererMI",2); + AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kFALSE); + tpcRecoParam->SetTimeInterval(60,940); + Double_t sysError[5]={0.3,1, 0.3/150., 1./150.,0.3/(150*150.)}; + tpcRecoParam->SetSystematicError(sysError); + tpcRecoParam->SetMinMaxCutAbs(4.); + tpcRecoParam->SetMinLeftRightCutAbs(6.); + tpcRecoParam->SetMinUpDownCutAbs(6.); + // tpcRecoParam->Dump(); AliTPCReconstructor::SetRecoParam(tpcRecoParam); - AliTPCReconstructor::SetStreamLevel(0); - reco.SetRunReconstruction("ITS TPC TRD TOF HMPID PHOS EMCAL MUON T0 VZERO FMD PMD ZDC"); - //Reco.SetInput("raw.root") ; - //AliPHOSRecoParam* recEmc = new AliPHOSRecoParamEmc(); - // recEmc->SetSubtractPedestals(kFALSE); - // AliPHOSReconstructor::SetRecoParamEmc(recEmc); - if (! reco.SetRunQA("ITS:ALL") ) { - printf("SetRunQA ERROR\n") ; - exit ; - } - reco.SetInLoopQA() ; - - //AliQA::SetQARefStorage(Form("%s%s/", AliQA::GetQARefDefaultStorage(), kYear)) ; - AliQA::SetQARefStorage("local://$ALICE_ROOT") ; - AliQA::SetQARefDataDirName(AliQA::kMONTECARLO) ; //RUN_TYPE + AliTPCReconstructor::SetStreamLevel(1); + + // TRD setting + // Settings for the TRD Raw Reader + AliTRDrawStreamBase::SetRawStreamVersion("TB"); + AliTRDrawStreamTB::SetNoErrorWarning(); + AliTRDrawStreamTB::AllowCorruptedData(); + AliTRDrawStreamTB::DisableStackNumberChecker(); + AliTRDrawStreamTB::DisableStackLinkNumberChecker(); + AliTRDrawStreamTB::SetSubtractBaseline(10); -// **** The field map settings must be the same as in Config.C ! - AliMagFMaps *field=new AliMagFMaps("Maps","Maps",2,1.,10.,AliMagFMaps::k5kG); - Bool_t uniform=kFALSE; - AliTracker::SetFieldMap(field,uniform); + // TRD reconstruction params + AliTRDrecoParam *fTRDrecoParam = AliTRDrecoParam::GetCosmicTestParam(); + AliTRDReconstructor::SetRecoParam(fTRDrecoParam); + AliTRDtrackerV1::SetNTimeBins(30); + + // PHOS settings + AliPHOSRecoParam* recEmc = new AliPHOSRecoParamEmc(); + recEmc->SetSubtractPedestals(kTRUE); + recEmc->SetMinE(0.05); + recEmc->SetClusteringThreshold(0.10); + AliPHOSReconstructor::SetRecoParamEmc(recEmc); + + // T0 settings + //AliLog::SetModuleDebugLevel("T0", 10); + + // MUON settings + //AliLog::SetClassDebugLevel("AliMUONRawStreamTracker",3); + AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam(); + muonRecoParam->CombineClusterTrackReco(kTRUE); + muonRecoParam->SetCalibrationMode("NOGAIN"); + //muonRecoParam->SetClusteringMode("PEAKFIT"); + //muonRecoParam->SetClusteringMode("PEAKCOG"); + muonRecoParam->Print("FULL"); + // AliRecoParam::Instance()->RegisterRecoParam(muonRecoParam); + + // Tracking settings + // **** The field map settings must be the same as in Config.C ! + AliMagFMaps *field=new AliMagFMaps("Maps", "Maps", 2, 1., 10., AliMagFMaps::k5kG); + Bool_t uniform = kFALSE; + AliTracker::SetFieldMap(field, uniform); + Double_t mostProbPt=0.35; + AliExternalTrackParam::SetMostProbablePt(mostProbPt); + + // AliReconstruction settings + AliReconstruction reco; + reco.SetUniformFieldTracking(uniform); + reco.SetWriteESDfriend(kTRUE); + reco.SetWriteAlignmentData(); + reco.SetInput("raw.root"); + reco.SetUseTrackingErrorsForAlignment("ITS"); + + // In case some detectors have to be switched off... + reco.SetRunReconstruction("ITS TPC TRD TOF HMPID PHOS MUON FMD PMD T0 VZERO ZDC ACORDE"); + + reco.SetRunVertexFinder(kTRUE); + + // all events in one single file + reco.SetNumberOfEventsPerFile(-1); + // switch off cleanESD + reco.SetCleanESD(kFALSE); + + reco.SetRunQA("ALL:ALL") ; + //AliQA::SetQARefStorage(Form("%s%s/", AliQA::GetQARefDefaultStorage(), kYear)) ; + AliQA::SetQARefStorage("local://$ALICE_ROOT") ; + + AliLog::Flush(); + TStopwatch timer; timer.Start(); - reco.Run("raw.root"); + reco.Run(); timer.Stop(); timer.Print(); } -- 2.39.3