#include <FP420TrackMain.h>

Public Member Functions
	FP420TrackMain (const edm::ParameterSet &conf)

void	run (edm::Handle< ClusterCollectionFP420 > &input, std::auto_ptr< TrackCollectionFP420 > &toutput)
	Runs the algorithm. More...

	~FP420TrackMain ()

Private Attributes
double	chiCutX_

double	chiCutY_

edm::ParameterSet	conf_

int	dn0

double	dXX_

double	dYY_

TrackProducerFP420 *	finderParameters_

double	gapBlade_

double	pitchX_

double	pitchXW_

double	pitchY_

double	pitchYW_

int	pn0_

int	rn0_

int	sn0_

std::string	trackMode_

bool	UseHalfPitchShiftInX_

bool	UseHalfPitchShiftInXW_

bool	UseHalfPitchShiftInY_

bool	UseHalfPitchShiftInYW_

bool	validTrackerizer_

int	verbosity

double	XsensorSize_

int	xytype_

double	YsensorSize_

double	z420_

double	zBlade_

double	zD2_

double	zD3_

double	ZGapLDet_

double	zinibeg_

double	ZSiDet_

double	ZSiPlane_

double	ZSiStep_

Detailed Description

Definition at line 21 of file FP420TrackMain.h.

Constructor & Destructor Documentation

FP420TrackMain::FP420TrackMain ( const edm::ParameterSet & conf )

Definition at line 23 of file FP420TrackMain.cc.

References chiCutX_, chiCutY_, conf_, gather_cfg::cout, dn0, dXX_, dYY_, finderParameters_, gapBlade_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), pitchX_, pitchXW_, pitchY_, pitchYW_, pn0_, rn0_, sn0_, trackMode_, UseHalfPitchShiftInX_, UseHalfPitchShiftInXW_, UseHalfPitchShiftInY_, UseHalfPitchShiftInYW_, validTrackerizer_, verbosity, XsensorSize_, xytype_, YsensorSize_, z420_, zBlade_, zD2_, zD3_, ZGapLDet_, zinibeg_, ZSiDet_, ZSiPlane_, and ZSiStep_.

                                                          :conf_(conf)  { 
   
   verbosity   = conf_.getUntrackedParameter<int>("VerbosityLevel");
   trackMode_  =  conf_.getParameter<std::string>("TrackModeFP420");
   dn0   = conf_.getParameter<int>("NumberFP420Detectors");
   sn0_ = conf_.getParameter<int>("NumberFP420Stations");
   pn0_ = conf_.getParameter<int>("NumberFP420SPlanes");
   rn0_ = 7;
   xytype_ = conf_.getParameter<int>("NumberFP420SPTypes");
   z420_           = conf_.getParameter<double>("z420");
   zD2_            = conf_.getParameter<double>("zD2");
   zD3_            = conf_.getParameter<double>("zD3");
   dXX_ = conf_.getParameter<double>("dXXFP420");
   dYY_ = conf_.getParameter<double>("dYYFP420");
   chiCutX_ = conf_.getParameter<double>("chiCutX420");
   chiCutY_ = conf_.getParameter<double>("chiCutY420");
   
   if (verbosity > 0) {
     std::cout << "FP420TrackMain constructor::" << std::endl;
     std::cout << "sn0=" << sn0_ << " pn0=" << pn0_ << " xytype=" << xytype_ << std::endl;
     std::cout << "trackMode = " << trackMode_ << std::endl;
     std::cout << "dXX=" << dXX_ << " dYY=" << dYY_ << std::endl;
     std::cout << "chiCutX=" << chiCutX_ << " chiCutY=" << chiCutY_ << std::endl;
   }
     // zD2_ = 1000.;  // dist between centers of 1st and 2nd stations
     // zD3_ = 8000.;  // dist between centers of 1st and 3rd stations
     
   UseHalfPitchShiftInX_= true;
   UseHalfPitchShiftInXW_= true;
   UseHalfPitchShiftInY_= true;
   UseHalfPitchShiftInYW_= true;
 
   pitchX_= 0.050;
   pitchY_= 0.050;// 
   pitchXW_= 0.400;
   pitchYW_= 0.400;// 
 
   XsensorSize_=8.0;
   YsensorSize_=7.2;
 
 //
   zBlade_ = 5.00;
   gapBlade_ = 1.6;
   double gapSupplane = 1.6;
   ZSiPlane_=2*zBlade_+gapBlade_+gapSupplane;
   
   double ZKapton = 0.1;
   ZSiStep_=ZSiPlane_+ZKapton;
   
   double ZBoundDet = 0.020;
   double ZSiElectr = 0.250;
   double ZCeramDet = 0.500;
 
   double eee1=11.;
   double eee2=12.;
   zinibeg_ = (eee1-eee2)/2.;
 //
   ZSiDet_ = 0.250;
 //
   ZGapLDet_= zBlade_/2-(ZSiDet_+ZSiElectr+ZBoundDet+ZCeramDet/2);
 //
     if (verbosity > 1) {
       std::cout << "FP420TrackMain constructor::" << std::endl;
       std::cout << " zD2=" << zD2_ << " zD3=" << zD3_ << " zinibeg =" << zinibeg_ << std::endl;
       std::cout << " UseHalfPitchShiftInX=" << UseHalfPitchShiftInX_ << " UseHalfPitchShiftInY=" << UseHalfPitchShiftInY_ << std::endl;
       std::cout << " UseHalfPitchShiftInXW=" << UseHalfPitchShiftInXW_ << " UseHalfPitchShiftInYW=" << UseHalfPitchShiftInYW_ << std::endl;
       std::cout << " pitchX=" << pitchX_ << " pitchY=" << pitchY_ << std::endl;
       std::cout << " pitchXW=" << pitchXW_ << " pitchYW=" << pitchYW_ << std::endl;
       std::cout << " zBlade_=" << zBlade_ << " gapBlade_=" << gapBlade_ << std::endl;
       std::cout << " ZKapton=" << ZKapton << " ZBoundDet=" << ZBoundDet << std::endl;
       std::cout << " ZSiElectr=" << ZSiElectr << " ZCeramDet=" << ZCeramDet << std::endl;
       std::cout << " ZSiDet=" << ZSiDet_ << " gapSupplane=" << gapSupplane << std::endl;
     }
 
 
 
       if ( trackMode_ == "TrackProducerSophisticatedFP420" ) {
   
   
   //trackMode_ == "TrackProducerVar1FP420" ||
   //trackMode_ == "TrackProducerVar2FP420" ||
 
      // if ( trackMode_ == "TrackProducerMaxAmplitudeFP420" ||
 //     trackMode_ == "TrackProducerMaxAmplitude2FP420"  ||
 //     trackMode_ == "TrackProducerSophisticatedFP420"  ||
 //     trackMode_ == "TrackProducer3DFP420" )  {
 
     finderParameters_ = new TrackProducerFP420(sn0_, pn0_, rn0_, xytype_, z420_, zD2_, zD3_,
                            pitchX_, pitchY_,
                            pitchXW_, pitchYW_,
                            ZGapLDet_, ZSiStep_,
                            ZSiPlane_, ZSiDet_,zBlade_,gapBlade_,
                            UseHalfPitchShiftInX_, UseHalfPitchShiftInY_,
                            UseHalfPitchShiftInXW_, UseHalfPitchShiftInYW_,
                            dXX_,dYY_,chiCutX_,chiCutY_,zinibeg_,verbosity,
                            XsensorSize_,YsensorSize_);
     validTrackerizer_ = true;
       } 
       else {
     std::cout << "ERROR:FP420TrackMain: No valid finder selected" << std::endl;
     validTrackerizer_ = false;
       }
 }

