| 7cc34f08 |
1 | |
| 2 | |
| 3 | |
| 4 | |
| 5 | // |
| 6 | // |
| 7 | |
| 8 | // This class estiamtes the material budget of the inner detectors in ALICE based |
| 9 | // on the "upper/lower track matching"-method of the ALICE TPC. |
| 10 | |
| 11 | // |
| 12 | // |
| 13 | // |
| 14 | |
| 15 | // |
| 16 | |
| 17 | |
| 18 | |
| 19 | |
| 20 | // ROOT includes |
| 21 | #include <TChain.h> |
| 22 | #include <TMath.h> |
| 23 | #include <TVectorD.h> |
| 24 | #include <TSystem.h> |
| 25 | #include <TFile.h> |
| 26 | #include "TGeoGlobalMagField.h" |
| 27 | |
| 28 | // ALIROOT includes |
| 29 | #include <TTreeStream.h> |
| 30 | #include <AliAnalysisManager.h> |
| 31 | #include <AliESDInputHandler.h> |
| 32 | #include "AliStack.h" |
| 33 | #include "AliMCEvent.h" |
| 34 | #include "AliMCEventHandler.h" |
| 35 | |
| 36 | #include <AliESD.h> |
| 37 | #include "AliMaterialBudget.h" |
| 38 | #include "AliGenInfoMaker.h" |
| 39 | #include "AliHelix.h" |
| 40 | |
| 41 | // |
| 42 | #include "AliMCInfo.h" |
| 7cc34f08 |
43 | #include "AliESDRecInfo.h" |
| 44 | #include "AliTPCParamSR.h" |
| 45 | #include "AliTracker.h" |
| 46 | #include "AliTPCseed.h" |
| 47 | |
| 48 | // STL includes |
| 49 | #include <iostream> |
| 50 | |
| 51 | using namespace std; |
| 52 | |
| 53 | ClassImp(AliMaterialBudget) |
| 54 | |
| 55 | //________________________________________________________________________ |
| 56 | AliMaterialBudget::AliMaterialBudget() : |
| 57 | AliAnalysisTask(), |
| 58 | fMCinfo(0), //! MC event handler |
| 59 | fESD(0), |
| 60 | fDebugStreamer(0), |
| 61 | fStreamLevel(0), |
| 62 | fDebugLevel(0), |
| 63 | fDebugOutputPath(), |
| 64 | fListHist(0), |
| 65 | fHistMult(0), |
| 66 | fCutMaxD(5), // maximal distance in rfi ditection |
| 67 | fCutMaxDz(40), // maximal distance in z ditection |
| 68 | fCutTheta(0.03), // maximal distan theta |
| 69 | fCutMinDir(-0.99) // direction vector products |
| 70 | { |
| 71 | // |
| 72 | // Default constructor (should not be used) |
| 73 | // |
| 74 | } |
| 75 | |
| 76 | AliMaterialBudget::AliMaterialBudget(const AliMaterialBudget& /*info*/) : |
| 77 | AliAnalysisTask(), |
| 78 | fMCinfo(0), //! MC event handler |
| 79 | fESD(0), |
| 80 | // |
| 81 | fDebugStreamer(0), |
| 82 | fStreamLevel(0), |
| 83 | fDebugLevel(), |
| 84 | fDebugOutputPath(), |
| 85 | fListHist(0), |
| 86 | fHistMult(0), |
| 87 | fCutMaxD(5), // maximal distance in rfi ditection |
| 88 | fCutMaxDz(40), // maximal distance in z ditection |
| 89 | fCutTheta(0.03), // maximal distan theta |
| 90 | fCutMinDir(-0.99) // direction vector products |
| 91 | { |
| 92 | // |
| 93 | // Default constructor |
| 94 | // |
| 95 | } |
| 96 | |
| 97 | |
| 98 | |
| 99 | //________________________________________________________________________ |
| 100 | AliMaterialBudget::AliMaterialBudget(const char *name) : |
| 101 | AliAnalysisTask(name, "AliMaterialBudget"), |
| 102 | fMCinfo(0), //! MC event handler |
| 103 | fESD(0), |
| 104 | fDebugStreamer(0), |
| 105 | fStreamLevel(0), |
| 106 | fDebugLevel(0), |
| 107 | fDebugOutputPath(), |
| 108 | fListHist(0), |
| 109 | fHistMult(0), |
| 110 | fCutMaxD(5), // maximal distance in rfi ditection |
| 111 | fCutMaxDz(40), // maximal distance in z ditection |
| 112 | fCutTheta(0.03), // maximal distan theta |
| 113 | fCutMinDir(-0.99) // direction vector products |
| 114 | { |
| 115 | // |
| 116 | // Normal constructor |
| 117 | // |
| 118 | // Input slot #0 works with a TChain |
| 119 | DefineInput(0, TChain::Class()); |
| 120 | // Output slot #0 writes into a TList |
| 121 | DefineOutput(0, TList::Class()); |
| 122 | // |
| 123 | // |
| 124 | } |
| 125 | |
| 126 | AliMaterialBudget::~AliMaterialBudget(){ |
| 127 | // |
| 128 | // |
| 129 | // |
| 130 | if (fDebugLevel>0) printf("AliMaterialBudget::~AliMaterialBudget\n"); |
| 131 | if (fDebugStreamer) delete fDebugStreamer; |
| 132 | fDebugStreamer=0; |
| 133 | } |
| 134 | |
| 135 | |
| 136 | //________________________________________________________________________ |
| 137 | void AliMaterialBudget::ConnectInputData(Option_t *) |
| 138 | { |
| 139 | // |
| 140 | // Connect the input data |
| 141 | // |
| 142 | if(fDebugLevel>3) |
| 143 | cout << "AnalysisTaskTPCCluster::ConnectInputData()" << endl; |
| 144 | |
| 145 | TTree* tree=dynamic_cast<TTree*>(GetInputData(0)); |
| 146 | if (!tree) { |
| 147 | //Printf("ERROR: Could not read chain from input slot 0"); |
| 148 | } |
| 149 | else { |
| 150 | AliESDInputHandler *esdH = dynamic_cast<AliESDInputHandler*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()); |
| 151 | if (!esdH) { |
| 152 | //Printf("ERROR: Could not get ESDInputHandler"); |
| 153 | } |
| 154 | else { |
| 155 | esdH->SetActiveBranches("ESDfriend"); |
| 156 | fESD = esdH->GetEvent(); |
| 157 | //Printf("*** CONNECTED NEW EVENT ****"); |
| 158 | } |
| 159 | } |
| 160 | AliMCEventHandler* mcinfo = (AliMCEventHandler*) (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()); |
| 161 | mcinfo->SetReadTR(kTRUE); |
| 162 | |
| 163 | fMCinfo = mcinfo->MCEvent(); |
| 164 | |
| 165 | |
| 166 | } |
| 167 | |
| 168 | |
| 169 | |
| 170 | |
| 171 | |
| 172 | |
| 173 | //________________________________________________________________________ |
| 174 | void AliMaterialBudget::CreateOutputObjects() |
| 175 | { |
| 176 | // |
| 177 | // Connect the output objects |
| 178 | // |
| 179 | if(fDebugLevel>3) |
| 180 | cout << "AnalysisTaskTPCCluster::CreateOutputObjects()" << endl; |
| 181 | // |
| 182 | fListHist = new TList(); |
| 183 | fListHist->SetOwner(); |
| 184 | // |
| 185 | fHistMult = new TH1F("HistMult", "Number of Tracks per Event; number of tracks per event; number of tracks",501,-0.5,500.5); |
| 186 | fListHist->Add(fHistMult); |
| 187 | |
| 188 | |
| 189 | } |
| 190 | |
| 191 | |
| 192 | //________________________________________________________________________ |
| 193 | void AliMaterialBudget::Exec(Option_t *) { |
| 194 | // |
| 195 | // Execute analysis for current event |
| 196 | // |
| 197 | |
| 198 | if(fDebugLevel>3) |
| 199 | cout << "AliMaterialBudget::Exec()" << endl; |
| 200 | |
| 201 | fHistMult->Fill(fESD->GetNumberOfTracks()); |
| 202 | //FindPairs(fESD); // nearly everything takes place in find pairs... |
| 203 | |
| 204 | // If MC has been connected |
| 205 | |
| 206 | if (!fMCinfo){ |
| 207 | cout << "Not MC info\n" << endl; |
| 208 | }else{ |
| 209 | ProcessMCInfo(); |
| 210 | //mcinfo->Print(); |
| 211 | //DumpInfo(); |
| 212 | } |
| 213 | // |
| 214 | PostData(0, fListHist); |
| 215 | } |
| 216 | |
| 217 | |
| 218 | |
| 219 | |
| 220 | //________________________________________________________________________ |
| 221 | void AliMaterialBudget::Terminate(Option_t *) { |
| 222 | // |
| 223 | // Terminate loop |
| 224 | // |
| 225 | if(fDebugLevel>3) |
| 226 | printf("AliMaterialBudget: Terminate() \n"); |
| 227 | // |
| 228 | if (fDebugLevel>0) printf("AliMCtrackingTestTask::Terminate\n"); |
| 229 | if (fDebugStreamer) delete fDebugStreamer; |
| 230 | fDebugStreamer = 0; |
| 231 | return; |
| 232 | } |
| 233 | |
| 234 | |
| 235 | |
| 236 | TTreeSRedirector *AliMaterialBudget::GetDebugStreamer(){ |
| 237 | // |
| 238 | // Get Debug streamer |
| 239 | // In case debug streamer not yet initialized and StreamLevel>0 create new one |
| 240 | // |
| 241 | if (fStreamLevel==0) return 0; |
| 242 | if (fDebugStreamer) return fDebugStreamer; |
| 243 | TString dsName; |
| 244 | dsName=GetName(); |
| 245 | dsName+="Debug.root"; |
| 246 | dsName.ReplaceAll(" ",""); |
| 247 | fDebugStreamer = new TTreeSRedirector(dsName.Data()); |
| 248 | return fDebugStreamer; |
| 249 | } |
| 250 | |
| 251 | |
| 252 | |
| 253 | |
| 254 | |
| 255 | |
| 256 | AliExternalTrackParam * AliMaterialBudget::MakeTrack(const AliTrackReference* ref, TParticle*part) |
| 257 | { |
| 258 | // |
| 259 | // Make track out of the track ref |
| 260 | // part - TParticle used to determine chargr |
| 261 | // the covariance matrix - equal 0 - starting from ideal MC position |
| 262 | if (!part->GetPDG()) return 0; |
| 263 | Double_t xyz[3]={ref->X(),ref->Y(),ref->Z()}; |
| 264 | Double_t pxyz[3]={ref->Px(),ref->Py(),ref->Pz()}; |
| c6d45264 |
265 | Int_t charge = TMath::Nint(part->GetPDG()->Charge()/3.); |
| 7cc34f08 |
266 | if (ref->X()*ref->Px()+ref->Y()*ref->Py() <0){ |
| 267 | pxyz[0]*=-1; |
| 268 | pxyz[1]*=-1; |
| 269 | pxyz[2]*=-1; |
| c6d45264 |
270 | charge*=-1; |
| 7cc34f08 |
271 | } |
| 272 | Double_t cv[21]; |
| 273 | for (Int_t i=0; i<21;i++) cv[i]=0; |
| 274 | AliExternalTrackParam * param = new AliExternalTrackParam(xyz,pxyz,cv,charge); |
| 275 | return param; |
| 276 | } |
| 277 | |
| 278 | Bool_t AliMaterialBudget::PropagateToPoint(AliExternalTrackParam *param, Double_t *xyz, Double_t mass, Float_t step){ |
| 279 | // |
| 280 | // Propagate track to point xyz using |
| 281 | // AliTracker::PropagateTo functionality |
| 282 | // |
| 283 | // param - track parameters |
| 284 | // xyz - position to propagate |
| 285 | // mass - particle mass |
| 286 | // step - step to be used |
| 287 | Double_t radius=TMath::Sqrt(xyz[0]*xyz[0]+xyz[1]*xyz[1]); |
| 288 | AliTracker::PropagateTrackToBxByBz(param, radius+step, mass, step, kTRUE,0.99); |
| 289 | AliTracker::PropagateTrackToBxByBz(param, radius+0.5, mass, step*0.1, kTRUE,0.99); |
| 290 | Double_t sxyz[3]={0,0,0}; |
| 291 | AliESDVertex vertex(xyz,sxyz); |
| 292 | Bool_t isOK = param->PropagateToDCA(&vertex,AliTracker::GetBz(),10); |
| 293 | return isOK; |
| 294 | } |
| 295 | |
| 296 | |
| 297 | void AliMaterialBudget::ProcessMCInfo(){ |
| 298 | // |
| 299 | // |
| 300 | // |
| 301 | // |
| 18e588e9 |
302 | |
| 7cc34f08 |
303 | const Double_t kPcut=0.05; |
| 304 | const Double_t kMinDrITS = 2.; // minimal distance between references |
| 305 | const Double_t kMinDrTRD = 8.; // minimal distance between references |
| 306 | const Double_t kMinDrTOF = 10.; // minimal distance between references |
| 307 | // |
| 308 | // |
| 309 | // Process tracks |
| 310 | // |
| 311 | Int_t npart = fMCinfo->GetNumberOfTracks(); |
| 312 | if (npart==0) return; |
| 313 | Double_t vertex[4]={0,0,0,0}; |
| 18e588e9 |
314 | |
| 315 | TClonesArray * trefs = new TClonesArray("AliTrackReference"); |
| 316 | TParticle * particle= new TParticle; |
| 317 | |
| 318 | if(particle && trefs) { |
| 319 | fMCinfo->GetParticleAndTR(0, particle, trefs); |
| 320 | } |
| 321 | |
| 322 | |
| 7cc34f08 |
323 | if (particle){ |
| 324 | vertex[0]=particle->Vx(); |
| 325 | vertex[1]=particle->Vy(); |
| 326 | vertex[2]=particle->Vz(); |
| 327 | vertex[3]=particle->R(); |
| 328 | } |
| 329 | // |
| 330 | // |
| 331 | AliTrackReference dummy,*pdummy= &dummy; |
| 332 | AliExternalTrackParam epdummy,*pepdummy= &epdummy; |
| 333 | Int_t nRefITS =0; |
| 334 | Int_t nRefTPC =0; |
| 335 | Int_t nRefTRD =0; |
| 336 | Int_t nRefTOF =0; |
| 337 | AliTrackReference * refITS0, *refITS1; |
| 338 | AliTrackReference * refTPC0, *refTPC1; |
| 339 | AliTrackReference * refTPCIn0, *refTPCIn1; |
| 340 | AliTrackReference * refTRD0, *refTRD1; |
| 341 | AliTrackReference * refTOF0, *refTOF1; |
| 342 | AliTrackReference *refMinR; |
| 343 | // |
| 18e588e9 |
344 | if(!particle) return; |
| 7cc34f08 |
345 | for (Int_t ipart=0;ipart<npart;ipart++){ |
| 346 | //Int_t status = fMCinfo->GetParticleAndTR(ipart, particle, trefs); |
| 347 | AliMCParticle * pp = (AliMCParticle*) fMCinfo->GetTrack(ipart); |
| 348 | if (!pp) continue; |
| 349 | if (particle->P()<kPcut) continue; |
| 350 | Double_t mass = particle->GetMass(); |
| 351 | |
| 352 | // RESET |
| 353 | nRefITS =0; |
| 354 | nRefTPC =0; |
| 355 | nRefTRD =0; |
| 356 | nRefTOF =0; |
| 357 | refITS0=pdummy; refITS1=pdummy; |
| 358 | refTPC0=pdummy; refTPC1=pdummy; |
| 359 | refTPCIn0=pdummy; refTPCIn1=pdummy; |
| 360 | refTRD0=pdummy; refTRD1=pdummy; |
| 361 | refTOF0=pdummy; refTOF1=pdummy; |
| 362 | refMinR = pdummy; |
| 363 | // |
| 364 | Int_t nref = pp->GetNumberOfTrackReferences(); |
| 365 | if (nref==0) continue; |
| 366 | for (Int_t iref = 0; iref < nref; iref++) { |
| 367 | AliTrackReference *ref = pp->GetTrackReference(iref); |
| 368 | if (ref->DetectorId()==AliTrackReference::kDisappeared) continue; |
| 369 | // if (ref.Px()*particle.Px()+ref.Py()*particle.Py()<0) break; // returning track |
| 370 | // |
| 371 | if (ref->DetectorId()==AliTrackReference::kITS){ |
| 372 | if (TMath::Abs(ref->R()-refITS1->R())>kMinDrITS) { |
| 373 | refITS1 = ref; |
| 374 | nRefITS++; |
| 375 | } |
| 376 | if (refITS0==pdummy) refITS0=ref; |
| 377 | } |
| 378 | if (ref->DetectorId()==AliTrackReference::kTPC){ |
| 379 | nRefTPC++; |
| 380 | refTPC1 = ref; |
| 381 | if (refTPC0==pdummy) refTPC0=ref; |
| 382 | } |
| 383 | if (ref->DetectorId()==AliTrackReference::kTRD){ |
| 384 | if (TMath::Abs(ref->R()-refTRD1->R())>kMinDrTRD) { |
| 385 | refTRD1 = ref; |
| 386 | nRefTRD++; |
| 387 | } |
| 388 | if (refTRD0==pdummy) refTRD0=ref; |
| 389 | } |
| 390 | if (ref->DetectorId()==AliTrackReference::kTOF){ |
| 391 | if (TMath::Abs(ref->X()-refTOF1->X()) + TMath::Abs(ref->Y()-refTOF1->Y()) + TMath::Abs(ref->Z()-refTOF1->Z())>kMinDrTOF) { |
| 392 | refTOF1 = ref; |
| 393 | nRefTOF++; |
| 394 | } |
| 395 | if (refTOF0==pdummy) refTOF0=ref; |
| 396 | } |
| 397 | // |
| 398 | // "find inner track ref" |
| 399 | if (ref->DetectorId()==AliTrackReference::kTPC){ |
| 400 | if (ref->Px()*ref->X()+ref->Py()*ref->Y()<0){ |
| 401 | // track in |
| 402 | if (refTPCIn0 == pdummy) refTPCIn0=ref; |
| 403 | if (refTPCIn0 != pdummy && refTPCIn0->R()>ref->R()) |
| 404 | refTPCIn0=ref; |
| 405 | } |
| 406 | if (ref->Px()*ref->X()+ref->Py()*ref->Y()>0){ |
| 407 | // track in |
| 408 | if (refTPCIn1 == pdummy) refTPCIn1=ref; |
| 409 | if (refTPCIn1 != pdummy && refTPCIn1->R()>ref->R()) |
| 410 | refTPCIn1=ref; |
| 411 | } |
| 412 | } |
| 413 | |
| 414 | |
| 415 | if (refMinR==pdummy && ref->P()>0 ){ |
| 416 | refMinR=ref; |
| 417 | } |
| 418 | if (refMinR->R()>ref->R() && ref->P()>0 ) refMinR=ref; |
| 419 | |
| 420 | } |
| 421 | // |
| 422 | AliExternalTrackParam * trackMC = pepdummy; |
| 423 | //track0->GetDZ(0,0,0,bz,dvertex0) |
| 424 | Float_t dist[2]={0,0}; |
| 425 | AliMagF* field = (AliMagF*) TGeoGlobalMagField::Instance()->GetField(); |
| 426 | Double_t esdfield= fESD->GetMagneticField(); |
| 427 | Double_t xyz[3]={0,0,0}; |
| 428 | Double_t bxyz[3]={0,0,0}; |
| 429 | field->Field(xyz,bxyz); |
| 430 | if (refMinR->P()>0) { |
| 431 | trackMC = MakeTrack(refMinR,particle); |
| 432 | trackMC->GetDZ(0,0,0,bxyz[2],dist); |
| 433 | } |
| 434 | Double_t alphaTOF0 = TMath::ATan2(refTOF0->Y(),refTOF0->X()); |
| 435 | Double_t alphaTOF1 = TMath::ATan2(refTOF1->Y(),refTOF1->X()); |
| 436 | Int_t dsecTOF = TMath::Nint(180*(alphaTOF0-alphaTOF1)/(TMath::Pi()*20.)-0.5); |
| 437 | // |
| 438 | // make the two different TPC tracks and propagate them to their DCA |
| 439 | // |
| 440 | Double_t dP = 0; |
| 441 | Bool_t statusProp = kFALSE; |
| 442 | Double_t dY = 0; |
| 443 | Double_t dZ = 0; |
| 444 | AliExternalTrackParam * track0 = pepdummy; |
| 445 | AliExternalTrackParam * track1 = pepdummy; |
| 446 | AliExternalTrackParam * otrack0 = pepdummy; |
| 447 | AliExternalTrackParam * otrack1 = pepdummy; |
| 448 | if (refTPCIn0!=pdummy && refTPCIn1!=pdummy) { |
| 449 | track0 = MakeTrack(refTPCIn0,particle); |
| 450 | track1 = MakeTrack(refTPCIn1,particle); |
| 451 | otrack0 = MakeTrack(refTPCIn0,particle); |
| 452 | otrack1 = MakeTrack(refTPCIn1,particle); |
| 453 | dP = track0->P() - track1->P(); // momentum loss |
| 454 | statusProp = AliMaterialBudget::PropagateCosmicToDCA(track0,track1,mass); |
| 455 | if (statusProp) { |
| 456 | dY = track0->GetY() - track1->GetY(); |
| 457 | dZ = track0->GetZ() - track1->GetZ(); |
| 458 | } |
| 459 | } |
| 460 | // |
| 461 | TTreeSRedirector *pcstream = GetDebugStreamer(); |
| 462 | if (pcstream){ |
| 463 | char name[100]; |
| 464 | for (Int_t id=0; id<3;id++){ |
| 465 | |
| 466 | // if (id==0) sprintf(name,"mcAll"); // all tracks: inconvenient to cut if on is only interest in tracks which reach the TPC |
| 467 | if (id==0) continue; // require TPC |
| 18e588e9 |
468 | if (id==1) snprintf(name,100,"mcITS"); |
| 469 | if (id==2) snprintf(name,100,"mcTPC"); |
| 7cc34f08 |
470 | if (id==1&& nRefITS==0) continue; |
| 471 | if (id==2&& nRefTPC==0) continue; |
| 472 | |
| 473 | (*pcstream)<<name<< |
| 474 | "ipart="<<ipart<< |
| 475 | "p.="<<particle<< |
| 476 | "mass="<<mass<< |
| 477 | "tbfield="<<bxyz[2]<< |
| 478 | "esdbfield="<<esdfield<< |
| 479 | // counter |
| 480 | "nref="<<nref<< |
| 481 | "nRefITS="<<nRefITS<< |
| 482 | "nRefTPC="<<nRefTPC<< |
| 483 | "nRefTRD="<<nRefTRD<< |
| 484 | "nRefTOF="<<nRefTOF<< |
| 485 | //references |
| 486 | "refMinR.="<<refMinR<< |
| 487 | "refITS0.="<<refITS0<< |
| 488 | "refITS1.="<<refITS1<< |
| 489 | "refTPC0.="<<refTPC0<< |
| 490 | "refTPC1.="<<refTPC1<< |
| 491 | "refTPCIn0.="<<refTPCIn0<< |
| 492 | "refTPCIn1.="<<refTPCIn1<< |
| 493 | "refTRD0.="<<refTRD0<< |
| 494 | "refTRD1.="<<refTRD1<< |
| 495 | "refTOF0.="<<refTOF0<< |
| 496 | "refTOF1.="<<refTOF1<< |
| 497 | //trigger variables |
| 498 | "dsecTOF="<<dsecTOF<< // delta TOF sectors |
| 499 | "aTOF0="<<alphaTOF0<< |
| 500 | "aTOF1="<<alphaTOF1<< |
| 501 | // |
| 502 | // track |
| 503 | "dr="<<dist[0]<< |
| 504 | "dz="<<dist[1]<< |
| 505 | // |
| 506 | // "two" TPC tracks |
| 507 | "status="<<statusProp<< |
| 508 | "dP="<<dP<< |
| 509 | "otrack0.="<<otrack0<< |
| 510 | "otrack1.="<<otrack1<< |
| 511 | "track0.="<<track0<< |
| 512 | "track1.="<<track1<< |
| 513 | "dY="<<dY<< |
| 514 | "dZ="<<dZ<< |
| 515 | "\n"; |
| 516 | } |
| 517 | } |
| 518 | } |
| 519 | } |
| 520 | |
| 521 | |
| 522 | |
| 523 | void AliMaterialBudget::ProcessRefTracker(AliTrackReference* refIn, AliTrackReference* refOut, TParticle*part,Int_t type){ |
| 524 | // |
| 525 | // Test propagation from In to out |
| 526 | // |
| 527 | AliExternalTrackParam *param = 0; |
| 528 | AliExternalTrackParam *paramMC = 0; |
| 529 | Double_t xyzIn[3]={refIn->X(),refIn->Y(), refIn->Z()}; |
| 530 | Double_t mass = part->GetMass(); |
| 531 | Double_t step=1; |
| 532 | // |
| 533 | param=MakeTrack(refOut,part); |
| 534 | paramMC=MakeTrack(refOut,part); |
| 535 | if (!param) return; |
| 536 | if (type<3) PropagateToPoint(param,xyzIn, mass, step); |
| 537 | if (type==3) { |
| 538 | AliTPCseed seed; |
| 539 | seed.Set(param->GetX(),param->GetAlpha(),param->GetParameter(),param->GetCovariance()); |
| 540 | Float_t alpha= TMath::ATan2(refIn->Y(),refIn->X()); |
| 26a35193 |
541 | if(seed.Rotate(alpha-seed.GetAlpha()) == kFALSE) return; |
| 7cc34f08 |
542 | seed.SetMass(mass); |
| 543 | for (Float_t xlayer= seed.GetX(); xlayer>refIn->R(); xlayer-=step){ |
| 544 | seed.PropagateTo(xlayer); |
| 545 | } |
| 546 | seed.PropagateTo(refIn->R()); |
| 547 | param->Set(seed.GetX(),seed.GetAlpha(),seed.GetParameter(),seed.GetCovariance()); |
| 548 | } |
| 549 | TTreeSRedirector *pcstream = GetDebugStreamer(); |
| 550 | TVectorD gpos(3); |
| 551 | TVectorD gmom(3); |
| 552 | param->GetXYZ(gpos.GetMatrixArray()); |
| 553 | param->GetPxPyPz(gmom.GetMatrixArray()); |
| 554 | if (pcstream){ |
| 555 | (*pcstream)<<"MC"<< |
| 556 | "type="<<type<< |
| 557 | "step="<<step<< |
| 558 | "refIn.="<<refIn<< |
| 559 | "refOut.="<<refOut<< |
| 560 | "p.="<<part<< |
| 561 | "par.="<<param<< |
| 562 | "parMC.="<<paramMC<< |
| 563 | "gpos.="<<&gpos<< |
| 564 | "gmom.="<<&gmom<< |
| 565 | "\n"; |
| 566 | } |
| 567 | } |
| 568 | |
| 569 | |
| 570 | void AliMaterialBudget::FitTrackRefs(TParticle * part, TClonesArray * trefs){ |
| 571 | // |
| 572 | // |
| 573 | // |
| 574 | // |
| 575 | const Int_t kMinRefs=6; |
| 576 | Int_t nrefs = trefs->GetEntries(); |
| 577 | if (nrefs<kMinRefs) return; // we should have enough references |
| 578 | Int_t iref0 =-1; |
| 579 | Int_t iref1 =-1; |
| 580 | |
| 581 | for (Int_t iref=0; iref<nrefs; iref++){ |
| 582 | AliTrackReference * ref = (AliTrackReference*)trefs->At(iref); |
| 583 | if (!ref) continue; |
| 584 | Float_t dir = ref->X()*ref->Px()+ref->Y()*ref->Py(); |
| 585 | if (dir<0) break; |
| 586 | if (ref->DetectorId()!=AliTrackReference::kTPC) continue; |
| 587 | if (iref0<0) iref0 = iref; |
| 588 | iref1 = iref; |
| 589 | } |
| 590 | if (iref1-iref0<kMinRefs) return; |
| c6d45264 |
591 | Double_t covar[15]; |
| 592 | for (Int_t icov=0; icov<15; icov++) covar[icov]=0; |
| 7cc34f08 |
593 | covar[0]=1; |
| 594 | covar[2]=1; |
| 595 | covar[5]=1; |
| 596 | covar[9]=1; |
| 597 | covar[14]=1; |
| 598 | |
| 599 | AliTrackReference * refIn = (AliTrackReference*)trefs->At(iref0); |
| 600 | AliTrackReference * refOut = (AliTrackReference*)trefs->At(iref1); |
| 601 | AliExternalTrackParam *paramPropagate= MakeTrack(refIn,part); |
| 602 | AliExternalTrackParam *paramUpdate = MakeTrack(refIn,part); |
| 603 | paramUpdate->AddCovariance(covar); |
| 604 | Double_t mass = part->GetMass(); |
| 605 | Double_t charge = part->GetPDG()->Charge()/3.; |
| 606 | /* |
| 607 | Float_t alphaIn= TMath::ATan2(refIn->Y(),refIn->X()); |
| 608 | Float_t radiusIn= refIn->R(); |
| 609 | Float_t alphaOut= TMath::ATan2(refOut->Y(),refOut->X()); |
| 610 | Float_t radiusOut= refOut->R(); |
| 611 | */ |
| 612 | Bool_t isOKP=kTRUE; |
| 613 | Bool_t isOKU=kTRUE; |
| 614 | AliMagF * field = (AliMagF*) TGeoGlobalMagField::Instance()->GetField(); |
| 615 | for (Int_t iref = iref0; iref<=iref1; iref++){ |
| 616 | AliTrackReference * ref = (AliTrackReference*)trefs->At(iref); |
| 617 | Float_t alphaC= TMath::ATan2(ref->Y(),ref->X()); |
| 618 | Double_t pos[3] = {ref->X(), ref->Y(), ref->Z()}; |
| 619 | Double_t mag[3]; |
| 620 | field->Field(pos,mag); |
| 621 | isOKP&=paramPropagate->Rotate(alphaC); |
| 622 | isOKU&=paramUpdate->Rotate(alphaC); |
| 623 | for (Float_t xref= paramPropagate->GetX(); xref<ref->R(); xref++){ |
| 624 | isOKP&=paramPropagate->PropagateTo(xref, mag[2]); |
| 625 | isOKU&=paramUpdate->PropagateTo(xref, mag[2]); |
| 626 | } |
| 627 | isOKP&=paramPropagate->PropagateTo(ref->R(), mag[2]); |
| 628 | isOKU&=paramUpdate->PropagateTo(ref->R(), mag[2]); |
| 629 | Double_t clpos[2] = {0, ref->Z()}; |
| 630 | Double_t clcov[3] = { 0.005,0,0.005}; |
| 631 | isOKU&= paramUpdate->Update(clpos, clcov); |
| 632 | } |
| 633 | TTreeSRedirector *pcstream = GetDebugStreamer(); |
| 634 | if (pcstream){ |
| 635 | TVectorD gposU(3); |
| 636 | TVectorD gmomU(3); |
| 637 | TVectorD gposP(3); |
| 638 | TVectorD gmomP(3); |
| 639 | paramUpdate->GetXYZ(gposU.GetMatrixArray()); |
| 640 | paramUpdate->GetPxPyPz(gmomU.GetMatrixArray()); |
| 641 | paramPropagate->GetXYZ(gposP.GetMatrixArray()); |
| 642 | paramPropagate->GetPxPyPz(gmomP.GetMatrixArray()); |
| 643 | |
| 644 | (*pcstream)<<"MCupdate"<< |
| 645 | "isOKU="<<isOKU<< |
| 646 | "isOKP="<<isOKP<< |
| 647 | "m="<<mass<< |
| 648 | "q="<<charge<< |
| 649 | "part.="<<part<< |
| 650 | "refIn.="<<refIn<< |
| 651 | "refOut.="<<refOut<< |
| 652 | "pP.="<<paramPropagate<< |
| 653 | "pU.="<<paramUpdate<< |
| 654 | "gposU.="<<&gposU<< |
| 655 | "gmomU.="<<&gmomU<< |
| 656 | "gposP.="<<&gposP<< |
| 657 | "gmomP.="<<&gmomP<< |
| 658 | "\n"; |
| 659 | } |
| 660 | } |
| 661 | |
| 662 | |
| 663 | |
| 664 | void AliMaterialBudget::FindPairs(AliESDEvent * event) { |
| 665 | // |
| 666 | // This function matches the cosmic tracks and calculates the energy loss in the material. |
| 667 | // If accessible the "true" energy loss is determined with the MC track references. |
| 668 | // |
| 669 | // |
| 670 | // Find cosmic pairs |
| 671 | // |
| 672 | // Track0 is choosen in upper TPC part |
| 673 | // Track1 is choosen in lower TPC part |
| 674 | // |
| 675 | if(fDebugLevel>3) |
| 676 | cout << "AliMaterialBudget::FindPairs()" << endl; |
| 677 | |
| 678 | |
| 679 | AliESDfriend *ESDfriend=static_cast<AliESDfriend*>(event->FindListObject("AliESDfriend")); |
| 680 | Int_t ntracks=event->GetNumberOfTracks(); |
| 681 | TObjArray tpcSeeds(ntracks); |
| 682 | if (ntracks==0) return; |
| 683 | Double_t vtxx[3]={0,0,0}; |
| 684 | Double_t svtxx[3]={0.000001,0.000001,100.}; |
| 685 | AliESDVertex vtx(vtxx,svtxx); |
| 686 | // |
| 687 | //track loop |
| 688 | // |
| 689 | for (Int_t i=0;i<ntracks;++i) { |
| 690 | AliESDtrack *track = event->GetTrack(i); |
| 691 | const AliExternalTrackParam * trackIn = track->GetInnerParam(); |
| 692 | const AliExternalTrackParam * trackOut = track->GetOuterParam(); |
| 693 | if (!trackIn) continue; |
| 694 | if (!trackOut) continue; |
| 695 | AliESDfriendTrack *friendTrack = ESDfriend->GetTrack(i); |
| ff7abee2 |
696 | if (!friendTrack) continue; |
| 7cc34f08 |
697 | TObject *calibObject; |
| 698 | AliTPCseed *seed = 0; |
| 699 | for (Int_t l=0;(calibObject=friendTrack->GetCalibObject(l));++l) { |
| 700 | if ((seed=dynamic_cast<AliTPCseed*>(calibObject))) break; |
| 701 | } |
| 702 | if (seed) tpcSeeds.AddAt(seed,i); |
| 703 | } |
| 704 | |
| 705 | if (ntracks<2) return; |
| 706 | // |
| 707 | // Find pairs |
| 708 | // |
| 709 | for (Int_t i=0;i<ntracks;++i) { |
| 710 | AliESDtrack *track0 = event->GetTrack(i); |
| 711 | // track0 - choosen upper part |
| 712 | if (!track0) continue; |
| 713 | if (!track0->GetOuterParam()) continue; |
| 714 | if (track0->GetOuterParam()->GetAlpha()<0) continue; |
| 715 | Double_t dir0[3]; |
| 716 | track0->GetDirection(dir0); |
| 717 | for (Int_t j=0;j<ntracks;++j) { |
| 718 | if (i==j) continue; |
| 719 | AliESDtrack *track1 = event->GetTrack(j); |
| 720 | //track 1 lower part |
| 721 | if (!track1) continue; |
| 722 | if (!track1->GetOuterParam()) continue; |
| 723 | if (track1->GetOuterParam()->GetAlpha()>0) continue; |
| 724 | // |
| 725 | Double_t dir1[3]; |
| 726 | track1->GetDirection(dir1); |
| 727 | |
| 728 | AliTPCseed * seed0 = (AliTPCseed*) tpcSeeds.At(i); |
| 729 | AliTPCseed * seed1 = (AliTPCseed*) tpcSeeds.At(j); |
| 730 | if (! seed0) continue; |
| 731 | if (! seed1) continue; |
| 732 | // |
| 733 | Float_t dir = (dir0[0]*dir1[0] + dir0[1]*dir1[1] + dir0[2]*dir1[2]); |
| 734 | Float_t d0 = track0->GetLinearD(0,0); |
| 735 | Float_t d1 = track1->GetLinearD(0,0); |
| 736 | // |
| 737 | // conservative cuts - convergence to be guarantied |
| 738 | // applying before track propagation |
| 739 | if (TMath::Abs(d0+d1)>fCutMaxD) continue; // distance to the 0,0 |
| 740 | if (dir>fCutMinDir) continue; // direction vector product |
| 741 | Float_t bz = AliTracker::GetBz(); |
| 742 | Float_t dvertex0[2]; //distance to 0,0 |
| 743 | Float_t dvertex1[2]; //distance to 0,0 |
| 744 | track0->GetDZ(0,0,0,bz,dvertex0); |
| 745 | track1->GetDZ(0,0,0,bz,dvertex1); |
| 746 | if (TMath::Abs(dvertex0[1])>250) continue; |
| 747 | if (TMath::Abs(dvertex1[1])>250) continue; |
| 748 | // |
| 749 | // |
| 750 | // |
| 751 | Float_t dmax = TMath::Max(TMath::Abs(d0),TMath::Abs(d1)); |
| 752 | AliExternalTrackParam param0(*track0); |
| 753 | AliExternalTrackParam param1(*track1); |
| 754 | // |
| 755 | // Propagate using Magnetic field and correct fo material budget |
| 756 | // |
| 757 | AliTracker::PropagateTrackToBxByBz(¶m0,dmax+1,0.0005,3,kTRUE); |
| 758 | AliTracker::PropagateTrackToBxByBz(¶m1,dmax+1,0.0005,3,kTRUE); |
| 759 | // |
| 760 | // Propagate rest to the 0,0 DCA - z should be ignored |
| 761 | // |
| 762 | Bool_t b0 = param0.PropagateToDCA(&vtx,bz,1000); |
| 763 | Bool_t b1 = param1.PropagateToDCA(&vtx,bz,1000); |
| 764 | // |
| 765 | param0.GetDZ(0,0,0,bz,dvertex0); |
| 766 | param1.GetDZ(0,0,0,bz,dvertex1); |
| 767 | if (TMath::Abs(param0.GetZ()-param1.GetZ())>fCutMaxDz) continue; |
| 768 | // |
| 769 | Double_t xyz0[3];//,pxyz0[3]; |
| 770 | Double_t xyz1[3];//,pxyz1[3]; |
| 771 | param0.GetXYZ(xyz0); |
| 772 | param1.GetXYZ(xyz1); |
| 773 | Bool_t isPair = IsPair(¶m0,¶m1); |
| 774 | // |
| 775 | // HERE WE WILL PUT THE ACCESS TO THE MC TRACKS AND MATCH THESE !!!! |
| 776 | // |
| 777 | Int_t label0 = TMath::Abs(track0->GetLabel()); |
| 778 | AliMCParticle *mcParticle0 = (AliMCParticle*) fMCinfo->GetTrack(label0); |
| 779 | TParticle *particle0 = mcParticle0->Particle(); |
| 780 | AliTrackReference *ref0 = GetFirstTPCTrackRef(mcParticle0); // get the first TPC track reference |
| 781 | if (!ref0) continue; |
| 782 | AliExternalTrackParam *paramMC0 = 0; |
| 783 | paramMC0 = MakeTrack(ref0, particle0); |
| 784 | // |
| 785 | Int_t label1 = TMath::Abs(track1->GetLabel()); |
| 786 | AliMCParticle *mcParticle1 = (AliMCParticle*) fMCinfo->GetTrack(label1); |
| 787 | TParticle *particle1 = mcParticle1->Particle(); |
| 788 | AliTrackReference *ref1 = GetFirstTPCTrackRef(mcParticle1); // get the first TPC track reference |
| 789 | if (!ref1) continue; |
| 790 | AliExternalTrackParam *paramMC1 = 0; |
| 791 | paramMC1 = MakeTrack(ref1, particle1); |
| 792 | // |
| 793 | // ACCESS TOF INFORMATION |
| 794 | Int_t nTrackRefTOF0 = 0; |
| 795 | Int_t nTrackRefITS0 = 0; |
| 796 | AliTrackReference * refLastTOF0 = 0; |
| 797 | AliTrackReference * refFirstTOF0 = GetAllTOFinfo(mcParticle0, nTrackRefTOF0, nTrackRefITS0); |
| 798 | Float_t alphaTOF0 = 0; |
| 799 | if (refFirstTOF0) alphaTOF0 = refFirstTOF0->Alpha(); |
| 800 | // |
| 801 | Int_t nTrackRefTOF1 = 0; |
| 802 | Int_t nTrackRefITS1 = 0; |
| 803 | AliTrackReference * refLastTOF1 = 0; |
| 804 | AliTrackReference * refFirstTOF1 =GetAllTOFinfo(mcParticle1, nTrackRefTOF1, nTrackRefITS1); |
| 805 | Float_t alphaTOF1 = 0; |
| 806 | if (refFirstTOF1) alphaTOF1 = refFirstTOF1->Alpha(); |
| 807 | //cout << " STATUS "<<nTrackRefTOF0<<" "<<refFirstTOF0<<" "<<refLastTOF0<<" " <<refFirstTOF1<<" "<<refLastTOF1<<endl; |
| 808 | // |
| 809 | // |
| 810 | // |
| 811 | if (fStreamLevel>0){ |
| 812 | TTreeSRedirector * cstream = GetDebugStreamer(); |
| 813 | AliExternalTrackParam *ip0 = (AliExternalTrackParam *)track0->GetInnerParam(); |
| 814 | AliExternalTrackParam *ip1 = (AliExternalTrackParam *)track1->GetInnerParam(); |
| 815 | AliExternalTrackParam *op0 = (AliExternalTrackParam *)track0->GetOuterParam(); |
| 816 | AliExternalTrackParam *op1 = (AliExternalTrackParam *)track1->GetOuterParam(); |
| 817 | // |
| 818 | // |
| 819 | // |
| 820 | if (cstream) { |
| 821 | (*cstream) << "Track0" << |
| 822 | "dir="<<dir<< // direction |
| 823 | "OK="<<isPair<< // will be accepted |
| 824 | "b0="<<b0<< // propagate status |
| 825 | "b1="<<b1<< // propagate status |
| 826 | // |
| 827 | "Particle.="<<particle0<< // TParticle with generated momentum |
| 828 | "NTrackRefTOF0="<<nTrackRefTOF0<< // Number of TOF track references upper |
| 829 | "NTrackRefTOF1="<<nTrackRefTOF1<< // Number of TOF track references lower |
| 830 | "NTrackRefITS0="<<nTrackRefITS0<< // Number of ITS track references upper |
| 831 | "NTrackRefITS1="<<nTrackRefITS1<< // Number of ITS track references lower |
| 832 | "Alpha0="<<alphaTOF0<< // alpha upper |
| 833 | "Alpha1="<<alphaTOF1<< // alpha lower |
| 834 | "RefFirstTOF0.="<<refFirstTOF0<< // first tof reference upper |
| 835 | "RefLastTOF0.="<<refLastTOF0<< // last tof reference upper |
| 836 | "RefFirstTOF1.="<<refFirstTOF1<< // first tof reference lower |
| 837 | "RefLastTOF1.="<<refLastTOF1<< // last tof reference lower |
| 838 | // |
| 839 | "Orig0.=" << track0 << // original track 0 |
| 840 | "Orig1.=" << track1 << // original track 1 |
| 841 | "Tr0.="<<¶m0<< // track propagated to the DCA 0,0 |
| 842 | "Tr1.="<<¶m1<< // track propagated to the DCA 0,0 |
| 843 | "Ip0.="<<ip0<< // inner param - upper |
| 844 | "Ip1.="<<ip1<< // inner param - lower |
| 845 | "Op0.="<<op0<< // outer param - upper |
| 846 | "Op1.="<<op1<< // outer param - lower |
| 847 | // |
| 848 | "paramTrackRef0.="<<paramMC0<< // "ideal" MC true track parameters from track references - upper |
| 849 | "paramTrackRef1.="<<paramMC1<< // "ideal" MC true track parameters from track references - lower |
| 850 | // |
| 851 | "v00="<<dvertex0[0]<< // distance using kalman |
| 852 | "v01="<<dvertex0[1]<< // |
| 853 | "v10="<<dvertex1[0]<< // |
| 854 | "v11="<<dvertex1[1]<< // |
| 855 | "d0="<<d0<< // linear distance to 0,0 |
| 856 | "d1="<<d1<< // linear distance to 0,0 |
| 857 | // |
| 858 | // |
| 859 | // |
| 860 | "x00="<<xyz0[0]<< // global position close to vertex |
| 861 | "x01="<<xyz0[1]<< |
| 862 | "x02="<<xyz0[2]<< |
| 863 | // |
| 864 | "x10="<<xyz1[0]<< // global position close to vertex |
| 865 | "x11="<<xyz1[1]<< |
| 866 | "x12="<<xyz1[2]<< |
| 867 | // |
| 868 | "dir00="<<dir0[0]<< // direction upper |
| 869 | "dir01="<<dir0[1]<< |
| 870 | "dir02="<<dir0[2]<< |
| 871 | // |
| 872 | "dir10="<<dir1[0]<< // direction lower |
| 873 | "dir11="<<dir1[1]<< |
| 874 | "dir12="<<dir1[2]<< |
| 875 | // |
| 876 | // |
| 877 | "Seed0.=" << seed0 << // original seed 0 |
| 878 | "Seed1.=" << seed1 << // original seed 1 |
| 879 | // |
| 880 | "\n"; |
| 881 | } |
| 882 | } |
| 883 | delete paramMC0; |
| 884 | delete paramMC1; |
| 885 | } |
| 886 | } |
| 887 | |
| 888 | return; |
| 889 | |
| 890 | |
| 891 | } |
| 892 | |
| 893 | Bool_t AliMaterialBudget::IsPair(AliExternalTrackParam *tr0, AliExternalTrackParam *tr1){ |
| 894 | // |
| 895 | // |
| 896 | /* |
| 897 | // 0. Same direction - OPOSITE - cutDir +cutT |
| 898 | TCut cutDir("cutDir","dir<-0.99") |
| 899 | // 1. |
| 900 | TCut cutT("cutT","abs(Tr1.fP[3]+Tr0.fP[3])<0.03") |
| 901 | // |
| 902 | // 2. The same rphi |
| 903 | TCut cutD("cutD","abs(Tr0.fP[0]+Tr1.fP[0])<5") |
| 904 | // |
| 905 | // |
| 906 | // |
| 907 | TCut cutPt("cutPt","abs(Tr1.fP[4]+Tr0.fP[4])<1&&abs(Tr0.fP[4])+abs(Tr1.fP[4])<10"); |
| 908 | // 1/Pt diff cut |
| 909 | */ |
| 910 | const Double_t *p0 = tr0->GetParameter(); |
| 911 | const Double_t *p1 = tr1->GetParameter(); |
| 912 | if (TMath::Abs(p0[3]+p1[3])>fCutTheta) return kFALSE; |
| 913 | if (TMath::Abs(p0[1]-p1[1])>fCutMaxDz) return kFALSE; |
| 914 | if (TMath::Abs(p0[0]+p1[0])>fCutMaxD) return kFALSE; |
| 915 | |
| 916 | Double_t d0[3], d1[3]; |
| 917 | tr0->GetDirection(d0); |
| 918 | tr1->GetDirection(d1); |
| 919 | if (d0[0]*d1[0] + d0[1]*d1[1] + d0[2]*d1[2] >fCutMinDir) return kFALSE; |
| 920 | // |
| 921 | return kTRUE; |
| 922 | } |
| 923 | |
| 924 | |
| 925 | AliTrackReference * AliMaterialBudget::GetFirstTPCTrackRef(AliMCParticle *mcParticle) |
| 926 | { |
| 927 | // return first TPC track reference |
| 928 | if(!mcParticle) return 0; |
| 929 | |
| 930 | // find first track reference |
| 931 | // check direction to select proper reference point for looping tracks |
| 932 | Int_t nTrackRef = mcParticle->GetNumberOfTrackReferences(); |
| 933 | AliTrackReference *ref = 0; |
| 934 | AliTrackReference *refIn = 0; |
| 935 | for (Int_t iref = 0; iref < nTrackRef; iref++) { |
| 936 | ref = mcParticle->GetTrackReference(iref); |
| 937 | if(ref && (ref->DetectorId()==AliTrackReference::kTPC)) |
| 938 | { |
| 939 | //Float_t dir = ref->X()*ref->Px()+ref->Y()*ref->Py(); |
| 940 | //if(dir < 0.) break; |
| 941 | |
| 942 | refIn = ref; |
| 943 | break; |
| 944 | } |
| 945 | } |
| 946 | |
| 947 | return refIn; |
| 948 | } |
| 949 | |
| 950 | |
| 951 | AliTrackReference * AliMaterialBudget::GetAllTOFinfo(AliMCParticle *mcParticle, Int_t &nTrackRef, Int_t &nTrackRefITS, Int_t retValue) { |
| 952 | // |
| 953 | // |
| 954 | // |
| 955 | |
| 956 | if(!mcParticle) return 0; |
| 957 | Int_t counter = 0; |
| 958 | nTrackRef = 0; |
| 959 | nTrackRefITS = 0; |
| 960 | AliTrackReference *ref = 0; |
| 961 | for (Int_t iref = 0; iref < mcParticle->GetNumberOfTrackReferences(); iref++) { |
| 962 | ref = mcParticle->GetTrackReference(iref); |
| 963 | if(ref && (ref->DetectorId()==AliTrackReference::kTOF)) { |
| 964 | counter = iref; |
| 965 | nTrackRef++; |
| 966 | } |
| 967 | if(ref && (ref->DetectorId()==AliTrackReference::kITS)) nTrackRefITS++; |
| 968 | } |
| 969 | if (nTrackRef ==0) return 0; |
| 970 | if (retValue == 0) return mcParticle->GetTrackReference(counter - nTrackRef +1); |
| 971 | return mcParticle->GetTrackReference(counter); |
| 972 | |
| 973 | } |
| 974 | |
| 975 | |
| 976 | void AliMaterialBudget::FinishTaskOutput() |
| 977 | { |
| 978 | // |
| 979 | // According description in AliAnalisysTask this method is call |
| 980 | // on the slaves before sending data |
| 981 | // |
| 982 | Terminate("slave"); |
| 983 | gSystem->Exec("pwd"); |
| 984 | RegisterDebugOutput(); |
| 985 | |
| 986 | } |
| 987 | |
| 988 | |
| 989 | void AliMaterialBudget::RegisterDebugOutput(){ |
| 990 | // |
| 991 | // |
| 992 | // |
| 993 | // |
| 994 | // store - copy debug output to the destination position |
| 995 | // currently ONLY for local copy |
| 996 | TString dsName; |
| 997 | dsName=GetName(); |
| 998 | dsName+="Debug.root"; |
| 999 | dsName.ReplaceAll(" ",""); |
| 1000 | TString dsName2=fDebugOutputPath.Data(); |
| 1001 | gSystem->MakeDirectory(dsName2.Data()); |
| 1002 | dsName2+=gSystem->HostName(); |
| 1003 | gSystem->MakeDirectory(dsName2.Data()); |
| 1004 | dsName2+="/"; |
| 1005 | dsName2+=gSystem->BaseName(gSystem->pwd()); |
| 1006 | dsName2+="/"; |
| 1007 | gSystem->MakeDirectory(dsName2.Data()); |
| 1008 | dsName2+=dsName; |
| 1009 | AliInfo(Form("copy %s\t%s\n",dsName.Data(),dsName2.Data())); |
| 1010 | printf("copy %s\t%s\n",dsName.Data(),dsName2.Data()); |
| 1011 | TFile::Cp(dsName.Data(),dsName2.Data()); |
| 1012 | } |
| 1013 | |
| 1014 | |
| 1015 | Bool_t AliMaterialBudget::PropagateCosmicToDCA(AliExternalTrackParam *param0, AliExternalTrackParam *param1, Double_t mass){ |
| 1016 | // |
| 1017 | // param0 - upper part of cosmic track |
| 1018 | // param1 - lower part of cosmic track |
| 1019 | // |
| 1020 | // 0. Propagate both tracks to DCA to (0,0,0) |
| 1021 | // 1. After propagation to DCA rotate track param1 to corrdinate system of track1 <-> rotate param0 to coordinate system of param 1 ???? |
| 1022 | // 2. Propagate track 1 to refernce x from track0 |
| 1023 | // |
| 1024 | |
| 1025 | // step 0. |
| 1026 | |
| 1027 | Float_t d0 = param0->GetLinearD(0,0); |
| 1028 | Float_t d1 = param1->GetLinearD(0,0); |
| 1029 | Float_t dmax = TMath::Max(TMath::Abs(d0),TMath::Abs(d1)); |
| 1030 | // |
| 1031 | // propagate in the beginning taking all material into account |
| 1032 | // |
| 1033 | AliTracker::PropagateTrackToBxByBz(param0,dmax+1.,mass,0.5,kTRUE,0.99,-1.); |
| 1034 | AliTracker::PropagateTrackToBxByBz(param1,dmax+1.,mass,0.5,kTRUE,0.99,1.); |
| 1035 | // |
| 1036 | Double_t vtxx[3]={0,0,0}; |
| 1037 | Double_t svtxx[3]={0.000001,0.000001,100.}; |
| 1038 | AliESDVertex vtx(vtxx,svtxx); |
| 1039 | // |
| 1040 | Bool_t b0 = param0->PropagateToDCA(&vtx,AliTracker::GetBz(),1000); |
| 1041 | Bool_t b1 = param1->PropagateToDCA(&vtx,AliTracker::GetBz(),1000); |
| 1042 | |
| 1043 | if (!(b0 && b1)) return 0; |
| 1044 | |
| 1045 | // step 1. |
| 1046 | |
| 1047 | Float_t dAlpha = param0->GetAlpha(); |
| 1048 | param1->Rotate(dAlpha); |
| 1049 | |
| 1050 | // step 2. |
| 1051 | |
| 1052 | Float_t refX = param0->GetX(); |
| 1053 | param1->PropagateTo(refX,AliTracker::GetBz()); |
| 1054 | |
| 1055 | return kTRUE; |
| 1056 | |
| 1057 | |
| 1058 | } |
| 1059 | |
| 1060 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame