ROOT logo
/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/


///////////////////////////////////////////////////////////////////////////
/*

Origin: marian.ivanov@cern.ch
Container classes with MC infomation

*/

#if !defined(__CINT__) || defined(__MAKECINT__)
#include <stdio.h>
#include <string.h>
//ROOT includes
#include "TROOT.h"
#include "Rtypes.h"
#include "TFile.h"
#include "TTree.h"
#include "TChain.h"
#include "TCut.h"
#include "TString.h"
#include "TStopwatch.h"
#include "TParticle.h"
#include "TDatabasePDG.h"
#include "TSystem.h"
#include "TCanvas.h"
#include "TPolyLine3D.h"

//ALIROOT includes
#include "AliRun.h"
#include "AliStack.h"
#include "AliSimDigits.h"
#include "AliTPCParam.h"
#include "AliTPC.h"
#include "AliTPCLoader.h"
#include "AliDetector.h"
#include "AliTrackReference.h"
#include "AliTPCParamSR.h"
#include "AliTracker.h"
#include "AliMagF.h"
#include "AliHelix.h"
#include "AliTrackPointArray.h"

#endif
#include "AliGenKinkInfo.h"

//
// 

ClassImp(AliGenKinkInfo)



/////////////////////////////////////////////////////////////////////////////////

AliGenKinkInfo::AliGenKinkInfo():
  fMCd(),       //info about daughter particle - second particle for V0
  fMCm(),       //info about mother particle   - first particle for V0
  fMCDist1(0),  //info about closest distance according closest MC - linear DCA
  fMCDist2(0),  //info about closest distance parabolic DCA
  fMCRr(0),     // rec position of the vertex 
  fMCR(0)       //exact r position of the vertex
{
  //
  // default constructor
  //
  for (Int_t i=0;i<3;i++){
    fMCPdr[i]=0;
    fMCPd[i]=0;
    fMCX[i]=0;
    fMCPm[i]=0;
    fMCAngle[i]=0;
    fMCXr[i] = 0;
  }
  for (Int_t i=0; i<2; i++) {
    fPdg[i]= 0;
    fLab[i]=0;
  }
}

