PndTrackingQualityBarrelTaskNewLinks.cxx File Reference

#include <iostream>
#include <functional>
#include "TROOT.h"
#include "TClonesArray.h"
#include "THStack.h"
#include "FairRootManager.h"
#include "FairRun.h"
#include "FairRuntimeDb.h"
#include "FairHit.h"
#include "FairMultiLinkedData.h"
#include "PndMCTrack.h"
#include "PndTrack.h"
#include "PndMCEntry.h"
#include "PndSttHit.h"
#include "PndSttTube.h"
#include "PndSttMapCreator.h"
#include "RhoHistogram/RhoTuple.h"
#include "PndTrackingQualityBarrelTaskNewLinks.h"
#include "PndTrackingQualityBarrelAnalysisNewLinks.h"
#include "PndSdsHit.h"
#include "TGeoManager.h"
#include "PndGeoHandling.h"

Go to the source code of this file.

Functions
Bool_t	IsBarrelMVD (const PndSdsHit *hit)
	ClassImp (PndTrackingQualityBarrelTaskNewLinks)

Function Documentation

ClassImp

(

PndTrackingQualityBarrelTaskNewLinks

)

Bool_t IsBarrelMVD

(

const PndSdsHit *

hit

)

Definition at line 424 of file PndTrackingQualityBarrelTaskNewLinks.cxx.

References Double_t, exit(), fabs(), PndSdsHit::GetSensorID(), gGeoManager, and PndGeoHandling::Instance().

{
  if(!hit) return 0;  
  Int_t sensorID = hit->GetSensorID();
  gGeoManager->cd(PndGeoHandling::Instance()->GetPath(sensorID));
  TGeoHMatrix* transMat = gGeoManager->GetCurrentMatrix();
  if( !transMat ){
    std::cout<<"\n\nIsBarrelMVD():: NO transition Matrix for MVD sensor "<<sensorID<<"\n\n!!!"<<std::endl;
    return 0;
  }
  Double_t *mmm = transMat->GetRotationMatrix();
  if( !mmm ){
    std::cout<<"\n\nIsBarrelMVD():: Can not extract transition Matrix for MVD sensor "<<sensorID<<", something is completely wrong\n\n!!!"<<std::endl;
    exit(0);//SG!!
    return 0;
  }
  bool forward = ( fabs(mmm[6]) < 0.999 && fabs(mmm[7]) < 0.999 );
  return ( !forward );
}