PndTrkCTFindTrackInXY2.h

Go to the documentation of this file.

#ifndef PndTrkCTFindTrackInXY2_H
#define PndTrkCTFindTrackInXY2_H 1

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


 struct FindTrackInXYProjection2_InputOutputData {
        // inputs;
        Double_t apotemastrawdetectormin;
        Double_t *Cosine;       // for the Legiandre fit;
        Short_t legiandre_nthetadiv;
        Short_t legiandre_nradiusdiv;
        bool* InclusionListStt;
        Double_t (*info)[7];
        Short_t maxhitsinfit;
        Short_t maxmvdpixelhitsintrack;
        Short_t maxmvdstriphitsintrack;
        bool * Mvdhits;
        Short_t nMvdPixelHit;
        Short_t nMvdStripHit;
        Short_t number_straws; // this is the number of all straws in the Stt detector;
        Short_t * nParContiguous;
        Short_t (*ListParContiguous)[6]; // this is the list of contiguous TubeID ;
        Double_t r_stt_inner_par_max;   // r_stt_inner_par_max ==> radius of the circumscribed
                                        // circumference to the
                                        //; outer hexagon defining THE INNER axial Stt straw region;
        Double_t *Sinus;        // for the Legiandre fit;
        Short_t *StrawCode;
        Short_t *StrawCode2;
        Short_t *TubeID ;  // TubeID[i] =  Tube ID corresponding to i-th Stt hit;
        Short_t *SttStrawOn;//  SttStrawOn[i] >= 0 --> it is the Stt hit number corresponding to Stt
                                // i-th Tube ID; tSttStrawOn[i] == -1 --> i-th Stt straw NOT hit;
        Double_t thetamax;      // input, maximum of Theta range (usually 2PI radians); 
        Double_t thetamin;      // input, minimum of Theta range (usually 0 radians);
        Double_t * XMvdPixel;
        Double_t * XMvdStrip;
        Double_t * YMvdPixel;
        Double_t * YMvdStrip;

        // outputs;
        Double_t* ALFA;
        Double_t* BETA;
        Short_t* Charge;
        Double_t *Fi_initial_helix_referenceframe;
        Double_t *Fi_low_limit;
        Double_t *Fi_up_limit;
        Double_t* GAMMA;
        Short_t *ListHitsinTrack;
        Short_t *ListMvdPixelHitsinTrack;
        Short_t *ListMvdStripHitsinTrack;
        Short_t *nMvdPixelHitsinTrack;
        Short_t *nMvdStripHitsinTrack;
        Short_t maxstthitsintrack;
        Short_t minimumhitspertrack;
        Short_t *nHitsinTrack;
        Double_t* Oxx;
        Double_t* Oyy;
        Double_t* Rr;
        Double_t* xTube;
        Double_t* yTube;
        Double_t* zTube;
        Double_t* xxyyTube;


//---


        Double_t apotemamaxskewstraw;
        Double_t deltanr;
        Double_t dimensionscitil;
        Short_t *FiConformalIndex;
        Double_t *Fi_final_helix_referenceframe;
        Short_t* HitsinBoxConf;
        int icounter; // counter useful for debugging;
        bool* InclusionListSciTil;
        Double_t (*infoparalConformal)[5];
        Short_t* ListSciTilHitsinTrack;
        Short_t* ListSttParHits;
        Short_t maxscitilhitsintrack;
        Short_t maxstthits;
        Short_t minouterhitspertrack;
        Short_t* nBoxConf;
        Short_t nFicell; // Fi cell of the seed hit;
        Short_t nfidivconformal;
        Short_t nRcell;  // R cell of the seed hit;
        Short_t nrdivconformal;
        Short_t nSciTilHits;
        Short_t* nSciTilHitsinTrack;
        Int_t nsttparhit;
        void * posizSciT;
        Double_t posizSciTilx;
        Double_t posizSciTily;
        Double_t* radiaConf;
        Short_t *RConformalIndex;
        Double_t rstrawdetectormax;
        Double_t strawradius;
        Double_t* S_SciTilHitsinTrack;
        Double_t *trajectory_vertex;
        bool* TypeConf;
        Double_t *U;
        Double_t *V;
        bool YesSciTil;
        };


class PndTrkCTFindTrackInXY2 : public TObject
{


 public:

  PndTrkCTFindTrackInXY2(){};


