1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 ///////////////////////////////////////////////////////////////////////////////
20 // This class creates and fills the monitor histograms for the ITS //
22 ///////////////////////////////////////////////////////////////////////////////
30 #include "AliITSRawStreamSDD.h"
31 #include "AliITSRawStreamSSD.h"
32 #include "AliITSclusterV2.h"
33 #include "AliITSgeom.h"
34 #include "AliMonitorHisto.h"
35 #include "AliMonitorITS.h"
36 #include "AliRawReader.h"
37 #include "AliRunLoader.h"
39 ClassImp(AliMonitorITS)
42 //_____________________________________________________________________________
43 AliMonitorITS::AliMonitorITS(AliITSgeom* geom)
45 // create a ITS monitor object with the given geometry
51 //_____________________________________________________________________________
52 AliMonitorITS::AliMonitorITS(const AliMonitorITS& monitor) :
55 AliFatal("copy constructor not implemented");
58 //_____________________________________________________________________________
59 AliMonitorITS& AliMonitorITS::operator = (const AliMonitorITS& /*monitor*/)
61 AliFatal("assignment operator not implemented");
65 //_____________________________________________________________________________
66 void AliMonitorITS::CreateHistos(TFolder* folder)
68 // create the ITS monitor histograms
70 fFolder = folder->AddFolder("ITS", "ITS");
72 fSDDDigitsCharge = CreateHisto1("SDDDigitsCharge",
73 "charge distribution of SDD digits",
74 128, 0, 256, "charge", "#Delta N/N",
75 AliMonitorHisto::kNormEvents);
77 fSSDDigitsCharge = CreateHisto1("SSDDigitsCharge",
78 "charge distribution of SSD digits",
79 100, 0, 200, "charge", "#Delta N/N",
80 AliMonitorHisto::kNormEvents);
82 fSDDClustersCharge = CreateHisto1("SDDClustersCharge",
83 "charge distribution of SDD clusters",
84 100, 0, 200, "charge", "#Delta N/N",
85 AliMonitorHisto::kNormEvents);
87 fSSDClustersCharge = CreateHisto1("SSDClustersCharge",
88 "charge distribution of SSD clusters",
89 100, 0, 400, "charge", "#Delta N/N",
90 AliMonitorHisto::kNormEvents);
92 Int_t nModulesSPD = fGeom->GetLastSPD() - fGeom->GetStartSPD() + 1;
93 fSPDNClustersVsModule = CreateHisto1("SPDNClustersVsModule",
94 "mean number of clusters per SPD module",
95 nModulesSPD, -0.5, nModulesSPD-0.5,
96 "module", "<N_{clusters}>",
97 AliMonitorHisto::kNormEvents);
99 Int_t nModulesSDD = fGeom->GetLastSDD() - fGeom->GetStartSDD() + 1;
100 fSDDNClustersVsModule = CreateHisto1("SDDNClustersVsModule",
101 "mean number of clusters per SDD module",
102 nModulesSDD, -0.5, nModulesSDD-0.5,
103 "module", "<N_{clusters}>",
104 AliMonitorHisto::kNormEvents);
106 Int_t nModulesSSD = fGeom->GetLastSSD() - fGeom->GetStartSSD() + 1;
107 fSSDNClustersVsModule = CreateHisto1("SSDNClustersVsModule",
108 "mean number of clusters per SSD module",
109 nModulesSSD, -0.5, nModulesSSD-0.5,
110 "module", "<N_{clusters}>",
111 AliMonitorHisto::kNormEvents);
113 Int_t nLayer = fGeom->GetNlayers();
114 fNClustersVsLayer = CreateHisto1("NClustersVsLayer",
115 "mean number of clusters per layer",
116 nLayer, 0.5, nLayer+0.5,
117 "layer", "<N_{clusters}>",
118 AliMonitorHisto::kNormEvents);
120 fNTracks = CreateHisto1("NTracks", "number of tracks per event",
121 300, 0, 30000, "N_{tracks}", "#Delta N/N",
122 AliMonitorHisto::kNormEvents);
124 fNTracksITSTPC = CreateHisto2("NTracksTPCITS",
125 "correlation of number of ITS and TPC tracks per event",
126 30, 0, 30000, 30, 0, 30000,
127 "N_{tracks}(ITS)", "N_{tracks}(TPC)",
128 "#Delta N/N", AliMonitorHisto::kNormEvents);
130 fTrackPt = CreateHisto1("TrackPt", "pt distribution of tracks",
131 90, 0, 3, "p_{t} [GeV/c]", "#Delta N/N",
132 AliMonitorHisto::kNormEntries);
134 fTrackEta = CreateHisto1("TrackEta", "eta distribution of tracks",
135 100, -2, 2, "#eta", "#Delta N/N",
136 AliMonitorHisto::kNormEntries);
138 fTrackPhi = CreateHisto1("TrackPhi", "phi distribution of tracks",
139 120, -180, 180, "#phi [#circ]", "#Delta N/N",
140 AliMonitorHisto::kNormEntries);
142 fTrackDEdxVsP = CreateHisto2("TrackDEdxVsP", "dE/dx of tracks",
143 100, 0, 3, 100, 0, 150,
144 "p [GeV/c]", "dE/dx", "#Delta N/N",
145 AliMonitorHisto::kNormEntries);
149 //_____________________________________________________________________________
150 void AliMonitorITS::FillHistos(AliRunLoader* runLoader,
151 AliRawReader* rawReader, AliESD* esd)
153 // fill the ITS monitor histogrms
156 AliITSRawStreamSDD inputSDD(rawReader);
157 while (inputSDD.Next()) {
158 fSDDDigitsCharge->Fill(inputSDD.GetSignal());
162 AliITSRawStreamSSD inputSSD(rawReader);
163 while (inputSSD.Next()) {
164 fSSDDigitsCharge->Fill(inputSSD.GetSignal());
167 AliLoader* itsLoader = runLoader->GetLoader("ITSLoader");
168 if (!itsLoader) return;
171 itsLoader->LoadRecPoints();
172 TTree* clustersTree = itsLoader->TreeR();
173 if (!clustersTree) return;
174 TClonesArray* clusters = new TClonesArray("AliITSclusterV2");
175 clustersTree->SetBranchAddress("Clusters", &clusters);
177 for (Int_t iModule = 0; iModule < clustersTree->GetEntries(); iModule++) {
178 clustersTree->GetEntry(iModule);
179 Int_t iLayer, iLadder, iDetector;
180 fGeom->GetModuleId(iModule, iLayer, iLadder, iDetector);
182 for (Int_t j = 0; j < clusters->GetEntriesFast(); j++) {
183 AliITSclusterV2* cluster = (AliITSclusterV2*) clusters->At(j);
186 fSPDNClustersVsModule->Fill(iModule - fGeom->GetStartSPD());
190 fSDDClustersCharge->Fill(cluster->GetQ());
191 fSDDNClustersVsModule->Fill(iModule - fGeom->GetStartSDD());
195 fSSDClustersCharge->Fill(cluster->GetQ());
196 fSSDNClustersVsModule->Fill(iModule - fGeom->GetStartSSD());
200 fNClustersVsLayer->Fill(iLayer);
205 itsLoader->UnloadRecPoints();
209 Int_t nTPCTracks = 0;
210 for (Int_t i = 0; i < esd->GetNumberOfTracks(); i++) {
211 AliESDtrack* track = esd->GetTrack(i);
212 if (!track) continue;
213 if ((track->GetStatus() | AliESDtrack::kTPCin) != 0) nTPCTracks++;
214 if ((track->GetStatus() | AliESDtrack::kITSin) == 0) continue;
218 track->GetPxPyPz(pxyz);
219 TVector3 pTrack(pxyz);
220 fTrackPt->Fill(pTrack.Pt());
221 fTrackEta->Fill(pTrack.Eta());
222 fTrackPhi->Fill(pTrack.Phi() * TMath::RadToDeg());
223 fTrackDEdxVsP->Fill(pTrack.Mag(), track->GetITSsignal());
225 fNTracks->Fill(nTracks);
226 fNTracksITSTPC->Fill(nTracks, nTPCTracks);