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 **************************************************************************/
16 //====================================================================================================================================================
18 // Description of an ALICE muon forward track pair, i.e. a pair of AliMuonForwardTrack objects
20 // Contact author: antonio.uras@cern.ch
22 //====================================================================================================================================================
25 #include "AliMUONTrackParam.h"
26 #include "TParticle.h"
27 #include "AliMuonForwardTrack.h"
28 #include "AliMUONTrackExtrap.h"
29 #include "TClonesArray.h"
30 #include "TDatabasePDG.h"
31 #include "TLorentzVector.h"
33 #include "AliMuonForwardTrackPair.h"
35 ClassImp(AliMuonForwardTrackPair)
37 //====================================================================================================================================================
39 AliMuonForwardTrackPair::AliMuonForwardTrackPair():
41 fMuonForwardTracks(0),
45 // default constructor
47 fMuonForwardTracks = new TClonesArray("AliMuonForwardTrack", 2);
51 //====================================================================================================================================================
53 AliMuonForwardTrackPair::AliMuonForwardTrackPair(AliMuonForwardTrack *track0, AliMuonForwardTrack *track1):
55 fMuonForwardTracks(0),
59 fMuonForwardTracks = new TClonesArray("AliMuonForwardTrack", 2);
61 new ((*fMuonForwardTracks)[0]) AliMuonForwardTrack(*track0);
62 new ((*fMuonForwardTracks)[1]) AliMuonForwardTrack(*track1);
68 //====================================================================================================================================================
70 AliMuonForwardTrackPair::AliMuonForwardTrackPair(const AliMuonForwardTrackPair& trackPair):
72 fMuonForwardTracks(trackPair.fMuonForwardTracks),
73 fKinemMC(trackPair.fKinemMC)
80 //====================================================================================================================================================
82 AliMuonForwardTrackPair& AliMuonForwardTrackPair::operator=(const AliMuonForwardTrackPair& trackPair) {
86 // check assignement to self
87 if (this == &trackPair) return *this;
89 // base class assignement
90 AliMuonForwardTrackPair::operator=(trackPair);
95 fMuonForwardTracks = trackPair.fMuonForwardTracks;
101 //====================================================================================================================================================
103 void AliMuonForwardTrackPair::SetTrack(Int_t iTrack, AliMuonForwardTrack *track) {
105 if (iTrack==0 || iTrack==1) {
106 new ((*fMuonForwardTracks)[iTrack]) AliMuonForwardTrack(*track);
111 //====================================================================================================================================================
113 Double_t AliMuonForwardTrackPair::GetWeightedOffset(Double_t x, Double_t y, Double_t z) {
115 Double_t weightedOffset[2]={0};
117 weightedOffset[0] = ((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetWeightedOffset(x, y, z);
118 weightedOffset[1] = ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetWeightedOffset(x, y, z);
120 Double_t weightedOffsetDimuon = TMath::Sqrt(0.5 * (weightedOffset[0]*weightedOffset[0] + weightedOffset[1]*weightedOffset[1]));
122 return weightedOffsetDimuon;
126 //====================================================================================================================================================
128 Double_t AliMuonForwardTrackPair::GetMass(Double_t z, Int_t nClusters) {
130 Int_t idCluster[2] = {0};
132 idCluster[0] = ((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetNMFTClusters() - nClusters;
133 idCluster[1] = ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetNMFTClusters() - nClusters;
135 if (idCluster[0]<0) idCluster[0] = 0;
136 if (idCluster[1]<0) idCluster[1] = 0;
138 Double_t momentum[2] = {0};
140 AliMUONTrackParam *param0 = ((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetTrackParamAtMFTCluster(idCluster[0]);
141 AliMUONTrackParam *param1 = ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetTrackParamAtMFTCluster(idCluster[1]);
143 AliDebug(2, Form("MFT before extrap: 1st muon = (%f, %f, %f) 2nd muon = (%f, %f, %f)",
144 param0->Px(), param0->Py(), param0->Pz(),
145 param1->Px(), param1->Py(), param1->Pz()));
147 if (TMath::Abs(z)<1e6) {
148 AliDebug(2, Form("Extrapolating 1st muon from z = %f to z = %f", param0->GetZ(), z));
149 AliMUONTrackExtrap::ExtrapToZCov(param0, z);
150 AliDebug(2, Form("Extrapolating 2nd muon from z = %f to z = %f", param1->GetZ(), z));
151 AliMUONTrackExtrap::ExtrapToZCov(param1, z);
154 AliDebug(2, Form("MFT after extrap: 1st muon = (%f, %f, %f) 2nd muon = (%f, %f, %f)",
155 param0->Px(), param0->Py(), param0->Pz(),
156 param1->Px(), param1->Py(), param1->Pz()));
158 momentum[0] = (param0->P());
159 momentum[1] = (param1->P());
161 Double_t mMu = TDatabasePDG::Instance()->GetParticle("mu-")->Mass();
165 dimu.SetE(TMath::Sqrt(mMu*mMu + momentum[0]*momentum[0]) + TMath::Sqrt(mMu*mMu + momentum[1]*momentum[1]));
167 dimu.SetPx(param0->Px() + param1->Px());
168 dimu.SetPy(param0->Py() + param1->Py());
169 dimu.SetPz(param0->Pz() + param1->Pz());
175 //====================================================================================================================================================
177 Double_t AliMuonForwardTrackPair::GetMassWithoutMFT(Double_t x, Double_t y, Double_t z, Int_t nClusters) {
179 Int_t idCluster[2] = {0};
181 idCluster[0] = ((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetNMUONClusters() - nClusters;
182 idCluster[1] = ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetNMUONClusters() - nClusters;
184 if (idCluster[0]<0) idCluster[0] = 0;
185 if (idCluster[1]<0) idCluster[1] = 0;
187 AliMUONTrackParam *param0 = ((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetTrackParamAtMUONCluster(idCluster[0]);
188 AliMUONTrackParam *param1 = ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetTrackParamAtMUONCluster(idCluster[1]);
190 AliDebug(2, Form("MUON before extrap: 1st muon = (%f, %f, %f) 2nd muon = (%f, %f, %f)",
191 param0->Px(), param0->Py(), param0->Pz(),
192 param1->Px(), param1->Py(), param1->Pz()));
194 AliDebug(2, Form("Extrapolating 1st muon from z = %f to z = %f", param0->GetZ(), z));
195 AliMUONTrackExtrap::ExtrapToVertex(param0, x, y, z, 0., 0.); // this should reproduce what is done in AliMUONESDInterface::MUONToESD(...)
196 AliDebug(2, Form("Extrapolating 2nd muon from z = %f to z = %f", param1->GetZ(), z));
197 AliMUONTrackExtrap::ExtrapToVertex(param1, x, y, z, 0., 0.); // this should reproduce what is done in AliMUONESDInterface::MUONToESD(...)
199 AliDebug(2, Form("MUON after extrap: 1st muon = (%f, %f, %f) 2nd muon = (%f, %f, %f)",
200 param0->Px(), param0->Py(), param0->Pz(),
201 param1->Px(), param1->Py(), param1->Pz()));
203 Double_t momentum[2] = {0};
205 momentum[0] = param0->P();
206 momentum[1] = param1->P();
208 Double_t mMu = TDatabasePDG::Instance()->GetParticle("mu-")->Mass();
212 dimu.SetE(TMath::Sqrt(mMu*mMu + momentum[0]*momentum[0]) + TMath::Sqrt(mMu*mMu + momentum[1]*momentum[1]));
214 dimu.SetPx(param0->Px() + param1->Px());
215 dimu.SetPy(param0->Py() + param1->Py());
216 dimu.SetPz(param0->Pz() + param1->Pz());
222 //====================================================================================================================================================
224 void AliMuonForwardTrackPair::SetKinemMC() {
226 AliDebug(2, Form("MC: 1st muon = (%f, %f, %f) 2nd muon = (%f, %f, %f)",
227 ((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetMCTrackRef()->Px(),
228 ((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetMCTrackRef()->Py(),
229 ((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetMCTrackRef()->Pz(),
230 ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetMCTrackRef()->Px(),
231 ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetMCTrackRef()->Py(),
232 ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetMCTrackRef()->Pz()));
234 fKinemMC.SetE(((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetMCTrackRef()->Energy() +
235 ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetMCTrackRef()->Energy());
237 fKinemMC.SetPx(((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetMCTrackRef()->Px() +
238 ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetMCTrackRef()->Px());
240 fKinemMC.SetPy(((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetMCTrackRef()->Py() +
241 ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetMCTrackRef()->Py());
243 fKinemMC.SetPz(((AliMuonForwardTrack*) fMuonForwardTracks->At(0))->GetMCTrackRef()->Pz() +
244 ((AliMuonForwardTrack*) fMuonForwardTracks->At(1))->GetMCTrackRef()->Pz());
248 //====================================================================================================================================================
250 Bool_t AliMuonForwardTrackPair::IsResonance() {
252 Bool_t result = kFALSE;
254 Int_t labelMC[2] = {0};
255 Int_t codePDG[2] = {0};
257 for (Int_t iTrack=0; iTrack<2; iTrack++) {
258 labelMC[iTrack] = ((AliMuonForwardTrack*) fMuonForwardTracks->At(iTrack))->GetParentMCLabel(0);
259 codePDG[iTrack] = ((AliMuonForwardTrack*) fMuonForwardTracks->At(iTrack))->GetParentPDGCode(0);
262 AliDebug(1, Form("Muons' mothers: (%d, %d)", labelMC[0], labelMC[1]));
264 if (labelMC[0]==labelMC[1] && codePDG[0]==codePDG[1] && (codePDG[0]== 113 ||
268 codePDG[0]==100443 ||
270 codePDG[0]==100553 ) ) result = kTRUE;
272 if (result) AliDebug(1, Form("Pair is a resonance %d", codePDG[0]));
278 //====================================================================================================================================================