void AliGenKinkInfo::Update()
{
  //
  // Update information
  //  some redundancy - faster acces to the values in analysis code
  //
  fMCPd[0] = fMCd.GetParticle().Px();
  fMCPd[1] = fMCd.GetParticle().Py();
  fMCPd[2] = fMCd.GetParticle().Pz();
  fMCPd[3] = fMCd.GetParticle().P();
  //
  fMCX[0]  = fMCd.GetParticle().Vx();
  fMCX[1]  = fMCd.GetParticle().Vy();
  fMCX[2]  = fMCd.GetParticle().Vz();
  fMCR       = TMath::Sqrt( fMCX[0]*fMCX[0]+fMCX[1]*fMCX[1]);
  //
  fPdg[0]    = fMCd.GetParticle().GetPdgCode();
  fPdg[1]    = fMCm.GetParticle().GetPdgCode();
  //
  fLab[0]    = fMCd.GetParticle().GetUniqueID();
  fLab[1]    = fMCm.GetParticle().GetUniqueID();
  //
  //
  //
  Double_t x1[3],p1[3];
  TParticle& pdaughter = fMCd.GetParticle();
  x1[0] = pdaughter.Vx();      
  x1[1] = pdaughter.Vy();
  x1[2] = pdaughter.Vz();
  p1[0] = pdaughter.Px();
  p1[1] = pdaughter.Py();
  p1[2] = pdaughter.Pz();
  Double_t sign = 0;
  sign = (pdaughter.GetPDG()->Charge()>0)? -1:1;
  AliHelix dhelix1(x1,p1,sign);
  //
  //
  Double_t x2[3],p2[3];            
  //
  TParticle& pmother = fMCm.GetParticle();
  x2[0] = pmother.Vx();      
  x2[1] = pmother.Vy();      
  x2[2] = pmother.Vz();      
  p2[0] = pmother.Px();
  p2[1] = pmother.Py();
  p2[2] = pmother.Pz();
  //
  const AliTrackReference & pdecay = fMCm.GetTRdecay();
  x2[0] = pdecay.X();      
  x2[1] = pdecay.Y();      
  x2[2] = pdecay.Z();      
  p2[0] = pdecay.Px();
  p2[1] = pdecay.Py();
  p2[2] = pdecay.Pz();
  //
  sign = (pmother.GetPDG()->Charge()>0) ? -1:1;
  AliHelix mhelix(x2,p2,sign);
  
  //
  //
  //
  //find intersection linear
  //
  Double_t distance1, distance2;
  Double_t phase[2][2] = { {0,0},{0,0} };
  Double_t radius[2] = {0};
  Int_t  points = dhelix1.GetRPHIintersections(mhelix, phase, radius);
  Double_t delta1=10000,delta2=10000;    
  if (points>0){
    dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
    dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
    dhelix1.LinearDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
  }
  if (points==2){    
    dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
    dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
    dhelix1.LinearDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
  }
  distance1 = TMath::Min(delta1,delta2);
  //
  //find intersection parabolic
  //
  points = dhelix1.GetRPHIintersections(mhelix, phase, radius);
  delta1=10000,delta2=10000;  
  
  if (points>0){
    dhelix1.ParabolicDCA(mhelix,phase[0][0],phase[0][1],radius[0],delta1);
  }
  if (points==2){    
    dhelix1.ParabolicDCA(mhelix,phase[1][0],phase[1][1],radius[1],delta2);
  }
  
  distance2 = TMath::Min(delta1,delta2);
  //
  if (delta1<delta2){
    //get V0 info
    dhelix1.Evaluate(phase[0][0],fMCXr);
    dhelix1.GetMomentum(phase[0][0],fMCPdr);
    mhelix.GetMomentum(phase[0][1],fMCPm);
    dhelix1.GetAngle(phase[0][0],mhelix,phase[0][1],fMCAngle);
    fMCRr = TMath::Sqrt(radius[0]);
  }
  else{
    dhelix1.Evaluate(phase[1][0],fMCXr);
    dhelix1.GetMomentum(phase[1][0],fMCPdr);
    mhelix.GetMomentum(phase[1][1],fMCPm);
    dhelix1.GetAngle(phase[1][0],mhelix,phase[1][1],fMCAngle);
    fMCRr = TMath::Sqrt(radius[1]);
  }     
  //            
  //
  fMCDist1 = TMath::Sqrt(distance1);
  fMCDist2 = TMath::Sqrt(distance2);      
    
}


