PndTrkComparisonMCtruth.h

Go to the documentation of this file.

#ifndef PndTrkComparisonMCtruth_H
#define PndTrkComparisonMCtruth_H 1

// #include "FairRootManager.h"
#include "PndTrkVectors.h"
#include "PndMCTrack.h"
// Root includes
#include "TROOT.h"
#include "TClonesArray.h"


struct PndTrkComparisonMCtruth_io_Data{
        Double_t Bfield;
        Short_t *Charge;
        Double_t Cvel;
        Short_t *daTrackFoundaTrackMC;
        Double_t DIMENSIONSciTil ;
        Double_t Errorsqpixel;
        Double_t Errorsqstrip;
        Double_t *FI0;
        TClonesArray *fMCTrackArray;
        TClonesArray *fMvdMCPointArray;
        Short_t fSciTilMaxNumber;
        TClonesArray *fSciTHitArray;
        TClonesArray *fSciTPointArray;
        TClonesArray *fSttPointArray;
        FILE     *HANDLE;
        FILE     *HANDLE2;
        Double_t *info;
        int       istampa;
        int       IVOLTE;
        Double_t *KAPPA;
        bool     *keepit;
        bool     *InclusionListStt;
        Short_t  *ListMvdPixelHitsinTrack;
        Short_t  *ListMvdStripHitsinTrack;
        Short_t  *ListSciTilHitsinTrack;
        Short_t  *ListSttParHitsinTrack;
        Short_t  *ListSttSkewHitsinTrack;
        Short_t  *ListTrackCandHit;
        Short_t  *ListTrackCandHitType;
        int       MAXMCTRACKS;
        int       MAXMVDPIXELHITS;
        int       MAXMVDPIXELHITSINTRACK;
        int       Maxmvdmcpoints;
        int       MAXMVDSTRIPHITS;
        int       MAXMVDSTRIPHITSINTRACK;
        int       MAXSCITILHITS;
        int       MAXSCITILHITSINTRACK;
        int       MAXSTTHITS;
        int       maxstthitsintrack;
        int       MAXTRACKSPEREVENT;
        Short_t  *MCMvdPixelAloneList;
        Short_t  *MCMvdStripAloneList;
        Short_t  *MCParalAloneList;
        Short_t  *MCSkewAloneList;
        Double_t *MCSkewAloneX; // dimension : [MAXSTTHITS]
        Double_t *MCSkewAloneY; // dimension : [MAXSTTHITS]
        Short_t  *MvdPixelCommonList;
        Short_t  *MvdPixelSpuriList;
        Short_t  *MCSciTilAloneList; // equivalent to a matrix
                        //  [nTotalCandidates][nSciTilHits];
        Short_t  *MvdStripCommonList;
        Short_t  *MvdStripSpuriList;
        Short_t  *nHitsInMCTrack;
        Short_t  *nHitsInSciTile;
        Short_t  *nMCMvdPixelAlone;
        Short_t  *nMCMvdStripAlone;
        Short_t  *nMCParalAlone;
        Short_t  *nMCSciTilAlone;
        Short_t  *nMCSkewAlone;
        Short_t  *nMvdPixelCommon;
        Short_t  *nMvdPixelHitsinTrack;
        Short_t  *nMvdStripHitsinTrack;
        Short_t   nMvdPixelHit;
        Short_t  *nMvdPixelSpuriinTrack;
        Short_t  *nMvdStripCommon;
        Short_t   nMvdStripHit;
        Short_t  *nMvdStripSpuriinTrack;
        Short_t  *nParalCommon;
        Short_t  *nSciTilCommon;
        Short_t   nSciTilHits;
        Short_t  *nSciTilHitsinTrack;
        Short_t  *nSciTilSpuriinTrack;
        Short_t  *nSkewCommon;
        Short_t  *nSkewHitsInMCTrack;
        Short_t  *nSpuriParinTrack;
        Short_t  *nSpuriSkewinTrack;
        Int_t     nSttHit;
        Short_t  *nSttParHitsinTrack;
        Short_t  *nSttSkewHitsinTrack;
        Short_t   nTotalCandidates;
        Short_t  *OriginalSciTilList;
        Double_t *Ox;
        Double_t *Oy;
        Short_t  *ParalCommonList;
        Short_t  *ParSpuriList;
        Double_t *R;
        Double_t *refindexMvdPixel;
        Double_t *refindexMvdStrip;
        Short_t  *resultFitSZagain;
        Short_t  *SciTilCommonList; // equivalent to a matrix
                        // [nTotalCandidates][MAXSCITILHITSINTRACK];
        Short_t  *SciTilSpuriList; // equivalent to a matrix
                        //  [nTotalCandidates][MAXSCITILHITSINTRACK];
        Short_t  *SkewCommonList;
        Short_t  *SkewSpuriList;
        bool     *SttSZfit;
        Double_t *XMvdPixel;
        Double_t *XMvdStrip;
        Double_t *XSciTilCenter;
        Double_t *YMvdPixel;
        Double_t *YMvdStrip;
        Double_t *YSciTilCenter;
        Double_t *ZMvdPixel;
        Double_t *ZMvdStrip;
        Double_t *ZSciTilCenter;
        };

 // inizio cambio_in_perl

