DiagMaker Class Reference

Inheritance diagram for DiagMaker:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)
void	beginJob ()
	DiagMaker (const PSet &)
void	endJob ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Public Attributes
int	delayCount
int	evCntr

Private Attributes
TCanvas *	c1
TCanvas *	c2
WordHist *	wordBend
WordHist *	wordBx
WordHist *	wordBx0
WordHist *	wordCscID
WordHist *	wordKeywire
WordHist *	wordMpclink
WordHist *	wordPattern
WordHist *	wordPHD1
WordHist *	wordPHD2
WordHist *	wordQuality
WordHist *	wordRank
WordHist *	wordStrip
WordHist *	wordSyncErr
WordHist *	wordTHD1
WordHist *	wordTHD2
WordHist *	wordTrknmb
WordHist *	wordTrkPhi
WordHist *	wordTrkTheta
WordHist *	wordValid

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 106 of file DiagMaker.cc.

Constructor & Destructor Documentation

DiagMaker::DiagMaker

(

const PSet &

p

)

Definition at line 144 of file DiagMaker.cc.

{
evCntr = 0;   

}

Member Function Documentation

void DiagMaker::analyze ( const edm::Event &  ev, const edm::EventSetup &  es  )

virtual

Implements edm::EDAnalyzer.

Definition at line 151 of file DiagMaker.cc.

References gather_cfg::cout, edm::Event::getByLabel(), and L1TMuon::InternalTrack::kCSC.