Float_t AliGenKinkInfo::GetQt(){
  //
  //
  Float_t momentumd = TMath::Sqrt(fMCPd[0]*fMCPd[0]+fMCPd[1]*fMCPd[1]+fMCPd[2]*fMCPd[2]);
  return TMath::Sin(fMCAngle[2])*momentumd;
}


 AliGenKinkInfo.cxx:1
 AliGenKinkInfo.cxx:2
 AliGenKinkInfo.cxx:3
 AliGenKinkInfo.cxx:4
 AliGenKinkInfo.cxx:5
 AliGenKinkInfo.cxx:6
 AliGenKinkInfo.cxx:7
 AliGenKinkInfo.cxx:8
 AliGenKinkInfo.cxx:9
 AliGenKinkInfo.cxx:10
 AliGenKinkInfo.cxx:11
 AliGenKinkInfo.cxx:12
 AliGenKinkInfo.cxx:13
 AliGenKinkInfo.cxx:14
 AliGenKinkInfo.cxx:15
 AliGenKinkInfo.cxx:16
 AliGenKinkInfo.cxx:17
 AliGenKinkInfo.cxx:18
 AliGenKinkInfo.cxx:19
 AliGenKinkInfo.cxx:20
 AliGenKinkInfo.cxx:21
 AliGenKinkInfo.cxx:22
 AliGenKinkInfo.cxx:23
 AliGenKinkInfo.cxx:24
 AliGenKinkInfo.cxx:25
 AliGenKinkInfo.cxx:26
 AliGenKinkInfo.cxx:27
 AliGenKinkInfo.cxx:28
 AliGenKinkInfo.cxx:29
 AliGenKinkInfo.cxx:30
 AliGenKinkInfo.cxx:31
 AliGenKinkInfo.cxx:32
 AliGenKinkInfo.cxx:33
 AliGenKinkInfo.cxx:34
 AliGenKinkInfo.cxx:35
 AliGenKinkInfo.cxx:36
 AliGenKinkInfo.cxx:37
 AliGenKinkInfo.cxx:38
 AliGenKinkInfo.cxx:39
 AliGenKinkInfo.cxx:40
 AliGenKinkInfo.cxx:41
 AliGenKinkInfo.cxx:42
 AliGenKinkInfo.cxx:43
 AliGenKinkInfo.cxx:44
 AliGenKinkInfo.cxx:45
 AliGenKinkInfo.cxx:46
 AliGenKinkInfo.cxx:47
 AliGenKinkInfo.cxx:48
 AliGenKinkInfo.cxx:49
 AliGenKinkInfo.cxx:50
 AliGenKinkInfo.cxx:51
 AliGenKinkInfo.cxx:52
 AliGenKinkInfo.cxx:53
 AliGenKinkInfo.cxx:54
 AliGenKinkInfo.cxx:55
 AliGenKinkInfo.cxx:56
 AliGenKinkInfo.cxx:57
 AliGenKinkInfo.cxx:58
 AliGenKinkInfo.cxx:59
 AliGenKinkInfo.cxx:60
 AliGenKinkInfo.cxx:61
 AliGenKinkInfo.cxx:62
 AliGenKinkInfo.cxx:63
 AliGenKinkInfo.cxx:64
 AliGenKinkInfo.cxx:65
 AliGenKinkInfo.cxx:66
 AliGenKinkInfo.cxx:67
 AliGenKinkInfo.cxx:68
 AliGenKinkInfo.cxx:69
 AliGenKinkInfo.cxx:70
 AliGenKinkInfo.cxx:71
 AliGenKinkInfo.cxx:72
 AliGenKinkInfo.cxx:73
 AliGenKinkInfo.cxx:74
 AliGenKinkInfo.cxx:75
 AliGenKinkInfo.cxx:76
 AliGenKinkInfo.cxx:77
 AliGenKinkInfo.cxx:78
 AliGenKinkInfo.cxx:79
 AliGenKinkInfo.cxx:80
 AliGenKinkInfo.cxx:81
 AliGenKinkInfo.cxx:82
 AliGenKinkInfo.cxx:83
 AliGenKinkInfo.cxx:84
 AliGenKinkInfo.cxx:85
 AliGenKinkInfo.cxx:86
 AliGenKinkInfo.cxx:87
 AliGenKinkInfo.cxx:88
 AliGenKinkInfo.cxx:89
 AliGenKinkInfo.cxx:90
 AliGenKinkInfo.cxx:91
 AliGenKinkInfo.cxx:92
 AliGenKinkInfo.cxx:93
 AliGenKinkInfo.cxx:94
 AliGenKinkInfo.cxx:95
 AliGenKinkInfo.cxx:96
 AliGenKinkInfo.cxx:97
 AliGenKinkInfo.cxx:98
 AliGenKinkInfo.cxx:99
 AliGenKinkInfo.cxx:100
 AliGenKinkInfo.cxx:101
 AliGenKinkInfo.cxx:102
 AliGenKinkInfo.cxx:103
 AliGenKinkInfo.cxx:104
 AliGenKinkInfo.cxx:105
 AliGenKinkInfo.cxx:106
 AliGenKinkInfo.cxx:107
 AliGenKinkInfo.cxx:108
 AliGenKinkInfo.cxx:109
 AliGenKinkInfo.cxx:110
 AliGenKinkInfo.cxx:111
 AliGenKinkInfo.cxx:112
 AliGenKinkInfo.cxx:113
 AliGenKinkInfo.cxx:114
 AliGenKinkInfo.cxx:115
 AliGenKinkInfo.cxx:116
 AliGenKinkInfo.cxx:117
 AliGenKinkInfo.cxx:118
 AliGenKinkInfo.cxx:119
 AliGenKinkInfo.cxx:120
 AliGenKinkInfo.cxx:121
 AliGenKinkInfo.cxx:122
 AliGenKinkInfo.cxx:123
 AliGenKinkInfo.cxx:124
 AliGenKinkInfo.cxx:125
 AliGenKinkInfo.cxx:126
 AliGenKinkInfo.cxx:127
 AliGenKinkInfo.cxx:128
 AliGenKinkInfo.cxx:129
 AliGenKinkInfo.cxx:130
 AliGenKinkInfo.cxx:131
 AliGenKinkInfo.cxx:132
 AliGenKinkInfo.cxx:133
 AliGenKinkInfo.cxx:134
 AliGenKinkInfo.cxx:135
 AliGenKinkInfo.cxx:136
 AliGenKinkInfo.cxx:137
 AliGenKinkInfo.cxx:138
 AliGenKinkInfo.cxx:139
 AliGenKinkInfo.cxx:140
 AliGenKinkInfo.cxx:141
 AliGenKinkInfo.cxx:142
 AliGenKinkInfo.cxx:143
 AliGenKinkInfo.cxx:144
 AliGenKinkInfo.cxx:145
 AliGenKinkInfo.cxx:146
 AliGenKinkInfo.cxx:147
 AliGenKinkInfo.cxx:148
 AliGenKinkInfo.cxx:149
 AliGenKinkInfo.cxx:150
 AliGenKinkInfo.cxx:151
 AliGenKinkInfo.cxx:152
 AliGenKinkInfo.cxx:153
 AliGenKinkInfo.cxx:154
 AliGenKinkInfo.cxx:155
 AliGenKinkInfo.cxx:156
 AliGenKinkInfo.cxx:157
 AliGenKinkInfo.cxx:158
 AliGenKinkInfo.cxx:159
 AliGenKinkInfo.cxx:160
 AliGenKinkInfo.cxx:161
 AliGenKinkInfo.cxx:162
 AliGenKinkInfo.cxx:163
 AliGenKinkInfo.cxx:164
 AliGenKinkInfo.cxx:165
 AliGenKinkInfo.cxx:166
 AliGenKinkInfo.cxx:167
 AliGenKinkInfo.cxx:168
 AliGenKinkInfo.cxx:169
 AliGenKinkInfo.cxx:170
 AliGenKinkInfo.cxx:171
 AliGenKinkInfo.cxx:172
 AliGenKinkInfo.cxx:173
 AliGenKinkInfo.cxx:174
 AliGenKinkInfo.cxx:175
 AliGenKinkInfo.cxx:176
 AliGenKinkInfo.cxx:177
 AliGenKinkInfo.cxx:178
 AliGenKinkInfo.cxx:179
 AliGenKinkInfo.cxx:180
 AliGenKinkInfo.cxx:181
 AliGenKinkInfo.cxx:182
 AliGenKinkInfo.cxx:183
 AliGenKinkInfo.cxx:184
 AliGenKinkInfo.cxx:185
 AliGenKinkInfo.cxx:186
 AliGenKinkInfo.cxx:187
 AliGenKinkInfo.cxx:188
 AliGenKinkInfo.cxx:189
 AliGenKinkInfo.cxx:190
 AliGenKinkInfo.cxx:191
 AliGenKinkInfo.cxx:192
 AliGenKinkInfo.cxx:193
 AliGenKinkInfo.cxx:194
 AliGenKinkInfo.cxx:195
 AliGenKinkInfo.cxx:196
 AliGenKinkInfo.cxx:197
 AliGenKinkInfo.cxx:198
 AliGenKinkInfo.cxx:199
 AliGenKinkInfo.cxx:200
 AliGenKinkInfo.cxx:201
 AliGenKinkInfo.cxx:202
 AliGenKinkInfo.cxx:203
 AliGenKinkInfo.cxx:204
 AliGenKinkInfo.cxx:205
 AliGenKinkInfo.cxx:206
 AliGenKinkInfo.cxx:207
 AliGenKinkInfo.cxx:208
 AliGenKinkInfo.cxx:209
 AliGenKinkInfo.cxx:210
 AliGenKinkInfo.cxx:211
 AliGenKinkInfo.cxx:212
 AliGenKinkInfo.cxx:213
 AliGenKinkInfo.cxx:214
 AliGenKinkInfo.cxx:215
 AliGenKinkInfo.cxx:216
 AliGenKinkInfo.cxx:217
 AliGenKinkInfo.cxx:218
 AliGenKinkInfo.cxx:219