class PndTrkComparisonMCtruth
{


 public:

 PndTrkComparisonMCtruth(){};
 virtual ~PndTrkComparisonMCtruth(){};

 void AssociateFoundTrackstoMCquater(
        Double_t BFIELD,
        Double_t CVEL,
        Vec <Short_t> *daTrackFoundaTrackMC,
        TClonesArray *fMCTrackArray,
        Vec <Int_t>   *FromPixeltoMCTrack,
        Vec <Int_t>   *FromStriptoMCTrack,
        Vec <int> *FromSciTiltoMCTrackList,
                //  dimension   nSciTilHits*nMCTracks, equivalent to a matrix
                //  of dimension  [nSciTilHits][nMCTracks].
        Vec <bool> *keepit,
        Vec <Double_t> *info,
        Vec <Short_t> *ListSttParHitsinTrack,
        Vec <Short_t> *ListMvdPixelHitsinTrack,
        Vec <Short_t> *ListSciTilHitsinTrack,
        Vec <Short_t> *ListSttSkewHitsinTrack,
        Vec <Short_t> *ListMvdStripHitsinTrack,
        int MAXMVDPIXELHITSINTRACK,
        int MAXMVDSTRIPHITSINTRACK,
        int MAXSCITILHITSINTRACK,
        int MAXSTTHITSINTRACK,
        Vec <int> *nFromSciTiltoMCTrack,
        Vec <Short_t> *nSttParHitsinTrack,
        int nMCTracks,
        Vec <Short_t> *nMvdPixelHitsinTrack,
        Short_t nSciTilHits,
        Vec <Short_t> *nSciTilHitsinTrack,
                // dimension [MAXTRACKSPEREVENT][MAXSCITILHITSINTRACK]
        Vec <Short_t> *nSttSkewHitsinTrack,
        Vec <Short_t> *nMvdStripHitsinTrack,
        Short_t  nTracksFoundSoFar,
        Vec <Double_t> *Ox,
        Vec <Double_t> *Oy,
        Vec <Double_t> *R,
        Vec <Double_t> *X1,
        Vec <Double_t> *Y1,
        Vec <Double_t> *X2,
        Vec <Double_t> *Y2,
        Vec <Double_t> *X3,
        Vec <Double_t> *Y3,
        int IVOLTE

                                        );

 int ComparisonwithMC(
        PndTrkComparisonMCtruth_io_Data ioData
                );

 Double_t FindDistance(
        Double_t Oxx,   //  center from wich distance is calculated
        Double_t Oyy,   //  center from wich distance is calculated
        Double_t Rr,
        Double_t tanlow,
        Double_t tanmid,
        Double_t tanup,
        Double_t alfa,  //  intersection circumference parameter
        Double_t beta,  //  intersection circumference parameter
        Double_t gamma  //  intersection circumference parameter
                        );
void getMCInfo(
        Double_t BFIELD,
        Double_t CVEL,
        Double_t * Cx,
        Double_t * Cy,
        TClonesArray *fMCTrackArray,
        Int_t MCTrack,
        Double_t * Rr
        );

void MvdMatchedSpurioustoTrackCand(
        Vec <Short_t> *daTrackFoundaTrackMC,
        Vec <Int_t> *FromPixeltoMCTrack,
        Vec <Int_t> *FromStriptoMCTrack,
        Vec <bool> *keepit,
        Vec <Short_t> *ListMvdPixelHitsinTrack,
        Vec <Short_t> *ListMvdStripHitsinTrack,
        int MAXMVDPIXELHITSINTRACK,
        int MAXMVDSTRIPHITSINTRACK,
        Short_t nMvdPixelHit,
        Short_t nMvdStripHit,
        Vec <Short_t> *nMvdPixelHitsinTrack,
        Vec <Short_t> *nMvdStripHitsinTrack,
        Short_t nSttTrackCand,          // input

        Vec <Short_t> *nMvdPixelCommon,
        Vec <Short_t> *MvdPixelCommonList,
        Vec <Short_t> *nMvdPixelSpuriinTrack,
        Vec <Short_t> *MvdPixelSpuriList,
        Vec <Short_t> *nMCMvdPixelAlone,
        Vec <Short_t> *MCMvdPixelAloneList,

        Vec <Short_t> *nMvdStripCommon,
        Vec <Short_t> *MvdStripCommonList,
        Vec <Short_t> *nMvdStripSpuriinTrack,
        Vec <Short_t> *MvdStripSpuriList,
        Vec <Short_t> *nMCMvdStripAlone,
        Vec <Short_t> *MCMvdStripAloneList      
                                );

void MvdMatchtoMC(
        Double_t ERRORSQPIXEL,
        Double_t ERRORSQSTRIP,
        TClonesArray *fMvdMCPointArray,
        Short_t nMvdMCPoint,
        int istampa,
        int IVOLTE,
        Short_t nMvdPixelHit,
        Short_t nMvdStripHit,
        Vec <Double_t> *refindexMvdPixel,
        Vec <Double_t> *refindexMvdStrip,
        Vec <Double_t> *XMvdPixel,
        Vec <Double_t> *XMvdStrip,
        Vec <Double_t> *YMvdPixel,
        Vec <Double_t> *YMvdStrip,
        Vec <Double_t> *ZMvdPixel,
        Vec <Double_t> *ZMvdStrip,
        Vec <Int_t> *FromPixeltoMCTrack,
        Vec <Int_t> *FromStriptoMCTrack 
        );

