PndTrkPlotMacros2.h

Go to the documentation of this file.

#ifndef PndTrkPlotMacros2_H
#define PndTrkPlotMacros2_H 1

#include "PndTrkVectors.h"

// Root includes
#include "TClonesArray.h"
#include "TROOT.h"


struct PndTrkPlotMacros2_InputData{
        Double_t apotemastrawdetectormin;
        Double_t apotemamaxinnerparstraw;
        Double_t apotemamaxskewstraw;
        Double_t apotemaminouterparstraw;
        Double_t apotemaminskewstraw;
        Double_t bfield;
        Short_t *Charge;
        Double_t cvel;
        Short_t *daTrackFoundaTrackMC;
        Double_t dimensionscitil;
        bool doMcComparison;
        Double_t *FI0;
        TClonesArray *fMCTrackArray;
        TClonesArray *fSttPointArray;
        Double_t *info;
        int IVOLTE;
        Double_t *KAPPA;
        bool *keepit;
        bool *InclusionListSciTil;
        bool *InclusionListStt;
        int istampa;
        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 MAXMVDSTRIPHITS;
        int MAXMVDSTRIPHITSINTRACK;
        int MAXSCITILHITSINTRACK;
        int MAXSCITILHITS;
        int MAXSTTHITS;
        int maxstthitsintrack;
        int MAXTRACKSPEREVENT;
        Short_t *MCMvdPixelAloneList;
        Short_t *MCMvdStripAloneList;
        Short_t *MCParalAloneList;

        Short_t *MCSciTilAloneList;// equivalent to a matrix
                //  [nTotalCandidates][nSciTilHits]

        Short_t *MCSkewAloneList;
        Double_t *MCSkewAloneX;
        Double_t *MCSkewAloneY;
        Short_t *MvdPixelCommonList;
        Short_t *MvdPixelSpuriList;
        Short_t *MvdStripCommonList;
        Short_t *MvdStripSpuriList;
        Short_t NFIDIVCONFORMAL;
        Short_t *nMCMvdPixelAlone;
        Short_t *nMCMvdStripAlone;
        Short_t *nMCParalAlone;

        Short_t *nMCSciTilAlone;

        Short_t *nMCSkewAlone;
        Short_t nMCTracks;
        Short_t *nMvdPixelCommon;
        Short_t  nMvdPixelHit;
        Short_t *nMvdPixelHitsinTrack;
        Short_t *nMvdPixelSpuriinTrack;
        Short_t *nMvdStripCommon;
        Short_t nMvdStripHit;
        Short_t *nMvdStripHitsinTrack;
        Short_t *nMvdStripSpuriinTrack;
        Short_t *nParalCommon;
        Short_t NRDIVCONFORMAL;
        Short_t *nSciTilCommon;

        Short_t nSciTilHits;
        Short_t *nSciTilHitsinTrack;

        Short_t *nSciTilSpuriinTrack;

        Short_t *nSkewCommon;
        Short_t *nSpuriParinTrack;
        Int_t    nSttHit;
        Int_t    nSttParHit;
        Short_t *nSttParHitsinTrack;
        Int_t    nSttSkewHit;
        Short_t *nSttSkewHitsinTrack;
        Short_t  nTotalCandidates;
        Short_t *nTrackCandHit;
        Short_t number_straws;
        Double_t *Ox;
        Double_t *Oy;
        Short_t *ParalCommonList;
        Short_t *ParSpuriList;
        Double_t *posizSciTil;
        Double_t *R;
        Double_t *radiaConf;
        Double_t rstrawdetectormax;
        Double_t *SchosenSkew;


        Short_t *SciTilCommonList ;// equivalent to a matrix
                        // [nTotalCandidates][MAXSCITILHITSINTRACK]
        Short_t *SciTilSpuriList ;// equivalent to a matrix
                        // [nTotalCandidates][MAXSCITILHITSINTRACK]


        Double_t *sigmaXMvdPixel;
        Double_t *sigmaXMvdStrip;
        Double_t *sigmaYMvdPixel;
        Double_t *sigmaYMvdStrip;
        Short_t *SkewCommonList;
        TClonesArray *SttTubeArray;
        Short_t *StrawCode;
        Short_t *StrawCode2;
        Double_t verticalgap;
        Double_t *XMvdPixel;
        Double_t *XMvdStrip;
        Double_t *YMvdPixel;
        Double_t *YMvdStrip;
        Double_t *WDX;
        Double_t *WDY;
        Double_t *WDZ;
        Double_t *ZMvdPixel;
        Double_t *ZMvdStrip;
        } ;

//inizio cambio_in_perl

class PndTrkPlotMacros2 : public TObject
{


 public:


 PndTrkPlotMacros2(){};
 ~PndTrkPlotMacros2(){};

void disegnaAssiXY(
        FILE * MACRO,
        double xmin,
        double xmax,
        double ymin,
        double ymax
        );

void disegnaSciTilHit(
        int colorcode,  // goes in the SetColor function of root;
        Double_t DIMENSIONSCITIL,
        FILE * MACRO,
        double posx,
        double posy,
        int ScitilHit,
        int tipo// if 0 then SciTil draw in XY; if 1 then SciTil draw in SZ; else
                // SciTil draw in UV.
        );

void DrawBiHexagonInMacro(
        Double_t vgap,
        FILE * MACRO,
        Double_t Ami,
        Double_t Ama,
        Short_t color,
        char *name
        );


void DrawHexagonCircleInMacro(
        Double_t GAP,
        FILE * MACRO,
        Double_t ApotemaMin,
        Double_t Rma,
        Short_t color,
        char *name
        );


 void SttInfoXYZParal(
        Vec <Double_t> *info,
        Short_t infopar,
        Double_t Oxx,
        Double_t Oyy,
        Double_t Rr,
        Double_t KAPPA,
        Double_t FI0,
        Short_t Charge,
        Double_t *Posiz
        );

 void  WriteAllMacros(
        PndTrkPlotMacros2_InputData In_Put
        );



 void WriteMacroAllHitsRestanti(
        Double_t APOTEMAMAXINNERPARSTRAW,
        Double_t APOTEMAMAXSKEWSTRAW,
        Double_t APOTEMAMINOUTERPARSTRAW,
        Double_t APOTEMAMINSKEWSTRAW,
        Vec <bool> *InclusionListSciTil,
        Vec <bool> *InclusionListStt,
        Vec <Double_t> *info,
        int IVOLTE,
        Vec <bool> *keepit,
        Vec <Short_t> *ListTrackCandHit,
        Vec <Short_t> *ListTrackCandHitType,
        int MAXMVDPIXELHITS,
        int MAXMVDPIXELHITSINTRACK,
        int MAXMVDSTRIPHITS,
        int MAXMVDSTRIPHITSINTRACK,
        int MAXSCITILHITSINTRACK,
        int MAXSTTHITS,
        int MAXSTTHITSINTRACK,
        int MAXTRACKSPEREVENT,
        Short_t nMvdPixelHit,
        Short_t nMvdStripHit,
        Short_t nSciTilHit,
        Short_t nSttHit,
        Short_t nSttParHit,
        Short_t nSttSkewHit,
        Short_t nSttTrackCand,
        Vec <Short_t> *nTrackCandHit,
        Vec <Double_t> *posizSciTil,
        Double_t RSTRAWDETECTORMAX,
        Double_t APOTEMASTRAWDETECTORMIN,
        Double_t VERTICALGAP,
        Vec <Double_t> *XMvdPixel,
        Vec <Double_t> *XMvdStrip,
        Vec <Double_t> *YMvdPixel,
        Vec <Double_t> *YMvdStrip

        );

 void WriteMacroParallelHitsGeneral(
        Double_t APOTEMAMAXINNERPARSTRAW,
        Double_t APOTEMAMAXSKEWSTRAW,
        Double_t APOTEMAMINOUTERPARSTRAW,
        Double_t APOTEMAMINSKEWSTRAW,
        Double_t BFIELD,
        Double_t CVEL,
        Double_t DIMENSIONSCITIL,
        bool doMcComparison,
        TClonesArray *fMCTrackArray,
        Int_t Nhits,
        Vec <Double_t> *info,
        PndTrkPlotMacros2_InputData In_Put,
        int IVOLTE,
        Short_t nMCTracks,
        Short_t nMvdPixelHit,
        Short_t nMvdStripHit,
        Short_t nSciTilHits,
        Short_t nTracksFoundSoFar,
        Vec <bool> *keepit,
        Vec <Double_t> *FI0,
        Vec <Double_t> *Ox,
        Vec <Double_t> *Oy,
        Vec <Double_t> *posizSciTil,
        Vec <Double_t> *primoangolo,
        Vec <Double_t> *R,
        Double_t RSTRAWDETECTORMAX,
        Double_t APOTEMASTRAWDETECTORMIN,
        Vec <Double_t> *sigmaXMvdPixel,
        Vec <Double_t> *sigmaXMvdStrip,
        Vec <Double_t> *sigmaYMvdPixel,
        Vec <Double_t> *sigmaYMvdStrip,
        Vec <Double_t> *ultimoangolo,
        Double_t VERTICALGAP,
        Vec <Double_t> *XMvdPixel,
        Vec <Double_t> *XMvdStrip,
        Vec <Double_t> *YMvdPixel,
        Vec <Double_t> *YMvdStrip
        );


