]>
Commit | Line | Data |
---|---|---|
e3990982 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | ////////////////////////////////////////////////////////// | |
17 | // Calls derived from AliAnalysisTaskESDfilter | |
18 | // Filter the ESD Events to AODs, only those events with | |
19 | // some signal in EMCAL, righ now at least a | |
20 | // cluster of high energy | |
21 | // | |
22 | // Author: Gustavo Conesa Balbastre (LPSC - Grenoble) | |
cd231d42 | 23 | // |
24 | // $Id$ | |
25 | // | |
e3990982 | 26 | ////////////////////////////////////////////////////////// |
27 | ||
28 | #include "AliESDCaloCluster.h" | |
29 | ||
30 | #include "AliAnalysisTaskESDfilterEMCALEventSelect.h" | |
31 | ||
32 | ClassImp(AliAnalysisTaskESDfilterEMCALEventSelect) | |
33 | ||
34 | //____________________________________________________________________________________ | |
35 | AliAnalysisTaskESDfilterEMCALEventSelect::AliAnalysisTaskESDfilterEMCALEventSelect() : | |
36 | AliAnalysisTaskESDfilter("ESD Filte : EMCAL selected events"), | |
37 | fEnergyCut(10), fNcellsCut (2), | |
38 | fRecoUtils(0x0), | |
39 | fGeometry(0), fGeoName("EMCAL_COMPLETE12SMV1") | |
40 | { | |
41 | // Default constructor | |
42 | ||
43 | fRecoUtils = new AliEMCALRecoUtils; | |
44 | ||
45 | } | |
46 | ||
47 | //____________________________________________________________________________________________________ | |
48 | AliAnalysisTaskESDfilterEMCALEventSelect::AliAnalysisTaskESDfilterEMCALEventSelect(const char *name) : | |
49 | AliAnalysisTaskESDfilter(name), | |
50 | fEnergyCut(10), fNcellsCut (2), | |
51 | fRecoUtils(0x0), | |
52 | fGeometry(0), fGeoName("EMCAL_COMPLETE12SMV1") | |
53 | { | |
54 | // Constructor | |
55 | ||
56 | fRecoUtils = new AliEMCALRecoUtils; | |
57 | ||
58 | } | |
59 | ||
60 | //_________________________________________________________________ | |
61 | Bool_t AliAnalysisTaskESDfilterEMCALEventSelect::AcceptEventEMCAL() | |
62 | { | |
63 | // Accept event given there is a cluster with enough energy | |
64 | ||
65 | if(!fGeometry) fGeometry = AliEMCALGeometry::GetInstance("EMCAL_COMPLETE12SMV1"); | |
66 | ||
67 | Int_t nCluster = InputEvent() -> GetNumberOfCaloClusters(); | |
68 | AliVCaloCells * caloCell = InputEvent() -> GetEMCALCells(); | |
69 | Int_t bc = InputEvent() -> GetBunchCrossNumber(); | |
70 | ||
71 | for(Int_t icalo = 0; icalo < nCluster; icalo++) | |
72 | { | |
73 | AliESDCaloCluster *clus = (AliESDCaloCluster*) (InputEvent()->GetCaloCluster(icalo)); | |
74 | ||
75 | if( ( clus->IsEMCAL() ) && ( clus->GetNCells() > fNcellsCut ) && ( clus->E() > fEnergyCut ) && | |
76 | fRecoUtils->IsGoodCluster(clus,fGeometry,caloCell,bc)) | |
77 | { | |
78 | ||
79 | // printf("Accept event %d, E %2.2f > %2.2f, nCells %d > %d \n", | |
80 | // (Int_t) Entry(),clus->E(), fEnergyCut, clus->GetNCells(), fNcellsCut); | |
81 | ||
82 | return kTRUE; | |
83 | } | |
84 | ||
85 | }// loop | |
86 | ||
87 | return kFALSE; | |
88 | ||
89 | } | |
90 | ||
91 | //_________________________________________________________________ | |
92 | void AliAnalysisTaskESDfilterEMCALEventSelect::UserExec(Option_t *) | |
93 | { | |
94 | // Main loop | |
95 | ||
96 | // Check if the events contains what we want in EMCAL, if not, | |
97 | // do not copy the ESD into AOD | |
98 | ||
99 | if(!AcceptEventEMCAL()) return ; | |
100 | ||
101 | // Continue the processing in the same way as in the ESD filter | |
102 | ||
103 | AliAnalysisTaskESDfilter::UserExec(""); | |
104 | ||
105 | } |