FP420TrackMain::~FP420TrackMain ( )

Definition at line 129 of file FP420TrackMain.cc.

References finderParameters_.

                                 {
   if ( finderParameters_ != 0 ) {
     delete finderParameters_;
   }
 }

Member Function Documentation

void FP420TrackMain::run	(	edm::Handle< ClusterCollectionFP420 > &	input,
		std::auto_ptr< TrackCollectionFP420 > &	toutput
	)

Runs the algorithm.

Definition at line 137 of file FP420TrackMain.cc.

References TrackFP420::ax(), TrackFP420::ay(), TrackFP420::bx(), TrackFP420::by(), TrackFP420::chi2x(), TrackFP420::chi2y(), gather_cfg::cout, dn0, finderParameters_, first, estimatePileup::inputRange, TrackFP420::nclusterx(), TrackFP420::nclustery(), TrackProducerFP420::trackFinderSophisticated(), trackMode_, validTrackerizer_, and verbosity.

Referenced by cms::TrackerizerFP420::produce().

 {
   
   if ( validTrackerizer_ ) {
     
     int number_detunits          = 0;
     int number_localelectroderechits = 0;
     /*
       for (int sector=1; sector<sn0_; sector++) {
       for (int zmodule=1; zmodule<pn0_; zmodule++) {
       for (int zside=1; zside<rn0_; zside++) {
       int sScale = 2*(pn0-1);
       //      int index = FP420NumberingScheme::packFP420Index(det, zside, sector, zmodule);
       // intindex is a continues numbering of FP420
       int zScale=2;  unsigned int detID = sScale*(sector - 1)+zScale*(zmodule - 1)+zside;
       ClusterMap.clear();
       ClusterCollectionFP420::Range clusterRange;
       clusterRange = input.get(detID);
       ClusterCollectionFP420::ContainerIterator clusterRangeIteratorBegin = clusterRange.first;
       ClusterCollectionFP420::ContainerIterator clusterRangeIteratorEnd   = clusterRange.second;
       for ( ;sort_begin != sort_end; ++sort_begin ) {
       ClusterMap.push_back(*sort_begin);
       } // for
       
       }//for
       }//for
       }//for
     */
     // get vector of detunit ids
     //    const std::vector<unsigned int> detIDs = input->detIDs();
     
     // to be used in put (besause of 0 in track collection for: 1) 1st track and 2) case of no track)
     // ignore 0, but to save info for 1st track record it second time on place 1   .
     
     bool first = true;
     // loop over detunits
     for (int det=1; det<dn0; det++) {
       ++number_detunits;
       int StID = 1111;
       if(det==2) StID = 2222;
       std::vector<TrackFP420> collector;
       //        std::vector<TrackFP420> collector;
       collector.clear();
       
       // if ( trackMode_ == "TrackProducerMaxAmplitudeFP420") {
       //     collector = finderParameters_->trackFinderMaxAmplitude(input); //std::vector<TrackFP420> collector;
       // }// if ( trackMode
       // else if (trackMode_ == "TrackProducerMaxAmplitude2FP420" ) {
       //     collector = finderParameters_->trackFinderMaxAmplitude2(input); //
       //  }// if ( trackMode
       /*
     else if (trackMode_ == "TrackProducerVar1FP420" ) {
     collector = finderParameters_->trackFinderVar1(input); //
     }// if ( trackMode
     else if (trackMode_ == "TrackProducerVar2FP420" ) {
     collector = finderParameters_->trackFinderVar2(input); //
     }// if ( trackMode
       */
       if (trackMode_ == "TrackProducerSophisticatedFP420" ) {
     collector = finderParameters_->trackFinderSophisticated(input,det); //
       }// if ( trackMode
       
       
       //  else if (trackMode_ == "TrackProducer3DFP420" ) {
       //     collector = finderParameters_->trackFinder3D(input); //
       // }// if ( trackMode
       
       if (collector.size()>0){
     TrackCollectionFP420::Range inputRange;
     inputRange.first = collector.begin();
     inputRange.second = collector.end();
     
     if ( first ) {
       // use it only if TrackCollectionFP420 is the TrackCollection of one event, otherwise, do not use (loose 1st cl. of 1st event only)
       first = false;
       unsigned int StID0 = 0;
       toutput->put(inputRange,StID0); // !!! put into adress 0 for detID which will not be used never
     } //if ( first ) 
     
     // !!! put                                        !!! put
     toutput->put(inputRange,StID);
     
     number_localelectroderechits += collector.size();
       } // if collector.size
     }//for det loop
     
     
     if (verbosity > 0) {
       std::cout << "FP420TrackMain: execution in mode " << trackMode_ << " generating " << number_localelectroderechits << " tracks in  " << number_detunits << " detectors" << std::endl; 
     }
     
     
     if (verbosity ==-29) {
       //     check of access to the collector:
       // loop over detunits
       for (int det=1; det<dn0; det++) {
     int StID = 1111;
     if(det==2) StID = 2222;
     std::vector<TrackFP420> collector;
     collector.clear();
     TrackCollectionFP420::Range outputRange;
     outputRange = toutput->get(StID);
     // fill output in collector vector (for may be sorting? or other checks)
     TrackCollectionFP420::ContainerIterator sort_begin = outputRange.first;
     TrackCollectionFP420::ContainerIterator sort_end = outputRange.second;
     for ( ;sort_begin != sort_end; ++sort_begin ) {
       collector.push_back(*sort_begin);
     } // for
     std::cout <<" ===" << std::endl;
     std::cout <<" ===" << std::endl;
     std::cout <<"=======FP420TrackMain:check size = " << collector.size() << "  det = " << det << std::endl;
     std::cout <<" ===" << std::endl;
     std::cout <<" ===" << std::endl;
     vector<TrackFP420>::const_iterator simHitIter = collector.begin();
     vector<TrackFP420>::const_iterator simHitIterEnd = collector.end();
     // loop in #tracks
     for (;simHitIter != simHitIterEnd; ++simHitIter) {
       const TrackFP420 itrack = *simHitIter;
       
       std::cout << "FP420TrackMain:check: nclusterx = " << itrack.nclusterx() << "  nclustery = " << itrack.nclustery() << std::endl;
       std::cout << "  ax = " << itrack.ax() << "  bx = " << itrack.bx() << std::endl;
       std::cout << "  ay = " << itrack.ay() << "  by = " << itrack.by() << std::endl;
       std::cout << " chi2x= " << itrack.chi2x() << " chi2y= " << itrack.chi2y() << std::endl;
       std::cout <<" ===" << std::endl;
       std::cout <<" ===" << std::endl;
       std::cout <<" =======================" << std::endl;
     }
     
     //==================================
     
     //     end of check of access to the strip collection
     std::cout <<"=======            FP420TrackMain:                    end of check     " << std::endl;
     
       }//for det
     }// if verbosity
     
     
     
     
   }// if ( validTrackerizer_
   
   
   
 }