 void WriteMacroParallelHitsGeneralConformalwithMC(
        Double_t APOTEMAMAXINNERPARSTRAW,
        Double_t APOTEMAMAXSKEWSTRAW,
        Double_t APOTEMAMINOUTERPARSTRAW,
        Double_t APOTEMAMINSKEWSTRAW,
        Double_t BFIELD,
        Double_t CVEL,
        Double_t DIMENSIONSCITIL,
        bool doMcComparison,
        TClonesArray *fMCTrackArray,
        Int_t Nhits,
        Vec <Double_t> *info,
        PndTrkPlotMacros2_InputData In_Put,
        int IVOLTE,
        Short_t nMCTracks,
        Short_t nMvdPixelHit,
        Short_t nMvdStripHit,
        Short_t nSciTilHits,
        Short_t nTracksFoundSoFar,
        Vec <bool> *keepit,
        Vec <Double_t> *FI0,
        Vec <Double_t> *Ox,
        Vec <Double_t> *Oy,
        Vec <Double_t> *posizSciTil,
        Vec <Double_t> *primoangolo,
        Vec <Double_t> *R,
        Double_t RSTRAWDETECTORMAX,
        Double_t APOTEMASTRAWDETECTORMIN,
        Vec <Double_t> *sigmaXMvdPixel,
        Vec <Double_t> *sigmaXMvdStrip,
        Vec <Double_t> *sigmaYMvdPixel,
        Vec <Double_t> *sigmaYMvdStrip,
        Vec <Double_t> *ultimoangolo,
        Double_t VERTICALGAP,
        Vec <Double_t> *XMvdPixel,
        Vec <Double_t> *XMvdStrip,
        Vec <Double_t> *YMvdPixel,
        Vec <Double_t> *YMvdStrip
        );


 void WriteMacroParallel_MvdHitsGeneralConformalwithMC(
        Double_t APOTEMAMAXINNERPARSTRAW,
        Double_t APOTEMAMAXSKEWSTRAW,
        Double_t APOTEMAMINOUTERPARSTRAW,
        Double_t APOTEMAMINSKEWSTRAW,
        Double_t BFIELD,
        Double_t CVEL,
        Double_t DIMENSIONSCITIL,
        bool doMcComparison,
        TClonesArray *fMCTrackArray,
        Int_t Nhits,
        Vec <Double_t> *info,
        PndTrkPlotMacros2_InputData In_Put,
        int IVOLTE,
        Short_t nMCTracks,
        Short_t nMvdPixelHit,
        Short_t nMvdStripHit,
        Short_t nSciTilHits,
        Short_t nTracksFoundSoFar,
        Vec <bool> *keepit,
        Vec <Double_t> *FI0,
        Vec <Double_t> *Ox,
        Vec <Double_t> *Oy,
        Vec <Double_t> *posizSciTil,
        Vec <Double_t> *primoangolo,
        Vec <Double_t> *R,
        Double_t RSTRAWDETECTORMAX,
        Double_t APOTEMASTRAWDETECTORMIN,
        Vec <Double_t> *sigmaXMvdPixel,
        Vec <Double_t> *sigmaXMvdStrip,
        Vec <Double_t> *sigmaYMvdPixel,
        Vec <Double_t> *sigmaYMvdStrip,
        Vec <Double_t> *ultimoangolo,
        Double_t VERTICALGAP,
        Vec <Double_t> *XMvdPixel,
        Vec <Double_t> *XMvdStrip,
        Vec <Double_t> *YMvdPixel,
        Vec <Double_t> *YMvdStrip
        );

 void WriteMacroSkewAssociatedHitswithMC(
        Double_t *ESSE,
        Double_t *ESSEalone,
        PndTrkPlotMacros2_InputData In_Put,
        int iNome, // questo e' per il nome delle Macro solamente.
        int iTrack
        );


 void WriteMacroSkewAssociatedHitswithMC_Degree(
        Double_t *ESSE,
        Double_t *ESSEalone,
        PndTrkPlotMacros2_InputData In_Put,
        int iNome, // questo e' per il nome delle Macro solamente.
        int iTrack
        );


 void WriteMacroSttParallelAssociatedHitsandMvdwithMC(
        PndTrkPlotMacros2_InputData In_Put,
        Double_t Oxx,
        Double_t Oyy,
        Double_t Rr,
        Double_t primoangolo,
        Double_t ultimoangolo,
        Short_t Nhits,
        int iTrack,
        int iNome,
        Short_t daSttTrackaMCTrack,
        Short_t nMvdPixelHitsAssociatedToSttTra,
        Short_t nMvdStripHitsAssociatedToSttTra,
        Short_t nSkewHitsinTrack
        );
// fine cambio_in_perl  ;

 void WriteMacroSttParallel(
        PndTrkPlotMacros2_InputData In_Put
        );



 void WriteMacroSttParallelExternal(
        PndTrkPlotMacros2_InputData In_Put
        );

  ClassDef(PndTrkPlotMacros2,1);

};

#endif