{


    cout<<"Begin DiagMaker Analyzer:::::::::::::::::::::::::::::::::::::\n:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::\n\n";

  cout << "Reached analyze part of DiagMaker.\n";

  edm::Handle<InternalTrackCollection> emuITC;
  edm::Handle<InternalTrackCollection> dataITC;

  ev.getByLabel( "sptf" , "EmuITC"  , emuITC  );
  ev.getByLabel( "L1TMuonTrkFinder" , "DataITC"  , dataITC );

  int trknmbData, validData, qualityData, keywireData, stripData, patternData; 
  int bendData, bxData, mpclinkData, bx0Data, syncErrData, cscIDData;
  int trknmbEmu, validEmu, qualityEmu, keywireEmu, stripEmu, patternEmu;
  int bendEmu, bxEmu, mpclinkEmu, bx0Emu, syncErrEmu, cscIDEmu;

  auto itEmu     = emuITC->cbegin();
  auto itendEmu  = emuITC->cend();

  auto itData    = dataITC->cbegin();
  auto itendData = dataITC->cend(); 
  int trkCntr = 0;
  // Loop over all tracks in the collection
  for( ; (itEmu != itendEmu) && (itData != itendData) ; ++itEmu, ++itData)
  {
      trkCntr++;
      cout << "TRACK(Emu) "  << trkCntr << " Phi:  " << itEmu->phi  << ",   Theta:  " << itEmu ->theta << endl;
      cout << "TRACK(Data) " << trkCntr << " Phi:  " << itData->phi << ",   Theta:  " << itData->theta << endl;
      cout << "Track(Emu)  " << trkCntr << " Rank: " << itEmu->rank << endl;
      cout << "Track(Data) " << trkCntr << " Rank: " << itData->rank<< endl;
      cout << "Emu:\n dp1- "<<itEmu->deltas[0][0]<<" dp2- "<<itEmu->deltas[0][1]<<" dt1- "<<itEmu->deltas[1][0]<<" dt2- "<<itEmu->deltas[1][1]<<endl;
      cout << "Data:\n dp1- "<<itData->deltas[0][0]<<" dp2- "<<itData->deltas[0][1]<<" dt1- "<<itData->deltas[1][0]<<" dt2- "<<itData->deltas[1][1]<<endl;

      wordTrkPhi   -> fillWord( itData->phi          , itEmu->phi         );
      wordTrkTheta -> fillWord( itData->theta        , itEmu->theta       );
      wordRank     -> fillWord( itData->rank         , itEmu->rank        );
      wordPHD1     -> fillWord( fabs(itData->deltas[0][0]) , fabs(itEmu->deltas[0][0]));
      wordPHD2     -> fillWord( fabs(itData->deltas[0][1]) , fabs(itEmu->deltas[0][1]));
      wordTHD1     -> fillWord( fabs(itData->deltas[1][0]) , fabs(itEmu->deltas[1][0]));
      wordTHD2     -> fillWord( fabs(itData->deltas[1][1]) , fabs(itEmu->deltas[1][1]));
    
      if(itData->phi != itEmu->phi){std::cout<<"Track Phi Mismatch\n\n";}
      if(itData->theta != itEmu->theta){std::cout<<"Track Theta Mismatch\n\n";}
      if(itData->rank != itEmu->rank){std::cout<<"Track Rank Mismatch\n\n";}
      if(fabs(itData->deltas[0][0]) != fabs(itEmu->deltas[0][0])){std::cout<<"00 D Mismatch\n";}
      if(fabs(itData->deltas[0][1]) != fabs(itEmu->deltas[0][1])){std::cout<<"01 D Mismatch\n";}
      if(fabs(itData->deltas[1][0]) != fabs(itEmu->deltas[1][0])){std::cout<<"10 D Mismatch\n";}
      if(fabs(itData->deltas[1][1]) != fabs(itEmu->deltas[1][1])){std::cout<<"11 D Mismatch\n";}
      
      
      


      // Creating station maps for this track
      TriggerPrimitiveStationMap tpsmData = itData->getStubs();
      TriggerPrimitiveStationMap tpsmEmu  = itEmu->getStubs();

      // Getting the unique station ID number for ME1
      const unsigned id = 4*L1TMuon::InternalTrack::kCSC;

      // Looping over all four stations
      for(unsigned meNum=id; meNum<(id+4); meNum++)
      { 
          // Getting the trig prim lists for this station
          //TriggerPrimitiveList tplData = tpsmData[meNum];
          //TriggerPrimitiveList tplEmu  = tpsmEmu[meNum];
          std::vector<TriggerPrimitive> tplData = tpsmData[meNum];
          std::vector<TriggerPrimitive> tplEmu  = tpsmEmu[meNum];

          cout << "ME " << meNum-id+1 << " -  # Trig Prims = " << tplData.size() << endl;

          // Looping over all the trigger primitives in the lists
          for(unsigned tpNum = 0; (tpNum < tplData.size()) && (tpNum < tplEmu.size()) ; tpNum++)
          {
              cout << " ----- tp #" << tpNum << endl; 
              // Creating references to the trig prim info
              TriggerPrimitive tprData = tplData.at(tpNum);
              TriggerPrimitive tprEmu  = tplEmu.at(tpNum);

              trknmbData     = (tprData).getCSCData().trknmb;
              validData      = (tprData).getCSCData().valid;
              qualityData    = (tprData).getCSCData().quality;
              keywireData    = (tprData).getCSCData().keywire;
              stripData      = (tprData).getCSCData().strip;
              patternData    = (tprData).getCSCData().pattern;
              bendData       = (tprData).getCSCData().bend;
              bxData         = (tprData).getCSCData().bx;
              mpclinkData    = (tprData).getCSCData().mpclink;
              bx0Data        = (tprData).getCSCData().bx0;
              syncErrData    = (tprData).getCSCData().syncErr;
              cscIDData      = (tprData).getCSCData().cscID;
             
              trknmbEmu     = (tprEmu).getCSCData().trknmb;
              validEmu      = (tprEmu).getCSCData().valid;
              qualityEmu    = (tprEmu).getCSCData().quality;
              keywireEmu    = (tprEmu).getCSCData().keywire;
              stripEmu      = (tprEmu).getCSCData().strip;
              patternEmu    = (tprEmu).getCSCData().pattern;
              bendEmu       = (tprEmu).getCSCData().bend;
              bxEmu         = (tprEmu).getCSCData().bx;
              mpclinkEmu    = (tprEmu).getCSCData().mpclink;
              bx0Emu        = (tprEmu).getCSCData().bx0;
              syncErrEmu    = (tprEmu).getCSCData().syncErr;
              cscIDEmu      = (tprEmu).getCSCData().cscID;

              cout << "trknmb:   " << "data=" << trknmbData  << ",  emu=" << trknmbEmu  << endl;
              cout << "valid:    " << "data=" << validData   << ",  emu=" << validEmu   << endl;
              cout << "quality:  " << "data=" << qualityData << ",  emu=" << qualityEmu << endl;
              cout << "keywire:  " << "data=" << keywireData << ",  emu=" << keywireEmu << endl;
              cout << "strip:    " << "data=" << stripData   << ",  emu=" << stripEmu   << endl;
              cout << "pattern:  " << "data=" << patternData << ",  emu=" << patternEmu << endl;
              cout << "bend:     " << "data=" << bendData    << ",  emu=" << bendEmu    << endl;
              cout << "bx:       " << "data=" << bxData      << ",  emu=" << bxEmu      << endl;
              cout << "mpclink:  " << "data=" << mpclinkData << ",  emu=" << mpclinkEmu << endl;
              cout << "bx0:      " << "data=" << bx0Data     << ",  emu=" << bx0Emu     << endl;
              cout << "syncErr:  " << "data=" << syncErrData << ",  emu=" << syncErrEmu << endl;
              cout << "cscID:    " << "data=" << cscIDData   << ",  emu=" << cscIDEmu   << endl;
              cout << "--------------------------------\n";

              // Filling the word histograms
              wordTrknmb   -> fillWord( trknmbData  , trknmbEmu  );
              wordValid    -> fillWord( validData   , validEmu   );
              wordQuality  -> fillWord( qualityData , qualityEmu );
              wordKeywire  -> fillWord( keywireData , keywireEmu );
              wordStrip    -> fillWord( stripData   , stripEmu   );
              wordPattern  -> fillWord( patternData , patternEmu );
              wordBend     -> fillWord( bendData    , bendEmu    );
              wordBx       -> fillWord( bxData      , bxEmu      );
              wordMpclink  -> fillWord( mpclinkData , mpclinkEmu );
              wordBx0      -> fillWord( bx0Data     , bx0Emu     );
              wordSyncErr  -> fillWord( syncErrData , syncErrEmu );
              wordCscID    -> fillWord( cscIDData   , cscIDEmu   );
          
          if(keywireData != keywireEmu){std::cout<<"Wire Mismatch\n\n";}
          if(stripData != stripEmu){std::cout<<"Strip Mismatch\n\n";}
          

          }
      }
  }
  
  cout<<"End DiagMaker Analyzer:::::::::::::::::::::::::::::::::::::\n:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::\n\n";
  
}

