From bb72d99069b8b4c4efdba019a11c03262de3d834 Mon Sep 17 00:00:00 2001 From: cvetan Date: Fri, 24 Jul 2009 12:27:05 +0000 Subject: [PATCH 1/1] Correct shell scipt and macro to be used with 2009 data. --- test/cosmic/rec.C | 127 +++++---------------------------------------- test/cosmic/rec.sh | 4 +- 2 files changed, 15 insertions(+), 116 deletions(-) diff --git a/test/cosmic/rec.C b/test/cosmic/rec.C index 8d54b56278d..b2fc5c0e76d 100644 --- a/test/cosmic/rec.C +++ b/test/cosmic/rec.C @@ -1,14 +1,11 @@ -void rec(const char *filename="raw.root", const Int_t mfield=1) +void rec(const char *filename="raw.root") { ///////////////////////////////////////////////////////////////////////////////////////// // - // Second version of the reconstruction - // script for the 2008 cosmic data (LHC08b) + // First version of the reconstruction + // script for the 2009 data (LHC09b) // ///////////////////////////////////////////////////////////////////////////////////////// - //AliLog::SetGlobalLogLevel(AliLog::kWarning); - AliLog::SetGlobalLogLevel(AliLog::kError); - gSystem->Load("libRAliEn.so"); gSystem->Load("libNet.so"); gSystem->Load("libMonaLisa.so"); @@ -17,130 +14,32 @@ void rec(const char *filename="raw.root", const Int_t mfield=1) // Set the CDB storage location AliCDBManager * man = AliCDBManager::Instance(); -// man->SetDefaultStorage("local://$ALICE_ROOT/OCDB"); - man->SetDefaultStorage("alien://folder=/alice/data/2008/LHC08d/OCDB/"); - man->SetSpecificStorage("ITS/Calib/*","local://$ALICE_ROOT/OCDB"); + man->SetDefaultStorage("raw://"); - // Example in case a specific CDB storage is needed - // man->SetSpecificStorage("ITS/Calib/MapsAnodeSDD","local://$ALICE_ROOT/OCDB"); - - AliQA::SetQARefStorage("local://$ALICE_ROOT/QAref") ; - // Reconstruction settings AliReconstruction rec; - // 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(0.1); // [cm] - itsRecoParam->SetSAUseAllClusters(); - rec.SetRecoParam("ITS",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(); - rec.SetRecoParam("TPC",tpcRecoParam); - AliTPCReconstructor::SetStreamLevel(1); - - // TRD setting - // Settings for the TRD Raw Reader - AliTRDrawStreamBase::SetRawStreamVersion("TB"); - - // TRD reconstruction params - AliTRDrecoParam *fTRDrecoParam = AliTRDrecoParam::GetCosmicTestParam(); - rec.SetRecoParam("TRD",fTRDrecoParam); - AliTRDtrackerV1::SetNTimeBins(30); - - // PHOS settings - AliPHOSRecoParam* recPHOS = new AliPHOSRecoParam(); - recPHOS->SetEMCSubtractPedestals(kTRUE); - recPHOS->SetEMCMinE(0.05); - recPHOS->SetEMCClusteringThreshold(0.10); - rec.SetRecoParam("PHOS",recPHOS); - - // T0 settings - //AliLog::SetModuleDebugLevel("T0", 10); + // QA options + rec.SetRunQA("ALL:ALL") ; + rec.SetQARefDefaultStorage("local://$ALICE_ROOT/QAref") ; - // 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"); - rec.SetRecoParam("MUON",muonRecoParam); - - // Tracking settings - AliMagF* field; - if (mfield) - field = new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG); - else - field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k2kG); - - Double_t mostProbPt=0.35; - AliExternalTrackParam::SetMostProbablePt(mostProbPt); + AliGRPRecoParam *grpRecoParam = AliGRPRecoParam::GetLowFluxParam(); + grpRecoParam->SetVertexerTracksConstraintITS(kFALSE); + grpRecoParam->SetVertexerTracksConstraintTPC(kFALSE); + grpRecoParam->SetMostProbablePt(3.0); + rec.SetRecoParam("GRP",grpRecoParam); // AliReconstruction settings rec.SetWriteESDfriend(kTRUE); rec.SetWriteAlignmentData(); rec.SetInput(filename); - // rec.SetRunReconstruction("ALL"); + rec.SetRunReconstruction("ALL"); rec.SetUseTrackingErrorsForAlignment("ITS"); - // In case some detectors have to be switched off... - rec.SetRunReconstruction("ITS TPC TRD TOF HMPID PHOS MUON FMD PMD T0 VZERO ZDC ACORDE"); - - // Enable vertex finder - it is needed for cosmic track reco - rec.SetRunVertexFinder(kTRUE); - - // To be enabled if some equipment IDs are not set correctly by DAQ - // rec.SetEquipmentIdMap("EquipmentIdMap.data"); - - // Detector options if any - rec.SetOption("ITS","cosmics"); - rec.SetOption("MUON","SAVEDIGITS"); - rec.SetOption("T0","cosmic"); - - // Enabled when CTP readout starts - rec.SetFillTriggerESD(kTRUE); - - // all events in one single file - rec.SetNumberOfEventsPerFile(-1); - // switch off cleanESD rec.SetCleanESD(kFALSE); - //rec.SetEventRange(0,15); - - rec.SetRunQA(":"); - rec.SetRunGlobalQA(kTRUE); AliLog::Flush(); rec.Run(); - //cout << "-----------------------------------------------------------------" << endl; - //cout << "-----------------------------------------------------------------" << endl; - //cout << "--------- Reconstruction Completed. Start merging QAs -----------" << endl; - //cout << "-----------------------------------------------------------------" << endl; - //cout << "-----------------------------------------------------------------" << endl; - //AliQADataMakerSteer qas; - //qas.Merge(); } diff --git a/test/cosmic/rec.sh b/test/cosmic/rec.sh index df03e180eac..287be52be68 100755 --- a/test/cosmic/rec.sh +++ b/test/cosmic/rec.sh @@ -12,7 +12,7 @@ # SET THE FOLLOWING PARAMETERS IF NEEDED: # --------------------------------------- -YEAR=08 +YEAR=09 DIALOG=/usr/bin/dialog # --------------------------------------- @@ -33,7 +33,7 @@ TITLE="Standalone reconstruction of Grid rawdata chunks. v$VERSION" # Retrieve the list of chunks from AliEn....... BASEDIR="/alice/data/20"$YEAR PATTERN="/raw/"$YEAR"0000"$RUNNUM"*0.root" -aliensh -c "gbbox find $BASEDIR $PATTERN" | head --lines=-1 > collection.tmp +gbbox find $BASEDIR $PATTERN | head --lines=-1 > collection.tmp [ $(stat -c%s collection.tmp) -eq 0 ] && { echo "No chunks found for the given run"; exit 1; } rm -f collection.tmp2 -- 2.39.3