Inheritance diagram for SiStripGainFromCalibTree:

Classes
class	isEqual

Public Member Functions
	SiStripGainFromCalibTree (const edm::ParameterSet &)

	~SiStripGainFromCalibTree ()

Public Member Functions inherited from ConditionDBWriter< SiStripApvGain >
	ConditionDBWriter (const edm::ParameterSet &iConfig)

virtual	~ConditionDBWriter ()

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 (const std::string &iProcessName, std::vector< const char * > &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 ()

Private Member Functions
virtual void	algoAnalyze (const edm::Event &, const edm::EventSetup &) override

void	algoAnalyzeTheTree ()

virtual void	algoBeginRun (const edm::Run &run, const edm::EventSetup &iSetup) override

void	algoComputeMPVandGain ()

virtual void	algoEndJob () override

SiStripApvGain *	getNewObject () override

void	getPeakOfLandau (TH1 InputHisto, double FitResults, double LowRange=50, double HighRange=5400)

bool	IsFarFromBorder (TrajectoryStateOnSurface trajState, const uint32_t detid, const edm::EventSetup iSetup)

bool	IsGoodLandauFit (double *FitResults)

void	MakeCalibrationMap ()

bool	produceTagFilter ()

void	storeOnTree ()

Private Attributes
std::string	AlgoMode

bool	AllowSaturation

__gnu_cxx::hash_map< unsigned int, stAPVGain *, __gnu_cxx::hash< unsigned int >, isEqual >	APVsColl

std::vector< stAPVGain * >	APVsCollOrdered

unsigned int	BAD

int	CalibrationLevel

MonitorElement *	Charge_Vs_Index

MonitorElement *	Charge_Vs_Index_Absolute

MonitorElement *	Charge_Vs_PathlengthTECM1

MonitorElement *	Charge_Vs_PathlengthTECM2

MonitorElement *	Charge_Vs_PathlengthTECP1

MonitorElement *	Charge_Vs_PathlengthTECP2

MonitorElement *	Charge_Vs_PathlengthTIB

MonitorElement *	Charge_Vs_PathlengthTIDM

MonitorElement *	Charge_Vs_PathlengthTIDP

MonitorElement *	Charge_Vs_PathlengthTOB

DQMStore *	dbe

unsigned int	ERun

bool	FirstSetOfConstants

unsigned int	GOOD

bool	harvestingMode

string	m_calibrationPath

unsigned int	MASKED

double	MaxChi2OverNDF

double	MaxMPVError

unsigned int	MaxNrStrips

double	MaxTrackChiOverNdf

double	MaxTrackEta

double	MaxTrackMomentum

double	MinNrEntries

double	MinTrackEta

unsigned int	MinTrackHits

double	MinTrackMomentum

unsigned int	NCluster

unsigned int	NEvent

unsigned int	NTrack

bool	OldGainRemoving

std::string	OutputGains

unsigned int	SRun

double	tagCondition_GoodFrac

double	tagCondition_NClusters

TFileService *	tfs

edm::InputTag	theTracksLabel

bool	useCalibration

bool	Validation

vector< string >	VInputFiles

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 ConditionDBWriter< SiStripApvGain >
void	setDoStore (const bool doStore)
	When set to false the payload will not be written to the db. More...

void	storeOnDbNow ()

cond::Time_t	timeOfLastIOV ()

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 112 of file SiStripGainFromCalibTree.cc.

Constructor & Destructor Documentation

SiStripGainFromCalibTree::SiStripGainFromCalibTree ( const edm::ParameterSet & iConfig )

explicit

Definition at line 202 of file SiStripGainFromCalibTree.cc.

References AlgoMode, AllowSaturation, CalibrationLevel, dbe, FirstSetOfConstants, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), harvestingMode, m_calibrationPath, MaxChi2OverNDF, MaxMPVError, MaxNrStrips, MaxTrackChiOverNdf, MaxTrackEta, MaxTrackMomentum, MinNrEntries, MinTrackEta, MinTrackHits, MinTrackMomentum, OldGainRemoving, cppFunctionSkipper::operator, OutputGains, AlCaHLTBitMon_QueryRunRegistry::string, tagCondition_GoodFrac, tagCondition_NClusters, theTracksLabel, useCalibration, Validation, and VInputFiles.

                                                                                  : ConditionDBWriter<SiStripApvGain>(iConfig)
 {
    OutputGains         = iConfig.getParameter<std::string>("OutputGains");
    theTracksLabel      = iConfig.getUntrackedParameter<edm::InputTag>("Tracks");
 
    AlgoMode            = iConfig.getUntrackedParameter<std::string>("AlgoMode", "CalibTree");
    MinNrEntries        = iConfig.getUntrackedParameter<double>  ("minNrEntries"       ,  20);
    MaxMPVError         = iConfig.getUntrackedParameter<double>  ("maxMPVError"        ,  500.0);
    MaxChi2OverNDF      = iConfig.getUntrackedParameter<double>  ("maxChi2OverNDF"     ,  5.0);
    MinTrackMomentum    = iConfig.getUntrackedParameter<double>  ("minTrackMomentum"   ,  3.0);
    MaxTrackMomentum    = iConfig.getUntrackedParameter<double>  ("maxTrackMomentum"   ,  99999.0);
    MinTrackEta         = iConfig.getUntrackedParameter<double>  ("minTrackEta"        , -5.0);
    MaxTrackEta         = iConfig.getUntrackedParameter<double>  ("maxTrackEta"        ,  5.0);
    MaxNrStrips         = iConfig.getUntrackedParameter<unsigned>("maxNrStrips"        ,  2);
    MinTrackHits        = iConfig.getUntrackedParameter<unsigned>("MinTrackHits"       ,  8);
    MaxTrackChiOverNdf  = iConfig.getUntrackedParameter<double>  ("MaxTrackChiOverNdf" ,  3);
    AllowSaturation     = iConfig.getUntrackedParameter<bool>    ("AllowSaturation"    ,  false);
    FirstSetOfConstants = iConfig.getUntrackedParameter<bool>    ("FirstSetOfConstants",  true);
    Validation          = iConfig.getUntrackedParameter<bool>    ("Validation"         ,  false);
    OldGainRemoving     = iConfig.getUntrackedParameter<bool>    ("OldGainRemoving"    ,  false);
 
    CalibrationLevel    = iConfig.getUntrackedParameter<int>     ("CalibrationLevel"   ,  0);
    VInputFiles         = iConfig.getParameter<vector<string> >  ("InputFiles");
 
 
    useCalibration      = iConfig.getUntrackedParameter<bool>    ("UseCalibration"     , false);
    m_calibrationPath   = iConfig.getUntrackedParameter<string>  ("calibrationPath");
    harvestingMode      = iConfig.getUntrackedParameter<bool>    ("harvestingMode"     , false);
 
    tagCondition_NClusters  = iConfig.getUntrackedParameter<double>    ("NClustersForTagProd"     , 2E8);
    tagCondition_GoodFrac   = iConfig.getUntrackedParameter<double>    ("GoodFracForTagProd"     , 0.95);
 
    dbe = edm::Service<DQMStore>().operator->();
    dbe->setVerbose(10);
 }

SiStripGainFromCalibTree::~SiStripGainFromCalibTree ( )

Definition at line 778 of file SiStripGainFromCalibTree.cc.

779 {

780 }

Member Function Documentation

