HLT/MUON/AliHLTMUONDecision.cxx

   1 /**************************************************************************
   2  * This file is property of and copyright by the ALICE HLT Project        *
   3  * All rights reserved.                                                   *
   4  *                                                                        *
   5  * Primary Authors:                                                       *
   6  *   Artur Szostak <artursz@iafrica.com>                                  *
   7  *                                                                        *
   8  * Permission to use, copy, modify and distribute this software and its   *
   9  * documentation strictly for non-commercial purposes is hereby granted   *
  10  * without fee, provided that the above copyright notice appears in all   *
  11  * copies and that both the copyright notice and this permission notice   *
  12  * appear in the supporting documentation. The authors make no claims     *
  13  * about the suitability of this software for any purpose. It is          *
  14  * provided "as is" without express or implied warranty.                  *
  15  **************************************************************************/
  16
  17 /* $Id: $ */
  18
  19 ///
  20 /// @file   AliHLTMUONDecision.cxx
  21 /// @author Artur Szostak <artursz@iafrica.com>
  22 /// @date   12 May 2008
  23 /// @brief  Implementation of the AliHLTMUONDecision class.
  24 ///
  25 /// The class is used to store the dHLT decision in a more convenient ROOT
  26 /// object format for testing and analysis.
  27 ///
  28
  29 #include "AliHLTMUONDecision.h"
  30 #include "AliHLTMUONMansoTrack.h"
  31 #include "AliLog.h"
  32 #include <cstring>
  33 #include <iostream>
  34 #include <iomanip>
  35
  36 ClassImp(AliHLTMUONDecision);
  37 ClassImp(AliHLTMUONDecision::AliTrackDecision);
  38 ClassImp(AliHLTMUONDecision::AliPairDecision);
  39
  40
  41 std::ostream& operator << (
  42                 std::ostream& stream,
  43                 const AliHLTMUONDecision& decision
  44         )
  45 {
  46 /// Stream operator for std::ostream classes.
  47 /// \param stream  The output stream object being written to.
  48 /// \param track  The dHLT decision object to print to the stream.
  49 /// \returns  Returns 'stream'.
  50
  51         stream  << "No. low/high pT: [" << decision.fNlowPt
  52                 << ", " << decision.fNhighPt
  53                 << "]; No. any/low/high unlike: [" << decision.fNunlikeAnyPt
  54                 << ", " << decision.fNunlikeLowPt
  55                 << ", " << decision.fNunlikeHighPt
  56                 << "]; No. any/low/high like = " << decision.fNlikeAnyPt
  57                 << ", " << decision.fNlikeLowPt
  58                 << ", " << decision.fNlikeHighPt
  59                 << "]; No. any/low/high mass = " << decision.fNmassAny
  60                 << ", " << decision.fNmassLow
  61                 << ", " << decision.fNmassHigh
  62                 << "]";
  63         return stream;
  64 }
  65
  66
  67 AliHLTMUONDecision::AliHLTMUONDecision(
  68                 UInt_t nLowPt, UInt_t nHiPt,
  69                 UInt_t nUnlikeAnyPt, UInt_t nUnlikeLowPt, UInt_t nUnlikeHighPt,
  70                 UInt_t nLikeAnyPt, UInt_t nLikeLowPt, UInt_t nLikeHighPt,
  71                 UInt_t nMassAny, UInt_t nMassLow, UInt_t nMassHigh
  72         ) :
  73         TObject(),
  74         fNlowPt(nLowPt),
  75         fNhighPt(nHiPt),
  76         fNunlikeAnyPt(nUnlikeAnyPt),
  77         fNunlikeLowPt(nUnlikeLowPt),
  78         fNunlikeHighPt(nUnlikeHighPt),
  79         fNlikeAnyPt(nLikeAnyPt),
  80         fNlikeLowPt(nLikeLowPt),
  81         fNlikeHighPt(nLikeHighPt),
  82         fNmassAny(nMassAny),
  83         fNmassLow(nMassLow),
  84         fNmassHigh(nMassHigh),
  85         fTrackDecisions("AliHLTMUONDecision::AliTrackDecision"),
  86         fPairDecisions("AliHLTMUONDecision::AliPairDecision")
  87 {
  88 /// Constructor for creating a dHLT decision object.
  89 /// \param nLowPt  Number of tracks above low pT cut.
  90 /// \param nHiPt   Number of tracks above high pT cut.
  91 /// \param nUnlikeAnyPt  Number of track pairs with unlike sign.
  92 /// \param nUnlikeLowPt  Number of unlike sign track pairs with pT > low cut.
  93 /// \param nUnlikeHighPt Number of unlike sign track pairs with pT > high cut.
  94 /// \param nLikeAnyPt   Number of track pairs with like sign.
  95 /// \param nLikeLowPt   Number of like sign track pairs with pT > low cut.
  96 /// \param nLikeHighPt  Number of like sign track pairs with pT > high cut.
  97 /// \param nMassAny   Number of unlike sign track pairs with invariant mass > low cut.
  98 /// \param nMassLow   Number of unlike sign track pairs with invariant mass > low mass cut and pT > low pT cut.
  99 /// \param nMassHigh  Number of unlike sign track pairs with invariant mass > high mass cut and pT > high pT cut.
 100 }
 101
 102
 103 void AliHLTMUONDecision::AddDecision(
 104                 Float_t pt, Bool_t passedLowCut, Bool_t passedHighCut,
 105                 const AliHLTMUONMansoTrack* track
 106         )
 107 {
 108 /// Add a single track decision to the dHLT trigger.
 109 /// \param pt  The calculated pT value used for the trigger decision.
 110 /// \param passedLowCut  Flag indicating if the track passed the low pT cut.
 111 /// \param passedHighCut  Flag indicating if the track passed the high pT cut.
 112 /// \param track  Pointer to the associated track object.
 113
 114         new (fTrackDecisions[fTrackDecisions.GetEntriesFast()])
 115                 AliTrackDecision(
 116                         pt, passedLowCut, passedHighCut, track
 117                 );
 118 }
 119
 120
 121 void AliHLTMUONDecision::AddDecision(const AliTrackDecision* decision)
 122 {
 123 /// Add a single track decision to the dHLT trigger.
 124
 125         new (fTrackDecisions[fTrackDecisions.GetEntriesFast()])
 126                 AliTrackDecision(*decision);
 127 }
 128
 129
 130 void AliHLTMUONDecision::AddDecision(
 131                 Float_t mass, Bool_t passedLowCut,
 132                 Bool_t passedHighCut, Bool_t unlike,
 133                 UChar_t lowPtCount, UChar_t highPtCount,
 134                 const AliHLTMUONMansoTrack* trackA,
 135                 const AliHLTMUONMansoTrack* trackB
 136         )
 137 {
 138 /// Add a track pair decision to the dHLT trigger.
 139 /// \param mass   The invariant mass of the track pair.
 140 /// \param passedLowCut  Indicates if the pair passed the low mass cut.
 141 /// \param passedHighCut  Indicates if the pair passed the high mass cut.
 142 /// \param unlike   Indicates if the tracks have opposite sign.
 143 /// \param lowPtCount  The number of tracks in the pair that passed the low pT cut.
 144 ///            Should be in the range [0..2].
 145 /// \param highPtCount  The number of tracks in the pair that passed the high pT cut.
 146 ///            Should be in the range [0..2].
 147 /// \param trackA  Pointer to the first associated track object.
 148 /// \param trackB  Pointer to the second associated track object.
 149
 150         new (fPairDecisions[fPairDecisions.GetEntriesFast()])
 151                 AliPairDecision(
 152                         mass, passedLowCut, passedHighCut, unlike,
 153                         lowPtCount, highPtCount, trackA, trackB
 154                 );
 155 }
 156
 157
 158 void AliHLTMUONDecision::AddDecision(const AliPairDecision* decision)
 159 {
 160 /// Add a track pair decision to the dHLT trigger.
 161
 162         new (fPairDecisions[fPairDecisions.GetEntriesFast()])
 163                 AliPairDecision(*decision);
 164 }
 165
 166
 167 void AliHLTMUONDecision::Print(Option_t* option) const
 168 {
 169 /// Prints the trigger decision information to standard output (screen).
 170 /// \param option  Can be one of the following:
 171 ///      - "compact" - prints in a compact format.
 172 ///      - "detail" - prints trigger information in a more detailed format.
 173 ///      - "all" - prints a full dump of the trigger object.
 174
 175         using namespace std;
 176
 177         if (option == NULL or strcmp(option, "") == 0 or
 178             strcmp(option, "compact") == 0
 179            )
 180         {
 181                 cout << *this << endl;
 182         }
 183         else if (strcmp(option, "detail") == 0)
 184         {
 185                 cout << "dHLT trigger decision scalars:" << endl;
 186                 cout << "No. tracks passed pT cut," << endl;
 187                 cout << "  low = " << fNlowPt << endl;
 188                 cout << " high = " << fNhighPt << endl;
 189                 cout << "No. unlike sign pairs," << endl;
 190                 cout << "          total = " << fNunlikeAnyPt << endl;
 191                 cout << "   pT > low cut = " << fNunlikeLowPt << endl;
 192                 cout << "  pT > high cut = " << fNunlikeHighPt << endl;
 193                 cout << "No. like sign pairs," << endl;
 194                 cout << "          total = " << fNlikeAnyPt << endl;
 195                 cout << "   pT > low cut = " << fNlikeLowPt << endl;
 196                 cout << "  pT > high cut = " << fNlikeHighPt << endl;
 197                 cout << "No. pairs with," << endl;
 198                 cout << "          invariant mass > low cut = " << fNmassAny << endl;
 199                 cout << "   invariant mass and pT > low cut = " << fNmassLow << endl;
 200                 cout << "  invariant mass and pT > high cut = " << fNmassHigh << endl;
 201
 202                 streamsize w = cout.width();
 203                 ios::fmtflags f = cout.flags();
 204                 cout << "Triggers for single tracks:" << endl;
 205                 cout << setw(10) << "" << setw(12) << "pT  " << setw(15) << "Passed pT cut" << endl;
 206                 cout    << setw(10) << "Track" << setw(12) << "(GeV/c)" << setw(6) << "low "
 207                         << setw(3) << " | " << setw(6) << "high" << endl;
 208                 for (Int_t i = 0; i < NumberOfTracks(); i++)
 209                 {
 210                         const AliTrackDecision* decision = SingleTrackDecision(i);
 211                         if (decision == NULL) continue;
 212
 213                         if (decision->Track() == NULL)
 214                         {
 215                                 cout << setw(10) << "-";
 216                         }
 217                         else
 218                         {
 219                                 cout << setw(10) << decision->Track()->Id();
 220                         }
 221
 222                         cout    << setw(12) << decision->Pt()
 223                                 << setw(6) << (decision->PassedLowPtCut() ? "yes" : "no")
 224                                 << setw(3) << "   "
 225                                 << setw(6) << (decision->PassedHighPtCut() ? "yes" : "no")
 226                                 << endl;
 227                 }
 228                 cout << "Triggers for track pairs:" << endl;
 229                 cout    << setw(20) << "Track pair" << setw(6) << "Like" << setw(12) << "mass  "
 230                         << setw(17) << "Passed mass cut"
 231                         << setw(20) << "No. with pT > than" << endl;
 232                 cout    << setw(10) << "track A" << setw(10) << "track B"
 233                         << setw(6) << "sign" << setw(12) << "(GeV/c^2)"
 234                         << setw(8) << "low " << setw(3) << " | " << setw(6) << "high"
 235                         << setw(11) << "low " << setw(3) << " | " << setw(6) << "high" << endl;
 236                 for (Int_t j = 0; j < NumberOfPairs(); j++)
 237                 {
 238                         const AliPairDecision* decision = TrackPairDecision(j);
 239                         if (decision == NULL) continue;
 240
 241                         if (decision->TrackA() == NULL)
 242                         {
 243                                 cout << setw(10) << "-";
 244                         }
 245                         else
 246                         {
 247                                 cout << setw(10) << decision->TrackA()->Id();
 248                         }
 249
 250                         if (decision->TrackB() == NULL)
 251                         {
 252                                 cout << setw(10) << "-";
 253                         }
 254                         else
 255                         {
 256                                 cout << setw(10) << decision->TrackB()->Id();
 257                         }
 258
 259                         cout    << setw(6) << (decision->LikeSign() ? "yes" : "no")
 260                                 << setw(12) << decision->Mass()
 261                                 << setw(8) << (decision->PassedLowMassCut() ? "yes" : "no")
 262                                 << setw(3) << "   "
 263                                 << setw(6) << (decision->PassedHighMassCut() ? "yes" : "no")
 264                                 << setw(11) << Int_t(decision->NumberPassedLowPtCut())
 265                                 << setw(3) << "   "
 266                                 << setw(6) << Int_t(decision->NumberPassedHighPtCut())
 267                                 << endl;
 268                 }
 269                 cout.width(w); // reset the field width to previous value.
 270                 cout.flags(f); // reset the flags to previous values.
 271         }
 272         else if (strcmp(option, "all") == 0)
 273         {
 274                 cout << "dHLT trigger decision scalars:" << endl;
 275                 cout << "No. tracks passed pT cut," << endl;
 276                 cout << "  low = " << fNlowPt << endl;
 277                 cout << " high = " << fNhighPt << endl;
 278                 cout << "No. unlike sign pairs," << endl;
 279                 cout << "          total = " << fNunlikeAnyPt << endl;
 280                 cout << "   pT > low cut = " << fNunlikeLowPt << endl;
 281                 cout << "  pT > high cut = " << fNunlikeHighPt << endl;
 282                 cout << "No. like sign pairs," << endl;
 283                 cout << "          total = " << fNlikeAnyPt << endl;
 284                 cout << "   pT > low cut = " << fNlikeLowPt << endl;
 285                 cout << "  pT > high cut = " << fNlikeHighPt << endl;
 286                 cout << "No. pairs with," << endl;
 287                 cout << "          invariant mass > low cut = " << fNmassAny << endl;
 288                 cout << "   invariant mass and pT > low cut = " << fNmassLow << endl;
 289                 cout << "  invariant mass and pT > high cut = " << fNmassHigh << endl;
 290
 291                 cout << "===============================================" << endl;
 292                 cout << "========== Triggers for single tracks: ========" << endl;
 293                 for (Int_t i = 0; i < NumberOfTracks(); i++)
 294                 {
 295                         const AliTrackDecision* decision = SingleTrackDecision(i);
 296                         if (decision == NULL) continue;
 297                         decision->Print("all");
 298                 }
 299                 cout << "===============================================" << endl;
 300                 cout << "========== Triggers for track pairs: ==========" << endl;
 301                 for (Int_t j = 0; j < NumberOfPairs(); j++)
 302                 {
 303                         const AliPairDecision* decision = TrackPairDecision(j);
 304                         if (decision == NULL) continue;
 305                         decision->Print("all");
 306                 }
 307         }
 308         else
 309         {
 310                 AliError("Unknown option specified. Can only be one of 'compact',"
 311                         " 'detail' or 'all'."
 312                 );
 313         }
 314 }
 315
 316
 317 Int_t AliHLTMUONDecision::Compare(const TObject* obj) const
 318 {
 319 /// We compare this object with 'obj' first by the trigger scalars and then
 320 /// by the signle track and track pair decision lists.
 321 /// \param obj  This is the object to compare to. It must be of type AliHLTMUONDecision.
 322 /// \returns  -1 if 'this' is smaller than 'obj', 1 if greater and zero if both
 323 ///      objects are the same.
 324
 325         if (obj->IsA() == AliHLTMUONDecision::Class())
 326         {
 327                 const AliHLTMUONDecision* d =
 328                         static_cast<const AliHLTMUONDecision*>(obj);
 329                 if (fNlowPt < d->fNlowPt) return -1;
 330                 if (fNlowPt > d->fNlowPt) return 1;
 331                 if (fNhighPt < d->fNhighPt) return -1;
 332                 if (fNhighPt > d->fNhighPt) return 1;
 333                 if (fNunlikeAnyPt < d->fNunlikeAnyPt) return -1;
 334                 if (fNunlikeAnyPt > d->fNunlikeAnyPt) return 1;
 335                 if (fNunlikeLowPt < d->fNunlikeLowPt) return -1;
 336                 if (fNunlikeLowPt > d->fNunlikeLowPt) return 1;
 337                 if (fNunlikeHighPt < d->fNunlikeHighPt) return -1;
 338                 if (fNunlikeHighPt > d->fNunlikeHighPt) return 1;
 339                 if (fNlikeAnyPt < d->fNlikeAnyPt) return -1;
 340                 if (fNlikeAnyPt > d->fNlikeAnyPt) return 1;
 341                 if (fNlikeLowPt < d->fNlikeLowPt) return -1;
 342                 if (fNlikeLowPt > d->fNlikeLowPt) return 1;
 343                 if (fNlikeHighPt < d->fNlikeHighPt) return -1;
 344                 if (fNlikeHighPt > d->fNlikeHighPt) return 1;
 345                 if (fNmassAny < d->fNmassAny) return -1;
 346                 if (fNmassAny > d->fNmassAny) return 1;
 347                 if (fNmassLow < d->fNmassLow) return -1;
 348                 if (fNmassLow > d->fNmassLow) return 1;
 349                 if (fNmassHigh < d->fNmassHigh) return -1;
 350                 if (fNmassHigh > d->fNmassHigh) return 1;
 351
 352                 // Now check the track decision arrays.
 353                 if (NumberOfTracks() < d->NumberOfTracks()) return -1;
 354                 if (NumberOfTracks() > d->NumberOfTracks()) return 1;
 355                 for (Int_t i = 0; i < NumberOfTracks(); i++)
 356                 {
 357                         Int_t result = SingleTrackDecision(i)->Compare( d->SingleTrackDecision(i) );
 358                         if (result != 0) return result;
 359                 }
 360                 if (NumberOfPairs() < d->NumberOfPairs()) return -1;
 361                 if (NumberOfPairs() > d->NumberOfPairs()) return 1;
 362                 for (Int_t j = 0; j < NumberOfPairs(); j++)
 363                 {
 364                         Int_t result = TrackPairDecision(j)->Compare( d->TrackPairDecision(j) );
 365                         if (result != 0) return result;
 366                 }
 367
 368                 // At this point everything was equal so return 0 to indicate this fact.
 369                 return 0;
 370         }
 371         else
 372         {
 373                 AliError(Form("Do not know how to compare %s to %s.",
 374                         this->ClassName(),
 375                         obj->ClassName()
 376                 ));
 377                 return -999;
 378         }
 379 }
 380
 381
 382 bool AliHLTMUONDecision::operator == (const AliHLTMUONDecision& decision) const
 383 {
 384 /// Comparison operator just compares if the scalars are the same.
 385 /// \param decision  The trigger decision object to compare to.
 386 /// \returns  true if 'this' object has the same scalars as 'decision', else false.
 387
 388         return  fNlowPt == decision.fNlowPt
 389                 and fNhighPt == decision.fNhighPt
 390                 and fNunlikeAnyPt == decision.fNunlikeAnyPt
 391                 and fNunlikeLowPt == decision.fNunlikeLowPt
 392                 and fNunlikeHighPt == decision.fNunlikeHighPt
 393                 and fNlikeAnyPt == decision.fNlikeAnyPt
 394                 and fNlikeLowPt == decision.fNlikeLowPt
 395                 and fNlikeHighPt == decision.fNlikeHighPt
 396                 and fNmassAny == decision.fNmassAny
 397                 and fNmassLow == decision.fNmassLow
 398                 and fNmassHigh == decision.fNmassHigh;
 399 }
 400
 401
 402 std::ostream& operator << (
 403                 std::ostream& stream,
 404                 const AliHLTMUONDecision::AliTrackDecision& decision
 405         )
 406 {
 407 /// Stream operator for std::ostream classes.
 408 /// \param stream  The output stream object being written to.
 409 /// \param track  The dHLT decision object to print to the stream.
 410 /// \returns  Returns 'stream'.
 411
 412         stream  << "Passed low/high pT cut: [" << (decision.fPassedLowCut ? "yes" : "no")
 413                 << ", " << (decision.fPassedHighCut ? "yes" : "no")
 414                 << "]; with pT = " << decision.fPt;
 415         return stream;
 416 }
 417
 418
 419 void AliHLTMUONDecision::AliTrackDecision::Print(Option_t* option) const
 420 {
 421 /// Prints the trigger decision to standard output (screen).
 422 /// \param option  Can be one of the following:
 423 ///      - "compact" - prints in a compact format.
 424 ///      - "detail" - prints trigger information in a more detailed format.
 425 ///      - "all" - prints a full dump of the trigger object.
 426
 427         using namespace std;
 428
 429         if (option == NULL or strcmp(option, "") == 0 or
 430             strcmp(option, "compact") == 0
 431            )
 432         {
 433                 cout << *this << endl;
 434         }
 435         else if (strcmp(option, "detail") == 0)
 436         {
 437                 cout << "Trigger decision for track: " << fTrack->Id() << endl;
 438                 cout << "pT = " << fPt << " GeV/c" << endl;
 439                 cout << "pT cut | passed" << endl;
 440                 cout << "-------+--------" << endl;
 441                 cout << "  low  |  " << (fPassedLowCut ? "yes" : "no") << endl;
 442                 cout << " high  |  " << (fPassedHighCut ? "yes" : "no") << endl;
 443         }
 444         else if (strcmp(option, "all") == 0)
 445         {
 446                 cout << "Trigger decision for track: " << fTrack->Id() << endl;
 447                 cout << "pT = " << fPt << " GeV/c" << endl;
 448                 cout << "pT cut | passed" << endl;
 449                 cout << "-------+--------" << endl;
 450                 cout << "  low  |  " << (fPassedLowCut ? "yes" : "no") << endl;
 451                 cout << " high  |  " << (fPassedHighCut ? "yes" : "no") << endl;
 452                 cout << "===== Track details =====" << endl;
 453                 fTrack->Print("all");
 454         }
 455         else
 456         {
 457                 AliError("Unknown option specified. Can only be one of 'compact',"
 458                         " 'detail' or 'all'."
 459                 );
 460         }
 461 }
 462
 463
 464 Int_t AliHLTMUONDecision::AliTrackDecision::Compare(const TObject* obj) const
 465 {
 466 /// We compare this object with 'obj' field by field.
 467 /// \param obj  This is the object to compare to. It must be of type
 468 ///      AliHLTMUONDecision::AliTrackDecision.
 469 /// \returns  -1 if 'this' is smaller than 'obj', 1 if greater and zero if both
 470 ///      objects are the same.
 471
 472         if (obj->IsA() == AliHLTMUONDecision::Class())
 473         {
 474                 const AliHLTMUONDecision::AliTrackDecision* d =
 475                         static_cast<const AliHLTMUONDecision::AliTrackDecision*>(obj);
 476                 if (fPt < d->fPt) return -1;
 477                 if (fPt > d->fPt) return 1;
 478                 if (fPassedLowCut < d->fPassedLowCut) return -1;
 479                 if (fPassedLowCut > d->fPassedLowCut) return 1;
 480                 if (fPassedHighCut < d->fPassedHighCut) return -1;
 481                 if (fPassedHighCut > d->fPassedHighCut) return 1;
 482                 return fTrack->Compare(d->fTrack);
 483         }
 484         else
 485         {
 486                 AliError(Form("Do not know how to compare %s to %s.",
 487                         this->ClassName(),
 488                         obj->ClassName()
 489                 ));
 490                 return -999;
 491         }
 492 }
 493
 494
 495 std::ostream& operator << (
 496                 std::ostream& stream,
 497                 const AliHLTMUONDecision::AliPairDecision& decision
 498         )
 499 {
 500 /// Stream operator for std::ostream classes.
 501 /// \param stream  The output stream object being written to.
 502 /// \param track  The dHLT decision object to print to the stream.
 503 /// \returns  Returns 'stream'.
 504
 505         stream  << (decision.fUnlike ? "Unlike" : "Like")
 506                 << " sign pair passed low/high mass cut: [" << (decision.fPassedLowCut ? "yes" : "no")
 507                 << ", " << (decision.fPassedHighCut ? "yes" : "no")
 508                 << "]; with mass = " << decision.fMass;
 509         return stream;
 510 }
 511
 512
 513 void AliHLTMUONDecision::AliPairDecision::Print(Option_t* option) const
 514 {
 515 /// Prints the trigger decision to standard output (screen).
 516 /// \param option  Can be one of the following:
 517 ///      - "compact" - prints in a compact format.
 518 ///      - "detail" - prints trigger information in a more detailed format.
 519 ///      - "all" - prints a full dump of the trigger object.
 520
 521         using namespace std;
 522
 523         if (option == NULL or strcmp(option, "") == 0 or
 524             strcmp(option, "compact") == 0
 525            )
 526         {
 527                 cout << *this << endl;
 528         }
 529         else if (strcmp(option, "detail") == 0)
 530         {
 531                 cout    << "Trigger decision for track pair: {" << fTrackA->Id()
 532                         << ", " << fTrackB->Id() << "}" << endl;
 533                 cout << "Invariant mass = " << fMass << " GeV/c^2" << endl;
 534                 cout << "mass cut | passed" << endl;
 535                 cout << "---------+--------" << endl;
 536                 cout << "   low   |  " << (fPassedLowCut ? "yes" : "no") << endl;
 537                 cout << "  high   |  " << (fPassedHighCut ? "yes" : "no") << endl;
 538                 cout << "Number of tracks in pair that passed," << endl;
 539                 cout << "   low pT cut = " << Int_t(fLowPtCount) << endl;
 540                 cout << "  high pT cut = " << Int_t(fHighPtCount) << endl;
 541         }
 542         else if (strcmp(option, "all") == 0)
 543         {
 544                 cout    << "Trigger decision for track pair: {" << fTrackA->Id()
 545                         << ", " << fTrackB->Id() << "}" << endl;
 546                 cout << "Invariant mass = " << fMass << " GeV/c^2" << endl;
 547                 cout << "mass cut | passed" << endl;
 548                 cout << "---------+--------" << endl;
 549                 cout << "   low   |  " << (fPassedLowCut ? "yes" : "no") << endl;
 550                 cout << "  high   |  " << (fPassedHighCut ? "yes" : "no") << endl;
 551                 cout << "Number of tracks in pair that passed," << endl;
 552                 cout << "   low pT cut = " << Int_t(fLowPtCount) << endl;
 553                 cout << "  high pT cut = " << Int_t(fHighPtCount) << endl;
 554                 cout << "===== First track details =====" << endl;
 555                 fTrackA->Print("all");
 556                 cout << "===== Second track details =====" << endl;
 557                 fTrackB->Print("all");
 558         }
 559         else
 560         {
 561                 AliError("Unknown option specified. Can only be one of 'compact',"
 562                         " 'detail' or 'all'."
 563                 );
 564         }
 565 }
 566
 567
 568 Int_t AliHLTMUONDecision::AliPairDecision::Compare(const TObject* obj) const
 569 {
 570 /// We compare this object with 'obj' field by field.
 571 /// \param obj  This is the object to compare to. It must be of type
 572 ///      AliHLTMUONDecision::AliPairDecision.
 573 /// \returns  -1 if 'this' is smaller than 'obj', 1 if greater and zero if both
 574 ///      objects are the same.
 575
 576         if (obj->IsA() == AliHLTMUONDecision::Class())
 577         {
 578                 const AliHLTMUONDecision::AliPairDecision* d =
 579                         static_cast<const AliHLTMUONDecision::AliPairDecision*>(obj);
 580                 if (fMass < d->fMass) return -1;
 581                 if (fMass > d->fMass) return 1;
 582                 if (fPassedLowCut < d->fPassedLowCut) return -1;
 583                 if (fPassedLowCut > d->fPassedLowCut) return 1;
 584                 if (fPassedHighCut < d->fPassedHighCut) return -1;
 585                 if (fPassedHighCut > d->fPassedHighCut) return 1;
 586                 if (fUnlike < d->fUnlike) return -1;
 587                 if (fUnlike > d->fUnlike) return 1;
 588                 if (fLowPtCount < d->fLowPtCount) return -1;
 589                 if (fLowPtCount > d->fLowPtCount) return 1;
 590                 if (fHighPtCount < d->fHighPtCount) return -1;
 591                 if (fHighPtCount > d->fHighPtCount) return 1;
 592                 Int_t result = fTrackA->Compare(d->fTrackA);
 593                 if (result != 0) return result;
 594                 return fTrackB->Compare(d->fTrackB);
 595         }
 596         else
 597         {
 598                 AliError(Form("Do not know how to compare %s to %s.",
 599                         this->ClassName(),
 600                         obj->ClassName()
 601                 ));
 602                 return -999;
 603         }
 604 }
 605