void DiagMaker::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 300 of file DiagMaker.cc.

References alignmentValidation::c1.

{


    

    // Constructing all the word objects  
    c1 = new TCanvas();
    c1 -> cd();
    
    wordTrkPhi      = new WordHist( "Phi"     , 12 , 0   , 0   );
    c1 -> cd();
    wordTrkTheta    = new WordHist( "Theta"   ,  6 , 0   , 0.2 );
    c1 -> cd();
    wordTrknmb      = new WordHist( "Phhit"   ,  6 , 0   , 0.4 );
    c1 -> cd();
    wordValid       = new WordHist( "Valid"   , 12 , 0   , 0.6 );
    c1 -> cd();
    wordQuality     = new WordHist( "Quality" ,  7 , 0   , 0.8 );
    c1 -> cd();
    wordKeywire     = new WordHist( "Keywire" , 12 , 0.5 , 0   );
    c1 -> cd();
    wordStrip       = new WordHist( "Strip"   ,  7 , 0.5 , 0.2 );
    c1 -> cd();
    wordPattern     = new WordHist( "Pattern" , 12 , 0.5 , 0.4 );
    c1 -> cd();
    wordBend        = new WordHist( "Bend"    , 12 , 0.5 , 0.6 );
    c1 -> cd();
    wordBx          = new WordHist( "Bx"      , 12 , 0.5 , 0.8 );

    c2 = new TCanvas();
    c2 -> cd();
    wordMpclink     = new WordHist( "Mpclink" , 12 , 0   , 0   );
    c2 -> cd();
    wordBx0         = new WordHist( "Bx"      , 12 , 0   , 0.2 );
    c2 -> cd();
    wordSyncErr     = new WordHist( "SyncErr" , 12 , 0   , 0.4 );
    c2 -> cd();
    wordCscID       = new WordHist( "CscID"   , 12 , 0   , 0.6 );
    c2 -> cd();
    wordRank        = new WordHist( "Rank"    , 12 , 0.5 , 0   );
    c2 -> cd();
    wordPHD1        = new WordHist( "PHD1"    , 12 , 0.5 , 0.2 );
    c2 ->cd();
    wordPHD2        = new WordHist( "PHD2"    , 12 , 0.5 , 0.4 );
    c2 ->cd();
    wordTHD1        = new WordHist( "THD1"    , 12 , 0.5 , 0.6 );
    c2 ->cd();
    wordTHD2        = new WordHist( "THD2"    , 12 , 0.5 , 0.8 );

}

void DiagMaker::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 352 of file DiagMaker.cc.

References alignmentValidation::c1.