  ~PndTrkCTFindTrackInXY2(){};




//----------------

void AddMvdHitsToSttTracks(
        Double_t delta,                 // input;
        Double_t highqualitycut,        // input;
        Double_t FiRangeMvdLow,         // input;
        Double_t FiRangeMvdUp,          // input;
        Short_t maxmvdpixelhitsintrack, // input;
        Short_t maxmvdstriphitsintrack, // input;
        Short_t nMvdPixelHit,           // input;
        Short_t nMvdStripHit,           // input;
        Double_t  Ox,                   // input;
        Double_t  Oy,                   // input;
        Double_t  R,                    // input;
        Double_t * XMvdPixel,           // input;
        Double_t * XMvdStrip,           // input;
        Double_t * YMvdPixel,           // input;
        Double_t * YMvdStrip,           // input;

        Short_t &nPixelHitsinTrack, // output
        Short_t *ListPixelHitsinTrack, // output; dimensionality : [MAXMVDSTRIPHITSINTRACK];
        Short_t &nStripHitsinTrack, // output
        Short_t *ListStripHitsinTrack // output; dimensionality : [MAXMVDSTRIPHITSINTRACK];
        );


  Short_t AssociateSciTilHit(
        Double_t dimensionscitil,
        Double_t *esse, // output, list of  S of the SciTil hits associated. 
        bool* InclusionListSciTil,
        Short_t *List, // output, list of SciTil hits associated (max. 2);
        Short_t maxscitilhitsintrack,
        Short_t nSciTilHits,
        Double_t Oxx,
        Double_t Oyy,
        Double_t posizSciTil[][3],
        Double_t Rr
        );

 void DecideWhichAngularRangeAndCharge(
        Double_t fiCenter,              // input, fi of the (0,0) in the Helix reference frame;
        Double_t fi_low_limit[2],       // input; fi low limit of the Stt detector in the Helix frame;
        Double_t fi_up_limit[2],        // input; fi up limit of the Stt detector in the Helix frame;
        Double_t (*info)[7],            // input
        Short_t *ListHitsinTrack,       // input
        Short_t  nHitsinTrack,          // input
        Double_t Oxx,                   // input
        Double_t Oyy,                   // input

        Short_t & charge,               // output; charge of the particle;
        Double_t & FiRangeMvdLow,       // output; Fi range of possible Mvd hits; FiRangeMvdLow always
                                        // between 0 and 2PI; and always  FiRangeMvdLow<FiRangeMvdUp;
        Double_t & FiRangeMvdUp,        // output; Fi range of possible Mvd hits;
        Double_t & Fi_low_limit,        // output;
        Double_t & Fi_up_limit          // output;

        );


  bool FindTrackInXYProjection(
        struct FindTrackInXYProjection2_InputOutputData* InOut,
        int istampa,
        int IVOLTE
        );


  void  OrderingUsingConformal(
        Short_t  Charge,  // input;
        Double_t info[][7],     // input;
        Int_t nHits,  // input;
        Double_t oX,  // input;
        Double_t oY,  // input;
        Short_t * ListHits      // input and output (ordered);
        );


  void  OrderingUsingFi(
        Short_t  Charge,  // input;
        Double_t info[][7],     // input;
        Int_t nHits,  // input;
        Double_t oX,  // input;
        Double_t oY,  // input;
        Short_t * ListHits      // input and output (ordered);
        );

  void  OrderingUsingR(
        Double_t info[][7],     // input;
        Int_t nHits,  // input;
        Short_t * ListHits      // input and output (ordered);
        );




  Short_t TrkAssociatedParallelHitsToHelix5(
        Short_t *auxListHitsinTrack,
        bool *InclusionListStt,
        Double_t Fi_low,
        Double_t Fi_up,
        Double_t info[][7],
        Short_t *ListSttParHits,
        Int_t NhitsParallel,
        Double_t Oxx,
        Double_t Oyy,
        Double_t Rr,
        Double_t strawradius
        );



  Short_t TrkAssociatedParallelHitsToHelix6(
        Short_t *auxListHitsinTrack,
        bool *InclusionListStt,
        Double_t Fi_low,
        Double_t Fi_up,
        Double_t info[][7],
        Short_t *ListSttParHits,
        Int_t NhitsParallel,
        Double_t Oxx,
        Double_t Oyy,
        Double_t Rr,
        Double_t maximum_distance
        );
  ClassDef(PndTrkCTFindTrackInXY2,1);

};

#endif