 void SciTilMatchtoMC(
        Double_t BFIELD,
        Double_t CVEL,
        Double_t DIMENSIONSCITIL,
        TClonesArray *fMCTrackArray,
        Vec <int> *FromSciTiltoMCTrackList,
        TClonesArray *fSciTHitArray,
        Short_t fSciTilMaxNumber,
        TClonesArray *fSciTPointArray,
        Vec <int> *nFromSciTiltoMCTrack,
        Short_t *nHitsInSciTile,
        int nMCTracks,
        Short_t nSciTilHits,
        Short_t *OriginalSciTilList,
        Vec <Double_t> *XSciTilCenter,
        Vec <Double_t> *YSciTilCenter,
        Vec <Double_t> *ZSciTilCenter

                );

 void SciTilMatchedSpurioustoTrackCand(
        Vec <Short_t> *daTrackFoundaTrackMC,
        Vec <int> *FromSciTiltoMCTrackList,
                                        //  of dimension  [nSciTilHits][nMCTracks]
        Vec <bool> *keepit,
        Vec <Short_t> *ListSciTilHitsinTrack,
                                        // [MAXTRACKSPEREVENT][MAXSCITILHITSINTRACK]
        int MAXSCITILHITSINTRACK,       // input
        Short_t *MCSciTilAloneList,     // output; equivalent to a matrix of dimension
                                        // [MAXTRACKSPEREVENT][MAXSCITILHITSINTRACK]
        Vec <int> *nFromSciTiltoMCTrack,
        Short_t *nMCSciTilAlone,        // output
        int nMCTracks,                  // input
        Short_t nSciTilHits,            // input
        Vec <Short_t>  *nSciTilHitsinTrack,
        Short_t *nSciTilCommon,         // output
        Short_t *nSciTilSpuriinTrack,   // output
        Short_t nSttTrackCand,          // input
        Short_t *SciTilCommonList,      // output; equivalent to a matrix of dimension
                                        // [MAXTRACKSPEREVENT][MAXSCITILHITSINTRACK]
        Short_t *SciTilSpuriList        // output; equivalent to a matrix of dimension
                                        // [MAXTRACKSPEREVENT][MAXSCITILHITSINTRACK].
                        );

 void stampaMCTracks(
        Double_t BFIELD,
        Double_t CVEL,
        TClonesArray* fMCTrackArray,
        int nMCTracks
                        );

void SttMatchedSpurious(
        Vec <Short_t> *daTrackFoundaTrackMC,
        Vec <bool> *InclusionListStt,
        Vec <Double_t> *info,
        Vec <bool> *keepit,
        int MAXSTTHITS,
        int MAXSTTHITSINTRACK,
        int MAXTRACKSPEREVENT,
        Vec <Short_t> *ListSttParHitsinTrack,
        Vec <Short_t> *ListSttSkewHitsinTrack,
        Vec <Short_t> *MCParalAloneList,
        Vec <Short_t> *MCSkewAloneList,
        Vec <Short_t> *nHitsInMCTrack,
        Vec <Short_t> *nSttParHitsinTrack,
        Vec <Short_t> *nMCParalAlone,
        Vec <Short_t> *nMCSkewAlone,
        Vec <Short_t> *nParalCommon,
        Vec <Short_t> *nSkewCommon,
        Vec <Short_t> *nSkewHitsInMCTrack,
        Vec <Short_t> *nSttSkewHitsinTrack,
        Vec <Short_t> *nSpuriParinTrack,
        Vec <Short_t> *nSpuriSkewinTrack,
        Short_t  nSttHits,
        Short_t  nTracksFoundSoFar, // those found by PR
        Vec <Short_t> *ParalCommonList,
        Vec <Short_t> *ParSpuriList,
        Vec <Short_t> *SkewCommonList,
        Vec <Short_t> *SkewSpuriList

        );

 // fine cambio_in_perl

  ClassDef(PndTrkComparisonMCtruth,1);

};

#endif