/data/refman/pasoursint/CMSSW_6_1_2_SLHC2/src/Geometry/CSCGeometry/interface/CSCChamberSpecs.h

Go to the documentation of this file.

#ifndef Geometry_CSCGeometry_CSCChamberSpecs_H
#define Geometry_CSCGeometry_CSCChamberSpecs_H

#include "Geometry/CommonDetUnit/interface/GeomDetType.h"

#include <cmath>
#include <string>
#include <vector>

class CSCGeometry;
class CSCLayerGeometry;
class CSCWireGroupPackage;
class Topology; //@@ Can go once GeomDetType has it removed.
class TrapezoidalPlaneBounds;

class CSCChamberSpecs : public GeomDetType { 

public:

  typedef std::vector<float> CSCSpecsParcel;

  CSCChamberSpecs();

  CSCChamberSpecs( const CSCGeometry* geom, int iChamberType, 
                   const TrapezoidalPlaneBounds& mediaShape,
                   const CSCSpecsParcel& fupar,
                   const CSCWireGroupPackage& wg 
                  );

  ~CSCChamberSpecs();

  bool operator!=( const CSCChamberSpecs& specs ) const;
  bool operator==( const CSCChamberSpecs& specs ) const;

  //@@ Topology() will be removed from GeomDetType interface (good, so we're no longer forced to pick something)
  //@@ But still there as of Aug-2007. So much for design.

  virtual const Topology& topology() const;

  const CSCLayerGeometry* oddLayerGeometry( int iendcap ) const 
   { return (iendcap==1? poszOddLayerGeometry:negzOddLayerGeometry);}
  const CSCLayerGeometry* evenLayerGeometry( int iendcap ) const 
   { return (iendcap==1? poszEvenLayerGeometry:negzEvenLayerGeometry);}

  int chamberType() const; 

  std::string chamberTypeName() const;


  // CHAMBER 'PARAMETER' FUNCTIONS
  //@@ FIXME these must all be sorted out in a world of real conditions & calibration

  float stripResolution()     const {return specsValue(0);}
  float wireResolution()      const {return specsValue(1);}
  float efficiency()          const {return specsValue(2);}
  float timeWindow()          const {return specsValue(3);}

  float neutronRate()         const {return specsValue(4);}

  int nStrips()               const {return nstrips;}

  int   nNodes()              const {return int( specsValue(6) );}
  int   nWiresPerGroup()      const {return int( specsValue(7) );}
  int   nFloatingStrips()     const {return int( specsValue(8) );}

  float stripPhiPitch() const {return 1.E-03 * stripDeltaPhi;}

  float ctiOffset() const {return centreToIntersectionOffset;}

  float wireSpacing() const;
 
  float anodeCathodeSpacing() const {return specsValue(9);}

  float gasGain()             const;

  float voltage()             const {return specsValue(11);}

  float calibrationError()    const {return specsValue(13);}

  float electronAttraction()  const {return specsValue(14);}

  float fractionQS()          const {return specsValue(15);}

  float chargePerCount() const;

  float wireRadius()          const {return specsValue(17);}

  float shaperPeakingTime()   const {return specsValue(18);}
  
  float constantNoise() const {return specsValue(22);}

  float e_pF() const {return specsValue(23);}

  //@@ Why 1/sqrt(time)?!?
  float wireNoise(float timeInterval) const {
    const float WIRE_CAPACITANCE = 8.; // in pF
    return (constantNoise() + nWiresPerGroup()*WIRE_CAPACITANCE)
           / sqrt(timeInterval/100.);
  }
 
  float stripNoise(float timeInterval) const;

  //@@ The following is nonsense to be fixed at some stage
  //  float adcThreshold()        const {return 9.99;}  

   bool gangedStrips() const { return gangedStrips_; }

  // STATIC FUNCTION

  static int whatChamberType( int istation, int iring );

 private:

  float specsValue( int index ) const {
    return theSpecsValues[ index ];  
  }

  // A ChamberSpecs has 4 associated LayerGeometry's
  CSCLayerGeometry* poszOddLayerGeometry;
  CSCLayerGeometry* poszEvenLayerGeometry;
  CSCLayerGeometry* negzOddLayerGeometry;
  CSCLayerGeometry* negzEvenLayerGeometry;

  //  theChamberType is a unique integer 1-10 for a station, ring pair.

  //  The type value is defined as <br>
  //        1           for S = 1  and R=A=4 split strips in ME11 <br>
  //      2,3,4 =  R+1  for S = 1  and R = 1,2,3 <br>
  //      5-10  = 2*S+R for S = 2,3,4 and R = 1,2 <br>
  
  int theChamberType;

  // Careful that the following order of data members is required to ensure
  // correct initialization in constructor initialization list!

  // Store for specs parameter values
  CSCSpecsParcel theSpecsValues;

  int nstrips; // no. of strips per layer
  float stripDeltaPhi;   // Delta-phi width of strip in this chamber type (in mrad)
  float centreToIntersectionOffset; // Possible correction to whereStripsMeet

  bool gangedStrips_;

  // Names of chamber types
  static const std::string theName[10];

  // Name of this class 
  static const std::string myName;

};

#endif // Geometry_CSCGeometry_CSCChamberSpecs_H