TrackProducerFP420.h

Go to the documentation of this file.

#ifndef TrackProducerFP420_h
#define TrackProducerFP420_h

#include <string>
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "DataFormats/FP420Cluster/interface/TrackFP420.h"
#include "DataFormats/FP420Cluster/interface/ClusterFP420.h"
#include "DataFormats/FP420Cluster/interface/ClusterCollectionFP420.h"
#include "SimG4CMS/FP420/interface/FP420NumberingScheme.h"

#include <vector>
#include <algorithm>
#include <cmath>


class TrackProducerFP420 {
public:

  typedef std::vector<ClusterFP420>::const_iterator           ClusterFP420Iter;

  TrackProducerFP420(int, int, int, int, double, double, double, double, double, double, double, double, double, double, double, double, double, bool, bool, bool, bool, double, double, float, float, double, int, double, double);

  //  std::vector<TrackFP420> trackFinderMaxAmplitude(ClusterCollectionFP420 input);
  //  std::vector<TrackFP420> trackFinderMaxAmplitude2(ClusterCollectionFP420 input);

  //  // std::vector<TrackFP420> trackFinderVar1(ClusterCollectionFP420 input);
  //  //std::vector<TrackFP420> trackFinderVar2(ClusterCollectionFP420 input);

    //    std::vector<TrackFP420> trackFinderSophisticated(ClusterCollectionFP420 input);
    std::vector<TrackFP420> trackFinderSophisticated(edm::Handle<ClusterCollectionFP420> input, int det);

  //  std::vector<TrackFP420> trackFinder3D(ClusterCollectionFP420 input);

private:
  ClusterCollectionFP420 soutput;

  std::vector<TrackFP420> rhits; 

 // Number of Stations:
 int sn0;
 // Number of superplanes:
 int pn0;
 // Number of sensors
 int rn0;
 // Number of planes:
 int xytype;

 // shift of planes:
    bool UseHalfPitchShiftInX;
    bool UseHalfPitchShiftInY;
    bool UseHalfPitchShiftInXW;
    bool UseHalfPitchShiftInYW;

    //double zUnit; 
    double z420; 
    double zD2; 
    double zD3; 
    double pitchX;
    double pitchY;
    double pitchXW;
    double pitchYW;
        double ZGapLDet;
    //double ZBoundDet;
    double ZSiStep;
    double ZSiPlane;
    double ZSiDet ;
    double zBlade;
    double gapBlade;

    double dXX;
    double dYY;
    float chiCutX;
    float chiCutY;

    double zinibeg;

    int verbos;

    double XsensorSize;
    double YsensorSize;

};



#endif