]>
Commit | Line | Data |
---|---|---|
807998f3 | 1 | //============================================================================= |
2 | // | |
3 | // *** AddTaskTRDCalib | |
4 | // | |
5 | // This macros setup the TRD calibration task | |
6 | // | |
7 | //============================================================================= | |
8 | ||
9 | AliAnalysisTask *AddTaskTRDCalib(Int_t runNumber) | |
10 | { | |
11 | gSystem->Load("libTRDcalib"); | |
12 | // pointer to the analysis manager | |
13 | AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager(); | |
14 | if (!mgr) { | |
15 | Error("AddTaskTRDCalib", "No analysis manager to connect to."); | |
16 | return NULL; | |
17 | } | |
18 | ||
19 | // check the input handler | |
20 | if (!mgr->GetInputEventHandler()) { | |
21 | ::Error("AddTask", "This task requires an input event handler"); | |
22 | return NULL; | |
23 | } | |
24 | ||
c1ca05ad | 25 | // check the presence of the detectors |
26 | AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/GRP/Data"); | |
27 | AliGRPObject* grpData = dynamic_cast<AliGRPObject*>(entry->GetObject()); | |
28 | if (!grpData) {printf("Failed to get GRP data for run",runNumber); return;} | |
29 | Int_t activeDetectors = grpData->GetDetectorMask(); | |
30 | TString detStr = AliDAQ::ListOfTriggeredDetectors(activeDetectors); | |
31 | ||
32 | ||
807998f3 | 33 | //////////////////////////////////////////// |
34 | // Number of timebins | |
35 | /////////////////////////////////////////// | |
36 | AliTRDcalibDB *calib = AliTRDcalibDB::Instance(); | |
37 | Int_t nbOfTimeBins = calib->GetNumberOfTimeBinsDCS(); | |
38 | if(nbOfTimeBins < 0) nbOfTimeBins = 27; | |
39 | //////////////////////////////////////////// | |
40 | // | |
41 | ///////////////////////////////////////////// | |
42 | Int_t versiongain, subversiongain, versionvdrift, subversionvdrift; | |
43 | ||
c1ca05ad | 44 | |
807998f3 | 45 | ///////////////////////// |
46 | // The TRD calib Task | |
47 | ///////////////////////// | |
48 | AliTRDCalibTask *calibTask = new AliTRDCalibTask(); | |
49 | ||
50 | // Disabling TRD CPAss0 as per https://savannah.cern.ch/bugs/?88813 | |
51 | //calibTask->SetMaxEvent(-1); | |
52 | ||
a88f27f8 | 53 | calibTask->SetMaxNbTracks(999999999); |
807998f3 | 54 | calibTask->SetHisto2d(kTRUE); |
55 | calibTask->SetVector2d(kFALSE); | |
56 | calibTask->SetVdriftLinear(kTRUE); | |
57 | calibTask->SetExbAlt(kFALSE); | |
58 | calibTask->SetNz(0,0); | |
59 | calibTask->SetNrphi(0,0); | |
60 | calibTask->SetNz(0,1); | |
61 | calibTask->SetNrphi(0,1); | |
62 | calibTask->SetNz(0,2); | |
63 | calibTask->SetNrphi(0,2); | |
64 | calibTask->SetLow(0); | |
65 | calibTask->SetHigh(30); | |
66 | calibTask->SetFillZero(kFALSE); | |
67 | // now | |
68 | calibTask->AddSelectedTriggerClass("C0OB0-ABCE-NOPF-ALL"); | |
69 | calibTask->AddSelectedTriggerClass("CTRDCO2-ABCE-NOPF-CENT"); | |
70 | calibTask->AddSelectedTriggerClass("CTRDCO2-ABCE-NOPF-TRD"); | |
71 | calibTask->AddSelectedTriggerClass("CTRDCO2-ABCE-NOPF-ALL"); | |
72 | calibTask->SetReject(kTRUE); | |
73 | // before | |
74 | //calibTask->AddSelectedTriggerClass("CINT1B-ABCE-NOPF-ALL"); | |
75 | //calibTask->AddSelectedTriggerClass("CINT1WU-B-NOPF-ALL"); | |
76 | //calibTask->AddSelectedTriggerClass("CINT7WU-B-NOPF-ALL"); | |
77 | //calibTask->AddSelectedTriggerClass("CINT7WU-I-NOPF-ALL"); | |
78 | //calibTask->SetReject(kFALSE); | |
79 | //calibTask->SetDebug(2); | |
80 | calibTask->SetNbTimeBins(nbOfTimeBins); | |
5fdbc802 | 81 | calibTask->SetNumberBinCharge(100); |
807998f3 | 82 | //calibTask->SetMaxEvent(10); |
83 | //calibTask->SetThresholdP(1.0); | |
84 | calibTask->SetRequirePrimaryVertex(kTRUE); | |
85 | calibTask->SetMinNbOfContributors(1); | |
86 | calibTask->SetMaxCluster(100.0); | |
87 | calibTask->SetNbMaxCluster(2); | |
c1ca05ad | 88 | |
89 | if ( detStr.Contains("ITSSPD") && (!detStr.Contains("ITSSDD") || !detStr.Contains("ITSSSD"))) calibTask->SetUseSPDVertex(); | |
90 | ||
807998f3 | 91 | //calibTask->SetLimitChargeIntegration(kTRUE); |
92 | ||
93 | ||
94 | ///////////////////////////// | |
95 | // Track cuts | |
96 | ///////////////////////////// | |
97 | AliESDtrackCuts *trackCuts = new AliESDtrackCuts("trackcuts","trackcuts"); | |
98 | trackCuts->SetMinNClustersTPC(50); | |
99 | trackCuts->SetMaxChi2PerClusterTPC(3.5); | |
100 | //trackCuts->SetMaxCovDiagonalElements(2,2,0.5,0.5,2); | |
101 | trackCuts->SetRequireTPCRefit(kTRUE); | |
c04132ba | 102 | // trackCuts->SetRequireITSRefit(kTRUE); //removed by BD after discussion with Raphaelle, was introduced to handle pile-up |
a88f27f8 | 103 | trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); |
807998f3 | 104 | //trackCuts->SetMinNsigmaToVertex(10); |
105 | trackCuts->SetRequireSigmaToVertex(kFALSE); | |
106 | trackCuts->SetAcceptKinkDaughters(kFALSE); | |
107 | trackCuts->SetMaxDCAToVertexZ(30.0); | |
108 | trackCuts->SetMaxDCAToVertexXY(3.0); | |
109 | trackCuts->SetDCAToVertex2D(kFALSE); | |
110 | ||
111 | calibTask->SetESDtrackCuts(trackCuts); | |
112 | ||
113 | mgr->AddTask(calibTask); | |
114 | ||
115 | AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer(); | |
116 | //AliAnalysisDataContainer *cinput = mgr->GetCommonOutputContainer(); | |
117 | ||
118 | if (!cinput) cinput = mgr->CreateContainer("cchain",TChain::Class(), | |
119 | AliAnalysisManager::kInputContainer); | |
120 | ||
121 | AliAnalysisDataContainer *coutput =mgr->CreateContainer("TRDCalib",TList::Class(), AliAnalysisManager::kOutputContainer, "AliESDfriends_v1.root"); | |
122 | ||
123 | ||
124 | mgr->ConnectInput(calibTask,0,cinput); | |
125 | mgr->ConnectOutput(calibTask,1,coutput); | |
126 | return calibTask; | |
127 | ||
128 | } |