MUON/AliMUONTriggerTrackToTrackerClusters.cxx

   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 // $Id$
  17
  18 #include "AliMUONTriggerTrackToTrackerClusters.h"
  19
  20 ///\class AliMUONTriggerTrackToTrackerClusters
  21 ///
  22 /// Class to convert trigger tracks into "fake" clusters in stations 4 and 5
  23 ///
  24 /// Only intent is to be able to reconstruct data where stations 4 and 5 were
  25 /// not functionning, typically early cosmic runs
  26 ///
  27 ///\author Laurent Aphecetche, Subatech
  28
  29 #include "AliCodeTimer.h"
  30 #include "AliLog.h"
  31 #include "AliMUONConstants.h"
  32 #include "AliMUONGeometryTransformer.h"
  33 #include "AliMUONTrack.h"
  34 #include "AliMUONTrackExtrap.h"
  35 #include "AliMUONTriggerTrack.h"
  36 #include "AliMUONVCluster.h"
  37 #include "AliMUONVClusterStore.h"
  38 #include "AliMUONVTriggerTrackStore.h"
  39 #include "AliMUONTrackParam.h"
  40 #include "AliMpArea.h"
  41 #include "AliMpDEManager.h"
  42 #include <TMath.h>
  43
  44 ///\cond CLASSIMP
  45 ClassImp(AliMUONTriggerTrackToTrackerClusters)
  46 ///\endcond
  47
  48 //_____________________________________________________________________________
  49 AliMUONTriggerTrackToTrackerClusters::AliMUONTriggerTrackToTrackerClusters(const AliMUONGeometryTransformer& transformer,
  50                                                                            AliMUONVTriggerTrackStore* trackStore)
  51 : TObject(), fkTransformer(transformer), fTriggerTrackStore(trackStore)
  52 {
  53   /// ctor. We do not take ownership of the trigger track store.
  54 }
  55
  56 //_____________________________________________________________________________
  57 AliMUONTriggerTrackToTrackerClusters::~AliMUONTriggerTrackToTrackerClusters()
  58 {
  59   /// dtor
  60 }
  61
  62 //_____________________________________________________________________________
  63 Int_t
  64 AliMUONTriggerTrackToTrackerClusters::DetElemId(Int_t chamber, Double_t x, Double_t y,
  65                                                 Double_t ex, Double_t ey,
  66                                                 Double_t& z) const
  67 {
  68   /// Find in which detection element (x,y) (global) position is.
  69
  70   AliMpDEIterator it;
  71
  72   AliMpArea a(x, y, ex, ey);
  73
  74   it.First(chamber);
  75
  76   while ( !it.IsDone() )
  77   {
  78     Int_t detElemId = it.CurrentDEId();
  79
  80     AliMpArea* area = fkTransformer.GetDEArea(detElemId);
  81
  82     if ( area->Overlap(a) )
  83     {
  84       // get z of the center of that DE.
  85       Double_t dummyx, dummyy;
  86       fkTransformer.Local2Global(detElemId,0,0,0,dummyx,dummyy,z);
  87       return detElemId;
  88     }
  89     it.Next();
  90   }
  91
  92   return -1;
  93 }
  94
  95 //_____________________________________________________________________________
  96 Int_t
  97 AliMUONTriggerTrackToTrackerClusters::GenerateClusters(Int_t iChamber,
  98                                                        AliMUONVClusterStore& clusterStore) const
  99 {
 100   /// Generate clusters in given chamber
 101   /// Return the number of clusters added to the clusterStore
 102
 103   AliCodeTimerAuto(Form("Chamber %d",iChamber));
 104
 105   TIter next(fTriggerTrackStore->CreateIterator());
 106
 107   AliMUONTriggerTrack* track;
 108   Int_t nadded(0);
 109
 110   while ( ( track = static_cast<AliMUONTriggerTrack*>(next()) ) )
 111   {
 112     nadded += GenerateClusters(iChamber,*track,clusterStore);
 113   }
 114   return nadded;
 115 }
 116
 117 //_____________________________________________________________________________
 118 Int_t
 119 AliMUONTriggerTrackToTrackerClusters::GenerateClusters(Int_t iChamber,
 120                                                        const AliMUONTriggerTrack& track,
 121                                                        AliMUONVClusterStore& clusterStore) const
 122 {
 123   /// From a trigger track, generate 1 cluster in given chamber
 124
 125   /// Get a (rough) guestimate of the track momentum
 126
 127   Int_t nadded(0);
 128
 129   Double_t z = AliMUONConstants::DefaultChamberZ(10);
 130
 131   Double_t bendingCoord = track.GetY11();
 132   Double_t bendingSlope = TMath::Tan(track.GetThetay());
 133
 134   Double_t bendingImpact = bendingCoord - z * bendingSlope;
 135
 136   AliDebug(1,Form("TriggerTrack impact parameter=%e",bendingImpact));
 137
 138   //  StdoutToAliDebug(1,track.Print());
 139
 140   Double_t inverseBendingMomentum = 1. / AliMUONTrackExtrap::GetBendingMomentumFromImpactParam(bendingImpact);
 141
 142   // Construct an AliMUONTrackParam from the trigger track, in order to be able to extrapolate it
 143   // to chambers 6..9 planes.
 144
 145   AliMUONTrackParam trackParam;
 146
 147   trackParam.SetZ(z);
 148   trackParam.SetNonBendingCoor(track.GetX11());
 149   trackParam.SetNonBendingSlope(TMath::Tan(track.GetThetax()));
 150   trackParam.SetBendingCoor(bendingCoord);
 151   trackParam.SetBendingSlope(bendingSlope);
 152   trackParam.SetInverseBendingMomentum(inverseBendingMomentum);
 153
 154   Double_t dZ = TMath::Abs(AliMUONConstants::DefaultChamberZ(12) - AliMUONConstants::DefaultChamberZ(10));
 155
 156   Double_t sigmaX = AliMUONConstants::TriggerNonBendingReso();
 157
 158   Double_t sigmaY = AliMUONConstants::TriggerBendingReso();
 159
 160   // Compute and set track parameters covariances
 161   TMatrixD paramCov(5,5);
 162   paramCov.Zero();
 163
 164   // Non bending plane
 165   paramCov(0,0) = sigmaX*sigmaX;
 166   paramCov(0,1) = -sigmaX/dZ;
 167   paramCov(1,0) = paramCov(0,1);
 168
 169   paramCov(1,1) = 2.0*sigmaX/dZ/dZ;
 170
 171   // Bending plane
 172   paramCov(2,2) = sigmaY*sigmaY;
 173   paramCov(2,3) = -sigmaY/dZ;
 174   paramCov(3,2) = paramCov(2,3);
 175
 176   paramCov(3,3) = 2.0*sigmaY/dZ/dZ;
 177
 178   // Inverse bending momentum (50% error)
 179   paramCov(4,4) = 0.5*inverseBendingMomentum * 0.5*inverseBendingMomentum;
 180
 181   // Set covariances
 182   trackParam.SetCovariances(paramCov);
 183
 184   // Now we extrapolate this trackParam to chambers 6 -> 9
 185
 186   const Float_t kFilterThickness = TMath::Abs(AliMUONConstants::MuonFilterZEnd()-AliMUONConstants::MuonFilterZBeg()); // cm
 187
 188   Int_t nclusters = clusterStore.GetSize();
 189
 190   AliMUONTrackParam tp(trackParam);
 191
 192   Double_t zg = AliMUONConstants::DefaultChamberZ(iChamber);
 193   AliMUONTrackExtrap::ExtrapToZCov(&tp, zg); // Extrap to iChamber
 194
 195   AliMUONTrackExtrap::AddMCSEffect(&tp, kFilterThickness, AliMUONConstants::MuonFilterX0()); // Add MCS effects
 196
 197   AliDebug(1,Form("iChamber=%d",iChamber));
 198
 199   StdoutToAliDebug(1,tp.Print("FULLCOV"););
 200
 201   Double_t x = tp.GetNonBendingCoor();
 202   Double_t y = tp.GetBendingCoor();
 203   const TMatrixD& cov = tp.GetCovariances();
 204   Double_t ex = TMath::Sqrt(cov(0,0));
 205   Double_t ey = TMath::Sqrt(cov(2,2));
 206
 207   Double_t zde;
 208
 209   Int_t detElemId = DetElemId(iChamber,x,y,ex,ey,zde);
 210
 211   AliDebug(1,Form("zg = %e zde = %e",zg,zde));
 212
 213   if ( AliMpDEManager::IsValidDetElemId(detElemId) )
 214   {
 215     AliMUONVCluster* rawCluster = clusterStore.Add(AliMpDEManager::GetChamberId(detElemId), detElemId, nclusters);
 216
 217     ++nclusters;
 218     ++nadded;
 219
 220     rawCluster->SetCharge(100.0);
 221     rawCluster->SetXYZ(x, y, zg);
 222     rawCluster->SetErrXY(ex,ey);
 223   }
 224   else
 225   {
 226     AliWarning(Form("No DE found at xg=%e yg=%e",detElemId,x,y));
 227   }
 228
 229   return nadded;
 230 }