void SiStripGainFromCalibTree::algoAnalyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 804 of file SiStripGainFromCalibTree.cc.

References a, AlgoMode, AllowSaturation, SiStripCluster::amplitudes(), sistrip::APV, APVsColl, stAPVGain::CalibGain, Charge_Vs_Index, Charge_Vs_Index_Absolute, Charge_Vs_PathlengthTECM1, Charge_Vs_PathlengthTECM2, Charge_Vs_PathlengthTECP1, Charge_Vs_PathlengthTECP2, Charge_Vs_PathlengthTIB, Charge_Vs_PathlengthTIDM, Charge_Vs_PathlengthTIDP, Charge_Vs_PathlengthTOB, reco::TrackBase::chi2(), SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), edm::AssociationMap< edm::OneToOne< std::vector< Trajectory >, reco::TrackCollection, unsigned short > >::const_iterator, ERun, stAPVGain::Eta, reco::TrackBase::eta(), edm::false, MonitorElement::Fill(), FirstSetOfConstants, SiStripCluster::firstStrip(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByLabel(), h, harvestingMode, edm::EventBase::id(), stAPVGain::Index, IsFarFromBorder(), edm::ESHandleBase::isValid(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localDirection(), PV3DBase< T, PVType, FrameType >::mag(), MaxNrStrips, MaxTrackChiOverNdf, MaxTrackEta, MaxTrackMomentum, Trajectory::measurements(), MinTrackEta, MinTrackHits, MinTrackMomentum, SiStripMatchedRecHit2D::monoCluster(), SiStripMatchedRecHit2D::monoId(), NCluster, reco::TrackBase::ndof(), NEvent, NTrack, NULL, reco::TrackBase::numberOfValidHits(), OldGainRemoving, reco::TrackBase::p(), CSCOverlapsAlignmentAlgorithm_cff::Path, stAPVGain::PreviousGain, DetId::rawId(), edm::EventID::run(), alignCSCRings::s, SRun, SiStripMatchedRecHit2D::stereoCluster(), SiStripMatchedRecHit2D::stereoId(), stAPVGain::SubDet, StripSubdetector::TEC, theTracksLabel, stAPVGain::Thickness, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, testEve_cfg::tracks, useCalibration, Validation, and PV3DBase< T, PVType, FrameType >::z().

 {
   // in AlCaHarvesting mode we just need to run the logic in the endJob step
   if(harvestingMode) return;
   
    if(AlgoMode=="CalibTree")return;
 
    if(NEvent==0){
       SRun          = iEvent.id().run();
    }
    ERun             = iEvent.id().run();
    NEvent++;
 
   //FROM SHALLOW GAIN
   edm::ESHandle<TrackerGeometry> theTrackerGeometry;         iSetup.get<TrackerDigiGeometryRecord>().get( theTrackerGeometry );  
   edm::ESHandle<SiStripGain> gainHandle;                     iSetup.get<SiStripGainRcd>().get(gainHandle);
   edm::Handle<edm::View<reco::Track> > tracks;               iEvent.getByLabel(theTracksLabel, tracks);    
   edm::Handle<TrajTrackAssociationCollection> associations;  iEvent.getByLabel(theTracksLabel, associations);
 
   for( TrajTrackAssociationCollection::const_iterator association = associations->begin(); association != associations->end(); association++) {
        const Trajectory*  traj  = association->key.get();
        const reco::Track* track = association->val.get();
 
        //clean on the tracks
        if(fabs(track->eta())          < MinTrackEta        )continue;
        if(fabs(track->eta())          > MaxTrackEta        )continue;
        if(track->p()                  < MinTrackMomentum   )continue;
        if(track->p()                  > MaxTrackMomentum   )continue;
        if(track->numberOfValidHits()  < MinTrackHits       )continue;
        if(track->chi2()/track->ndof() > MaxTrackChiOverNdf )continue;
        NTrack++;
 
        vector<TrajectoryMeasurement> measurements = traj->measurements();
        for(vector<TrajectoryMeasurement>::const_iterator measurement_it = measurements.begin(); measurement_it!=measurements.end(); measurement_it++){
           TrajectoryStateOnSurface trajState = measurement_it->updatedState();
           if( !trajState.isValid() ) continue;     
 
           const TrackingRecHit*         hit                 = (*measurement_it->recHit()).hit();
           const SiStripRecHit1D*        sistripsimple1dhit  = dynamic_cast<const SiStripRecHit1D*>(hit);
           const SiStripRecHit2D*        sistripsimplehit    = dynamic_cast<const SiStripRecHit2D*>(hit);
           const SiStripMatchedRecHit2D* sistripmatchedhit   = dynamic_cast<const SiStripMatchedRecHit2D*>(hit);
 
           const SiStripCluster*   Cluster = NULL;
           uint32_t                DetId = 0;
 
           for(unsigned int h=0;h<2;h++){
             if(!sistripmatchedhit && h==1){
                continue;
             }else if(sistripmatchedhit  && h==0){
                Cluster = &sistripmatchedhit->monoCluster();
                DetId = sistripmatchedhit->monoId();
             }else if(sistripmatchedhit  && h==1){
                Cluster = &sistripmatchedhit->stereoCluster();;
                DetId = sistripmatchedhit->stereoId();
             }else if(sistripsimplehit){
                Cluster = (sistripsimplehit->cluster()).get();
                DetId = sistripsimplehit->geographicalId().rawId();
             }else if(sistripsimple1dhit){
                Cluster = (sistripsimple1dhit->cluster()).get();
                DetId = sistripsimple1dhit->geographicalId().rawId();
             }else{
                continue;
             }
 
             LocalVector             trackDirection = trajState.localDirection();
             double                  cosine         = trackDirection.z()/trackDirection.mag();
             const auto&             Ampls          = Cluster->amplitudes();
             int                     FirstStrip     = Cluster->firstStrip();
             int                     APVId          = FirstStrip/128;
             bool                    Saturation     = false;
             bool                    Overlapping    = false;
             unsigned int            charge         = 0;
             double                  PrevGain       = -1;
             stAPVGain* APV = APVsColl[(DetId<<3) | (FirstStrip/128)];
             double                  Path           = (10.0*APV->Thickness)/fabs(cosine);
 
 
             if(gainHandle.isValid()){ 
                SiStripApvGain::Range detGainRange = gainHandle->getRange(DetId);
                PrevGain = *(detGainRange.first + APVId);
             }
 
             for(unsigned int a=0;a<Ampls.size();a++){               
                charge+=Ampls[a];
                if(Ampls[a] >=254)Saturation =true;
             }
 
             if(FirstStrip==0                                  )Overlapping=true;
             if(FirstStrip==128                                )Overlapping=true;
             if(FirstStrip==256                                )Overlapping=true;
             if(FirstStrip==384                                )Overlapping=true;
             if(FirstStrip==512                                )Overlapping=true;
             if(FirstStrip==640                                )Overlapping=true;
 
             if(FirstStrip<=127 && FirstStrip+Ampls.size()>127)Overlapping=true;
             if(FirstStrip<=255 && FirstStrip+Ampls.size()>255)Overlapping=true;
             if(FirstStrip<=383 && FirstStrip+Ampls.size()>383)Overlapping=true;
             if(FirstStrip<=511 && FirstStrip+Ampls.size()>511)Overlapping=true;
             if(FirstStrip<=639 && FirstStrip+Ampls.size()>639)Overlapping=true;
 
             if(FirstStrip+Ampls.size()==127                   )Overlapping=true;
             if(FirstStrip+Ampls.size()==255                   )Overlapping=true;
             if(FirstStrip+Ampls.size()==383                   )Overlapping=true;
             if(FirstStrip+Ampls.size()==511                   )Overlapping=true;
             if(FirstStrip+Ampls.size()==639                   )Overlapping=true;
             if(FirstStrip+Ampls.size()==767                   )Overlapping=true;
 
             //cleaning on the cluster
             if(IsFarFromBorder(&trajState,DetId, &iSetup)  == false           )continue;
             if(Overlapping                                 == true            )continue;
             if(Saturation  && !AllowSaturation                                )continue;
             if(Ampls.size()                              > MaxNrStrips        )continue;
             NCluster++;
 
 
             int Charge = 0;
             if(useCalibration || !FirstSetOfConstants){
                bool Saturation = false;
                for(unsigned int s=0;s<Ampls.size();s++){
                   int StripCharge =  Ampls[s];
                   if(useCalibration && !FirstSetOfConstants){ StripCharge=(int)(StripCharge*(APV->PreviousGain/APV->CalibGain));
                   }else if(useCalibration){                   StripCharge=(int)(StripCharge/APV->CalibGain);
                   }else if(!FirstSetOfConstants){             StripCharge=(int)(StripCharge*APV->PreviousGain);}
                   if(StripCharge>1024){
                      StripCharge = 255;
                      Saturation = true;
                   }else if(StripCharge>254){
                      StripCharge = 254;
                      Saturation = true;
                   }
                   Charge += StripCharge;
                }
                if(Saturation && !AllowSaturation)continue;
             }else{
                Charge = charge;
             }
 
             //printf("ChargeDifference = %i Vs %i with Gain = %f\n",(*charge)[i],Charge,APV->CalibGain);
 
             double ClusterChargeOverPath   =  ( (double) Charge )/Path ;
             if(Validation)     {ClusterChargeOverPath/=PrevGain;}
             if(OldGainRemoving){ClusterChargeOverPath*=PrevGain;}
             Charge_Vs_Index_Absolute->Fill(APV->Index,Charge);   
             Charge_Vs_Index         ->Fill(APV->Index,ClusterChargeOverPath);
 
 
                   if(APV->SubDet==StripSubdetector::TIB){ Charge_Vs_PathlengthTIB  ->Fill(Path,Charge);
             }else if(APV->SubDet==StripSubdetector::TOB){ Charge_Vs_PathlengthTOB  ->Fill(Path,Charge);
             }else if(APV->SubDet==StripSubdetector::TID){ 
                      if(APV->Eta<0){                      Charge_Vs_PathlengthTIDM ->Fill(Path,Charge);
                }else if(APV->Eta>0){                      Charge_Vs_PathlengthTIDP ->Fill(Path,Charge);
                }
             }else if(APV->SubDet==StripSubdetector::TEC){
                      if(APV->Eta<0){
                         if(APV->Thickness<0.04){          Charge_Vs_PathlengthTECM1->Fill(Path,Charge);
                   }else if(APV->Thickness>0.04){          Charge_Vs_PathlengthTECM2->Fill(Path,Charge);
                   }
                }else if(APV->Eta>0){
                         if(APV->Thickness<0.04){          Charge_Vs_PathlengthTECP1->Fill(Path,Charge);
                   }else if(APV->Thickness>0.04){          Charge_Vs_PathlengthTECP2->Fill(Path,Charge);
                   }
                }
             }
 
           }//loop on  clusters
        }//loop on measurements
   }//loop on tracks
 
 }

void SiStripGainFromCalibTree::algoAnalyzeTheTree ( )

private

Definition at line 413 of file SiStripGainFromCalibTree.cc.

Referenced by algoEndJob().

 {
    for(unsigned int i=0;i<VInputFiles.size();i++){
       printf("Openning file %3i/%3i --> %s\n",i+1, (int)VInputFiles.size(), (char*)(VInputFiles[i].c_str())); fflush(stdout);
       TChain* tree = new TChain("gainCalibrationTree/tree");
       tree->Add(VInputFiles[i].c_str());
 
       TString EventPrefix("");
       TString EventSuffix("");
 
       TString TrackPrefix("track");
       TString TrackSuffix("");
 
       TString CalibPrefix("GainCalibration");
       TString CalibSuffix("");
 
       unsigned int                 eventnumber    = 0;    tree->SetBranchAddress(EventPrefix + "event"          + EventSuffix, &eventnumber   , NULL);
       unsigned int                 runnumber      = 0;    tree->SetBranchAddress(EventPrefix + "run"            + EventSuffix, &runnumber     , NULL);
       std::vector<bool>*           TrigTech    = 0;    tree->SetBranchAddress(EventPrefix + "TrigTech"    + EventSuffix, &TrigTech   , NULL);
 
       std::vector<double>*         trackchi2ndof  = 0;    tree->SetBranchAddress(TrackPrefix + "chi2ndof"       + TrackSuffix, &trackchi2ndof , NULL);
       std::vector<float>*          trackp         = 0;    tree->SetBranchAddress(TrackPrefix + "momentum"       + TrackSuffix, &trackp        , NULL);
       std::vector<float>*          trackpt        = 0;    tree->SetBranchAddress(TrackPrefix + "pt"             + TrackSuffix, &trackpt       , NULL);
       std::vector<double>*         tracketa       = 0;    tree->SetBranchAddress(TrackPrefix + "eta"            + TrackSuffix, &tracketa      , NULL);
       std::vector<double>*         trackphi       = 0;    tree->SetBranchAddress(TrackPrefix + "phi"            + TrackSuffix, &trackphi      , NULL);
       std::vector<unsigned int>*   trackhitsvalid = 0;    tree->SetBranchAddress(TrackPrefix + "hitsvalid"      + TrackSuffix, &trackhitsvalid, NULL);
 
       std::vector<int>*            trackindex     = 0;    tree->SetBranchAddress(CalibPrefix + "trackindex"     + CalibSuffix, &trackindex    , NULL);
       std::vector<unsigned int>*   rawid          = 0;    tree->SetBranchAddress(CalibPrefix + "rawid"          + CalibSuffix, &rawid         , NULL);
       std::vector<float>*          localdirx      = 0;    tree->SetBranchAddress(CalibPrefix + "localdirx"      + CalibSuffix, &localdirx     , NULL);
       std::vector<float>*          localdiry      = 0;    tree->SetBranchAddress(CalibPrefix + "localdiry"      + CalibSuffix, &localdiry     , NULL);
       std::vector<float>*          localdirz      = 0;    tree->SetBranchAddress(CalibPrefix + "localdirz"      + CalibSuffix, &localdirz     , NULL);
       std::vector<unsigned short>* firststrip     = 0;    tree->SetBranchAddress(CalibPrefix + "firststrip"     + CalibSuffix, &firststrip    , NULL);
       std::vector<unsigned short>* nstrips        = 0;    tree->SetBranchAddress(CalibPrefix + "nstrips"        + CalibSuffix, &nstrips       , NULL);
       std::vector<bool>*           saturation     = 0;    tree->SetBranchAddress(CalibPrefix + "saturation"     + CalibSuffix, &saturation    , NULL);
       std::vector<bool>*           overlapping    = 0;    tree->SetBranchAddress(CalibPrefix + "overlapping"    + CalibSuffix, &overlapping   , NULL);
       std::vector<bool>*           farfromedge    = 0;    tree->SetBranchAddress(CalibPrefix + "farfromedge"    + CalibSuffix, &farfromedge   , NULL);
       std::vector<unsigned int>*   charge         = 0;    tree->SetBranchAddress(CalibPrefix + "charge"         + CalibSuffix, &charge        , NULL);
       std::vector<float>*          path           = 0;    tree->SetBranchAddress(CalibPrefix + "path"           + CalibSuffix, &path          , NULL);
       std::vector<float>*          chargeoverpath = 0;    tree->SetBranchAddress(CalibPrefix + "chargeoverpath" + CalibSuffix, &chargeoverpath, NULL);
       std::vector<unsigned char>*  amplitude      = 0;    tree->SetBranchAddress(CalibPrefix + "amplitude"      + CalibSuffix, &amplitude     , NULL);
       std::vector<double>*         gainused       = 0;    tree->SetBranchAddress(CalibPrefix + "gainused"       + CalibSuffix, &gainused      , NULL);
 
 
       printf("Number of Events = %i + %i = %i\n",NEvent,(unsigned int)tree->GetEntries(),(unsigned int)(NEvent+tree->GetEntries()));
       printf("Progressing Bar              :0%%       20%%       40%%       60%%       80%%       100%%\n");
       printf("Looping on the Tree          :");
       int TreeStep = tree->GetEntries()/50;if(TreeStep<=1)TreeStep=1;
       for (unsigned int ientry = 0; ientry < tree->GetEntries(); ientry++) {
       if(ientry%TreeStep==0){printf(".");fflush(stdout);}
          tree->GetEntry(ientry);
 
          if(runnumber<SRun)SRun=runnumber;
          if(runnumber>ERun)ERun=runnumber;
 
          NEvent++;
          NTrack+=(*trackp).size();
 
      unsigned int FirstAmplitude=0;
          for(unsigned int i=0;i<(*chargeoverpath).size();i++){
             FirstAmplitude+=(*nstrips)[i];
                  int    TI = (*trackindex)[i];
             if((*tracketa      )[TI]  < MinTrackEta        )continue;
             if((*tracketa      )[TI]  > MaxTrackEta        )continue;
             if((*trackp        )[TI]  < MinTrackMomentum   )continue;
             if((*trackp        )[TI]  > MaxTrackMomentum   )continue;
             if((*trackhitsvalid)[TI]  < MinTrackHits       )continue;
             if((*trackchi2ndof )[TI]  > MaxTrackChiOverNdf )continue;
             if((*farfromedge)[i]        == false           )continue;
             if((*overlapping)[i]        == true            )continue;
             if((*saturation )[i]        && !AllowSaturation)continue;
             if((*nstrips    )[i]      > MaxNrStrips        )continue;
             NCluster++;
 
             stAPVGain* APV = APVsColl[((*rawid)[i]<<3) | ((*firststrip)[i]/128)];
             //printf("detId=%7i run=%7i event=%9i charge=%5i cs=%3i\n",(*rawid)[i],runnumber,eventnumber,(*charge)[i],(*nstrips)[i]);
 
             //double trans = atan2((*localdiry)[i],(*localdirx)[i])*(180/3.14159265);
             //double alpha = acos ((*localdirx)[i] / sqrt( pow((*localdirx)[i],2) +  pow((*localdirz)[i],2) ) ) * (180/3.14159265);
             //double beta  = acos ((*localdiry)[i] / sqrt( pow((*localdirx)[i],2) +  pow((*localdirz)[i],2) ) ) * (180/3.14159265);
       
         //printf("NStrip = %i : Charge = %i --> Path = %f  --> ChargeOverPath=%f\n",(*nstrips)[i],(*charge)[i],(*path)[i],(*chargeoverpath)[i]);
             //printf("Amplitudes: ");
             //for(unsigned int a=0;a<(*nstrips)[i];a++){printf("%i ",(*amplitude)[FirstAmplitude+a]);}
             //printf("\n");
 
             
             int Charge = 0;
             if(useCalibration || !FirstSetOfConstants){
                bool Saturation = false;
                for(unsigned int s=0;s<(*nstrips)[i];s++){
                   int StripCharge =  (*amplitude)[FirstAmplitude-(*nstrips)[i]+s];
                   if(useCalibration && !FirstSetOfConstants){ StripCharge=(int)(StripCharge*(APV->PreviousGain/APV->CalibGain));
                   }else if(useCalibration){                   StripCharge=(int)(StripCharge/APV->CalibGain);
                   }else if(!FirstSetOfConstants){             StripCharge=(int)(StripCharge*APV->PreviousGain);}
                   if(StripCharge>1024){
                      StripCharge = 255;
                      Saturation = true;
                   }else if(StripCharge>254){
                      StripCharge = 254;
                      Saturation = true;
                   }
                   Charge += StripCharge;
                }
                if(Saturation && !AllowSaturation)continue;
             }else{
                Charge = (*charge)[i];
             }
 
             //printf("ChargeDifference = %i Vs %i with Gain = %f\n",(*charge)[i],Charge,APV->CalibGain);
 
             double ClusterChargeOverPath   =  ( (double) Charge )/(*path)[i] ;
             if(Validation)     {ClusterChargeOverPath/=(*gainused)[i];}
             if(OldGainRemoving){ClusterChargeOverPath*=(*gainused)[i];}
             Charge_Vs_Index_Absolute->Fill(APV->Index,Charge);   
             Charge_Vs_Index         ->Fill(APV->Index,ClusterChargeOverPath);
 
 
               if(APV->SubDet==StripSubdetector::TIB){ Charge_Vs_PathlengthTIB  ->Fill((*path)[i],Charge); 
             }else if(APV->SubDet==StripSubdetector::TOB){ Charge_Vs_PathlengthTOB  ->Fill((*path)[i],Charge);
             }else if(APV->SubDet==StripSubdetector::TID){
                      if(APV->Eta<0){              Charge_Vs_PathlengthTIDM ->Fill((*path)[i],Charge);
                }else if(APV->Eta>0){                      Charge_Vs_PathlengthTIDP ->Fill((*path)[i],Charge);
                }
             }else if(APV->SubDet==StripSubdetector::TEC){
                      if(APV->Eta<0){
                         if(APV->Thickness<0.04){          Charge_Vs_PathlengthTECM1->Fill((*path)[i],Charge);
                   }else if(APV->Thickness>0.04){          Charge_Vs_PathlengthTECM2->Fill((*path)[i],Charge);
                   }
                }else if(APV->Eta>0){
                         if(APV->Thickness<0.04){          Charge_Vs_PathlengthTECP1->Fill((*path)[i],Charge);
                   }else if(APV->Thickness>0.04){          Charge_Vs_PathlengthTECP2->Fill((*path)[i],Charge);
                   }
                }
             }
 
          }// END OF ON-CLUSTER LOOP
       }printf("\n");// END OF EVENT LOOP
 
    }
 }

void SiStripGainFromCalibTree::algoBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 238 of file SiStripGainFromCalibTree.cc.

 {
   cout << "algoBeginRun start" << endl;
   if(!harvestingMode){
      cout << "   booking start" << endl;
      dbe->setCurrentFolder("AlCaReco/SiStripGains/");
      Charge_Vs_Index           = dbe->book2D("Charge_Vs_Index"          , "Charge_Vs_Index"          , 72785, 0   , 72784,1000,0,2000);
      Charge_Vs_Index_Absolute  = dbe->book2D("Charge_Vs_Index_Absolute" , "Charge_Vs_Index_Absolute" , 72785, 0   , 72784, 500,0,2000);
      Charge_Vs_PathlengthTIB   = dbe->book2D("Charge_Vs_PathlengthTIB"  , "Charge_Vs_PathlengthTIB"  , 20   , 0.3 , 1.3  , 250,0,2000);
      Charge_Vs_PathlengthTOB   = dbe->book2D("Charge_Vs_PathlengthTOB"  , "Charge_Vs_PathlengthTOB"  , 20   , 0.3 , 1.3  , 250,0,2000);
      Charge_Vs_PathlengthTIDP  = dbe->book2D("Charge_Vs_PathlengthTIDP" , "Charge_Vs_PathlengthTIDP" , 20   , 0.3 , 1.3  , 250,0,2000);
      Charge_Vs_PathlengthTIDM  = dbe->book2D("Charge_Vs_PathlengthTIDM" , "Charge_Vs_PathlengthTIDM" , 20   , 0.3 , 1.3  , 250,0,2000);
      Charge_Vs_PathlengthTECP1 = dbe->book2D("Charge_Vs_PathlengthTECP1", "Charge_Vs_PathlengthTECP1", 20   , 0.3 , 1.3  , 250,0,2000);
      Charge_Vs_PathlengthTECP2 = dbe->book2D("Charge_Vs_PathlengthTECP2", "Charge_Vs_PathlengthTECP2", 20   , 0.3 , 1.3  , 250,0,2000);
      Charge_Vs_PathlengthTECM1 = dbe->book2D("Charge_Vs_PathlengthTECM1", "Charge_Vs_PathlengthTECM1", 20   , 0.3 , 1.3  , 250,0,2000);
      Charge_Vs_PathlengthTECM2 = dbe->book2D("Charge_Vs_PathlengthTECM2", "Charge_Vs_PathlengthTECM2", 20   , 0.3 , 1.3  , 250,0,2000);
    } 
    edm::ESHandle<TrackerGeometry> tkGeom;
    iSetup.get<TrackerDigiGeometryRecord>().get( tkGeom );
    auto const & Det = tkGeom->dets();
 
    edm::ESHandle<SiStripGain> gainHandle;
    iSetup.get<SiStripGainRcd>().get(gainHandle);
    if(!gainHandle.isValid()){printf("\n#####################\n\nERROR --> gainHandle is not valid\n\n#####################\n\n");exit(0);}
  
    for(unsigned int i=0;i<gainHandle->getNumberOfTags();i++){
       printf("Reccord %i --> Rcd Name = %s    Label Name = %s\n",i,gainHandle->getRcdName(i).c_str(), gainHandle->getLabelName(i).c_str());
    }
  
    edm::ESHandle<SiStripQuality> SiStripQuality_;
    iSetup.get<SiStripQualityRcd>().get(SiStripQuality_);
 
    unsigned int Index=0;
    for(unsigned int i=0;i<Det.size();i++){
       DetId  Detid  = Det[i]->geographicalId(); 
       int    SubDet = Detid.subdetId();
 
       if( SubDet == StripSubdetector::TIB ||  SubDet == StripSubdetector::TID ||
           SubDet == StripSubdetector::TOB ||  SubDet == StripSubdetector::TEC  ){
 
           auto DetUnit     = dynamic_cast<const StripGeomDetUnit*> (Det[i]);
       if(!DetUnit)continue;
 
           const StripTopology& Topo     = DetUnit->specificTopology();  
           unsigned int         NAPV     = Topo.nstrips()/128;
           for(unsigned int j=0;j<NAPV;j++){
                 stAPVGain* APV = new stAPVGain;
                 APV->Index         = Index;
                 APV->Bin           = -1;
                 APV->DetId         = Detid.rawId();
                 APV->APVId         = j;
                 APV->SubDet        = SubDet;
                 APV->FitMPV        = -1;
                 APV->FitMPVErr     = -1;
                 APV->FitWidth      = -1;
                 APV->FitWidthErr   = -1;
                 APV->FitChi2       = -1;
                 APV->Gain          = -1;
                 APV->PreviousGain  = 1;
                 APV->x             = DetUnit->position().basicVector().x();
                 APV->y             = DetUnit->position().basicVector().y();
                 APV->z             = DetUnit->position().basicVector().z();
                 APV->Eta           = DetUnit->position().basicVector().eta();
                 APV->Phi           = DetUnit->position().basicVector().phi();
                 APV->R             = DetUnit->position().basicVector().transverse();
                 APV->Thickness     = DetUnit->surface().bounds().thickness();
         APV->NEntries      = 0;
                 APV->isMasked      = SiStripQuality_->IsApvBad(Detid.rawId(),j);
 
         if(!FirstSetOfConstants){
            if(gainHandle->getNumberOfTags()!=2){printf("ERROR: NUMBER OF GAIN TAG IS EXPECTED TO BE 2\n");fflush(stdout);exit(0);};        
            APV->PreviousGain  = gainHandle->getApvGain(APV->APVId,gainHandle->getRange(APV->DetId, 1),1);
                    //printf("DETID = %7i APVID=%1i Previous Gain=%8.4f\n",APV->DetId,APV->APVId,APV->PreviousGain);
         }
 
                 APVsCollOrdered.push_back(APV);
         APVsColl[(APV->DetId<<3) | APV->APVId] = APV;
                 Index++;
           }
       }
    }
 
    MakeCalibrationMap();
 
    NEvent     = 0;
    NTrack     = 0;
    NCluster   = 0;
    SRun       = 1<<31;
    ERun       = 0;
    GOOD       = 0;
    BAD        = 0;
    MASKED     = 0;
    
    cout << "algoBeginRun end" << endl;
 }

void SiStripGainFromCalibTree::algoComputeMPVandGain ( )

private

Definition at line 557 of file SiStripGainFromCalibTree.cc.

References sistrip::APV, stAPVGain::APVId, APVsColl, BAD, stAPVGain::Bin, CalibrationLevel, Charge_Vs_Index, stAPVGain::DetId, stAPVGain::FitChi2, stAPVGain::FitMPV, stAPVGain::FitMPVErr, stAPVGain::FitWidth, stAPVGain::FitWidthErr, stAPVGain::Gain, getPeakOfLandau(), MonitorElement::getTH2F(), GOOD, Exhume::I, i, stAPVGain::Index, IsGoodLandauFit(), stAPVGain::isMasked, MASKED, stAPVGain::NEntries, NULL, and stAPVGain::PreviousGain.

Referenced by algoEndJob().

                                                      {
    unsigned int I=0;
    TH1F* Proj = NULL;
    double FitResults[5];
    double MPVmean = 300;
 
    TH2F *chvsidx = Charge_Vs_Index->getTH2F();
 
 
    printf("Progressing Bar              :0%%       20%%       40%%       60%%       80%%       100%%\n");
    printf("Fitting Charge Distribution  :");
    int TreeStep = APVsColl.size()/50;
    for(__gnu_cxx::hash_map<unsigned int, stAPVGain*,  __gnu_cxx::hash<unsigned int>, isEqual >::iterator it = APVsColl.begin();it!=APVsColl.end();it++,I++){
    if(I%TreeStep==0){printf(".");fflush(stdout);}
       stAPVGain* APV = it->second;
       if(APV->Bin<0) APV->Bin = chvsidx->GetXaxis()->FindBin(APV->Index);
 
       if(APV->isMasked){MASKED++; continue;}
 
       Proj = (TH1F*)(chvsidx->ProjectionY("",chvsidx->GetXaxis()->FindBin(APV->Index),chvsidx->GetXaxis()->FindBin(APV->Index),"e"));
       if(!Proj)continue;
 
       if(CalibrationLevel==0){
       }else if(CalibrationLevel==1){
          int SecondAPVId = APV->APVId;
          if(SecondAPVId%2==0){    SecondAPVId = SecondAPVId+1; }else{ SecondAPVId = SecondAPVId-1; }
      stAPVGain* APV2 = APVsColl[(APV->DetId<<3) | SecondAPVId];
          if(APV2->Bin<0) APV2->Bin = chvsidx->GetXaxis()->FindBin(APV2->Index);
          TH1F* Proj2 = (TH1F*)(chvsidx->ProjectionY("",APV2->Bin,APV2->Bin,"e"));
          if(Proj2){Proj->Add(Proj2,1);delete Proj2;}
       }else if(CalibrationLevel==2){
           for(unsigned int i=0;i<6;i++){
             __gnu_cxx::hash_map<unsigned int, stAPVGain*,  __gnu_cxx::hash<unsigned int>, isEqual >::iterator tmpit;
             tmpit = APVsColl.find((APV->DetId<<3) | i);
             if(tmpit==APVsColl.end())continue;
             stAPVGain* APV2 = tmpit->second;
         if(APV2->DetId != APV->DetId || APV2->APVId == APV->APVId)continue;            
             if(APV2->Bin<0) APV2->Bin = chvsidx->GetXaxis()->FindBin(APV2->Index);
             TH1F* Proj2 = (TH1F*)(chvsidx->ProjectionY("",APV2->Bin,APV2->Bin,"e"));
             if(Proj2){Proj->Add(Proj2,1);delete Proj2;}
           }          
       }else{
          CalibrationLevel = 0;
          printf("Unknown Calibration Level, will assume %i\n",CalibrationLevel);
       }
 
       getPeakOfLandau(Proj,FitResults);
       APV->FitMPV      = FitResults[0];
       APV->FitMPVErr   = FitResults[1];
       APV->FitWidth    = FitResults[2];
       APV->FitWidthErr = FitResults[3];
       APV->FitChi2     = FitResults[4];
       APV->NEntries    = Proj->GetEntries();
 
       if(IsGoodLandauFit(FitResults)){
           APV->Gain = APV->FitMPV / MPVmean;
           GOOD++;
       }else{
           APV->Gain = APV->PreviousGain;
           BAD++;
       }
       if(APV->Gain<=0)           APV->Gain  = 1;
 
       //printf("%5i/%5i:  %6i - %1i  %5E Entries --> MPV = %f +- %f\n",I,APVsColl.size(),APV->DetId, APV->APVId, Proj->GetEntries(), FitResults[0], FitResults[1]);fflush(stdout);
       delete Proj;
    }printf("\n");
 }

void SiStripGainFromCalibTree::algoEndJob ( )

overrideprivatevirtual

Reimplemented from ConditionDBWriter< SiStripApvGain >.

Definition at line 335 of file SiStripGainFromCalibTree.cc.

References algoAnalyzeTheTree(), algoComputeMPVandGain(), AlgoMode, Charge_Vs_Index, Charge_Vs_Index_Absolute, Charge_Vs_PathlengthTECM1, Charge_Vs_PathlengthTECM2, Charge_Vs_PathlengthTECP1, Charge_Vs_PathlengthTECP2, Charge_Vs_PathlengthTIB, Charge_Vs_PathlengthTIDM, Charge_Vs_PathlengthTIDP, Charge_Vs_PathlengthTOB, gather_cfg::cout, dbe, MonitorElement::getTH2F(), harvestingMode, TFileService::make(), cppFunctionSkipper::operator, storeOnTree(), and tfs.

                                      {
    if(AlgoMode == "PCL" && !harvestingMode)return;//nothing to do in that case
 
    if(AlgoMode=="CalibTree"){
       // Loop on calibTrees to fill the 2D histograms
       algoAnalyzeTheTree();
    }else if(harvestingMode){
      // Load the 2D histograms from the DQM objects
      // When running in AlCaHarvesting mode the histos are already booked and should be just retrieved from
      // DQMStore so that they can be used in the fit
      cout << "   retrieving from DQMStore" << endl;
      Charge_Vs_Index           = dbe->get("AlCaReco/SiStripGains/Charge_Vs_Index");
      Charge_Vs_Index_Absolute  = dbe->get("AlCaReco/SiStripGains/Charge_Vs_Index_Absolute");
      Charge_Vs_PathlengthTIB   = dbe->get("AlCaReco/SiStripGains/Charge_Vs_PathlengthTIB");
      Charge_Vs_PathlengthTOB   = dbe->get("AlCaReco/SiStripGains/Charge_Vs_PathlengthTOB");
      Charge_Vs_PathlengthTIDP  = dbe->get("AlCaReco/SiStripGains/Charge_Vs_PathlengthTIDP");
      Charge_Vs_PathlengthTIDM  = dbe->get("AlCaReco/SiStripGains/Charge_Vs_PathlengthTIDM");
      Charge_Vs_PathlengthTECP1 = dbe->get("AlCaReco/SiStripGains/Charge_Vs_PathlengthTECP1");
      Charge_Vs_PathlengthTECP2 = dbe->get("AlCaReco/SiStripGains/Charge_Vs_PathlengthTECP2");
      Charge_Vs_PathlengthTECM1 = dbe->get("AlCaReco/SiStripGains/Charge_Vs_PathlengthTECM1");
      Charge_Vs_PathlengthTECM2 = dbe->get("AlCaReco/SiStripGains/Charge_Vs_PathlengthTECM2");
    }
 
    // Now that we have the full statistics we can extract the information of the 2D histograms
    algoComputeMPVandGain();
    
    //FIXME: for the moment the tree is disabled in PCL, among other reasons the fact that the TFileSevice is not available @ Tier0
    if(AlgoMode != "PCL"){
       storeOnTree();
       //also save the 2D monitor elements to this file as TH2D tfs
       tfs = edm::Service<TFileService>().operator->();
       tfs->make<TH2F> (*Charge_Vs_Index->getTH2F());
       tfs->make<TH2F> (*Charge_Vs_Index_Absolute->getTH2F());
       tfs->make<TH2F> (*Charge_Vs_PathlengthTIB->getTH2F());
       tfs->make<TH2F> (*Charge_Vs_PathlengthTOB->getTH2F());
       tfs->make<TH2F> (*Charge_Vs_PathlengthTIDP->getTH2F());
       tfs->make<TH2F> (*Charge_Vs_PathlengthTIDM->getTH2F());
       tfs->make<TH2F> (*Charge_Vs_PathlengthTECP1->getTH2F());
       tfs->make<TH2F> (*Charge_Vs_PathlengthTECP2->getTH2F());
       tfs->make<TH2F> (*Charge_Vs_PathlengthTECM1->getTH2F());
       tfs->make<TH2F> (*Charge_Vs_PathlengthTECM2->getTH2F());
    }
 }

SiStripApvGain * SiStripGainFromCalibTree::getNewObject ( )

overrideprivatevirtual

Implements ConditionDBWriter< SiStripApvGain >.

Definition at line 741 of file SiStripGainFromCalibTree.cc.

References a, sistrip::APV, APVsCollOrdered, gather_cfg::cout, stAPVGain::DetId, stAPVGain::Gain, harvestingMode, NULL, getGTfromDQMFile::obj, produceTagFilter(), SiStripApvGain::put(), and ConditionDBWriter< SiStripApvGain >::setDoStore().

 {
    SiStripApvGain* obj = new SiStripApvGain();
    if(!harvestingMode) return obj;
 
    if(!produceTagFilter()){
        cout << "[SiStripGainFromCalibTree] getNewObject -> will not produce a paylaod because produceTagFilter returned false " << endl;       
        setDoStore(false); 
        return obj;
    }
 
 
    std::vector<float>* theSiStripVector = NULL;
    unsigned int PreviousDetId = 0; 
    for(unsigned int a=0;a<APVsCollOrdered.size();a++)
    {
       stAPVGain* APV = APVsCollOrdered[a];
       if(APV==NULL){ printf("Bug\n"); continue; }
       if(APV->DetId != PreviousDetId){
          if(theSiStripVector!=NULL){
         SiStripApvGain::Range range(theSiStripVector->begin(),theSiStripVector->end());
         if ( !obj->put(PreviousDetId,range) )  printf("Bug to put detId = %i\n",PreviousDetId);
      }
      theSiStripVector = new std::vector<float>;
          PreviousDetId = APV->DetId;
       }
       theSiStripVector->push_back(APV->Gain);
    }
    if(theSiStripVector!=NULL){
       SiStripApvGain::Range range(theSiStripVector->begin(),theSiStripVector->end());
       if ( !obj->put(PreviousDetId,range) )  printf("Bug to put detId = %i\n",PreviousDetId);
    }
 
    return obj;
 }

void SiStripGainFromCalibTree::getPeakOfLandau	(	TH1 *	InputHisto,
		double *	FitResults,
		double	LowRange = `50`,
		double	HighRange = `5400`
	)

private

Definition at line 380 of file SiStripGainFromCalibTree.cc.

References MinNrEntries.

Referenced by algoComputeMPVandGain().

 { 
    FitResults[0]         = -0.5;  //MPV
    FitResults[1]         =  0;    //MPV error
    FitResults[2]         = -0.5;  //Width
    FitResults[3]         =  0;    //Width error
    FitResults[4]         = -0.5;  //Fit Chi2/NDF
 
    if( InputHisto->GetEntries() < MinNrEntries)return;
 
    // perform fit with standard landau
    TF1* MyLandau = new TF1("MyLandau","landau",LowRange, HighRange);
    MyLandau->SetParameter(1,300);
    InputHisto->Fit(MyLandau,"0QR WW");
 
    // MPV is parameter 1 (0=constant, 1=MPV, 2=Sigma)
    FitResults[0]         = MyLandau->GetParameter(1);  //MPV
    FitResults[1]         = MyLandau->GetParError(1);   //MPV error
    FitResults[2]         = MyLandau->GetParameter(2);  //Width
    FitResults[3]         = MyLandau->GetParError(2);   //Width error
    FitResults[4]         = MyLandau->GetChisquare() / MyLandau->GetNDF();  //Fit Chi2/NDF
 
    delete MyLandau;
 }

bool SiStripGainFromCalibTree::IsFarFromBorder	(	TrajectoryStateOnSurface *	trajState,
		const uint32_t	detid,
		const edm::EventSetup *	iSetup
	)

private

Definition at line 975 of file SiStripGainFromCalibTree.cc.

References Surface::bounds(), gather_cfg::cout, edm::EventSetup::get(), Bounds::length(), TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localPosition(), Parameters::parameters, LocalTrajectoryError::positionError(), GeomDet::surface(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

Referenced by algoAnalyze().

 { 
   edm::ESHandle<TrackerGeometry> tkGeom; iSetup->get<TrackerDigiGeometryRecord>().get( tkGeom );
 
   LocalPoint  HitLocalPos   = trajState->localPosition();
   LocalError  HitLocalError = trajState->localError().positionError() ;
 
   const GeomDetUnit* it = tkGeom->idToDetUnit(DetId(detid));
   if (dynamic_cast<const StripGeomDetUnit*>(it)==0 && dynamic_cast<const PixelGeomDetUnit*>(it)==0) {
      std::cout << "this detID doesn't seem to belong to the Tracker" << std::endl;
      return false;
   }
 
   const BoundPlane plane = it->surface();
   const TrapezoidalPlaneBounds* trapezoidalBounds( dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
   const RectangularPlaneBounds* rectangularBounds( dynamic_cast<const RectangularPlaneBounds*>(&(plane.bounds())));
 
   double DistFromBorder = 1.0;    
   double HalfLength     = it->surface().bounds().length() /2.0;
 
   if(trapezoidalBounds)
   {
      std::array<const float, 4> const & parameters = (*trapezoidalBounds).parameters();
      //std::vector<float> const & parameters = (*trapezoidalBounds).parameters();
      HalfLength     = parameters[3];
   }else if(rectangularBounds){
      HalfLength     = it->surface().bounds().length() /2.0;
   }else{return false;}
 
   if (fabs(HitLocalPos.y())+HitLocalError.yy() >= (HalfLength - DistFromBorder) ) return false;
 
   return true;
 }

bool SiStripGainFromCalibTree::IsGoodLandauFit ( double * FitResults )

private

Definition at line 405 of file SiStripGainFromCalibTree.cc.

Referenced by algoComputeMPVandGain().

                                                                 {
    if(FitResults[0] <= 0             )return false;
 //   if(FitResults[1] > MaxMPVError   )return false;
 //   if(FitResults[4] > MaxChi2OverNDF)return false;
    return true;   
 }

void SiStripGainFromCalibTree::MakeCalibrationMap ( )

private

Definition at line 782 of file SiStripGainFromCalibTree.cc.

References sistrip::APV, APVsColl, stAPVGain::CalibGain, m_calibrationPath, and useCalibration.

Referenced by algoBeginRun().

                                                  {
    if(!useCalibration)return;
   
    TChain* t1 = new TChain("SiStripCalib/APVGain");
    t1->Add(m_calibrationPath.c_str());
 
    unsigned int  tree_DetId;
    unsigned char tree_APVId;
    double        tree_Gain;
 
    t1->SetBranchAddress("DetId"             ,&tree_DetId      );
    t1->SetBranchAddress("APVId"             ,&tree_APVId      );
    t1->SetBranchAddress("Gain"              ,&tree_Gain       );
 
    for (unsigned int ientry = 0; ientry < t1->GetEntries(); ientry++) {
       t1->GetEntry(ientry);
        stAPVGain* APV = APVsColl[(tree_DetId<<3) | (unsigned int)tree_APVId];
        APV->CalibGain = tree_Gain;
    }
 }

bool SiStripGainFromCalibTree::produceTagFilter ( )

private

Definition at line 727 of file SiStripGainFromCalibTree.cc.

References BAD, Charge_Vs_Index, gather_cfg::cout, MonitorElement::getTH2F(), GOOD, tagCondition_GoodFrac, and tagCondition_NClusters.

Referenced by getNewObject().

                                                {
   
    // The goal of this function is to check wether or not there is enough statistics to produce a meaningful tag for the DB or not 
   if(Charge_Vs_Index->getTH2F()->Integral() < tagCondition_NClusters) {
     cout << "[SiStripGainFromCalibTree] produceTagFilter -> Return false: Statistics is too low : " << Charge_Vs_Index->getTH2F()->Integral() << endl;
     return false;
   }
   if((1.0 * GOOD) / (GOOD+BAD) < tagCondition_GoodFrac) {
     cout << "[SiStripGainFromCalibTree] produceTagFilter ->  Return false: ratio of GOOD/TOTAL is too low: " << (1.0 * GOOD) / (GOOD+BAD) << endl;
     return false;
   }
   return true; 
 }

void SiStripGainFromCalibTree::storeOnTree ( )

private

Definition at line 626 of file SiStripGainFromCalibTree.cc.

References a, AlgoMode, sistrip::APV, stAPVGain::APVId, APVsColl, APVsCollOrdered, BAD, Charge_Vs_Index, stAPVGain::DetId, stAPVGain::Eta, stAPVGain::FitChi2, stAPVGain::FitMPV, stAPVGain::FitMPVErr, stAPVGain::FitWidth, stAPVGain::FitWidthErr, stAPVGain::Gain, MonitorElement::getTH2F(), GOOD, stAPVGain::Index, stAPVGain::isMasked, TFileService::make(), NCluster, stAPVGain::NEntries, NEvent, NTrack, NULL, cppFunctionSkipper::operator, OutputGains, stAPVGain::Phi, stAPVGain::PreviousGain, stAPVGain::R, stAPVGain::SubDet, tfs, stAPVGain::Thickness, stAPVGain::x, stAPVGain::y, and stAPVGain::z.

Referenced by algoEndJob().

 {
   tfs = edm::Service<TFileService>().operator->();
 
    unsigned int  tree_Index;
    unsigned int  tree_Bin;
    unsigned int  tree_DetId;
    unsigned char tree_APVId;
    unsigned char tree_SubDet;
    float         tree_x;
    float         tree_y;
    float         tree_z;
    float         tree_Eta;
    float         tree_R;
    float         tree_Phi;
    float         tree_Thickness;
    float         tree_FitMPV;
    float         tree_FitMPVErr;
    float         tree_FitWidth;
    float         tree_FitWidthErr;
    float         tree_FitChi2NDF;
    double        tree_Gain;
    double        tree_PrevGain;
    double        tree_NEntries;
    bool          tree_isMasked;
 
    TTree*         MyTree;
    MyTree = tfs->make<TTree> ("APVGain","APVGain");
    MyTree->Branch("Index"             ,&tree_Index      ,"Index/i");
    MyTree->Branch("Bin"               ,&tree_Bin        ,"Bin/i");
    MyTree->Branch("DetId"             ,&tree_DetId      ,"DetId/i");
    MyTree->Branch("APVId"             ,&tree_APVId      ,"APVId/b");
    MyTree->Branch("SubDet"            ,&tree_SubDet     ,"SubDet/b");
    MyTree->Branch("x"                 ,&tree_x          ,"x/F"); 
    MyTree->Branch("y"                 ,&tree_y          ,"y/F");   
    MyTree->Branch("z"                 ,&tree_z          ,"z/F");   
    MyTree->Branch("Eta"               ,&tree_Eta        ,"Eta/F");
    MyTree->Branch("R"                 ,&tree_R          ,"R/F");
    MyTree->Branch("Phi"               ,&tree_Phi        ,"Phi/F");
    MyTree->Branch("Thickness"         ,&tree_Thickness  ,"Thickness/F");
    MyTree->Branch("FitMPV"            ,&tree_FitMPV     ,"FitMPV/F");
    MyTree->Branch("FitMPVErr"         ,&tree_FitMPVErr  ,"FitMPVErr/F");
    MyTree->Branch("FitWidth"          ,&tree_FitWidth   ,"FitWidth/F");
    MyTree->Branch("FitWidthErr"       ,&tree_FitWidthErr,"FitWidthErr/F");
    MyTree->Branch("FitChi2NDF"        ,&tree_FitChi2NDF ,"FitChi2NDF/F");
    MyTree->Branch("Gain"              ,&tree_Gain       ,"Gain/D");
    MyTree->Branch("PrevGain"          ,&tree_PrevGain   ,"PrevGain/D");
    MyTree->Branch("NEntries"          ,&tree_NEntries   ,"NEntries/D");
    MyTree->Branch("isMasked"          ,&tree_isMasked   ,"isMasked/O");
 
 
    FILE* Gains = stdout;
    if(AlgoMode!="PCL"){
    fprintf(Gains,"NEvents   = %i\n",NEvent);
    fprintf(Gains,"NTracks   = %i\n",NTrack);
    fprintf(Gains,"NClusters = %i\n",NCluster);
    fprintf(Gains,"Number of APVs = %lu\n",static_cast<unsigned long>(APVsColl.size()));
    fprintf(Gains,"GoodFits = %i BadFits = %i ratio = %f\n",GOOD,BAD,(100.0*GOOD)/(GOOD+BAD));
    Gains=fopen(OutputGains.c_str(),"w");
    }
    fprintf(Gains,"NEvents   = %i\n",NEvent);
    fprintf(Gains,"NTracks   = %i\n",NTrack);
    fprintf(Gains,"NClusters = %i\n",NCluster);
    fprintf(Gains,"Number of APVs = %lu\n",static_cast<unsigned long>(APVsColl.size()));
    fprintf(Gains,"GoodFits = %i BadFits = %i ratio = %f\n",GOOD,BAD,(100.0*GOOD)/(GOOD+BAD));
 
    for(unsigned int a=0;a<APVsCollOrdered.size();a++){
       stAPVGain* APV = APVsCollOrdered[a];
       if(APV==NULL)continue;
 //     printf(      "%i | %i | PreviousGain = %7.5f NewGain = %7.5f (#clusters=%8.0f)\n", APV->DetId,APV->APVId,APV->PreviousGain,APV->Gain, APV->NEntries);
       fprintf(Gains,"%i | %i | PreviousGain = %7.5f NewGain = %7.5f (#clusters=%8.0f)\n", APV->DetId,APV->APVId,APV->PreviousGain,APV->Gain, APV->NEntries);
 
       tree_Index      = APV->Index;
       tree_Bin        = Charge_Vs_Index->getTH2F()->GetXaxis()->FindBin(APV->Index);
       tree_DetId      = APV->DetId;
       tree_APVId      = APV->APVId;
       tree_SubDet     = APV->SubDet;
       tree_x          = APV->x;
       tree_y          = APV->y;
       tree_z          = APV->z;
       tree_Eta        = APV->Eta;
       tree_R          = APV->R;
       tree_Phi        = APV->Phi;
       tree_Thickness  = APV->Thickness;
       tree_FitMPV     = APV->FitMPV;
       tree_FitMPVErr  = APV->FitMPVErr;
       tree_FitWidth   = APV->FitWidth;
       tree_FitWidthErr= APV->FitWidthErr;
       tree_FitChi2NDF = APV->FitChi2;
       tree_Gain       = APV->Gain;
       tree_PrevGain   = APV->PreviousGain;
       tree_NEntries   = APV->NEntries;
       tree_isMasked   = APV->isMasked;
 
       MyTree->Fill();
    }
    if(AlgoMode!="PCL")fclose(Gains);
 
 
 }