{
  
  c1 -> cd();

  wordTrkPhi   -> pad -> cd();
  wordTrkPhi   -> formatWordHist();
  wordTrkTheta -> pad -> cd();
  wordTrkTheta -> formatWordHist();
  wordTrknmb   -> pad -> cd();
  wordTrknmb   -> formatWordHist();
  wordValid    -> pad -> cd();
  wordValid    -> formatWordHist();
  wordQuality  -> pad -> cd();
  wordQuality  -> formatWordHist();
  wordKeywire  -> pad -> cd();
  wordKeywire  -> formatWordHist();
  wordStrip    -> pad -> cd();
  wordStrip    -> formatWordHist();
  wordPattern  -> pad -> cd();
  wordPattern  -> formatWordHist();
  wordBend     -> pad -> cd();
  wordBend     -> formatWordHist();
  wordBx       -> pad -> cd();
  wordBx       -> formatWordHist();

  gSystem->ProcessEvents();
  TImage *img1 = TImage::Create();
  img1 -> FromPad(c1);
  img1 -> WriteImage("wordHist1.png");
  delete c1;
  delete img1;

  c2 -> cd();

  wordMpclink  -> pad -> cd();
  wordMpclink  -> formatWordHist();
  wordBx0      -> pad -> cd();
  wordBx0      -> formatWordHist();
  wordSyncErr  -> pad -> cd();
  wordSyncErr  -> formatWordHist();
  wordCscID    -> pad -> cd();
  wordCscID    -> formatWordHist();
  wordRank     -> pad -> cd();
  wordRank     -> formatWordHist();
  wordPHD1     -> pad -> cd();
  wordPHD1     -> formatWordHist();
  wordPHD2     -> pad -> cd();
  wordPHD2     -> formatWordHist();
  wordTHD1     -> pad -> cd();
  wordTHD1     -> formatWordHist();
  wordTHD2     -> pad -> cd();
  wordTHD2     -> formatWordHist();

  gSystem->ProcessEvents();
  TImage *img2 = TImage::Create();
  img2->FromPad(c2);
  img2->WriteImage("wordHist2.png");
  delete c2;
  delete img2;
  
  

}

Member Data Documentation

TCanvas* DiagMaker::c1

private

Definition at line 116 of file DiagMaker.cc.

TCanvas * DiagMaker::c2

private

Definition at line 116 of file DiagMaker.cc.

int DiagMaker::delayCount

Definition at line 113 of file DiagMaker.cc.

int DiagMaker::evCntr

Definition at line 114 of file DiagMaker.cc.

WordHist* DiagMaker::wordBend

private

Definition at line 126 of file DiagMaker.cc.

WordHist* DiagMaker::wordBx

private

Definition at line 127 of file DiagMaker.cc.

WordHist* DiagMaker::wordBx0

private

Definition at line 129 of file DiagMaker.cc.

WordHist* DiagMaker::wordCscID

private

Definition at line 131 of file DiagMaker.cc.

WordHist* DiagMaker::wordKeywire

private

Definition at line 123 of file DiagMaker.cc.

WordHist* DiagMaker::wordMpclink

private

Definition at line 128 of file DiagMaker.cc.

WordHist* DiagMaker::wordPattern

private

Definition at line 125 of file DiagMaker.cc.

WordHist* DiagMaker::wordPHD1

private

Definition at line 133 of file DiagMaker.cc.

WordHist * DiagMaker::wordPHD2

private

Definition at line 133 of file DiagMaker.cc.

WordHist* DiagMaker::wordQuality

private

Definition at line 122 of file DiagMaker.cc.

WordHist* DiagMaker::wordRank

private

Definition at line 132 of file DiagMaker.cc.

WordHist* DiagMaker::wordStrip

private

Definition at line 124 of file DiagMaker.cc.

WordHist* DiagMaker::wordSyncErr

private

Definition at line 130 of file DiagMaker.cc.

WordHist* DiagMaker::wordTHD1

private

Definition at line 134 of file DiagMaker.cc.

WordHist * DiagMaker::wordTHD2

private

Definition at line 134 of file DiagMaker.cc.

WordHist* DiagMaker::wordTrknmb

private

Definition at line 120 of file DiagMaker.cc.

WordHist* DiagMaker::wordTrkPhi

private

Definition at line 119 of file DiagMaker.cc.

WordHist* DiagMaker::wordTrkTheta

private

Definition at line 118 of file DiagMaker.cc.

WordHist* DiagMaker::wordValid

private

Definition at line 121 of file DiagMaker.cc.